“…Thus, many approximations are adopted such as averaging out environment degrees of freedom 18 and classical description of the system and/or environment. 19,20 In the past three decades, great progress has been made in the development of theoretical approaches for open quantum systems. These approaches include the perturbative Redfield equation, 21 the exact Nakajima-Zwanzig formalism 22,23 and its kernel-based expansions, [24][25][26] the quantum-classical and fully classical approaches, 20,[27][28][29][30] Green's function formalism, 31 the transfer tensor method 32 and its extension, 33 the multi-configuration time-dependent Hartree (MCTDH), 34,35 the pseudo-mode approach, 36,37 the reaction coordinate (RC) approach, 38,39 the quantum Monte Carlo (QMC), 40,41 the time-dependent numerical renormalization group (NRG), 42 the density matrix renormalization group (tDMRG), 43 the polaron transformation, 24,44 the time evolving density matrix using the orthogonal polynomial algorithm (TEDOPA), 45,46 the quasiadiabatic propagator path integral (QuAPI), 47,48 the numerical variational method (NVM), 49 the automated compression of environment (ACE) method, 50 the hierarchy equations of motion (HEOM), [51][52][53][54]…”