The near-edge structure in energy-loss spectroscopy: many-electron and magnetic effects in transition metal nitrides and carbides

We investigate the ability of the local density approximation (LDA) in density functional theory to predict the near-edge structure in electron energy-loss spectroscopy in the dipole approximation. We include screening of the core hole within the LDA using Slater's transition state theory. We find that anion K-edge threshold energies are systematically overestimated by 4.22 +/- 0.44 eV in twelve transition metal carbides and nitrides in the rock-salt (B1) structure. When we apply this 'universal' many-electron correction to energy-loss spectra calculated within the transition state approximation to LDA, we find quantitative agreement with experiment to within one or two eV for TiC, TiN and VN. We compare our calculations to a simpler approach using a projected Mulliken density which honours the dipole selection rule, in place of the dipole matrix element itself. We find remarkably close agreement between these two approaches. Finally, we show an anomaly in the near-edge structure in CrN to be due to magnetic structure. In particular, we find that the N K edge in fact probes the magnetic moments and alignments of ther sublattice.

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Hydrogen bonding and collective proton modes in clusters and periodic layers of squaric acid: A density functional study

Hydrogen bonding in clusters and extended layers of squaric acid molecules has been investigated by density functional computations. Equilibrium geometries, harmonic vibrational frequencies, and energy barriers for proton transfer along hydrogen bonds have been determined using the Car-Parrinello method. The results provide crucial parameters for a first principles modeling of the potential energy surface, and highlight the role of collective modes in the low-energy proton dynamics. The importance of quantum effects in condensed squaric acid systems has been investigated, and shown to be negligible for the lowest-energy collective proton modes. This information provides a quantitative basis for improved atomistic models of the order-disorder and displacive transitions undergone by squaric acid crystals as a function of temperature and pressure. (C) 2001 American Institute of Physics.

Local and semilocal density functional computations for crystals of 1-alkyl-3-methyl-imidazolium salts

The accuracy and reliability of popular density functional approximations for the compounds giving origin to room temperature ionic liquids have been assessed by computing the T=0 K crystal structure of several 1-alkyl-3-methyl-imidazolium salts. Two prototypical exchange-correlation approximations have been considered, i.e., the local density approximation (LDA) and one gradient corrected scheme [PBE-GGA, Phys. Rev. Lett. 77, 3865 (1996)]. Comparison with low-temperature x-ray diffraction data shows that the equilibrium volume predicted by either approximations is affected by large errors, nearly equal in magnitude (~10%), and of opposite sign. In both cases the error can be traced to a poor description of the intermolecular interactions, while the intramolecular structure is fairly well reproduced by LDA and PBE-GGA. The PBE-GGA optimization of atomic positions within the experimental unit cell provides results in good agreement with the x-ray structure. The correct system volume can also be restored by supplementing PBE-GGA with empirical dispersion terms reproducing the r-6 attractive tail of the van der Waals interactions.

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.

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 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.