where ρ ρ is the density operator describing the state of the working medium and L h/c is the generator of dissipative dynamics which includes the heat transport terms from the baths. Using this construction, the total change in energy of the sub-system becomes: