“…Since the spectral poses challenge to many existing method, it arouse much attention recently. Till now, it has been studied by many different methods including the quasi-adiabatic propagator path integral (QUAPI) [22,27,28], the Van Vleck perturbation theory together with a Born-Markov master equation [29], the flow equation renormalization [30,31,32], the non-interacting blip approximation (NIBA) [31,33], and generalized polaron transformation method [34,35,36]. Again, most of the works have not considered the effect of temperature.…”