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

Many-Body Perturbation Theory for Ab Initio Nuclear Structure

A. Tichai

Published Jan 1, 2017• 1 citation
Publisher

Abstract

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.

DOI
10.26083/tuprints-00006942

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

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

Ab initioBogoliubov coupled cluster theory for open-shell nuclei

A. Signoracci, T. Duguet, G. Hagen +1 more

Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass $A\ensuremath{\approx}130$ on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking.Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory.Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in $m$ scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation.Results: Proof-of-principle calculations in an ${N}_{\text{max}}=6$ spherical harmonic oscillator basis for $^{16,18}\mathrm{O}$ and $^{18}\mathrm{Ne}$ in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while $^{20}\mathrm{O}$ and $^{20}\mathrm{Mg}$ display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is relatively constant for all five nuclei, in both the Hartree-Fock-Bogoliubov and BCCD approximations.Conclusions: The newly developed many-body formalism increases the potential span of ab initio calculations based on single-reference coupled cluster techniques tremendously, i.e., potentially to reach several hundred additional midmass nuclei. The new formalism offers a wealth of potential applications and further extensions dedicated to the description of ground and excited states of open-shell nuclei. Short-term goals include the implementation of three-nucleon forces at the normal-ordered two-body level. Midterm extensions include the approximate treatment of triples corrections and the development of the equation-of-motion methodology to treat both excited states and odd nuclei. Long-term extensions include exact restoration of U(1) and SU(2) symmetries.

arXiv:1412.2696•Physical Review C•Jun 29, 2015•96 citations

Bogoliubov many-body perturbation theory for open-shell nuclei

A. Tichai, P. Arthuis, T. Duguet +3 more

A Rayleigh–Schrödinger many-body perturbation theory (MBPT) approach is introduced by making use of a particle-number-breaking Bogoliubov reference state to tackle (near-)degenerate open-shell fermionic systems. By choosing a reference state that solves the Hartree–Fock–Bogoliubov variational problem, the approach reduces to the well-tested Møller–Plesset, i.e., Hartree–Fock based, MBPT when applied to closed-shell systems. Due to its algorithmic simplicity, the newly developed framework provides a computationally simple yet accurate alternative to state-of-the-art non-perturbative many-body approaches. At the price of working in the quasi-particle basis associated with a single-particle basis of sufficient size, the computational scaling of the method is independent of the particle number. This paper presents the first realistic applications of the method ranging from the oxygen to the nickel isotopic chains on the basis of a modern nuclear Hamiltonian derived from chiral effective field theory.

arXiv:1806.10931•Physics Letters B•Sep 25, 2018•96 citations

Multi-reference many-body perturbation theory for nuclei I -- Novel PGCM-PT formalism

Mikaël Frosini, T. Duguet, J.-P. Ebran +1 more

Perturbative and non-perturbative expansion methods already constitute a tool of choice to perform ab initio calculations over a significant part of the nuclear chart. In this context, the categories of accessible nuclei directly reflect the class of unperturbed state employed in the formulation of the expansion. The present work generalizes to the nuclear many-body context the versatile method of Ref. \cite{burton20a} by formulating a perturbative expansion on top of a multi-reference unperturbed state mixing deformed non-orthogonal Bogoliubov vacua, i.e. a state obtained from the projected generator coordinate method (PGCM). Particular attention is paid to the part of the mixing taking care of the symmetry restoration, showing that it can be exactly contracted throughout the expansion, thus reducing significantly the dimensionality of the linear problem to be solved to extract perturbative corrections. While the novel expansion method, coined as PGCM-PT, reduces to the PGCM at lowest order, it reduces to single-reference perturbation theories in appropriate limits. Based on a PGCM unperturbed state capturing (strong) static correlations in a versatile and efficient fashion, PGCM-PT is indistinctly applicable to doubly closed-shell, singly open-shell and doubly open-shell nuclei. The remaining (weak) dynamical correlations are brought consistently through perturbative corrections. This symmetry-conserving multi-reference perturbation theory is state-specific and applies to both ground and excited PGCM unperturbed states, thus correcting each state belonging to the low-lying spectrum of the system under study.

arXiv:2110.15737•Oct 29, 2021•0 citations

Combining symmetry breaking and restoration with configuration interaction: A highly accurate many-body scheme applied to the pairing Hamiltonian

J. Ripoche, Denis Lacroix A, D. Gambacurta +2 more

Background: Ab initio many-body methods have been developed over the past ten years to address mid-mass nuclei. In their best current level of implementation, their accuracy is of the order of a few percent error on the ground-state correlation energy. Recently implemented variants of these methods are operating a breakthrough in the description of medium-mass open-shell nuclei at a polynomial computational cost while putting state-of-the-art models of internucleon interactions to the test.Purpose: As progress in the design of internucleon interactions is made, and as questions one wishes to answer are refined in connection with increasingly available experimental data, further efforts must be made to tailor many-body methods that can reach an even higher precision for an even larger number of observable quantum states or nuclei. The objective of the present work is to contribute to such a quest by designing and testing a new many-body scheme.Methods: We formulate a truncated configuration-interaction method that consists of diagonalizing the Hamiltonian in a highly truncated subspace of the total $N$-body Hilbert space. The reduced Hilbert space is generated via the particle-number projected BCS state along with projected seniority-zero two- and four-quasiparticle excitations. Furthermore, the extent by which the underlying BCS state breaks $\text{U}(1)$ symmetry is optimized in the presence of the projected two- and four-quasiparticle excitations. This constitutes an extension of the so-called restricted variation after projection method in use within the frame of multireference energy density functional calculations. The quality of the newly designed method is tested against exact solutions of the so-called attractive pairing Hamiltonian problem.Results: By construction, the method reproduces exact results for $N=2$ and $N=4$. For $N=(8,\phantom{\rule{0.16em}{0ex}}16,\phantom{\rule{0.16em}{0ex}}20)$, the error in the ground-state correlation energy is less than (0.006%, 0.1%, 0.15%) across the entire range of internucleon coupling defining the pairing Hamiltonian and driving the normal-to-superfluid quantum phase transition. The presently proposed method offers the advantage of automatic access to the low-lying spectroscopy, which it does with high accuracy.Conclusions: The numerical cost of the newly designed variational method is polynomial (${N}^{6}$) in system size. This method achieves unprecedented accuracy for the ground-state correlation energy, effective pairing gap, and one-body entropy as well as for the excitation energy of low-lying states of the attractive pairing Hamiltonian. This constitutes a sufficiently strong motivation to envision its application to realistic nuclear Hamiltonians in view of providing a complementary, accurate, and versatile ab initio description of mid-mass open-shell nuclei in the future.

arXiv:1610.04063•Physical Review C•Jan 24, 2017•21 citations
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.
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.
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.
Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass $A\ensuremath{\approx}130$ on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking.Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory.Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in mmm scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation.Results: Proof-of-principle calculations in an Nmax=6{N}_{\text{max}}=6Nmax​=6 spherical harmonic oscillator basis for 16,18O^{16,18}\mathrm{O}16,18O and 18Ne^{18}\mathrm{Ne}18Ne in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while 20O^{20}\mathrm{O}20O and 20Mg^{20}\mathrm{Mg}20Mg display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is relatively constant for all five nuclei, in both the Hartree-Fock-Bogoliubov and BCCD approximations.Conclusions: The newly developed many-body formalism increases the potential span of ab initio calculations based on single-reference coupled cluster techniques tremendously, i.e., potentially to reach several hundred additional midmass nuclei. The new formalism offers a wealth of potential applications and further extensions dedicated to the description of ground and excited states of open-shell nuclei. Short-term goals include the implementation of three-nucleon forces at the normal-ordered two-body level. Midterm extensions include the approximate treatment of triples corrections and the development of the equation-of-motion methodology to treat both excited states and odd nuclei. Long-term extensions include exact restoration of U(1) and SU(2) symmetries.
A Rayleigh–Schrödinger many-body perturbation theory (MBPT) approach is introduced by making use of a particle-number-breaking Bogoliubov reference state to tackle (near-)degenerate open-shell fermionic systems. By choosing a reference state that solves the Hartree–Fock–Bogoliubov variational problem, the approach reduces to the well-tested Møller–Plesset, i.e., Hartree–Fock based, MBPT when applied to closed-shell systems. Due to its algorithmic simplicity, the newly developed framework provides a computationally simple yet accurate alternative to state-of-the-art non-perturbative many-body approaches. At the price of working in the quasi-particle basis associated with a single-particle basis of sufficient size, the computational scaling of the method is independent of the particle number. This paper presents the first realistic applications of the method ranging from the oxygen to the nickel isotopic chains on the basis of a modern nuclear Hamiltonian derived from chiral effective field theory.
Perturbative and non-perturbative expansion methods already constitute a tool of choice to perform ab initio calculations over a significant part of the nuclear chart. In this context, the categories of accessible nuclei directly reflect the class of unperturbed state employed in the formulation of the expansion. The present work generalizes to the nuclear many-body context the versatile method of Ref. \cite{burton20a} by formulating a perturbative expansion on top of a multi-reference unperturbed state mixing deformed non-orthogonal Bogoliubov vacua, i.e. a state obtained from the projected generator coordinate method (PGCM). Particular attention is paid to the part of the mixing taking care of the symmetry restoration, showing that it can be exactly contracted throughout the expansion, thus reducing significantly the dimensionality of the linear problem to be solved to extract perturbative corrections. While the novel expansion method, coined as PGCM-PT, reduces to the PGCM at lowest order, it reduces to single-reference perturbation theories in appropriate limits. Based on a PGCM unperturbed state capturing (strong) static correlations in a versatile and efficient fashion, PGCM-PT is indistinctly applicable to doubly closed-shell, singly open-shell and doubly open-shell nuclei. The remaining (weak) dynamical correlations are brought consistently through perturbative corrections. This symmetry-conserving multi-reference perturbation theory is state-specific and applies to both ground and excited PGCM unperturbed states, thus correcting each state belonging to the low-lying spectrum of the system under study.
Background: Ab initio many-body methods have been developed over the past ten years to address mid-mass nuclei. In their best current level of implementation, their accuracy is of the order of a few percent error on the ground-state correlation energy. Recently implemented variants of these methods are operating a breakthrough in the description of medium-mass open-shell nuclei at a polynomial computational cost while putting state-of-the-art models of internucleon interactions to the test.Purpose: As progress in the design of internucleon interactions is made, and as questions one wishes to answer are refined in connection with increasingly available experimental data, further efforts must be made to tailor many-body methods that can reach an even higher precision for an even larger number of observable quantum states or nuclei. The objective of the present work is to contribute to such a quest by designing and testing a new many-body scheme.Methods: We formulate a truncated configuration-interaction method that consists of diagonalizing the Hamiltonian in a highly truncated subspace of the total NNN-body Hilbert space. The reduced Hilbert space is generated via the particle-number projected BCS state along with projected seniority-zero two- and four-quasiparticle excitations. Furthermore, the extent by which the underlying BCS state breaks U(1)\text{U}(1)U(1) symmetry is optimized in the presence of the projected two- and four-quasiparticle excitations. This constitutes an extension of the so-called restricted variation after projection method in use within the frame of multireference energy density functional calculations. The quality of the newly designed method is tested against exact solutions of the so-called attractive pairing Hamiltonian problem.Results: By construction, the method reproduces exact results for N=2N=2N=2 and N=4N=4N=4. For N=(8,16,20)N=(8,\phantom{\rule{0.16em}{0ex}}16,\phantom{\rule{0.16em}{0ex}}20)N=(8,16,20), the error in the ground-state correlation energy is less than (0.006%, 0.1%, 0.15%) across the entire range of internucleon coupling defining the pairing Hamiltonian and driving the normal-to-superfluid quantum phase transition. The presently proposed method offers the advantage of automatic access to the low-lying spectroscopy, which it does with high accuracy.Conclusions: The numerical cost of the newly designed variational method is polynomial (N6{N}^{6}N6) in system size. This method achieves unprecedented accuracy for the ground-state correlation energy, effective pairing gap, and one-body entropy as well as for the excitation energy of low-lying states of the attractive pairing Hamiltonian. This constitutes a sufficiently strong motivation to envision its application to realistic nuclear Hamiltonians in view of providing a complementary, accurate, and versatile ab initio description of mid-mass open-shell nuclei in the future.