Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integration

Hendrik Eshuis, Julian Yarkony, Filipp Furche

Research output: Contribution to journalArticleResearchpeer-review

136 Citations (Scopus)

Abstract

The random phase approximation (RPA) is an increasingly popular post-Kohn-Sham correlation method, but its high computational cost has limited molecular applications to systems with few atoms. Here we present an efficient implementation of RPA correlation energies based on a combination of resolution of the identity (RI) and imaginary frequency integration techniques. We show that the RI approximation to four-index electron repulsion integrals leads to a variational upper bound to the exact RPA correlation energy if the Coulomb metric is used. Auxiliary basis sets optimized for second-order Møller-Plesset (MP2) calculations are well suitable for RPA, as is demonstrated for the HEAT [A. Tajti, J. Chem. Phys. 121, 11599 (2004)] and MOLEKEL [F. Weigend, Chem. Phys. Lett. 294, 143 (1998)] benchmark sets. Using imaginary frequency integration rather than diagonalization to compute the matrix square root necessary for RPA, evaluation of the RPA correlation energy requires O (N4 log N) operations and O (N3) storage only; the price for this dramatic improvement over existing algorithms is a numerical quadrature. We propose a numerical integration scheme that is exact in the two-orbital case and converges exponentially with the number of grid points. For most systems, 30-40 grid points yield μH accuracy in triple zeta basis sets, but much larger grids are necessary for small gap systems. The lowest-order approximation to the present method is a post-Kohn-Sham frequency-domain version of opposite-spin Laplace-transform RI-MP2 [J. Jung, Phys. Rev. B 70, 205107 (2004)]. Timings for polyacenes with up to 30 atoms show speed-ups of two orders of magnitude over previous implementations. The present approach makes it possible to routinely compute RPA correlation energies of systems well beyond 100 atoms, as is demonstrated for the octapeptide angiotensin II.

Original languageEnglish
Article number234114
JournalJournal of Chemical Physics
Volume132
Issue number23
DOIs
StatePublished - 21 Jun 2010

Fingerprint

Atoms
approximation
Aromatic polymers
Correlation methods
Laplace transforms
Angiotensin II
energy
grids
Electrons
angiotensins
atoms
Costs
yield point
numerical integration
quadratures
time measurement
costs
orbitals
evaluation
matrices

Cite this

@article{ff784833f62248d2a5515bd4b289ce37,
title = "Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integration",
abstract = "The random phase approximation (RPA) is an increasingly popular post-Kohn-Sham correlation method, but its high computational cost has limited molecular applications to systems with few atoms. Here we present an efficient implementation of RPA correlation energies based on a combination of resolution of the identity (RI) and imaginary frequency integration techniques. We show that the RI approximation to four-index electron repulsion integrals leads to a variational upper bound to the exact RPA correlation energy if the Coulomb metric is used. Auxiliary basis sets optimized for second-order M{\o}ller-Plesset (MP2) calculations are well suitable for RPA, as is demonstrated for the HEAT [A. Tajti, J. Chem. Phys. 121, 11599 (2004)] and MOLEKEL [F. Weigend, Chem. Phys. Lett. 294, 143 (1998)] benchmark sets. Using imaginary frequency integration rather than diagonalization to compute the matrix square root necessary for RPA, evaluation of the RPA correlation energy requires O (N4 log N) operations and O (N3) storage only; the price for this dramatic improvement over existing algorithms is a numerical quadrature. We propose a numerical integration scheme that is exact in the two-orbital case and converges exponentially with the number of grid points. For most systems, 30-40 grid points yield μH accuracy in triple zeta basis sets, but much larger grids are necessary for small gap systems. The lowest-order approximation to the present method is a post-Kohn-Sham frequency-domain version of opposite-spin Laplace-transform RI-MP2 [J. Jung, Phys. Rev. B 70, 205107 (2004)]. Timings for polyacenes with up to 30 atoms show speed-ups of two orders of magnitude over previous implementations. The present approach makes it possible to routinely compute RPA correlation energies of systems well beyond 100 atoms, as is demonstrated for the octapeptide angiotensin II.",
author = "Hendrik Eshuis and Julian Yarkony and Filipp Furche",
year = "2010",
month = "6",
day = "21",
doi = "10.1063/1.3442749",
language = "English",
volume = "132",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics",
number = "23",

}

Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integration. / Eshuis, Hendrik; Yarkony, Julian; Furche, Filipp.

In: Journal of Chemical Physics, Vol. 132, No. 23, 234114, 21.06.2010.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integration

AU - Eshuis, Hendrik

AU - Yarkony, Julian

AU - Furche, Filipp

PY - 2010/6/21

Y1 - 2010/6/21

N2 - The random phase approximation (RPA) is an increasingly popular post-Kohn-Sham correlation method, but its high computational cost has limited molecular applications to systems with few atoms. Here we present an efficient implementation of RPA correlation energies based on a combination of resolution of the identity (RI) and imaginary frequency integration techniques. We show that the RI approximation to four-index electron repulsion integrals leads to a variational upper bound to the exact RPA correlation energy if the Coulomb metric is used. Auxiliary basis sets optimized for second-order Møller-Plesset (MP2) calculations are well suitable for RPA, as is demonstrated for the HEAT [A. Tajti, J. Chem. Phys. 121, 11599 (2004)] and MOLEKEL [F. Weigend, Chem. Phys. Lett. 294, 143 (1998)] benchmark sets. Using imaginary frequency integration rather than diagonalization to compute the matrix square root necessary for RPA, evaluation of the RPA correlation energy requires O (N4 log N) operations and O (N3) storage only; the price for this dramatic improvement over existing algorithms is a numerical quadrature. We propose a numerical integration scheme that is exact in the two-orbital case and converges exponentially with the number of grid points. For most systems, 30-40 grid points yield μH accuracy in triple zeta basis sets, but much larger grids are necessary for small gap systems. The lowest-order approximation to the present method is a post-Kohn-Sham frequency-domain version of opposite-spin Laplace-transform RI-MP2 [J. Jung, Phys. Rev. B 70, 205107 (2004)]. Timings for polyacenes with up to 30 atoms show speed-ups of two orders of magnitude over previous implementations. The present approach makes it possible to routinely compute RPA correlation energies of systems well beyond 100 atoms, as is demonstrated for the octapeptide angiotensin II.

AB - The random phase approximation (RPA) is an increasingly popular post-Kohn-Sham correlation method, but its high computational cost has limited molecular applications to systems with few atoms. Here we present an efficient implementation of RPA correlation energies based on a combination of resolution of the identity (RI) and imaginary frequency integration techniques. We show that the RI approximation to four-index electron repulsion integrals leads to a variational upper bound to the exact RPA correlation energy if the Coulomb metric is used. Auxiliary basis sets optimized for second-order Møller-Plesset (MP2) calculations are well suitable for RPA, as is demonstrated for the HEAT [A. Tajti, J. Chem. Phys. 121, 11599 (2004)] and MOLEKEL [F. Weigend, Chem. Phys. Lett. 294, 143 (1998)] benchmark sets. Using imaginary frequency integration rather than diagonalization to compute the matrix square root necessary for RPA, evaluation of the RPA correlation energy requires O (N4 log N) operations and O (N3) storage only; the price for this dramatic improvement over existing algorithms is a numerical quadrature. We propose a numerical integration scheme that is exact in the two-orbital case and converges exponentially with the number of grid points. For most systems, 30-40 grid points yield μH accuracy in triple zeta basis sets, but much larger grids are necessary for small gap systems. The lowest-order approximation to the present method is a post-Kohn-Sham frequency-domain version of opposite-spin Laplace-transform RI-MP2 [J. Jung, Phys. Rev. B 70, 205107 (2004)]. Timings for polyacenes with up to 30 atoms show speed-ups of two orders of magnitude over previous implementations. The present approach makes it possible to routinely compute RPA correlation energies of systems well beyond 100 atoms, as is demonstrated for the octapeptide angiotensin II.

UR - http://www.scopus.com/inward/record.url?scp=77953864922&partnerID=8YFLogxK

U2 - 10.1063/1.3442749

DO - 10.1063/1.3442749

M3 - Article

VL - 132

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 23

M1 - 234114

ER -