Here, t 0 is some initial time, where the total state of the system and bath is assumed to be factorized, and we have introduced the bath correlation function C m n ( t ) = tr ( E m , I ( t ) E n ρ ρ env,eq ) in terms of the density operator of the environment in thermal equilibrium, ρ ρ env,eq .