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.
We derive the coupling-constant expression for the kinetic energy (KE) functional of the pair-density-functional theory. Such the expression can be effectively used in developing the approximate KE functional. Indeed, on the basis of this expression, we present an approximate form of the KE functional by means of the perturbation theory. Furthermore, as the first step of an alternative way of obtaining the approximate KE functional, we derive several relations and bounds for the KE functional by investigating how the KE functional changes with the nonuniform scaling, rotation, and translation of electron coordinates. Such relations and bounds can be utilized as restrictive conditions on the approximate form. Thus, the coupling-constant expression, relations and bounds obtained here provide useful tools for devising various types of the approximate KE functional that is indispensable for advancing the pair-density-functional theory.
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
In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
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.
A new approach for the calculation of ground states of many-electron systems is proposed via an integral formulation of the Hohenberg-Kohn-Sham density-functional theory. Only equations for the total electron density are involved; orbitals are not employed. Exchange and correlation effects are incorporated. In place of the set of single-electron equations, the total electron density is explicitly expressed in terms of the Kohn-Sham effective local potential through multidimensional integrations. The development is based on the first-order density matrix as obtained from the one-body Green's function in polygonal and Fourier path-integral representations. This might open up the possibility of ab initio calculations for molecules with very many electrons. It also provides explicit solutions to two long-standing problems: electron kinetic energy and momentum density as functionals of the total electron density. The formulation can also be used in calculations for general fermions.