Nonlinear response of flow harmonics in Gubser flow with participant-reaction planes mismatch

This paper extends perturbative Gubser flow solutions to derive analytic nonlinear response relations for flow harmonics $v_2$ and $v_4$, revealing that a mismatch between participant and reaction planes introduces a critical angular factor that can significantly alter the magnitude and even the sign of effective nonlinear response coefficients.

The Gist

Imagine two heavy, lumpy balls (atomic nuclei) smashing into each other at nearly the speed of light. When they collide, they don't just shatter; they melt into a tiny, super-hot drop of "perfect fluid" called Quark-Gluon Plasma (QGP). This fluid expands and cools down incredibly fast, eventually turning back into a spray of particles that fly out in all directions.

Physicists want to understand the shape of this explosion. They look at the particles flying out and ask: "Did they fly out evenly in a circle, or did they prefer certain directions?"

The Main Characters: Eccentricity and Flow

1. The Shape of the Collision (Eccentricity):
   Imagine the two nuclei are like lumpy potatoes. When they collide, they don't always hit perfectly head-on. Sometimes they graze each other, creating a football-shaped collision zone. Sometimes they are more oval, or even have weird bumps.

   • $\epsilon_2$ (Epsilon-2): Think of this as how "oval" the collision is (like a rugby ball).
   • $\epsilon_4$ (Epsilon-4): Think of this as how "four-leaf clover" or "squared-off" the shape is.

2. The Flow of Particles ($v_2$ and $v_4$):
   As the hot fluid expands, the pressure pushes harder in the short direction of the oval than the long direction. This pushes the particles out faster in that direction.

   • $v_2$: How much the particles prefer to fly out in an oval pattern.
   • $v_4$: How much they prefer a four-leaf clover pattern.

The Big Question: How are they connected?

For a long time, scientists thought the relationship was simple:

• If you have a little oval shape ($\epsilon_2$), you get a little oval flow ($v_2$).
• If you have a four-leaf shape ($\epsilon_4$), you get a four-leaf flow ($v_4$).

But, it turns out it's more complicated. The fluid is so dynamic that the oval shape ($\epsilon_2$) can actually create a four-leaf flow ($v_4$) on its own, even if the initial shape wasn't a four-leaf clover! It's like if you squeeze a round balloon into an oval, the air inside might swirl in a way that creates a four-pointed star pattern on the surface.

This paper tries to write down the exact math for this "mixing" effect.

The Twist: The "Misaligned Compass"

Here is the most interesting part of this paper.

Imagine you are trying to measure the shape of a spinning top.

• The Participant Plane: This is the actual shape of the collision at the moment of impact (where the lumps of the potatoes hit).
• The Reaction Plane: This is the direction the particles actually fly out in the end.

In a perfect world, these two would line up perfectly. But in reality, because the collision is chaotic and happens in a split second, the "actual shape" and the "final flight direction" are slightly rotated relative to each other. It's like trying to draw a map of a city based on a photo taken from a helicopter that is slightly tilted.

The Paper's Discovery:
The authors found that this tiny rotation (the mismatch) changes the math significantly.

1. The "Volume Knob" Effect: The mismatch acts like a volume knob for the nonlinear effect. Depending on the angle, it can turn the signal up, turn it down, or even flip the sign (make a positive effect negative).
2. The "Noise" Myth: Previously, scientists thought this rotation was just random "static noise" that would average out if you looked at enough collisions. This paper says: No, it's not just noise. It's a fundamental part of the geometry.
3. The Nuclear Shape Detector: Because the rotation depends on the initial shape of the colliding nuclei (whether they are round, football-shaped, or have weird bumps), measuring this flow actually tells us about the internal structure of the atomic nuclei themselves.

A Simple Analogy: The Pizza Dough

Imagine two people throwing pizza dough into the air.

• The Eccentricity ($\epsilon$): How the dough is shaped when it leaves their hands (maybe it's a bit of a square, or a triangle).
• The Flow ($v$): How the dough spins and stretches as it flies.

If the dough is a perfect circle, it spins evenly. If it's a square, it stretches into a cross. But if the dough is a slightly tilted square, the way it spins might look like a diamond or a different shape entirely.

This paper says: "Hey, if you want to know exactly how the dough was shaped when it left the hands, you can't just look at how it spins. You have to account for the fact that the dough was tilted when it was thrown. If you ignore the tilt, you'll get the wrong recipe for the dough's shape."

Why Does This Matter?

1. Understanding the "Perfect Fluid": It helps us understand how this super-hot soup of particles behaves.
2. X-Raying the Nucleus: It gives physicists a new, precise tool to see the internal shape of atomic nuclei. Just like an MRI scan, but using high-speed particle collisions instead of magnets.
3. Fixing the Math: It corrects a long-held assumption that the "tilt" between the start and end of the collision was just random noise. It turns out that tilt holds the key to understanding the initial geometry.

In a nutshell: This paper uses advanced math to show that the way particles fly out of a nuclear collision depends not just on the shape of the collision, but on how that shape is rotated relative to the final explosion. By understanding this rotation, we can learn the secret shapes of the atomic nuclei themselves.

Technical Summary

1. Problem Statement

In relativistic heavy-ion collisions, the initial spatial anisotropy of the Quark-Gluon Plasma (QGP), characterized by eccentricities ($\epsilon_n$), drives the development of anisotropic flow harmonics ($v_n$). While the linear response of $v_2$ and $v_3$ to $\epsilon_2$ and $\epsilon_3$ is well-established, higher-order harmonics like $v_4$ exhibit significant nonlinear responses, primarily driven by the coupling of lower-order modes (e.g., $v_4 \sim v_2^2$).

Two critical gaps in the theoretical understanding of this phenomenon were addressed in this work:

1. Lack of Analytic Derivation: Most studies rely on numerical simulations. There is a scarcity of analytic solutions within relativistic hydrodynamics that explicitly map initial eccentricities ($\epsilon_2, \epsilon_4$) to nonlinear flow harmonics ($v_2, v_4$).
2. Geometric Phase Mismatch: Experimentally, flow harmonics are measured relative to the reaction plane ($\Psi_n$), while initial eccentricities are defined relative to the participant plane ($\psi_n$). These planes often misalign due to event-by-event fluctuations. Standard treatments often treat this mismatch as statistical noise that averages out. However, recent phenomenological studies suggest the nonlinear response coefficient is sensitive to initial nuclear deformation (specifically hexadecapole deformation $\beta_4$), implying the geometric phase mismatch is an intrinsic, non-removable feature of the hydrodynamic mapping that has not been rigorously derived analytically.

2. Methodology

The authors employed the Gubser flow framework, an exact analytic solution to ideal relativistic hydrodynamics, extended to include transverse perturbations.

• Framework: They utilized the Gubser solution, which maps Minkowski spacetime to de-Sitter spacetime via a Weyl transformation. To introduce anisotropy, they applied perturbative solutions around the azimuthally symmetric Gubser background, following the approach of Refs. [18–21].
• Perturbation Setup: The initial entropy density perturbation $S(\theta, \phi)$ was parameterized to include both second-order ($\epsilon_2$) and fourth-order ($\epsilon_4$) modes. Crucially, they introduced distinct participant plane angles ($\psi_2, \psi_4$) into the perturbation function:
  $$S(\theta, \phi) \propto a_2 \cos[2(\phi - \psi_2)] + a_4 \cos[4(\phi - \psi_4)]$$
• Calculation Steps:
  1. Derived the perturbed fluid velocity and temperature profiles in the early-time limit ($\tau \ll L$).
  2. Calculated the eccentricities ($\epsilon_n$) from the perturbed entropy density to establish the relationship between perturbation coefficients ($a_n$) and $\epsilon_n$.
  3. Applied the Cooper-Frye formula on an isothermal freeze-out surface to compute the particle distribution.
  4. Extracted the flow harmonics $v_2$ and $v_4$ by projecting the distribution onto the reaction plane ($\Psi_n$).
  5. Analyzed the results in both small and large transverse momentum ($p_T$) limits.

3. Key Contributions and Results

A. Analytic Nonlinear Response Relations

The study derived explicit analytic relations connecting initial eccentricities to flow harmonics:

• Linear Response: $v_2 = \chi_{22} \epsilon_2$
• Nonlinear Response: $v_4 = \chi_{44} \epsilon_4 + \chi_{42} \epsilon_2^2$
  • Here, $\chi_{ij}$ are response coefficients dependent on the freeze-out temperature ($T_f$), particle mass ($m_0$), and transverse momentum ($p_T$).
  • The term $\chi_{42} \epsilon_2^2$ represents the nonlinear contribution from the $v_2^2$ mode to $v_4$.

B. Universal High-$p_T$ Limit

The authors confirmed the well-known universal behavior for ideal fluids in the large $p_T$ limit:
$$\frac{v_4}{v_2^2} \to \frac{1}{2}$$
This result provides an analytical confirmation of previous numerical findings within the Gubser model.

C. The Impact of Participant-Reaction Plane Mismatch (Main Contribution)

The most significant finding is the rigorous derivation of how the angle difference between the participant plane ($\psi_n$) and the reaction plane ($\Psi_n$) modifies the observed nonlinear response.

• Modified Coefficient: The effective nonlinear response coefficient is not just $\chi_{42}$, but is modulated by a geometric factor:
  $$\chi_{\text{eff}} = \chi_{42} \cos[4(\psi_2 - \psi_4)]$$
• Phase Angle $\xi^*$: The mismatch induces a phase shift $\xi^*$ between the participant plane and the measured reaction plane for the 4th harmonic.
• Sign Reversal: Depending on the value of $4(\psi_2 - \psi_4)$, the effective coefficient can:
  • Remain positive (if $\cos > 0$).
  • Vanish (if $\cos = 0$, i.e., $4(\psi_2 - \psi_4) = \pi/2$).
  • Become negative (if $\cos < 0$).
• Implication: The standard experimental extraction of the nonlinear coefficient (often denoted $\chi_N$) is a convolution of the true medium response and this geometric projection factor. It is not a pure measure of QGP transport properties (like shear viscosity) but is intrinsically sensitive to the initial nuclear geometry.

4. Significance and Implications

1. Reinterpretation of Experimental Data: The paper challenges the conventional view that the nonlinear response coefficient is solely a probe of QGP viscosity. It demonstrates that the coefficient is highly sensitive to the initial spatial configuration of the colliding nuclei.
2. Nuclear Structure Probing: The sensitivity to the angle mismatch $\cos[4(\psi_2 - \psi_4)]$ provides a theoretical basis for using higher-order flow correlations to image nuclear deformation (specifically hexadecapole deformation $\beta_4$). In collisions of deformed nuclei, the correlation between $\psi_2$ and $\psi_4$ is driven by the nuclear shape, meaning the "noise" previously assumed to average out actually carries macroscopic structural information.
3. Analytic Foundation: By providing an exact analytic solution within the Gubser framework, the authors isolate geometric factors from dissipative effects (viscosity), offering a transparent mathematical foundation for understanding collective phenomena.
4. Future Directions: The work suggests that future phenomenological analyses must account for the event-by-event geometric phase mismatch to accurately extract transport coefficients or nuclear structure parameters.

In summary, this paper bridges the gap between analytic hydrodynamics and experimental observables, revealing that the "nonlinear response" is a complex interplay between medium dynamics and the geometric alignment of the initial collision state.