Abstract:The evolution of a relativistic heavy-ion collision is typically understood as a process that transmutes the initial geometry of the system into the final momentum distribution of observed hadrons, which can be described via a cumulant expansion of the initial distribution of energy density and is represented at leading order as the well-known eccentricity scaling of anisotropic flow. We extend this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. We confirm the validity of the framework in state-of-the-art hydrodynamic simulations of large and small systems. With this new framework, it is possible to separate the effects of early-time dynamics from those of final-state evolution, even in the case when the distribution of energy does not fully determine subsequent evolution, as for example, in small systems. Specifically, we answer the question of when and how azimuthal correlations from the initial state survive to the final state. In very small systems such as $p$-$p$, for example, initial momentum degrees of freedom dominate over energy. Thus, even if the system forms a quark-gluon plasma that is well described by hydrodynamics, the usual hydrodynamic picture of the transmutation of initial geometry to final momentum anisotropy is broken. Nevertheless, we show that the hydrodynamic response to the full energy-momentum tensor can be well understood in a similar manner as larger systems. Additionally, this framework elucidates the generic features of the system's evolution that are responsible for the impressive success of hydrodynamic simulations, but which may still hold even in cases where hydrodynamics is not applicable.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: How to more comprehensively understand the conversion process from the initial - state energy - momentum distribution to the final - particle momentum distribution in the process of heavy - ion collisions. Specifically, the existing frameworks mainly focus on the influence of the initial geometry on the final - particle momentum distribution, but ignore the role of other initial conditions (such as momentum - space properties). This paper proposes a new framework, which extends the traditional description method to include other components of the energy - momentum tensor \(T_{\mu\nu}\) in the initial state, so as to more accurately predict the final anisotropic flow.
### Specific problems
1. **Limitations of the existing framework**:
- The existing frameworks mainly predict the final - particle momentum distribution through the initial energy - density distribution, especially by describing it through the eccentricity - scaling relation of the anisotropic flow \(V_2=\kappa_2\epsilon_2\).
- This method performs well in large systems, but in small systems (such as proton - proton collisions), the initial - momentum degrees of freedom may be dominant, causing the traditional framework to fail.
2. **Requirement for an extended framework**:
- A more general framework is needed, which can consider both the energy density and other momentum - space properties (such as momentum density and stress tensor) in the initial state simultaneously.
- This extension helps to better understand the separation between the early - dynamic effects and the final - state evolution, especially in small systems where the energy distribution is not sufficient to fully determine the subsequent evolution.
### Solutions
1. **New generating function**:
- A new generating function \(\rho(\vec{x}_\perp)\) is proposed, which includes not only the energy density \(T_{\tau\tau}\), but also the contributions of the momentum density \(T_{\tau i}\) and the stress tensor \(T_{ij}\):
\[
\rho(\vec{x}_\perp)=T_{\tau\tau}(\vec{x}_\perp)-\alpha\partial_i T_{\tau i}(\vec{x}_\perp)+\beta\partial_i\partial_j T_{ij}(\vec{x}_\perp)
\]
- Here, \(\alpha\) and \(\beta\) are response coefficients, representing the importance of the momentum density and the initial transverse stress relative to the energy density respectively.
2. **Cumulant expansion**:
- Use the new generating function to perform cumulant expansion and construct the expressions of estimators \(V_n\), which can be systematically improved to improve the prediction accuracy.
- For example, for the elliptic flow \(V_2\) and the triangular flow \(V_3\), new lowest - order estimators are proposed:
\[
V_2^{(\text{est})}=-\kappa_2\frac{\langle r^2 e^{i2\phi}\rangle_E - 2\alpha\langle r e^{i\phi}\rangle_U - 2\beta\langle 1\rangle_C}{\langle|r e^{i\phi}-\langle r e^{i\phi}\rangle_E|^2\rangle_E}
\]
\[
V_3^{(\text{est})}=-\kappa_3\frac{\langle r^3 e^{i3\phi}\rangle_E - 3\alpha\langle r^2 e^{i2\phi}\rangle_U - 6\beta\langle r e^{i\phi}\rangle_C}{\langle|r e^{i\phi}-\langle r e^{i\phi}\rangle_E|^3\rangle_E}
\]