Eigenfields.m

Computes resonance fields. For a Hamiltonian Hc+b*Hz, returns all magnetic fields b for which the difference between two eigenvalues of Hc+b*Hz is equal to the frequency provided, and the transition moment across the specified operator is significant.

Syntax

    [tf,tm]=eigenfields(spin_system,parameters,Hc,Hz,Hmw)

Arguments

    Hc     -  laboratory frame Hamiltonian operator (Hilbert 
              space) or commutation superoperator (Liouville
              space, containing all spin-spin couplings, but
              no Zeeman terms

    Hz     -  laboratory frame Hamiltonian operator (Hilbert 
              space) or commutation superoperator (Liouville
              space, containing only Zeeman terms at 1 Tesla

    Hmw    -  microwave irradiation operator, without the am-
              plitude prefactor

    parameters.window   -  magnet field window, Tesla

    parameters.mw_freq  -  microwave frequency, Hz

    parameters.tm_tol   -  relative transition moment 
                           tolerance

    parameters.pp_tol   -  peak position tolerance in Tesla,
                           this should be much smaller than
                           the typical line width

Outputs

    tf     -  vector of transition fields in Tesla

    tm     -  vector of transition moments

Notes

In Hilbert space, the very efficient Schweiger-Stoll method is used. In Liouville space, the very general but rather slow generalised eigensolver supplied with Matlab is used.

See also

fieldsweep.m

Version 2.6, authors: Ilya Kuprov