  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.
