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

Irreducible 3-body forces contributions to the self-energy

Raimondi, F., C. Barbieri

Published Jan 27, 2017• 0 citations
arXivHTMLPDF

Abstract

The inclusion of the three-nucleon forces (3NFs) in \textit{ab initio} many-body approaches is a formidable task, due to the computational load implied by the treatment of their matrix elements. For this reason, practical applications have mostly been limited to contributions where 3NFs enter as effective two-nucleon interactions. In this contribution, we derive the algebraic diagrammatic construction (ADC) working equations for a specific Feynman diagram of the self-energy that contains a fully irreducible three-nucleon force. This diagram is expected to be the most important among those previously neglected, because it connects dominant excited intermediate state configurations.

DOI
10.48550/arxiv.1701.08127

Related papers

Algebraic diagrammatic construction formalism with three-body interactions

Francesco Raimondi, C. Barbieri

Background: Self-consistent Green's function theory has recently been extended to the basic formalism needed to account for three-body interactions [Carbone, Cipollone, Barbieri, Rios, and Polls, Phys. Rev. C 88, 054326 (2013)]. The contribution of three-nucleon forces has so far been included in ab initio calculations on nuclear matter and finite nuclei only as averaged two-nucleon forces.Purpose: We derive the working equations for all possible two- and three-nucleon terms that enter the expansion of the self-energy up to the third order, thus including the interaction-irreducible (i.e., not averaged) diagrams with three-nucleon forces that have been previously neglected.Methods: We employ the algebraic diagrammatic construction up to the third order as an organization scheme for generating a nonperturbative self-energy, in which ring (particle-hole) and ladder (particle-particle) diagrams are resummed to all orders.Results: We derive expressions of the static and dynamic self-energy up to the third order, by taking into account the set of diagrams required when either the skeleton or nonskeleton expansions of the single-particle propagator are assumed. A hierarchy of importance among different diagrams is revealed, and a particular emphasis is given to a third-order diagram [see Fig. 2(c)] that is expected to play a significant role among those featuring an interaction-irreducible three-nucleon force.Conclusion: A consistent formalism to resum at infinite order correlations induced by three-nucleon forces in the self-consistent Green's function theory is now available and ready to be implemented in the many-body solvers.

arXiv:1709.04330•Physical Review C•May 8, 2018•37 citations

Gorkov algebraic diagrammatic construction formalism at third order

C. Barbieri, T. Duguet, V. Somà

Background: The Gorkov approach to self-consistent Green's function theory has been formulated by Som\`a, Duguet, and Barbieri in [Phys. Rev. C 84, 064317 (2011)]. Over the past decade, it has become a method of reference for first-principles computations of semimagic nuclear isotopes. The currently available implementation is limited to a second-order self-energy and neglects particle-number nonconserving terms arising from contracting three-particle forces with anomalous propagators. For nuclear physics applications, this is sufficient to address first-order energy differences (i.e., two neutron separation energies, excitation energies of states dominating the one-nucleon spectral function), ground-state radii and moments on an accurate enough basis. However, addressing absolute binding energies, fine spectroscopic details of $N\ifmmode\pm\else\textpm\fi{}1$ particle systems or delicate quantities such as second-order energy differences associated with pairing gaps, requires going to higher truncation orders.Purpose: The formalism is extended to third order in the algebraic diagrammatic construction (ADC) expansion with two-body Hamiltonians.Methods: The expansion of Gorkov propagators in Feynman diagrams is combined with the algebraic diagrammatic construction up to the third order as an organization scheme to generate the Gorkov self-energy.Results: Algebraic expressions for the static and dynamic contributions to the self-energy, along with equations for the matrix elements of the Gorkov eigenvalue problem, are derived. It is first done for a general basis before specifying the set of equations to the case of spherical systems displaying rotational symmetry. Workable approximations to the full self-consistency problem are also elaborated on. The formalism at third order it thus complete for a general two-body Hamiltonian.Conclusions: Working equations for the full Gorkov-ADC(3) are now available for numerical implementation.

arXiv:2112.08322•Physical Review C•Apr 28, 2022•21 citations

Exact calculation of three-body contact interaction to second order1

Norbert Kaiser

For a system of fermions with a three-body contact interaction the second-order contributions to the energy per particle $\ensuremath \bar E(k_f)$ are calculated exactly. The three-particle scattering amplitude in the medium is derived in closed analytical form from the corresponding two-loop rescattering diagram. We compare the (genuine) second-order three-body contribution to $\ensuremath \bar E(k_f)\sim k_f^{10}$ with the second-order term due to the density-dependent effective two-body interaction, and find that the latter term dominates. The results of the present study are of interest for nuclear many-body calculations where chiral three-nucleon forces are treated beyond leading order via a density-dependent effective two-body interaction.

•Aug 16, 2016•11 citations

Closed-form intermediate representations of many-body propagators and resolvent matrices

J. Schirmer

A new type of ``intermediate'' state, \ensuremath{\Vert}\ensuremath{\Psi}${\mathrm{\ifmmode \tilde{}\else \~{}\fi{}}}_{\mathit{J}}$〉, is introduced in the representation of many-body propagators and related resolvent matrix elements. A simple algebraic procedure is presented for constructing the unitary transformation matrix Q that relates the intermediate states and the exact energy eigenstates, \ensuremath{\Vert}${\mathrm{\ensuremath{\Psi}}}_{\mathit{n}}$〉, of the interacting system. Here the starting point is the matrix X of generalized spectroscopic amplitudes 〈${\mathrm{\ensuremath{\Psi}}}_{\mathit{n}}$\ensuremath{\Vert}C${\mathrm{^}}_{\mathit{J}}$\ensuremath{\Vert}${\mathrm{\ensuremath{\Psi}}}_{0}^{\mathit{N}}$〉, where C${\mathrm{^}}_{\mathit{J}}$ denotes a physical excitation operator and \ensuremath{\Vert}${\mathrm{\ensuremath{\Psi}}}_{0}^{\mathit{N}}$〉 is the N-electron ground state. A block QR decomposition [G. H. Golub and C. F. von Loan, MatrixB Computations (Johns Hopkins University, Baltimore, 1989)] of X according to the equation X=Q F allows one to determine explicit expressions for the subblocks of Q and the intermediate propagator representations constituted by a nondiagonal effective interaction matrix C and an effective spectroscopic matrix f. These effective quantities f and C are formulated entirely in terms of ground-state density-matrix elements and related energy expectation values. The relevance of the intermediate representations as a means for deriving computational schemes is based on the regularity and compactness of the perturbation expansions for the effective matrices f and C. These basic properties also establish that the intermediate representations are closed-form versions of the algebraic-diagrammatic construction approximation schemes derived previously as a reformulation of the original diagrammatic propagator perturbation series. The intermediate representations represent a link between algebraic and diagrammatic approaches in the field of propagator methods and are expected to be useful also in the development of nonperturbative approximations.

Physical Review A•May 1, 1991•227 citations

Coupled-cluster theory for three-body Hamiltonians

G. Hagen, T. Papenbrock, David J. Dean +4 more

We derive coupled-cluster equations for three-body Hamiltonians. The equations for the one- and two-body cluster amplitudes are presented in a factorized form that leads to an efficient numerical implementation. We employ low-momentum two- and three-nucleon interactions and calculate the binding energy of $^{4}\mathrm{He}$. The results show that the main contribution of the three-nucleon interaction stems from its density-dependent zero-, one-, and two-body terms that result from the normal ordering of the Hamiltonian in coupled-cluster theory. The residual three-body terms that remain after normal ordering can be neglected.

arXiv:0704.2854•Physical Review C•Sep 4, 2007•191 citations

Exactness of the normal-ordered two-body truncation of three-nucleon forces

Maxwell Rothman, Ben Johnson-Toth, Francesca Bonaiti +3 more

Reference-state-based many-body methods start from Hamiltonians that are normal ordered with respect to the reference state. In low-energy nuclear physics applications, normal-ordered Hamiltonians consisting of two- and three-nucleon forces are usually truncated at the two-body rank with residual three-nucleon operators being discarded. Benchmark computations have shown that this truncation is accurate, but we lack an understanding about why it works. We show that the normal-ordered two-body truncation is exact for zero-range three-body forces when nuclei are computed using the coupled cluster with singles and doubles method. As the nuclear three-nucleon force is short ranged and a three-body contact is a leading term in effective field theories of quantum chromodynamics, our result provides an analytical basis for the popular normal-ordered two-body approximation.

arXiv:2508.01507•Physical Review C•Sep 30, 2025•6 citations
The inclusion of the three-nucleon forces (3NFs) in \textit{ab initio} many-body approaches is a formidable task, due to the computational load implied by the treatment of their matrix elements. For this reason, practical applications have mostly been limited to contributions where 3NFs enter as effective two-nucleon interactions. In this contribution, we derive the algebraic diagrammatic construction (ADC) working equations for a specific Feynman diagram of the self-energy that contains a fully irreducible three-nucleon force. This diagram is expected to be the most important among those previously neglected, because it connects dominant excited intermediate state configurations.
Background: Self-consistent Green's function theory has recently been extended to the basic formalism needed to account for three-body interactions [Carbone, Cipollone, Barbieri, Rios, and Polls, Phys. Rev. C 88, 054326 (2013)]. The contribution of three-nucleon forces has so far been included in ab initio calculations on nuclear matter and finite nuclei only as averaged two-nucleon forces.Purpose: We derive the working equations for all possible two- and three-nucleon terms that enter the expansion of the self-energy up to the third order, thus including the interaction-irreducible (i.e., not averaged) diagrams with three-nucleon forces that have been previously neglected.Methods: We employ the algebraic diagrammatic construction up to the third order as an organization scheme for generating a nonperturbative self-energy, in which ring (particle-hole) and ladder (particle-particle) diagrams are resummed to all orders.Results: We derive expressions of the static and dynamic self-energy up to the third order, by taking into account the set of diagrams required when either the skeleton or nonskeleton expansions of the single-particle propagator are assumed. A hierarchy of importance among different diagrams is revealed, and a particular emphasis is given to a third-order diagram [see Fig. 2(c)] that is expected to play a significant role among those featuring an interaction-irreducible three-nucleon force.Conclusion: A consistent formalism to resum at infinite order correlations induced by three-nucleon forces in the self-consistent Green's function theory is now available and ready to be implemented in the many-body solvers.
Background: The Gorkov approach to self-consistent Green's function theory has been formulated by Som\`a, Duguet, and Barbieri in [Phys. Rev. C 84, 064317 (2011)]. Over the past decade, it has become a method of reference for first-principles computations of semimagic nuclear isotopes. The currently available implementation is limited to a second-order self-energy and neglects particle-number nonconserving terms arising from contracting three-particle forces with anomalous propagators. For nuclear physics applications, this is sufficient to address first-order energy differences (i.e., two neutron separation energies, excitation energies of states dominating the one-nucleon spectral function), ground-state radii and moments on an accurate enough basis. However, addressing absolute binding energies, fine spectroscopic details of $N\ifmmode\pm\else\textpm\fi{}1$ particle systems or delicate quantities such as second-order energy differences associated with pairing gaps, requires going to higher truncation orders.Purpose: The formalism is extended to third order in the algebraic diagrammatic construction (ADC) expansion with two-body Hamiltonians.Methods: The expansion of Gorkov propagators in Feynman diagrams is combined with the algebraic diagrammatic construction up to the third order as an organization scheme to generate the Gorkov self-energy.Results: Algebraic expressions for the static and dynamic contributions to the self-energy, along with equations for the matrix elements of the Gorkov eigenvalue problem, are derived. It is first done for a general basis before specifying the set of equations to the case of spherical systems displaying rotational symmetry. Workable approximations to the full self-consistency problem are also elaborated on. The formalism at third order it thus complete for a general two-body Hamiltonian.Conclusions: Working equations for the full Gorkov-ADC(3) are now available for numerical implementation.
For a system of fermions with a three-body contact interaction the second-order contributions to the energy per particle $\ensuremath \bar E(k_f)$ are calculated exactly. The three-particle scattering amplitude in the medium is derived in closed analytical form from the corresponding two-loop rescattering diagram. We compare the (genuine) second-order three-body contribution to $\ensuremath \bar E(k_f)\sim k_f^{10}$ with the second-order term due to the density-dependent effective two-body interaction, and find that the latter term dominates. The results of the present study are of interest for nuclear many-body calculations where chiral three-nucleon forces are treated beyond leading order via a density-dependent effective two-body interaction.
A new type of ``intermediate'' state, \ensuremath{\Vert}\ensuremath{\Psi}${\mathrm{\ifmmode \tilde{}\else \~{}\fi{}}}_{\mathit{J}}$〉, is introduced in the representation of many-body propagators and related resolvent matrix elements. A simple algebraic procedure is presented for constructing the unitary transformation matrix Q that relates the intermediate states and the exact energy eigenstates, \ensuremath{\Vert}${\mathrm{\ensuremath{\Psi}}}_{\mathit{n}}$〉, of the interacting system. Here the starting point is the matrix X of generalized spectroscopic amplitudes 〈${\mathrm{\ensuremath{\Psi}}}_{\mathit{n}}$\ensuremath{\Vert}C${\mathrm{^}}_{\mathit{J}}$\ensuremath{\Vert}${\mathrm{\ensuremath{\Psi}}}_{0}^{\mathit{N}}$〉, where C${\mathrm{^}}_{\mathit{J}}$ denotes a physical excitation operator and \ensuremath{\Vert}${\mathrm{\ensuremath{\Psi}}}_{0}^{\mathit{N}}$〉 is the N-electron ground state. A block QR decomposition [G. H. Golub and C. F. von Loan, MatrixB Computations (Johns Hopkins University, Baltimore, 1989)] of X according to the equation X=Q F allows one to determine explicit expressions for the subblocks of Q and the intermediate propagator representations constituted by a nondiagonal effective interaction matrix C and an effective spectroscopic matrix f. These effective quantities f and C are formulated entirely in terms of ground-state density-matrix elements and related energy expectation values. The relevance of the intermediate representations as a means for deriving computational schemes is based on the regularity and compactness of the perturbation expansions for the effective matrices f and C. These basic properties also establish that the intermediate representations are closed-form versions of the algebraic-diagrammatic construction approximation schemes derived previously as a reformulation of the original diagrammatic propagator perturbation series. The intermediate representations represent a link between algebraic and diagrammatic approaches in the field of propagator methods and are expected to be useful also in the development of nonperturbative approximations.
We derive coupled-cluster equations for three-body Hamiltonians. The equations for the one- and two-body cluster amplitudes are presented in a factorized form that leads to an efficient numerical implementation. We employ low-momentum two- and three-nucleon interactions and calculate the binding energy of 4He^{4}\mathrm{He}4He. The results show that the main contribution of the three-nucleon interaction stems from its density-dependent zero-, one-, and two-body terms that result from the normal ordering of the Hamiltonian in coupled-cluster theory. The residual three-body terms that remain after normal ordering can be neglected.
Reference-state-based many-body methods start from Hamiltonians that are normal ordered with respect to the reference state. In low-energy nuclear physics applications, normal-ordered Hamiltonians consisting of two- and three-nucleon forces are usually truncated at the two-body rank with residual three-nucleon operators being discarded. Benchmark computations have shown that this truncation is accurate, but we lack an understanding about why it works. We show that the normal-ordered two-body truncation is exact for zero-range three-body forces when nuclei are computed using the coupled cluster with singles and doubles method. As the nuclear three-nucleon force is short ranged and a three-body contact is a leading term in effective field theories of quantum chromodynamics, our result provides an analytical basis for the popular normal-ordered two-body approximation.