Solving the ultranonlocality problem in time-dependent spin-density-functional theory

It has been known for some time that the exchange-correlation potential in time-dependent density-functional theory is an intrinsically nonlocal functional of the density as soon as one goes beyond the adiabatic approximation. In this paper we show that a much more severe nonlocality problem, with a completely different physical origin, plagues the exchange-correlation potentials in time-dependent spin-density functional theory. We show how the use of the spin current density as a basic variable solves this problem, and we provide an explicit local expression for the exchange-correlation fields as functionals of the spin currents.     

The exchange-correlation energy of a metallic surface

The exchange-correlation energy of a metal surface is analyzed in terms of the wavelength of the fluctuations which contribute to it. A scheme is proposed to interpolate between the short-wavelength region properly described by the local density functional approximation and the long-wavelength region for which an exact limiting form is established. This interpolation scheme is tested and verified for the soluble infinite barrier model and then applied to realistic density profiles.     

Post-Hartree-Fock methods and dynamic correlation in atoms and molecules

The methods of calculation of the matrix of the exchange-correlation interaction are considered within the framework of one post-Hartree-Fock one-electron method of investigation of the properties of many-electron systems. Such post-Hartree-Fock methods are based on two-step variational self-consistent calculations of the spin orbitals and superposition coefficients of configurations in the multiconfiguration approximation. The post-Hartree-Fock method used involves an approach related to the extended Koopmans’ theorem, which, in turn, proves to be a high-energy approximation for quantum Green’s functions. Obvious application areas of the calculations of the exchange-correlation interaction within the framework of the method proposed are the multiparticle perturbation theory, the parameterization of the energy representation as a functional of the single-particle density matrix, and the theory of Green’s functions in the multiconfiguration approximation. A relativistic generalization of the method with the aim of calculating the radiative corrections for many-electron atoms and for problems of interaction with an external field in the nonstationary Floquet theory is possible.     

Inclusion of the exchange-correlation effects in Ab initio methods for calculating the plasmon dispersion and line width in metals

A self-consistent analog of the Hubbard local-field factor accounting for a real band structure of the metal under investigation is proposed for the calculation of the exchange-correlation corrections to the randomphase approximation. The plasmon energy and the plasmon line width for potassium are calculated as a function of the wave vector. The calculations are performed within the random-phase approximation and with corrections for the local field. It is shown that the results of the calculations are in good agreement with the experimental data. This indicates that, first, the corrections to the random-phase approximation should be taken into account for metals with a low electron density and, second, the local-field factor approximation is an appropriate approach.     

Optimization of effective atom centered potentials for london dispersion forces in density functional theory

We add an effective atom-centered nonlocal term to the exchange-correlation potential in order to cure the lack of London dispersion forces in standard density functional theory. Calibration of this long-range correction is performed using density functional perturbation theory and an arbitrary reference. Without any prior assignment of types and structures of molecular fragments, our corrected generalized gradient approximation density functional theory calculations yield correct equilibrium geometries and dissociation energies of argon-argon, benzene-benzene, graphite-graphite, and argon-benzene complexes.     

Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids

The electron density, its gradient, and the Kohn-Sham orbital kinetic energy density are the local ingredients of a meta-generalized gradient approximation (meta-GGA). We construct a meta-GGA density functional for the exchange-correlation energy that satisfies exact constraints without empirical parameters. The exchange and correlation terms respect two paradigms: one- or two-electron densities and slowly varying densities, and so describe both molecules and solids with high accuracy, as shown by extensive numerical tests. This functional completes the third rung of "Jacob's ladder" of approximations, above the local spin density and GGA rungs.     

Non-local exchange correlation energy in quasi-two-dimensional electron layers

We introduce a simple, non-local approximation for the exchange-correlation energy of an inhomogeneous electron gas. This approximation is shown to be quantitatively accurate in several important limiting cases. The method is used to calculate the non-local exchange-correlation energy for density distributions representing quasi-two-dimensional electron layers. These energies are compared to the widely used three-dimensional local-density approximation and suggest that the latter approximation may contain errors on the order of 10–30% for layer profiles typical of real systems.     

Nonuniform scaling applied to surface energies of transition metals

We construct a generalized gradient approximation of the exchange-correlation energy that satisfies the nonuniform scaling in one dimension and is accurate in the whole quasi-two-dimensional (Q2D) regime. Using spatial and energetic analyses of metal (111) surfaces, we show that the Q2D behavior is important at the surface of most transition metals, and that the here proposed Q2D-generalized gradient approximation functional predicts for these metals accurate surface energies as well as bulk properties.     

Non-Local Density Functional Description of Poly-Para-Phenylene Vinylene

A fully non-local exchange-correlation formalism the weighted density approximation （WDA）, has within the framework of density functional theory, known as been applied to the conjugated polymer poly-para-phenylene vinylene （PPV） and is shown to lead to a marked improvement in the agreement of theory and experiment for the electronic band structure of the conjugated polymer. In particular, some new model WDA functions are developed, which substantially increase the electronic band gap of the polymer relative to those obtained with the local density approximation and generalized gradient approximation. The calculated band gap of PPV is quantitatively or atleast semiquantitatively in agreement with the experimental data.     

Density functional theory of transition metal phthalocyanines,I: electronic structure of NiPc and CoPc—self-interaction effects

We present a two-part systematic density functional theory study of the electronic structure of selected transition metal phthalocyanines. We use a semi-local generalized gradient approximation (GGA) functional, as well as several hybrid exchange-correlation functionals, and compare the results to experimental photoemission data. Here, we study the low-spin systems NiPc and CoPc. We show that hybrid functionals provide computed photoemission spectra in excellent agreement with experimental data, whereas the GGA functional fails qualitatively. This failure is primarily because of under-binding of localized orbitals due to self-interaction errors.     

Meta-generalized gradient approximation: non-empirical construction and performance of a density functional

The local ingredients of a meta-generalized gradient approximation (meta-GGA) include the electron density, its gradient, and the Kohn–Sham orbital kinetic energy density. We discuss the strategy of constructing a successful meta-GGA density functional for the exchange-correlation energy that satisfies exact constraints without empirical parameters. The new feature of this functional is that it simultaneously respects the two paradigms of electronic structure theory: one- or two-electron densities and slowly-varying densities, and so is uniformly accurate for atoms, molecules and solids. Results of extensive numerical tests of the new functional are summarized and evaluated.     

Multiresolution quantum chemistry in multiwavelet bases: time-dependent density functional theory with asymptotically corrected potentials in local density and generalized gradient approximations

A multiresolution solver for fully numerical linear response calculations of excitation states via the time-dependent Hartree–Fock and density functional theory (TD-HF/DFT) is presented. The linear response method Yanai et al. previously reported [J. Chem. Phys., submitted] was limited to the Tamm–Dancoff approximation and could only use the Hartree–Fock exchange and the local-spin density approximation (LSDA) with a crude asymptotic correction. The present development enables us to perform full TD-HF/DFT calculations employing generalized gradient approximation (GGA) exchange-correlation potentials as well as hybrid ones. The linear response of TD-HF/DFT is computed by means of iteratively solving the coupled integral equations with the Green's functions. In this study, Tozer and Handy's asymptotic correction (AC) is applied to existing DFT exchange-correlations, and is found numerically stable and efficient. Furthermore, the new hybrid exchange-correlation functional CAM-B3LYP, which was recently proposed by Yanai et al. [Chem. Phys. Lett. 393, 51 (2004)], is implemented. The implementation requires a new separated representation of the integral kernel for the Coulomb-attenuated potential. We demonstrate linear response calculations free of basis set error for the excited states of Be, N₂, C₂H₄ and C₆H₆ using LSDA, HCTH, CAM-B3LYP and PBE0 exchange-correlation functionals. The mean absolute errors of the C₆H₆ calculations with HCTH and CAM-B3LYP are 0.12 and 0.18?eV, respectively. The second derivative of exchange-correlation functionals is represented fully numerically at O(N) computation cost.     

First principles calculations of structural,electronic and optical properties of various phases of CaS

First principles calculations have been performed within the framework of density functional theory to investigate the structural, electronic and optical properties of all four possible B1, B2, B3 and B4 phases of CaS. Apart from the standard local density approximation (LDA) and GGA (PBE), a more accurate nonempirical density functional generalized gradient approximation (GGA), as proposed by Wu and Cohen [Phys. Rev. B 73, 235116 (2006)] for the exchange-correlation energy, E_XC, has been attempted in these calculations. Calculated electronic structure and the density of states are analyzed in terms of the contribution of Ca d states and S s and p states in determining the nature of the fundamental band gap in various phases. Reflectivity, R (ω), the real and imaginary part of the dielectric functions, ε(ω), have been calculated for all the phases and the results have been discussed and compared with the existing experimental data.     

Dynamical corrections to the DFT-LDA electron conductance in nanoscale systems

Using time-dependent current-density functional theory, we derive analytically the dynamical exchange-correlation correction to the dc conductance of nanoscale junctions. The correction pertains to the conductance calculated in the zero-frequency limit of time-dependent density functional theory within the adiabatic local-density approximation. In particular, we show that in linear response, the correction depends nonlinearly on the gradient of the electron density; thus, it is more pronounced for molecular junctions than for quantum point contacts. We provide specific numerical examples to illustrate these findings.     

高压下氧化镉弹性性质、电子结构和光学性质的第一性原理研究

 运用基于密度泛函理论的平面波赝势方法（PWP）,计算研究了氧化镉NaCl结构（B1结构）和CsCl结构（B2结构）在不同压力条件下的几何结构、弹性性质、电子结构和光学性质。交换关联能分别采用广义梯度近似（GGA）和局域密度近似（LDA）。通过比较计算和实验得到的晶格常数和体模量不难发现,LDA的计算结果更符合实验值。在高压的作用下,两种结构的导带能级有向高能级移动的趋势,而价带能级有向低能级移动的趋势,因此直接带隙变大。同时,对照态密度分布图及高压下能级的移动情况,分析了CdO两种结构在高压作用下的光学性质。     

Density functional theory of transition metal phthalocyanines,II: electronic structure of MnPc and FePc—symmetry and symmetry breaking

We present a two-part systematic density functional theory (DFT) study of the electronic structure of selected transition metal phthalocyanines. We use a semi-local generalized gradient approximation (GGA) functional, as well as several hybrid exchange-correlation functionals, and compare the results to experimental photoemission data. Here, we study the intermediate spin systems MnPc and FePc. We show that DFT calculations of these systems are extremely sensitive to the choice of functional and basis set with respect to the obtained electronic configuration and to symmetry breaking. Interestingly, all simulated spectra are in good agreement with experiment despite the differences in the underlying electronic configurations.     

Calculation of the nonlinear susceptibilities of interacting homogeneous electron gas by the density functional method

The density functional method is used for the calculation of nonlinear susceptibilities of an interacting electron gas. It is shown that this approach leads to the expression for nonlinear susceptibilities which satisfy exactly the strict sum rules derived previously. By using the local densities approximation we obtain the expression for the exchange-correlation contribution in the first non-linear susceptibility of interacting electron gas. It is shown that such a contribution may have a noticeable effect on the binding energy and phonon spectrum of polivalent simple metals.     

Local density functional approximation applied to Compton profiles

We have applied a self-consistent psuedopotential approach in the density functional approximation to study the dependence of the Compton profile of silicon on the choice of the exchange correlation potential. We find that the Compton profile depends strongly on the value of the exchange-correlation potential.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 1, pp. 76–79, January, 1996.     

Localization and delocalization errors in density functional theory and implications for band-gap prediction

The band-gap problem and other systematic failures of approximate exchange-correlation functionals are explained from an analysis of total energy for fractional charges. The deviation from the correct intrinsic linear behavior in finite systems leads to delocalization and localization errors in large and bulk systems. Functionals whose energy is convex for fractional charges such as the local density approximation display an incorrect apparent linearity in the bulk limit, due to the delocalization error. Concave functionals also have an incorrect apparent linearity in the bulk calculation, due to the localization error and imposed symmetry. This resolves an apparent paradox and identifies the physical nature of the error to be addressed to obtain accurate band gaps from density functional theory.