⚛️ nuclear experiments

Extraction of Neutron-Skin Parameters in Zr+Zr and Ru+Ru Collisions at $\sqrt{s_{\text{NN}}}$ = 200 GeV

By analyzing STAR experimental data on net-charge yield differences in 200 GeV Ru+Ru and Zr+Zr collisions using the UrQMD model, this study extracts a relative neutron-skin thickness difference of $0.278 \pm 0.018$ fm and demonstrates that initial-state nuclear geometry significantly influences baryon transport observables without requiring a baryon junction hypothesis.

Original authors: Junjian Zhou, Dongsheng Li, Fan Si, Xiujun Li, Yifei Zhang

Published 2026-08-17

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Junjian Zhou, Dongsheng Li, Fan Si, Xiujun Li, Yifei Zhang

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

Imagine the atomic nucleus not as a static marble, but as a bustling, chaotic city made of tiny particles called protons and neutrons. For decades, scientists have tried to map this city using low-energy tools, like taking a long-exposure photograph of a busy street. The result is a blurry image where the fast-moving details get smoothed out, hiding the rapid shifts and unique shapes of the city's layout. But recently, physicists discovered a way to take a "snapshot" with a super-fast camera: smashing two atomic nuclei together at nearly the speed of light. In this split-second collision, the nuclei don't just bounce; they melt into a hot, dense soup of fundamental particles called the Quark-Gluon Plasma (QGP). The way this soup expands and cools depends entirely on the exact shape and density of the two nuclei right before they crash. By studying the debris from these crashes, scientists can reverse-engineer the initial shapes of the nuclei, revealing secrets about how matter is built that low-energy experiments simply can't see.

This is the story of a recent investigation into two very similar atomic "cities": Zirconium-96 (Zr) and Ruthenium-96 (Ru). These two are "isobars," meaning they have the same total number of particles (96), but Zr has 40 protons while Ru has 44. Because they are so similar in size, scientists thought smashing them together would be the perfect control experiment to study the hot soup they create. However, the STAR collaboration at a massive particle collider noticed something strange. When they compared the debris from Zr crashes versus Ru crashes, the numbers didn't match the predictions of standard physics models. Specifically, the way protons and neutrons seemed to "transport" or move through the collision was off, leading to a ratio of baryons (heavy particles) to charge that was higher than anyone expected. This anomaly sparked a wild theory: perhaps a mysterious, invisible structure called a "baryon junction" was carrying the extra weight. But before jumping to exotic conclusions, a team of researchers asked a simpler question: Could the answer just be that the two cities had slightly different shapes all along?

In this study, the researchers used a powerful computer simulation called UrQMD to play detective with the shapes of these nuclei. They focused on a specific feature called the "neutron skin." Imagine the nucleus as a ball of mixed protons and neutrons; usually, the neutrons stick out a little further than the protons, forming a fuzzy "skin." The team suspected that the Zr nucleus might have a significantly thicker neutron skin than the Ru nucleus, a difference that would subtly change how the two cities overlap when they collide. To test this, they didn't just guess; they systematically adjusted the "fuzziness" of the Zr surface in their simulation, running millions of virtual collisions for each setting. They then compared the simulated debris against the real data from the STAR experiment, looking for the perfect match.

The results were striking. The standard models, which assumed both nuclei had identical shapes, failed to reproduce the experimental data. However, when the researchers adjusted the simulation to give the Zr nucleus a thicker neutron skin, the virtual debris suddenly lined up perfectly with the real-world measurements. They found that the best fit occurred when the difference in the neutron-skin thickness between the two systems was 0.278 ± 0.018 fm (femtometers). This small adjustment in the initial geometry was enough to explain the strange "net-charge yield difference" observed in the experiment.

Perhaps even more exciting is what this means for the mysterious "baryon transport" puzzle. The study showed that simply having a thicker neutron skin on the Zr nucleus naturally boosted the ratio of baryons to charge in the simulation, pushing it above the value of 1, just like the real data showed. This suggests that the "baryon junction" hypothesis might not be the only, or even the primary, explanation for the anomaly. Instead, the strange behavior could be a result of conventional nuclear physics—specifically, the subtle differences in how the protons and neutrons are arranged on the surface of the nucleus.

The authors emphasize that while their findings don't completely rule out the existence of exotic new physics, they prove that we must first understand the "baseline" of nuclear structure before claiming to have discovered something new. It's a reminder that in the high-stakes game of particle physics, sometimes the most revolutionary discovery is realizing that the map of the territory was just slightly off. By using high-energy collisions as a precise tool to measure these tiny nuclear shapes, this work bridges the gap between the chaotic world of particle collisions and the quiet, structured world of nuclear physics, offering a new, independent way to understand the fundamental forces that hold our universe together.

Technical Summary: Extraction of Neutron-Skin Parameters in Zr+Zr and Ru+Ru Collisions at $\sqrt{s_{NN}} = 200$ GeV

Problem Statement
Recent measurements by the STAR collaboration in isobaric collisions ( $^{96}$ Ru+ $^{96}$ Ru and $^{96}$ Zr+ $^{96}$ Zr) at $\sqrt{s_{NN}} = 200$ GeV revealed two significant anomalies: a net-charge yield difference per participant ( $\Delta Q/\langle N_{\text{part}} \rangle$ ) and a baryon-to-charge transport ratio ( $R_{B/Q} = \langle B \rangle / \Delta Q \times \Delta Z/A$ ) that exceeds unity ( $R_{B/Q} > 1$ ) in central collisions. Standard transport models utilizing isospin-independent Woods-Saxon (WS) distributions failed to reproduce the magnitude and centrality dependence of these observables. This discrepancy led to the hypothesis of a "baryon junction"—a novel non-perturbative gluonic configuration carrying baryon number—as the primary explanation. However, these models may have overlooked the sensitivity of these observables to the initial nuclear geometry, specifically the potential difference in neutron-skin thickness between the two isobars.

Methodology
The authors employed the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model (version 3.4) to simulate the collisions. To investigate the role of nuclear structure, they modified the standard initialization by incorporating an isospin-dependent WS distribution.

Parameter Scanning: While the $^{96}$ Ru nucleus was initialized with default parameters, the proton surface diffuseness parameter ( $a_{p,\text{Zr}}$ ) of the $^{96}$ Zr nucleus was systematically scanned from 0.523 fm down to 0.203 fm.
Correlation Mapping: The study established a clear, strong correlation between the input parameter $a_{p,\text{Zr}}$ and the relative neutron-skin thickness difference ( $\Delta r_{np}$ ) between the two systems.
$\chi^2$ Analysis: The model results for $\Delta Q/\langle N_{\text{part}} \rangle$ were compared against STAR experimental data across various centrality bins. To ensure robustness against finite model statistics, the authors utilized a specific $\chi^2$ definition that weights deviations solely by experimental uncertainty, while explicitly calculating the uncertainty of the $\chi^2$ value itself via error propagation.
Observables: The study focused on the net-charge yield difference and the derived baryon-to-charge transport ratio, applying kinematic cuts consistent with the STAR experiment ( $|y| < 0.5$ , $0.2 < p_T < 7.0$ GeV/c).

Key Results

Extraction of Neutron-Skin Thickness: By minimizing the $\chi^2$ between the model and experimental data, the authors extracted a relative neutron-skin thickness difference of $0.278 \pm 0.018$ fm between the $^{96}$ Zr and $^{96}$ Ru systems. The optimal fit was achieved with a proton surface diffuseness of $a_{p,\text{Zr}} \approx 0.323$ fm.
Resolution of $\Delta Q$ Discrepancy: The inclusion of this neutron-skin thickness difference allowed the UrQMD model to quantitatively reproduce the STAR experimental data for $\Delta Q/\langle N_{\text{part}} \rangle$ across all centrality bins, resolving the previous failure of the isospin-independent baseline.
Impact on Baryon Transport: The study demonstrated that the neutron-skin effect substantially modifies the baryon-to-charge transport ratio. Incorporating the best-fit neutron-skin thickness ( $\Delta r_{np} \approx 0.274$ fm) significantly enhanced the calculated $R_{B/Q}$ values, pushing them above unity in central to mid-central collisions (0–40%). This adjustment brought the model into much closer agreement with the observed $R_{B/Q} > 1$ behavior, particularly in the 0–10% centrality bin.
Centrality Dependence: The sensitivity of the observables to the neutron-skin thickness was found to be highest in the most central (0–10%) and most peripheral (60–80%) collisions, indicating that the conversion efficiency from initial geometric structure to final-state yields is maximized in these regimes.

Significance and Claims
The paper argues that nuclear structure effects, specifically the neutron-skin thickness difference, contribute non-negligibly to baryon transport observables without invoking the baryon junction hypothesis. The authors claim that:

Conventional Baseline: The observed anomalous $R_{B/Q} > 1$ behavior is partially driven by conventional nuclear structure effects. Therefore, any definitive identification of novel QCD phenomena (such as the baryon junction) must first rigorously disentangle these structural baseline contributions.
High-Energy "Snapshots": The work underscores the utility of relativistic heavy-ion collisions as a precision tool for nuclear structure physics. Unlike low-energy probes that provide time-averaged "blurred images," high-energy collisions offer instantaneous "snapshots" of nucleon distributions, allowing for the extraction of fine details like neutron-skin thickness.
Constraints on EoS: The results provide an independent constraint on the nuclear symmetry energy and the equation of state (EoS) of neutron-rich matter, effectively bridging the gap between relativistic heavy-ion collisions and low-energy nuclear structure physics.

The authors conclude that while their results do not rule out the existence of exotic mechanisms, they emphasize that precise initial-state geometry constraints are a critical prerequisite for exploring novel QCD phenomena in relativistic heavy-ion collisions.

More like this