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

Multi-reference many-body perturbation theory for nuclei III -- Ab initio calculations at second order in PGCM-PT

Mikaël Frosini, T. Duguet, J.-P. Ebran, B. Bally, H. Hergert, Tomás R. Rodríguez, Robert Roth, J. M. Yao, V. Somà

Published Nov 2, 2021• 4 citations
arXivHTMLPDF

Abstract

In spite of missing dynamical correlations, the projected generator coordinate method (PGCM) was recently shown to be a suitable method to tackle the low-lying spectroscopy of complex nuclei. Still, describing absolute binding energies and reaching high accuracy eventually requires the inclusion of dynamical correlations on top of the PGCM. In this context, the present work discusses the first realistic results of a novel multi-reference perturbation theory (PGCM-PT) that can do so within a symmetry-conserving scheme for both ground and low-lying excited states. First, proof-of-principle calculations in a small ($e_{\mathrm{max}}=4$) model space demonstrate that exact binding energies of closed- (\nucl{O}{16}) and open-shell (\nucl{O}{18}, \nucl{Ne}{20}) nuclei are reproduced within $0.5-1.5\%$ at second order, i.e. through PGCM-PT(2). Moreover, profiting from the pre-processing of the Hamiltonian via multi-reference in-medium similarity renormalization group transformations, PGCM-PT(2) can reach converged values within smaller model spaces than with an unevolved Hamiltonian. Doing so, dynamical correlations captured by PGCM-PT(2) are shown to bring essential corrections to low-lying excitation energies that become too dilated at leading order, i.e., at the strict PGCM level. The present work is laying the foundations for a better understanding of the optimal way to grasp static and dynamical correlations in a consistent fashion, with the aim of accurately describing ground and excited states of complex nuclei via ab initio many-body methods.

DOI
10.48550/arxiv.2111.01461

Related papers

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

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

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

Emulation of PGCM calculations using the Eigenvector continuation method

A. Y. Le Roux

An atomic nucleus is a quantum system of interacting nucleons and constitutes a problem difficult to solve exactly. For this reason, a diversity of approximate resolution methods has been designed, and Projected Generator Coordinate Method (PGCM) is one of them. The strong point of PGCM is to construct a physically inspired small dimensional space, in which an approximate solution of the nuclear many-body problem is easily found. However the numerical cost of PGCM space computation make this method inadapted for sensibility analysis of nuclear observables with restect to parametrisation of the interaction model, this analysis requiring an huge number of PGCM computations. In order to make this type of study possible, this thesis explore the concept of PGCM emulator. In this work, a combination of PGCM with Eigenvector Continuation (EC) is constructed and studied. This combination (the PGCM-EC emulator) takes advantage of mathematical similarities between PGCM and EC, and above all of the decomposition of the hamiltonian as a linear combination of parameter-independent terms. The latter property is used to concentrate the heavier numerical cost in the computation of parameter-independent quantities (the elementary kernels), and open the feasability of massive PGCM emulations, the price being having first-handedly computed the costly elementary kernels. Limits of the emulator are also explored, by introducing the concept of over-training, which is exactly a consequence of the aproximativeness of a PGCM computation. Eventually this thesis demonstrates the possibility to emulate millions of PGCM computations with an error on collective spectroscopy less than 3%, and with a low numerical cost fraction of 1% of the million PGCM calculations cost.

•Oct 15, 2024•0 citations

Nuclear Structure from the In-Medium Similarity Renormalization Group

H. Hergert, J. M. Yao, Titus Morris +3 more

Efforts to describe nuclear structure and dynamics from first principles have advanced significantly in recent years. Exact methods for light nuclei are now able to include continuum degrees of freedom and treat structure and reactions on the same footing, and multiple approximate, computationally efficient many-body methods have been developed that can be routinely applied for medium-mass nuclei. This has made it possible to confront modern nuclear interactions from Chiral Effective Field Theory, that are rooted in Quantum Chromodynamics with a wealth of experimental data. Here, we discuss one of these efficient new many-body methods, the In-Medium Similarity Renormalization Group (IMSRG), and its applications in modern nuclear structure theory. The IMSRG evolves the nuclear many-body Hamiltonian in second-quantized form through continuous unitary transformations that can be implemented with polynomial computational effort. Through suitably chosen generators, we drive the matrix representation of the Hamiltonian in configuration space to specific shapes, e.g., to implement a decoupling of low-and high-energy scales, or to extract energy eigenvalues for a given nucleus. We present selected results from Multireference IMSRG (MR-IMSRG) calculations of open-shell nuclei, as well as proof-of-principle applications for intrinsically deformed medium-mass nuclei. We discuss the successes and prospects of merging the (MR-)IMSRG with many-body methods ranging from Configuration Interaction to the Density Matrix Renormalization Group, with the goal of achieving an efficient simultaneous description of dynamic and static correlationsin atomic nuclei.

arXiv:1805.09221•Journal of Physics Conference Series•Jun 1, 2018•19 citations

Nuclear Structure from the In-Medium Similarity Renormalization Group

H. Hergert, Morris T.D., Yao Jm +3 more

Efforts to describe nuclear structure and dynamics from first principles have advanced significantly in recent years. Exact methods for light nuclei are now able to include continuum degrees of freedom and treat structure and reactions on the same footing, and multiple approximate, computationally efficient many-body methods have been developed that can be routinely applied for medium-mass nuclei. This has made it possible to confront modern nuclear interactions from Chiral Effective Field Theory, that are rooted in Quantum Chromodynamics with a wealth of experimental data. Here, we discuss one of these efficient new many-body methods, the In-Medium Similarity Renormalization Group (IMSRG), and its applications in modern nuclear structure theory. The IMSRG evolves the nuclear many-body Hamiltonian in second-quantized form through continuous unitary transformations that can be implemented with polynomial computational effort. Through suitably chosen generators, we drive the matrix representation of the Hamiltonian in configuration space to specific shapes, e.g., to implement a decoupling of low-and high-energy scales, or to extract energy eigenvalues for a given nucleus. We present selected results from Multireference IMSRG (MR-IMSRG) calculations of open-shell nuclei, as well as proof-of-principle applications for intrinsically deformed medium-mass nuclei. We discuss the successes and prospects of merging the (MR-)IMSRG with many-body methods ranging from Configuration Interaction to the Density Matrix Renormalization Group, with the goal of achieving an efficient simultaneous description of dynamic and static correlationsin atomic nuclei.

•Apr 3, 2026•0 citations
In spite of missing dynamical correlations, the projected generator coordinate method (PGCM) was recently shown to be a suitable method to tackle the low-lying spectroscopy of complex nuclei. Still, describing absolute binding energies and reaching high accuracy eventually requires the inclusion of dynamical correlations on top of the PGCM. In this context, the present work discusses the first realistic results of a novel multi-reference perturbation theory (PGCM-PT) that can do so within a symmetry-conserving scheme for both ground and low-lying excited states. First, proof-of-principle calculations in a small (emax=4e_{\mathrm{max}}=4emax​=4) model space demonstrate that exact binding energies of closed- (\nucl{O}{16}) and open-shell (\nucl{O}{18}, \nucl{Ne}{20}) nuclei are reproduced within 0.5−1.5%0.5-1.5\%0.5−1.5% at second order, i.e. through PGCM-PT(2). Moreover, profiting from the pre-processing of the Hamiltonian via multi-reference in-medium similarity renormalization group transformations, PGCM-PT(2) can reach converged values within smaller model spaces than with an unevolved Hamiltonian. Doing so, dynamical correlations captured by PGCM-PT(2) are shown to bring essential corrections to low-lying excitation energies that become too dilated at leading order, i.e., at the strict PGCM level. The present work is laying the foundations for a better understanding of the optimal way to grasp static and dynamical correlations in a consistent fashion, with the aim of accurately describing ground and excited states of complex nuclei via ab initio many-body methods.
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.
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.
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.
An atomic nucleus is a quantum system of interacting nucleons and constitutes a problem difficult to solve exactly. For this reason, a diversity of approximate resolution methods has been designed, and Projected Generator Coordinate Method (PGCM) is one of them. The strong point of PGCM is to construct a physically inspired small dimensional space, in which an approximate solution of the nuclear many-body problem is easily found. However the numerical cost of PGCM space computation make this method inadapted for sensibility analysis of nuclear observables with restect to parametrisation of the interaction model, this analysis requiring an huge number of PGCM computations. In order to make this type of study possible, this thesis explore the concept of PGCM emulator. In this work, a combination of PGCM with Eigenvector Continuation (EC) is constructed and studied. This combination (the PGCM-EC emulator) takes advantage of mathematical similarities between PGCM and EC, and above all of the decomposition of the hamiltonian as a linear combination of parameter-independent terms. The latter property is used to concentrate the heavier numerical cost in the computation of parameter-independent quantities (the elementary kernels), and open the feasability of massive PGCM emulations, the price being having first-handedly computed the costly elementary kernels. Limits of the emulator are also explored, by introducing the concept of over-training, which is exactly a consequence of the aproximativeness of a PGCM computation. Eventually this thesis demonstrates the possibility to emulate millions of PGCM computations with an error on collective spectroscopy less than 3%, and with a low numerical cost fraction of 1% of the million PGCM calculations cost.
Efforts to describe nuclear structure and dynamics from first principles have advanced significantly in recent years. Exact methods for light nuclei are now able to include continuum degrees of freedom and treat structure and reactions on the same footing, and multiple approximate, computationally efficient many-body methods have been developed that can be routinely applied for medium-mass nuclei. This has made it possible to confront modern nuclear interactions from Chiral Effective Field Theory, that are rooted in Quantum Chromodynamics with a wealth of experimental data. Here, we discuss one of these efficient new many-body methods, the In-Medium Similarity Renormalization Group (IMSRG), and its applications in modern nuclear structure theory. The IMSRG evolves the nuclear many-body Hamiltonian in second-quantized form through continuous unitary transformations that can be implemented with polynomial computational effort. Through suitably chosen generators, we drive the matrix representation of the Hamiltonian in configuration space to specific shapes, e.g., to implement a decoupling of low-and high-energy scales, or to extract energy eigenvalues for a given nucleus. We present selected results from Multireference IMSRG (MR-IMSRG) calculations of open-shell nuclei, as well as proof-of-principle applications for intrinsically deformed medium-mass nuclei. We discuss the successes and prospects of merging the (MR-)IMSRG with many-body methods ranging from Configuration Interaction to the Density Matrix Renormalization Group, with the goal of achieving an efficient simultaneous description of dynamic and static correlationsin atomic nuclei.
Efforts to describe nuclear structure and dynamics from first principles have advanced significantly in recent years. Exact methods for light nuclei are now able to include continuum degrees of freedom and treat structure and reactions on the same footing, and multiple approximate, computationally efficient many-body methods have been developed that can be routinely applied for medium-mass nuclei. This has made it possible to confront modern nuclear interactions from Chiral Effective Field Theory, that are rooted in Quantum Chromodynamics with a wealth of experimental data. Here, we discuss one of these efficient new many-body methods, the In-Medium Similarity Renormalization Group (IMSRG), and its applications in modern nuclear structure theory. The IMSRG evolves the nuclear many-body Hamiltonian in second-quantized form through continuous unitary transformations that can be implemented with polynomial computational effort. Through suitably chosen generators, we drive the matrix representation of the Hamiltonian in configuration space to specific shapes, e.g., to implement a decoupling of low-and high-energy scales, or to extract energy eigenvalues for a given nucleus. We present selected results from Multireference IMSRG (MR-IMSRG) calculations of open-shell nuclei, as well as proof-of-principle applications for intrinsically deformed medium-mass nuclei. We discuss the successes and prospects of merging the (MR-)IMSRG with many-body methods ranging from Configuration Interaction to the Density Matrix Renormalization Group, with the goal of achieving an efficient simultaneous description of dynamic and static correlationsin atomic nuclei.