Matter radii from interaction cross sections using microscopic nuclear densities

By integrating uncertainty-quantified nuclear densities from the Fayans energy density functional with a modernized Glauber reaction framework, this study reanalyzes calcium isotopic interaction cross sections to refute earlier claims of dramatic neutron swelling while establishing a robust methodology for extracting neutron radii across the nuclear chart.

The Gist

Imagine the atomic nucleus not as a tiny, solid marble, but as a bustling city made of two types of citizens: Protons (who carry a positive charge) and Neutrons (who are neutral).

For a long time, scientists have been very good at measuring the size of the "Proton City" because protons interact with electricity, making them easy to spot. However, measuring the size of the "Neutron City" is much harder. Neutrons are shy; they don't play well with electricity. To see them, scientists have to crash other particles into the nucleus and measure how much "space" the nucleus takes up during the collision. This is called an interaction cross-section.

The Problem: The "Swelling" Mystery

Recently, some scientists looked at Calcium isotopes (atoms with the same number of protons but different numbers of neutrons) and claimed something dramatic was happening. As they added more neutrons, the nucleus seemed to suddenly swell up like a balloon, creating a massive "neutron skin" on the outside.

This was a big deal because if true, it would change our understanding of how matter behaves inside neutron stars (which are basically giant neutron balls). But there was a catch: the way they calculated this size relied on a lot of guesswork and simplified models. It was like trying to guess the size of a house by looking at its shadow, but using a broken ruler and a cloudy day.

The New Approach: A High-Definition Camera

In this paper, the authors (Smith, Godbey, Hebborn, and colleagues) decided to build a better "camera" to take the measurement. They created a new, integrated pipeline that connects two things that are usually kept separate:

1. The Blueprint (Structure): They used a sophisticated computer model called Fayans EDF to predict exactly how protons and neutrons are arranged inside the nucleus. Think of this as a high-resolution 3D blueprint of the city.
2. The Crash Test (Reaction): They used a modern version of the Glauber model (a set of physics rules for particle collisions) to simulate what happens when these nuclei crash into a Carbon target.

The Secret Sauce:
Previous studies used a "one-size-fits-all" rulebook for how particles bounce off each other. This new team realized that rulebook was a bit sloppy. So, they recalibrated their rulebook using data from stable Calcium isotopes (the ones we know well) before applying it to the weird, unstable ones.

It's like a tailor measuring a person for a suit. Instead of using a generic size chart, they first measure a known person perfectly, adjust their measuring tape to be 100% accurate for that person, and then use that perfectly calibrated tape to measure the mystery person.

The Results: No Giant Swelling

When they ran their new, high-precision simulation on the Calcium isotopes, the results were surprising:

• The "Swelling" was an illusion: The dramatic, sudden explosion in size that earlier studies reported? It didn't happen.
• Growth is gradual: The neutron skin does grow as you add more neutrons, but it grows slowly and steadily, like a tree growing rings, not like a balloon popping.
• The Theory Matches Reality: Their new method showed that the theoretical models (the blueprints) actually agree with the experimental data once you account for all the errors and uncertainties.

Why Does This Matter?

Think of this like fixing a map.

• Old Map: Showed a massive, unexplained mountain range where there was actually just a gentle hill.
• New Map: Shows the gentle hill accurately.

This matters because:

1. Neutron Stars: If the "neutron skin" isn't as thick as we thought, it changes our calculations of how heavy and dense neutron stars can be before they collapse.
2. Future Experiments: Now that we have a reliable "calibrated tape measure," scientists at facilities like FRIB (where they make rare isotopes) can trust their new data. They can stop worrying about whether their tools are broken and start focusing on discovering new physics.

The Bottom Line

The authors didn't just find a new number; they built a better toolkit. They showed that by combining a precise theoretical blueprint with a carefully calibrated collision model, we can stop guessing and start knowing exactly how big these tiny nuclear cities really are. The "dramatic swelling" was likely just a measurement error, and the truth is a bit more modest, but much more reliable.

Technical Summary

1. Problem Statement

Understanding the evolution of nuclear size with proton and neutron numbers is critical for testing models of strongly interacting matter, particularly regarding isovector properties and neutron star structure. While nuclear charge radii are well-measured via electromagnetic probes, neutron radii and matter radii are difficult to measure directly.

• Current Limitation: The primary method for extracting neutron radii of short-lived, unstable nuclei is analyzing interaction cross sections ($\sigma_I$) using the Glauber model, specifically the Modified Optical Limit Approximation (MOL).
• The Discrepancy: Previous analyses (e.g., Tanaka et al., Phys. Rev. Lett. 124, 102501) suggested a "dramatic neutron swelling" (a rapid increase in neutron skin) in neutron-rich Calcium isotopes ($^{48-51}$Ca). However, these extractions relied on schematic two-parameter Fermi density models fitted to data, lacking a rigorous, end-to-end uncertainty quantification. This led to tensions between theoretical predictions (from Density Functional Theory) and experimental extractions, with unclear sources of error (theoretical model vs. reaction analysis).

2. Methodology

The authors developed an integrated, uncertainty-quantified pipeline linking microscopic nuclear structure calculations directly to reaction observables.

• Nuclear Structure (Input):

  • Used the Fayans Energy Density Functional (EDF), specifically the parametrization Fy(IVP3,0.9), which includes isospin-variable pairing (IVP) terms.
  • Instead of a single "best-fit" model, they generated an ensemble of 500 physically credible EDF samples (within 1$\sigma$, 2$\sigma$, and 3$\sigma$ confidence bands) based on the covariance matrix from the EDF calibration.
  • These samples provided self-consistent proton and neutron densities for $^{12}$C (target) and Calcium isotopes $^{42-51}$Ca (projectiles).

• Reaction Framework:

  • Employed the Modified Optical Limit Approximation (MOL) within the Glauber framework to calculate interaction cross sections.
  • Key Innovation (Recalibration): Unlike previous works that used fixed nucleon-nucleon (NN) profile functions from literature, the authors recalibrated the NN profile function parameters ($\alpha_{pp}, \alpha_{pn}, \alpha_{nn}$) for each EDF density sample.
  • The calibration was performed by minimizing a $\chi^2$ function against experimental cross-section data for stable/long-lived isotopes ($^{42-48}$Ca). This step accounts for in-medium effects and systematic deficiencies in the MOL approximation.

• Uncertainty Quantification (UQ):

  • Propagated three distinct sources of uncertainty to the final cross-section predictions:
    1. EDF Parametric Uncertainty: Variations in nuclear densities from the EDF ensemble.
    2. NN Cross-Section Uncertainty: Uncertainties in the input total NN cross sections used in the profile function.
    3. Profile Function Calibration Uncertainty: Uncertainties arising from the fitting of the profile function parameters (using covariance matrices).
  • The total uncertainty was derived from a massive ensemble of 250,000 calculations per isotope (500 EDF samples $\times$ 500 profile function samples).

3. Key Results

• Resolution of the "Neutron Swelling" Puzzle:

  • When applying the rigorous pipeline to the Calcium chain, the authors found no evidence for the dramatic neutron swelling reported in earlier analyses.
  • The theoretical predictions for interaction cross sections, when accounting for all uncertainties, are consistent with experimental data at the 2$\sigma$ and 3$\sigma$ confidence levels.
  • The apparent discrepancy in previous studies arose because those analyses ignored the large uncertainties associated with the reaction model inputs and the correlations between isotopes.

• Performance of the Fayans EDF:

  • The Fayans EDF successfully reproduces the odd-even staggering of cross sections and the complex pattern of charge radii along the Ca chain.
  • The predicted neutron skins (difference between neutron and proton radii) show a more modest growth across the chain compared to the steep rise suggested by previous data extractions.

• Uncertainty Analysis:

  • EDF Uncertainties: Negligible compared to experimental errors.
  • NN Cross-Section Uncertainties: Negligible due to the recalibration of profile functions, which compensates for inaccuracies in evaluated NN data.
  • Profile Function Uncertainties: The dominant theoretical uncertainty (approx. 3% of the cross section), but still significantly smaller than uncertainties in other reaction probes.
  • Correlations: The profile function parameters are strongly correlated across the isotopic chain; a single set of parameters shifts all cross sections up or down together. This correlation is crucial for the $\chi^2$ test.

• Model Robustness:

  • Sensitivity tests confirmed that using isospin-independent pairing (Fy($\Delta r$,HFB)) increases the discrepancy, validating the need for isospin-dependent pairing.
  • Results remained robust when changing the calibration dataset or using the standard Optical Limit Approximation (OLA) instead of MOL.

4. Significance and Impact

• Methodological Advancement: The paper establishes a new standard for analyzing interaction cross sections by creating a self-consistent, end-to-end framework that treats structure and reaction on equal footing with rigorous statistical error propagation.
• Correction of Previous Interpretations: It demonstrates that the "dramatic neutron swelling" in $^{48-51}$Ca was likely an artifact of insufficient uncertainty quantification in previous reaction analyses, rather than a failure of nuclear density functional theory.
• Future Applications: The modular pipeline is immediately applicable to rare isotope facilities (e.g., FRIB, RIKEN). It allows for the rapid, statistically rigorous extraction of matter and neutron radii from new high-statistics data, enabling direct constraints on EDF parametrizations and improving our understanding of isovector nuclear properties and neutron star equations of state.
• Call to Action: While the 2$\sigma$ tension for the most neutron-rich isotopes ($^{51}$Ca) remains, the framework provides a clear path forward: more precise measurements of $^{51}$Ca and heavier Calcium isotopes are needed to definitively constrain the models.