# thermalize.m

Modifies a symmetric relaxation superoperator to relax the system towards a user-specified state using the Levante-Ernst method ( https://doi.org/10.1016/0009-2614(95)00578-R ) which works by adding one-way cross-relaxation with the unit state.

## Syntax

    R=thermalize(spin_system,R,rho)


## Arguments

  R   - a relaxation superoperator that drives the
system to the zero state

rho - desired relaxation destination


## Outputs

  R   - the relaxation superoperator that drives the
system to the state requested


## Notes

Liouville space spherical tensor formalism only.