Convergence acceleration of Monte Carlo many-body perturbation methods by direct sampling

Alexander E Doran, So Hirata

Research output: Contribution to journalArticlepeer-review


In the Monte Carlo many-body perturbation (MC-MP) method, the conventional correlation-correction formula, which is a long sum of products of low-dimensional integrals, is first recast into a short sum of high-dimensional integrals over electron-pair and imaginary-time coordinates. These high-dimensional integrals are then evaluated by the Monte Carlo method with random coordinates generated by the Metropolis-Hasting algorithm according to a suitable distribution. The latter algorithm, while advantageous in its ability to sample nearly any distribution, introduces autocorrelation in sampled coordinates, which, in turn, increases the statistical uncertainty of the integrals and thus the computational cost. It also involves wasteful rejected moves and an initial "burn-in" step as well as displays hysteresis. Here, an algorithm is proposed that directly produces a random sequence of electron-pair coordinates for the same distribution used in the MC-MP method, which is free from autocorrelation, rejected moves, a burn-in step, or hysteresis. This direct-sampling algorithm is shown to accelerate second- and third-order Monte Carlo many-body perturbation calculations by up to 222% and 38%, respectively.

Original languageEnglish (US)
Article number104112
Pages (from-to)104112
JournalThe Journal of Chemical Physics
Issue number10
StatePublished - Sep 14 2020

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry


Dive into the research topics of 'Convergence acceleration of Monte Carlo many-body perturbation methods by direct sampling'. Together they form a unique fingerprint.

Cite this