Bottomonium suppression and elliptic flow in heavy-ion collisions

In this proceedings contribution I review recent progress concerning the suppression of bottomonium production in the quark-gluon plasma. Making use of open quantum system methods applied to potential non-relativistic quantum chromodynamics one can show that the dynamics of heavy-quarkonium bound states satisfying the scale hierarchy 1/a_0 >> pi T ~ m_D >> E obey a Lindblad equation whose solution provides the quantum evolution of the heavy-quarkonium reduced density matrix. To solve the resulting Lindblad equation we use a quantum trajectories algorithm which allows one to include all possible angular momentum states of the quark-antiquark probe in a scalable manner. We solve the Lindblad equation using a tuned 3+1D dissipative hydrodynamics code for the background temperature evolution. We then consider a large number of Monte-Carlo sampled bottomonium trajectories embedded in this background. This allows us to extract the centrality- and p_T-dependence of the nuclear suppression factor R_AA[Upsilon] and elliptic flow v_2[Upsilon]. We find good agreement between our model predictions and available sqrt(s_NN) = 5.02 TeV Pb-Pb collision experimental data from the ALICE, ATLAS, and CMS collaborations.