“…[62,63,77,78], the HEOM method has been extended to arbitrary spectral density function as well as finite temperature environment situation. Moreover, it has been reported that the HEOM method can be extended to simulate the dissipative dynamics of a few-level system embedded in a fermionic environment [53,65,79] or a spin environment [80,81]. It would be very interesting to extrapolate our study to these more general situations.…”