Dnp field scan.m

Magnetic field scan DNP experiment. Returns the equilibrium magnetization of all spins as a function of magnetic field. Syntax:

    answer=dnp_field_scan(spin_system,parameters,H,R,K)

where H is the Hamiltonian commutation superoperator, R is the relaxation superoperator and K is the chemical kinetics superoperator. The following parameters are required:

  parameters.mw_pwr      -   microwave power, Hz

  parameters.mw_frq      -   microwave frequency, Hz

  parameters.fields      -   a vector of magnetic field off-
                             sets from the B0 field, T

  parameters.rho_eq      -   equilibrium state at B0

  parameters.coil        -   coil state vector or a horizon-
                             tal stack thereof

  parameters.mw_oper     -   microwave irradiation operator

  parameters.ez_oper     -   Lz operator on the electrons

  parameters.method      -   'backslash' to use Matlab's
                             linear equation solver, 'gmres'
                             to use ILU preconditioned GMRES

Note: the relaxation superoperator should NOT be thermalized for this type of calculation.

Note: thermal equilibrium state is assumed to be the same at all fields in the sweep - do not use with broad magnetic field sweep experiments.