“…Specifically, the dynamic equation for the perturbation density matrix can be obtained as ∂ρ/∂t = [H, ρ]/ih−(dρ/dt) relax , where H is the Hamiltonian and the term (dρ/dt) relax describes the decay. Taking the rotating-wave, weak-field and steady-state approximations, we can find a general solution of the dynamic equation in which the third-and fifth-order density matrix elements are [15,16] …”