Many interesting fundamental and practical physical problems arise in open quantum systems [1]}, [2]} interacting with an environment acting as a thermal bath. One of the most widely used equations to study open quantum systems is the quantum master equation (QME), which describes the temporal evolution of the reduced density operator of the system. The reduced density operator is obtained by tracing over the variables describing the environment, for which an explicit description of dynamics is not usually feasible and often not of interest.
