Sources of longitudinal flow decorrelations in high-energy nuclear collisions

Introduction. High-energy nuclear collisions serve as laboratories for generating and examining the QGP under varied conditions Busza et al. (2018). Understanding the dynamics and properties of this exotic state requires precise control of its initial conditions (IC), including the distributions of colliding nucleons, their constituent partons, the deposition of energy, and the 3D geometric profile. Recent advancements showcased the potential to fine-tune the IC by comparing collision systems with similar masses yet distinct shapes, allowing for precise control over the QGP’s initial shape without altering its hydrodynamic response Giacalone et al. (2021); Jia and Zhang (2023); Nijs and van der Schee (2023); Giacalone (2023). Leveraging ratios of bulk observables, such as the elliptic flow, between these systems facilitates the extraction of valuable information about the IC Bally et al. (2022); Xu et al. (2023). Such experimental comparisons, e.g. between ²³⁸U+²³⁸U and ¹⁹⁷Au+¹⁹⁷Au STAR Collaboration (2024) or ⁹⁶Ru+⁹⁶Ru and ⁹⁶Zr+⁹⁶Zr Abdallah et al. (2022); Jia and Zhang (2023), have indicated strong impacts of nuclear structure on the geometrical properties of the IC, which in turn can be used to constrain the longitudinal profile Bhatta et al. (2023); Zhang et al. (2024).

Most investigations of the QGP focused on transverse $xy$ profile of the IC near mid-pseudorapidity ($\eta\approx 0$). However, the energy deposition process by the colliding nucleons is inherently non-boost-invariant and fluctuating along $\eta$, known as longitudinal fluctuations or decorrelations Bozek and Broniowski (2016); Pang et al. (2016); Li and Yan (2020); Franco and Luzum (2020). Consequently, the initial shape of the QGP and its transverse expansion vary with $\eta$ even within a single collision event. In popular QCD-inspired string models (PYTHIA Sjöstrand (2020), HIJING Gyulassy and Wang (1994), AMPT Lin et al. (2005)), where nucleons deposit energy in color flux tubes, longitudinal fluctuations arise from variations in the starting positions and lengths of the flux tubes along $\eta$ Bozek and Broniowski (2016); Pang et al. (2016). This gives rise to a global geometry that varies smoothly with $\eta$ as well as local stochastic fluctuations (see Fig. 1).

We illustrate this concept through an analysis of elliptic flow ($V_{2}=v_{2}e^{2i\Psi}$), which characterizes the elliptical distribution of particles in the $xy$-plane. The $V_{2}$ captures the hydrodynamic response to an elliptically-shaped overlap region, quantified by the eccentricity vector

$\displaystyle\mathcal{E}_{2}=\varepsilon_{2}e^{2i\Phi}\equiv-\frac{\int e(r,\phi)r^{2}e^{2i\phi}rdrd\phi}{\int e(r,\phi)r^{2}rdrd\phi}\;,$

(1)

defined from transverse energy density $e(r,\phi)$ at a given $\eta$. The longitudinal structure of $V_{2}$ is driven by the asymmetrical energy deposition by forward and backward moving participants, resulting in a $\mathcal{E}_{2}$ that interpolates between $\mathcal{E}_{2,\mathrm{F}}$ and $\mathcal{E}_{2,\mathrm{B}}$, the eccentricity vectors defined by the forward-going and backward-going colliding nucleons, respectively Jia and Huo (2014). In other words, the initially produced QGP has an elliptic shape that twists and deforms from $\mathcal{E}_{2,\mathrm{F}}$ to $\mathcal{E}_{2,\mathrm{B}}$, driving a final-state $V_{2}$ with a similar $\eta$ dependence (top of Fig. 1).

Experimentally, elliptic flow is accessed using the two-particle correlation (2PC) method, which quantifies the second moment of its event-by-event distribution,

$\displaystyle V_{2\Delta}(\eta_{1},\eta_{2})=\left\langle V_{2}(\eta_{1})V_{2}^{*}(\eta_{2})\right\rangle\;.$

(2)

The initial state longitudinal fluctuations result in a decrease of $V_{2\Delta}$ for pairs with increasing separation $\Delta\eta=\eta_{1}-\eta_{2}$. Specifically, “global” decorrelations arising from nucleon fluctuations and nuclear deformation lead to gradual variations of $V_{2\Delta}$ with $\Delta\eta$ (the two sources have similar $\eta$ dependence Zhang et al. (2024)), whereas “local” decorrelations in energy deposition results in variation of $V_{2\Delta}$ at small $|\Delta\eta|$ (bottom of Fig. 1) Bozek and Broniowski (2016). However, the local decorrelations and non-flow are difficult to distinguish, as both appear at relatively small $|\Delta\eta|$. An experimental separation between the initial state short-range correlations and final-state non-flow remains elusive.

In this paper, using correlations between $V_{2}$ and different estimates of $\mathcal{E}_{2}$, we identify sources of initial-state-driven decorrelations contributing to $V_{2\Delta}(\eta_{1},\eta_{2})$. By comparing U+U with Au+Au collisions, we discover a fingerprint of short-range correlation as a dilution of the impact of nuclear deformation, leading to an increase of the ratio of $V_{2\Delta,\mathrm{UU}}/V_{2\Delta,\mathrm{AuAu}}$ with $\Delta\eta$. We further develop a method to disentangle between long-range correlation, short-range correlation, and final-state non-flow.

Setup. The study is carried out using the transport model, AMPT Lin et al. (2005), which simulates the full space-time evolution of heavy-ion collisions event by event. The nucleons are sampled from the Woods-Saxon distribution in polar coordinate: $\rho(r,\theta,\phi)\propto[1+\exp((r-R_{0}(1+\beta Y_{2,0}(\theta,\phi))/a_{0})]^{-1}$, where the $\beta$ is the quadrupole deformation, $R_{0}$ is the nuclear radius, and $a_{0}$ is the skin thickness. Collisions for three systems are simulated, whose geometrical parameters, listed in Table 1, are taken from Ref. Ryssens et al. (2023). The collision events are generated at $\mbox{$\sqrt{s_{\mathrm{NN}}}$}=193$ GeV for U+U and $\mbox{$\sqrt{s_{\mathrm{NN}}}$}=200$ GeV for Au+Au, employing the string-melting mode with a partonic cross-section of 3 $m$b. Elliptic eccentricities are calculated for all participating nucleons, denoted as $\mathcal{E}_{2}$, as well as for the partons at mid-rapidity $|\eta|<0.4$ after string melting but before the partonic transport, $\mathcal{E}_{2,\mathrm{quark}}$. The partons are weighted by their total energy in the calculation according to Eq. 1, hence $\mathcal{E}_{2,\mathrm{quark}}$ captures the initial state local fluctuation associated with strings (Fig. 1) without being impacted by non-flow. The final state particles used for the analysis are chosen from the transverse momentum range $0.2<p_{\mathrm{T}}<3$ GeV/$c$ and pseudorapidity range $|\eta|<5$. Event centrality is defined by charged particle multiplicity in $|\eta|<0.5$ ($N_{\mathrm{ch}}$). This analysis focuses on the 0–5% most central collisions, where the impact of nuclear deformation is prominent.

Elliptic flow is calculated using two different approaches. First, the flow is obtained using the standard 2PC method, $V_{2\Delta}(\eta_{1},\eta_{2})$, in the pseudorapidity range $|\eta_{1}|,|\eta_{2}|<2$. This is the typical range covered by the detectors at RHIC and the Large Hadron Collider (LHC). The $V_{2\Delta}(\Delta\eta)=\int V_{2\Delta}(\eta,\eta-\Delta\eta)d\eta/\int d\eta$, covering the range $|\Delta\eta|<4$, captures all the initial-state-driven long-range and short-range correlations illustrated by Fig. 1, as well as final-state non-flow.

In the second approach, we calculate the projection of the flow along the direction of eccentricity vectors,

$\displaystyle v_{2,\varepsilon}(\eta)\equiv\frac{\left\langle V_{2}(\eta)\mathcal{E}_{2}^{*}\right\rangle}{\sqrt{\left\langle\mathcal{E}_{2}\mathcal{E}_{2}^{*}\right\rangle}}\;,\;v_{2,\varepsilon_{\mathrm{quark}}}(\eta)\equiv\frac{\left\langle V_{2}(\eta)\mathcal{E}_{2,{\mathrm{quark}}}^{*}\right\rangle}{\sqrt{\left\langle\mathcal{E}_{2,{\mathrm{quark}}}\mathcal{E}_{2,{\mathrm{quark}}}^{*}\right\rangle}}\;$

(3)

where “$\left\langle\right\rangle$” indicates an average over events in a specific centrality range. The $v_{2,\varepsilon}(\eta)$ captures initial-state-driven long-range correlations, while $v_{2,\varepsilon_{\mathrm{quark}}}(\eta)$ captures in addition the initial-state-driven short-range correlations. Both observables are free of non-flow. Therefore, comparing the shapes of $v_{2,\varepsilon}^{2}(\eta)$ and $v_{2,\varepsilon_{\mathrm{quark}}}^{2}(\eta)$ to the shape of $V_{2\Delta}(\Delta\eta)$ facilitates the separation of the initial- and final-state effects.

The values of $V_{2\Delta}$, displayed in Fig. 2(a), decrease gradually with $\Delta\eta$ for all three systems. This decreasing trend is characterized by a mild structure at $|\Delta\eta|<2$ and a broader one at $|\Delta\eta|>2$, which as we show later indicates non-flow and initial-state-driven short-range correlations, respectively. The differences

$\displaystyle V_{2\Delta,\beta}(\Delta\eta)=V_{2\Delta}\{\beta=0.28\}-V_{2\Delta}\{\beta=0\}\approx b(\Delta\eta)\beta^{2}\;,$

(4)

capture the deformation-driven contribution of Uranium, and has a quadratic dependence on deformation parameter Jia (2022); Jia and Zhang (2023). Based on this discussion, the 2PC flow has four components,

$\displaystyle V_{2\Delta}(\Delta\eta)=a(\Delta\eta)+b(\Delta\eta)\beta^{2}+\mathrm{src}(\Delta\eta)+\mathrm{nf}(\Delta\eta)\;.$

(5)

The $a(\Delta\eta)$ and $\mathrm{src}(\Delta\eta)$ capture the initial-state-driven long-range and short-range correlations in collisions of spherical nuclei, respectively, $b(\Delta\eta)$ has similar long-range structure as $a(\Delta\eta)$ Zhang et al. (2024), and $\mathrm{nf}(\Delta\eta)$ is associated with non-flow.

Figure 2(b) shows the ratio of $V_{2\Delta}$ between U+U with and without deformation, as well as between U+U and Au+Au. The former can be expressed as

$\displaystyle R_{V_{2\Delta}}(\Delta\eta)=\frac{a(\Delta\eta)+b(\Delta\eta)\beta^{2}+\mathrm{src}(\Delta\eta)+\mathrm{nf}(\Delta\eta)}{a(\Delta\eta)+\mathrm{src}(\Delta\eta)+\mathrm{nf}(\Delta\eta)}\;.$

(6)

Since $\mathrm{src}(\Delta\eta)$ and $\mathrm{nf}(\Delta\eta)$ are short-range, they dilute the impact of nuclear deformation at small $|\Delta\eta|$, leading to an increase of the ratio towards large $|\Delta\eta|$. The increase slows down at $|\Delta\eta|\gtrsim 1.5$, suggesting $\mathrm{src}(\Delta\eta)$ and $\mathrm{nf}(\Delta\eta)$ have different ranges in $\Delta\eta$.

Figure 2(c) and (e) display the results of $v_{2,\varepsilon}^{2}$ and $v_{2,\varepsilon_{\mathrm{quark}}}^{2}$ over the full $\eta$ range. The values are smaller than $V_{2\Delta}$ as expected since they are projections of flow signal and free of non-flow. Their shapes signify the expected $\eta$ dependence of the single particle flow.

The corresponding ratios between U+U with and without deformation are displayed in the bottom panels. The $R_{v_{2,\varepsilon}^{2}}$ is nearly independent $\eta$, which is expected since $v_{2,\varepsilon}$ only captures the long-range components $a$ and $b\beta^{2}$ in Eq. 6, both of which have a similar shape Zhang et al. (2024). On the other hand, $R_{v_{2,\varepsilon_{\mathrm{quark}}}^{2}}$ displays a gradual increase with $\eta$, with a shape that closely resembles Fig. 2(a). Our interpretation is that the local structures of the initial state at mid-rapidity, captured by $\mathcal{E}_{2,{\mathrm{quark}}}$, generate additional flow around mid-rapidity. This is supported by Fig. 2(e), which shows that $v_{2,\varepsilon_{\mathrm{quark}}}>v_{2,\varepsilon}$ at $|\Delta\eta|<2$ but they are much closer at $|\Delta\eta|>2$. The presence of initial state short-range correlation in $\varepsilon_{\mathrm{quark}}$, contributing to $V_{2}$, is responsible for the gradual increase of $R_{v_{2,\varepsilon_{\mathrm{quark}}}^{2}}$ in Fig. 2(f). The bottom panels of Fig. 2 also show ratios between U+U and Au+Au by red circles, which are smaller than ratios between U+U with and without deformation and exhibit weaker $\Delta\eta$ oder $\eta$ dependence as well. This suggests that the impacts of Uranium deformation are partially canceled by the moderate deformation of ¹⁹⁷Au.

Clearly, the $\Delta\eta$ dependence of the deformation-driven flow is a reflection of the short-range correlation and its rapidity structure. To quantify the sources of decorrelations in $V_{2\Delta}$, we perform a convolution of $v_{2,\varepsilon}(\eta)$ and $v_{2,\varepsilon_{\mathrm{quark}}}(\eta)$ within $|\eta|<2$:

$\displaystyle V_{2\Delta,\varepsilon}(\Delta\eta)=\frac{1}{4}\int_{-2}^{2}v_{2,\varepsilon}(\eta_{1})v_{2,\varepsilon}(\eta_{2})\delta(\eta_{1}-\eta_{2})d\eta_{1}d\eta_{2}$

(7)

The expression for $V_{2\Delta,\varepsilon_{\mathrm{quark}}}(\Delta\eta)$ is analogous. These quantities approximate the initial-state-driven flow decorrelations and can be directly compared to $V_{2\Delta}$.

This comparison is carried out in Fig. 3. The $V_{2\Delta,\varepsilon_{\mathrm{quark}}}$ has a shape that matches almost perfectly with $V_{2\Delta}$ at $|\Delta\eta|>2$, while the $V_{2\Delta,\varepsilon}$ is much flatter in $|\Delta\eta|$. The shape of $V_{2\Delta,\varepsilon}$ is also similar to the shape of $V_{2\Delta,\beta}$ in Eq. 4, suggesting a similarly weak $\Delta\eta$ dependence Zhang et al. (2024).

These insights enable us to separate components of longitudinal flow decorrelations (shaded areas). The long-range component, $V_{2\Delta,\varepsilon}$, dominates, followed by the short-range component induced by local fluctuations. This short-range component starts to decrease at $|\Delta\eta|>2$, reflecting the typical coherence length of the local fluctuations. The estimated non-flow component, between $V_{2\Delta}$ and $V_{2\Delta,\varepsilon_{\mathrm{quark}}}$, spans a range of $|\Delta\eta|<2$. Since $V_{2\Delta,\varepsilon_{\mathrm{quark}}}$ may not capture all the initial-state-driven short-range components, this method likely gives an upper-limit of non-flow. In the bottom panels, we compare this component to an estimate obtained by assuming non-flow is unmodified and scales inversely with $N_{\mathrm{ch}}$ from Ref. STAR Collaboration (2024)

$\displaystyle\mathrm{nf}_{s}(\Delta\eta)=V_{2\Delta,90-100\%}N_{\mathrm{ch},90-100\%}/N_{\mathrm{ch},0-5\%}\;.$

(8)

Note that 90–100% centrality events are selected based on number of participanting nucleons $N_{\mathrm{part}}$, instead of $N_{\mathrm{ch}}$ to avoid selection bias. When integrated over $|\Delta\eta|<2$, the non-flow from the scaling method is about 30–40% lower. The long tail in $\mathrm{nf}_{s}(\Delta\eta)$, arising mainly from the away-side jet fragmentation, is completely absent in the subtraction method. This is due to jet quenching in AMPT, which is very effective in erasing $\cos(2\Delta\phi)$ component of the away-side at low $p_{\mathrm{T}}$ ¹¹1Away-side non-flow shape at low $p_{\mathrm{T}}$ is very broad and dominated by $\cos(\Delta\phi)$ driven by momentum conservation, which is the only harmonics surviving after medium effects..

Equipped with this knowledge, we are ready to dissect the rising trend of the $R_{V_{2\Delta}}$ in Fig. 2(b). We calculate three ratios between U+U in Fig. 3(a) and Au+Au in Fig. 3(c), $R_{V_{2\Delta}}$, $R_{V_{2\Delta,\varepsilon_{\mathrm{quark}}}}$, and $R_{V_{2\Delta,\varepsilon}}$, and display the results in Fig. 4. According to Eq. 6, the ordering, $R_{V_{2\Delta}}\lesssim R_{V_{2\Delta,\varepsilon_{\mathrm{quark}}}}<R_{V_{2\Delta,\varepsilon}}$, signifies the dilution of the impact of Uranium deformation by non-flow and initial-state short-range correlations. About half of the increase of $R_{V_{2\Delta}}$ is due to non-flow concentrated at $|\Delta\eta|<2$, while the remaining increase over $2<|\Delta\eta|<4$ is due to short-range correlations captured by $R_{V_{2\Delta,\varepsilon_{\mathrm{quark}}}}$. These initial-state-driven short-range correlations have a fairly wide range not exhausted even at $\Delta\eta=4$ where $R_{V_{2\Delta}}$ and $R_{V_{2\Delta,\varepsilon_{\mathrm{quark}}}}$ is still below $R_{V_{2\Delta,\varepsilon}}$. We stress that the impact of short-range correlations and non-flow was revealed only because of the large $\beta_{\mathrm{U}}$. Therefore, ratios of isobar or isobar-like systems with different nuclear shapes are excellent tools to explore the rapidity structure of the QGP’s initial state.

Recently, STAR has measured the $R_{V_{2\Delta}}$ between U+U and Au+Au collisions in the range $|\Delta\eta|<2$, and compared with a hydrodynamic model to constrain the deformation of the Uranium nucleus STAR Collaboration (2024). The hydrodynamic model includes non-flow from resonance decays but not jets, and has no long-range decorrelation effects. The STAR measurement was performed within the two-subevent method, nearly equivalent to the 2PC method with a cut of $0.8<|\Delta\eta|<2$. For this choice, the estimated non-flow in AMPT model is 6.8% for Au+Au collisions (see Table 2), reducing to 4% for $R_{V_{2\Delta}}$ between U+U and Au+Au, consistent with the uncertainty quoted by the STAR measurement STAR Collaboration (2024). Previous estimations rely on multi-parameter fit of $R_{V_{2\Delta}}(\Delta\eta)$ assume that non-flow either has a double-gaussian shape Abdulhamid et al. (2024) or having no flow decorrealtions Abdelwahab et al. (2015). Due to the unconstrained nature of these assumptions, the extracted non-flow fraction may vary up to 20% in central Au+Au collisions. Our approach serves as a valuable tool to expose potential limitations of any non-flow estimation methods.

The presence of a significant initial-state-driven short-range correlation impacts the interpretation of the flow decorrelations based on the factorization ratio $r_{2}(\eta)=\left\langle V_{2}(-\eta)V_{2}^{*}(\eta_{\mathrm{ref}}\right\rangle/\left\langle V_{2}(\eta)V_{2}^{*}(\eta_{\mathrm{ref}})\right\rangle$ Khachatryan et al. (2015); Aaboud et al. (2018), where $\eta_{\mathrm{ref}}$ is typically chosen from forward region, e.g. $\eta_{\mathrm{ref}}>4$ with a gap of at least 2 units from $\eta$, e.g. $|\eta|<2$. However, the initial-state-driven short-range correlations may span a $|\Delta\eta|$ range larger than $|\eta_{\mathrm{ref}}-\eta|$, reducing the values of $r_{2}$. This component is local and stochastic in nature and therefore is expected to cause random fluctuations of amplitude and direction of $V_{2}$ Pang et al. (2016); Xu et al. (2022). In this sense, $r_{2}(\Delta\eta)$ can not be interpreted as reflecting just the decorrelations between $-\eta$ and $\eta$.

Summary. We show that initial-state-driven longitudinal flow decorrelations have several components with different ranges in $\Delta\eta$. In addition to the usual global long-range component that are nearly independent of $\Delta\eta$, the flow decorrelations has a significant short-range component that is broader than non-flow correlations. This short-range component was shown to dilute the impact of nuclear deformation out to $\Delta\eta=4$. Our approach provides a new tool to analyze the initial longitudinal structure and final state response in transport or hydrodynamic models and can be extended to other bulk observables such as triangular flow and $p_{\mathrm{T}}$ fluctuations.

Understanding the longitudinal dynamics of the QGP is the frontier of high-energy nuclear research. The ability to vary the nuclear shape in isobar or isobar-like collisions provides a means to expose these dynamics. Future research should explore such a possibility by leveraging the large acceptance detectors such as ALICE Phase2 upgrade Acharya et al. (2024) and fixed target program of the SMOG2 at LHCb LHCb Collaboration .

This work is supported by DOE Research Grant Number DE-SC0024602.