Difference between revisions of "Eigenfields.m"
(Update function See also links and function index membership) |
(sync with Spinach main f053e432: add the missing parameters.rspt_order documentation) |
||
| Line 34: | Line 34: | ||
parameters.fwhm - transition full width at half maximum, Tesla | parameters.fwhm - transition full width at half maximum, Tesla | ||
| + | |||
| + | parameters.rspt_order - perturbation theory order to use to account | ||
| + | for the off-diagonal part of the Hamiltonian, | ||
| + | Inf for exact diagonalisation | ||
==Outputs== | ==Outputs== | ||
Latest revision as of 09:36, 30 August 2026
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 Hmw is significant.
Syntax
tran=eigenfields(spin_system,parameters,Hz,Hc,Hmw)
Parameters
Hz - field-dependent part of the laboratory-frame Hamiltonian operator
in Hilbert space, or commutation superoperator in Liouville space,
normalised to 1 Tesla
Hc - field-independent part of the laboratory-frame Hamiltonian operator
in Hilbert space, or commutation superoperator in Liouville space,
containing couplings and offsets
Hmw - observable operator in Hilbert space, or observable vector in
Liouville space, without the amplitude prefactor
parameters.window - magnetic-field window, Tesla
parameters.mw_freq - microwave frequency, Hz
parameters.orientation - three Euler angles in radians specifying the
system orientation
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
parameters.fwhm - transition full width at half maximum, Tesla
parameters.rspt_order - perturbation theory order to use to account
for the off-diagonal part of the Hamiltonian,
Inf for exact diagonalisation
Outputs
tran.tf - vector of transition fields in Tesla tran.tm - vector of transition moments tran.tw - vector of transition FWHMs in Tesla tran.pd - vector of energy level population differences tran.ti - transition identity array, one row per transition tran.tj - vector of scaled field-sweep Jacobians
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, acomm.m, arnoldi.m, atranspose.m, aux_mat.m, binpack.m, cheap_norm.m, cheb_coeff.m, clean_up.m, dirdiff.m, expdrop.m, expmint.m, expmint2.m, fftdiff.m, fourdif.m, fourlap.m, frob_chop.m, gaussfun.m, hdot.m, herm_spline.m, jacobianest.m, keep_rank.m, krondelta.m, kronm_new.m, logfactorial.m, lorentzcon.m, lorentzfun.m, md5_hash.m, mprealloc.m, remncomm.m, remtrace.m, rspert.m, rspt_eig.m, snormpdf.m, svd_shrink.m, tikhoind.m, tikhonov.m, trapdiff.m, unit_oper.m, unit_state.m, vvpert.m, Kernel_utilities
Version 2.6, authors: Ilya Kuprov