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:2605.30527

Time-ordered Diagrammatic Monte Carlo for atomic nuclei

Stefano Brolli, C. Barbieri

Published May 28, 2026• 0 citations
arXivHTMLPDF

Abstract

Diagrammatic Monte Carlo provides a systematically improvable framework for stochastically resumming many-body expansions to high orders through direct sampling of diagram topologies. We advance our earlier work by introducing a novel time-ordered Diagrammatic Monte Carlo algorithm for the single-particle Green's function. The algorithm is tailored to finite nuclei, formulated in discrete model spaces and applicable to arbitrary two-body interactions. The new time-ordered diagrammatic Monte Carlo algorithm is based on the on-the-fly evaluation of time-ordered Goldstone diagrams, avoiding explicit diagram enumeration and expensive frequency integration. We show the algorithm by computing ${}^{16}$O up to fifth order in a reduced model space using optimized reference state orbitals and including effective three-body forces. Benchmarking against established truncation schemes in ab initio nuclear theory demonstrates its potential to overcome the limitations of current many-body approaches.

DOI
10.48550/arxiv.2605.30527

Related papers

Diagrammatic Monte Carlo for Finite Systems at Zero Temperature

Stefano Brolli, C. Barbieri, E. Vigezzi

High-order virtual excitations play an important role in microscopic models of nuclear reactions at intermediate energies. However, the factorial growth of their complexity has prevented their consistent inclusion in ab initio many-body calculations. For infinite systems at finite temperature, such drawbacks can be overcome using diagrammatic Monte Carlo techniques to resum entire series of Feynman diagrams. We present a diagrammatic Monte Carlo algorithm that can be applied to self-bound systems with discrete energy levels at zero temperature, and demonstrate its potential for the Richardson model of nuclear pairing. We show that sampling the topological space of diagrams allows the inclusion of high-order excitations that are neglected in state-of-the-art approximations used in nuclear physics and quantum chemistry. We propose that sampling the diagrammatic space can overcome the long-standing gap between our microscopic understanding of structure and reactions in nuclear physics.

arXiv:2501.02646•Physical Review Letters•May 9, 2025•1 citation

Stochastic, real-space, imaginary-time evaluation of third-order Feynman–Goldstone diagrams

Soohaeng Yoo Willow, So Hirata

A new, alternative set of interpretation rules of Feynman-Goldstone diagrams for many-body perturbation theory is proposed, which translates diagrams into algebraic expressions suitable for direct Monte Carlo integrations. A vertex of a diagram is associated with a Coulomb interaction (rather than a two-electron integral) and an edge with the trace of a Green's function in real space and imaginary time. With these, 12 diagrams of third-order many-body perturbation (MP3) theory are converted into 20-dimensional integrals, which are then evaluated by a Monte Carlo method. It uses redundant walkers for convergence acceleration and a weight function for importance sampling in conjunction with the Metropolis algorithm. The resulting Monte Carlo MP3 method has low-rank polynomial size dependence of the operation cost, a negligible memory cost, and a naturally parallel computational kernel, while reproducing the correct correlation energies of small molecules within a few mEh after 10(6) Monte Carlo steps.

The Journal of Chemical Physics•Jan 13, 2014•22 citations

Diagrammatic Methods for Strongly-Correlated Electrons in Realistic Systems

Jia Li

This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.

•Jan 1, 2021•0 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

Diagrammatic Monte Carlo Approach to Angular Momentum in Quantum Many-Particle Systems

Giacomo Bighin, Timur V. Tscherbul, Mikhail Lemeshko

We introduce a diagrammatic Monte Carlo approach to angular momentum properties of quantum many-particle systems possessing a macroscopic number of degrees of freedom. The treatment is based on a diagrammatic expansion that merges the usual Feynman diagrams with the angular momentum diagrams known from atomic and nuclear structure theory, thereby incorporating the non-Abelian algebra inherent to quantum rotations. Our approach is applicable at arbitrary coupling, is free of systematic errors and of finite-size effects, and naturally provides access to the impurity Green function. We exemplify the technique by obtaining an all-coupling solution of the angulon model; however, the method is quite general and can be applied to a broad variety of systems in which particles exchange quantum angular momentum with their many-body environment.

arXiv:1803.07990•Physical Review Letters•Oct 16, 2018•18 citations

Bogoliubov Many-Body Perturbation Theory for Nuclei : Systematic Generation and Evaluation of Diagrams and First ab initio Calculations

P. Arthuis

The last few decades in nuclear structure theory have seen a rapid expansion of ab initio theories, aiming at describing the properties of nuclei starting from the inter-nucleonic interaction. Such an expansion relied both on the tremendous growth of computing power and novel formal developments. This work focuses on the development of the recently proposed Bogoliubov Many-Body Perturbation Theory that relies on a particle-number-breaking reference state to tackle singly open-shell nuclei. The formalism is first described in details, and diagrammatic and algebraic contributions are derived up to second order. Its link to standard Many-Body Perturbation Theory is made explicit, as well as its connexion to Bogoliubov Coupled-Cluster theory. An automated extension to higher orders based on graph theory methods is then detailed, and the ADG numerical program generating and evaluating BMBPT diagrams at arbitrary order is introduced. Such a formal development carries implications that are not restricted to the present work, as the developed methods can be applied to other many-body methods. Finally, first numerical results obtained for oxygen, calcium and nickel isotopes are presented. They establish BMBPT as a method of interest for large-scale computations of isotopic or isotonic chains in the mid-mass sector of the nuclear chart.

•Sep 27, 2018•0 citations
Diagrammatic Monte Carlo provides a systematically improvable framework for stochastically resumming many-body expansions to high orders through direct sampling of diagram topologies. We advance our earlier work by introducing a novel time-ordered Diagrammatic Monte Carlo algorithm for the single-particle Green's function. The algorithm is tailored to finite nuclei, formulated in discrete model spaces and applicable to arbitrary two-body interactions. The new time-ordered diagrammatic Monte Carlo algorithm is based on the on-the-fly evaluation of time-ordered Goldstone diagrams, avoiding explicit diagram enumeration and expensive frequency integration. We show the algorithm by computing 16{}^{16}16O up to fifth order in a reduced model space using optimized reference state orbitals and including effective three-body forces. Benchmarking against established truncation schemes in ab initio nuclear theory demonstrates its potential to overcome the limitations of current many-body approaches.
High-order virtual excitations play an important role in microscopic models of nuclear reactions at intermediate energies. However, the factorial growth of their complexity has prevented their consistent inclusion in ab initio many-body calculations. For infinite systems at finite temperature, such drawbacks can be overcome using diagrammatic Monte Carlo techniques to resum entire series of Feynman diagrams. We present a diagrammatic Monte Carlo algorithm that can be applied to self-bound systems with discrete energy levels at zero temperature, and demonstrate its potential for the Richardson model of nuclear pairing. We show that sampling the topological space of diagrams allows the inclusion of high-order excitations that are neglected in state-of-the-art approximations used in nuclear physics and quantum chemistry. We propose that sampling the diagrammatic space can overcome the long-standing gap between our microscopic understanding of structure and reactions in nuclear physics.
A new, alternative set of interpretation rules of Feynman-Goldstone diagrams for many-body perturbation theory is proposed, which translates diagrams into algebraic expressions suitable for direct Monte Carlo integrations. A vertex of a diagram is associated with a Coulomb interaction (rather than a two-electron integral) and an edge with the trace of a Green's function in real space and imaginary time. With these, 12 diagrams of third-order many-body perturbation (MP3) theory are converted into 20-dimensional integrals, which are then evaluated by a Monte Carlo method. It uses redundant walkers for convergence acceleration and a weight function for importance sampling in conjunction with the Metropolis algorithm. The resulting Monte Carlo MP3 method has low-rank polynomial size dependence of the operation cost, a negligible memory cost, and a naturally parallel computational kernel, while reproducing the correct correlation energies of small molecules within a few mEh after 10(6) Monte Carlo steps.
This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.
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.
We introduce a diagrammatic Monte Carlo approach to angular momentum properties of quantum many-particle systems possessing a macroscopic number of degrees of freedom. The treatment is based on a diagrammatic expansion that merges the usual Feynman diagrams with the angular momentum diagrams known from atomic and nuclear structure theory, thereby incorporating the non-Abelian algebra inherent to quantum rotations. Our approach is applicable at arbitrary coupling, is free of systematic errors and of finite-size effects, and naturally provides access to the impurity Green function. We exemplify the technique by obtaining an all-coupling solution of the angulon model; however, the method is quite general and can be applied to a broad variety of systems in which particles exchange quantum angular momentum with their many-body environment.
The last few decades in nuclear structure theory have seen a rapid expansion of ab initio theories, aiming at describing the properties of nuclei starting from the inter-nucleonic interaction. Such an expansion relied both on the tremendous growth of computing power and novel formal developments. This work focuses on the development of the recently proposed Bogoliubov Many-Body Perturbation Theory that relies on a particle-number-breaking reference state to tackle singly open-shell nuclei. The formalism is first described in details, and diagrammatic and algebraic contributions are derived up to second order. Its link to standard Many-Body Perturbation Theory is made explicit, as well as its connexion to Bogoliubov Coupled-Cluster theory. An automated extension to higher orders based on graph theory methods is then detailed, and the ADG numerical program generating and evaluating BMBPT diagrams at arbitrary order is introduced. Such a formal development carries implications that are not restricted to the present work, as the developed methods can be applied to other many-body methods. Finally, first numerical results obtained for oxygen, calcium and nickel isotopes are presented. They establish BMBPT as a method of interest for large-scale computations of isotopic or isotonic chains in the mid-mass sector of the nuclear chart.