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

High-order perturbative calculations of nuclear ground states: Automated evaluation of many-body diagrams

Zhen Li, A. Tichai, A. Schwenk, N. A. Smirnova

Published Jun 30, 2026• 0 citations
arXivHTMLPDF

Abstract

We advance the many-body perturbation theory (MBPT) calculations of the ground-state energy and radius of closed-shell nuclei beyond third order. Using automated diagram generation and evaluation up to fifth order, we present ground-state properties of selected closed-shell nuclei up to $^{78}$Ni with two- and three-nucleon interactions derived from chiral effective field theory. A clear convergence trend is observed for the ground-state energy enabling calculations at improved accuracy. We further investigate in detail the decomposition of the fourth-order contributions. For the ground-state energy, the magnitude of the fourth-order contribution is typically less than half of the third order, and a typical cancellation among different classes of diagrams is observed. Finally, we perform a comprehensive comparison between MBPT and non-perturbative in-medium similarity renormalization group (IMSRG) calculations, with the goal to provide insight into many-body uncertainties associated with the IMSRG(2) truncation.

DOI
10.48550/arxiv.2607.00119

Related papers

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

Tensor-decomposition techniques for ab initio nuclear structure calculations: From chiral nuclear potentials to ground-state energies

A. Tichai, Roman Schutski, Gustavo E. Scuseria +1 more

Background: The computational resources needed to generate the ab initio solution of the nuclear many-body problem for increasing mass number and/or accuracy necessitates innovative developments to improve upon (i) the storage of many-body operators and (ii) the scaling of many-body methods used to evaluate nuclear observables. The storing and efficient handling of many-body operators with high particle ranks is currently one of the major bottlenecks limiting the applicability range of ab initio studies with respect to mass number and accuracy. Recently, the application of tensor decomposition techniques to many-body tensors has proven highly beneficial to reduce the computational cost of ab initio calculations in quantum chemistry and solid-state physics.Purpose: The impact of applying state-of-the-art tensor factorization techniques to modern nuclear Hamiltonians derived from chiral effective field theory is investigated. Subsequently, the error induced by the tensor decomposition of the input Hamiltonian on ground-state energies of closed-shell nuclei calculated via second-order many-body perturbation theory is benchmarked.Methods: The first proof-of-principles application of tensor-decomposition techniques to the nuclear Hamiltonian is performed. Two different tensor formats are investigated by systematically benchmarking the approximation error on matrix elements stored in various bases of interest. The analysis is achieved while including normal-ordered three-nucleon interactions that are nowadays used as input to the most advanced ab initio calculations in medium-mass nuclei. With the aid of the factorized Hamiltonian, the second-order perturbative correction to ground-state energies is decomposed and the scaling properties of the underlying tensor network are discussed.Results: The employed tensor formats are found to lead to an efficient data compression of two-body matrix elements of the nuclear Hamiltonian. In particular, the sophisticated tensor hypercontraction scheme yields low tensor ranks with respect to both harmonic-oscillator and Hartree-Fock single-particle bases. It is found that the tensor rank depends on the two-body total angular momentum $J$ for which one performs the decomposition, which is itself directly related to the sparsity of the corresponding tensor. Furthermore, including normal-ordered two-body contributions originating from three-body interactions does not compromise the efficient data compression. Ultimately, the use of factorized matrix elements authorizes controlled approximations of the exact second-order ground-state energy corrections. In particular, a small enough error is obtained from low-rank factorizations in $^{4}\mathrm{He},^{16}\mathrm{O}$, and $^{40}\mathrm{Ca}$.Conclusions: It is presently demonstrated that tensor-decomposition techniques can be efficiently applied to systematically approximate the nuclear many-body Hamiltonian in terms of lower-rank tensors. Beyond the input Hamiltonian, tensor-decomposition techniques can be envisioned to scale down the cost of state-of-the-art nonperturbative many-body methods in order to extend ab initio studies to (i) higher precisions, (ii) larger masses, and (iii) nuclei of doubly open-shell character.

arXiv:1810.08419•Physical Review C•Mar 18, 2019•22 citations

Ab initio Rayleigh-Schrödinger perturbation calculation including three-body force

B. S. Hu, T. Li, F. R. Xu

We first derive the Rayleigh-Schrödinger many-body perturbation theory up to third order (RSPT3) for Hamiltonians with three-body interaction. The structure of closed-shell nuclei in a wide mass range from 4He to 48Ca has been investigated by the RSPT3 with explicit NN+3N Hamiltonian. The RSPT3 calculations are performed within Hartree-Fock bases. The perturbative contribution of antisymmetrized Goldstone diagrams (diagrammatic expansion for RSPT) with normal-ordered interaction has been analyzed. We demonstrate that the normal-ordered two-body level (NO2B) approximation which neglects the residual three-body term can catch the main effect of full three-body force in RSPT calculation. We also present rigorous benchmarks for RSPT3 with non-perturbative coupled cluster and in-medium similarity renormalization group using the same chiral NNLOsat NO2B interaction. The three methods are in good agreement with each other for binding energies and radii. However, the RSPT3 provides an alternative that is computationally inexpensive but comparable in accuracy to state-of-the-art non-perturbative many-body methods.

arXiv:1810.08804•Oct 20, 2018•1 citation

Projected Bogoliubov Many-Body Perturbation Theory : Overcoming formal and technical challenges

Julien Ripoche

Given elementary inter-nucleon interactions,the resolution of the A-body Schrödinger equation gives access to the properties of the nuclear eigenstates. Many-body perturbation theory (MBPT) based on a Hartree-Fock mean field allows one to treat the dynamical correlations at play in doubly-closed-shell nuclei. Beyond shell closures,breaking U(1) symmetry associated with the conservation of neutron and proton numbers further permits to capture static correlations at the meanfield level and to formulate the well-behaved Bogoliubov many-body perturbation theory (BMBPT).Nevertheless, BMBPT results show contaminations associated with the broken symmetry that is only emergent on finite quantum system such as theatomic nucleus. Thus, the restoration of U(1) symmetry beyond the mean field is necessary for acorrect description and gives rise to the projected BMBPT method (PBMBPT). The goal is to implement PBMBPT in order to perform ab initio calculations of singly-open-shell mid-mass nuclei.The present work provides systematic solutions to the formal and technical problems arising during the implementation of PBMBPT: an automatic and safe method to generate PBMBPT diagrams and associated expressions, a framework to calculate other observables than the energy, an extension of the normal ordering approximation to symmetry breaking methods as well as an automatic and safe tool for spherical reduction of tensor networks.

•Sep 24, 2019•2 citations

Many-Body Perturbation Theory for Ab Initio Nuclear Structure

A. Tichai

The solution of the quantum many-body problem for medium-mass nuclei using realistic nuclear interactions poses a superbe challenge for nuclear structure research. Because an exact solution can only be provided for the lightest nuclei, one has to rely on approximate solutions when proceeding to heavier systems. Over the past years, tremendous progress has been made in the development and application of systematically improvable expansion methods and an accurate description of nuclear observables has become viable up to mass number $A \approx 100$. While closed-shell systems are consistently described via a plethora of different many-body methods, the extension to genuine open-shell systems still remains a major challenge and up to now there is no \textit{ab initio} many-body method which applies equally well to systems with even and odd mass numbers. The goal of this thesis is the development and implementation of innovative perturbative approaches with genuine open-shell capabilities. This requires the extension of well-known single-reference approaches to more general vacua. In this work we choose two complementary routes for the usage of generalized reference states. First, we derive a new \emph{ab initio} approach based on multi-configurational reference states that are conveniently derived from a prior no-core shell model calculation. Perturbative corrections are derived via second-order many-body perturbation theory, thus, merging configuration interaction and many-body perturbation theory. The generality of this ansatz enables for a treatment of medium-mass systems with arbitrary mass number, as well as the extension to low-lying excited states such that ground and excited states are treated on an equal footing. In a complementary approach, we use reference states that break a symmetry of the underlying Hamiltonian. In the simplest case this corresponds to the expansion around a particle-number-broken Hartree-Fock-Bogoliubov vacuum which is obtained from a mean-field calculation. Pairing correlations are already absorbed in the reference state. The mild scaling allows for investigations up to tin isotopic chains. All benchmarks were performed using state-of-the-art chiral two- plus three-nucleon interactions thus allowing for a universal description of nuclear observables in the medium-mass regime.

•Jan 1, 2017•1 citation

In-medium similarity renormalization group for closed and open-shell nuclei

H Hergert

Abstract We present a pedagogical introduction to the in-medium similarity renormalization group (IMSRG) framework for ab initio calculations of nuclei. The IMSRG performs continuous unitary transformations of the nuclear many-body Hamiltonian in second-quantized form, which can be implemented with polynomial computational effort. Through suitably chosen generators, it is possible to extract eigenvalues of the Hamiltonian in a given nucleus, or drive the Hamiltonian matrix in configuration space to specific structures, e.g., band- or block-diagonal form. Exploiting this flexibility, we describe two complementary approaches for the description of closed- and open-shell nuclei: the first is the multireference IMSRG (MR-IMSRG), which is designed for the efficient calculation of nuclear ground-state properties. The second is the derivation of non-empirical valence-space interactions that can be used as input for nuclear shell model (i.e., configuration interaction (CI)) calculations. This IMSRG+shell model approach provides immediate access to excitation spectra, transitions, etc, but is limited in applicability by the factorial cost of the CI calculations. We review applications of the MR-IMSRG and IMSRG+shell model approaches to the calculation of ground-state properties for the oxygen, calcium, and nickel isotopic chains or the spectroscopy of nuclei in the lower sd shell, respectively, and present selected new results, e.g., for the ground- and excited state properties of neon isotopes.

arXiv:1607.06882•Physica Scripta•Dec 19, 2016•96 citations
We advance the many-body perturbation theory (MBPT) calculations of the ground-state energy and radius of closed-shell nuclei beyond third order. Using automated diagram generation and evaluation up to fifth order, we present ground-state properties of selected closed-shell nuclei up to 78^{78}78Ni with two- and three-nucleon interactions derived from chiral effective field theory. A clear convergence trend is observed for the ground-state energy enabling calculations at improved accuracy. We further investigate in detail the decomposition of the fourth-order contributions. For the ground-state energy, the magnitude of the fourth-order contribution is typically less than half of the third order, and a typical cancellation among different classes of diagrams is observed. Finally, we perform a comprehensive comparison between MBPT and non-perturbative in-medium similarity renormalization group (IMSRG) calculations, with the goal to provide insight into many-body uncertainties associated with the IMSRG(2) truncation.
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.
Background: The computational resources needed to generate the ab initio solution of the nuclear many-body problem for increasing mass number and/or accuracy necessitates innovative developments to improve upon (i) the storage of many-body operators and (ii) the scaling of many-body methods used to evaluate nuclear observables. The storing and efficient handling of many-body operators with high particle ranks is currently one of the major bottlenecks limiting the applicability range of ab initio studies with respect to mass number and accuracy. Recently, the application of tensor decomposition techniques to many-body tensors has proven highly beneficial to reduce the computational cost of ab initio calculations in quantum chemistry and solid-state physics.Purpose: The impact of applying state-of-the-art tensor factorization techniques to modern nuclear Hamiltonians derived from chiral effective field theory is investigated. Subsequently, the error induced by the tensor decomposition of the input Hamiltonian on ground-state energies of closed-shell nuclei calculated via second-order many-body perturbation theory is benchmarked.Methods: The first proof-of-principles application of tensor-decomposition techniques to the nuclear Hamiltonian is performed. Two different tensor formats are investigated by systematically benchmarking the approximation error on matrix elements stored in various bases of interest. The analysis is achieved while including normal-ordered three-nucleon interactions that are nowadays used as input to the most advanced ab initio calculations in medium-mass nuclei. With the aid of the factorized Hamiltonian, the second-order perturbative correction to ground-state energies is decomposed and the scaling properties of the underlying tensor network are discussed.Results: The employed tensor formats are found to lead to an efficient data compression of two-body matrix elements of the nuclear Hamiltonian. In particular, the sophisticated tensor hypercontraction scheme yields low tensor ranks with respect to both harmonic-oscillator and Hartree-Fock single-particle bases. It is found that the tensor rank depends on the two-body total angular momentum JJJ for which one performs the decomposition, which is itself directly related to the sparsity of the corresponding tensor. Furthermore, including normal-ordered two-body contributions originating from three-body interactions does not compromise the efficient data compression. Ultimately, the use of factorized matrix elements authorizes controlled approximations of the exact second-order ground-state energy corrections. In particular, a small enough error is obtained from low-rank factorizations in 4He,16O^{4}\mathrm{He},^{16}\mathrm{O}4He,16O, and 40Ca^{40}\mathrm{Ca}40Ca.Conclusions: It is presently demonstrated that tensor-decomposition techniques can be efficiently applied to systematically approximate the nuclear many-body Hamiltonian in terms of lower-rank tensors. Beyond the input Hamiltonian, tensor-decomposition techniques can be envisioned to scale down the cost of state-of-the-art nonperturbative many-body methods in order to extend ab initio studies to (i) higher precisions, (ii) larger masses, and (iii) nuclei of doubly open-shell character.
We first derive the Rayleigh-Schrödinger many-body perturbation theory up to third order (RSPT3) for Hamiltonians with three-body interaction. The structure of closed-shell nuclei in a wide mass range from 4He to 48Ca has been investigated by the RSPT3 with explicit NN+3N Hamiltonian. The RSPT3 calculations are performed within Hartree-Fock bases. The perturbative contribution of antisymmetrized Goldstone diagrams (diagrammatic expansion for RSPT) with normal-ordered interaction has been analyzed. We demonstrate that the normal-ordered two-body level (NO2B) approximation which neglects the residual three-body term can catch the main effect of full three-body force in RSPT calculation. We also present rigorous benchmarks for RSPT3 with non-perturbative coupled cluster and in-medium similarity renormalization group using the same chiral NNLOsat NO2B interaction. The three methods are in good agreement with each other for binding energies and radii. However, the RSPT3 provides an alternative that is computationally inexpensive but comparable in accuracy to state-of-the-art non-perturbative many-body methods.
Given elementary inter-nucleon interactions,the resolution of the A-body Schrödinger equation gives access to the properties of the nuclear eigenstates. Many-body perturbation theory (MBPT) based on a Hartree-Fock mean field allows one to treat the dynamical correlations at play in doubly-closed-shell nuclei. Beyond shell closures,breaking U(1) symmetry associated with the conservation of neutron and proton numbers further permits to capture static correlations at the meanfield level and to formulate the well-behaved Bogoliubov many-body perturbation theory (BMBPT).Nevertheless, BMBPT results show contaminations associated with the broken symmetry that is only emergent on finite quantum system such as theatomic nucleus. Thus, the restoration of U(1) symmetry beyond the mean field is necessary for acorrect description and gives rise to the projected BMBPT method (PBMBPT). The goal is to implement PBMBPT in order to perform ab initio calculations of singly-open-shell mid-mass nuclei.The present work provides systematic solutions to the formal and technical problems arising during the implementation of PBMBPT: an automatic and safe method to generate PBMBPT diagrams and associated expressions, a framework to calculate other observables than the energy, an extension of the normal ordering approximation to symmetry breaking methods as well as an automatic and safe tool for spherical reduction of tensor networks.
The solution of the quantum many-body problem for medium-mass nuclei using realistic nuclear interactions poses a superbe challenge for nuclear structure research. Because an exact solution can only be provided for the lightest nuclei, one has to rely on approximate solutions when proceeding to heavier systems. Over the past years, tremendous progress has been made in the development and application of systematically improvable expansion methods and an accurate description of nuclear observables has become viable up to mass number A≈100A \approx 100A≈100. While closed-shell systems are consistently described via a plethora of different many-body methods, the extension to genuine open-shell systems still remains a major challenge and up to now there is no \textit{ab initio} many-body method which applies equally well to systems with even and odd mass numbers. The goal of this thesis is the development and implementation of innovative perturbative approaches with genuine open-shell capabilities. This requires the extension of well-known single-reference approaches to more general vacua. In this work we choose two complementary routes for the usage of generalized reference states. First, we derive a new \emph{ab initio} approach based on multi-configurational reference states that are conveniently derived from a prior no-core shell model calculation. Perturbative corrections are derived via second-order many-body perturbation theory, thus, merging configuration interaction and many-body perturbation theory. The generality of this ansatz enables for a treatment of medium-mass systems with arbitrary mass number, as well as the extension to low-lying excited states such that ground and excited states are treated on an equal footing. In a complementary approach, we use reference states that break a symmetry of the underlying Hamiltonian. In the simplest case this corresponds to the expansion around a particle-number-broken Hartree-Fock-Bogoliubov vacuum which is obtained from a mean-field calculation. Pairing correlations are already absorbed in the reference state. The mild scaling allows for investigations up to tin isotopic chains. All benchmarks were performed using state-of-the-art chiral two- plus three-nucleon interactions thus allowing for a universal description of nuclear observables in the medium-mass regime.
Abstract We present a pedagogical introduction to the in-medium similarity renormalization group (IMSRG) framework for ab initio calculations of nuclei. The IMSRG performs continuous unitary transformations of the nuclear many-body Hamiltonian in second-quantized form, which can be implemented with polynomial computational effort. Through suitably chosen generators, it is possible to extract eigenvalues of the Hamiltonian in a given nucleus, or drive the Hamiltonian matrix in configuration space to specific structures, e.g., band- or block-diagonal form. Exploiting this flexibility, we describe two complementary approaches for the description of closed- and open-shell nuclei: the first is the multireference IMSRG (MR-IMSRG), which is designed for the efficient calculation of nuclear ground-state properties. The second is the derivation of non-empirical valence-space interactions that can be used as input for nuclear shell model (i.e., configuration interaction (CI)) calculations. This IMSRG+shell model approach provides immediate access to excitation spectra, transitions, etc, but is limited in applicability by the factorial cost of the CI calculations. We review applications of the MR-IMSRG and IMSRG+shell model approaches to the calculation of ground-state properties for the oxygen, calcium, and nickel isotopic chains or the spectroscopy of nuclei in the lower sd shell, respectively, and present selected new results, e.g., for the ground- and excited state properties of neon isotopes.