At the initial moment, we assume that the system and the environment are factorized and the environment is in its own thermal equilibrium at an inverse temperature β = 1/(k B T ) with the Boltzmann constant k B and the temperature T . The initial density operator for the composite system-plusenvironment is given by