One fluid to rule them all: viscous hydrodynamic description of event-by-event central p+p, p+Pb and Pb+Pb collisions at sqrt{s}=5.02 TeV

The matter created in central p+p, p+Pb and Pb+Pb collisions at sqrt{s}=5.02 TeV is simulated event-by-event using the superSONIC model, which combines pre-equilibrium flow, viscous hydrodynamic evolution and late-stage hadronic rescatterings. Employing a generalization of the Monte Carlo Glauber model where each nucleon possesses three constituent quarks, superSONIC describes the experimentally measured elliptic and triangular flow at central rapidity in all systems using a single choice for the fluid parameters, such as shear and bulk viscosities. This suggests a common hydrodynamic origin of the experimentally observed flow patterns in all high energy nuclear collisions, including p+p.