Magnetic tight binding and the iron-chromium enthalpy anomaly

We describe a self-consistent magnetic tight-binding theory based in an expansion of the Hohenberg-Kohn density functional to second order, about a non-spin-polarized reference density. We show how a first order expansion about a density having a trial input magnetic moment leads to a fixed moment model. We employ a simple set of tight-binding parameters that accurately describes electronic structure and energetics, and show these to be transferable between first row transition metals and their alloys. We make a number of calculations of the electronic structure of dilute Cr impurities in Fe, which we compare with results using the local spin density approximation. The fixed moment model provides a powerful means for interpreting complex magnetic configurations in alloys; using this approach, we are able to advance a simple and readily understood explanation for the observed anomaly in the enthalpy of mixing.

Spontaneous spin polarization and electron localization in constrained geometries: The Wigner transition in nanowires

The Wigner transition in a jellium model of cylindrical nanowires has been investigated by density-functional computations using the local spin-density approximation. A wide range of background densities rho(b) has been explored from the nearly ideal metallic regime (r(s)=[3/4 pi rho(b)](1/3)=1) to the high correlation limit (r(s)=100). Computations have been performed using an unconstrained plane wave expansion for the Kohn-Sham orbitals and a large simulation cell with up to 480 electrons. The electron and spin distributions retain the cylindrical symmetry of the Hamiltonian at high density, while electron localization and spin polarization arise nearly simultaneously in low-density wires (r(s)similar to 30). At sufficiently low density (r(s)>= 40), the ground-state electron distribution is the superposition of well defined and nearly disjoint droplets, whose charge and spin densities integrate almost exactly to one electron and 1/2 mu(B), respectively. Droplets are arranged on radial shells and define a distorted lattice whose structure is intermediate between bcc and fcc. Dislocations and grain boundaries are apparent in the droplets' configuration found by our simulations. Our computations aim at modeling the behavior of experimental low-carried density systems made of lightly doped semiconductor nanostructures or conducting polymers.

Electronic structure of the nitrogen-vacancy center in diamond from first-principles theory

The nitrogen-vacancy (NV) center is a paramagnetic defect in diamond with applications as a qubit. Here, we investigate its electronic structure by using ab initio density functional theory for five different NV center models of two different cluster sizes. We describe the symmetry and energetics of the low-lying states and compare the optical frequencies obtained to experimental results. We compute the major transition of the negatively charged NV centers to within 25–100 meV accuracy and find that it is energetically favorable for substitutional nitrogens to donate an electron to NV0. The excited state of the major transition and the NV0 state with a neutral donor nitrogen are found to be close in energy.

The energetics of tetrahydrocarbazole aromatization over Pd(111): A computational analysis

The carbazole moiety is a component of many important pharmaceuticals including anticancer and anti-HIV agents and is commonly utilized in the production of modern polymeric materials with novel photophysical and electronic properties. Simple carbazoles are generally produced via the aromatization of the respective tetrahydrocarbazole (THCZ). In this work, density functional theory calculations are used to model the reaction pathway of tetrahydrocarbazole aromatization over Pd(111). The geometry of each of the intermediate surface species has been determined and how each structure interacts with the metal surface addressed. The reaction energies and barriers of each of the elementary surface reactions have also been calculated, and a detailed analysis of the energetic trends performed. Our calculations have shown that the surface intermediates remain fixed to the surface via the aromatic ring in a manner similar to that of THCZ. Moreover, the aliphatic ring becomes progressively more planer with the dissociation of each subsequent hydrogen atom. Analysis of the reaction energy profile has revealed that the trend in reaction barriers is determined by the two factors: (i) the strength of the dissociating ring-H bond and (ii) the subsequent gain in energy due to the geometric relaxation of the aliphatic ring. (c) 2008 American Institute of Physics.

Dramatic liquid-phase dehydrogenation rate enhancements using gas-phase hydrogen acceptors

The dehydrogenation of 1,2,3,4-tetrahydrocarbazole (THCZ) to form carbazole (CZ) over supported palladium catalysts was examined in the presence of hydrogen acceptors. As expected, liquid hydrogen acceptors increased the rate of reaction but, importantly, gaseous hydrogen acceptors also have been used. Ethene, propene, and but-1-ene showed up to a fivefold increase in the rate of dehydrogenation. Moreover, compared with the analogous liquid systems, the gaseous alternatives are a potentially more economic method of enhancing the activity and provide a simpler workup. The mechanism for the increase in rate was examined by density functional theory calculations, which showed that the propene hydrogenation competes effectively with the back-hydrogenation of the intermediates formed during the THCZ dehydrogenation, resulting in a shift in the equilibrium toward to the formation of CZ. (C) 2007 Elsevier Inc. All rights reserved.

A DFT study of the chain growth probability in Fischer-Tropsch synthesis

The chain growth probability (alpha value) is one of the most significant parameters in Fischer-Tropsch (FT) synthesis. To gain insight into the chain growth probability, we systematically studied the hydrogenation and C-C coupling reactions with different chain lengths on the stepped Co(0001) surface using density functional theory calculations. Our findings elucidate the relationship between the barriers of these elementary reactions and the chain length. Moreover, we derived a general expression of the chain growth probability and investigated the behavior of the alpha value observed experimentally. The high methane yield results from the lower chain growth rate for C-1 + C-1 coupling compared with the other coupling reactions. After C-1, the deviation of product distribution in FT synthesis from the Anderson-Schulz-Flory distribution is due to the chain length-dependent paraffin/olefin ratio. (C) 2008 Elsevier Inc. All rights reserved.

Insight into the adsorption competition and the relationship between dissociation and association reactions in ammonia synthesis

Ammonia synthesis on three metal surfaces (Zr, Ru, and Pd) is investigated using density functional theory calculations. In addition to N-2 dissociation, all the transition states of the hydrogenation reactions from N to NH3 are located and the reaction energy profiles at both low and high surface coverages are compared and analyzed. The following are found: (i) Surface coverage effect on dissociation reactions is more significant than that on association reactions. (ii) The difference between N and H chemisorption energies, the so-called chemisorption energy gap which is a measure of adsorption competition, is vital to the reactivity of the catalysts. (iii) The hydrogenation barriers can considerably affect the overall rate of ammonia synthesis. A simple model to describe the relationship between dissociation and association reactions is proposed. (c) 2007 American Institute of Physics.

Trends in C-O and C-N bond formations over transition metal surfaces: An insight into kinetic sensitivity in catalytic reactions

Transition metal catalyzed bond formation is a fundamental process in catalysis and is of general interest throughout chemistry. To date, however, the knowledge of association reactions is rather limited, relative to what is known about dissociative processes. For example, surprisingly little is known about how the bond-forming ability of a metal, in general, varies across the Periodic Table. In particular, the effect of reactant valency on such trends is poorly understood. Herein, the authors examine these key issues by using density functional theory calculations to study CO and CN formations over the 4d metals. The calculations reveal that the chemistries differ in a fundamental way. In the case of CO formation, the reaction enthalpies span a much greater range than those of CN formation. Moreover, CO formation is found to be kinetically sensitive to the metal; here the reaction barriers (E-a) are found to be influenced by the reaction enthalpy. CN formation, conversely, is found to be relatively kinetically insensitive to the metal, and there is no correlation found between the reaction barriers and the reaction enthalpy. Analysis has shown that at the final adsorbed state, the interaction between N and the surface is relatively greater than that of O. Furthermore, in comparison with O, relatively less bonding between the surface and N is observed to be lost during transition state formation. These greater interactions between N and the surface, which can be related to the larger valency of N, are found to be responsible for the relatively smaller enthalpy range and limited variation in E-a for CN formation. (C) 2007 American Institute of Physics.

The dissociation of molecularly adsorbed CO and CN over the 4d transition metals: A universal relationship between the reaction barriers and the reaction enthalpies

Density functional theory calculations are used to study the stability of molecularly adsorbed CO and CN over transition metal surfaces. The minimum energy reaction pathways, corresponding reaction barriers (E-a), and reaction enthalpies (Delta H) for the dissociation of CO and CN to atomic products over the 4d transition metals from Zr to Pd have been determined. CO is found to be the more stable adsorbate on the right hand side of the period (from Tc onwards), whereas CN is the more stable surface species on the early metals (Zr, Nb and Mo). A single linear relationship is found to exist that correlates the barriers of both reactions with their respective reaction enthalpies. (c) 2006 Elsevier B.V. All rights reserved.

Importance of electronegativity differences and surface structure in molecular dissociation reactions at transition metal surfaces

The dissociative adsorption of N-2 has been studied at both monatomic steps and flat regions on the surfaces of the 4d transition metals from Zr to Pd. Using density functional theory (DFT) calculations, we have determined and analyzed the trends in both straight reactivity and structure sensitivity across the periodic table. With regards to reactivity, we find that the trend in activation energy (Ea) is determined mainly by a charge transfer from the surface metal atoms to the N atoms during transition state formation, namely, the degree of ionicity of the N-surface bond at the transition state. Indeed, we find that the strength of the metal-N bond at the transition state (and therefore the trend in Ea) can be predicted by the difference in Mulliken electronegativity between the metal and N. Structure sensitivity is analyzed in terms of geometric and electronic effects. We find that the lowering of Ea due to steps is more pronounced on the right-hand side of the periodic table. It is found that for the early transition metals the geometric and electronic effects work in opposition when going from terrace to step active site. In the case of the late 4d metals, however, these effects work in combination, producing a more marked reduction in Ea.

Reactivity of the 4d transition metals toward N hydrogenation and NH dissociation: A DFT-based HSAB analysis

The density functional theory (DFT) based hard-soft acid-base (HSAB) reactivity indices, including the electrophilicity index, have been successfully applied to many areas of molecular chemistry. In this work we test the applicability of such an approach to fundamental surface chemistry. We have considered, as prototypical surface reactions, both the hydrogenation of atomic nitrogen and the dissociative adsorption of the NH molecular radical. By use of a DFT methodology, the minimum energy reaction pathways, and corresponding reaction barriers, of the above reactions over Zr(001), Nb(110), Mo(110), Tc(001), Ru(001), Rh(111), and Pd(111) have been determined. By consideration of the chemical potential and chemical hardness of the surface metal atoms, and the principle of electronegativity equalization, it is found that the charge transferred to the NH radical during the process of dissociative adsorption correlates very well with that determined by Mulliken population analysis. Furthermore, it is found that the stability of the NH/surface transition state complex relates directly to this charge transfer and that the trend in transition state stability predicted by a HSAB; treatment correlates very strongly with that determined by DFT calculations. With regards to N hydrogenation, we find that during the course of the reaction, H loses cohesion to the surface, as it must migrate from a 3-fold hollow site to either a bridge or top site, to react with N. Partial density of states (PDOS) and Mulliken population analysis reveal that this loss of bonding is accompanied by charge transfer from H to the surface metal atoms. Moreover, by simple modeling, we show that the reaction barriers are directly proportional to this mandatory charge transfer. Indeed, it is found that the reaction barriers correlate very well with the electrophilicity index of the metal atoms.

The importance of hydrogen's potential-energy surface and the strength of the forming R-H bond in surface hydrogenation reactions

An understanding of surface hydrogenation reactivity is a prevailing issue in chemistry and vital to the rational design of future catalysts. In this density-functional theory study, we address hydrogenation reactivity by examining the reaction pathways for N+H -> NH and NH+H -> NH2 over the close-packed surfaces of the 4d transition metals from Zr-Pd. It is found that the minimum-energy reaction pathway is dictated by the ease with which H can relocate between hollow-site and top-site adsorption geometries. A transition state where H is close to a top site reduces the instability associated with bond sharing of metal atoms by H and N (NH) (bonding competition). However, if the energy difference between hollow-site and top-site adsorption energies (Delta E-H) is large this type of transition state is unfavorable. Thus we have determined that hydrogenation reactivity is primarily controlled by the potential-energy surface of H on the metal, which is approximated by Delta E-H, and that the strength of N (NH) chemisorption energy is of less importance. Delta E-H has also enabled us to make predictions regarding the structure sensitivity of these reactions. Furthermore, we have found that the degree of bonding competition at the transition state is responsible for the trend in reaction barriers (E-a) across the transition series. When this effect is quantified a very good linear correlation is found with E-a. In addition, we find that when considering a particular type of reaction pathway, a good linear correlation is found between the destabilizing effects of bonding competition at the transition state and the strength of the forming N-H (HN-H) bond. (c) 2006 American Institute of Physics.

The catalytic role of water in CO oxidation

Water, one of the most popular species in our planet, can play a catalytic role in many reactions, including reactions in heterogeneous catalysis. In a recent experimental work, Bergeld, Kasemo, and Chakarov demonstrated that water is able to promote CO oxidation under low temperatures (similar to200 K). In this study, we choose CO oxidation on Pt(111) in the presence of water as a model system to address the catalytic role of water for surface reactions in general using density functional theory. Many elementary steps possibly involved in the CO oxidation on Pt(111) at low temperatures have been investigated. We find the following. First, in the presence of water, the CO oxidation barrier is reduced to 0.33 eV (without water the barrier is 0.80 eV). This barrier reduction is mainly due to the H-bonding between the H in the H2O and the O at the transition state (TS), which stabilizes the TS. Second, CO can readily react with OH with a barrier of 0.44 eV, while COOH dissociation to produce CO2 is not easy (the barrier is 1.02 eV). Third, in the H2O+OH mixed phase, CO can be easily converted into CO2. It occurs through two steps: CO reacts with OH, forming COOH; and COOH transfers the H to a nearby H2O and, at the same time, an H in the H2O transfers to a OH, leading to CO2 formation. The reaction barrier of this process is 0.60 eV under CO coverage of 1/6 ML and 0.33 eV under CO coverage of 1/3 ML. The mechanism of CO oxidation at low temperatures is discussed. On the basis of our calculations, we propose that the water promotion effect can in general be divided into two classes: (i) By H-bonding between the H of H2O and an electron negative species such as the O in the reaction of CO+O+H2O-->CO2+H2O, H2O can stabilize the TS of the reaction and hence reduce the barrier. (ii) H2O first dissociates into H and OH and then OH or H participates directly in the reaction to induce new reaction mechanism with more favorable routes, in which OH or H can act as an intermediate. (C) 2003 American Institute of Physics.

Resolution of an ancient surface science anomaly: Work function change induced by N adsorption on W{100}

For many decades it has been assumed that an adsorbate centered above a metal surface and with a net negative charge should increase the work function of the surface. However, despite their electronegativity, N adatoms on W{100} cause a significant work function decrease. Here we present a resolution of this anomaly. Using density functional theory, we demonstrate that while the N atom carries a negative charge, of overriding importance is a reduction in the surface overspill electron density into the vacuum, when that charge is engaged in bonding to the adatom. This novel interpretation is fundamentally important in the general understanding of work function changes induced by atomic adsorbates.

Mechanism of H-2 metabolism on Fe-only hydrogenases

The metabolism of hydrogen (H-2 2H(+) + 2e(-)) constitutes a central process in the global biological energy cycle. Among all the enzymes that can mediate this process, Fe-only hydrogenases are unique in their particular high reactivity. Recently, some important progresses have been achieved. Following our previous paper [Z.-P. Liu and P. Hu, J. Am. Chem. Soc. 124, 5175 (2002)] that characterizes the individual redox state of the active site of Fe-only hydrogenase, in this work we have determined the feasible reaction pathways and energetics for the H-2 metabolism on the active site of Fe-only hydrogenases, using density functional theory. We show that H-2 metabolism possesses very low reaction barriers and a proximal base from a nearby protein plays an important role in H-2 metabolism. (C) 2002 American Institute of Physics.