Difference between revisions of "Thermalize.m"
m (Correct function page syntax from current Spinach source) |
m (Update function page text from current Spinach source) |
||
| Line 1: | Line 1: | ||
{{DISPLAYTITLE:thermalize.m}} | {{DISPLAYTITLE:thermalize.m}} | ||
| − | Modifies | + | |
| + | 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== | ==Syntax== | ||
| Line 9: | Line 10: | ||
==Arguments== | ==Arguments== | ||
| − | R | + | |
| − | + | 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== | ==Outputs== | ||
| − | R | + | |
| − | + | R - thermalised relaxation superoperator | |
==Notes== | ==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== | ==See also== | ||
Revision as of 15:19, 5 June 2026
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.
Contents
Syntax
R=thermalize(spin_system,R,HLSPS,T,rho_eq,method)
Arguments
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
