⚛️ nuclear experiments

Scaling approach to rigid and soft nuclear deformation through flow fluctuations in high-energy nuclear collisions

This paper proposes a scaling approach using triangular flow fluctuations in relativistic $^{238}$ U+ $^{238}$ U collisions to extract the mean and variance of octupole deformation, thereby distinguishing between static and vibrational nuclear deformation origins and refining quark-gluon plasma initial conditions.

Original authors: Lumeng Liu, Chunjian Zhang, Jinhui Chen, Jiangyong Jia, Xu-Guang Huang, Yu-Gang Ma

Published 2026-07-13

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Lumeng Liu, Chunjian Zhang, Jinhui Chen, Jiangyong Jia, Xu-Guang Huang, Yu-Gang Ma

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 atomic nuclei not as smooth, perfect spheres, but as squishy, shape-shifting blobs. For decades, physicists have been trying to figure out the exact "personality" of a heavy nucleus called Uranium-238. Specifically, they are looking at its "pear-shaped" wobble, known as octupole deformation. The big mystery is: Is this pear shape a permanent, rigid feature of the nucleus (like a hard plastic toy), or is it a soft, vibrating wobble that changes from moment to moment (like a jellyfish jiggling in water)?

Until now, the tools used to look at these nuclei were like trying to guess the shape of a spinning top just by looking at its shadow. They could tell the nucleus was deformed, but they couldn't tell if that deformation was a solid, unchanging fact or a fluctuating vibration.

The New "Flow" Microscope
In this paper, the authors propose a clever new way to solve this puzzle using high-energy nuclear collisions. They suggest smashing Uranium-238 nuclei together at nearly the speed of light. When these nuclei collide, they create a super-hot, super-dense soup called a quark-gluon plasma (QGP). As this soup expands and cools, it flows out like water from a squeezed balloon.

The shape of the initial collision determines the direction of this flow. If the nuclei are pear-shaped, the resulting flow will have a specific "triangular" pattern. The authors realized that by looking at the fluctuations in this flow—how much the triangular pattern changes from one collision to the next—they could peek inside the nucleus's shape secrets.

The "Cumulant" Detective Work
To do this, the team used a mathematical tool called "multi-particle cumulants." Think of this like a detective using different levels of magnification:

The Two-Particle Clue (The "Average"): They first looked at the standard triangular flow. Their simulations showed this clue is sensitive to the average amount of pear-shape. It's like measuring the average size of a crowd; it tells you how big the crowd is, but not how much the people are jiggling around.
The Four-Particle Clue (The "Wiggle"): Then, they looked at a more complex four-particle correlation. This is the real magic. Their simulations suggest that this specific measurement is sensitive to the variance or the "jiggle" of the shape. It's like measuring not just the crowd size, but how wildly the people are dancing.

The Simulation Results
The authors ran massive computer simulations (using a model called TRENTo) with 2 × 10⁸ (200 million) collision events for each scenario to get a clear picture. They tested two extreme cases:

The "Rigid" Case: A nucleus with a fixed, solid pear shape (mean deformation $\bar{\beta}_3 = 0.10$ , variance $\sigma_{\beta_3} = 0.000$ ).
The "Soft" Case: A nucleus with no average pear shape, but a lot of random, soft vibrations (mean deformation $\bar{\beta}_3 = 0.00$ , variance $\sigma_{\beta_3} = 0.100$ ).

Both cases produced the same "average" triangular flow signal, which is why previous experiments couldn't tell them apart. However, the four-particle signal was completely different. The simulations showed that the strength of this four-particle signal scales linearly with the fourth-order moment of the deformation, $\langle \beta_3^4 \rangle$ .

The Big Breakthrough (In Simulation)
By combining the "average" clue (two-particle) with the "jiggle" clue (four-particle), the authors demonstrate that you can mathematically separate the mean shape from the fluctuation.

If the nucleus is rigid, the four-particle signal will be low because the shape doesn't change.
If the nucleus is soft/vibrating, the four-particle signal will be high because the shape is constantly shifting.

They found that for Uranium-238, the existing data suggests a small but non-zero average octupole deformation, with a value around $\beta_{3,U} \sim 0.08 - 0.1$ . However, the paper explicitly states that the origin of this deformation (static vs. dynamic) remains an open question because previous measurements only saw the average.

What This Means
This work doesn't claim to have solved the mystery yet; rather, it provides a new, powerful tool to solve it. The authors suggest that by applying this scaling approach to real experimental data from the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), scientists can finally distinguish between a rigid nuclear pear and a soft, vibrating one.

They also note that this method assumes the vibrations follow a "Gaussian" (bell-curve) distribution. If the vibrations are weird or non-Gaussian (skewed or having heavy tails), they would need even more complex measurements (like six- or eight-particle correlations) to figure it out. But for now, this new scaling approach offers a clear path to refining our understanding of nuclear shapes and the initial conditions of the quark-gluon plasma, turning high-energy collisions into a high-resolution microscope for the atomic nucleus.

Technical Summary: Scaling Approach to Rigid and Soft Nuclear Deformation

Problem Statement
The fundamental nature of octupole deformation in atomic nuclei—specifically whether it manifests as a static, rigid deformation of the ground state or as soft, vibrational modes fluctuating around a mean shape—remains an open question in nuclear physics. While low-energy probes (e.g., excited state parities and transition rates) have provided evidence for non-zero octupole correlations in heavy nuclei like $^{238}\text{U}$ , they struggle to distinguish between these two scenarios because both can produce qualitatively similar spectra. In high-energy nuclear collisions, previous analyses of two- and three-particle correlations ( $\langle v_3^2 \rangle$ and $\langle v_3^2 \delta p_T \rangle$ ) in relativistic $^{238}\text{U}+^{238}\text{U}$ collisions successfully constrained the mean square octupole deformation $\langle \beta_3^2 \rangle$ . However, this observable is a sum of the mean deformation squared and its variance ( $\bar{\beta}_3^2 + \sigma_{\beta_3}^2$ ), rendering it insensitive to the critical distinction between a narrow distribution (static/rigid) and a wide distribution (dynamic/vibrational).

Methodology
The authors propose a scaling approach utilizing multi-particle cumulants of triangular flow fluctuations to probe the higher-order moments of the nuclear deformation distribution. The methodology involves:

Analytical Framework: Deriving relations connecting the second- and fourth-order cumulants of the anisotropic eccentricity ( $\varepsilon_n\{2\}$ $ε_{n} {2}$ and $\varepsilon_n\{4\}$ $ε_{n} {4}$ ) to the second- and fourth-order moments of the deformation parameter $\beta_n$ $β_{n}$ ( $\langle \beta_n^2 \rangle$ $⟨ β_{n}^{2} ⟩$ and $\langle \beta_n^4 \rangle$ $⟨ β_{n}^{4} ⟩$ ).
- The second-order cumulant $c_{n,\varepsilon}\{2\}$ is shown to depend linearly on $\langle \beta_n^2 \rangle$ .
- The fourth-order cumulant $c_{n,\varepsilon}\{4\}$ , which suppresses flow fluctuations, is derived to depend on a combination of $\langle \beta_n^4 \rangle$ and $\langle \beta_n^2 \rangle^2$ .
Scaling Hypothesis: The authors hypothesize that for a fixed second-order moment $\langle \beta_3^2 \rangle$ , the magnitude of the fourth-order cumulant $|c_{3,\varepsilon}\{4\}|$ and the ratio $|c_{3,\varepsilon}\{4\}/c_{3,\varepsilon}^2\{2\}|$ scale linearly with the fourth-order moment $\langle \beta_3^4 \rangle$ .
Model Simulations: These relations are validated using the TRENTo initial condition model for $^{238}\text{U}+^{238}\text{U}$ collisions at $\sqrt{s_{NN}} = 193$ GeV. The simulations systematically vary the mean deformation $\bar{\beta}_3$ and variance $\sigma_{\beta_3}$ while keeping $\langle \beta_3^2 \rangle$ fixed at 0.01. A reference system of $^{197}\text{Au}+^{197}\text{Au}$ is used to calculate double ratios ( $R_O$ ) to minimize final-state hydrodynamic response effects.
Statistical Requirements: To obtain stable estimates of the small four-particle cumulants in ultra-central collisions, the study generates $2 \times 10^8$ minimum-bias events per deformation setting.

Key Contributions and Results

Linear Scaling Relations: The simulations confirm that while the two-particle cumulant $c_{3,\varepsilon}\{2\}$ is insensitive to $\langle \beta_3^4 \rangle$ when $\langle \beta_3^2 \rangle$ is fixed, the four-particle cumulant $|c_{3,\varepsilon}\{4\}|$ exhibits a clear, approximately linear dependence on $\langle \beta_3^4 \rangle$ .
Robust Observables: The study demonstrates that the normalized ratio $|c_{3,\varepsilon}\{4\}/c_{3,\varepsilon}^2\{2\}|$ and the system ratio $|R_{c_{3,\varepsilon}\{4\}}|$ (comparing U+U to Au+Au) also scale linearly with $\langle \beta_3^4 \rangle$ . These ratios effectively cancel out linear response coefficients, isolating the initial-state geometric fluctuations.
Extraction of Parameters: By combining the established linear relation of $c_{3,\varepsilon}\{2\}$ $c_{3, ε} {2}$ to $\langle \beta_3^2 \rangle$ $⟨ β_{3}^{2} ⟩$ with the new linear relation of the fourth-order cumulant to $\langle \beta_3^4 \rangle$ $⟨ β_{3}^{4} ⟩$ , the authors show that both the mean octupole deformation ( $\bar{\beta}_3$ $\overset{ˉ}{β}_{3}$ ) and its variance ( $\sigma_{\beta_3}$ $σ_{β_{3}}$ ) can be simultaneously extracted.
- Under the assumption of Gaussian fluctuations, $\bar{\beta}_3$ and $\sigma_{\beta_3}$ are uniquely determined by $\langle \beta_3^2 \rangle$ and $\langle \beta_3^4 \rangle$ via analytical expressions derived in the paper.
Differentiation of Modes: The results indicate that "rigid" (static) and "soft" (vibrational) deformation scenarios, which yield identical $\langle \beta_3^2 \rangle$ values, produce distinct values for $\langle \beta_3^4 \rangle$ . Consequently, this approach provides a direct method to discriminate between static and dynamic origins of octupole collectivity.
Extension to Non-Gaussianity: The paper briefly outlines how higher-order cumulants (sixth and eighth order) could, in principle, provide sensitivity to non-Gaussian features (skewness and kurtosis) of the deformation distribution, though a full extraction would require additional constraints or assumptions regarding response coefficients.

Significance
The paper claims to establish a new tool for probing static and dynamic collective modes in high-energy nuclear collisions. By linking low-energy nuclear structure properties (octupole deformation nature) to high-energy flow fluctuations, the work advances the refinement of initial conditions for Quark-Gluon Plasma (QGP) studies. The authors assert that applying this scaling approach to existing high-energy $^{238}\text{U}$ collision data could provide the first direct experimental evidence to distinguish the nature of the octupole collective mode, thereby resolving a fundamental question in nuclear structure through a unique high-energy collision experiment. The study is presented as an initial-state sensitivity study, with the authors noting that full quantitative validation via event-by-event viscous hydrodynamic simulations remains a task for future work due to computational demands.

Technical Summary: Scaling Approach to Rigid and Soft Nuclear Deformation

More like this