Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1409.0354 · Physical Review B

Correlation energy within exact-exchange adiabatic connection fluctuation-dissipation theory: Systematic development and simple approximations

Nicola Colonna, Maria Hellgren, Stefano de Gironcoli

Published Sep 26, 2014• 47 citations
arXivHTMLPDF

Abstract

We have calculated the correlation energy of the homogeneous electron gas (HEG) and the dissociation energy curves of molecules with covalent bonds from an efficient implementation of the adiabatic connection fluctuation dissipation expression including the exact-exchange (EXX) kernel. The EXX kernel is defined from first-order perturbation theory and used in the Dyson equation of time-dependent density-functional theory. Within this approximation (RPAx), the correlation energies of the HEG are significantly improved with respect to the random phase approximation (RPA) up to densities of the order of ${r}_{s}\ensuremath{\approx}10$. However, beyond this value, the RPAx response function exhibits an unphysical divergence and the approximation breaks down. Total energies of molecules at equilibrium are also highly accurate, but we find a similar instability at stretched geometries. Staying within an exact first-order approximation to the response function, we use an alternative resummation of the higher-order terms. This slight redefinition of RPAx fixes the instability in total energy calculations without compromising the overall accuracy of the approach.

DOI
10.1103/physrevb.90.125150

Related papers

Tension between predicting accurate ground state correlation energies and excitation energies from adiabatic approximations in TDDFT

Lucas M. Everhart, Julio A. Derteano, Jefferson E. Bates

The connection between the adiabatic excitation energy of time-dependent density functional theory and the ground state correlation energy from the adiabatic connection fluctuation-dissipation theorem (ACFDT) is explored in the limiting case of one excited state. An exact expression is derived for any adiabatic Hartree-exchange-correlation kernel that connects the excitation energy and the potential contribution to correlation. The resulting formula is applied to the asymmetric Hubbard dimer, a system where this limit is exact. Results from a hierarchy of approximations to the kernel, including the random phase approximation (RPA) with and without exchange and the adiabatically exact (AE) approximation, are compared to the exact ones. At full coupling, the numerical results indicate a tension between predicting an accurate excitation energy and an accurate potential contribution to correlation. The AE approximation is capable of making accurate predictions of both quantities, but only in parts of the parameter space that classify as weakly correlated, while RPA tends to be unable to accurately predict these properties simultaneously anywhere. For a strongly correlated dimer, the AE approximation greatly overestimates the excitation energy yet continues to yield an accurate ground state correlation energy due to its accurate prediction of the adiabatic connection integrand. If similar trends hold for real systems, the development of correlation kernels will be important for applications of the ACFDT in systems with large potential contributions to correlation.

The Journal of Chemical Physics•Feb 25, 2022•0 citations

Accurate total energies from the adiabatic-connection fluctuation-dissipation theorem

Nick Woods, Mike Entwistle, Rex W. Godby

In the context of inhomogeneous one-dimensional finite systems, recent numerical advances [Phys. Rev. B 103, 125155 (2021)] allow us to compute the exact coupling-constant dependent exchange-correlation kernel ${f}_{\text{xc}}^{\ensuremath{\lambda}}(x,{x}^{\ensuremath{'}},\ensuremath{\omega})$ within linear response time-dependent density-functional theory. This permits an improved understanding of ground-state total energies derived from the adiabatic-connection fluctuation-dissipation theorem (ACFDT). We consider both one-shot and self-consistent ACFDT calculations, and demonstrate that chemical accuracy is reliably preserved when the frequency dependence in the exact functional ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ is neglected. This performance is understood on the grounds that the exact ${f}_{\text{xc}}[n]$ varies slowly over the most relevant $\ensuremath{\omega}$ range (but not in general), and hence the spatial structure in ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ is able to largely remedy the principal issue in the present context: self-interaction (examined from the perspective of the exchange-correlation hole). Moreover, we find that the implicit orbitals contained within a self-consistent ACFDT calculation utilizing the adiabatic exact kernel ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ are remarkably similar to the exact Kohn-Sham orbitals, thus further establishing that the majority of the physics required to capture the ground-state total energy resides in the spatial dependence of ${f}_{\text{xc}}[n]$ at $\ensuremath{\omega}=0$.

arXiv:2107.14434•Physical review. B./Physical review. B•Sep 17, 2021•4 citations

Adiabatic Connection without Coupling Constant Integration

Jefferson E. Bates, Niladri Sengupta, J. Sensenig +1 more

Using a second-order approximation to Random Phase Approximation renormalized (RPAr) many-body perturbation theory for the interacting density-density response function, we have developed a so-called higher-order terms (HOT) approximation for the correlation energy. In combination with the first-order RPAr correction, our new method faithfully captures the infinite-order correlation for a given exchange-correlation kernel, yielding errors of the total correlation energy on the order of 1% or less for most systems. For exchange-like kernels, our new method has the further benefit that the coupling-strength integration can be completely eliminated resulting in a modest reduction in computational cost compared to the traditional approach. When the correlation energy is accurately reproduced by the HOT approximation, structural properties and energy differences are also accurately reproduced, as we demonstrate for several periodic solids and some molecular systems. Energy differences involving fragmentation are challenging for the HOT method, however, due to errors that may not cancel between a composite system and its constituent pieces.

Journal of Chemical Theory and Computation•May 7, 2018•14 citations

Correlation beyond the random phase approximation: A consistent many-body perturbation theory approach

F. Bechstedt

The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function $L$. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential $W$ in the BSE for $L$ leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.

Physical review. B./Physical review. B•Jun 15, 2018•5 citations

Exchange-correlation potentials in the adiabatic connection fluctuation-dissipation framework

Yann‐Michel Niquet, Martin Fuchs, Xavier Gonze

We provide the expression of the exchange-correlation potential in the adiabatic connection fluctuation-dissipation (ACFD) framework, for arbitrary time-dependent (TD) kernels. We investigate the asymptotic behavior of the ACFD potential in three relevant approximations: the random-phase approximation, the exact-exchange kernel in two-electron systems, and the adiabatic local-density approximation. We show that these potentials have the expected $\ensuremath{-}{1/r+Q/r}^{3}\ensuremath{-}\ensuremath{\alpha}{/(2r}^{4})$ tail (in closed-shell systems with spherical symmetry), where Q and $\ensuremath{\alpha}$ depend on the TD kernel and reflect the physics included in each approximation. We also discuss approximate ACFD potentials that are much simpler to compute than the exact ones while being likely of reasonable accuracy.

Physical Review A•Sep 11, 2003•89 citations

Development of low-scaling methods to calculate ground state energies and analytical gradients based on the adiabatic-connection fluctuation-dissipation theorem

Matthias Beuerle

New methods to efficiently calculate energetics and first order-properties for mean-field and correlated electronic structure theories are presented. In linear-scaling short-range hybrid calculations new integral screening criteria exploiting the locality of the attenuated Coulomb operator are introduced that allow to significantly increase the performance of these density functionals. This enables short-range hybrid calculations at a similar cost as pure semi-local density functional theory (DFT) calculations with increased accuracy due to the admixture of exchange. At the level of correlated electronic structure theory, the realm of systems which can be calculated with a linear-scaling random-phase-approximation (RPA) method is extended with a novel multi-node parallel algorithm. Furthermore, more efficient quadrature schemes are used based on the equivalence of the employed integral transforms with the Fourier transform of the non-interacting polarizability. In combination with a new way to introduce Cholesky orbitals this results in a more effient and numerically more accurate linear-scaling RPA correlation energy method with improved memory requirements. Combining these techniques further with an approach to calculate analytical first order properties using quantities from many body perturbation theory, allows to present a low-scaling method that gives access to RPA gradients for molecular systems with several hundred atoms. Moreover linear- and low-scaling methods to calculate different beyond RPA correlation energies are devised, exploiting the locality of exchange type contractions. This enables to apply these more accurate RPA methods to significantly larger systems, which allows to demonstrate the gain in accuracy for large, dispersion dominated systems. The newly developed methods extend the scope of RPA and RPA with exchange correlation energies, and first order RPA properties significantly, while the accuracy is under full numerical control. Furthermore, there is no overhead to the respective canonical algorithm, making the presented methods competitive also for small system sizes. Finally, a new beyond RPA scheme is presented that combines the benefits of plain RPA and RPA with second-order screened exchange, leading to a more balanced, highly accurate post Kohn-Sham method.

•Jan 1, 2019•0 citations
We have calculated the correlation energy of the homogeneous electron gas (HEG) and the dissociation energy curves of molecules with covalent bonds from an efficient implementation of the adiabatic connection fluctuation dissipation expression including the exact-exchange (EXX) kernel. The EXX kernel is defined from first-order perturbation theory and used in the Dyson equation of time-dependent density-functional theory. Within this approximation (RPAx), the correlation energies of the HEG are significantly improved with respect to the random phase approximation (RPA) up to densities of the order of rs\ensuremath≈10{r}_{s}\ensuremath{\approx}10rs​\ensuremath≈10. However, beyond this value, the RPAx response function exhibits an unphysical divergence and the approximation breaks down. Total energies of molecules at equilibrium are also highly accurate, but we find a similar instability at stretched geometries. Staying within an exact first-order approximation to the response function, we use an alternative resummation of the higher-order terms. This slight redefinition of RPAx fixes the instability in total energy calculations without compromising the overall accuracy of the approach.
The connection between the adiabatic excitation energy of time-dependent density functional theory and the ground state correlation energy from the adiabatic connection fluctuation-dissipation theorem (ACFDT) is explored in the limiting case of one excited state. An exact expression is derived for any adiabatic Hartree-exchange-correlation kernel that connects the excitation energy and the potential contribution to correlation. The resulting formula is applied to the asymmetric Hubbard dimer, a system where this limit is exact. Results from a hierarchy of approximations to the kernel, including the random phase approximation (RPA) with and without exchange and the adiabatically exact (AE) approximation, are compared to the exact ones. At full coupling, the numerical results indicate a tension between predicting an accurate excitation energy and an accurate potential contribution to correlation. The AE approximation is capable of making accurate predictions of both quantities, but only in parts of the parameter space that classify as weakly correlated, while RPA tends to be unable to accurately predict these properties simultaneously anywhere. For a strongly correlated dimer, the AE approximation greatly overestimates the excitation energy yet continues to yield an accurate ground state correlation energy due to its accurate prediction of the adiabatic connection integrand. If similar trends hold for real systems, the development of correlation kernels will be important for applications of the ACFDT in systems with large potential contributions to correlation.
In the context of inhomogeneous one-dimensional finite systems, recent numerical advances [Phys. Rev. B 103, 125155 (2021)] allow us to compute the exact coupling-constant dependent exchange-correlation kernel ${f}_{\text{xc}}^{\ensuremath{\lambda}}(x,{x}^{\ensuremath{'}},\ensuremath{\omega})$ within linear response time-dependent density-functional theory. This permits an improved understanding of ground-state total energies derived from the adiabatic-connection fluctuation-dissipation theorem (ACFDT). We consider both one-shot and self-consistent ACFDT calculations, and demonstrate that chemical accuracy is reliably preserved when the frequency dependence in the exact functional ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ is neglected. This performance is understood on the grounds that the exact fxc[n]{f}_{\text{xc}}[n]fxc​[n] varies slowly over the most relevant $\ensuremath{\omega}$ range (but not in general), and hence the spatial structure in ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ is able to largely remedy the principal issue in the present context: self-interaction (examined from the perspective of the exchange-correlation hole). Moreover, we find that the implicit orbitals contained within a self-consistent ACFDT calculation utilizing the adiabatic exact kernel ${f}_{\text{xc}}[n](\ensuremath{\omega}=0)$ are remarkably similar to the exact Kohn-Sham orbitals, thus further establishing that the majority of the physics required to capture the ground-state total energy resides in the spatial dependence of fxc[n]{f}_{\text{xc}}[n]fxc​[n] at $\ensuremath{\omega}=0$.
Using a second-order approximation to Random Phase Approximation renormalized (RPAr) many-body perturbation theory for the interacting density-density response function, we have developed a so-called higher-order terms (HOT) approximation for the correlation energy. In combination with the first-order RPAr correction, our new method faithfully captures the infinite-order correlation for a given exchange-correlation kernel, yielding errors of the total correlation energy on the order of 1% or less for most systems. For exchange-like kernels, our new method has the further benefit that the coupling-strength integration can be completely eliminated resulting in a modest reduction in computational cost compared to the traditional approach. When the correlation energy is accurately reproduced by the HOT approximation, structural properties and energy differences are also accurately reproduced, as we demonstrate for several periodic solids and some molecular systems. Energy differences involving fragmentation are challenging for the HOT method, however, due to errors that may not cancel between a composite system and its constituent pieces.
The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function LLL. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential WWW in the BSE for LLL leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.
We provide the expression of the exchange-correlation potential in the adiabatic connection fluctuation-dissipation (ACFD) framework, for arbitrary time-dependent (TD) kernels. We investigate the asymptotic behavior of the ACFD potential in three relevant approximations: the random-phase approximation, the exact-exchange kernel in two-electron systems, and the adiabatic local-density approximation. We show that these potentials have the expected $\ensuremath{-}{1/r+Q/r}^{3}\ensuremath{-}\ensuremath{\alpha}{/(2r}^{4})$ tail (in closed-shell systems with spherical symmetry), where Q and $\ensuremath{\alpha}$ depend on the TD kernel and reflect the physics included in each approximation. We also discuss approximate ACFD potentials that are much simpler to compute than the exact ones while being likely of reasonable accuracy.
New methods to efficiently calculate energetics and first order-properties for mean-field and correlated electronic structure theories are presented. In linear-scaling short-range hybrid calculations new integral screening criteria exploiting the locality of the attenuated Coulomb operator are introduced that allow to significantly increase the performance of these density functionals. This enables short-range hybrid calculations at a similar cost as pure semi-local density functional theory (DFT) calculations with increased accuracy due to the admixture of exchange. At the level of correlated electronic structure theory, the realm of systems which can be calculated with a linear-scaling random-phase-approximation (RPA) method is extended with a novel multi-node parallel algorithm. Furthermore, more efficient quadrature schemes are used based on the equivalence of the employed integral transforms with the Fourier transform of the non-interacting polarizability. In combination with a new way to introduce Cholesky orbitals this results in a more effient and numerically more accurate linear-scaling RPA correlation energy method with improved memory requirements. Combining these techniques further with an approach to calculate analytical first order properties using quantities from many body perturbation theory, allows to present a low-scaling method that gives access to RPA gradients for molecular systems with several hundred atoms. Moreover linear- and low-scaling methods to calculate different beyond RPA correlation energies are devised, exploiting the locality of exchange type contractions. This enables to apply these more accurate RPA methods to significantly larger systems, which allows to demonstrate the gain in accuracy for large, dispersion dominated systems. The newly developed methods extend the scope of RPA and RPA with exchange correlation energies, and first order RPA properties significantly, while the accuracy is under full numerical control. Furthermore, there is no overhead to the respective canonical algorithm, making the presented methods competitive also for small system sizes. Finally, a new beyond RPA scheme is presented that combines the benefits of plain RPA and RPA with second-order screened exchange, leading to a more balanced, highly accurate post Kohn-Sham method.