Quarkonium production in ultrarelativistic heavy ions collisions is one of the best probes of the QGP formed in these collisions. Resorting to accurate methods to describe the $Q\bar{Q}$
evolution in a QGP is a prerequisite for the precise interpretation of experimental data. Among these methods, the quantum master equations (QME) derived within the formalism of open quantum systems are particularly relevant. We present exact numerical solutions in a 1D setting of previously derived quantum master equations (QME) in their quantum Brownian regime. Distinctive features of the in-medium bottomonia
evolution with the QME are presented; some phenomenological consequences are addressed by considering evolutions for a fixed as well as EPOS4 temperature profiles. Next, we investigate the accuracy of the semiclassical approximation (often used to describe charmonium production in URHIC) by benchmarking the corresponding evolutions on the exact solutions derived with the QME for the case of a $c\bar{c}$
pair.

