  In this proceedings contribution I review recent work which computes the
suppression of bottomonium production in heavy-ion collisions using open
quantum systems methods applied within the potential non-relativistic quantum
chromodynamics (pNRQCD) effective field theory. I discuss how the computation
of bottomonium suppression can be reduced to solving a
Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) quantum master equation for the
evolution of the b-bbar reduced density matrix. The open quantum systems
approach used allows one to take into account the non-equilibrium dynamics and
decoherence of bottomonium in the quark-gluon plasma. Finally, I present
comparisons of phenomenological predictions obtained using a recently obtained
next-to-leading-order GKSL equation with ALICE, ATLAS, and CMS experimental
data for bottomonium suppression and elliptic flow.
