We compute the quarkonium nuclear modification factor in a strongly coupled quark-gluon plasma for quarkonia that are $S$-wave Coulombic bound states.
We perform the analysis in a non-relativistic effective field theory framework that is accurate at leading-order in the heavy-quark density expansion and at next-to-leading order in the multipole expansion. We write and solve the Lindblad equation for the heavy quark-antiquark density.
Thermal mass shift, width and the Lindblad equation depend on only two non-perturbative parameters: the heavy-quark momentum diffusion coefficient and its dispersive counterpart. Finally, we provide numerical results for the nuclear modification factors of the $1S$ and $2S$ bottomonium states.