Semi-classical treatment of bottomonium suppression in p-Pb collisions

We present a detailed study of bottomonium (Υ(1S), Y(2S), and Y(3S)) suppression in minimum-bias proton-Lead (p+Pb) collisions at 5.02 and 8.16 TeV. This study incorporates both cold nuclear matter (CNM) effects (nuclear parton distribution function(nPDF) effect, coherent energy loss and momentum broadening), and hot nuclear matter (HNM) effect due to the hot Quark-Gluon-Plasma (QGP) medium. The nPDF effect is computed using EPPS21, while coherent energy loss and momentum broadening follow the approach developed by Arleo and Peigne. The QGP evolution is simulated using 3+1D anisotropic hydrodynamics (aHydro). To calculate in-medium bottomonium suppression, we employ a semi-classical kinetic rate equation using perturbative and non-perturbative bottomonium reaction rates [1], and compared our findings to those obtained from a next-to-leading order open quantum system framework formulated within potential nonrelativistic quantum chromodynamics (pNRQCD) [2]. Finally, combining the CNM and HNM effects, we compute the nuclear modification factor (R^Y_pA) of the bottomonium states as a function of its rapidity (y) and transverse momentum (p_T), and compare our predictions with experimental data from the ALICE, ATLAS, CMS, and LHCb Collaborations. Incorporating all these effects gives a good description of the experimental data, supporting the idea of hot and short-lived QGP formation in the min-bias p-Pb collisions at the LHC energies.