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:2602.08265 · Physical Review C

Stochastic many-body perturbation theory for high-order calculations

X. Zhen, R. Z. Hu, Junchen Pei, F. R. Xu

Published Apr 7, 2026• 0 citations
arXivHTMLPDF

Abstract

High-order perturbative $\textit{ab initio}$ calculations are challenging due to the rapidly growing configuration space and the difficulty of assessing convergence. In this letter, we introduce perturbation theory quantum Monte Carlo (PTQMC), a stochastic approach designed to compute high-order many-body perturbative corrections. By representing the perturbative wave function with random walkers in configuration space, PTQMC avoids the exponential scaling inherent to conventional constructions of high-rank excitation operators. Benchmark calculations for the Richardson pairing model demonstrate that PTQMC accurately reproduces exact many-body perturbation theory (MBPT) coefficients up to 16th order, even in strongly divergent regimes. We further show that combining PTQMC with series resummation techniques yields stable and precise energy estimates in cases where the straightforward perturbative series fails. Finally, we propose the effective number of configurations, $e^{S}$, as a global measure of perturbative wave-function complexity that can be directly extracted within PTQMC. We demonstrate that the saturation behavior of $e^{S}$ provides a more reliable indicator of the validity of perturbative expansions than energy convergence alone.

DOI
10.1103/q3vn-8y8s

Related papers

Stochastic many-body perturbation theory for electron correlation energies

Zhendong Li

Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone’s time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone’s time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.

arXiv:1909.00307•The Journal of Chemical Physics•Dec 27, 2019•12 citations

Full Configuration Interaction Quantum Monte Carlo for Accurate $\textit{Ab Initio}$ Nuclear Structure Calculations

Rongzhe Hu, Furong Xu, B. S. Hu +1 more

We introduce novel full configuration interaction quantum Monte Carlo (FCIQMC) as an accurate many-body solver for $\textit{ab initio}$ nuclear structure calculations. This stochastic approach directly samples the exact wave function in the full configuration space, enabling high-fidelity treatment of high-order many-body correlations in strongly interacting nuclear systems. Using interactions from chiral effective field theory, we have computed ground-state energies and charge radii of $^4$He, $^8$Be, $^{12}$C and $^{16}$O with sub-percent-level many-body uncertainties. These results establish FCIQMC as a stochastic full-configuration-space solver capable of treating systems beyond the reach of the conventional no-core shell model, and as an accurate benchmark for truncated many-body expansion methods.

arXiv:2606.22341•Jun 21, 2026•0 citations

Efficient and stochastic multireference perturbation theory for large active spaces within a full configuration interaction quantum Monte Carlo framework

Robert J. Anderson, Toru Shiozaki, George H. Booth

Full Configuration Interaction Quantum Monte Carlo (FCIQMC) has been effectively applied to very large configuration interaction (CI) problems and was recently adapted for use as an active space solver and combined with orbital optimization. In this work, we detail an approach within FCIQMC to allow for efficient sampling of fully internally contracted multireference perturbation theories within the same stochastic framework. Schemes are described to allow for the close control over the resolution of stochastic sampling of the effective higher-body intermediates within the active space. It is found that while complete active space second-order perturbation theory seems less amenable to a stochastic reformulation, strongly contracted N-Electron Valence second-order Perturbation Theory (NEVPT2) is far more stable, requiring a similar number of walkers to converge the sc-NEVPT2 expectation values as to converge the underlying CI problem. We demonstrate the application of the stochastic approach to the computation of sc-NEVPT2 within a (24, 24) active space in a biologically relevant system and show that small numbers of walkers are sufficient for a faithful sampling of the sc-NEVPT2 energy to chemical accuracy, despite the active space already exceeding the limits of practicality for traditional approaches. This raises prospects of an efficient stochastic solver for multireference chemical problems requiring large active spaces, with an accurate treatment of external orbitals.

arXiv:1911.11042•The Journal of Chemical Physics•Feb 3, 2020•41 citations

Multi-state effective Hamiltonian and size-consistency corrections in stochastic configuration interactions

Seiichiro L. Ten-no

Model space quantum Monte Carlo (MSQMC) is an extension of full configuration interaction QMC that allows us to calculate quasi-degenerate and excited electronic states by sampling the effective Hamiltonian in the model space. We introduce a novel algorithm based on the state-selective partitioning for the effective Hamiltonian using left eigenvectors to calculate several electronic states simultaneously at much less computational cost than the original MSQMC with the energy-dependent partitioning. The sampling of walkers in MSQMC is analyzed in the single reference limit using a stochastic algorithm for higher-order perturbation energies by the analogy of the deterministic case utilizing a full configuration interaction program. We further develop size-consistency corrections of the initiator adaptation (i-MSQMC) in three different ways, i.e., the coupled electron pair approximation, a posteriori, and second-order perturbative corrections. It is clearly demonstrated that most of the initiator error is originating from the deficiency of proper scaling of correlation energy due to its truncated CI nature of the initiator approximation and that the greater part of the error can be recovered by the size-consistency corrections developed in this work.

arXiv:1708.09498•The Journal of Chemical Physics•Dec 26, 2017•29 citations

Diagrammatic Monte Carlo for electronic correlation in molecules: High-order many-body perturbation theory with low scaling

Giacomo Bighin, Quoc P. Ho, M. Lemeshko +1 more

We present a low-scaling diagrammatic Monte Carlo approach to molecular correlation energies. Using combinatorial graph theory to encode many-body Hugenholtz diagrams, we sample the M\o{}ller-Plesset ($\mathrm{MP}n$) perturbation series, obtaining accurate correlation energies up to $n=5$, with quadratic scaling in the number of basis functions. Our technique reduces the computational complexity of the molecular many-fermion correlation problem, opening up the possibility of low-scaling, accurate stochastic computations for a wide class of many-body systems described by Hugenholtz diagrams.

arXiv:2203.12666•Physical review. B./Physical review. B•Jul 11, 2023•1 citation

Full configuration interaction quantum Monte Carlo for accurate $\textit{ab initio}$ nuclear structure calculations: algorithms and calculation details

Rongzhe Hu, Furong Xu, Baishan Hu +1 more

Full configuration interaction quantum Monte Carlo (FCIQMC) is a stochastic many-body solver that has been widely applied to electronic, molecular, and condensed-matter systems. In this work we apply FCIQMC to $\textit{ab initio}$ nuclear structure calculations using interactions derived from chiral effective field theory. We describe the algorithm in detail, including imaginary-time propagation, excitation generation, estimator choices, the initiator approximation with adaptive shift correction, and reduced-density-matrix (RDM) sampling. Benchmark calculations in small model spaces, where deterministic full configuration interaction (FCI) results are available, validate the stochastic calculation of energies, radii, and RDM-based pure estimators. For large model spaces, we analyze the residual finite-walker bias through systematic walker-number convergence and infinite-walker extrapolations. We also demonstrate that FCIQMC can be extended beyond ground-state calculations by computing the low-lying spectrum of $^6$Li.

arXiv:2607.05525•Jul 6, 2026•0 citations
High-order perturbative ab initio\textit{ab initio}ab initio calculations are challenging due to the rapidly growing configuration space and the difficulty of assessing convergence. In this letter, we introduce perturbation theory quantum Monte Carlo (PTQMC), a stochastic approach designed to compute high-order many-body perturbative corrections. By representing the perturbative wave function with random walkers in configuration space, PTQMC avoids the exponential scaling inherent to conventional constructions of high-rank excitation operators. Benchmark calculations for the Richardson pairing model demonstrate that PTQMC accurately reproduces exact many-body perturbation theory (MBPT) coefficients up to 16th order, even in strongly divergent regimes. We further show that combining PTQMC with series resummation techniques yields stable and precise energy estimates in cases where the straightforward perturbative series fails. Finally, we propose the effective number of configurations, eSe^{S}eS, as a global measure of perturbative wave-function complexity that can be directly extracted within PTQMC. We demonstrate that the saturation behavior of eSe^{S}eS provides a more reliable indicator of the validity of perturbative expansions than energy convergence alone.
Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone’s time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone’s time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.
Full Configuration Interaction Quantum Monte Carlo for Accurate Ab Initio\textit{Ab Initio}Ab Initio Nuclear Structure Calculations
We introduce novel full configuration interaction quantum Monte Carlo (FCIQMC) as an accurate many-body solver for ab initio\textit{ab initio}ab initio nuclear structure calculations. This stochastic approach directly samples the exact wave function in the full configuration space, enabling high-fidelity treatment of high-order many-body correlations in strongly interacting nuclear systems. Using interactions from chiral effective field theory, we have computed ground-state energies and charge radii of 4^44He, 8^88Be, 12^{12}12C and 16^{16}16O with sub-percent-level many-body uncertainties. These results establish FCIQMC as a stochastic full-configuration-space solver capable of treating systems beyond the reach of the conventional no-core shell model, and as an accurate benchmark for truncated many-body expansion methods.
Full Configuration Interaction Quantum Monte Carlo (FCIQMC) has been effectively applied to very large configuration interaction (CI) problems and was recently adapted for use as an active space solver and combined with orbital optimization. In this work, we detail an approach within FCIQMC to allow for efficient sampling of fully internally contracted multireference perturbation theories within the same stochastic framework. Schemes are described to allow for the close control over the resolution of stochastic sampling of the effective higher-body intermediates within the active space. It is found that while complete active space second-order perturbation theory seems less amenable to a stochastic reformulation, strongly contracted N-Electron Valence second-order Perturbation Theory (NEVPT2) is far more stable, requiring a similar number of walkers to converge the sc-NEVPT2 expectation values as to converge the underlying CI problem. We demonstrate the application of the stochastic approach to the computation of sc-NEVPT2 within a (24, 24) active space in a biologically relevant system and show that small numbers of walkers are sufficient for a faithful sampling of the sc-NEVPT2 energy to chemical accuracy, despite the active space already exceeding the limits of practicality for traditional approaches. This raises prospects of an efficient stochastic solver for multireference chemical problems requiring large active spaces, with an accurate treatment of external orbitals.
Model space quantum Monte Carlo (MSQMC) is an extension of full configuration interaction QMC that allows us to calculate quasi-degenerate and excited electronic states by sampling the effective Hamiltonian in the model space. We introduce a novel algorithm based on the state-selective partitioning for the effective Hamiltonian using left eigenvectors to calculate several electronic states simultaneously at much less computational cost than the original MSQMC with the energy-dependent partitioning. The sampling of walkers in MSQMC is analyzed in the single reference limit using a stochastic algorithm for higher-order perturbation energies by the analogy of the deterministic case utilizing a full configuration interaction program. We further develop size-consistency corrections of the initiator adaptation (i-MSQMC) in three different ways, i.e., the coupled electron pair approximation, a posteriori, and second-order perturbative corrections. It is clearly demonstrated that most of the initiator error is originating from the deficiency of proper scaling of correlation energy due to its truncated CI nature of the initiator approximation and that the greater part of the error can be recovered by the size-consistency corrections developed in this work.
We present a low-scaling diagrammatic Monte Carlo approach to molecular correlation energies. Using combinatorial graph theory to encode many-body Hugenholtz diagrams, we sample the M\o{}ller-Plesset (MPn\mathrm{MP}nMPn) perturbation series, obtaining accurate correlation energies up to n=5n=5n=5, with quadratic scaling in the number of basis functions. Our technique reduces the computational complexity of the molecular many-fermion correlation problem, opening up the possibility of low-scaling, accurate stochastic computations for a wide class of many-body systems described by Hugenholtz diagrams.
Full configuration interaction quantum Monte Carlo for accurate ab initio\textit{ab initio}ab initio nuclear structure calculations: algorithms and calculation details
Full configuration interaction quantum Monte Carlo (FCIQMC) is a stochastic many-body solver that has been widely applied to electronic, molecular, and condensed-matter systems. In this work we apply FCIQMC to ab initio\textit{ab initio}ab initio nuclear structure calculations using interactions derived from chiral effective field theory. We describe the algorithm in detail, including imaginary-time propagation, excitation generation, estimator choices, the initiator approximation with adaptive shift correction, and reduced-density-matrix (RDM) sampling. Benchmark calculations in small model spaces, where deterministic full configuration interaction (FCI) results are available, validate the stochastic calculation of energies, radii, and RDM-based pure estimators. For large model spaces, we analyze the residual finite-walker bias through systematic walker-number convergence and infinite-walker extrapolations. We also demonstrate that FCIQMC can be extended beyond ground-state calculations by computing the low-lying spectrum of 6^66Li.