Running the gamut of high energy nuclear collisions

Abstract: We present calculations of bulk properties and multiparticle correlations in a large variety of collision systems within a hybrid formalism consisting of IP-Glasma initial conditions, MUSIC viscous relativistic hydrodynamics, and UrQMD microscopic hadronic transport. In particular, we study heavy ion collisions at the Large Hadron Collider (LHC), including Pb+Pb, Xe+Xe, and O+O collisions, and Au+Au, U+U, Ru+Ru, Zr+Zr, and O+O collisions at the Relativistic Heavy Ion Collider (RHIC). We further study asymmetric systems, including p+Au, d+Au, $^3$He+Au, and p+Pb collisions at various energies as well as p+p collisions at 0.5 and 13 TeV. We describe experimental observables in all heavy ion systems well with one fixed set of parameters, validating the energy and system dependence of the framework. Many observables in the smaller systems are also well described, although they test the limits of the model. Calculations of O+O collisions provide predictions for potential future runs at RHIC and LHC.