Investigation of Nonlinear Collective Dynamics in Relativistic Heavy-Ion Collisions Using A Multi-Phase Transport Model

Using the AMPT model, this study demonstrates that while the nonlinear response coefficient $\chi_{4,22}$ grows dynamically during the evolution of relativistic heavy-ion collisions, the ratio of this coefficient between U+U and Au+Au systems remains stable across all stages, effectively isolating intrinsic initial-state geometric correlations to support experimental efforts in extracting high-order nuclear structure.

The Gist

Imagine you are trying to figure out the shape of a mysterious, invisible object by smashing two of them together at nearly the speed of light and watching how the debris flies apart. This is essentially what scientists do when they smash heavy atomic nuclei (like Uranium or Gold) together in particle accelerators to create a "soup" of fundamental particles called Quark-Gluon Plasma (QGP).

This paper is a detective story about how scientists can use the patterns of that flying debris to learn about the hidden shapes of the original nuclei, specifically looking for a specific type of squishiness called hexadecapole deformation.

Here is the breakdown of the paper using simple analogies:

1. The Setup: Smashing Squishy Balls

Think of the atomic nuclei as two giant, squishy clay balls.

• Gold (Au) is roughly a perfect sphere, like a smooth marble.
• Uranium (U) is not a perfect sphere; it's shaped like a football (or a rugby ball) with some extra bumps and curves.

When these balls collide, they don't just bounce off; they melt into a super-hot, super-dense fluid (the QGP) that expands outward like an exploding balloon. The way this "balloon" expands depends on the shape of the original clay balls.

2. The Problem: The "Echo" is Too Loud

Scientists measure how the debris flies out in different directions. They look for a specific pattern called flow.

• Linear Flow: If you squish a round ball, it flies out in a simple oval shape. This is easy to predict.
• Non-Linear Flow: If you squish a weirdly shaped ball, the debris doesn't just follow the simple oval; it creates complex, wavy patterns.

The paper focuses on a specific measurement called $\chi_{4,22}$. Think of this as a "complexity meter." It measures how much the debris pattern is influenced by the shape of the collision versus just the random chaos of the smash.

The Confusion:
For a long time, scientists thought this "complexity meter" only told them about how the fluid (the QGP) reacted to the smash. They thought it was like an echo: You hit the wall, the echo comes back. They assumed the echo didn't tell you anything about the wall's original texture, only about the air in the room.

3. The Discovery: The Echo Does Remember the Wall

The authors used a super-computer simulation (called the AMPT model) to watch the collision happen in slow motion, stage by stage. They looked at three moments:

1. The Partonic Phase: The moment the nuclei first melt into the fluid soup.
2. The Coalescence Phase: The moment the fluid starts cooling and clumping back together into particles.
3. The Final Phase: The moment the particles fly out to be detected.

What they found:

• The Absolute Value: As the collision progresses from soup to solid particles, the "complexity meter" ($\chi_{4,22}$) gets bigger and bigger. This confirms that the fluid does act like a dynamic amplifier, building up complexity as it expands.
• The Ratio Trick: Here is the clever part. The authors compared the Uranium smash (weird shape) to the Gold smash (round shape).
  • They calculated the ratio: (Complexity of Uranium) divided by (Complexity of Gold).
  • The Magic: Even though the "complexity meter" changed wildly as the collision evolved, the ratio between Uranium and Gold stayed exactly the same at every single stage.

4. The Analogy: The Amplifier and the Microphone

Imagine you have two microphones:

• Mic A is attached to a round drum (Gold).
• Mic B is attached to a weirdly shaped drum (Uranium).

You hit both drums. The sound gets louder as it travels through a giant, echoing hall (the expanding fluid).

• If you look at the volume of Mic B, it changes a lot depending on where you stand in the hall (the stage of evolution).
• But if you take the ratio of Mic B's volume to Mic A's volume, that number stays constant, no matter where you stand in the hall.

Why? Because the hall (the fluid) amplifies both sounds equally. The difference in volume between the two mics is purely due to the shape of the drums, not the hall.

5. The Conclusion: Why This Matters

This paper proves that by taking this ratio, scientists can cancel out all the messy, complicated physics of the fluid expansion (the "hall"). What is left is a clean, pure signal of the original shape of the Uranium nucleus.

• The Result: This gives scientists a powerful new tool to measure the "hexadecapole deformation" (a specific type of 16-sided squishiness) of atomic nuclei.
• The Impact: It's like finally being able to see the fingerprint of the nucleus clearly, even though it was hidden inside a chaotic explosion. This helps us understand the fundamental structure of matter in the universe.

In a nutshell: The authors showed that while the "soup" created in these collisions changes and grows, the difference in behavior between a round nucleus and a weird-shaped nucleus remains a constant, reliable fingerprint. By comparing the two, we can ignore the soup and read the shape of the nucleus directly.

Technical Summary

1. Problem Statement

Relativistic heavy-ion collisions create a Quark-Gluon Plasma (QGP), a state of deconfined matter. A primary challenge in this field is that experimental observables are measured only in the final state, while the QGP exists only in the early stages of the collision ($\sim 0.5$ fm/$c$). Understanding the QGP requires mapping initial-state geometric configurations to final-state momentum anisotropies.

While linear relationships between initial spatial eccentricities ($\epsilon_n$) and final flow harmonics ($v_n$) hold well for $n=2,3$, higher-order harmonics (e.g., $v_4$) involve significant nonlinear response mechanisms. Specifically, $v_4$ receives large contributions from the square of the elliptic flow ($v_2^2$). The nonlinear response coefficient, $\chi_{4,22}$, quantifies this coupling.

• The Conflict: Traditionally, $\chi_{4,22}$ was viewed solely as a medium response property, independent of the initial state. However, recent studies suggest it is sensitive to the hexadecapole deformation ($\beta_4$) of colliding nuclei (e.g., $^{238}$U).
• The Challenge: Extracting intrinsic nuclear structure (like $\beta_4$) is complicated by theoretical uncertainties arising from the complex dynamical evolution of the medium (viscosity, hadronic rescattering, etc.). It is unclear how much of the observed $\chi_{4,22}$ is generated dynamically versus how much reflects the initial geometry.

2. Methodology

The authors utilize the A Multi-Phase Transport (AMPT) model to simulate $^{238}$U+$^{238}$U and $^{197}$Au+$^{197}$Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The study employs a "stage-by-stage" analysis to disentangle the contributions of different evolutionary phases.

• Initial Geometry: Nuclei are modeled using deformed Woods-Saxon distributions.
  • Uranium ($^{238}$U): Includes quadrupole ($\beta_2 = 0.286$) and hexadecapole ($\beta_4 = 0.100$) deformations.
  • Gold ($^{197}$Au): Treated as approximately spherical with small deformations ($\beta_2 = -0.131, \beta_4 = 0.031$).
• Simulation Stages: The AMPT model (string-melting version) tracks the system through three distinct snapshots:
  1. Partonic Phase: After initial parton scatterings (ZPC model) but before hadronization.
  2. Hadronization: Immediately after quark coalescence, effectively switching off subsequent hadronic rescattering ($t_{max} = 0.6$ fm/$c$).
  3. Full Evolution: After final-state hadronic rescattering and kinetic freeze-out ($t_{max} = 30$ fm/$c$).
• Observables:
  • Nonlinear Response Coefficient ($\chi_{4,22}$): Defined as $\chi_{4,22} \equiv v_4\{\Phi_2\} / \langle v_2^4 \rangle^{1/2}$, where $v_4\{\Phi_2\}$ is the projection of $v_4$ onto the second-order event plane.
  • Ratio Observable ($R(X)$): The ratio of an observable $X$ between U+U and Au+Au systems ($R(X) = X_{UU} / X_{AuAu}$) is used to cancel out common medium evolution effects.

3. Key Contributions

• Microscopic Origin Mapping: The paper provides the first systematic, stage-resolved investigation of how $\chi_{4,22}$ is generated and modified from the partonic phase through hadronization to the final hadronic rescattering phase.
• Validation of Ratio Method: It theoretically validates the use of the U+U/Au+Au ratio as a robust experimental tool to isolate initial-state nuclear deformations ($\beta_4$) by canceling out complex medium amplification dynamics.
• Decoupling Dynamics from Geometry: The study demonstrates that while the absolute magnitude of flow observables is heavily dependent on the medium's evolution, the relative geometric correlations between different nuclei are "locked in" at the initial state.

4. Key Results

• Dynamic Generation of $\chi_{4,22}$:

  • The absolute magnitude of $\chi_{4,22}$ increases monotonically as the system evolves from the partonic phase to the final kinetic freeze-out.
  • This confirms that $\chi_{4,22}$ is a dynamically generated medium response, accumulating continuously through collective expansion.
  • The constituent components (asymmetric cumulant $ac_2\{3\}$ and elliptic flow $v_2$) also increase during evolution, while the four-particle cumulant $\langle v_2^4 \rangle$ remains relatively stable during the hadronization phase.

• Stability of the Ratio $R(\chi_{4,22})$:

  • Despite the significant evolution of absolute values, the ratio $R(\chi_{4,22}) = \chi_{4,22}^{UU} / \chi_{4,22}^{AuAu}$ remains stable across all evolutionary stages within statistical uncertainties.
  • This indicates that the ratio effectively cancels the "medium amplification efficiency," isolating the intrinsic geometric differences (specifically the hexadecapole deformation $\beta_4$) present at the initial state.

• Flow Harmonic Behavior:

  • Hadronization: The elliptic flow ($v_2$) remains stable, while hexadecapole flow ($v_4$) undergoes modification.
  • Hadronic Rescattering: $v_4$ remains largely unchanged, while $v_2$ increases substantially.
  • Event Plane Correlation: The angular correlation $\langle \cos 4(\Phi_4 - \Phi_2) \rangle$ increases monotonically throughout the evolution, driving the growth of $\chi_{4,22}$.

5. Significance

• Experimental Guidance: The findings provide strong theoretical support for recent experimental efforts (e.g., by the STAR collaboration) to extract high-order nuclear structure parameters, such as the hexadecapole deformation of Uranium, using nonlinear flow observables.
• Theoretical Insight: The study clarifies that while nonlinear flow is a complex dynamical process, the relative differences between collision systems are robust against medium uncertainties. This justifies the use of comparative ratio observables in heavy-ion physics to probe nuclear structure.
• Future Applications: The methodology suggests that applying this comparative ratio approach to high-statistics experimental data, while carefully mitigating non-flow effects, will allow for precise resolution of the high-order spatial structure of deformed nuclei.