Adds pumping terms to the relxation superoperator to enable approximate simulation of CIDNP type effects.
R - relaxation superoperator, from relaxation.m rho - the state to be pumped, from state.m rate - pumping rate, Hz
R - modified relaxation superoperator
See cidnp_pumping_1.m and cidnp_pumping_2.m files in examples/spin_chemistry directory.
1. For the pumping to work correctly, the unit state population (first element) in the state vector that R will be acting on must be set to 1.
2. This function is only available in sphten-liouv formalism, and may be called repeatedly if multiple states are pumped.
Version 2.1, authors: Ilya Kuprov