Total energy evaluation in the Strutinsky shell correction method

We analyze the total energy evaluation in the Strutinsky shell correction method (SCM) of Ullmo et al. [Phys. Rev. B 63, 125339 (2001)], where a series expansion of the total energy is developed based on perturbation theory. In agreement with Yannouleas and Landman [Phys. Rev. B 48, 8376 (1993)], we also identify the first-order SCM result to be the Harris functional [Phys. Rev. B 31, 1770 (1985)]. Further, we find that the second-order correction of the SCM turns out to be the second-order error of the Harris functional, which involves the a priori unknown exact Kohn-Sham (KS) density, rho(KS)(r). Interestingly, the approximation of rho(KS)(r) by rho(out)(r), the output density of the SCM calculation, in the evaluation of the second-order correction leads to the Hohenberg-Kohn-Sham functional. By invoking an auxiliary system in the framework of orbital-free density functional theory, Ullmo et al. designed a scheme to approximate rho(KS)(r), but with several drawbacks. An alternative is designed to utilize the optimal density from a high-quality density mixing method to approximate rho(KS)(r). Our new scheme allows more accurate and complex kinetic energy density functionals and nonlocal pseudopotentials to be employed in the SCM. The efficiency of our new scheme is demonstrated in atomistic calculations on the cubic diamond Si and face-centered-cubic Ag systems.     

Comparing ab initio density-functional and wave function theories: the impact of correlation on the electronic density and the role of the correlation potential

The framework of ab initio density-functional theory (DFT) has been introduced as a way to provide a seamless connection between the Kohn-Sham (KS) formulation of DFT and wave-function based ab initio approaches [R. J. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005)]. Recently, an analysis of the impact of dynamical correlation effects on the density of the neon atom was presented [K. Jankowski, K. Nowakowski, I. Grabowski, and J. Wasilewski, J. Chem. Phys. 130, 164102 (2009)], contrasting the behaviour for a variety of standard density functionals with that of ab initio approaches based on second-order M?ller-Plesset (MP2) and coupled cluster theories at the singles-doubles (CCSD) and singles-doubles perturbative triples [CCSD(T)] levels. In the present work, we consider ab initio density functionals based on second-order many-body perturbation theory and coupled cluster perturbation theory in a similar manner, for a range of small atomic and molecular systems. For comparison, we also consider results obtained from MP2, CCSD, and CCSD(T) calculations. In addition to this density based analysis, we determine the KS correlation potentials corresponding to these densities and compare them with those obtained for a range of ab initio density functionals via the optimized effective potential method. The correlation energies, densities, and potentials calculated using ab initio DFT display a similar systematic behaviour to those derived from electronic densities calculated using ab initio wave function theories. In contrast, typical explicit density functionals for the correlation energy, such as VWN5 and LYP, do not show behaviour consistent with this picture of dynamical correlation, although they may provide some degree of correction for already erroneous explicitly density-dependent exchange-only functionals. The results presented here using orbital dependent ab initio density functionals show that they provide a treatment of exchange and correlation contributions within the KS framework that is more consistent with traditional ab initio wave function based methods.     

Improving density functional theory for crystal polymorph energetics

We show that the quality of density functional theory (DFT) predictions for the relative stabilities of polymorphs of crystalline para-diiodobenzene (PDIB) is dramatically improved through a simple two-body correction using wavefunction-based electronic structure theory. PDIB has two stable polymorphs under ambient conditions, and like Hongo et al. [J. Phys. Chem. Lett., 1, 1789 (2010)] we find that DFT makes wildly variable predictions of the relative stabilities, depending on the approximate functional used. The two-body corrected scheme, using Grimme's spin-scaled variant of second-order M?ller-Plesset perturbation theory and any of the tested density functionals, predicts the α-polymorph to be more stable, consistent with experiment, and produces a relative stability that agrees with the benchmark quantum Monte-Carlo results of Hongo et al. within statistical uncertainty.     

Semiempirical double-hybrid density functional with improved description of long-range correlation

The recently proposed new family of "double-hybrid" density functionals [Grimme, S. J. Chem. Phys. 2006, 124, 34108] replaces a fraction of the semi-local correlation energy by a non-local correlation energy expression that employs the Kohn-Sham orbitals in second-order many-body perturbation theory. These functionals have provided results of high accuracy over a wide range of properties but fail to accurately describe long-range van der Waals interactions. In this work, a distance-dependent scaling factor for the non-local correlation energy is introduced to address this problem, and two new double-hybrid density functionals are proposed. The new functionals are optimized with the finite cc-pVTZ basis on training sets of atomization energies and intermolecular interaction energies. They are compared against (scaled) second-order M?ller-Plesset perturbation theories and popular density functionals including the hybrid-GGA functional B3-LYP and the first double-hybrid functional (B2-PLYP). Tests are performed on an extensive set including reaction energies, barrier heights, weakly interacting complexes, transition-metal systems, molecular geometries, and harmonic vibrational frequencies. Within the cc-pVTZ atomic orbital basis, we have demonstrated the ability to find a parametrization scheme which is simultaneously able to describe thermochemistry and weakly bound systems with a satisfactory degree of accuracy.     

XYG3型双杂化密度泛函方法新进展：从能量到能量的解析梯度

本文回顾了现代密度泛函理论的基础,着重评述了XYG3型双杂化（XYG3typeofdoublyhybrid,xDH）泛函的最新进展,解析能量梯度的实现．XYG3是首个依照绝热途径理论建立的双杂化泛函,在具体实现上具有独特的构架．该类型泛函利用常用泛函（如B3LYP或PBE0等）作母泛函来进行自洽计算,以期获得更好的密度和轨道,然后将所得到的轨道和密度信息带入到xDH泛函中以得到最终能量．由于自洽泛函和最终能量泛函不同,因而在计算解析能量梯度时需要求解耦合微扰Kohn—Sham方程．在此基础上,还评述了xDH泛函在能量,尤其是构型优化方面的具体表现．测试的构型集包括以共价键键合的分子和非键相互作用体系的平衡结构,以及反应过渡态结构．结果表明,xDH双杂化泛函总体上给出了比母泛函更好的能量和几何构型．     

Koopmans-like approximation in the Kohn-Sham method and the impact of the frozen core approximation on the computation of the reactivity parameters of the density functional theory

A Koopmans-like approximation is introduced in the spin-polarized version of the Kohn-Sham (KS) density functional theory to obtain a relation between KS orbital energies and vertical ionization potential and electron affinity. Expressions for reactivity indexes (like electronegativity, hardness, electrophilicity, and excitation energies) include KS frontier orbital energies and additional contributions associated with the self-interaction correction. Those reactivity parameters were computed with different exchange-correlation functionals to test the approach for a set of small molecules. The results show that the present approximation provides a better way to estimate hardness, electronegativity, and electrophilicity than just the use of frontier orbital energy values. However KS HOMO and LUMO energy gap gives a better agreement with excitation energies.     

Orbital-dependent correlation energy in density-functional theory based on a second-order perturbation approach: success and failure

We have developed a second-order perturbation theory (PT) energy functional within density-functional theory (DFT). Based on PT with the Kohn-Sham (KS) determinant as a reference, this new ab initio exchange-correlation functional includes an exact exchange (EXX) energy in the first order and a correlation energy including all single and double excitations from the KS reference in the second order. The explicit dependence of the exchange and correlation energy on the KS orbitals in the functional fits well into our direct minimization approach for the optimized effective potential, which is a very efficient method to perform fully self-consistent calculations for any orbital-dependent functionals. To investigate the quality of the correlation functional, we have applied the method to selected atoms and molecules. For two-electron atoms and small molecules described with small basis sets, this new method provides excellent results, improving both second-order Moller-Plesset expression and any conventional DFT results significantly. For larger systems, however, it performs poorly, converging to very low unphysical total energies. The failure of PT based energy functionals is analyzed, and its origin is traced back to near degeneracy problems due to the orbital- and eigenvalue-dependent algebraic structure of the correlation functional. The failure emerges in the self-consistent approach but not in perturbative post-EXX calculations, emphasizing the crucial importance of self-consistency in testing new orbital-dependent energy functionals.     

Density functional theory of complex transition densities

We present an extension of Hohenberg-Kohn-Sham density functional theory to the domain of complex local potentials and complex electron densities. The approach is applicable to resonance (Siegert) [Phys. Rev. 56, 750 (1939)] states and other scattering and transport problems that can be described by a normalized state of a Hamiltonian containing a complex local potential. Such Hamiltonians are non-Hermitian and their eigenvalues are in general complex, the imaginary part being inversely proportional to the lifetime of the system. The one-to-one correspondence between complex local potentials nu and complex electron densities rho is established provided that the complex variables are sufficiently close to real local potentials and densities of nondegenerate ground states. We show that the exchange-correlation functionals, contributing to the complex energy, are determined through analytic continuation of their ground-state-theory counterparts. This implies that the exchange-correlation effects on the lifetime of a resonance are, under appropriate conditions, already determined by the functionals of the ground-state theory.     

Exact Kohn-Sham versus Hartree-Fock in momentum space: Examples of two-fermion systems

The question of how density functional theory (DFT) compares with Hartree-Fock (HF) for the computation of momentum-space properties is addressed in relation to systems for which (near) exact Kohn-Sham (KS) and HF one-electron matrices are known. This makes it possible to objectively compare HF and exact KS and hence to assess the potential of DFT for momentum-space studies. The systems considered are the Moshinsky [Am. J. Phys. 36, 52 (1968)] atom, Hooke's atom, and light two-electron ions, for which expressions for correlated density matrices or momentum densities have been derived in closed form. The results obtained show that it is necessary to make a distinction between true and approximate DFTs.     

Density functional study of double ionization energies

In this paper, double ionization energies (DIEs) of gas-phase atoms and molecules are calculated by energy difference method with density functional theory. To determine the best functional for double ionization energies, we first study 24 main group atoms in the second, third, and fourth periods. An approximation is used in which the electron density is first obtained from a density functional computation with the exchange-correlation potential V xc known as statistical average of orbital potentials, after which the energy is computed from that density with 59 different exchange-correlation energy functionals E xc. For the 24 atoms, the two best E xc functional providing DIEs with average absolute deviation (AAD) of only 0.25 eV are the Perdew-Burke-Ernzerhof functional modified by Hammer et al. [Phys. Rev. B 59, 6413 (1999)] and one known as the Krieger-Chen-Iafrate-Savin functional modified by Krieger et al. (unpublished). Surprisingly, none of the 20 available hybrid functionals is among the top 15 functionals for the DIEs of the 24 atoms. A similar procedure is then applied to molecules, with opposite results: Only hybrid functionals are among the top 15 functionals for a selection of 29 molecules. The best E xc functional for the 29 molecules is found to be the Becke 1997 functional modified by Wilson et al. [J. Chem. Phys. 115, 9233 (2001)]. With that functional, the AAD from experiment for DIEs of 29 molecules is just under 0.5 eV. If the two suspected values for C2H2 and Fe(CO)5 are excluded, the AAD improves to 0.3(2) eV. Many other hybrid functionals perform almost as well.     

Assessment of the efficiency of long-range corrected functionals for some properties of large compounds

Using the long-range correction (LC) density functional theory (DFT) scheme introduced by Iikura et al. [J. Chem. Phys. 115, 3540 (2001)] and the Coulomb-attenuating model (CAM-B3LYP) of Yanai et al. [Chem. Phys. Lett. 393, 51 (2004)], we have calculated a series of properties that are known to be poorly reproduced by standard functionals: Bond length alternation of pi-conjugated polymers, polarizabilities of delocalized chains, and electronic spectra of extended dyes. For each of these properties, we present cases in which traditional hybrid functionals do provide accurate results and cases in which they fail to reproduce the correct trends. The quality of the results is assessed with regard to experimental values and/or data arising from electron-correlated wave function approaches. It turns out that (i) both LC-DFT and CAM-B3LYP provide an accurate bond length alternation for polyacetylene and polymethineimine, although for the latter they decrease slightly too rapidly with chain length. (ii) The LC generalized gradient approximation and MP2 polarizabilities of long polyphosphazene and polymethineimine oligomers agree almost perfectly. In the same way, CAM-B3LYP corrects the major part of the B3LYP faults. (iii) LC and CAM techniques do not help in correcting the nonrealistic evolution with chain length of the absorption wavelengths of cyanine derivatives. In addition, though both schemes significantly overestimate the ground to excited state transition energy of substituted anthraquinone dyes, they provide a more consistent picture once a statistical treatment is performed than do traditional hybrid functionals.     

Second-order Møller-Plesset perturbation energy obtained from divide-and-conquer Hartree-Fock density matrix

The density matrix (DM) obtained from Yang's [Phys. Rev. Lett. 66, 1438 (1991)] divide-and-conquer (DC) Hartree-Fock (HF) calculation is applied to the explicit second-order M?ller-Plesset perturbation (MP2) energy functional of the HF DM, which was firstly mentioned by Ayala and Scuseria [J. Chem. Phys. 110, 3660 (1999)] and was improved by Surján [Chem. Phys. Lett. 406, 318 (2005)] as DM-Laplace MP2. This procedure, termed DC-DM MP2, requires the HF DM of holes, for which we propose two evaluation schemes in DC manner. Numerical studies reveal that the DC-DM MP2 energy deviation from canonical MP2 is the same order of magnitude as DC-HF energy deviation from conventional HF whichever type of hole DM is adopted. It is also confirmed that the central processing unit time of DC-DM MP2 is less than that of DM-Laplace MP2 because the DC-HF DM is sparser than conventional DM.     

Regional self-interaction correction of density functional theory

We propose a new simple scheme for self-interaction correction (SIC) of exchange functionals in the density functional theory. In the new scheme, exchange energies are corrected by substituting exchange self-interactions for exchange functionals in regions of self-interaction. To classify the regions of self-interaction, we take advantage of the property of the total kinetic energy density approaching the Weizs?cker density in the case of electrons in isolated orbitals. The scheme differs from conventional SIC methods in that it produces optimized molecular structures. Applying the scheme to the calculation of reaction energy barriers showed that it provides a clear improvement in cases where the barriers are underestimated by conventional "pure" functionals. In particular, we found that this scheme even reproduces a transition state that is not given by pure functionals.     

Performance of the density matrix functional theory in the quantum theory of atoms in molecules

The generalization to arbitrary molecular geometries of the energetic partitioning provided by the atomic virial theorem of the quantum theory of atoms in molecules (QTAIM) leads to an exact and chemically intuitive energy partitioning scheme, the interacting quantum atoms (IQA) approach, that depends on the availability of second-order reduced density matrices (2-RDMs). This work explores the performance of this approach in particular and of the QTAIM in general with approximate 2-RDMs obtained from the density matrix functional theory (DMFT), which rests on the natural expansion (natural orbitals and their corresponding occupation numbers) of the first-order reduced density matrix (1-RDM). A number of these functionals have been implemented in the promolden code and used to perform QTAIM and IQA analyses on several representative molecules and model chemical reactions. Total energies, covalent intra- and interbasin exchange-correlation interactions, as well as localization and delocalization indices have been determined with these functionals from 1-RDMs obtained at different levels of theory. Results are compared to the values computed from the exact 2-RDMs, whenever possible.     

Translational and rotational invariance requisites for density functional derivatives

Translational and rotational invariance of functionals lead to hierarchies of equations between successive derivatives. These hierarchies allow alternating series expansions of some density functionals in terms of functional derivatives and charge density. Translational and rotational invariance also give rise to integrodifferential equations that link derivatives of all orders. From the minimal properties of the kinetic energy functional T_s[ρ] and the functional F[ρ] = min_Ψ→ρ <Ψ|T + V_ee|Ψ>, it follows that $int dˆ3 r dˆ3 rˆprimef({bold r}) {delˆ2 T―s[rho]}over{delrho({bold r})delrho({bold rˆprime})}f({bold rˆprime}) geq 0hspace{1cm}hbox{and}hspace{1cm}int dˆ3 r dˆ3 rˆprimef({bold r}) {delˆ2 F[rho]}over{delrho({bold r})delrho({bold rˆprime})}f({bold rˆprime}) geq 0$ for all ∫ d³ r d³ f′ f(r) = 0. This property combined with constraints on functionals due to translational invariance lead to inequalities satisfied by first derivatives of selected density functionals. © 1997 John Wiley & Sons, Inc.     

Universal mathematical identities in density functional theory: results from three different spin-resolved representations

This paper supersedes previous theoretical approaches to conceptual DFT because it provides a unified and systematic approach to all of the commonly considered formulations of conceptual DFT, and even provides the essential mathematical framework for new formulations. Global, local, and nonlocal chemical reactivity indicators associated with the "closed-system representation" ([N(alpha),N(beta),nu(alpha)(r),nu(beta)(r)]) of spin-polarized density functional theory (SP-DFT) are derived. The links between these indicators and the ones associated with the "open-system representation" ([mu(alpha),mu(beta),nu(alpha)(r),nu(beta)(r)]) are derived, including the spin-resolved Berkowitz-Parr identity. The Legendre transform to the "density representation" ([rho(alpha)(r),rho(beta)(r)]) is performed, and the spin-resolved Harbola-Chattaraj-Cedillo-Parr identities linking the density representation to the closed-system and open-system representations are derived. Taken together, these results provide the framework for understanding chemical reactions from both the electron-following perspective (using either the closed-system or the open-system representation) and electron-preceding perspective (density representation). A powerful matrix-vector notation is developed; with this notation, identities in conceptual DFT become universal. Specifically, this notation allows the fundamental identities in conventional (spin-free) conceptual DFT, the [N(alpha),N(beta)] representation, and the [N=N(alpha)+N(beta),N(S)=N(alpha)-N(beta)] representation to be written in exactly the same forms. In cases where spin transfer and electron transfer are coupled (e.g., radical+molecule reactions), we believe that the [N(alpha),N(beta)] representation may be more useful than the more common [N,N(S)] representation.     

Time-dependent density-functional studies of the D2 coulomb explosion

Real-time first principle simulations are presented of the D(2) Coulomb explosion dynamics detonated by exposure to very intense few-cycle laser pulse. Three approximate functionals within the time-dependent density functional theory (TDDFT) functionals are examined for describing the electron dynamics, including time-dependent Hartree-Fock theory. Nuclei are treated classically with quantum corrections. The calculated results are sensitive to the underlying electronic structure theory, showing too narrow kinetic energy distribution peaked at too high kinetic energy when compared with recent experimental results (Phys. Rev. Lett. 2003, 91, 093002). Experiment also shows a low energy peak which is not seen in the present calculation. We conclude that while Ehrenfest-adiabatic-TDDFT can qualitatively account for the dynamics, it requires further development, probably beyond the adiabatic approximation, to be quantitative.