We can integrate this equation on condition that the interaction picture of the reduced density matrix ρ ρ I ( t ) changes slowly in time (which is true if R is small), then ρ ρ a b , I ( t ) ≈ ≈ ρ ρ a b , I ( 0 ) , getting