pith. sign in

arxiv: 2602.08265 · v3 · submitted 2026-02-09 · ⚛️ nucl-th

Stochastic many-body perturbation theory for high-order calculations

Xin Zhen , Rongzhe Hu , Junchen Pei , Furong Xu This is my paper

Pith reviewed 2026-05-16 06:23 UTC · model grok-4.3

classification ⚛️ nucl-th
keywords many-body perturbation theoryquantum Monte Carlostochastic samplingRichardson pairing modelperturbative series resummationwave-function complexityhigh-order corrections
0
0 comments X

The pith

Stochastic walkers compute exact many-body perturbation coefficients to 16th order.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

High-order many-body perturbation theory faces exponential growth in the number of configurations, which blocks direct calculation of corrections beyond low orders. PTQMC represents the perturbative wave function by an ensemble of random walkers that move through configuration space, sidestepping the need to construct explicit high-rank excitation operators. Benchmarks on the Richardson pairing model show that the sampled coefficients match exact MBPT results through order 16, including regimes where the series diverges. The method also extracts the effective number of configurations e^S directly from the walker distribution; its saturation signals the breakdown of perturbation theory more reliably than monitoring energy alone. When the sampled series is resummed, stable and accurate energies are recovered where plain perturbation fails.

Core claim

PTQMC reproduces exact many-body perturbation theory coefficients up to sixteenth order by sampling the perturbative wave function with random walkers in configuration space, as demonstrated in the Richardson pairing model even in divergent regimes, while the extracted e^S provides a global diagnostic whose saturation indicates when the perturbative expansion remains valid.

What carries the argument

PTQMC, the representation of the perturbative wave function by an ensemble of random walkers that propagate in many-body configuration space to evaluate high-order corrections without explicit operator construction.

If this is right

  • High-order corrections become feasible for systems whose configuration spaces are too large for conventional determinant expansions.
  • Resummation applied to the PTQMC series produces stable energies in cases where the unresummed series diverges.
  • The e^S diagnostic extracted from walker statistics supplies a practical criterion for choosing the maximum perturbative order in production calculations.
  • Nuclear many-body calculations can incorporate higher-order perturbative effects without exponential cost growth.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach could be applied to realistic nuclear Hamiltonians to determine the perturbative order at which convergence fails for medium-mass nuclei.
  • Integration with existing quantum Monte Carlo codes would allow stochastic evaluation of response functions or transition amplitudes at high orders.
  • Variance analysis of the walker population on larger pairing models would reveal the practical scaling limits before bias appears.

Load-bearing premise

The random-walker sampling must faithfully represent the high-order perturbative wave function without uncontrolled statistical bias or variance that would distort the extracted coefficients.

What would settle it

A systematic deviation between PTQMC coefficients and independently computed exact MBPT coefficients for the Richardson model at order 16 or higher would show that the sampling introduces bias.

Figures

Figures reproduced from arXiv: 2602.08265 by Furong Xu, Junchen Pei, Rongzhe Hu, Xin Zhen.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: shows the effective number of configurations, e S, as a function of the coupling strength g for different perturbative orders. The dark-shaded regions indicate coupling intervals where e S exhibits poor convergence up to tenth order, quantified by the criterion [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper introduces perturbation theory quantum Monte Carlo (PTQMC), a stochastic method representing the perturbative wave function via random walkers in configuration space to compute high-order many-body perturbation theory (MBPT) corrections without exponential scaling. Benchmark calculations on the Richardson pairing model are reported to reproduce exact MBPT coefficients up to 16th order even in divergent regimes; the method is combined with series resummation for stable energies, and the effective number of configurations e^S is proposed as a diagnostic for perturbative wave-function complexity whose saturation indicates validity better than energy convergence alone.

Significance. If the stochastic representation is unbiased and variance remains controllable, PTQMC could enable routine high-order MBPT in nuclear many-body systems where configuration spaces grow rapidly and series diverge, providing a practical route to precise energies via resummation. The e^S diagnostic offers a potentially useful global complexity measure extractable directly from the sampling.

major comments (2)
  1. [Abstract and benchmark calculations] Abstract and benchmark calculations: The headline claim that PTQMC reproduces exact MBPT coefficients up to order 16 requires that the random-walker sampling yields unbiased estimators with controllable variance. No information is supplied on walker population size, time-step errors, or variance-control techniques, which are load-bearing for validating the numerical agreement in strongly divergent regimes.
  2. [Section introducing e^S] Section introducing e^S: The effective number of configurations e^S is presented as a directly extractable diagnostic whose saturation reliably signals the breakdown of perturbation theory. Its precise definition from the walker distribution, including any normalization or averaging procedure, must be specified to confirm it is independent of sampling bias.
minor comments (2)
  1. [Method section] Notation for the perturbative wave function and walker representation should be introduced with explicit equations early in the manuscript to clarify how the stochastic sampling maps onto the standard MBPT expansion.
  2. [Figures] Figure captions for the Richardson-model benchmarks should include the specific walker parameters and error bars used in each panel to allow direct assessment of statistical precision.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have revised the manuscript to supply the requested technical details.

read point-by-point responses

  1. Referee: [Abstract and benchmark calculations] Abstract and benchmark calculations: The headline claim that PTQMC reproduces exact MBPT coefficients up to order 16 requires that the random-walker sampling yields unbiased estimators with controllable variance. No information is supplied on walker population size, time-step errors, or variance-control techniques, which are load-bearing for validating the numerical agreement in strongly divergent regimes.

    Authors: We agree that these implementation parameters are necessary to substantiate the benchmark claims. In the revised manuscript we have added a dedicated paragraph in the Methods section specifying a walker population of 10^6, imaginary-time step Δτ = 0.005, and the use of branching with population control. The PTQMC propagator is constructed to satisfy detailed balance for the perturbative expansion, guaranteeing unbiased estimators in the infinite-walker limit; residual time-step and population-control errors are shown to lie below the reported statistical uncertainties (typically < 0.5 % relative error) up to 16th order. Multiple independent runs with error bars are now included in the revised figures to demonstrate variance control even in divergent regimes. revision: yes

  2. Referee: [Section introducing e^S] Section introducing e^S: The effective number of configurations e^S is presented as a directly extractable diagnostic whose saturation reliably signals the breakdown of perturbation theory. Its precise definition from the walker distribution, including any normalization or averaging procedure, must be specified to confirm it is independent of sampling bias.

    Authors: We appreciate the request for an explicit definition. In the revised manuscript we now state that e^S is obtained from the equilibrated walker weights via the Shannon entropy S = −∑_c p_c ln p_c, where p_c = w_c / ∑ w_j is the normalized probability of configuration c. The sum is performed over all visited configurations and S is averaged over the final 20 % of the Monte Carlo steps after equilibration. Because the underlying stochastic process samples the exact perturbative weights (unbiased by construction), the resulting e^S inherits the same lack of bias; we have added a short derivation and a convergence check with respect to walker number to confirm this independence. revision: yes

Circularity Check

0 steps flagged

No significant circularity; method validated against independent exact benchmarks

full rationale

The paper introduces PTQMC as a stochastic walker representation of the perturbative wave function and validates it by direct numerical reproduction of independently known exact MBPT coefficients for the Richardson pairing model up to order 16. No load-bearing step reduces the reported coefficients, the e^S diagnostic, or the resummation results to a fit performed on the same data or to a self-citation chain. The benchmark comparison is external and falsifiable; the sampling construction itself does not presuppose the target coefficients. This is the normal case of a self-contained computational method with external verification.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claims rest on standard quantum-mechanical assumptions about the many-body Hilbert space and the validity of the perturbative expansion itself; no new free parameters or invented entities are introduced beyond the stochastic representation and the e^S definition.

axioms (2)
  • domain assumption The many-body Hilbert space admits a perturbative expansion in powers of the interaction Hamiltonian.
    Invoked implicitly when the paper states that PTQMC computes MBPT coefficients.
  • domain assumption Random-walker sampling in configuration space can be made unbiased for the perturbative wave-function components.
    Required for the claim that PTQMC reproduces exact coefficients.
invented entities (1)
  • effective number of configurations e^S no independent evidence
    purpose: Global scalar measure of perturbative wave-function complexity extracted from the walker population.
    New diagnostic proposed in the paper; no independent evidence outside the PTQMC runs is supplied.

pith-pipeline@v0.9.0 · 5486 in / 1386 out tokens · 48606 ms · 2026-05-16T06:23:37.180219+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages · 1 internal anchor

  1. [1]

    The Padé approximation represents a perturbative series by a ratio of two polynomials and is known to cap- ture nonanalytic structures such as singularities

    As a simple and widely used resummation scheme, the Padé approximation was employed [ 41, 42], with sta- tistical uncertainties propagated from the PTQMC input data. The Padé approximation represents a perturbative series by a ratio of two polynomials and is known to cap- ture nonanalytic structures such as singularities. It has therefore been used to ext...

    work page

  2. [2]

    Machleidt and F

    R. Machleidt and F. Sammarruca, Recent advances in chiral EFT based nuclear forces and their applications , Prog. Part. Nucl. Phys. 137, 104117 (2024)

    work page 2024

  3. [3]

    Hergert, A guided tour of ab initio nuclear many-body theory, Front

    H. Hergert, A guided tour of ab initio nuclear many-body theory, Front. Phys. 8, 379 (2020)

    work page 2020

  4. [4]

    Ekström, et al

    A. Ekström, et al. , What is ab initio in nuclear theory? , Front. Phys. 11, 1129094 (2023)

    work page 2023

  5. [5]

    Gandolfi, D

    S. Gandolfi, D. Lonardoni, A. Lovato, and M. Piarulli, Atomic nuclei from quantum Monte Carlo calculations with chiral EFT interactions , Front. Phys. 8, 117 (2020)

    work page 2020

  6. [6]

    Drischler, Nuclear matter from chiral effective field theory, Ph.D

    C. Drischler, Nuclear matter from chiral effective field theory, Ph.D. thesis, Technische Univ. Darmstadt (Ger- many) (2017)

    work page 2017

  7. [7]

    Drischler, K

    C. Drischler, K. Hebeler, and A. Schwenk, Chiral inter- actions up to next-to-next-to-next-to-leading order and nuclear saturation, Phys. Rev. Lett. 122, 042501 (2019)

    work page 2019

  8. [8]

    Drischler, J

    C. Drischler, J. W. Holt, and C. Wellenhofer, Chiral ef- fective field theory and the high-density nuclear equation of state , Annu. Rev. Nucl. Part. Sci. 71, 403 (2021)

    work page 2021

  9. [9]

    Drischler, R

    C. Drischler, R. J. Furnstahl, J. A. Melendez, and D. R. Phillips, How well do we know the neutron-matter equa- tion of state at the densities inside neutron stars? A Bayesian approach with correlated uncertainties , Phys. Rev. Lett. 125, 202702 (2020) . 7

    work page 2020

  10. [10]

    Drischler, A

    C. Drischler, A. Carbone, K. Hebeler, and A. Schwenk, Neutron matter from chiral two- and three-nucleon cal- culations up to N 3LO, Phys. Rev. C 94, 054307 (2016)

    work page 2016

  11. [11]

    Tichai, R

    A. Tichai, R. Roth, and T. Duguet, Many-body perturba- tion theories for finite nuclei , Front. Phys. 8, 164 (2020)

    work page 2020

  12. [12]

    W. G. Jiang, C. Forssén, T. Djärv, and G. Hagen, Nuclear-matter saturation and symmetry energy within ∆-full chiral effective field theory , Phys. Rev. C 109, L061302 (2024)

    work page 2024

  13. [13]

    Baardsen, A

    G. Baardsen, A. Ekström, G. Hagen, and M. Hjorth- Jensen, Coupled-cluster studies of infinite nuclear mat- ter, Phys. Rev. C 88, 054312 (2013)

    work page 2013

  14. [14]

    Hagen, et al

    G. Hagen, et al. , Coupled-cluster calculations of nucle- onic matter , Phys. Rev. C 89, 014319 (2014)

    work page 2014

  15. [15]

    Xu, et al

    X.-Y. Xu, et al. , Progress in ab initio in-medium similar- ity renormalization group and coupled-channel method with coupling to the continuum , Nucl. Sci. Tech. 35, 215 (2024)

    work page 2024

  16. [16]

    Zhen, et al

    X. Zhen, et al. , Non-perturbative calculations of nu- clear matter using in-medium similarity renormalization group, Phys. Lett. B 862, 139350 (2025)

    work page 2025

  17. [18]

    J. G. Lietz, S. Novario, G. R. Jansen, G. Hagen, and M. Hjorth-Jensen, Computational Nuclear Physics and Post Hartree-Fock Methods , in An Advanced Course in Computational Nuclear Physics: Bridging the Scales from Quarks to Neutron Stars , Lecture Notes in Physics, Vol. 936, edited by M. Hjorth-Jensen, M. P. Lombardo, and U. van Kolck (Springer, Cham, 201...

    work page 2017

  18. [20]

    Tichai, J

    A. Tichai, J. Langhammer, S. Binder, and R. Roth, Hartree–Fock many-body perturbation theory for nuclear ground-states, Phys. Lett. B 756, 283 (2016)

    work page 2016

  19. [21]

    R. J. Bartlett and M. Musiał, Coupled-cluster theory in quantum chemistry, Rev. Mod. Phys. 79, 291 (2007)

    work page 2007

  20. [22]

    Hagen, T

    G. Hagen, T. Papenbrock, M. Hjorth-Jensen, and D. J. Dean, Coupled-cluster computations of atomic nuclei , Rep. Prog. Phys. 77, 096302 (2014)

    work page 2014

  21. [23]

    Hergert, S

    H. Hergert, S. K. Bogner, T. D. Morris, A. Schwenk, and K. Tsukiyama, The in-medium similarity renormalization group: A novel ab initio method for nuclei , Phys. Rep. 621, 165 (2016)

    work page 2016

  22. [24]

    Li, Advances in many-body perturbation theory for closed-and open-shell nuclei , Ph.D

    Z. Li, Advances in many-body perturbation theory for closed-and open-shell nuclei , Ph.D. thesis, Université de Bordeaux (2023)

    work page 2023

  23. [25]

    M. L. Leininger, W. D. Allen, H. F. Schaefer, III, and C. D. Sherrill, Is Møller–Plesset perturbation theory a convergent ab initio method?, J. Chem. Phys. 112, 9213 (2000)

    work page 2000

  24. [26]

    Kvaal, E

    S. Kvaal, E. Jarlebring, and W. Michiels, Computing sin- gularities of perturbation series , Phys. Rev. A 83, 032505 (2011)

    work page 2011

  25. [27]

    Brolli, C

    S. Brolli, C. Barbieri, and E. Vigezzi, Diagrammatic Monte Carlo for Finite Systems at Zero Temperature , Phys. Rev. Lett. 134, 182502 (2025)

    work page 2025

  26. [28]

    Faribault, O

    A. Faribault, O. El Araby, C. Sträter, and V. Gritsev, Gaudin model solver based on the Bethe ansatz/ordinary differential equations correspondence , Phys. Rev. B 83, 235124 (2011)

    work page 2011

  27. [29]

    Snyman and H

    I. Snyman and H. B. Geyer, Strong-coupling limit of the Richardson Hamiltonian analyzed using Dyson mapping , Phys. Rev. B 73, 144516 (2006)

    work page 2006

  28. [30]

    P. W. Claeys, S. De Baerdemacker, M. Van Raemdonck, and D. Van Neck, Eigenvalue-based method and form- factor determinant representations for integrable XXZ Richardson-Gaudin models , Phys. Rev. B 91, 155102 (2015)

    work page 2015

  29. [31]

    G. H. Booth, A. J. W. Thom, and A. Alavi, Fermion Monte Carlo without fixed nodes: A game of life, death, and annihilation in Slater determinant space , J. Chem. Phys. 131, 054106 (2009)

    work page 2009

  30. [32]

    J. J. Shepherd, G. Booth, A. Grüneis, and A. Alavi, Full configuration interaction perspective on the homo- geneous electron gas , Phys. Rev. B 85, 081103 (2012)

    work page 2012

  31. [33]

    G. H. Booth, A. Grüneis, G. Kresse, and A. Alavi, To- wards an exact description of electronic wavefunctions in real solids , Nature 493, 365 (2013)

    work page 2013

  32. [34]

    Lu, et al

    B.-N. Lu, et al. , Perturbative Quantum Monte Carlo Method for Nuclear Physics , Phys. Rev. Lett. 128, 242501 (2022)

    work page 2022

  33. [35]

    Curry, J

    R. Curry, J. E. Lynn, K. E. Schmidt, and A. Gezerlis, Second-order perturbation theory in continuum quantum Monte Carlo calculations , Phys. Rev. Res. 5, L042021 (2023)

    work page 2023

  34. [36]

    Wlazłowski, J

    G. Wlazłowski, J. W. Holt, S. Moroz, A. Bulgac, and K. J. Roche, Auxiliary-field quantum Monte Carlo simu- lations of neutron matter in chiral effective field theory , Phys. Rev. Lett. 113, 182503 (2014)

    work page 2014

  35. [37]

    J. E. Lynn, et al. , Chiral three-nucleon interactions in light nuclei, neutron- α scattering, and neutron matter , Phys. Rev. Lett. 116, 062501 (2016)

    work page 2016

  36. [38]

    Roggero, A

    A. Roggero, A. Mukherjee, and F. Pederiva, Quantum Monte Carlo calculations of neutron matter with nonlocal chiral interactions, Phys. Rev. Lett. 112, 221103 (2014)

    work page 2014

  37. [39]

    Arthuis, C

    P. Arthuis, C. Barbieri, F. Pederiva, and A. Rog- gero, Quantum Monte Carlo calculations in configura- tion space with three-nucleon forces , Phys. Rev. C 107, 044303 (2023)

    work page 2023

  38. [40]

    Carlson, et al

    J. Carlson, et al. , Quantum Monte Carlo methods for nuclear physics, Rev. Mod. Phys. 87, 1067 (2015)

    work page 2015

  39. [41]

    [ ? ? ? ]

    See Supplemental Material for a detailed description of the PTQMC algorithm, which includes Refs. [ ? ? ? ]

    work page

  40. [42]

    Padé, Sur la représentation approchée d’une fonction par des fractions rationnelles , Annales scientifiques de l’École Normale Supérieure Série 3, 9, 3 (1892)

    H. Padé, Sur la représentation approchée d’une fonction par des fractions rationnelles , Annales scientifiques de l’École Normale Supérieure Série 3, 9, 3 (1892)

    work page

  41. [43]

    Bender and S

    C. Bender and S. Orszag, Advanced Mathematical Meth- ods for Scientists and Engineers: Asymptotic Methods and Perturbation Theory , Vol. 1 (Springer, 1999)

    work page 1999

  42. [44]

    Roth and J

    R. Roth and J. Langhammer, Padé-resummed high-order perturbation theory for nuclear structure calculations , Phys. Lett. B 683, 272 (2010)

    work page 2010

  43. [45]

    Nuclear and neutron matter equations of state from high-quality potentials up to fifth order of the chiral expansion

    F. Sammarruca, et al. , Nuclear and neutron mat- ter equations of state from high-quality potentials up to fifth order of the chiral expansion , arXiv preprint arXiv:1807.06640 (2018)

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  44. [46]

    Langhammer, R

    J. Langhammer, R. Roth, and C. Stumpf, Spectra of open-shell nuclei with Padé-resummed degenerate per- turbation theory , Phys. Rev. C 86, 054315 (2012)

    work page 2012

  45. [47]

    Heinz, Analysis of Three-Body Effects in the In- Medium Similarity Renormalization Group, Master’s the- sis, Technische Univ

    M. Heinz, Analysis of Three-Body Effects in the In- Medium Similarity Renormalization Group, Master’s the- sis, Technische Univ. Darmstadt (Germany) (2020)

    work page 2020

  46. [48]

    Marino, W

    F. Marino, W. G. Jiang, and S. J. Novario, Diagrammatic ab initio methods for infinite nuclear matter with modern chiral interactions, Phys. Rev. C 110, 054322 (2024) . 8

    work page 2024

  47. [49]

    R. J. Bartlett and J. Noga, The expectation value coupled-cluster method and analytical energy deriva- tives, Chem. Phys. Lett. 150, 29 (1988)

    work page 1988

  48. [50]

    C. E. Shannon, A mathematical theory of communica- tion, Bell Syst. Tech. J. 27, 379 (1948)

    work page 1948

  49. [51]

    V. V. Sokolov, B. A. Brown, and V. Zelevinsky, Invariant correlational entropy and complexity of quantum states , Phys. Rev. E 58, 56 (1998)

    work page 1998

  50. [52]

    V. V. Flambaum and F. M. Izrailev, Entropy production and wave packet dynamics in the Fock space of closed chaotic many-body systems , Phys. Rev. E 64, 036220 (2001)

    work page 2001

  51. [53]

    J. R. Armstrong, S. Åberg, S. M. Reimann, and V. G. Zelevinsky, Complexity of quantum states in the two- dimensional pairing model , Phys. Rev. E 86, 066204 (2012)

    work page 2012

  52. [54]

    V. K. B. Kota and R. Sahu, Information entropy and number of principal components in shell model transition strength distributions , Phys. Lett. B 429, 1 (1998)

    work page 1998

  53. [55]

    X. Zhen, R. Hu, J. Pei, and F. Xu, PTQMC data repos- itory, https://github.com/AlbusZhen/PTQMC_Pairing (2026). Supplemental Material for: Stochastic many-body perturbation theory for high-order calculations X. Zhen ,1 R. Z. Hu ,1 J. C. Pei ,1, ∗ and F. R. Xu 1, 2, † 1School of Physics, and State Key Laboratory of Nuclear Physics and Technology, Peking Uni...

    work page 2026

  54. [56]

    We initialize the walker population to represent the first-order amplitudes w(1) (which enter the second-order energy estimator E(2) in Eq

    Initialization. We initialize the walker population to represent the first-order amplitudes w(1) (which enter the second-order energy estimator E(2) in Eq. (8) of the main text). To control the population size, we monitor the total walker number Nw = X I ˆNI , (1) and normalize it to a target value N0. Because both the walker population and the number of ...

    work page

  55. [57]

    For each walker on |ΦJ ⟩, we randomly choose a connected determinant |ΦI ⟩ with proposal probability pgen(I|J) using the excitation-generation algorithm [ 6]

    Spawning. For each walker on |ΦJ ⟩, we randomly choose a connected determinant |ΦI ⟩ with proposal probability pgen(I|J) using the excitation-generation algorithm [ 6]. A spawn from J to I is attempted with probability ps(I|J) = 1 ∆I0 |HIJ | pgen(I|J) , (2) and the phase of the spawned walker on |ΦI ⟩ is given by −ˆhIJ ˆnJ , where ˆhIJ = HIJ /|HIJ |. The ...

    work page

  56. [58]

    Annihilation and evaluation. After spawning, all contributions on the same determinant |ΦI ⟩ are combined (annihilated) by summing their phases, yielding the updated complex population ˆNI = X a∈ΦI ˆna = NI ˆnI . (3) Once the redistribution is completed, the MBPT correlation energy is evaluated using Eq. (10) of the main text. A normalization factor has t...

    work page

  57. [59]

    Shavitt and R

    I. Shavitt and R. J. Bartlett, Many-body methods in chemistry and physics: MBPT and coupled-cluster theory (Cambridge University Press, 2009)

    work page 2009

  58. [60]

    Tichai, R

    A. Tichai, R. Roth, and T. Duguet, Many-body perturbation theories for finite nuclei , Front. Phys. 8, 164 (2020)

    work page 2020

  59. [61]

    J. C. Slater, The Theory of Complex Spectra , Phys. Rev. 34, 1293 (1929)

    work page 1929

  60. [62]

    E. U. Condon, The Theory of Complex Spectra , Phys. Rev. 36, 1121 (1930)

    work page 1930

  61. [63]

    Arthuis, T

    P. Arthuis, T. Duguet, A. Tichai, R.-D. Lasseri, and J.-P. Ebran, ADG: Automated generation and evaluation of many-body diagrams I. Bogoliubov many-body perturbation theory , Comput. Phys. Commun. 240, 202 (2019)

    work page 2019

  62. [64]

    Hu, et al

    R. Hu, et al. , Ab initio exact calculation of strongly correlated nucleonic matter (2026), arXiv:2508.09252 [nucl-th]

    work page arXiv 2026