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:2508.01507 · Physical Review C

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

Maxwell Rothman, Ben Johnson-Toth, Francesca Bonaiti, G. Hagen, M. Heinz, T. Papenbrock

Published Sep 30, 2025• 6 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/psnh-gc5k

Related papers

Normal-ordered k-body approximation in particle-number-breaking theories

Ripoche, J., A. Tichai, T. Duguet

The reach of ab initio many-body theories is rapidly extending over the nuclear chart. However, dealing fully with three-nucleon, possibly four-nucleon, interactions makes the solving of the A-body Schrodinger equation particularly cumbersome, if not impossible beyond a certain nuclear mass. Consequently, ab initio calculations of mid-mass nuclei are typically performed on the basis of the so-called normal-ordered two-body (NO2B) approximation that captures dominant effects of three-nucleon forces while effectively working with two-nucleon operators. A powerful idea currently employed to extend ab initio calculations to open-shell nuclei consists of expanding the exact solution of the A-body Schrodinger equation while authorizing the approximate solution to break symmetries of the Hamiltonian. In this context, operators are normal ordered with respect to a symmetry-breaking reference state such that proceeding to a naive truncation may lead to symmetry-breaking approximate operators. The purpose of the present work is to design a normal-ordering approximation of operators that is consistent with the symmetries of the Hamiltonian while working in the context of symmetry broken (and potentially restored) methods. Focusing on many-body formalisms in which U(1) global-gauge symmetry associated with particle-number conservation is broken (and potentially restored), a particle-number-conserving normal-orderedk-body (PNOkB) approximation of an arbitrary N-body operator is designed on the basis of Bogoliubov reference states. A numerical test based on particle-number-projected Hartree-Fock-Bogoliubov calculations permits to check the particle-number conserving/violating character of a given approximation to a particle-number conserving operator. The PNOkB approximation of an arbitrary N-body operator is formulated. Based on this systematic approach, it is demonstrated that naive extensions of the normal-ordered two-body (NO2B) approximation employed so far on the basis of symmetry-conserving reference states lead to particle non-conserving operators. Alternatively, the PNOkB procedure is now available to generate particle-number-conserving approximate operators. The formal analysis is validated numerically. Using the presently proposed PNOkB approximation, ab initio calculations based on symmetry-breaking and restored formalisms can be safely performed. The future formulation of an angular-momentum-conserving normal-ordered k-body approximation based on deformed Slater determinant or Bogoliubov reference states is envisioned.

•Feb 1, 2020•20 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

Normal-ordering approximations and translational (non)invariance

T. Djärv, A. Ekström, C. Forssén +1 more

Normal ordering provides an approach to approximate three-body forces as effective two-body operators and it is therefore an important tool in many-body calculations with realistic nuclear interactions. The corresponding neglect of certain three-body terms in the normal-ordered Hamiltonian is known to influence translational invariance, although the magnitude of this effect has not yet been systematically quantified. In this paper we study in particular the normal-ordering two-body approximation applied to a single harmonic-oscillator reference state. We explicate the breaking of translational invariance and demonstrate the magnitude of the approximation error as a function of model space parameters for $^{4}\mathrm{He}$ and $^{16}\mathrm{O}$ by performing full no-core shell-model calculations with and without three-nucleon forces. We combine two different diagnostics to better monitor the breaking of translational invariance. While the center-of-mass effect is shown to become potentially very large for $^{4}\mathrm{He}$, it is also shown to be much smaller for $^{16}\mathrm{O}$ although full convergence is not reached. These tools can be easily implemented in studies using other many-body frameworks and bases.

arXiv:2104.14820•Physical Review C•Aug 23, 2021•11 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

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

Short-Ranged Central and Tensor Correlations in Nuclear Many-Body Systems

Thomas Neff

In this work a unitary correlation operator is presented that explicitly describes the short-ranged central and tensor correlations in the nuclear many-body system. These short-ranged correlations are induced by the repulsive core and the pronounced tensor force of realistic nucleon-nucleon interactions and cannot be described by the simple many-body states of a mean-field or shell model approach. The unitary correlation operator is discussed for the Argonne V18 and the Bonn-A interactions. Applying the correlation operator onto the Hamiltonian a common effective interaction for low momenta is obtained. Calculations for 4He using the one- and two-body part of the correlated Hamiltonian compare favorably with exact many-body methods. Calculations for 16O and 40Ca which are not possible with exact methods are performed using harmonic oscillator shell model states. The observed deviation from the experimental binding energies and radii is attributed to the missing three-body forces. The correlated interaction in a basis-free representation is used in the Fermionic Molecular Dynamics model to calculate the nuclei of the p- and the sd-shell. In this model the uncorrelated many-particle state is given as a Slaterdeterminant of Gaussian wave-packets with spin and isospin degrees of freedom which allows a consistent description of spherically symmetric, intrinsically deformed and clustered nuclei.

•Jan 1, 2008•5 citations
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.
The reach of ab initio many-body theories is rapidly extending over the nuclear chart. However, dealing fully with three-nucleon, possibly four-nucleon, interactions makes the solving of the A-body Schrodinger equation particularly cumbersome, if not impossible beyond a certain nuclear mass. Consequently, ab initio calculations of mid-mass nuclei are typically performed on the basis of the so-called normal-ordered two-body (NO2B) approximation that captures dominant effects of three-nucleon forces while effectively working with two-nucleon operators. A powerful idea currently employed to extend ab initio calculations to open-shell nuclei consists of expanding the exact solution of the A-body Schrodinger equation while authorizing the approximate solution to break symmetries of the Hamiltonian. In this context, operators are normal ordered with respect to a symmetry-breaking reference state such that proceeding to a naive truncation may lead to symmetry-breaking approximate operators. The purpose of the present work is to design a normal-ordering approximation of operators that is consistent with the symmetries of the Hamiltonian while working in the context of symmetry broken (and potentially restored) methods. Focusing on many-body formalisms in which U(1) global-gauge symmetry associated with particle-number conservation is broken (and potentially restored), a particle-number-conserving normal-orderedk-body (PNOkB) approximation of an arbitrary N-body operator is designed on the basis of Bogoliubov reference states. A numerical test based on particle-number-projected Hartree-Fock-Bogoliubov calculations permits to check the particle-number conserving/violating character of a given approximation to a particle-number conserving operator. The PNOkB approximation of an arbitrary N-body operator is formulated. Based on this systematic approach, it is demonstrated that naive extensions of the normal-ordered two-body (NO2B) approximation employed so far on the basis of symmetry-conserving reference states lead to particle non-conserving operators. Alternatively, the PNOkB procedure is now available to generate particle-number-conserving approximate operators. The formal analysis is validated numerically. Using the presently proposed PNOkB approximation, ab initio calculations based on symmetry-breaking and restored formalisms can be safely performed. The future formulation of an angular-momentum-conserving normal-ordered k-body approximation based on deformed Slater determinant or Bogoliubov reference states is envisioned.
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.
Normal ordering provides an approach to approximate three-body forces as effective two-body operators and it is therefore an important tool in many-body calculations with realistic nuclear interactions. The corresponding neglect of certain three-body terms in the normal-ordered Hamiltonian is known to influence translational invariance, although the magnitude of this effect has not yet been systematically quantified. In this paper we study in particular the normal-ordering two-body approximation applied to a single harmonic-oscillator reference state. We explicate the breaking of translational invariance and demonstrate the magnitude of the approximation error as a function of model space parameters for 4He^{4}\mathrm{He}4He and 16O^{16}\mathrm{O}16O by performing full no-core shell-model calculations with and without three-nucleon forces. We combine two different diagnostics to better monitor the breaking of translational invariance. While the center-of-mass effect is shown to become potentially very large for 4He^{4}\mathrm{He}4He, it is also shown to be much smaller for 16O^{16}\mathrm{O}16O although full convergence is not reached. These tools can be easily implemented in studies using other many-body frameworks and bases.
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.
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.
In this work a unitary correlation operator is presented that explicitly describes the short-ranged central and tensor correlations in the nuclear many-body system. These short-ranged correlations are induced by the repulsive core and the pronounced tensor force of realistic nucleon-nucleon interactions and cannot be described by the simple many-body states of a mean-field or shell model approach. The unitary correlation operator is discussed for the Argonne V18 and the Bonn-A interactions. Applying the correlation operator onto the Hamiltonian a common effective interaction for low momenta is obtained. Calculations for 4He using the one- and two-body part of the correlated Hamiltonian compare favorably with exact many-body methods. Calculations for 16O and 40Ca which are not possible with exact methods are performed using harmonic oscillator shell model states. The observed deviation from the experimental binding energies and radii is attributed to the missing three-body forces. The correlated interaction in a basis-free representation is used in the Fermionic Molecular Dynamics model to calculate the nuclei of the p- and the sd-shell. In this model the uncorrelated many-particle state is given as a Slaterdeterminant of Gaussian wave-packets with spin and isospin degrees of freedom which allows a consistent description of spherically symmetric, intrinsically deformed and clustered nuclei.