Search papers, labs, and topics across Lattice.
This paper introduces a diagrammatic Monte Carlo method to evaluate the ladder series contributions to the correlation potential of a positron interacting with a molecule, specifically using the $GW$@TDHF and positron-hole Goldstone ladder series. The approach employs stochastic sampling within the Tamm-Dancoff approximation and utilizes Cesáro-Riesz resummation to achieve an efficient infinite-order extrapolation. Benchmark results for lithium hydride reveal that this method significantly reduces memory requirements while maintaining quantitative accuracy compared to exact diagonalization techniques, demonstrating its effectiveness in many-body theory.
Stochastic sampling of positron-molecule interactions achieves a dramatic reduction in memory usage while accurately capturing complex many-body effects.
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{á}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.