Curing the numerical Cherenkov instability in 3+1D Glasma simulations

4 Nov 2019, 17:40
Wanda Han Show Theatre & Wanda Reign Wuhan Hotel

Wanda Han Show Theatre & Wanda Reign Wuhan Hotel

Poster Presentation Initial state and approach to equilibrium Poster Session


Andreas Ipp (TU Wien)


We present our results on genuinely 3+1D simulations of the Glasma in heavy ion collisions. By incorporating nuclei with finite thickness along the beam axis, we naturally break boost invariance and obtain Gaussian-like rapidity profiles of energy densities already at the tree level [1]. The profiles resemble strong coupling results and agree surprisingly well with experimental data of pion multiplicities as obtained at RHIC. The applicability of these simulations is limited by the numerical Cherenkov instability. We construct a new lattice action which allows us to cure this instability along one lattice direction [2]. This paves the way for simulations at higher energies and better resolution.

[1] AI, D. Müller, Phys.Lett. B771 (2017) 74-79, arXiv:1703.00017
[2] AI, D. Müller, Eur.Phys.J. C78 (2018) no.11, 884, arXiv:1804.01995

Primary author

Andreas Ipp (TU Wien)


David Mueller (Vienna University of Technology)

Presentation Materials