Macroscopic limit of time-dependent density-functional theory for adiabatic local approximations of the exchange-correlation kernel

Time-dependent density-functional theory is a rather accurate and efficient way to compute electronic excitations for finite systems. However, in the macroscopic limit (systems of increasing size), for the usual adiabatic random-phase, local-density, or generalized-gradient approximations, one recovers the Kohn-Sham independent-particle picture, and thus the incorrect band gap. To clarify this trend, we investigate the macroscopic limit of the exchange-correlation kernel in such approximations by means of an algebraical analysis complemented with numerical studies of a one-dimensional tight-binding model. We link the failure to shift the Kohn-Sham spectrum of these approximate kernels to the fact that the corresponding operators in the transition space act only on a finite subspace.

Many-body theory of positron-atom interactions

A many-body theory approach is developed for the problem of positron-atom scattering and annihilation. Strong electron- positron correlations are included nonperturbatively through the calculation of the electron-positron vertex function. It corresponds to the sum of an infinite series of ladder diagrams, and describes the physical effect of virtual positronium formation. The vertex function is used to calculate the positron-atom correlation potential and nonlocal corrections to the electron-positron annihilation vertex. Numerically, we make use of B-spline basis sets, which ensures rapid convergence of the sums over intermediate states. We have also devised an extrapolation procedure that allows one to achieve convergence with respect to the number of intermediate- state orbital angular momenta included in the calculations. As a test, the present formalism is applied to positron scattering and annihilation on hydrogen, where it is exact. Our results agree with those of accurate variational calculations. We also examine in detail the properties of the large correlation corrections to the annihilation vertex.

C-H bond activation over metal oxides: A new insight into the dissociation kinetics from density functional theory

The C-H activation on metal oxides is a fundamental process in chemistry. In this paper, we report a density functional theory study on the process of the C-H activation of CH4 on Pd(111), Pt(111), Ru(0001), Tc(0001), Cu(111), PdO(001), PdO(110), and PdO(100). A linear relationship between the C-H activation barrier and the chemisorption in the dissociation final state on the metal surfaces is obtained, which is consistent with the work in the literature. However, the relationship is poor on the metal oxide surfaces. Instead, a strong linear correlation between the barrier and the lattice O-H bond strength is found on the oxides. The new linear relationship is analyzed and the physical origin is identified. (c) 2008 American Institute of Physics.

Many-body theory of gamma spectra from positron-atom annihilation

A many-body theory approach to the calculation of gamma spectra of positron annihilation on many-electron atoms is developed. We evaluate the first-order correlation correction to the annihilation vertex and perform numerical calculations for the noble gas atoms. Extrapolation with respect to the maximal orbital momentum of the intermediate electron and positron states is used to achieve convergence. The inclusion of correlation corrections improves agreement with experimental gamma spectra.

Many-body theory calculations of positron binding to negative ions

A many-body theory approach developed by the authors [Phys. Rev. A 70, 032720 (2004)] is applied to positron bound states and annihilation rates in atomic systems. Within the formalism, full account of virtual positronium (Ps) formation is made by summing the electron-positron ladder diagram series, thus enabling the theory to include all important many-body correlation effects in the positron problem. Numerical calculations have been performed for positron bound states with the hydrogen and halogen negative ions, also known as Ps hydride and Ps halides. The Ps binding energies of 1.118, 2.718, 2.245, 1.873 and 1.393 eV and annihilation rates of 2.544, 2.482, 1.984, 1.913 and 1.809 ns^{-1}, have been obtained for PsH, PsF, PsCl, PsBr and PsI, respectively.

Density functionals from many-body perturbation theory: The band gap for semiconductors and insulators

Theoretically the Kohn-Sham band gap differs from the exact quasiparticle energy gap by the derivative discontinuity of the exchange-correlation functional. In practice for semiconductors and insulators the band gap calculated within any local or semilocal density approximations underestimates severely the experimental energy gap. On the other hand, calculations with an "exact" exchange potential derived from many-body perturbation theory via the optimized effective potential suggest that improving the exchange-correlation potential approximation can yield a reasonable agreement between the Kohn-Sham band gap and the experimental gap. The results in this work show that this is not the case. In fact, we add to the exact exchange the correlation that corresponds to the dynamical (random phase approximation) screening in the GW approximation. This accurate exchange-correlation potential provides band structures similar to the local density approximation with the corresponding derivative discontinuity that contributes 30%-50% to the energy gap. Our self-consistent results confirm substantially the results for Si and other semiconductors obtained perturbatively [R. W. Godby , Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid. (c) 2006 American Institute of Physics.

Shape corrections to exchange-correlation potentials by gradient-regulated seamless connection of model potentials for inner and outer region

Shape corrections to the standard approximate Kohn-Sham exchange-correlation (xc) potentials are considered with the aim to improve the excitation energies (especially for higher excitations) calculated with time-dependent density functional perturbation theory. A scheme of gradient-regulated connection (GRAC) of inner to outer parts of a model potential is developed. Asymptotic corrections based either on the potential of Fermi and Amaldi or van Leeuwen and Baerends (LB) are seamlessly connected to the (shifted) xc potential of Becke and Perdew (BP) with the GRAC procedure, and are employed to calculate the vertical excitation energies of the prototype molecules N-2, CO, CH2O, C2H4, C5NH5, C6H6, Li-2, Na-2, K-2. The results are compared with those of the alternative interpolation scheme of Tozer and Handy as well as with the results of the potential obtained with the statistical averaging of (model) orbital potentials. Various asymptotically corrected potentials produce high quality excitation energies, which in quite a few cases approach the benchmark accuracy of 0.1 eV for the electronic spectra. Based on these results, the potential BP-GRAC-LB is proposed for molecular response calculations, which is a smooth potential and a genuine "local" density functional with an analytical representation. (C) 2001 American Institute of Physics.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Density functional theory study on the activation of molecular oxygen on a stepped gold surface in an aqueous environment: a new approach for simulating reactions in solution

The activation of oxygen molecules is an important issue in the gold-catalyzed partial oxidation of alcohols in aqueous solution. The complexity of the solution arising from a large number of solvent molecules makes it difficult to study the reaction in the system. In this work, O-2 activation on an Au catalyst is investigated using an effective approach to estimate the reaction barriers in the presence of solvent. Our calculations show that O-2 can be activated, undergoing OOH* in the presence of water molecules. The OOH* can readily be formed on Au(211) via four possible pathways with almost equivalent free energy barriers at the aqueous-solid interface: the direct or indirect activation of O-2 by surface hydrogen or the hydrolysis of O-2 following a Langmuir-Hinshelwood mechanism or an Eley-Rideal mechanism. Among them, the Eley-Rideal mechanism may be slightly more favorable due to the restriction of the low coverage of surface H on Au(211) in the other mechanisms. The results shed light on the importance of water molecules on the activation of oxygen in gold-catalyzed systems. Solvent is found to facilitate the oxygen activation process mainly by offering extra electrons and stabilizing the transition states. A correlation between the energy barrier and the negative charge of the reaction center is found. The activation barrier is substantially reduced by the aqueous environment, in which the first solvation shell plays the most important role in the barrier reduction. Our approach may be useful for estimating the reaction barriers in aqueous systems.

An Energy Descriptor To Quantify Methane Selectivity in Fischer-Tropsch Synthesis:A Density Functional Theory Study

Selectivity is a fundamental issue in heterogeneous catalysis. In this study, the CH(4) selectivity in Fischer-Tropsch synthesis is chosen to be investigated: CH4 selectivity on Rh, Co, Ru, Fe, and Re surfaces is computed by first-principles methods. In conjunction with kinetic analyses, we are able to derive the effective barrier difference between methane formation and chain growth (Delta E(eff)) to quantify the CH(4) selectivity. By using this energy descriptor, the ranking of methane selectivity predicted from density functional theory (DFT) calculations is consistent with experimental work. Moreover, a linear correlation between Delta E(eff) and the chemisorption energy of C + 4H (Delta H) is found. This fundamental finding possesses the following significance: (i) it shows that the selectivity, which appears to have kinetic characteristics, is largely determined by thermodynamic properties; and (ii) it suggests that an increase of the binding strength of C + 4H will suppress methane selectivity.

Dielectrics in a time-dependent electric field: a real-time approach based on density-polarization functional theory

In the presence of a (time-dependent) macroscopic electric field the electron dynamics of dielectrics cannot be described by the time-dependent density only. We present a real-time formalism that has the density and the macroscopic polarization P as key quantities. We show that a simple local function of P already captures long-range correlation in linear and nonlinear optical response functions. Specifically, after detailing the numerical implementation, we examine the optical absorption, the second- and third-harmonic generation of bulk Si, GaAs, AlAs and CdTe at different level of approximation. We highlight links with ultranonlocal exchange-correlation functional approximations proposed within linear response time-dependent density functional theory framework.