“…The spinboson model has attracted considerable attention in past decades because it provides a universal model for numerous physical and chemical processes. The reduced system dynamics of spin-boson model has been studied by various analytical and numerical methods, for example, the polaron [5,6] or generalized Silbey-Harris [7][8][9][10][11][12] transformation approach, the time-dependent numerical renormalization group method [13][14][15][16], the quasi adiabatic propagator path integral, [9,17,18] and hierarchical equations of motion (HEOM) [19][20][21][22][23][24][25][26]. Each method has its own regimes of validity depending on the systembath coupling strength, the bath temperature and the bath spectral density function.…”