Difference between revisions of "Eigenfields.m"

From Spinach Documentation Wiki
Jump to: navigation, search
(Created page with "{{DISPLAYTITLE:eigenfields.m}} __NOTOC__ Computes resonance fields. For a Hamiltonian Hc+b*Hz, returns all magnetic fields b for which the difference between any two eigenvalu...")
 
Line 1: Line 1:
 
{{DISPLAYTITLE:eigenfields.m}} __NOTOC__
 
{{DISPLAYTITLE:eigenfields.m}} __NOTOC__
−
Computes resonance fields. For a Hamiltonian Hc+b*Hz, returns all magnetic fields b for which the difference between any two eigenvalues of Hc+b*Hz is equal to omega.
+
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==
 
==Syntax==
  
−
     [b,uv]=eigenfields(Hc,Hz,omega,window)
+
     [tf,tm]=eigenfields(spin_system,parameters,Hc,Hz,Hmw)
  
 
==Arguments==
 
==Arguments==
  
−
     Hc    -  laboratory frame Hamiltonian commutation  
+
     Hc    -  laboratory frame Hamiltonian operator (Hilbert
−
               superoperator containing all spin-spin
+
              space) or commutation superoperator (Liouville
−
              couplings, but no Zeeman interactions
+
               space, containing all spin-spin couplings, but
 +
              no Zeeman terms
 
   
 
   
−
     Hz    -  laboratory frame Hamiltonian commutation  
+
     Hz    -  laboratory frame Hamiltonian operator (Hilbert
−
               superoperator containing Zeeman interac-
+
              space) or commutation superoperator (Liouville
−
              tions at 1 Tesla
+
               space, containing only Zeeman terms at 1 Tesla
 
   
 
   
−
     omega - resonance frequency in rad/s
+
     Hmw    - microwave irradiation operator, without the am-
 +
              plitude prefactor
 
   
 
   
−
     window -  frequency window to consider, a vector
+
     parameters.window   -  magnet field window, Tesla
−
              with two elements, rad/s
+
 +
    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==
 
==Outputs==
  
−
     b      -  vector of magnetic fields in Tesla
+
     tf    -  vector of transition fields in Tesla
−
+
 
−
     uv     -  a matrix of |u><v| dyadics where |u>
+
     tm     -  vector of transition moments
−
              and |v> are separated by omega under
 
−
              Hc+b*Hz
 
  
 
==Notes==
 
==Notes==
−
A very expensive and barely stable generalised eigensolver supplied with Matlab is used.
+
In Liouville space, a very expensive and barely stable generalised eigensolver supplied with Matlab is used.
  
 
==See also==
 
==See also==
−
[[eigentracks.m]], [[fieldsweep.m]]
+
[[fieldsweep.m]]
  
  
−
''Version 2.3, authors: [[Ilya Kuprov]]''
+
''Version 2.6, authors: [[Ilya Kuprov]]''

Revision as of 14:04, 1 July 2021

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 Liouville space, a very expensive and barely stable generalised eigensolver supplied with Matlab is used.

See also

fieldsweep.m


Version 2.6, authors: Ilya Kuprov