steady.m

Steady state under the repeated action by the same dissipative evolution propagator.

Syntax

   rho=steady(spin_system,P,rho,method)

Parameters

   spin_system - Spinach spin system object

   P           - propagator, an exponential of the Liouvillian that
                 contains a thermalised relaxation superoperator, or
                 a product thereof, for example, from a repeating block
                 of a pulse or a pulse sequence

   rho         - optional initial guess for the steady state; a good
                 guess can significantly accelerate this function. Leave
                 empty otherwise; the state must have unit trace, which in
                 sphten-liouv means a first element equal to 1

   method      - steady state solver method: 'newton' by default, or
                 'squaring' for propagator squaring

Outputs

   rho         - steady state under the repeated application of the propagator P

Notes

Available for sphten-liouv and zeeman-liouv formalisms; the Newton-Raphson solver pins the first state vector element in sphten-liouv and the density matrix trace in zeeman-liouv.

See also

evolution.m, isergen.m, iserstep.m, krylov.m, step.m, Kernel_functions

Version 2.11, authors: Ilya Kuprov