3+1D hydrodynamic simulation of relativistic heavy-ion collisions

Abstract: We present MUSIC, an implementation of the Kurganov-Tadmor algorithm for relativistic 3+1 dimensional fluid dynamics in heavy-ion collision scenarios. This Riemann-solver-free, second-order, high-resolution scheme is characterized by a very small numerical viscosity and its ability to treat shocks and discontinuities very well. We also incorporate a sophisticated algorithm for the determination of the freeze-out surface using a three dimensional triangulation of the hyper-surface. Implementing a recent lattice based equation of state, we compute p_T-spectra and pseudorapidity distributions for Au+Au collisions at root s = 200 GeV and present results for the anisotropic flow coefficients v_2 and v_4 as a function of both p_T and pseudorapidity. We were able to determine v_4 with high numerical precision, finding that it does not strongly depend on the choice of initial condition or equation of state.

• The paper introduces a 3+1D hydrodynamic simulation using the KT scheme to accurately model the evolution of QGP and hadron gas.
• It utilizes a lattice-based EOS and a triangulated freeze-out surface to precisely calculate observables like elliptic (v2) and higher-order (v4) flow coefficients.
• The simulation results reliably replicate experimental pT-spectra and pseudorapidity distributions, paving the way for future viscous and jet quenching studies.

Overview of "3+1D Hydrodynamic Simulation of Relativistic Heavy-Ion Collisions"

This paper presents a detailed implementation of a 3+1 dimensional hydrodynamic simulation model for relativistic heavy-ion collisions, specifically using the Kurganov-Tadmor (KT) scheme. The focus is on simulating the evolution of a Quark-Gluon Plasma (QGP) and hadron gas formed in these collisions, within an ideal hydrodynamic framework. The authors, Schenke, Jeon, and Gale, introduce a computational tool called music, which employs a high-resolution central difference scheme, allowing efficient calculation of flow properties with minimal numerical viscosity.

Key Contributions and Methods

The primary contribution of the paper is the implementation of the KT algorithm for relativistic fluid dynamics within the specific context of heavy-ion collisions. The KT method is notable for its ability to effectively handle shocks and discontinuities, crucial for accurate simulations of high-energy collision dynamics. This approach circumvents the usual Riemann solver techniques, which simplifies calculations and enhances computational efficiency.

The authors incorporate a lattice-based equation of state (EOS) model that matches the transition temperature—derived from QCD lattice simulations—allowing the exploration of thermal features of the evolving QGP and hadron gas. This is critical for producing realistic simulations that match experimental results.

Another significant technical advancement is the triangulation method for determining freeze-out surfaces. This involves discretizing the freeze-out hyper-surface into tetrahedra, which are used to calculate particle spectra and flow coefficients with high precision. This development is particularly relevant for accurately simulating experimental observables such as the anisotropic flow coefficients, $v_2$ (elliptic flow) and $v_4$, across varying transverse momentum $p_T$ and pseudorapidity $\eta$.

Results and Implications

The paper reports calculated $p_T$-spectra and pseudorapidity distributions for Au+Au collisions at $\sqrt{s} = 200 \, \text{GeV}$. The simulations effectively match experimental data from STAR and PHENIX collaborations, validating the efficacy of the KT method and the EOS models in reproducing key collision characteristics.

A critical result is the detailed prediction of the $v_2$ and $v_4$ flow coefficients. The authors find that $v_4$ does not show strong sensitivity to the initial conditions or to variations in the EOS, offering new insights into the flow dynamics. This conclusion contradicts some earlier expectations in the field and suggests that other factors, such as viscosity, might play a more significant role in determining higher harmonic flows.

Future Directions

The authors plan to extend the current simulations to include viscous hydrodynamic effects, which remain beyond the current ideal model and will likely refine the predictive capability of the simulation, especially concerning higher momentum and rapidity regions of the collision. Moreover, coupling the model with high $p_T$ jet physics simulations, as intended, could offer deeper insights into jet quenching phenomena and energy loss mechanisms in QGP.

This work sets a robust foundation for future research by demonstrating the viability of using sophisticated numerical schemes such as the KT method for complex hydrodynamic simulations in nuclear physics. It provides a potential pathway for integrating more complex physical effects and achieving a more comprehensive theoretical understanding of heavy-ion collisions.