Thermalize.m

Modifies the relaxation superoperator to drive the system to the user-specified target state using the inhomogeneous master equation formalism, or to the equilibrium state of the lab-frame Hamiltonian at the temperature provided by the user using the DiBari-Levitt formalism.

Syntax

    R=thermalize(spin_system,R,HLSPS,T,rho_eq,method)

Parameters

  R       - symmetric negative definite relaxation superoperator that
            drives the system towards the zero state vector; this may
            be obtained from relaxation.m if inter.equilibrium is 'zero'

  HLSPS   - lab-frame Hamiltonian left side product superoperator,
            available from hamiltonian.m; also call orientation.m if
            necessary. This is not required for IME formalism, pass []

  T       - absolute temperature, not required for IME formalism,
            pass []

  rho_eq  - thermal equilibrium state, not required for DiBari-Levitt
            formalism, pass []

  method  - 'dibari' for DiBari-Levitt thermalisation, 'IME' for the
            inhomogeneous master equation

Outputs

  R       - thermalised relaxation superoperator

Notes

IME requires the population of the unit state in the state vector to be exactly 1; Spinach cannot check or enforce this requirement. DiBari-Levitt thermalisation is computationally expensive, but tends to work better than IME, particularly in exotic regimes.

See also

magpump.m, relaxation.m, kinetics.m, lindbladian.m

Version 2.2, authors: Ilya Kuprov