Assuming that the population remains small, one can truncate at given values of n1 and n2 to yield a finite set of differential equations. For an initial condition, we take the population to be zero, that is, ρ1,1,1,1 = 1 with all other elements of the density matrix set to zero. Eq. A.6 can then be solved numerically (we used the Adams-Moulton-Bashforth procedure) until a steady state is reach...