Bottomonium suppression in 5.02 and 8.16 TeV $p$-Pb collisions

We compute the suppression of $\Upsilon \left(1S\right)$, $\Upsilon \left(2S\right)$, and $\Upsilon \left(3S\right)$ states in $p$-Pb collisions relative to $pp$ collisions, including nuclear parton distribution function (nPDF) effects, coherent energy loss, momentum broadening, and final-state interactions in the quark-gluon plasma. We employ the EPPS21 nPDFs and calculate the uncertainty resulting from variation over the associated error sets. To compute coherent energy loss and momentum broadening, we follow the approach of Arleo, Peigne, and collaborators. The $3+1D$ viscous hydrodynamical background evolution of the quark-gluon plasma is generated by anisotropic hydrodynamics. The in-medium suppression of bottomonium in the quark-gluon plasma is computed using a next-to-leading-order open quantum system framework formulated within potential nonrelativistic quantum chromodynamics. We find that inclusion of all these effects provides a reasonable description of experimental data from the ALICE, ATLAS, CMS, and LHCb Collaborations for the suppression of $\Upsilon \left(1S\right)$, $\Upsilon \left(2S\right)$, and $\Upsilon \left(3S\right)$ as a function of both transverse momentum and rapidity.