Semianalytical calculation of the trajectory of relativistic nuclear collisions in the QCD phase diagram

Abstract

We extend a semianalytical model that includes the finite nuclear thickness to calculate the energy density $\ensuremath{\epsilon}(t)$ and conserved-charge densities including the net-baryon density ${n}_{\mathrm{B}}(t)$ produced at mid-spacetime-rapidity in central Au + Au collisions. Assuming the formation of a quark-gluon plasma with an ideal gas equation of state of either quantum or Boltzmann statistics or with a lattice QCD-based equation of state, we extract the temperature $T$ and chemical potentials ${\ensuremath{\mu}}_{\mathrm{B}}, {\ensuremath{\mu}}_{\mathrm{Q}}$, and ${\ensuremath{\mu}}_{\mathrm{S}}$ as functions of time. This then allows us to semianalytically calculate the $T\text{\ensuremath{-}}{\ensuremath{\mu}}_{\mathrm{B}}$ trajectory of relativistic nuclear collisions in the QCD phase diagram, which should benefit the studies of high density physics including the search for the critical endpoint. This model is also useful for exploring the trajectories in the more general $T\text{\ensuremath{-}}{\ensuremath{\mu}}_{\mathrm{B}}\text{\ensuremath{-}}{\ensuremath{\mu}}_{\mathrm{Q}}\text{\ensuremath{-}}{\ensuremath{\mu}}_{\mathrm{S}}$ QCD phase space.