Monte Carlo explicitly correlated many-body Green's function theory

Cole M. Johnson, Alexander E. Doran, Seiichiro L. Ten-No, So Hirata

Research output: Contribution to journalArticlepeer-review


A highly scalable stochastic algorithm is proposed and implemented for computing the basis-set-incompleteness correction to the diagonal, frequency-independent self-energy of the second-order many-body Green's function (GF2) theory within the explicitly correlated (F12) formalism. The 6-, 9-, 12-, and 15-dimensional integrals comprising the F12 correction are directly evaluated by the Monte Carlo method using appropriate weight functions for importance sampling. The method is naturally and easily parallelized, involves minimal memory space and no disk I/O, and can use virtually any mathematical form of a correlation factor. Its computational cost to correct all ionization energies (IEs) is observed to increase as the fourth power of system size, as opposed to the fifth power in the case of the deterministic counterparts. The GF2 calculations and their F12 corrections for the first IEs of C60 and C70 were executed on 128 graphical processing units (GF2) and 896 central processing units (F12), respectively, to reach the results with statistical errors of 0.04 eV or less. They showed that the basis-set-incompleteness (from aug-cc-pVDZ) accounts for only 50%-60% of the deviations from experiments, suggesting the significance of higher-order perturbation corrections.

Original languageEnglish (US)
Article number174112
JournalJournal of Chemical Physics
Issue number17
StatePublished - Nov 7 2018

ASJC Scopus subject areas

  • General Physics and Astronomy
  • Physical and Theoretical Chemistry


Dive into the research topics of 'Monte Carlo explicitly correlated many-body Green's function theory'. Together they form a unique fingerprint.

Cite this