pith. sign in

arxiv: 2606.02549 · v1 · pith:JX4FVEFMnew · submitted 2026-06-01 · ⚛️ physics.chem-ph · physics.atom-ph· physics.comp-ph

Diagrammatic Monte Carlo for positron-molecule many-body theory

T. A. Scott , S. K. Gregg , D. G. Green This is my paper

Pith reviewed 2026-06-28 11:49 UTC · model grok-4.3

classification ⚛️ physics.chem-ph physics.atom-phphysics.comp-ph
keywords diagrammatic Monte Carlopositron self-energyladder seriesTamm-Dancoff approximationmany-body perturbation theorydensity fittingvirtual positroniummolecular calculations
0
0 comments X

The pith

Diagrammatic Monte Carlo samples infinite ladder series for positron self-energy with memory reduced by the number of orbitals.

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

This paper introduces a Monte Carlo method to evaluate the ladder diagram contributions to the self-energy of a positron in a molecular field by sampling them stochastically order by order. The approach uses the Tamm-Dancoff approximation, which is exact for the virtual-positronium and positron-hole series, and applies Cesàro-Riesz resummation to reach infinite order from finite samples. It employs Gaussian bases and density fitting so that only three-center integrals need storage, cutting memory use by a factor of roughly the number of molecular orbitals compared to solving the full Bethe-Salpeter equations deterministically. Benchmarks on lithium hydride reproduce exact diagonalization results, confirming the method can sum the virtual-positronium series accurately.

Core claim

The stochastic diagrammatic Monte Carlo evaluation samples the GW@TDHF, virtual-positronium T-matrix, and positron-hole Goldstone ladder series to the positron correlation potential within the Tamm-Dancoff approximation and extrapolates to infinite order via Cesàro-Riesz resummation, achieving quantitative agreement with exact diagonalisation for lithium hydride while reducing the memory footprint of the largest arrays by a factor on the order of the number of molecular orbitals N ~ 10^2-10^3.

What carries the argument

Stochastic order-by-order sampling of ladder diagrams combined with Cesàro-Riesz resummation and density fitting of Coulomb integrals.

If this is right

  • Calculations become feasible for molecules requiring basis sets too large for full Bethe-Salpeter storage.
  • The successful summation of the infinite virtual-positronium ladder series opens accurate treatment of electron-positron correlations.
  • Gaussian basis sets with density fitting allow seamless integration into existing molecular codes.
  • Memory scaling improves from quadratic or worse in N to linear or better in practice for the sampled quantities.

Where Pith is reading between the lines

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

  • This technique may generalize to sampling other diagram classes in molecular many-body calculations.
  • Reduced memory could allow studies of positron binding in larger biomolecules or clusters.
  • Hybrid stochastic-deterministic schemes might further optimize different parts of the self-energy.
  • The method's accuracy on LiH suggests it could benchmark against experimental positron annihilation rates in molecules.

Load-bearing premise

The Cesàro-Riesz resummation procedure extrapolates the stochastically sampled finite-order results to infinite order without introducing uncontrolled bias.

What would settle it

Disagreement between the extrapolated Monte Carlo results and exact diagonalization benchmarks for the lithium hydride virtual-positronium contribution would indicate either sampling errors or bias in the resummation.

Figures

Figures reproduced from arXiv: 2606.02549 by D. G. Green, S. K. Gregg, T. A. Scott.

Figure 1
Figure 1. Figure 1: FIG. 1. The three infinite classes of diagram contributions to [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Positron binding energy of LiH at the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Positron binding energy of LiH at the combined [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
read the original abstract

A diagrammatic Monte Carlo evaluation of the ladder series contributions to the correlation potential (self energy) of a positron in the field of a molecule is presented. The $GW$@TDHF, virtual-positronium ($T$-matrix), and positron-hole Goldstone ladder series contributions are stochastically sampled order-by-order within the Tamm-Dancoff approximation, which is exact for the latter two classes, with Ces{\'a}ro-Riesz resummation used to extrapolate to infinite order. Gaussian bases are employed and Coulomb matrix elements are represented via density fitting, with the three centre integrals the largest arrays required to be stored in memory. The stochastic approach thus realizes a reduction in memory of the largest arrays required on the order of the number of molecular orbitals in the basis $N\sim$10$^2$--10$^3$ compared to the exact deterministic solution of Bethe-Salpeter equations [J. Hofierka, B. Cunningham, C. M. Rawlins, C. H. Patterson and D. G. Green, Nature {\bf 606}, {688} (2022)]. Benchmark results for lithium hydride show quantitative agreement with exact diagonalisation, notably demonstrating the successful stochastic summation of the virtual-positronium infinite electron-positron ladder series.

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

3 major / 1 minor

Summary. The paper presents a diagrammatic Monte Carlo method for stochastically sampling the GW@TDHF, virtual-positronium (T-matrix), and positron-hole Goldstone ladder series contributions to the positron self-energy in molecules. Sampling is performed order-by-order within the Tamm-Dancoff approximation (asserted exact for the latter two classes), with Cesàro-Riesz resummation to extrapolate to infinite order. Coulomb integrals use density fitting, yielding a memory reduction of order N (number of molecular orbitals, ~10^2-10^3) relative to deterministic Bethe-Salpeter equation solutions. Benchmarks on LiH show quantitative agreement with exact diagonalization, particularly for the virtual-positronium ladder series.

Significance. If the resummation procedure is free of uncontrolled bias, the approach would enable inclusion of infinite ladder diagrams in positron-molecule calculations at substantially reduced memory cost, facilitating larger basis sets than deterministic methods allow. The reported quantitative agreement with exact diagonalization on LiH provides supporting evidence for the stochastic sampling strategy on this system.

major comments (3)
  1. [Abstract] Abstract: The central claim of 'successful stochastic summation of the virtual-positronium infinite electron-positron ladder series' rests on Cesàro-Riesz extrapolation of finite-order stochastic samples. No discussion or test is provided of possible systematic bias in this resummation for sign-alternating or slowly convergent series, nor comparison against analytic resummations or higher-order exact results; this directly impacts the reliability of the infinite-order results and the claimed memory-reduction advantage over exact BSE.
  2. [Abstract] Abstract (benchmark paragraph): The statement of 'quantitative agreement with exact diagonalisation' for LiH is given without reported statistical uncertainties, sampling variances, number of Monte Carlo samples, maximum sampled order, or extrapolation parameters. This omission prevents assessment of whether the agreement is robust or could mask extrapolation errors that appear only at higher orders or for other molecules.
  3. [Abstract] Abstract: The assertion that the Tamm-Dancoff approximation 'is exact for the latter two classes' (virtual-positronium and positron-hole) is stated without derivation, proof sketch, or reference. Since this underpins the validity of restricting the stochastic sampling, a concise justification is needed to support the method's scope.
minor comments (1)
  1. [Abstract] The abstract could usefully specify the Gaussian basis sets and density-fitting details employed for the LiH benchmarks to allow reproducibility assessment.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and constructive comments. We address each major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim of 'successful stochastic summation of the virtual-positronium infinite electron-positron ladder series' rests on Cesàro-Riesz extrapolation of finite-order stochastic samples. No discussion or test is provided of possible systematic bias in this resummation for sign-alternating or slowly convergent series, nor comparison against analytic resummations or higher-order exact results; this directly impacts the reliability of the infinite-order results and the claimed memory-reduction advantage over exact BSE.

    Authors: We agree that the abstract would benefit from a concise note on the resummation. In revision we will add a sentence referencing the known convergence properties of Cesàro-Riesz summation for diagrammatic series and noting that the LiH benchmark (quantitative agreement within statistical error) provides empirical validation; a fuller discussion of bias tests appears in the main text. revision: yes

  2. Referee: [Abstract] Abstract (benchmark paragraph): The statement of 'quantitative agreement with exact diagonalisation' for LiH is given without reported statistical uncertainties, sampling variances, number of Monte Carlo samples, maximum sampled order, or extrapolation parameters. This omission prevents assessment of whether the agreement is robust or could mask extrapolation errors that appear only at higher orders or for other molecules.

    Authors: The abstract is space-limited, but we will revise the benchmark sentence to include or explicitly reference the statistical uncertainties, sample counts, maximum order, and extrapolation parameters that are already reported in the main text and supplementary material. revision: yes

  3. Referee: [Abstract] Abstract: The assertion that the Tamm-Dancoff approximation 'is exact for the latter two classes' (virtual-positronium and positron-hole) is stated without derivation, proof sketch, or reference. Since this underpins the validity of restricting the stochastic sampling, a concise justification is needed to support the method's scope.

    Authors: The exactness follows because the virtual-positronium T-matrix and positron-hole ladder diagrams contain no anti-resonant contributions that the TDA neglects; this is a standard property of these channels when the positron is treated as distinguishable. We will insert a brief parenthetical justification and a reference to the relevant diagrammatic analysis in the revised abstract. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained with external validation

full rationale

The paper applies standard diagrammatic Monte Carlo sampling to ladder series within the Tamm-Dancoff approximation (asserted exact for virtual-positronium and positron-hole classes) followed by Cesàro-Riesz extrapolation, with memory reduction relative to deterministic BSE. The benchmark agreement with exact diagonalization on LiH is an independent external check, not a self-referential fit or prediction. The cited Nature paper provides the deterministic reference but is not used to justify the stochastic method itself by construction. No quoted step reduces a claimed result to its own inputs or to a self-citation chain; the central claims rest on perturbation theory and numerical validation outside the present work.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach rests on standard many-body perturbation theory assumptions plus the stated exactness of Tamm-Dancoff for two series classes; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption Tamm-Dancoff approximation is exact for virtual-positronium and positron-hole ladder series
    Explicitly stated in the abstract as the basis for the stochastic sampling.
  • domain assumption Cesàro-Riesz resummation converges to the correct infinite-order result from finite-order stochastic samples
    Used to extrapolate; no independent justification supplied in abstract.

pith-pipeline@v0.9.1-grok · 5775 in / 1378 out tokens · 24037 ms · 2026-06-28T11:49:36.238602+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

51 extracted references

  1. [1]

    sider the positron-holeΛladder series

    The black horizontal dashed line marks the reference EXCITON+Bethe-Salpeter equation solution via exact diago- nalisation. sider the positron-holeΛladder series. The convergence behaviour of this is qualitatively similar to that of the RPA@TDA and TDHF@TDA series. Cesàro–Riesz re- summation produces stable extrapolations to1/N→0 across the full rangeδ= 0–...

  2. [2]

    Hofierka, B

    J. Hofierka, B. Cunningham, C. M. Rawlins, C. H. Pat- terson, and D. G. Green, Many-body theory of positron bindingtopolyatomicmolecules,Nature606,688(2022)

    2022

  3. [3]

    Tuomisto and I

    F. Tuomisto and I. Makkonen, Defect identification in semiconductors with positron annihilation: Experiment and theory, Rev. Mod. Phys.85, 1583 (2013)

    2013

  4. [4]

    Hugenschmidt, Positrons in surface physics, Surf

    C. Hugenschmidt, Positrons in surface physics, Surf. Sci. Rep.71, 547 (2016)

    2016

  5. [5]

    J. R. Danielson, D. H. E. Dubin, R. G. Greaves, and C. M. Surko, Plasma and trap-based techniques for sci- ence with positrons, Rev. Mod. Phys.87, 247 (2015)

    2015

  6. [6]

    G.B.Saha,Basics of PET imaging in physics, chemistry, and regulations(Springer, New York, 2005)

    2005

  7. [7]

    R. L. Wahal,Principles and Practice of Positron Emis- sion Tomography(Lippincott, Williams and Wilkins, Philadelphia, 2008)

    2008

  8. [8]

    Moskal, J

    P. Moskal, J. Baran,et al., Positronium image of the human brain in vivo, Science Advances10, eadp2840 (2024)

    2024

  9. [9]

    Moskal, A

    P. Moskal, A. Bilewicz, M. Das, B. Huang, A. Khreptak, 6 S. Parzych, J. Qi, A. Rominger, R. Seifert, S. Sharma, K. Shi, W. M. Steinberger, R. Walczak, and E. Stępień, Positronium imaging: History, current status, and fu- ture perspectives, IEEE Transactions on Radiation and Plasma Medical Sciences9, 981 (2025)

    2025

  10. [10]

    R. J. Drachman, Why positron physics is fun, AIP Con- ference Proceedings360, 369 (1996)

    1996

  11. [11]

    Prantzos, C

    N. Prantzos, C. Boehm, A. M. Bykov, R. Diehl, K. Fer- rière, N. Guessoum, P. Jean, J. Knoedlseder, A. Mar- cowith, I. V. Moskalenko, A. Strong, and G. Weidens- pointner, The 511 keV emission from positron annihila- tion in the Galaxy, Rev. Mod. Phys.83, 1001 (2011), publisher: American Physical Society

    2011

  12. [12]

    G. M. Fuller, A. Kusenko, D. Radice, and V. Takhistov, Positrons and 511 kev radiation as tracers of recent bi- nary neutron star mergers, Phys. Rev. Lett.122, 121101 (2019)

    2019

  13. [13]

    V. V. Flambaum and I. B. Samsonov, Radiation from matter-antimatter annihilation in the quark nugget model of dark matter, Phys. Rev. D104, 063042 (2021)

    2021

  14. [14]

    G. F. Gribakin, J. A. Young, and C. M. Surko, Positron- molecule interactions: Resonant attachment, annihila- tion, and bound states, Rev. Mod. Phys.82, 2557 (2010)

    2010

  15. [15]

    G. F. Gribakin, J. F. Stanton, J. R. Danielson, M. R. Natisin, and C. M. Surko, Mode coupling and mul- tiquantum vibrational excitations in feshbach-resonant positron annihilation in molecules, Phys. Rev. A96, 062709 (2017)

    2017

  16. [16]

    M. Y. Amusia, N. A. Cherepkov, L. V. Chernysheva, and S. G. Shapiro, Elastic scattering of slow positrons by he- lium, J. Phys. B: Atom. Mol. Phys.9, L531 (1976)

    1976

  17. [17]

    V. A. Dzuba, V. V. Flambaum, W. A. King, B. N. Miller, and O. P. Sushkov, Interaction between slow positrons and atoms, Phys. Scr.T46, 248 (1993)

    1993

  18. [18]

    V. A. Dzuba, V. V. Flambaum, G. F. Gribakin, and W.A.King,Boundstatesofpositronsandneutralatoms, Phys. Rev. A52, 4541 (1995)

    1995

  19. [19]

    G. F. Gribakin and J. Ludlow, Many-body theory of positron-atom interactions, Phys. Rev. A70, 032720 (2004)

    2004

  20. [20]

    Harabati, V

    C. Harabati, V. Dzuba, and V. Flambaum, Identifica- tion of atoms that can bind positrons, Phys. Rev. A.89, 022517 (2014)

    2014

  21. [21]

    Müller and L

    M. Müller and L. S. Cederbaum, Many-body theory of composite electronic-positronic systems, Phys. Rev. A 42, 170 (1990)

    1990

  22. [22]

    D. G. Green and G. F. Gribakin, Positron scattering and annihilation in hydrogenlike ions, Phys. Rev. A88, 032708 (2013)

    2013

  23. [23]

    D. G. Green, J. A. Ludlow, and G. F. Gribakin, Positron scattering and annihilation on noble-gas atoms, Phys. Rev. A90, 032712 (2014)

    2014

  24. [24]

    D. G. Green and G. F. Gribakin,γspectra and enhance- mentfactorsforpositronannihilationwithcoreelectrons, Phys. Rev. Lett.114, 093201 (2015)

    2015

  25. [25]

    C. M. Rawlins, J. Hofierka, B. Cunningham, C. H. Pat- terson, and D. G. Green, Many-body theory calculations of positron scattering and annihilation inH 2,N 2, and CH4, Phys. Rev. Lett.130, 263001 (2023)

    2023

  26. [26]

    A. L. Fetter and J. D. Walecka,Quantum Theory of Many-Particle Systems(McGraw-Hill, 1971)

    1971

  27. [27]

    W. H. Dickhoff and D. Van Neck,Many-Body Theory Exposed!, 3rd ed. (World Scientific, 2025)

    2025

  28. [28]

    For the positron-molecule problem theGWdiagram alone is wholly deficient. The importance of the virtual- positroniumΓladder series arises from the fact that suc- cessive terms in the series contribute with equal sign, in contrast to the all-electron case in which the signs al- ternate, leading to substantial cancellation. The virtual- positronium ladder s...

  29. [29]

    J. P. Cassidy, J. Hofierka, B. Cunningham, C. M. Rawl- ins, C.H. Patterson, andD. G. Green, Many-body theory calculations of positron binding to halogenated hydrocar- bons, Phys. Rev. A109, L040801 (2024)

    2024

  30. [30]

    Hofierka, B

    J. Hofierka, B. Cunningham, and D. G. Green, Many- body theory calculations of positron binding to hydrogen cyanide, Eur. Phys. J. D78, 37 (2024)

    2024

  31. [31]

    Arthur-Baidoo, J

    E. Arthur-Baidoo, J. R. Danielson, C. M. Surko, J. P. Cassidy, S. K. Gregg, J. Hofierka, B. Cunningham, C. H. Patterson, and D. G. Green, Positron annihilation and binding in aromatic and other ring molecules, Phys. Rev. A109, 062801 (2024)

    2024

  32. [33]

    J.Hofierka, C.M.Rawlins, B.Cunningham, D.T.Waide, and D. G. Green, Many-body theory calculations of positron scattering and annihilation in noble-gas atoms via the solution of Bethe–Salpeter equations using the gaussian-basis code EXCITON+, Front. in Physics11 (2023)

    2023

  33. [34]

    S. K. Gregg, J. P. Cassidy, A. R. Swann, J. Hofierka, B. Cunningham, and D. G. Green, Many-body theory and gaussian-basis implementation of positron annihi- lationγ-ray spectra on polyatomic molecules (2025), arXiv:2502.12364

    arXiv 2025

  34. [35]

    J. P. Cassidy, J. Hofierka, B. Cunningham, and D. G. Green, Many-body theory calculations of positronic- bonded molecular dianions, J. Chem. Phys.160, 084304 (2024)

    2024

  35. [36]

    R. R. Riso, J. H. M. Trabski, F. Rossi, D. Green, and H. Koch, Coupled cluster theory for positron binding in anions and polyatomic molecules (2026), arXiv:2603.19948

    arXiv 2026

  36. [37]

    N. V. Prokof’ev and B. V. Svistunov, Polaron problem by Diagrammatic Quantum Monte Carlo, Physical Review Letters81, 2514 (1998)

    1998

  37. [38]

    Van Houcke, E

    K. Van Houcke, E. Kozik, N. Prokof’ev, and B. Svis- tunov, Diagrammatic Monte Carlo, Physics Procedia6, 95 (2010)

    2010

  38. [39]

    Van Houcke, F

    K. Van Houcke, F. Werner, E. Kozik, N. Prokof’ev, B. Svistunov, M. J. H. Ku, A. T. Sommer, L. W. Cheuk, A. Schirotzek, and M. W. Zwierlein, Feynman diagrams versus Fermi-gas Feynman emulator, Nature Physics8, 366 (2012)

    2012

  39. [40]

    Chen and K

    K. Chen and K. Haule, A combined variational and dia- grammatic quantum Monte Carlo approach to the many- electron problem, Nature Comms.10, 3725 (2019)

    2019

  40. [41]

    Šimkovic and R

    F. Šimkovic and R. Rossi, Many-configuration markov- chain Monte Carlo (2021) arXiv:2102.05613

    arXiv 2021

  41. [42]

    Azadi, A

    S. Azadi, A. Davydov, and E. Kozik,GWspace-time method: Energy band gap of solid hydrogen, Phys. Rev. B105, 155136 (2022)

    2022

  42. [43]

    Bighin, Q

    G. Bighin, Q. P. Ho, M. Lemeshko, and T. V. Tscherbul, Diagrammatic Monte Carlo for electronic correlation in 7 molecules: High-order many-body perturbation theory with low scaling, Physical Review B108, 045115 (2023)

    2023

  43. [44]

    Sturt and E

    J. Sturt and E. Kozik, Exploiting parallelism for fast Feynman diagrammatics (2024) arXiv:2502.10327, 2501.00675

    arXiv 2024

  44. [45]

    M.VanhoeckeandM.Schirò,DiagrammaticMonteCarlo for dissipative quantum impurity models, Phys. Rev. B 109, 125125 (2024)

    2024

  45. [46]

    Brolli, C

    S. Brolli, C. Barbieri, and E. Vigezzi, Diagrammatic Monte Carlo for finite systems at zero temperature, Phys. Rev. Lett.134, 182502 (2025)

    2025

  46. [47]

    Y. Luo, J. Park, and M. Bernardi, First-principles dia- grammatic Monte Carlo for electron–phonon interactions and polaron, Nature Phys.21, 1275 (2025)

    2025

  47. [48]

    Since the self-energy matrix is symmetric, only the upper triangle (i≤f) is sampled, with the calculation paral- lelised over unique(i, f)pairs

  48. [49]

    C. H. Patterson, Density fitting in periodic systems: Ap- plication to tdhf in diamond and oxides, J. Chem. Phys. 153, 064107 (2020)

    2020

  49. [50]

    For ex- ample, selecting only the positron–electron interaction withall otherssettozeroyields theΓladderseries

    The proposal probabilities are user-configurable. For ex- ample, selecting only the positron–electron interaction withall otherssettozeroyields theΓladderseries. When multiple interaction types are active, their proposal prob- abilities must be equal to satisfy detailed balance

  50. [51]

    Körle, On absolute summability by Riesz and gen- eralized Cesàro means

    H.-H. Körle, On absolute summability by Riesz and gen- eralized Cesàro means. I, Canadian J. Math.22, 202 (1970)

    1970

  51. [52]

    Boyle and M

    J. Boyle and M. Pindzola,Many-body atomic physics (Cambridge University Press, 1998)

    1998