Strong electron correlation remains difficult to describe within conventional density functional approximations, particularly when natural-orbital occupations become strongly fractional. We introduce a systematic particle-hole expansion that generalizes density functional theory (DFT) through an explicit functional of the one-electron reduced density matrix. The previously derived cumulant-based correction emerges as the leading term and the thermal (entropic) correction as an infinite-order member whose coefficients are fixed by a single parameter. The expansion can be formulated as a moment problem solved by semidefinite programming with computational scaling comparable to DFT. We show that both the first- and second-order corrections yield closed-form natural-orbital occupations. Applications to periacenes and Fe-porphyrin demonstrate that the second-order correction substantially reduces excessive orbital fractionalization, yielding natural-orbital occupations and fundamental gaps in improved agreement with correlated calculations.
Strong electron correlation remains difficult to describe within conventional density functional approximations, particularly when natural-orbital occupations become strongly fractional. We introduce a systematic particle-hole expansion that generalizes density functional theory (DFT) through an explicit functional of the one-electron reduced density matrix. The previously derived cumulant-based correction emerges as the leading term and the thermal (entropic) correction as an infinite-order member whose coefficients are fixed by a single parameter. The expansion can be formulated as a moment problem solved by semidefinite programming with computational scaling comparable to DFT. We show that both the first- and second-order corrections yield closed-form natural-orbital occupations. Applications to periacenes and Fe-porphyrin demonstrate that the second-order correction substantially reduces excessive orbital fractionalization, yielding natural-orbital occupations and fundamental gaps in improved agreement with correlated calculations.
Density functional theory (DFT) is the de facto method for the electronic structure of weakly correlated systems. But for strongly correlated materials, common density functional approximations break down. Here, we derive a many-electron expansion (MEE) in DFT that accounts for successive one-, two-, three-, ... particle interactions within the system. To compute the correction terms, the density is first decomposed into a sum of localized, nodeless one-electron densities (${\ensuremath{\rho}}_{i}$). These one-electron densities are used to construct relevant two- (${\ensuremath{\rho}}_{i}+{\ensuremath{\rho}}_{j}$), three- (${\ensuremath{\rho}}_{i}+{\ensuremath{\rho}}_{j}+{\ensuremath{\rho}}_{k}$), ... electron densities. Numerically exact results for these few-particle densities can then be used to correct an approximate density functional via any of several many-body expansions. We show that the resulting hierarchy gives accurate results for several important model systems: the Hubbard and Peierls-Hubbard models in 1D and the pure Hubbard model in 2D. We further show that the method is numerically convergent for strongly correlated systems: applying successively higher order corrections leads to systematic improvement of the results. MEE thus provides a hierarchy of density functional approximations that applies to both weakly and strongly correlated systems.
We introduce a novel energy functional for ground-state electronic-structure calculations. Its fundamental variables are the natural spin-orbitals of the implied singlet many-body wave function and their joint occupation probabilities. The functional derives from a sequence of controlled approximations to the two-particle density matrix. Algebraic scaling of computational cost with electron number is obtainable in general, and Hartree-Fock scaling in the seniority-zero version of the theory. Results obtained with the latter version for saturated small molecular systems are compared with those of highly-accurate quantum-chemical computations. The numerical results are variational, capturing most of the correlation energy from equilibrium to dissociation. Their accuracy is considerably greater than that obtainable with current density-functional theory approximations and with current functionals of the one-particle density matrix only.
Density-Functional Theory (DFT), although widely used and very successful in the calculation of several observables, fails to correctly describe strongly correlated materials. In the first part of this work we, therefore, introduce reduced-density-matrix-functional theory (RDMFT) which is one possible way to treat electron correlation beyond DFT. Within this theory the one-body reduced density matrix (1-RDM) is used as the basic variable. Our main interest is the calculation of the fundamental gap which proves very problematic within DFT. In order to calculate the fundamental gap we generalize RDMFT to fractional particle numbers M by describing the system as an ensemble of an N and an N+1 particle system (with N≤ M≤N+1). For each fixed particle number, M, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to the total energy as a function of M. The derivative of this function with respect to the particle number has a discontinuity at integer particle number which is identical to the gap. In addition, we investigate the necessary and sufficient conditions for the 1-RDM of a system with fractional particle number to be N-representable. Numerical results are presented for alkali atoms, small molecules, and periodic systems. Another problem within DFT is the description of non-relativistic many- electron systems in the presence of magnetic fields. It requires the paramagnetic current density and the spin magnetization to be used as basic variables besides the electron density. However, electron-gas-based functionals of current-spin-density-functional Theory (CSDFT) exhibit derivative discontinuities as a function of the magnetic field whenever a new Landau level is occupied, which makes them difficult to use in practice. Since the appearance of Landau levels is, intrinsically, an orbital effect it is appealing to use orbital-dependent functionals. We have developed a CSDFT version of the optimized effective potential (OEP) method which allows for the use of explicitly orbital-dependent functionals in the context of CSDFT. We present the derivation of corresponding equations and show results for a quantum dot in external magnetic fields in the second part of this thesis.
Accurately describing chemical bond breaking requires capturing strong/static electron correlation, a major challenge for single-reference self-consistent-field methods such as Hartree--Fock and Kohn--Sham density-functional theory (DFT). Fractional occupations substantially improve this situation, as demonstrated by information-density-matrix-functional theory (i-DMFT) and thermally-assisted-occupation DFT (TAO-DFT). Although useful, they are formulated within rather different theoretical frameworks, are somewhat semi-empirical, and their relationship remains unclear. Here we introduce finite-temperature adiabatic-connection reduced density-matrix functional theory (FT-AC-RDMFT), in which i-DMFT and TAO-DFT emerge as different reductions of a unified framework. The effective eigenvalue equation for natural orbitals and the Fermi--Dirac occupation numbers follow directly from the variational principle, while a thermal adiabatic connection relates the empirical entropy coefficient of i-DMFT to a characteristic temperature of the auxiliary non-interacting system. Within this framework, i-DMFT is the exact-exchange-based reduction, whereas TAO-DFT results when exact exchange is replaced by a density-based exchange-correlation functional. Numerical calculations confirm that i-DMFT, apart from a constant, coincides with the exact-exchange realization of the framework, and show that the improved dissociation behavior from finite-temperature occupations extends beyond exact exchange to hybrid and semilocal density functionals. For strongly correlated V2O3, finite-temperature occupations substantially reduce the severely overestimated non-magnetic LDA bulk modulus, bringing it close to experiment. These results establish FT-AC-RDMFT as a common theoretical foundation for finite-temperature occupation-based electronic-structure methods in molecular and solid-state systems.
In many cases, density-functional theory (DFT) with current standard approximate functionals offers a relatively accurate and computationally cheap description of the short-range dynamic electron correlation effects. However, in general, standard DFT does not treat the dispersion interaction effects adequately which, on the other hand, can be described by many-body perturbation theory (MBPT). It is therefore of interest to develop a hybrid model which combines the best of both the MBPT and DFT approaches. This can be achieved by splitting the two-electron interaction into long-range and short-range parts; the long-range part is then treated by MBPT and the short-range part by DFT. This work deals with the formulation of a general MBPT-DFT model (i.e., valid for any type of zeroth-order Hamiltonian) based on such a range separation. Applying the Rayleigh-Schr\"odinger formalism in this context, one finds that the generalized Bloch equation becomes self-consistent at each order of perturbation. This complication arises because the short-range part of the energy is a functional of the exact electron density, which is expanded in a perturbation series. In order to address this ``self-consistency problem'' and provide computable orbital-based expressions for any order of perturbation, a general one-electron reduced-density-matrix-based formalism is proposed. Two applications of our general formalism are presented: The derivation of a hybrid second-order M\o{}ller-Plesset-DFT model and the formulation of a range-separated optimized effective potential method based on a hybrid G\"orling-Levy-type MBPT-DFT model.
The cumulant expansion is a powerful approach for including correlation effects in electronic structure calculations beyond the GW approximation. However, the expansion is not generally valid, as current implementations ignore terms that mix particle and hole states and lead to partial occupation numbers of one-electron states. These limitations are corrected here using a cumulant expansion of the retarded one-electron Green's function that includes both particle and hole contributions. The approach provides a consistent framework that improves on the GW approximation to the spectral function without additional computational effort. The method is illustrated with results for the homogeneous electron gas and comparisons to experiment and other methods.