Constraining the neutron skin of $^{208}$Pb with anisotropic flow in Pb+Pb collisions at the LHC

This study utilizes an improved multi-phase transport model to analyze Pb+Pb collisions at the LHC, finding that while anisotropic flow measurements can exclude large neutron skins in $^{208}$Pb, they exhibit a geometric degeneracy that limits their ability to precisely distinguish between zero and moderate neutron skin values.

The Gist

Imagine two giant, fluffy marshmallows smashing into each other at nearly the speed of light. These aren't just any marshmallows; they are Lead-208 nuclei, the heavy atoms used in particle accelerators like the Large Hadron Collider (LHC).

When these "marshmallows" collide, they don't just bounce off; they melt into a super-hot, super-dense soup called Quark-Gluon Plasma (QGP). As this soup expands and cools, it flows like a fluid. Scientists study how this fluid flows to learn about the universe's earliest moments.

But here's the twist: This paper asks a different question. Instead of just looking at the soup, the scientists are trying to figure out the shape of the marshmallows before they smashed.

The Mystery: The "Neutron Skin"

Think of an atom like an orange. The center is the core (protons and neutrons mixed together), and the rind is the surface. In heavy atoms like Lead, the neutrons (the neutral particles) tend to puff out a little further than the protons (the positively charged ones).

This extra layer of neutrons sticking out is called the Neutron Skin.

• The Problem: Scientists have been arguing about how thick this skin is. One experiment (PREX-II) said it's quite thick (like a thick orange rind). Another (CREX) said it's very thin (like a paper-thin peel). They can't agree, and it's a big puzzle in physics.

The Experiment: Smashing to Measure

The authors of this paper decided to use the LHC collisions as a giant measuring tape. They asked: "If we change the thickness of the neutron skin in our computer simulations, does it change how the resulting fluid flows?"

They used a sophisticated computer model (a "multi-phase transport model") to simulate millions of collisions. They created five different versions of the Lead nucleus:

1. Negative Skin: Neutrons are actually inside the protons (a weird, compact version).
2. Zero Skin: Neutrons and protons are the same size.
3. Moderate Skin: A standard, thin layer.
4. Thick Skin: A very puffy layer.
5. Super Thick Skin: An enormous layer.

The Analogy: The "Eccentricity" of the Crash

When two round objects collide, they don't always hit perfectly head-on. They often graze each other, creating an oval-shaped overlap zone.

• The "Skin" Effect: If the neutron skin is thick and puffy, the overlap zone looks different than if the skin is thin and tight.
• The Flow: Just like water rushing through a narrow, oval channel flows differently than water in a round channel, the "soup" created in the collision flows differently depending on the shape of the initial crash.

The scientists looked at something called Anisotropic Flow. Think of this as the "directionality" of the splash.

• If the crash is very oval, the fluid shoots out more in one direction (like a squashed water balloon).
• If the crash is rounder, the fluid shoots out more evenly.

The Findings: What Did They Discover?

1. The "Skin" Survives the Crash
The most exciting finding is that the "memory" of the neutron skin survives the entire chaotic explosion. Even though the collision is violent and messy, the initial shape of the nucleus leaves a fingerprint on the final flow of particles. The computer model showed a clear link: Thicker skin = Different flow pattern.

2. The "Goldilocks" Zone
The scientists compared their computer simulations to real data collected by the ALICE experiment at the LHC.

• The "Too Big" and "Too Small" Skins: The simulations with a very thick skin or a negative skin (where neutrons are inside) produced flow patterns that did not match the real data at all. These scenarios were ruled out.
• The "Just Right" Skins: The real data matched best with a zero skin or a moderate skin (around 0.16 fm).

3. The "Geometric Degeneracy" (The Tricky Part)
Here is the catch. While they could rule out the extreme cases, they couldn't perfectly tell the difference between "Zero Skin" and "Moderate Skin."

• The Analogy: Imagine trying to guess the exact thickness of a tire by looking at a car driving in a straight line. If the tire is slightly thin or slightly thick, the car drives almost exactly the same way. You can tell if the tire is flat or giant, but you can't easily tell the difference between a standard tire and a slightly wider one just by watching the car drive.
• The Result: The flow measurements are sensitive to the overall size of the collision, but they aren't sensitive enough to see the tiny, fine details of the surface profile. The "Zero" and "Moderate" skins create such similar collision shapes that the fluid flows almost identically.

The Conclusion

This paper is a bit of a "Yes, but..." story.

• Yes: Heavy-ion collisions at the LHC can tell us about the neutron skin. They successfully ruled out the extreme theories (the super-thick and the negative skins).
• But: They hit a wall of "geometric degeneracy." The current data isn't sharp enough to distinguish between a thin skin and a moderate skin because the overall shape of the collision is dominated by the size of the nucleus, not the tiny surface details.

In simple terms: The LHC is a powerful microscope, but for this specific question, it's like looking at a blurry photo. We can clearly see that the subject isn't a giant or a dwarf, but we can't quite tell if they are wearing a thin hat or a medium hat. To solve the puzzle, we need either sharper photos (better data) or a different angle (different types of collisions).

Technical Summary

Here is a detailed technical summary of the paper "Constraining the neutron skin of 208Pb with anisotropic flow in Pb+Pb collisions at the LHC."

1. Problem Statement

The thickness of the neutron skin in heavy nuclei, specifically $^{208}\text{Pb}$ (defined as $\Delta r_{np} = \sqrt{\langle r_n^2 \rangle} - \sqrt{\langle r_p^2 \rangle}$), is a critical parameter linking nuclear structure to the density dependence of the nuclear symmetry energy and the equation of state of neutron-rich matter. However, current experimental determinations are inconsistent:

• PREX-II (parity-violating electron scattering) suggests a large neutron skin ($\Delta r_{np} \approx 0.283$ fm).
• CREX and theoretical constraints suggest a much thinner skin ($\Delta r_{np} \approx 0.12\text{--}0.19$ fm).

This discrepancy, known as the PREX-CREX puzzle, necessitates independent probes. The authors investigate whether relativistic heavy-ion collisions (specifically Pb+Pb at the LHC) can constrain the neutron skin thickness by examining how the initial nuclear density distribution influences the final-state anisotropic flow ($v_n$) of the Quark-Gluon Plasma (QGP). A key challenge addressed is whether neutron-skin effects survive the complex dynamical evolution of the collision or are washed out by geometric degeneracies.

2. Methodology

The study employs a microscopic Multi-Phase Transport (AMPT) model, specifically an improved string-melting version (AMPT-SM).

• Model Improvements: The version used incorporates modern parton distribution functions, impact-parameter-dependent nuclear shadowing, an improved quark-coalescence mechanism, and updated heavy-flavor production treatments. This ensures a consistent description of bulk particle production and flow across different collision systems.
• Neutron Skin Configurations: The authors systematically varied the neutron density distribution of $^{208}\text{Pb}$ while keeping the proton distribution fixed. They utilized the Woods-Saxon parameterization to generate five distinct configurations with neutron skin thicknesses ($\Delta r_{np}$) of:
  • $-0.19$ fm (compact)
  • $0$ fm (reference)
  • $0.16$ fm (moderate)
  • $0.37$ fm
  • $0.74$ fm (very diffuse)
• Simulation Process:
  1. Initial nucleon positions are sampled based on the specific Woods-Saxon parameters.
  2. The full collision evolution is simulated at $\sqrt{s_{NN}} = 5.02$ TeV, covering initial conditions, partonic scatterings (ZPC), hadronization, and hadronic rescattering (ART).
  3. Observables are calculated and compared against ALICE experimental data.
• Observables Analyzed:
  • Initial State: Eccentricities ($\varepsilon_2, \varepsilon_3$).
  • Final State: Charged-particle multiplicity ($\langle N_{ch} \rangle$), mean transverse momentum ($\langle p_T \rangle$), and anisotropic flow coefficients ($v_2\{2\}, v_3\{2\}, v_2\{4\}$).
• Statistical Analysis: A $\chi^2$ analysis was performed to quantify the agreement between model predictions and ALICE data, focusing on the mid-central centrality range (20–60%).

3. Key Results

A. Initial Geometry and Evolution

• Survival of Effects: Neutron-skin effects are robustly preserved throughout the entire transport evolution. The ordering of initial geometric eccentricities is maintained in the final-state flow.
• Ellipticity ($\varepsilon_2$ and $v_2$):
  • More compact nuclei (smaller/negative $\Delta r_{np}$) yield larger initial elliptic eccentricities ($\varepsilon_2$) and higher elliptic flow ($v_2$).
  • More diffuse nuclei (larger $\Delta r_{np}$) suppress $\varepsilon_2$ and $v_2$.
  • $v_2\{4\}$ (four-particle cumulant) shows a pronounced, systematic dependence on $\Delta r_{np}$, indicating that neutron-skin modifications to the participant geometry survive the partonic and hadronic stages.
• Triangularity ($\varepsilon_3$ and $v_3$):
  • Conversely, $\varepsilon_3$ and $v_3$ increase monotonically with increasing neutron skin thickness. Thicker skins enhance event-by-event spatial irregularities due to the role of peripheral nucleons.
• Multiplicity and $p_T$:
  • $\langle N_{ch} \rangle$ decreases monotonically with increasing $\Delta r_{np}$ due to the dilution of the nuclear overlap region.
  • $\langle p_T \rangle$ shows a similar but weaker decreasing trend.

B. Constraints on Neutron Skin Thickness

• $\chi^2$ Analysis: The comparison with ALICE data reveals a clear selectivity:
  • Excluded: Configurations with extreme negative ($\Delta r_{np} = -0.19$ fm) or very large positive ($\Delta r_{np} = 0.74$ fm) neutron skins yield very high $\chi^2$ values and are strongly disfavored.
  • Preferred: The reference case ($\Delta r_{np} = 0$) and the moderate case ($\Delta r_{np} = 0.16$ fm) provide the best agreement with data, with $\chi^2$ values of 3.43 and 3.67, respectively.
• Geometric Degeneracy: A critical finding is the geometric degeneracy between $\Delta r_{np} = 0$ and $\Delta r_{np} = 0.16$ fm. These two configurations have nearly identical root-mean-square nuclear radii, resulting in indistinguishable transverse density profiles and overlap geometries within current experimental precision. Consequently, anisotropic flow cannot currently distinguish between a zero and a moderate neutron skin.

4. Significance and Conclusions

• Validation of Heavy-Ion Collisions as a Probe: The study confirms that relativistic heavy-ion collisions are a powerful, complementary tool for probing nuclear structure, capable of ruling out extreme neutron-skin scenarios that contradict other experimental data.
• Limitations of Current Observables: The research highlights a fundamental limitation: in large collision systems like Pb+Pb, anisotropic flow is primarily driven by the overall collision geometry and size. It lacks the extreme sensitivity required to resolve fine details of the nuclear surface profile when different configurations yield similar global radii.
• Resolution of the PREX-CREX Tension: While the study cannot definitively resolve the tension between PREX-II and CREX (as it cannot distinguish between 0 and 0.16 fm), it successfully excludes the large neutron skin scenario ($\approx 0.28$ fm) favored by PREX-II within the AMPT framework.
• Future Directions: To break the geometric degeneracy and achieve precise constraints, the authors suggest the need for:
  1. Higher-statistics experimental data.
  2. Observables more sensitive to the fine structure of the nuclear surface.
  3. Studies in smaller or asymmetric collision systems (e.g., isobaric collisions) where geometric degeneracies are less dominant.

In summary, the paper establishes that while anisotropic flow in Pb+Pb collisions provides a robust mechanism to exclude unphysical or extreme neutron-skin configurations, the current precision is insufficient to uniquely determine the exact thickness of the neutron skin due to geometric degeneracies in large systems.