Bottomonium suppression and elliptic flow using Heavy Quarkonium Quantum Dynamics

Abstract

Abstract We introduce a framework called Heavy Quarkonium Quantum Dynamics (HQQD) which can be used to compute the dynamical suppression of heavy quarkonia propagating in the quark-gluon plasma using real-time in-medium quantum evolution. Using HQQD we compute large sets of real-time solutions to the Schrödinger equation using a realistic in-medium complex-valued potential. We sample 2 million quarkonia wave packet trajectories and evolve them through the QGP using HQQD to obtain their survival probabilities. The computation is performed using three different HQQD model parameter sets in order to estimate our systematic uncertainty. After taking into account final state feed down we compare our results to existing experimental data for the suppression and elliptic flow of bottomonium states and find that HQQD predictions are good agreement with available data for RAA as a function of Npart and pT collected at $$ \sqrt{s_{\mathrm{NN}}} $$ s NN = 5.02 TeV. In the case of v2 for the various states, we find that the path-length dependence of ϒ(1s) suppression results in quite small v2 for ϒ(1s). Our prediction for the integrated elliptic flow for ϒ(1s) in the 10−90% centrality class, which now includes an estimate of the systematic error, is v2[ϒ(1s)] = 0.003 ± 0.0007 ±$$ {}_{0.0013}^{0.0006} $$ 0.0013 0.0006 . We also find that, due to their increased suppression, excited bottomonium states have a larger elliptic flow. Based on this observation we make predictions for v2[ϒ(2s)] and v2[ϒ(3s)] as a function of centrality and transverse momentum.