Imprints of fluctuating proton shapes on flow in proton-lead collisions at the LHC

Results for particle production in root-s = 5.02 TeV p+Pb collisions at the Large Hadron Collider within a combined classical Yang-Mills and relativistic viscous hydrodynamic calculation are presented. We emphasize the importance of sub-nucleon scale fluctuations in the proton projectile to describe the experimentally observed azimuthal harmonic coefficients, demonstrating their sensitivity to the proton shape. We stress that the proton shape and its fluctuations are not free parameters in our calculations. Instead, they have been constrained using experimental data from HERA on exclusive vector meson production. Including temperature dependent shear and bulk viscosities, as well as UrQMD for the low temperature regime, we present results for mean transverse momenta, harmonic flow coefficients for charged hadrons and identified particles, as well as Hanbury-Brown-Twiss radii.