Ab initio calculations of nuclear charge radii across and beyond ${}^{132}$Sn: Putting chiral EFT nuclear interactions to the test

This study employs ab initio Bogoliubov coupled cluster calculations to demonstrate that current chiral effective field theory interactions fail to simultaneously reproduce the absolute charge radii, parabolic isotopic shifts, and the kink at ${}^{132}$Sn in the Tin chain, thereby highlighting the need for improved theoretical models and new experimental data to better constrain nuclear forces.

The Gist

Imagine the atomic nucleus as a tiny, bustling city made of protons and neutrons. The "charge radius" of this city is simply a measure of how big the city is. Scientists have long tried to build a perfect map of these cities using a set of rules called Chiral Effective Field Theory (χEFT). Think of these rules as a "Grand Blueprint" that tells the protons and neutrons how to behave, how to stick together, and how far apart they should stand.

For a long time, this Blueprint worked okay for small cities (light elements), but when scientists tried to use it for the massive, heavy city of Tin (Sn), things got messy. Specifically, the Blueprint struggled to predict the size of the city as you added more and more neutrons (the "citizens" that don't have an electric charge).

Here is a simple breakdown of what this paper did and what it found:

1. The Experiment: Building a Better Map

The researchers used a super-advanced computer simulation (called Bogoliubov Coupled Cluster) to build a 3D model of Tin isotopes (different versions of the Tin city with varying numbers of neutrons). They tested three different versions of the "Grand Blueprint" (three different sets of interaction rules) to see which one could best predict the actual size of these cities.

They focused on a specific stretch of the periodic table: from Tin-96 up to Tin-150. This range is special because it crosses two "magic numbers" (50 and 82), which are like major city borders where the population structure changes dramatically.

2. The "Kink" in the Road

When you look at how the size of the Tin city changes as you add neutrons, it usually follows a smooth, curved path (like a parabola). However, at a specific point—Tin-132—the road suddenly bends sharply. This is called a "kink."

• The Challenge: The scientists wanted to see if their Blueprints could predict this sharp bend.
• The Result:
  • Blueprint A (1.8/2.0 EM): This one was too conservative. It predicted the city would be too small and missed the sharp bend entirely.
  • Blueprint B (ΔNNLOGO): This was better at getting the general size right, but it still missed the sharp bend at Tin-132.
  • Blueprint C (1.8/2.0 EM7.5): This one was a "miracle worker" for Tin-132! It predicted the sharp bend perfectly, matching real-world measurements.

3. The Twist: A "Fix" That Broke Something Else

Here is the plot twist. While Blueprint C got the bend at Tin-132 right, it did so for the wrong reasons.

Imagine you are trying to predict the traffic flow in a city. Blueprint C predicted the traffic jam at the border correctly, but only because it assumed the streets were arranged in a completely different way than they actually are.

• The Problem: Beyond Tin-132 (specifically around Tin-142), Blueprint C started predicting weird, inverted shapes in the city's growth that don't make sense physically. It predicted a "reverse kink" that likely won't happen in reality.
• The Lesson: Just because a model gets one specific number right doesn't mean the model is perfect. Blueprint C got the size of Tin-132 right by accidentally compensating for other errors, but it failed to describe the underlying structure of the city correctly.

4. The "Magic" of Neutron Shells

To understand why this happened, the scientists looked at the "neighborhoods" inside the nucleus (called shells).

• In the real world, neutrons fill up specific neighborhoods in a specific order.
• The successful Blueprint (C) assumed neutrons were moving into a "spacious" neighborhood (the $1h_{9/2}$ shell) right after the city border at 132. This spacious neighborhood pulled the city walls outward, creating the correct "kink."
• However, other evidence suggests neutrons should actually be moving into a "tighter" neighborhood. If they move into the tight one, the city shouldn't grow as fast, and the kink wouldn't be as sharp.

5. What's Next?

The paper concludes that while we have made great progress, our "Grand Blueprint" isn't finished yet.

• We need more data: We need to measure the sizes of Tin cities even further out (beyond Tin-134) and closer to the beginning (near Tin-100) to see which Blueprint is truly correct.
• We need better math: The current computer models are very good, but they are missing some "fine details" (called "triples" corrections) that might be necessary to get the physics 100% right.

The Big Picture Analogy

Think of the atomic nucleus like a balloon.

• The Blueprints are instructions on how much air to pump in as you add more people (neutrons) inside.
• The Kink is a sudden, sharp bulge in the balloon at a specific point.
• The Study found that one set of instructions (Blueprint C) made the balloon bulge at the right spot, but it did so by stretching the rubber in a way that would cause the balloon to pop or twist strangely a few inches later.
• The Goal: We need to find the set of instructions that makes the balloon bulge correctly and keeps the rest of the balloon smooth and logical.

In short: This paper is a stress test for our best theories of nuclear physics. It shows that while we are getting closer to the truth, our current "rules" still have flaws, and the heavy, neutron-rich versions of Tin are the perfect testing ground to fix them.

Technical Summary

1. Problem Statement

The paper addresses the persistent challenge in nuclear physics of accurately reproducing nuclear charge radii using ab initio methods based on Chiral Effective Field Theory ($\chi$EFT). While modern interactions can successfully predict binding energies, they often fail to simultaneously reproduce absolute charge radii and their evolution across isotopic chains.

• Specific Context: The study focuses on the Tin (Sn) isotopic chain ($Z=50$), a semi-magic system spanning from neutron-deficient $^{96}$Sn to neutron-rich $^{150}$Sn.
• Key Phenomena: The study targets specific structural features that serve as stringent tests for nuclear forces:
  • The inverted parabolic behavior of isotopic shifts between magic numbers $N=50$ and $N=82$.
  • The pronounced "kink" (sharp increase) in charge radii at the doubly magic nucleus $^{132}$Sn.
  • The behavior of radii beyond $N=82$ (e.g., $^{134}$Sn and $^{142}$Sn), where recent experimental data has revealed significant rises.
• Goal: To evaluate whether state-of-the-art fine-tuned $\chi$EFT Hamiltonians can simultaneously reproduce binding energies, absolute radii, and the complex evolution of isotopic shifts across this heavy mass region.

2. Methodology

The authors employ Bogoliubov Coupled Cluster (BCC) theory, an ab initio many-body method suitable for open-shell nuclei.

• Many-Body Framework:

  • Reference State: A particle-number-breaking Bogoliubov vacuum ($|\Phi\rangle$) obtained from Hartree-Fock-Bogoliubov (HFB) mean-field equations, enforcing rotational invariance ($J=0$).
  • Cluster Operator: The ground state is parameterized as $|\Psi_0\rangle = e^T |\Phi\rangle$. The calculations are truncated at the Singles and Doubles (BCCSD) level ($T = T_1 + T_2$).
  • Observable Evaluation: Expectation values are computed using a linear-response framework with a de-excitation operator $\Lambda$. In this work, $\Lambda$ is approximated to first order ($\Lambda \approx T^\dagger$).
  • Justification: While BCCSD has a residual error of $\sim$10% for binding energies (due to missing triples), it is expected to be highly accurate for charge radii, as correlation effects are less dominant for radii than for energies.

• Hamiltonians Tested: Three distinct fine-tuned $\chi$EFT-based interactions were employed:

  1. 1.8/2.0 (EM): Accurate for binding energies but underpredicts absolute radii by $\sim$5%.
  2. $\Delta$NNLOGO: Includes $\Delta$-isobar degrees of freedom; adjusted to better reproduce absolute radii via saturation properties.
  3. 1.8/2.0 (EM7.5): Features a refitted three-body low-energy constant ($c_D = 7.5$) to match $^{16}$O data, drastically improving radii in doubly magic nuclei.

• Computational Details:

  • Basis: Spherical harmonic oscillator (HO) with 13 major shells ($e_{max}=12$).
  • Uncertainty Quantification: Model-space uncertainties were gauged by varying the HO frequency ($\hbar\omega \in [10, 14]$ MeV) and the shell truncation ($e_{max}$).
  • Operator: The mean-square charge radius includes point-proton radius, proton/neutron finite size corrections, spin-orbit corrections, and the Darwin-Foldy term.

3. Key Contributions

• Systematic Benchmarking: This is one of the first comprehensive ab initio studies of charge radii across the entire Sn chain ($A=96-150$) using BCC theory, extending beyond the traditional doubly magic limits.
• Decomposition of Isotopic Shifts: The authors decompose the isotopic shifts into linear and quasi-parabolic components, allowing for a detailed comparison of how different Hamiltonians capture the underlying nuclear structure evolution.
• Identification of "Wrong Reason" Success: The study reveals that while the 1.8/2.0 (EM7.5) Hamiltonian successfully reproduces the experimental kink at $^{132}$Sn, it does so due to an incorrect prediction of the single-particle shell ordering beyond $N=82$.
• Prediction of Future Challenges: The paper highlights that the region beyond $^{132}$Sn (specifically $^{142}$Sn) acts as a "stress test," where different Hamiltonians diverge significantly, exposing critical deficiencies in current $\chi$EFT interactions.

4. Key Results

A. Absolute Charge Radii

• 1.8/2.0 (EM): Significantly underestimates experimental radii (RMS error $\sim$5.6%).
• $\Delta$NNLOGO & 1.8/2.0 (EM7.5): Both significantly improve the agreement, reducing RMS errors to 1.2% and 0.7%, respectively.
• Correlation Effects: Many-body correlations (BCCSD vs. HFB) generally improve agreement with experiment. They increase radii when the mean-field baseline is too low and decrease them when the baseline is too high.

B. Isotopic Shifts and the $^{132}$Sn Kink

• General Trend: All three Hamiltonians reproduce the general parabolic trend between $N=50$ and $N=82$ reasonably well.
• The $^{132}$Sn Kink:
  • $\Delta$NNLOGO: Underestimates the rise from $^{132}$Sn to $^{134}$Sn (predicts $\sim$0.11 fm$^2$ vs. experimental 0.22 fm$^2$).
  • 1.8/2.0 (EM7.5): Reproduces the experimental rise accurately (0.23 fm$^2$).
  • Mechanism: The success of EM7.5 is attributed to neutrons occupying the spatially extended $1h_{9/2}$ orbital immediately after $N=82$. This pulls protons outward, increasing the radius. However, this shell ordering contradicts experimental spin data for odd-even neighbors (e.g., $^{133}$Sn), suggesting the success is accidental.

C. The $^{142}$Sn Anomaly

• Inverted Kink: The 1.8/2.0 (EM7.5) Hamiltonian predicts a sharp "inverted kink" (decrease) in the isotopic shift at $^{142}$Sn, correlated with a large, unphysical rise in two-neutron separation energies ($S_{2n}$).
• Shell Structure: This behavior arises because EM7.5 places the $2f_{7/2}$ orbital (less spatially extended) as the valence shell beyond $N=92$.
• Contrast: The 1.8/2.0 (EM) and $\Delta$NNLOGO Hamiltonians do not show this inverted kink, as they predict a different shell ordering where the $2f_{7/2}$ shell is filled earlier.

D. Linear and Parabolic Components

• The linear slope of the isotopic shifts is well captured by $\Delta$NNLOGO and EM7.5.
• The quasi-parabolic component is well reproduced by 1.8/2.0 (EM) and $\Delta$NNLOGO but significantly underestimated by EM7.5.
• Local irregularities (sub-shell closures) are visible in BCCSD results but smoothed out in experimental data, suggesting missing collective correlations in the theory.

5. Significance and Outlook

• Validation of $\chi$EFT: The study demonstrates that while current $\chi$EFT interactions can be tuned to fit specific observables (like the $^{132}$Sn kink), they often fail to capture the underlying single-particle structure correctly. This highlights the need for interactions that are robust across different mass regions and observables.
• Future Directions:
  • Experimental: New measurements of isotopic shifts are urgently needed for nuclei near $^{100}$Sn and beyond $^{134}$Sn (up to $^{142}$Sn) to constrain the Hamiltonians.
  • Theoretical:
    • Inclusion of triples corrections ($T_3$) to BCCSD to improve binding energy accuracy and potentially refine radii.
    • Inclusion of two-body charge density corrections to the operator.
    • Development of methods to capture truly collective correlations that are difficult to access via single-reference expansions.
• Conclusion: The Sn isotopic chain, particularly the region beyond $N=82$, serves as a unique "playground" to pin down critical attributes of $\chi$EFT inter-nucleon interactions. The inability of current models to consistently predict both binding energies and charge radii in this region remains a major challenge for ab initio nuclear structure theory.