Higher-Order Flow Harmonics in QGP Dynamics

• Higher-order flow harmonics are Fourier coefficients (vₙ for n ≥ 3) that capture multi-particle azimuthal anisotropies, providing clear probes of initial geometry fluctuations and QGP dynamics.
• They are extracted using multi-particle cumulant methods and event-shape engineering to disentangle linear hydrodynamic responses from nonlinear mode-coupling effects.
• Measurement of these harmonics constrains key QGP properties such as specific shear viscosity and validates initial-state modeling in precision QGP tomography.

Higher-order flow harmonics are Fourier coefficients ($v_n$ with $n \geq 3$) quantifying multi-particle azimuthal anisotropies in the momentum distributions of final-state hadrons produced in relativistic heavy-ion collisions. These observables, extending beyond elliptic flow ($v_2$), are sensitive probes of both the event-by-event fluctuating initial geometry and the transport properties of the QCD medium formed in such collisions. Higher-order harmonics are generated through a combination of linear hydrodynamic response to initial spatial eccentricities ($\epsilon_n$) and nonlinear, mode-coupling contributions mediated by the development of collective flow. Their measurement and decomposition provide essential constraints for precision modeling of quark-gluon plasma (QGP) dynamics and the extraction of its specific shear viscosity ($\eta/s$) (Prasad et al., 2 May 2025).

1. Theoretical Foundations and Response Formalism

The azimuthal distribution of produced particles can be decomposed as

$$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$

where $v_n$ are the flow coefficients and $\Psi_n$ are the symmetry (event-plane) angles of order $n$.

• Initial State Geometry: The initial transverse spatial distribution is quantified by participant eccentricities $\epsilon_n$, defined event-by-event as

$n \geq 3$0

with $n \geq 3$1 the initial energy density (Qian et al., 2013).

• Linear and Nonlinear Response: For $n \geq 3$2 (elliptic) and $n \geq 3$3 (triangular) flow, the final coefficients scale linearly with the corresponding eccentricities ($n \geq 3$4, $n \geq 3$5). For $n \geq 3$6, both linear and nonlinear hydrodynamic response terms are relevant:

$n \geq 3$7

where $n \geq 3$8 are linear response coefficients and $n \geq 3$9 encapsulate nonlinear mode-mixing (Prasad et al., 2 May 2025, Qian et al., 2013, Giacalone et al., 2018).

Higher harmonics are thus sensitive to both the initial geometry—including its fluctuations—and the collective nonlinear response of the medium (Prasad et al., 2 May 2025, Qian et al., 2013, Collaboration et al., 2020, Giacalone et al., 2018).

2. Experimental Determination and Multi-particle Methods

The extraction of $v_2$0 relies on suppressing nonflow correlations and resolving event-by-event fluctuations:

• Multi-particle Q-cumulant Method: In each event,

$v_2$1

with $v_2$2 the event multiplicity. Two- and four-particle cumulants are constructed as

$v_2$3

with $v_2$4, $v_2$5 (when $v_2$6) (Prasad et al., 2 May 2025, Sorensen, 2011).

• Suppression of Nonflow: Large pseudorapidity ($v_2$7) gaps and partitioning into subevents further reduce short-range nonflow effects (Prasad et al., 2 May 2025, 2002.04583).
• Mixed-Harmonic Correlators: Three- and four-particle mixed correlators,

$v_2$8

and similar constructs enable separation of linear and nonlinear components (Collaboration et al., 2020, Magdy, 2021, Giacalone et al., 2018).

3. Event Shape Engineering and Spherocity

Event-shape engineering with transverse spherocity ($v_2$9) is a key technique for modulating the relative strength of nonlinear mode-mixing:

• Definition: Spherocity is defined for each event as

$\epsilon_n$0

with $\epsilon_n$1. Low-$\epsilon_n$2 ("jetty") events correspond to pencil-like topologies with large $\epsilon_n$3, while high-$\epsilon_n$4 ("isotropic") events have small $\epsilon_n$5 (Prasad et al., 2 May 2025).

• Isolation of Nonlinear/Linear Response: Selection on $\epsilon_n$6 allows systematic variation of $\epsilon_n$7 at fixed centrality. Nonlinear contributions to higher $\epsilon_n$8 (e.g., $\epsilon_n$9, $\eta/s$0) are maximized in low-$\eta/s$1 classes and suppressed in high-$\eta/s$2, isolating the linear response to initial $\eta/s$3 (Prasad et al., 2 May 2025).

This methodology enables the clean experimental disentanglement of linear and nonlinear contributions in higher-order harmonics (Prasad et al., 2 May 2025, 2002.04583).

4. Scaling Properties, Centrality Dependence, and Mode-Coupling

Higher-order flow harmonics exhibit distinctive scaling and centrality dependence:

• Ordering and Centrality Trends:
  • $\eta/s$4 for integrated events in all but ultra-central collisions (Prasad et al., 2 May 2025, Mohapatra, 2011).
  • $\eta/s$5 increase from central to mid-central events, peaking typically around $\eta/s$6–$\eta/s$7 centrality.
  • The ratios $\eta/s$8 and $\eta/s$9 are sensitive to the nonlinear hydrodynamic response and show strong variation with centrality and event shape (Prasad et al., 2 May 2025, Pandey et al., 6 May 2025).
• Mode-mixing and Nonlinear Response Coefficients:
  • Nonlinear terms such as $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$0 and $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$1 dominate the higher harmonics for $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$2.
  • Mode-coupling coefficients $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$3, extracted via multi-particle correlators, are found to be nearly centrality- and $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$4-independent and insensitive to variations in partonic cross-section or viscosity, supporting their origin in initial-state geometry (Collaboration et al., 2020, Magdy, 2021, Pandey et al., 6 May 2025).
  • In hydrodynamic and transport models, $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$550% of $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$6 in mid-central collisions can originate from nonlinear mode mixing (Qian et al., 2013, Magdy, 2021).
• Empirical Scaling: Ratios $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$7 are empirically flat as functions of $$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$8 and grow gently with system size ($$\frac{dN}{d\phi} = \frac{N}{2\pi} \left[ 1 + 2\sum_{n=1}^{\infty} v_n \cos n(\phi-\Psi_n) \right]$$9), consistent with "acoustic" viscous damping (Lacey et al., 2011, Pandey et al., 6 May 2025).

5. Physical Interpretation and Implications for QGP Properties

The measured higher-order harmonics and their decompositions provide uniquely stringent constraints on initial conditions and QCD transport coefficients:

• Sensitivity to Viscosity:
  • Higher $v_n$0 harmonics are suppressed by shear viscosity with $v_n$1, $v_n$2.
  • Extraction of $v_n$3 from joint fits to $v_n$4, $v_n$5, and $v_n$6 yields values close to the conjectured lower bound, $v_n$7 (Prasad et al., 2 May 2025, Lacey et al., 2011, Lacey, 2011).
  • Nonlinear contributions are particularly sensitive to event-by-event $v_n$8 fluctuations: Spherocity-selected samples with larger $v_n$9 have up to 30–50% higher medium response for $\Psi_n$0 (Prasad et al., 2 May 2025).
• Discrimination of Initial-State Models:
  • The scaling and centrality trends of higher-order $\Psi_n$1, as well as ratios such as $\Psi_n$2 and $\Psi_n$3, serve as powerful discriminants between initial condition models (e.g., Monte Carlo Glauber vs. MC-KLN) (Lacey et al., 2010, Lacey et al., 2011).
  • Accurate extraction of $\Psi_n$4 demands unambiguous identification of the initial geometry, achievable through systematic study of these scaling features (Lacey et al., 2010).
• Event-by-Event Fluctuations and Symmetric Cumulants:
  • Measurements of symmetric cumulants, $\Psi_n$5, and their $\Psi_n$6 dependence, provide information on the interplay of different harmonic fluctuations and constraints on temperature-dependent transport properties (Collaboration, 2017).

6. Model and Data Systematics

Simulations and comparisons with hydrodynamic and transport models reinforce and contextualize experimental findings:

• Hydrodynamics with Fluctuating IC: State-of-the-art event-by-event viscous hydrodynamics employing fluctuating initial conditions quantitatively reproduce $\Psi_n$7 spectra, mode-coupling coefficients, and event-plane correlations across RHIC and LHC energies (Collaboration et al., 2020, Magdy, 2021, Giacalone et al., 2018).
• AMPT and String Melting: Multi-phase transport approaches, e.g., AMPT with string melting, match experimental trends when tuned to reproduce charged particle multiplicity and low-$\Psi_n$8 spectra (Prasad et al., 2 May 2025, Magdy, 2021).
• Role of Event Shape Selection: Event-shape engineering via $\Psi_n$9 or reduced flow vector ($n$0) classes isolates nonlinear contributions, supporting the essentially uncorrelated nature of linear and nonlinear $n$1 components—validating the separation of physics encoded in $n$2 and $n$3 (Prasad et al., 2 May 2025, 2002.04583).

These developments have turned higher-order flow harmonics into a cornerstone of QGP tomography, providing both diagnostic tools for the initial-state geometry and precision handles on the QCD medium's transport properties.

7. Outlook and Precision QGP Tomography

Higher-order flow harmonics have evolved into precision observables underpinning the tomographic capabilities of ultrarelativistic heavy-ion collision experiments. Key avenues include:

• Extension to $n$4, $n$5, $n$6: Nonlinear hydrodynamic frameworks and multi-particle correlator measurements now allow full decomposition of hexagonal and higher harmonics, leveraging the entire hierarchy of mode-coupling terms (Giacalone et al., 2018, Pandey et al., 6 May 2025).
• Differential Observables: $n$7-dependent, species-dependent, and event-class-resolved analyses continue to refine the extraction of QGP properties, such as the temperature dependence of $n$8 and the role of bulk viscosity.
• Unified Description Across Systems: The empirical and model trends are robust from low to top RHIC energy, including baryon-rich ($n$9 few GeV) to highly transparent ($\epsilon_n$0 TeV) regimes (collaboration et al., 2020, Pandit, 2012).
• Constraints on Initial Fluctuations and Collectivity Origin: Mass ordering, number-of-constituent-quark (NCQ) scaling in higher-order $\epsilon_n$1, and the scaling properties of ratios such as $\epsilon_n$2, systematically probe the onset of collectivity, the precision of hydrodynamic modeling, and the character of initial-state fluctuations (Xu et al., 2016, Devi et al., 2023, collaboration et al., 2020).

Collectively, higher-order flow harmonics have established themselves as indispensable messengers of QGP dynamics, uniquely sensitive both to the fluctuating space-time geometry of heavy-ion collisions and to the microscopic properties—such as viscosity—of strongly interacting matter (Prasad et al., 2 May 2025, Lacey et al., 2011, Mohapatra, 2011, Giacalone et al., 2018).