⚛️ nuclear theory

Neutron skins: A perspective from dispersive optical models

This paper reviews neutron skin predictions from a dispersive optical model that successfully describes both ${}^{208}$ Pb and ${}^{48}$ Ca, highlighting the model's unique ability to treat structure and reaction data consistently while underscoring the ongoing puzzle of reconciling a large lead skin with the unexpectedly thin calcium skin observed in CREX.

Original authors: M. C. Atkinson, W. H. Dickhoff

Published 2026-07-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: M. C. Atkinson, W. H. Dickhoff

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

The Big Mystery: The "Neutron Skin"

Imagine an atom's nucleus as a crowded dance floor. The dancers are protons (positively charged) and neutrons (neutral). Usually, they dance together in a tight mix. But in some heavy atoms, like Calcium-48 and Lead-208, there are way more neutrons than protons.

The big question physicists have is: Where do these extra neutrons go?

Do they mix evenly throughout the dance floor?
Or do they crowd to the very edge, forming a "skin" around the core?

The thickness of this "neutron skin" is a crucial clue. It tells us how nuclear matter behaves under extreme pressure, which helps us understand everything from the heaviest atoms to neutron stars (which are basically giant atomic nuclei floating in space).

The Conflict: Two Experiments, Two Different Answers

Recently, two major experiments tried to measure this skin thickness, and they got very different results:

PREX-2 (Lead-208): Found a thick skin.
CREX (Calcium-48): Found a thin skin.

This created a puzzle. Most computer models (theories) say that if Lead has a thick skin, Calcium should also have a thick skin. They are "cousins" in the nuclear family, so they should behave similarly. The fact that one is thick and the other is thin has left scientists scratching their heads.

The New Detective: The Dispersive Optical Model (DOM)

The authors of this paper used a special tool called the Dispersive Optical Model (DOM).

The Analogy: The Master Architect
Imagine you are trying to understand a house.

Old Methods (Mean-Field Theory): These are like looking at a blueprint. They guess how the house is built based on general rules of architecture. They often get the general shape right but miss the specific details of how the rooms connect.
The DOM Method: This is like the Master Architect who has a complete set of blueprints and a record of every time someone walked through the front door, bounced off a wall, or sat in a chair.

The DOM is unique because it treats scattering data (how particles bounce off the nucleus, like throwing a ball at a wall) and bound-state data (how particles sit inside the nucleus, like furniture in a room) as part of the same story. It uses a mathematical "bridge" (called a dispersion relation) to connect what happens outside the nucleus to what happens inside. If you know how the ball bounces, you can deduce exactly how the furniture is arranged inside.

What the DOM Found

The authors used this "Master Architect" approach to analyze Calcium-48 and Lead-208 using all available data.

For Lead-208: The DOM predicted a thick skin (about 0.25 fm). This matched the PREX-2 experiment perfectly.
For Calcium-48: The DOM also predicted a thick skin (about 0.25 fm).

The Problem: The DOM's prediction for Calcium-48 clashes with the CREX experiment, which found a very thin skin.

Why the Discrepancy? (The Missing Puzzle Piece)

The paper argues that the DOM is very good at predicting things when it has enough data.

For Lead-208, they had a mountain of data (how protons and neutrons bounce off it at many different speeds). The model was confident.
For Calcium-48, the data is much scarcer. They have plenty of data on how protons bounce off it, but very little data on how neutrons bounce off it.

The Analogy: The Blindfolded Detective
Imagine trying to solve a crime in a dark room.

In the Lead case, the lights are on; the detective sees everything clearly.
In the Calcium case, the detective is wearing a blindfold for half the room (the neutron side). The DOM tried to guess the layout based on the protons and the little bit of neutron data it had, and it guessed "Thick Skin."

The authors suggest that the "thin skin" result from CREX might be correct, but the DOM couldn't see it yet because it lacked enough neutron "bouncing" data to force the model to change its mind.

The "What If" Scenario

The paper explores what would happen if we forced the model to accept the "thin skin" result from CREX.

The Consequence: If the skin is thin, the extra neutrons must be hiding deeper inside the nucleus, closer to the center.
The Ripple Effect: According to the laws of physics (specifically the Heisenberg Uncertainty Principle), if you squeeze neutrons into a smaller space, they have to move faster. This would mean Calcium-48 should have a lot more "fast-moving" (high-momentum) neutrons than protons.
The Conflict: This contradicts what the DOM currently sees (and what other experiments suggest), which is that there are actually more fast-moving protons than neutrons.

The Conclusion

The paper concludes that the "Neutron Skin Puzzle" isn't necessarily a failure of the theory, but a sign that we need more data.

The DOM is a powerful tool that links different types of experiments together.
It successfully predicted the thick skin in Lead.
It predicts a thick skin in Calcium, but this prediction is likely shaky because we don't have enough neutron scattering data for Calcium yet.
To solve the mystery, scientists need to measure how neutrons bounce off Calcium-48 much more precisely. Once they do that, the DOM should be able to tell us the true thickness of the skin and help us understand the rules that govern neutron stars.

In short: The model says "Thick Skin," the experiment says "Thin Skin," and the authors say, "We need to look closer at the Calcium data before we can be sure who is right."

Technical Summary: Neutron Skins: A Perspective from Dispersive Optical Models

Problem Statement
The distribution of neutrons relative to protons in neutron-rich nuclei is a fundamental question in nuclear physics, quantified by the neutron-skin thickness ( $R_{skin} = R_n - R_p$ ). This observable constrains the density dependence of the nuclear symmetry energy, which has profound implications for nuclear masses, drip lines, and astrophysical phenomena such as neutron star structure. A significant puzzle has emerged following recent experimental results: the PREX-2 experiment measured a thick neutron skin in $^{208}\text{Pb}$ ( $R_{skin} \approx 0.283$ fm), while the CREX experiment measured a significantly thinner skin in $^{48}\text{Ca}$ ( $R_{skin} \approx 0.121$ fm). Standard theoretical approaches, including mean-field and ab initio methods, generally predict a strong positive correlation between the skins of these two nuclei, failing to simultaneously reproduce a thick skin in $^{208}\text{Pb}$ and a thin skin in $^{48}\text{Ca}$ . Furthermore, most experimental determinations of neutron radii are model-dependent, with the exception of parity-violating electron scattering.

Methodology
The authors employ an empirical nonlocal Dispersive Optical Model (DOM) to address this discrepancy. Unlike mean-field or ab initio approaches that often treat structure and reactions separately, the DOM unifies these domains by linking bound-state and scattering data through a subtracted dispersion relation.

Theoretical Framework: The method utilizes the Dyson equation to determine the single-particle propagator (Green's function), $G_{\ell j}(r, r'; E)$ , from the nucleon self-energy, $\Sigma_{\ell j}$ . The self-energy acts as a complex, nonlocal, and energy-dependent potential.
Dispersion Relation: The real part of the self-energy is constrained by the imaginary part via a subtracted dispersion relation, using the Fermi energy ( $\varepsilon_F$ ) as the subtraction point. This ensures consistency between the potential at negative energies (bound states) and positive energies (scattering).
Nonlocality: The model explicitly includes nonlocality in both the static Hartree-Fock (HF) term and the imaginary part of the self-energy. This is critical for accurately calculating particle numbers, charge densities, and spectral functions, which local, energy-dependent potentials fail to reproduce correctly.
Data Constraints: The model parameters are constrained via weighted $\chi^2$ $χ^{2}$ minimization against a comprehensive dataset including:
- Elastic scattering angular distributions (protons and neutrons).
- Analyzing powers.
- Total and reaction cross-sections.
- Charge densities (from electron scattering).
- Single-particle energy levels.
- Ground-state binding energies and particle numbers.
Specifics for $^{48}\text{Ca}$ and $^{208}\text{Pb}$ : For $^{48}\text{Ca}$ , the isoscalar part is constrained by the nearby $N=Z$ $^{40}\text{Ca}$ nucleus, with only asymmetric parameters fitted. For $^{208}\text{Pb}$ , both isoscalar and isovector parameters were varied. The authors note that while charge-exchange data could further constrain the isovector potential, it was not included in the fits reviewed here.

Key Results

Predicted Neutron Skins: The DOM analysis predicts a thick neutron skin for both nuclei:
- $^{48}\text{Ca}$ : $R_{skin} = 0.25 \pm 0.023$ fm.
- $^{208}\text{Pb}$ : $R_{skin} = 0.25 \pm 0.05$ fm.
  The prediction for $^{208}\text{Pb}$ is consistent within $1\sigma$ of the PREX-2 measurement. However, the prediction for $^{48}\text{Ca}$ is over $2\sigma$ larger than the thin skin reported by CREX.
Data Consistency: The model successfully reproduces the full range of experimental data used for fitting, including elastic scattering cross-sections, analyzing powers, charge densities, and single-particle energy levels. The inclusion of the dynamic part of the self-energy is shown to be necessary to reduce the particle-hole gap predicted by static mean-field models.
Sensitivity Analysis: The authors investigated the sensitivity of the $^{48}\text{Ca}$ skin to the input data. They found that the total neutron cross-section data (scattering domain) provides a strong constraint on the skin thickness. When the skin was artificially forced to the CREX value ($0.132$ fm), the fit to the total cross-section data degraded significantly, suggesting that the current neutron scattering data favors a thicker skin.
Momentum Distributions and Spectroscopic Factors: The model predicts significant high-momentum content in both nuclei due to short-range correlations (SRC), with a higher fraction of high-momentum protons than neutrons. This "np-dominance" is consistent with CLAS experiments. The spectroscopic factors derived from the model are consistent with $(e, e'p)$ reaction data.

Significance and Limitations
The paper argues that the DOM offers a unique perspective on the "neutron skin puzzle" because it treats structure and reaction data on the same footing through the dispersion relation. The authors claim that their results for $^{208}\text{Pb}$ are robust given the available data. However, regarding $^{48}\text{Ca}$ , they explicitly state that the discrepancy with CREX likely arises from an insufficient dataset for neutrons in the scattering domain.

Data Deficit: Unlike protons in $^{48}\text{Ca}$ , which are well-constrained by multiple scattering datasets, neutron data is sparse (limited to a few elastic scattering energies and total cross-sections). The authors suggest that the thin skin reported by CREX implies a redistribution of neutrons toward the nuclear interior, which would increase high-momentum neutrons and potentially invert the current proton-neutron high-momentum hierarchy.
Future Requirements: The authors conclude that to resolve the discrepancy, more neutron data must be measured in both the scattering and bound-state domains for $^{48}\text{Ca}$ . Specifically, including the CREX result as a constraint in a future DOM fit is necessary to bring the neutron dataset closer to "completeness."
Theoretical Implications: The paper posits that the size of $^{48}\text{Ca}$ may be too small to apply bulk nuclear matter properties (like the Equation of State) directly, and that shell effects (specifically the $f_{7/2}$ closure) might play a more dominant role in skin formation than macroscopic asymmetry alone.

The authors maintain that parity-violating electron scattering remains the most unambiguous approach to determine neutron skins, and that the DOM provides a rigorous framework to link these observables to the underlying nuclear self-energy and Equation of State.