The neutron skin effect in Pb+Pb collisions at 2.76A TeV at the LHC

This study investigates how the neutron skin thickness in lead nuclei influences the space-time evolution of the quark-gluon plasma in 2.76A TeV Pb+Pb collisions at the LHC, finding that it significantly enhances initial spatial anisotropy and elliptic flow, particularly in peripheral collisions and at lower beam energies.

The Gist

Imagine a heavy nucleus, like the lead atom used in these experiments, as a giant, fuzzy ball of marbles. Inside this ball, there are two types of marbles: protons (which have a positive charge) and neutrons (which are neutral).

Usually, we think of these marbles as being mixed together perfectly, like blue and red jellybeans in a jar. But in reality, especially in heavy atoms like lead, the neutral neutrons tend to push a little further out toward the edge than the charged protons. This creates a thin, fuzzy layer of extra neutrons on the outside, like a skin or a hairy coat on the ball. Scientists call this the "neutron skin."

This paper asks a simple but profound question: Does this "hairy coat" change how the ball behaves when it smashes into another ball at nearly the speed of light?

The Big Smash: The "Fireball"

At the Large Hadron Collider (LHC), scientists smash lead nuclei together at incredible speeds. When they collide, they don't just bounce off; they melt into a super-hot, super-dense soup of particles called Quark-Gluon Plasma (QGP). Think of this as a tiny, fleeting "fireball" that exists for a fraction of a second before cooling down and turning back into ordinary particles.

The shape and behavior of this fireball depend entirely on how the two lead balls hit each other.

• Head-on collision: The balls smash flat. The fireball is round.
• Glancing blow: The balls clip each other. The fireball is shaped like a football (oval).

The Experiment: Two Scenarios

The researchers ran a computer simulation to see what happens in two different scenarios:

1. The "Smooth Ball" Scenario: They assumed the protons and neutrons are mixed perfectly, with no extra skin.
2. The "Hairy Ball" Scenario: They added the "neutron skin," making the edge of the nucleus slightly fuzzier and larger.

What They Found

Here is the breakdown of their findings using everyday analogies:

1. The Shape of the Crash (Spatial Anisotropy)
When two "hairy" balls collide, especially in a glancing blow (peripheral collision), that extra layer of neutrons changes the geometry of the crash. It's like two fuzzy tennis balls hitting each other; the fuzziness makes the contact point slightly different than if they were smooth rubber balls.

• Result: The fireball created in the "hairy" scenario is slightly more oval-shaped (more eccentric) than in the "smooth" scenario. This difference is most noticeable when the collision isn't a direct head-on hit.

2. The Flow of the Soup (Elliptic Flow)
Once the fireball is formed, it expands outward like a balloon popping. Because the fireball is oval-shaped, it doesn't expand equally in all directions; it squirts out faster along the short axis of the oval. This is called "elliptic flow."

• Result: Because the "hairy" collision created a more oval fireball, the resulting flow of particles was stronger. The "hairy" skin made the soup squirt out with more force in specific directions.

3. The Messengers: Photons vs. Hadrons
The paper looked at two types of "messengers" coming out of the fireball:

• Hadrons (like pions): These are like particles that get stuck in the soup and only escape when the soup cools down and freezes. They mostly tell us about the surface of the fireball.
  • Finding: The neutron skin didn't change much about these particles. They looked almost the same in both scenarios.
• Thermal Photons (Light): These are like tiny flashes of light emitted throughout the entire life of the fireball, from the very hot center to the cooling edges. They carry information from the inside out.
  • Finding: The photons were very sensitive to the neutron skin. Because the skin changed the initial shape of the fireball, the pattern of light emitted changed significantly. The "elliptic flow" of the light was much stronger when the neutron skin was included.

The "Why" It Matters

Think of the neutron skin like a secret ingredient in a recipe.

• If you are just tasting the crust of the cake (the hadrons), you might not notice the secret ingredient.
• But if you taste the very center of the cake (the photons), the secret ingredient changes the flavor completely.

The researchers found that if you ignore the neutron skin in your calculations, your predictions for how the light (photons) flows will be wrong, especially for glancing collisions. By adding the "neutron skin" into the model, their predictions got closer to what the actual experiments at the LHC are seeing.

The Bottom Line

This paper tells us that the "fuzziness" of the lead nucleus (the neutron skin) is not just a tiny detail to be ignored. It acts like a subtle steering wheel that changes the shape of the initial crash. This change ripples through the entire explosion, affecting how the fireball expands and how the light it emits flows.

For scientists trying to understand the fundamental laws of the universe, knowing exactly how "hairy" these atomic balls are helps them build better maps of the Quark-Gluon Plasma, the state of matter that existed just moments after the Big Bang.

Technical Summary

1. Problem Statement

In neutron-rich heavy nuclei like Lead-208 ($^{208}\text{Pb}$), the spatial distribution of neutrons extends further than that of protons, creating a "neutron skin." The thickness of this skin ($\Delta_{np}$), defined as the difference between the root-mean-square radii of neutron and proton distributions, is a critical observable for constraining the nuclear symmetry energy and the equation of state (EoS) of nuclear matter.

While the PREX collaboration recently measured this thickness ($\approx 0.28$ fm), its impact on relativistic heavy-ion collisions remains a subject of investigation. Standard models often assume a single nucleon density distribution, ignoring the distinct proton and neutron profiles. This paper addresses the gap in understanding how the neutron skin effect (NSE) influences the initial geometry, space-time evolution of the Quark-Gluon Plasma (QGP), and final-state observables in Pb+Pb collisions at LHC energies (2.76 A TeV).

2. Methodology

The authors employ a hybrid theoretical framework combining the Glauber model for initial state generation and (2+1)D ideal relativistic hydrodynamics for the evolution of the medium.

• Initial State Modeling (Glauber Model):

  • Case A (w/o NSE): Uses a standard single Woods-Saxon distribution with common radius ($d_{Pb} = 6.624$ fm) and skin thickness ($a_{Pb} = 0.549$ fm) for the total nuclear density.
  • Case B (w NSE): Constructs the initial state using separate Woods-Saxon distributions for protons and neutrons, incorporating distinct radii and surface diffuseness parameters based on PREX data:
    • Protons: $d_p = 6.680$ fm, $a_p = 0.447$ fm.
    • Neutrons: $d_n = 6.70 \pm 0.03$ fm, $a_n = 0.55 \pm 0.03$ fm.
  • The initial energy density is a linear combination of wounded nucleons ($n_{WN}$) and binary collisions ($n_{BC}$), with a weighting factor $\alpha = 0.05$.

• Hydrodynamic Evolution:

  • The MUSIC code is used to simulate the (2+1)D ideal hydrodynamic evolution.
  • Parameters: Initial formation time $\tau_0 = 0.14$ fm/c, initial central energy density $576$ GeV/fm$^3$, and a lattice-based EoS.
  • Freeze-out: Hadrons are produced via the Cooper-Frye formalism at a freeze-out temperature of 145 MeV.
  • Probes: The study calculates charged pion spectra, average transverse momentum ($\langle p_T \rangle$), anisotropic flow ($v_2$), and thermal photon emission (integrated over the entire space-time history).

• Comparison: Results from the "with NSE" and "without NSE" scenarios are compared across different impact parameters (centrality classes) and beam energies (2.76 A TeV at LHC and projected 39 A TeV at FCC).

3. Key Contributions

• Quantification of Geometric Impact: The study demonstrates that while the total number of participants ($N_{WN}$) and binary collisions ($N_{BC}$) are only marginally affected by the neutron skin, the initial spatial eccentricity ($\varepsilon_2$) is significantly altered.
• Centrality Dependence: It establishes that the neutron skin effect is most pronounced in peripheral collisions, where the surface structure of the nucleus plays a dominant role, compared to central collisions.
• Differential Sensitivity of Observables: The paper distinguishes between observables sensitive to the initial geometry (flow) versus those sensitive to the bulk medium properties (spectra), showing that anisotropic flow is a much more sensitive probe of the neutron skin than particle spectra.
• Photon vs. Hadron Sensitivity: It highlights that thermal photons, being emitted throughout the collision evolution, exhibit a stronger sensitivity to the initial spatial anisotropy induced by the neutron skin compared to hadrons, which freeze out at the surface.

4. Key Results

• Initial Spatial Eccentricity ($\varepsilon_2$):

  • The inclusion of the neutron skin leads to a significant increase in initial spatial eccentricity, particularly for peripheral collisions (e.g., impact parameter $b=12$ fm).
  • This difference is largest in the early stages of the fireball evolution (first few fm/c).
  • At higher energies (39 A TeV at FCC), the relative increase in $\varepsilon_2$ is smaller, suggesting reduced sensitivity to NSE at higher beam energies.

• Hydrodynamic Evolution:

  • Temperature: The time evolution of the average temperature $\langle T \rangle$ is essentially insensitive to the neutron skin.
  • Flow Velocity: The average transverse flow velocity $\langle v_T \rangle$ shows a slight enhancement when NSE is included, particularly in peripheral collisions.
  • Expansion: Energy density contours show slightly more expansion in the "without NSE" case, indicating that the neutron skin slightly modifies the transverse expansion dynamics.

• Final State Observables:

  • Spectra: The transverse momentum ($p_T$) spectra for charged pions and thermal photons are largely unaffected by the neutron skin.
  • Average $p_T$: A slight enhancement in $\langle p_T \rangle$ is observed for both pions and protons when NSE is included, more so in peripheral collisions.
  • Anisotropic Flow ($v_2$):
    • Hadrons: The elliptic flow ($v_2$) of charged pions shows a noticeable increase when NSE is included, mirroring the increase in initial eccentricity.
    • Photons: Thermal photons exhibit a substantial enhancement in elliptic flow in peripheral collisions when NSE is considered. This is attributed to the direct mapping of initial spatial anisotropy to photon flow without the dilution effects seen in hadronic freeze-out.

5. Significance and Conclusion

The paper concludes that the neutron skin effect is a non-negligible factor in modeling relativistic heavy-ion collisions, particularly for peripheral events and observables sensitive to initial geometry.

• Theoretical Implications: Ignoring the neutron skin in initial state models leads to an underestimation of initial spatial eccentricity and, consequently, the anisotropic flow of particles.
• Experimental Relevance: The inclusion of NSE helps partially reduce the discrepancy between theoretical hydrodynamic calculations and experimental data for photon elliptic flow in peripheral collisions. While current models still underestimate photon $v_2$, accounting for the neutron skin moves the predictions closer to experimental measurements.
• Future Directions: The authors suggest that a quantitative assessment requires incorporating event-by-event fluctuations and viscous effects. However, the qualitative conclusion—that neutron skin significantly shapes the medium evolution and final state flow—remains robust.

In summary, this work provides a crucial link between nuclear structure physics (neutron skin thickness) and high-energy nuclear collision phenomenology, advocating for the inclusion of separate proton and neutron density profiles in hydrodynamic simulations to achieve higher precision in interpreting LHC data.