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:nucl-th/9911039 · Physical Review Letters

Response Function Technique for Calculating the Random-Phase Approximation Correlation Energy

Y. R. Shimizu, Paola Donati, R.A. Broglia

Published Sep 11, 2000• 16 citations
arXivHTMLPDF

Abstract

We develop a scheme to exactly evaluate the correlation energy in the random-phase approximation, based on linear response theory [Y. R. Shimizu, J. D. Garrett, R. A. Broglia, M. Gallardo, and E. Vigezzi, Rev. Mod. Phys. 61, 131 (1989)]. It is demonstrated that our formula is equivalent to a contour integral representation recently proposed [F. Donau, D. Almehed, and R. G. Nazmitdinov, Phys. Rev. Lett. 83, 280 (1999)] being numerically more efficient for realistic calculations. Examples are presented for pairing correlations in rapidly rotating nuclei.

DOI
10.1103/physrevlett.85.2260

Related papers

Integral Representation of the Random-Phase Approximation Correlation Energy

F. Dönau, D. Almehed, R. G. Nazmitdinov

Using the spectral function ${\ensuremath{\varphi}}^{\ensuremath{'}}(z)/\ensuremath{\varphi}(z)$, the random-phase approximation correlation energy and other properties of a finite system can be written as a contour integral in a compact way. This yields a transparent expression and reduces drastically the numerical efforts for obtaining reliable values. For illustration the method is applied to pairing vibrations in rotating nuclei.

arXiv:nucl-th/9902073•Physical Review Letters•Jul 12, 1999•18 citations

Linear-response time-dependent density-functional theory with pairing fields

Degao Peng, Helen van Aggelen, Yang Yang +1 more

Recent development in particle-particle random phase approximation (pp-RPA) broadens the perspective on ground state correlation energies [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013), Y. Yang, H. van Aggelen, S. N. Steinmann, D. Peng, and W. Yang, J. Chem. Phys. 139, 174110 (2013); D. Peng, S. N. Steinmann, H. van Aggelen, and W. Yang, J. Chem. Phys. 139, 104112 (2013)] and N ± 2 excitation energies [Y. Yang, H. van Aggelen, and W. Yang, J. Chem. Phys. 139, 224105 (2013)]. So far Hartree-Fock and approximated density-functional orbitals have been utilized to evaluate the pp-RPA equation. In this paper, to further explore the fundamentals and the potential use of pairing matrix dependent functionals, we present the linear-response time-dependent density-functional theory with pairing fields with both adiabatic and frequency-dependent kernels. This theory is related to the density-functional theory and time-dependent density-functional theory for superconductors, but is applied to normal non-superconducting systems for our purpose. Due to the lack of the proof of the one-to-one mapping between the pairing matrix and the pairing field for time-dependent systems, the linear-response theory is established based on the representability assumption of the pairing matrix. The linear response theory justifies the use of approximated density-functionals in the pp-RPA equation. This work sets the fundamentals for future density-functional development to enhance the description of ground state correlation energies and N ± 2 excitation energies.

The Journal of Chemical Physics•Mar 14, 2014•39 citations

Correlation energy from random phase approximations: A reduced density matrices perspective

Katarzyna Pernal

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.

International Journal of Quantum Chemistry•Aug 28, 2017•32 citations

Computational schemes for the ground-state pair density

Katsuhiko Higuchi, Masahiko Higuchi

We reconfirm the performance of the initial scheme for calculating the ground-state pair density (Higuchi and Higuchi 2007 Physica B 387 117, 2008 Phys. Rev. B 78 125101) by using the alternative approximation of the correlating kinetic energy functional. It is shown that about 20% of the correlation energy can be reproduced by the initial scheme, irrespective of the approximate form of the correlating kinetic energy functional. On the basis of the initial scheme, various kinds of schemes that go beyond the initial one can be developed. We illustrate two kinds of computational schemes.

Journal of Physics Condensed Matter•Jan 20, 2009•23 citations

Random phase approximation vs exact shell-model correlation energies

Ionel Stetcu, Calvin W. Johnson

The random phase approximation (RPA) builds in correlations left out by mean-field theory. In full $0\ensuremath{\Elzxh}\ensuremath{\omega}$ shell-model spaces we calculate the Hartree-Fock $(\mathrm{HF})+\mathrm{RPA}$ binding energy, and compare it to exact diagonalization. We find that often, but not always, the $\mathrm{HF}+\mathrm{RPA}$ gives a good approximation to the ``exact'' ground state energy. In those cases where the RPA is less satisfactory, however, there is no obvious correlation with properties of the HF state, such as deformation or overlap with the exact ground state wave function. This weakens the reliability of the RPA for computing binding energies.

arXiv:nucl-th/0205029•Physical Review C•Sep 3, 2002•37 citations

Convergence behavior of the random phase approximation renormalized correlation energy

Jefferson E. Bates, J. Sensenig, Adrienn Ruzsinszky

Based on the random phase approximation (RPA), RPA renormalization [J. E. Bates and F. Furche, J. Chem. Phys. 139, 171103 (2013)] is a robust many-body perturbation theory that works for molecules and materials because it does not diverge as the Kohn-Sham gap approaches zero. Additionally, RPA renormalization enables the simultaneous calculation of RPA and beyond-RPA correlation energies since the total correlation energy is the sum of a series of independent contributions. The first-order approximation (RPAr1) yields the dominant beyond-RPA contribution to the correlation energy for a given exchange-correlation kernel, but systematically underestimates the total beyond-RPA correction. For both the homogeneous electron gas model and real systems, we demonstrate numerically that RPA renormalization beyond first order converges monotonically to the infinite-order beyond-RPA correlation energy for several model exchange-correlation kernels and that the rate of convergence is principally determined by the choice of the kernel and spin polarization of the ground state. The monotonic convergence is rationalized from an analysis of the RPA renormalized correlation energy corrections, assuming the exchange-correlation kernel and response functions satisfy some reasonable conditions. For spin-unpolarized atoms, molecules, and bulk solids, we find that RPA renormalization is typically converged to 1 meV error or less by fourth order regardless of the band gap or dimensionality. Most spin-polarized systems converge at a slightly slower rate, with errors on the order of 10 meV at fourth order and typically requiring up to sixth order to reach 1 meV error or less. Slowest to converge, however, open-shell atoms present the most challenging case and require many higher orders to converge.

Physical review. B./Physical review. B•May 26, 2017•17 citations
We develop a scheme to exactly evaluate the correlation energy in the random-phase approximation, based on linear response theory [Y. R. Shimizu, J. D. Garrett, R. A. Broglia, M. Gallardo, and E. Vigezzi, Rev. Mod. Phys. 61, 131 (1989)]. It is demonstrated that our formula is equivalent to a contour integral representation recently proposed [F. Donau, D. Almehed, and R. G. Nazmitdinov, Phys. Rev. Lett. 83, 280 (1999)] being numerically more efficient for realistic calculations. Examples are presented for pairing correlations in rapidly rotating nuclei.
Recent development in particle-particle random phase approximation (pp-RPA) broadens the perspective on ground state correlation energies [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013), Y. Yang, H. van Aggelen, S. N. Steinmann, D. Peng, and W. Yang, J. Chem. Phys. 139, 174110 (2013); D. Peng, S. N. Steinmann, H. van Aggelen, and W. Yang, J. Chem. Phys. 139, 104112 (2013)] and N ± 2 excitation energies [Y. Yang, H. van Aggelen, and W. Yang, J. Chem. Phys. 139, 224105 (2013)]. So far Hartree-Fock and approximated density-functional orbitals have been utilized to evaluate the pp-RPA equation. In this paper, to further explore the fundamentals and the potential use of pairing matrix dependent functionals, we present the linear-response time-dependent density-functional theory with pairing fields with both adiabatic and frequency-dependent kernels. This theory is related to the density-functional theory and time-dependent density-functional theory for superconductors, but is applied to normal non-superconducting systems for our purpose. Due to the lack of the proof of the one-to-one mapping between the pairing matrix and the pairing field for time-dependent systems, the linear-response theory is established based on the representability assumption of the pairing matrix. The linear response theory justifies the use of approximated density-functionals in the pp-RPA equation. This work sets the fundamentals for future density-functional development to enhance the description of ground state correlation energies and N ± 2 excitation energies.
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 reconfirm the performance of the initial scheme for calculating the ground-state pair density (Higuchi and Higuchi 2007 Physica B 387 117, 2008 Phys. Rev. B 78 125101) by using the alternative approximation of the correlating kinetic energy functional. It is shown that about 20% of the correlation energy can be reproduced by the initial scheme, irrespective of the approximate form of the correlating kinetic energy functional. On the basis of the initial scheme, various kinds of schemes that go beyond the initial one can be developed. We illustrate two kinds of computational schemes.
The random phase approximation (RPA) builds in correlations left out by mean-field theory. In full $0\ensuremath{\Elzxh}\ensuremath{\omega}$ shell-model spaces we calculate the Hartree-Fock (HF)+RPA(\mathrm{HF})+\mathrm{RPA}(HF)+RPA binding energy, and compare it to exact diagonalization. We find that often, but not always, the HF+RPA\mathrm{HF}+\mathrm{RPA}HF+RPA gives a good approximation to the ``exact'' ground state energy. In those cases where the RPA is less satisfactory, however, there is no obvious correlation with properties of the HF state, such as deformation or overlap with the exact ground state wave function. This weakens the reliability of the RPA for computing binding energies.
Based on the random phase approximation (RPA), RPA renormalization [J. E. Bates and F. Furche, J. Chem. Phys. 139, 171103 (2013)] is a robust many-body perturbation theory that works for molecules and materials because it does not diverge as the Kohn-Sham gap approaches zero. Additionally, RPA renormalization enables the simultaneous calculation of RPA and beyond-RPA correlation energies since the total correlation energy is the sum of a series of independent contributions. The first-order approximation (RPAr1) yields the dominant beyond-RPA contribution to the correlation energy for a given exchange-correlation kernel, but systematically underestimates the total beyond-RPA correction. For both the homogeneous electron gas model and real systems, we demonstrate numerically that RPA renormalization beyond first order converges monotonically to the infinite-order beyond-RPA correlation energy for several model exchange-correlation kernels and that the rate of convergence is principally determined by the choice of the kernel and spin polarization of the ground state. The monotonic convergence is rationalized from an analysis of the RPA renormalized correlation energy corrections, assuming the exchange-correlation kernel and response functions satisfy some reasonable conditions. For spin-unpolarized atoms, molecules, and bulk solids, we find that RPA renormalization is typically converged to 1 meV error or less by fourth order regardless of the band gap or dimensionality. Most spin-polarized systems converge at a slightly slower rate, with errors on the order of 10 meV at fourth order and typically requiring up to sixth order to reach 1 meV error or less. Slowest to converge, however, open-shell atoms present the most challenging case and require many higher orders to converge.