Disentangling nuclear structure through multiparticle azimuthal correlations in high-energy isobar collisions

This study demonstrates that multiparticle azimuthal correlations in central $^{96}$Ru+$^{96}$Ru and $^{96}$Zr+$^{96}$Zr collisions at 200 GeV are sensitive probes of nuclear deformation and neutron skin thickness, offering a robust method to disentangle nuclear structure effects with minimal dependence on shear viscosity.

The Gist

Imagine two massive, invisible billiard balls smashing into each other at nearly the speed of light. This isn't just a simple crash; it's a cosmic experiment designed to melt the balls into a super-hot, super-dense soup called Quark-Gluon Plasma (QGP). Scientists study this soup to understand how the universe behaved just moments after the Big Bang.

But here's the catch: to understand the soup, you first need to know exactly what the billiard balls looked like before they hit.

The "Isobar" Experiment: Twins with a Secret

The scientists in this paper are studying two specific types of atomic nuclei: Ruthenium-96 (Ru) and Zirconium-96 (Zr).

Think of these two nuclei as identical twins. They have the exact same weight (mass number) and the same number of protons and neutrons combined. If you put them on a scale, they weigh the same. However, they have a secret difference in their internal "DNA" (nuclear structure):

1. Shape: One is slightly squashed like a rugby ball (prolate), while the other is slightly stretched like a peanut (octupole deformation).
2. Skin: One has a slightly thicker "skin" of neutrons on the outside than the other.

The researchers wanted to know: If we smash these "twins" together, can we see the difference in their shapes and skins by looking at the debris?

The "Dance" of Particles

When these nuclei collide, they don't just explode randomly. The particles fly out in a specific pattern, like a dance formation. Physicists call this "flow."

• If the nuclei are round, the particles fly out in a circle.
• If the nuclei are squashed or weirdly shaped, the particles fly out in an oval or a star shape.

The paper focuses on multiparticle correlations. Instead of just watching one particle, they look at how groups of particles (3, 4, or more) dance together. It's like watching a group of people at a party:

• Do they all move in sync?
• Does the shape of the room (the initial collision) force them to move in a specific pattern?

The "Shadow" of the Nucleus

The researchers used a super-computer simulation (called the AMPT model) to crash these "twins" together millions of times. They were looking for shadows.

Just as a shadow reveals the shape of an object blocking the light, the way the particles flow reveals the shape of the nuclei that collided.

• The Finding: They discovered that the "dance patterns" of the particles are extremely sensitive to the shape (deformation) of the nuclei.
• The Sweet Spot: This sensitivity is strongest when the nuclei hit head-on (central collisions). It's like trying to see the texture of a pillow; it's easiest to see when you press right in the middle, rather than grazing the edge.
• The Skin: They also found that the thickness of the neutron "skin" leaves a mark, but it's a bit harder to spot because it often cancels itself out with other factors.

The "Viscosity" Check

One big worry for scientists is: "Does the stickiness (viscosity) of the plasma soup mess up our measurements?"

• Imagine trying to see a shape through honey versus water. If the honey is too thick, you can't see the shape clearly.
• The paper's great news: The "dance patterns" they measured are very stable. Even if the plasma is a bit stickier or a bit runnier, the signal of the nuclear shape remains clear. This means the experiment is robust and reliable.

Why Does This Matter?

This study is like a decoder ring for future experiments.

1. Proof of Concept: It proves that we can use high-energy collisions to measure the internal structure of atomic nuclei with incredible precision.
2. Future Guidance: It tells experimentalists at the STAR experiment (at the Relativistic Heavy Ion Collider) exactly what to look for. They can now say, "If we see this specific pattern in the particle dance, it means the nucleus was shaped like a peanut."
3. New Physics: It opens the door to studying other "twins" in the nuclear family, helping us understand how matter is built from the inside out.

In a nutshell: The paper shows that by watching how particles "dance" after a high-speed crash, we can deduce the hidden shapes and skins of the atomic nuclei that started the crash, even if those nuclei look identical on a scale. It's a new way to take an X-ray of the nucleus using the debris of a cosmic collision.

Technical Summary

1. Problem Statement

High-energy heavy-ion collisions create a Quark-Gluon Plasma (QGP), and the resulting anisotropic flow of particles serves as a primary probe of the QGP's transport properties. However, the initial state of these collisions is determined by the spatial distribution of nucleons within the colliding nuclei. While the "flow paradigm" is well-established, there is a need to systematically disentangle the specific effects of nuclear structure (such as quadrupole deformation $\beta_2$, octupole deformation $\beta_3$, neutron skin thickness, and nuclear radius) from the medium's hydrodynamic response.

The isobar collision program at the Relativistic Heavy Ion Collider (RHIC), specifically comparing $^{96}\text{Ru}+^{96}\text{Ru}$ and $^{96}\text{Zr}+^{96}\text{Zr}$ at $\sqrt{s_{NN}} = 200$ GeV, offers a unique opportunity. These nuclei have the same mass number ($A=96$) but different proton numbers ($Z=44$ vs. $Z=40$) and distinct nuclear structures. While previous studies focused on bulk observables, the sensitivity of multiparticle azimuthal correlations (higher-order flow fluctuations) to these specific nuclear structural differences had not been systematically assessed.

2. Methodology

The authors employed the AMPT (A Multi-Phase Transport) model to simulate the collisions. The study utilized a "string-melting" version of AMPT (v2.26t5) with the following key methodological components:

• Initial Conditions: The spatial distribution of nucleons was sampled using a deformed Woods-Saxon (WS) density profile. The model incorporated specific nuclear structure parameters for Ru and Zr, including:
  • Half-density radius ($R_0$)
  • Surface diffuseness/skin depth ($a_0$)
  • Quadrupole deformation ($\beta_2$)
  • Octupole deformation ($\beta_3$)
• Parameter Variation: Five distinct simulation cases were defined (Table I) to isolate the impact of individual parameters:
  • Case 1: Full $^{96}\text{Ru}$ parameters.
  • Case 5: Full $^{96}\text{Zr}$ parameters.
  • Cases 2-4: Systematic variations of $\beta_2$, $\beta_3$, $a_0$, and $R_0$ to isolate their specific contributions.
• Shear Viscosity Checks: The study varied the partonic cross-section ($\sigma$) to simulate different shear viscosity to entropy density ratios ($\eta/s$), ranging from 0.087 to 0.387, to test the robustness of the results against medium properties.
• Observables: The analysis focused on multiparticle cumulants to suppress non-flow effects (e.g., resonance decays, jets). Key observables included:
  • Asymmetric Cumulants: $asc_{nm, n+m}\{3\}$ (3-particle correlations).
  • Non-linear Coupling Coefficients: $\chi_n$ (e.g., $\chi_4, \chi_5$).
  • Symmetric Cumulants: $sc_{n,m}\{4\}$ and Normalized Symmetric Cumulants ($nsc_{n,m}\{4\}$).
  • Pearson Correlation Coefficients: $\rho(v_n^2, [p_T])$ between flow harmonics and mean transverse momentum.
• Statistical Power: Simulations generated over 300 million minimum-bias events per configuration to minimize statistical uncertainties in higher-order correlations.

3. Key Contributions

• Systematic Disentanglement: The paper provides the first systematic breakdown of how specific nuclear structure parameters ($\beta_2, \beta_3, a_0, R_0$) influence multiparticle flow correlations in isobar collisions.
• Isobar Ratio Analysis: By calculating ratios of observables between Ru and Zr ($O_{Ru}/O_{Zr}$), the study demonstrates how final-state hydrodynamic effects largely cancel out, leaving the differences attributable directly to initial-state nuclear structure.
• Viscosity Independence: The authors establish that multiparticle flow ratios are largely insensitive to the shear viscosity ($\eta/s$) of the QGP, validating their use as robust probes of nuclear structure.
• Prediction for STAR: The study offers quantitative predictions for the STAR experiment, specifically regarding the sensitivity of normalized cumulants and flow-momentum correlations to octupole and quadrupole deformations.

4. Key Results

• Sensitivity to Deformation:
  • Octupole Deformation ($\beta_3$): A significant difference in $\beta_3$ between Ru and Zr leads to a ~50% suppression in the symmetric cumulant $sc_{2,3}\{4\}$ in central collisions.
  • Quadrupole Deformation ($\beta_2$): The larger $\beta_2$ in Ru enhances the $sc_{2,4}\{4\}$ ratio in central collisions, while $\beta_3$ induces an opposing effect.
  • Pearson Correlations: The correlation $\rho(v_2^2, [p_T])$ decreases steadily in Ru (reflecting its large prolate deformation), while $\rho(v_3^2, [p_T])$ increases in Zr (reflecting its large octupole deformation).
• Centrality Dependence:
  • Multiparticle correlations are most sensitive to nuclear structure in the most central (0–10%) collision events.
  • The sensitivity diminishes in mid-central and peripheral collisions.
• Neutron Skin and Radius:
  • Differences in neutron skin thickness ($a_0$) and nuclear radius ($R_0$) generally have compensatory effects, often canceling each other out in the ratios, making deformation parameters the dominant drivers of the observed differences.
• Viscosity Robustness:
  • The ratios of multiparticle observables (e.g., $asc_{22|4}$, $asc_{23|5}$) remain nearly constant (deviating from unity by less than ~7%) across different $\eta/s$ values. This confirms that these ratios are insensitive to the QGP's transport properties and are primarily driven by initial-state geometry.

5. Significance

• New Probe for Nuclear Structure: The study confirms that multiparticle azimuthal correlations are a powerful, high-precision tool for measuring nuclear deformation parameters (specifically $\beta_2$ and $\beta_3$) and neutron skin thickness, complementing traditional low-energy nuclear physics methods.
• Experimental Guidance: The results provide a critical theoretical baseline for the STAR collaboration's ongoing analysis of isobar data. The prediction that normalized cumulants ($nsc$) are sensitive to nuclear structure offers a clear target for experimental verification.
• Decoupling Initial and Final State: By demonstrating that isobar ratios cancel out final-state medium effects, the paper validates a methodology to study initial-state fluctuations without the confounding variable of QGP viscosity.
• Future Applications: The framework established here can be extended to other isobar pairs (e.g., $^{16}\text{O}+^{16}\text{O}$, $^{20}\text{Ne}+^{20}\text{Ne}$) to further constrain nuclear parameters and test ab initio nuclear many-body theories.

In conclusion, this work successfully demonstrates that multiparticle flow observables in isobar collisions are highly sensitive to the intrinsic nuclear structure of the colliding nuclei, offering a novel pathway to disentangle initial-state geometry from the dynamical evolution of the QGP.