  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.
