⚛️ nuclear experiments

Chiral three-nucleon forces for the new local position-space two-nucleon potential in $\textit{ab initio}$ many-body calculations

This paper constructs a chiral three-nucleon force tailored to the new local Idaho position-space two-nucleon potential, demonstrating that a hybrid local and nonlocal regulator improves the description of binding energies and radii for nuclei up to $^{132}$ Sn when the force's low-energy constants are constrained by $^3$ H and $^{16}$ O.

Original authors: Rongzhe Hu, Jianguo Li, Siqin Fan, Furong Xu

Published 2026-06-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Rongzhe Hu, Jianguo Li, Siqin Fan, Furong Xu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the atomic nucleus as a tiny, crowded dance floor where protons and neutrons (collectively called nucleons) are constantly spinning and interacting. For decades, physicists have tried to write the "rules of the dance" (the laws of physics) to predict exactly how these particles move, how tightly they hold hands (binding energy), and how big the dance floor is (charge radius).

For a long time, the rules worked well for predicting how tightly the dancers held hands, but they kept getting the size of the dance floor wrong. They always predicted the nucleus was too small.

This paper introduces a new set of rules that finally gets both the "grip" and the "size" right, at least for a wide range of nuclei from tiny Helium-4 up to the heavy Tin-132.

Here is the breakdown of their discovery using simple analogies:

1. The Problem: The "Two-Person" vs. "Three-Person" Dance

In the past, physicists mostly focused on how two dancers interact with each other (the Two-Nucleon Force). They built a very good map for this, specifically a new "local" map created by a team at Idaho. Think of this map as a set of instructions that says, "If you are standing here, you feel a push or pull from your neighbor right next to you."

However, the paper argues that you can't just look at pairs. Sometimes, three dancers interact all at once in a way that isn't just the sum of two pairs. This is the Three-Nucleon Force (3NF).

The Issue: Most existing rules for these three-person interactions were written in a different language (momentum space) than the new Idaho map (position space). It's like trying to use a GPS designed for driving a car to navigate a hiking trail; it's close, but not quite right.
The Goal: The authors wanted to build a new "Three-Person Rulebook" that speaks the exact same language as the new Idaho map.

2. The Solution: The "Hybrid" Regulator

To make the math work without blowing up the computer, physicists use "regulators." You can think of a regulator as a noise-canceling filter. It tells the math, "Ignore the super-fast, high-energy jitters that we can't measure, and focus on the smooth, slow movements we care about."

There are two main types of filters:

Local Filter: Only looks at what's happening right here.
Non-local Filter: Looks at what's happening here and a bit over there simultaneously.

The authors found that using only the "Local" filter (which matched the Idaho map) made the nucleus too small. It was like the dancers were huddling too tightly in the center.

The Breakthrough: They created a Hybrid Filter (Local-Nonlocal).

The Analogy: Imagine the dancers are wearing a mix of stiff, tight shoes (Local) and flexible, stretchy socks (Non-local). The stiff shoes keep them grounded, but the stretchy socks allow them to spread out just enough.
The Result: By mixing these two types of filters, the nucleus expanded to the correct size. The "hybrid" approach softened the interactions just enough to fix the size problem without breaking the energy calculations.

3. Tuning the Knobs (The Constants)

Every new set of rules has a few "knobs" or dials that need to be turned to match reality. In this paper, these are called Low-Energy Constants ( $c_D$ and $c_E$ ).

The authors didn't guess these numbers. They used a "calibration" method.
Step 1: They tuned the knobs so the math perfectly matched the energy of a tiny nucleus called Tritium (3H).
Step 2: They checked if those same knobs worked for a slightly larger nucleus, Oxygen-16.
Once the knobs were set using these two "test cases," they applied the rules to the rest of the periodic table.

4. The Results: A Perfect Fit

When they tested their new "Hybrid" rules against real-world data:

Binding Energy (The Grip): The calculations matched the experimental data almost perfectly. The nuclei were the right weight.
Charge Radius (The Size): This was the big win. Previous methods (using only local filters) consistently underestimated the size of the nucleus. The new Hybrid method got the size right for everything from light elements up to heavy Tin-132.
Density (The Crowd): They even looked at how the dancers were distributed inside the nucleus. The old method made the center too crowded. The new Hybrid method spread the density out correctly, matching what scientists actually see in experiments.

Summary

The paper claims that by building a specific "Three-Person Force" that is compatible with a new "Two-Person Force" map, and by using a hybrid filter (mixing local and non-local rules), they solved a long-standing puzzle. They finally have a set of rules that can accurately predict both how heavy a nucleus is and how big it is, all the way up to heavy elements like Tin.

This doesn't just fix a math problem; it gives physicists a more reliable "first-principles" tool to understand the fundamental structure of matter.

Technical Summary: Chiral Three-Nucleon Forces for the New Local Position-Space Two-Nucleon Potential

Problem Statement
A primary objective in ab initio nuclear theory is to describe atomic nuclei and nuclear matter using fundamental degrees of freedom derived from chiral effective field theory ( $\chi$ EFT). While two-nucleon forces (2NF) and three-nucleon forces (3NF) naturally emerge within this framework, a persistent challenge remains: developing a unified chiral Hamiltonian that simultaneously reproduces the experimental ground-state energies and charge radii of medium-mass nuclei, as well as the saturation properties of nuclear matter. Previous studies utilizing purely local regulators for 3NFs have often resulted in a systematic underestimation of charge radii, while nonlocal momentum-space potentials have been difficult to apply to heavier nuclei due to computational constraints in Quantum Monte Carlo (QMC) methods. Furthermore, the Idaho group recently developed a new local position-space 2NF, but a companion chiral 3NF tailored specifically to this potential was lacking.

Methodology
The authors constructed a companion chiral 3NF at next-to-next-to-leading order (N2LO) specifically designed to work with the new Idaho local position-space 2NF (truncated at N3LO). The 2NF was evolved to a low-momentum scale ( $\lambda = 2.2 \text{ fm}^{-1}$ ) using the Similarity Renormalization Group (SRG) to accelerate convergence in many-body calculations.

The study focused on comparing two regularization schemes for the 3NF:

Purely Local Regulator: Using a local form function.
Hybrid Local-Nonlocal (lnl) Regulator: A product of local and nonlocal form functions.

The low-energy constants (LECs), $c_D$ and $c_E$ , which parameterize the intermediate-range one-pion exchange and short-range contact terms of the 3NF, were constrained using the ground-state energies of $^3$ H and $^{16}$ O. The $^3$ H energy was calculated using the No-Core Shell Model (NCSM) in Jacobi coordinates, while $^{16}$ O was calculated using the In-Medium Similarity Renormalization Group (IMSRG) in the Magnus formulation. The authors accounted for theoretical uncertainties from $\chi$ EFT truncation (estimated at $\sim$ 300 keV for $^3$ H) and many-body approximations (estimated at $\sim$ 2% for IMSRG).

The resulting Hamiltonians were applied to calculate binding energies and charge radii for a wide range of nuclei, from $^4$ He to $^{132}$ Sn. Closed-shell nuclei were treated with single-reference IMSRG, while open-shell isotopic chains (O, Mg, Si, S, Ar, Ca) were calculated using the valence-space IMSRG (VS-IMSRG) with ensemble normal ordering (ENO) to handle open-shell orbitals. Charge density distributions for $^{16}$ O and $^{40}$ Ca were also analyzed.

Key Contributions

Construction of a Companion 3NF: The paper provides the first chiral 3NF specifically tailored to the new Idaho local position-space 2NF.
Systematic Comparison of Regulators: It offers a systematic comparison between local and hybrid local-nonlocal (lnl) 3NF regulators within the context of the new local 2NF.
Optimization of LECs: The study determines optimal values for $c_D$ and $c_E$ for both regulator types, constrained by $^3$ H and $^{16}$ O, demonstrating that the optimal parameters differ significantly depending on the regulator scheme (e.g., $c_D = -2.0$ for local vs. $c_D = 6.0$ for lnl).

Results

Binding Energies: Both the local and lnl 3NF families, when combined with the new local 2NF, provide good descriptions of ground-state energies for medium-mass closed-shell nuclei and isotopic chains.
Charge Radii: A distinct difference was observed in the description of charge radii. The purely local 3NF combined with the local 2NF systematically underestimated charge radii. In contrast, the hybrid lnl 3NF regulator significantly improved the description of charge radii, bringing deviations from experimental data below 2% for nuclei ranging from $^4$ He to $^{132}$ Sn.
Charge Density Distributions: Calculations for $^{16}$ O and $^{40}$ Ca showed that the lnl regulator yields charge density distributions in excellent agreement with experimental data. The local 3NF resulted in densities that were significantly larger in the central region of the nucleus. The authors attribute this to the fact that local regulators are "harder," inducing stronger short-range correlations, whereas the nonlocal component in the lnl regulator softens these correlations, increasing the exterior density and the overall radius.
Magnitude of Improvement: The transition from a local to a local-nonlocal regulator resulted in approximately a 7% increase in the charge radius for medium-mass nuclei, correcting the systematic underestimation seen with purely local schemes.

Significance
The paper concludes that the regularization scheme of the chiral 3NF is a critical factor influencing the quality of ab initio many-body predictions, particularly for nuclear radii. The developed chiral Hamiltonian, combining the new local position-space 2NF with the hybrid lnl 3NF, offers a unified framework capable of accurately describing both binding energies and charge radii across a broad range of the nuclear chart. The authors assert that this Hamiltonian should be a valuable tool for future high-precision ab initio nuclear structure calculations.

1. The Problem: The "Two-Person" vs. "Three-Person" Dance

2. The Solution: The "Hybrid" Regulator

3. Tuning the Knobs (The Constants)

4. The Results: A Perfect Fit

Summary

More like this