“…In open quantum systems (OQSs), the time-evolution of a subsystem interacting with either a single or multiple environments, is modeled by a quantum master equation (QME) where ρ is the reduced density matrix, and F is the Liouvillian [1]. The simulation of open quantum many-body systems has lead to the development of various algorithms, e.g., renormalization group [2][3][4], meanfield methods [5][6][7], tensor networks [8][9][10][11][12], hierarchical equations of motion [13][14][15], Heisenberg equation of motion approaches [16,17], secular and non-secular Redfield theory [18,19], tensor transfer methods [20][21][22], and mixed quantum-classical methods [23][24][25].…”