Heavy quarkonium dynamics at next-to-leading order in the binding energy over temperature*

We present a systematic method for solving the Lindblad equation for heavy-quarkonium dynamics in the quark-gluon plasma which accounts for corrections that are next-to-leading order (NLO) in the ratio of the binding energy of the state and the temperature.  The method used relies on mapping the three-dimensional Lindblad evolution to the solution of the one-dimensional Schrodinger evolution with stochastically sampled quantum jumps.  We demonstrate how to achieve this dimensional reduction by writing the NLO complex Hamiltonian and jump operators in the spherical basis in which the operators act only on the radial part of the wave function.  As result, one can implement the quantum trajectories method to solve the NLO Lindblad equation.  Using the resulting NLO framework we can more reliably extend the calculation of heavy-quarkonium suppression to lower temperatures than is possible with the LO formalism.