Numerical generalised cumulant expansion route to the relaxation superoperator to second order, which is equivalent to Redfield theory.
H - a cell array of spin Hamiltonian commutation superoperators obtained from an MD trajectory
dt - time step of the MD trajectory, seconds
R - laboratory frame relaxation superoperator
The trajectory must be at least five times longer than the longest correlation time expected in the system.
Version 2.2, authors: Ilya Kuprov