Effect of spatial nonlocality on the density functional band gap

We show that for a large class of exchange-correlation functionals the local exchange-correlation potential obtained within an optimized effective potential severely underestimates the band gap. On the other hand, the corresponding nonlocal potential obtained from a generalized Kohn-Sham scheme provides a much better description of the band gap, in good agreement with experiments. These results strongly indicate that a local exchange-correlation potential, however good the exchange-correlation approximation, cannot capture the delicate interplay between correlation effects and spatial localization in the KS band structure, unless the (cumbersome) contribution from the derivative discontinuity of the exchange-correlation energy functional is considered.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Catalytic Water Formation on Platinum: A First-Principles Study

The study of catalytic behavior begins with one seemingly simple process, namely the hydrogenation of O to H2O on platinum. Despite the apparent simplicity its mechanism has been much debated. We have used density functional theory with,gradient corrections to examine microscopic reaction pathways for several elementary steps implicated in this fundamental catalytic process. We find that H2O formation from chemisorbed O and H atoms is a highly activated process. The largest barrier along this route, with a value of similar to1 eV, is the addition of the first H to O to produce OH. Once formed, however, OH groups are easily hydrogenated to H2O with a barrier of similar to0.2 eV. Disproportionation reactions with 1:1 and 2:1 stoichiometries of H2O and O have been examined as alternative routes for OH formation. Both stoichiometries of reaction produce OH groups with barriers that are much lower than that associated with the O + H reaction. H2O, therefore, acts as an autocatalyst in the overall H O formation process. Disproportionation with a 2:1 stoichiometry is thermodynamically and kinetically favored over disproportionation with a l:I stoichiometry. This highlights an additional (promotional) role of the second H2O molecule in this process. In support of our previous suggestion that the key intermediate in the low-temperature H2O formation reaction is a mixed OH and H2O overlayer we find that then is a very large barrier for the dissociation of the second H2O molecule in the 2:1 disproportionation process. We suggest that the proposed intermediate is then hydrogenated to H2O through a very facile proton transfer mechanism.

DFT and in situ EXAFS investigation of gold/ceria-zirconia low-temperature water gas shift catalysts: Identification of the nature of the active form of gold

A combined experimental and theoretical investigation of the nature of the active form of gold in oxide-supported gold catalysts for the water gas shift reaction has been performed. In situ extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES) experiments have shown that in the fresh catalysts the gold is in the form of highly dispersed gold ions. However, under water gas shift reaction conditions, even at temperatures as low as 100 degrees C, the evidence from EXAFS and XANES is only 14 consistent with rapid, and essentially complete, reduction of the gold to form metallic clusters containing about 50 atoms. The presence of Au-Ce distances in the EXAFS spectra, and the fact that about 15% of the gold atoms can be reoxidized after exposure to air at 150 degrees C, is indicative of a close interaction between a fraction (ca. 15%) of the gold atoms and the oxide support. Density functional theory (DFT) calculations are entirely consistent with this model and suggest that an important aspect of the active and stable form of gold under water gas shift reaction conditions is the location of a partially oxidized gold (Audelta+) species at a cerium cation vacancy in the surface of the oxide support. It is found that even with a low loading gold catalysts (0.2%) the fraction of ionic gold under water gas shift conditions is below the limit of detection by XANES (<5%). It is concluded that under water gas shift reaction conditions the active form of gold comprises small metallic gold clusters in intimate contact with the oxide support.

On the solvation of L-aspartic acid

We use molecular statics and dynamics to study the stability of L-aspartic acid both in vacuo and solvated by polar and non-polar molecules using density functional theory in the generalized gradient approximation. We find that structures stable in vacuo are unstable in aqueous solution and vice versa. From our simulations we are able to come to some conclusions about the mechanism of stabilisation of zwitterions by polar protic solvents, water and methanol.

Bandstructure approach to near edge structure

We review the current state of the art in EELS fingerprinting by computer simulation, focusing on the bandstructure approach to the problem. Currently calculations are made using a one electron theory, but we describe in principle the way to go beyond this to include final state effects. We include these effects within the one electron framework using the Slater transition state formula and assess the errors involved. Two examples are then given which illustrate the use of the one electron approximation within density functional theory. Our approach is to combine predicted atomic structure with predicted electronic structure to assist in fingerprinting of complex crystal structures.

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.

Fixed-node diffusion Monte Carlo computations for closed-shell jellium spheres

Fixed-node diffusion Monte Carlo computations are used to determine the ground state energy and electron density for jellium spheres with up to N = 106 electrons and background densities corresponding to the electron gas parameter 1 less than or equal to r(s)less than or equal to5.62. We analyze the density and size dependence of the surface energy, and we extrapolate our data to the thermodynamic limit. The results agree well with the predictions of density functional computations using the local density approximation. In the case of N = 20, we extend our computation to higher densities and identify a transition between atomic- and jelliumlike nodal structures occurring at the background density corresponding to r(s)=0.13. In this case the local density approximation is unable to reproduce the changes in the correlation energy due to the discontinuous transition in the ground state nodal structure. We discuss the relevance of our results for nonlocal approximations to density functional theory.

Thin Ta films: growth, stability, and diffusion studied by molecular-dynamics simulations

We report results of classical molecular-dynamics simulations of bcc and beta-Ta thin films. Thermal PVD film growth, surface roughness, argon ion bombardment, phase stability and transformation, vacancy and adatom diffusion, and thermal relaxation kinetics are discussed. Distinct differences between the two structures are observed, including a complex vacancy diffusion mechanism in beta-Ta. Embedded atom method potentials, which were fitted to bcc properties, have been used to model the Ta-Ta interactions. In order to verify the application of these potentials to the more complex beta-Ta structure, we have also performed density functional theory calculations. Results and implications of these calculations are discussed.

Catalytic Role of Metal Oxides in Gold-Based Catalysts: A First Principles Study of CO Oxidation on TiO2 Supported Au

CO oxidation on TiO2 supported Au has been studied using density functional theory calculations. Important catalytic roles of the oxide have been identified: (i) CO oxidation occurs at the interface between Au and the oxide with a very small barrier; and (ii) O-2 adsorption at the interface is the key step in the reaction. The physical origin of the oxide promotion effect has been further investigated: The oxide enhances electron transfer from the Au to the antibonding states of O-2, giving rise to (i) strong ionic bonding between the adsorbed O-2, Au, and the Ti cation; and (ii) a significant activation of O-2 towards CO oxidation.

Molecular effects in the ionization of N-2, O-2, and F-2 by intense laser fields

In this paper we study the response in time of N2, O2, and F2 to laser pulses having a wavelength of 390 nm. We find single-ionization suppression in O2 and its absence in F2, in accordance with experimental results at lambda= 800 nm. Within our framework of time-dependent density functional theory we are able to explain deviations from the predictions of intense-field many-body S-matrix theory (IMST). We confirm the connection of ionization suppression with destructive interference of outgoing electron waves from the ionized electron orbital. However, the prediction of ionization suppression, justified within the IMST approach through the symmetry of the highest occupied molecular orbital (HOMO), is not reliable since it turns out that—e.g., in the case of F2—the electronic response to the laser pulse is rather complicated and does not lead to dominant depletion of the HOMO. Therefore, the symmetry of the HOMO is not sufficient to predict ionization suppression. However, at least for F2, the symmetry of the dominantly ionized orbital is consistent with the nonsuppression of ionization.

Transport in nanoscale systems: the microcanonical versus grand-canonical picture

We analyse a picture of transport in which two large but finite charged electrodes discharge across a nanoscale junction. We identify a functional whose minimization, within the space of all bound many-body wavefunctions, defines an instantaneous steady state. We also discuss factors that favour the onset of steady-state conduction in such systems, make a connection with the notion of entropy, and suggest a novel source of steady-state noise. Finally, we prove that the true many-body total current in this closed system is given exactly by the one-electron total current, obtained from time-dependent density-functional theory.

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.