Abstract Random phase approximation (RPA) electron correlation methods have gained in popularity in the recent years. A number of RPA correlation energy variants emerged in the Kohn–Sham DFT framework and in the theory of strongly orthogonal geminals. Foundations of most RPA approaches trace back to an exact expression for two‐electron reduced density matrix written in terms of one‐electron density matrix and dynamic one‐electron response functions, originally presented in the seminal paper of McLachlan and Ball (Rev. Mod. Phys. 1964, 36, 844). The aim of this article is to give a pedagogical introduction to possible approaches for describing electron correlation energy based on the McLachlan and Ball formula. The focus of the presentation is to formulate electron‐interaction energy expressions as functions of reduced density matrices. On one hand, it provides a common umbrella for RPA approximations proposed for uncorrelated (Hartree–Fock or Kohn–Sham) as well as partially correlated (strongly orthogonal geminals) references. On the other hand, such presentation may stimulate new developments in density matrix functional theory.
Abstract Random phase approximation (RPA) electron correlation methods have gained in popularity in the recent years. A number of RPA correlation energy variants emerged in the Kohn–Sham DFT framework and in the theory of strongly orthogonal geminals. Foundations of most RPA approaches trace back to an exact expression for two‐electron reduced density matrix written in terms of one‐electron density matrix and dynamic one‐electron response functions, originally presented in the seminal paper of McLachlan and Ball (Rev. Mod. Phys. 1964, 36, 844). The aim of this article is to give a pedagogical introduction to possible approaches for describing electron correlation energy based on the McLachlan and Ball formula. The focus of the presentation is to formulate electron‐interaction energy expressions as functions of reduced density matrices. On one hand, it provides a common umbrella for RPA approximations proposed for uncorrelated (Hartree–Fock or Kohn–Sham) as well as partially correlated (strongly orthogonal geminals) references. On the other hand, such presentation may stimulate new developments in density matrix functional theory.
We calculate the interacting bandgap energy of a solid within the random-phase approximation (RPA) to density functional theory (DFT). The interacting bandgap energy is defined as ${E}_{g}={E}^{\mathrm{RPA}}(N+1)+{E}^{\mathrm{RPA}}(N\ensuremath{-}1)\ensuremath{-}2{E}^{\mathrm{RPA}}(N)$, where ERPA(N) is the total DFT-RPA energy of the N-electron system. We compare the interacting bandgap energy with the Kohn-Sham bandgap energy, which is the difference between the conduction and valence band edges in the Kohn-Sham band structure. We show that they differ by an unrenormalized ``G0W0'' self-energy correction (i.e., a GW self-energy correction computed using Kohn-Sham orbitals and energies as input). This provides a well-defined and meaningful interpretation to G0W0 quasiparticle bandgap calculations, but questions the physics behind the renormalization factors in the expression of the bandgap energy. We also separate the kinetic from the Coulomb contributions to the DFT-RPA bandgap energy, and discuss the related problem of the derivative discontinuity in the DFT-RPA functional. Last we discuss the applicability of our results to other functionals based on many-body perturbation theory.
The well-known and well-used Density Functional Theory (DFT) is a variational theory for the ground-state. There are no guarantees that it delivers good results for excited state properties such as quasiparticle energies or wave function. All that is guaranteed is that the total energy is exact, the chemical potential is correct, and that the eective electronic states (Kohn-Sham states) generate the right electron density at the variational minimum. Luttinger-Ward approaches to electronic structure are attractive as beyond-DFT schemes since, in principle, they oer a variational functional of the one-electron Green’s function that yields both the correct total energy and all one-particle properties at the variational extremum. However, compared to DFT, these methods face a number of additional challenges beyond the familiar problem of choosing an exchange-correlation functional: (a) they have a much greater computational cost compared to DFT, (b) one must choose some set of trial Green’s functions, and (c) one must locate an extremum (not necessarily a minimum). This talk aims both to be pedagogical and to provide a summary of our most recent work in this area. We focus mainly on the random-phase approximation (RPA) for the correlation functional which yields the GW self-energy for the Green’s function. We will describe ecient approaches for computing the RPA correlation energy; how various types of approximations to the correlation energy create a ladder of approximate self-energy operators and what type of physics is included in these approximations (exchange, screening, Hubbard U, etc.); and how the extremization is in fact highly problematic. We will end with some ideas for making practical progress with the extremization and their relation to the successful Quasiparticle Self-consistent GW (QSGW) scheme [1].
The correlation energy in the direct random phase approximation (dRPA) can be written, among other possibilities, either in terms of the interaction strength averaged correlation density matrix, or in terms of the coupled cluster doubles amplitudes obtained in the direct ring approximation (drCCD). Although the corresponding dRPA correlation density matrix on the one hand, and the drCCD amplitude matrix on the other hand, differ significantly, they yield identical energies. Similarly, the analogous RPA and rCCD correlation energies calculated from antisymmetrized two-electron integrals are identical to each other despite very different underlying working equations. In the present communication, a direct correspondence between amplitudes and densities is established and investigated with perturbation theory arguments. Our analysis also sheds some light on the properties of recently proposed RPA/rCCD variants which use antisymmetrized integrals in part of the equations and nonantisymmetrized integrals in others.
Abstract Recently a novel approach to find approximate exchange–correlation functionals in density-functional theory was presented (Mordovina et al 2019 J. Chem. Theory Comput. 15 5209), which relies on approximations to the interacting wave function using density-matrix embedding theory (DMET). This approximate interacting wave function is constructed by using a projection determined by an iterative procedure that makes parts of the reduced density matrix of an auxiliary system the same as the approximate interacting density matrix. If only the diagonal of both systems are connected this leads to an approximation of the interacting-to-non-interacting mapping of the Kohn–Sham approach to density-functional theory. Yet other choices are possible and allow to connect DMET with other density-functional theories such as kinetic-energy density functional theory or reduced density-matrix functional theory. In this work we give a detailed review of the basics of the DMET procedure from a density-functional perspective and show how both approaches can be used to supplement each other. We do not present a specific realization of combining density-functional methods with DMET but rather provide common grounds to facilitate future developments that encompass both approaches. We do so explicitly for the case of a one-dimensional lattice system, as this is the simplest setting where we can apply DMET and the one that was originally presented. Among others we highlight how the mappings of density-functional theories can be used to identify uniquely defined auxiliary systems and projections in DMET and how to construct approximations for different density-functional theories using DMET inspired projections. Such alternative approximation strategies become especially important for density-functional theories that are based on non-linearly coupled observables such as kinetic-energy density-functional theory, where the Kohn–Sham fields are no longer obtainable by functional differentiation of an energy expression, or for reduced density-matrix functional theories, where a straightforward Kohn–Sham construction is not feasible.
Local correlation methods rely on the assumption that electronic correlation is nearsighted. In this work, we develop a method to alleviate this assumption. The first step is to approximately decompose the electron correlation to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation by Langreth and Perdew. The nearsighted component is then calculated using the recently developed embedded cluster density approximation (ECDA) which is a local correlation method formulated in the context of density functional theory. The farsighted component is calculated based on the system’s Kohn-Sham orbitals. The accuracy of this new method depends on the quality of the decomposition. We examined the method’s accuracy by patching the random phase approximation (RPA) correlation energy in a H 24 chain in which the electron correlation is highly nonlocal. This new method predicts bond stretching energies, RPA correlation potential, and Kohn-Sham eigenvalues in good agreement with the benchmarks. Our results demonstrate the importance of including the farsighted part of electron correlation for studying systems having nonlocal correlations.
Abstract Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one‐dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short‐wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long‐wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn‐Sham orbitals. This new method is termed λ‐ECDA, where λ indicates the wavelength decomposition. The performance of λ‐ECDA is examined on a one‐dimensional model system: a H24 chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron‐electron and electron‐ion interactions, and slightly stronger nuclear charges (1.2e) are assigned to the pseudo‐H atoms. Bond stretching energies, RPA correlation potentials, and Kohn‐Sham eigenvalues predicted by λ‐ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy.