Solid effect.m

Solid effect DNP experiment. Computes either the time dependence of the longitudinal magnetization in the system or its asymptotic amplitude. Syntax:

              answer=solid_effect(spin_system,parameters)
 

Parameters:

        parameters.mw_power   - microwave power in rad/s

          parameters.theory   - level of theory. Set to 'exact' for 
                                the electron rotating frame calcula-
                                tion or to any of the following six 
                                options for the average Hamiltonian
                                theory calculation on top of the el-
                                ectron + nuclear rotating frame:
                                'ah_first_order', 'ah_second_order',
                                'ah_third_order', 'kb_first_order',
                                'kb_second_order', 'kb_third_order'.
                                See the header of average.m functi-
                                on for the meaning of these options.

     parameters.nuclear_frq   - nuclear Zeeman frequency in rad/s

       parameters.calc_type   - set to 'time_dependence' to get the 
                                time dependence of the longitudinal
                                magnetization and to 'steady_state'
                                to get the asymptotic longitudinal
                                magnetization.

       parameters.time_step   - if 'time_dependence' is set in the
                                calc_type parameter, sets the time
                                step, seconds.

         parameters.n_steps   - if 'time_dependence' is set in the
                                calc_type parameter, sets the num-
                                ber of time steps.

Note: the system must have exactly one electron.

Note: this function generates its own Liouvillian and should be called directly, without "liquid", "powder", or other wrappers.