Neutron skins probed in proton knockout from neutron-rich nuclei

This paper establishes a unified theoretical framework demonstrating that proton-induced knockout reactions, particularly two-proton removal, serve as sensitive probes of neutron-skin thickness and the nuclear symmetry energy by exhibiting systematic decreases in cross sections and momentum widths as neutron excess increases.

The Gist

Imagine an atomic nucleus not as a solid marble, but as a bustling city. In the center, you have the "downtown" (the core), packed tightly with both protons (positive charge) and neutrons (neutral). But on the outskirts, the "suburbs," things get interesting. In most stable atoms, the suburbs are a mix of both. However, in neutron-rich nuclei (atoms with too many neutrons), the suburbs become a "neutron-only" zone. This extra layer of neutrons sticking out is called a neutron skin.

Measuring how thick this skin is is a huge deal for physicists. It tells us about the "rules of the road" inside the nucleus and even helps us understand the structure of neutron stars (the super-dense remnants of exploded stars).

This paper proposes a new, clever way to measure the thickness of this neutron skin using a specific type of "collision experiment."

The Experiment: A High-Speed Pinball Game

Imagine firing a high-speed proton (a single particle) at a heavy atom.

• The Old Way (Electron Scattering): Scientists used to shoot electrons at atoms. Electrons are polite; they barely bump into anything and just glance off, giving a clean picture of the inside. But you can't do this easily with unstable, radioactive atoms found in space or rare labs.
• The New Way (Proton Knockout): This paper focuses on shooting protons at these rare atoms. It's like playing a game of high-speed pinball. You shoot a proton in, and it smashes into the atom, knocking other protons out.
  • The (p, 2p) Reaction: The incoming proton hits one proton inside, and both fly out.
  • The (p, 3p) Reaction: The incoming proton hits one proton, which then hits another proton, and all three fly out. This is a "domino effect" or a "chain reaction."

The Problem: The "Fog" of Neutrons

Here is the tricky part: When you shoot a proton into a neutron-rich atom, it has to travel through that "neutron-only" suburb (the skin) to get to the downtown core.

Think of the neutrons as a thick, sticky fog.

• Protons hate traveling through this fog. They get absorbed or slowed down easily because they interact strongly with neutrons.
• If the neutron skin is thick, the fog is dense. The incoming proton gets stopped before it can reach the core protons. The outgoing protons also get stuck trying to escape.
• Result: You see fewer "knocked out" protons than you expected. The "signal" is weaker.

The authors realized that by counting exactly how many protons get knocked out, they can measure how thick the fog (the neutron skin) is.

The Big Discovery: The "Double Knockout" is a Better Detective

The paper compares two scenarios:

1. Knocking out 1 proton: This gives you a general idea of the fog.
2. Knocking out 2 protons (the chain reaction): This is the paper's star player.

The Analogy:
Imagine trying to see how thick a forest is.

• If you throw a ball and it hits one tree and bounces out, you learn a little about the forest edge.
• If you throw a ball that hits a tree, which then hits a second tree, and then they both bounce out, that ball had to travel deeper into the forest to find that second tree.

The authors found that the two-proton knockout (p, 3p) is much more sensitive to the neutron skin than the one-proton version. Because it involves a chain reaction, it is more likely to be "filtered out" by the thick neutron fog. If the skin is thick, the double-knockout signal drops dramatically. This makes it a super-sensitive detector for the thickness of the neutron layer.

The "Momentum" Clue

The paper also looks at how fast the remaining piece of the atom (the "residue") flies away after the collision.

• The Goldhaber Model (The Old Map): This is like a statistical guess. It assumes the atom is a uniform bag of gas. It predicts a specific speed for the flying residue, but it doesn't care about the shape of the suburbs.
• The New Model (The GPS): The authors built a computer simulation that accounts for the fact that protons are more likely to be knocked out from the surface (the suburbs) than the deep core, especially when there's a thick neutron skin.

They found that the speed (momentum) of the flying residue changes depending on how thick the neutron skin is. The thicker the skin, the slower the residue flies. This is a new "fingerprint" that the old statistical models missed.

Why Does This Matter?

1. Mapping the Unknown: We can't easily measure the neutron skin of unstable, radioactive atoms (like those found in supernovas) with current tools. This method uses proton beams, which are easier to generate in modern labs.
2. Unlocking Neutron Stars: The "rules" that determine how thick a neutron skin is on Earth are the same rules that determine how big a neutron star is. If we know the skin thickness, we can predict the size of these cosmic giants.
3. Better Models: The paper shows that simple statistical guesses aren't enough. We need to understand the detailed "geometry" of the collision to get accurate answers.

Summary

In short, this paper says: "If you want to measure the invisible neutron skin on a rare atom, don't just knock out one proton. Knock out two in a chain reaction. It's like using a deeper probe; the signal will change much more dramatically if the skin is thick, giving us a clearer picture of the nuclear world and the secrets of neutron stars."

Technical Summary

1. Problem Statement

The paper addresses the challenge of probing the neutron-skin thickness ($\Delta R_{np}$) of neutron-rich nuclei and the density dependence of the nuclear symmetry energy. While electron scattering and parity-violating electron scattering provide clean probes of charge and weak charge distributions, hadronic probes like proton-induced knockout reactions offer complementary information but are complicated by strong final-state interactions (FSI) and absorption.

Specifically, the authors aim to:

• Develop a unified theoretical framework to calculate inclusive $(p,2p)$ (single-proton knockout) and sequential $(p,3p)$ (two-proton knockout) reaction cross sections and fragment momentum distributions.
• Quantify the sensitivity of these reactions to neutron-skin thickness along isotopic chains (specifically Tin isotopes).
• Determine whether $(p,3p)$ reactions offer enhanced sensitivity to isovector nuclear structure compared to $(p,2p)$ reactions.
• Distinguish between reaction-induced effects (attenuation, surface bias) and genuine many-body correlations (short-range correlations).

2. Methodology

The authors propose a probabilistic extension of Glauber multiple-scattering theory combined with microscopic nuclear densities.

• Theoretical Framework:

  • Eikonal Approximation: The reaction is modeled using eikonal waves for nucleon propagation between binary collisions. Unlike standard Glauber theory which assumes straight-line trajectories for the entire path, this model allows protons to scatter laterally after each collision.
  • Probabilistic Propagator: Instead of calculating complex quantum interference amplitudes, the authors use a reaction propagator $R(b, z_1, z_2)$ representing the survival probability of a nucleon. This simplifies the treatment of sequential collisions (e.g., a projectile proton hits a target proton, then the leading projectile proton hits a second target proton).
  • Reaction Dynamics: The $(p,3p)$ process is treated as a sequential two-step collision:
    1. Incoming proton collides with a bound proton (scattering angle $\Omega_1$).
    2. The leading proton collides with a second bound proton (scattering angle $\Omega_2$).
  • Cross Section Calculation: The total cross section is obtained by integrating over impact parameters ($b$), collision coordinates ($r_1, r_2$), and scattering angles, summing over all occupied single-particle orbitals.

• Nuclear Structure Inputs:

  • Densities: Microscopic ground-state densities are obtained from Hartree-Fock-Bogoliubov (HFB) calculations using 22 different Skyrme energy-density functionals (e.g., SIII, SLy4, SKM*, UNEDF). This allows the authors to assess model dependence and the correlation with symmetry energy parameters.
  • Single-Particle Wavefunctions: Calculated using a Woods-Saxon potential with spin-orbit coupling, adjusted to reproduce experimental separation energies.
  • Interaction: The elementary proton-nucleon cross sections are derived from the SAID model (partial-wave analysis). Medium modifications to the cross sections are included.

• Momentum Distribution Formalism:

  • A vectorized Monte Carlo (MC) abrasion model is used to calculate the longitudinal momentum dispersion ($\Delta_\parallel$) of the residual fragment.
  • The local Fermi momentum is sampled from the local density ($p_F(r) \propto \rho(r)^{1/3}$).
  • The recoil momentum of the fragment is the negative vector sum of the intrinsic momenta of the removed nucleons.

3. Key Contributions

1. Unified Formalism for $(p,2p)$ and $(p,3p)$: The paper provides a consistent probabilistic framework that handles both single and double proton removal, explicitly accounting for sequential collision dynamics and lateral scattering, which is crucial for $(p,3p)$ reactions.
2. Quantification of Neutron-Skin Sensitivity: The study systematically demonstrates that reaction cross sections and momentum widths are highly sensitive to the neutron-skin thickness, primarily due to the difference in proton-neutron ($\sigma_{pn}$) vs. proton-proton ($\sigma_{pp}$) cross sections.
3. Superiority of $(p,3p)$ as a Probe: The authors establish that $(p,3p)$ reactions are more sensitive to neutron-skin variations than $(p,2p)$ reactions because the double-removal process probes deeper into the nuclear interior and is more strongly affected by the surface-to-volume ratio changes induced by the skin.
4. Critique of Statistical Models: The work highlights the limitations of the Goldhaber model (a statistical model assuming a uniform Fermi gas), showing it fails to reproduce the mass dependence and neutron-skin sensitivity observed in the microscopic MC calculations.

4. Key Results

• Cross Section Trends:

  • Both $(p,2p)$ and $(p,3p)$ total cross sections exhibit a systematic decrease as neutron excess (and thus neutron-skin thickness) increases.
  • This reduction is caused by enhanced attenuation of the incoming and outgoing protons by the neutron-rich surface layer (since $\sigma_{pn} > \sigma_{pp}$ at intermediate energies).
  • The decrease is more pronounced for $(p,3p)$, with a sensitivity derivative estimated as $\frac{d \ln \sigma_{p,3p}}{d \Delta R_{np}} \sim -(0.2 - 0.4) \text{ fm}^{-1}$.
  • Calculations using 22 Skyrme functionals show a spread in predicted cross sections that correlates with the spread in predicted neutron-skin thicknesses.

• Momentum Distributions:

  • The longitudinal momentum dispersion ($\Delta_\parallel$) of the residual fragment decreases with increasing neutron skin thickness.
  • Physical Mechanism: A thicker neutron skin pushes the proton removal probability toward the nuclear surface, where the local density (and thus the local Fermi momentum) is lower. This results in a narrower momentum distribution.
  • Hierarchy: $\Delta_\parallel(p,3p) > \Delta_\parallel(p,2p)$, consistent with the recoil mechanism involving the sum of two independent momenta.
  • Comparison with Goldhaber: The Goldhaber model predicts a constant width (dependent only on mass number) and fails to capture the reduction in width caused by the neutron skin, whereas the microscopic MC model reproduces the trend.

• Model Dependence: The results confirm that the sensitivity to the symmetry energy slope parameter ($L$) is robust. Larger $L$ values predict thicker skins, which in turn predict smaller cross sections and narrower momentum widths.

5. Significance

• Complementary Probe: The paper establishes proton-induced knockout reactions as a powerful, complementary hadronic probe for neutron skins, offering constraints on the nuclear symmetry energy that are distinct from electroweak probes (like PREX/CREX).
• Isovector Structure: The enhanced sensitivity of $(p,3p)$ reactions suggests they are ideal for future experiments at facilities like FRIB and FAIR to constrain the isovector sector of the nuclear energy density functional.
• Astrophysical Implications: Since the neutron skin thickness is directly correlated with the radius and crust structure of neutron stars, improved constraints from $(p,3p)$ reactions can refine our understanding of the neutron star equation of state and gravitational wave signatures.
• Methodological Advance: By moving beyond the "straight-line" approximation and statistical models, the authors provide a more realistic description of reaction dynamics that separates geometric/attenuation effects from intrinsic nuclear correlations, aiding in the interpretation of spectroscopic factors.

In conclusion, the authors demonstrate that the systematic reduction of $(p,2p)$ and $(p,3p)$ cross sections and momentum widths in neutron-rich nuclei is a direct signature of the neutron skin, with $(p,3p)$ reactions offering a particularly sensitive window into the density dependence of the nuclear symmetry energy.