To check the validity of the recently proposed kinetic energy (KE) functional of the pair density (PD) functional theory, we perform the numerical calculation for the neutral Ne atom. This is the first computational trial for the PD functional theory which includes the KE functional that is approximately derived on the basis of the rigorous expression with the coupling-constant integration. In this work we adopt the correlated Thomas-Fermi functional as the approximate KE functional. Although it is shown that non-negligible errors exist in the resulting PD, the present results will be useful as reference data in cases of developing further the approximate KE functional.
We present a computational pair density (PD) functional theory using multiple Slater determinants in constructing the PD. Compared to the effective initial scheme [M. Higuchi and K. Higuchi, Phys. Rev. B 78, 125101 (2008)], the search region for the ground-state PD is substantially extended within the set of
N-representable PDs. The merit of the present scheme is to guarantee the resultant PD to be
N-representable and to describe correlation effects beyond the effective initial scheme, which enables us to tackle with the development of the approximate form of the kinetic energy functional. We derive two restrictive conditions on the kinetic energy functional. With the aid of these restrictive conditions, we also propose several types of approximate functionals of the kinetic energy. To check the abilities and limitations of the present scheme including such approximate functionals, actual calculations are performed for the neutral neon atom. The results show that the present scheme provides the
N-representable and correlated PDs, though it would have some limitations. What we need to do first for improvement of the computational PD functional theory is that the development of the approximate functional of the kinetic energy is further advanced along the present scheme.
Several explicit formulas for the kinetic energy of a many-electron system as a functional of the k-electron density are derived, with emphasis on the electron pair density. The emphasis is on general techniques for deriving approximate kinetic energy functionals and features generalized Weisacker bounds and methods using density-matrix reconstruction. Adapting results from statistical mechanics, a hierarchy of equations is derived that links electron pairs, triplets, quadruplets, etc.; this may be used to derive more accurate approximations. Several methods for defining the exact kinetic energy functional are presented, including the generalizations of the Levy and Lieb formulations of density-functional theory. Together with N-representability constraints on the k-density, this paper provides the basis for “generalized density functional theories” based on the electron pair density. There are also implications for conventional density-functional theory, notably regarding the development of more accurate density functionals for the kinetic energy.
In pair density-functional theory, the only unknown piece of the energy is the kinetic energy
T as a functional of the pair density
P(x1,x2). Although
T[P] has a simpler structure than the Hohenberg-Kohn functional of conventional density-functional theory, computational requirements are still moderate. In the present work, a particularly convenient model system to represent many-electron pair densities is introduced. This ``boson pair model'' (BPM) approximately treats electron pairs as noninteracting bosons. The resulting explicit model for the kinetic energy
T2[P] is shown to be exact for two-electron systems and a lower bound to
T[P] for more than two electrons. The one- and two-particle density matrices obtained from the BPM yield upper bounds for the corresponding many-electron quantities. This suggests a partitioning
T[P]=T2[P]+Teff[P], where only the remainder ${T}_{\mathrm{eff}}\phantom{\rule{0.1em}{0ex}}[P]\ensuremath{\geqslant}0$ needs to be approximated. If the BPM is constrained to yield the exact ground-state pair density, a two-electron Schr\"odinger equation with an effective local two-particle potential results; the latter is identified as a sum of the bare Coulomb interaction and the functional derivative of
Teff[P]. This self-consistent scheme to minimize the energy with respect to
P is more efficient than previous procedures. Further information on the functional derivative of
Teff[P] is derived from a contracted Schr\"odinger equation. Since
Teff[P] is explicitly known in the two-electron and noninteracting (Hartree-Fock) limits, the present method provides an alternative to density-matrix functional theories, which can be exact in the same limits and are similar in computational cost.
A perturbation theory is developed for the correlation energy
Ec[n], of a finite-density system, with respect to the coupling constant \ensuremath{\alpha} which multiplies the electron-electron repulsion operator in ${\mathit{H}}^{\mathrm{\ensuremath{\alpha}}}$=T^+\ensuremath{\alpha}V${\mathrm{^}}_{\mathit{e}\mathit{e}}$+
tsumi${\mathit{v}}_{\mathrm{\ensuremath{\alpha}}}$(
ri). The external potential ${\mathit{v}}_{\mathrm{\ensuremath{\alpha}}}$ is constrained to keep the gound-state density n fixed for all \ensuremath{\alpha}\ensuremath{\ge}0. ${\mathit{v}}_{\mathrm{\ensuremath{\alpha}}}$ is given completely in terms of functional derivatives at full charge (\ensuremath{\alpha}=1), from which
Ec[${\mathit{n}}_{\ensuremath{\lambda}}$]=
ec,2[n]+ ${\ensuremath{\lambda}}^{\mathrm{\ensuremath{-}}1}$
ec,3[n]+${\ensuremath{\lambda}}^{\mathrm{\ensuremath{-}}2}$
ec,4[n]+..., where each
ec,j[n] is expressed in terms of integrals involving Kohn-Sham determinants. Here, ${\mathit{n}}_{\ensuremath{\lambda}}$(x,y,x)=${\ensuremath{\lambda}}^{3}$n(\ensuremath{\lambda}x,\ensuremath{\lambda}y,\ensuremath{\lambda}z) and \ensuremath{\lambda}=${\mathrm{\ensuremath{\alpha}}}^{\mathrm{\ensuremath{-}}1}$. The identification of ${\mathrm{lim}}_{\ensuremath{\lambda}\ensuremath{\rightarrow}\mathrm{\ensuremath{\infty}}}$
Ec[${\mathit{n}}_{\ensuremath{\lambda}}$], which is a high-density limit, as the second-order energy
ec,2[n] allows one to compute bounds upon ${\mathrm{lim}}_{\ensuremath{\lambda}\ensuremath{\rightarrow}\mathrm{\ensuremath{\infty}}}$
Ec[${\mathit{n}}_{\ensuremath{\lambda}}$]; the bounds imply that ${\mathrm{lim}}_{\ensuremath{\lambda}\ensuremath{\rightarrow}\mathrm{\ensuremath{\infty}}}$
Ec[${\mathit{n}}_{\ensuremath{\lambda}}$]\ensuremath{\simeq}
Ec[n] for a large class of small atoms and molecules, and suggest that ${\mathrm{lim}}_{\ensuremath{\lambda}\ensuremath{\rightarrow}\mathrm{\ensuremath{\infty}}}$
Ec[${\mathit{n}}_{\ensuremath{\lambda}}$] should be of the same order of magnitude as
Ec[n] in finite insulators and semiconductors.Approximations to
Ec[n] should reflect all this. In contrast, perhaps the well-known overbinding of the local-density approximation (LDA) in molecules and solids is due, in part, to the fact that the LDA correlation energy is too sensitive to a coordinate scaling of n. Indeed, the LDA for
Ec[${\mathit{n}}_{\ensuremath{\lambda}}$] diverges when \ensuremath{\lambda}\ensuremath{\rightarrow}\ensuremath{\infty} because of the presence of the -ln(\ensuremath{\lambda}) term in the Gell-Mann and Brueckner high-density expression for the correlation energy, per particle, of a homogeneous density, which is infinite. In a sense, the derived perturbation expansion transforms the Gell-Mann and Brueckner expression into one that applies specifically to an inhomogeneous density which integrates to a finite number of electrons.
Abstract Exact equations for correlation energy and kinetic energy correction functionals are derived in the time‐dependent density functional theory. Hierarchies of equations are also obtained and approximate functionals are proposed by truncation of the hierarchies at different orders. Explicit relations for the correlation energy and its kinetic component along the adiabatic path are given in terms of powers of the moments of the electron density. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003