Difference between revisions of "Deer 3p soft deer.m"

From Spinach Documentation Wiki
Jump to: navigation, search
 
Line 1: Line 1:
Three-pulse DEER/PELDOR pulse sequence. The sequence uses soft  
+
{{DISPLAYTITLE:deer_3p_soft_deer.m}} __NOTOC__
pulses computed with the Fokker-Planck formalism. Syntax:
+
Three-pulse DEER/PELDOR pulse sequence. The sequence uses soft pulses computed with the Fokker-Planck formalism.
  
            deer=deer_3p(spin_system,parameters,H,R,K)
+
==Syntax==
  
where H is the Hamiltonian commutation superoperator, R is the
+
    echo_stack=deer_3p_soft_deer(spin_system,parameters,H,R,K)
relaxation superoperator and K is the chemical kinetics super-
 
operator. The following parameters are required:
 
 
 
      parameters.pulse_frq  - frequencies for the
 
                              three pulses, Hz
 
  
 +
==Arguments==
 +
 +
      parameters.pulse_frq  - frequencies for the three
 +
                              pulses, Hz
 +
 
       parameters.pulse_pwr  - power levels for the three
 
       parameters.pulse_pwr  - power levels for the three
                               pulses, Hz
+
                               pulses, rad/s
 
+
 
       parameters.pulse_dur  - durations for the three
 
       parameters.pulse_dur  - durations for the three
 
                               pulses, seconds
 
                               pulses, seconds
 
+
 
       parameters.pulse_phi  - initial phases for the three  
 
       parameters.pulse_phi  - initial phases for the three  
 
                               pulses, radians
 
                               pulses, radians
 
+
 
       parameters.pulse_rnk  - Fokker-Planck ranks for the
 
       parameters.pulse_rnk  - Fokker-Planck ranks for the
 
                               three pulses
 
                               three pulses
 
+
 
       parameters.p1_p3_gap  - time between the first and the
 
       parameters.p1_p3_gap  - time between the first and the
 
                               third pulses, seconds
 
                               third pulses, seconds
 
+
 
       parameters.p2_nsteps  - number of second pulse posi-
 
       parameters.p2_nsteps  - number of second pulse posi-
 
                               tions in the interval between
 
                               tions in the interval between
 
                               the first and the third pulse
 
                               the first and the third pulse
 
+
 
       parameters.echo_time  - time to sample around the ex-
 
       parameters.echo_time  - time to sample around the ex-
 
                               pected echo position
 
                               pected echo position
 
+
 
       parameters.echo_npts  - number of points in the echo
 
       parameters.echo_npts  - number of points in the echo
 
                               discretization
 
                               discretization
 
+
 
       parameters.rho0      - initial state
 
       parameters.rho0      - initial state
 
+
 
       parameters.coil      - detection state
 
       parameters.coil      - detection state
 
+
 
       parameters.method    - soft puse propagation method,
 
       parameters.method    - soft puse propagation method,
 
                               'expv' for Krylov propagation,
 
                               'expv' for Krylov propagation,
Line 45: Line 45:
 
                               gation, 'evolution' for Spin-
 
                               gation, 'evolution' for Spin-
 
                               ach evolution function
 
                               ach evolution function
 +
 +
    H  - Hamiltonian matrix, received from context function
 +
 +
    R  - relaxation superoperator, received from context function
 +
 +
    K  - kinetics superoperator, received from context function
 +
 +
==Outputs==
 +
 +
    deer  - DEER echo stack, a matrix with p2_nsteps echoes
 +
            with echo_npts points each
 +
 +
==Examples==
 +
A number of complete examples for two- and three-electron systems is given in examples/esr_solids folder. A good way to proceed (soft_3_pulse_deer_2e.m example file) is to look at how the pulses affect the system:
 +
 +
[[File:pulse_effect.png]]
 +
 +
then to take a look at the echo and its Fourier transform:
 +
 +
[[File:deer_echo.png]]
 +
 +
and finally at the real and imaginary parts of the DEER trace:
 +
 +
[[File:deer_trace.png]]
 +
 +
==Notes==
 +
# For the method, start with 'expm', change to 'expv' if the calculation runs out of memory, and use 'evolution' as the last resort.
 +
 +
# Simulated echoes tend to be sharp and hard to catch because simulation does not have distributions in experimental      parameters. Fourier transforming the echo prior to integration is recommended.
 +
 +
==See also==
 +
[[deer_3p_hard_echo.m]], [[deer_3p_hard_deer.m]], [[deer_analyt.m]], [[eseem.m]], [[oopeseem.m]], [[deer_3p_soft_diag.m]], [[deer_3p_soft_hole.m]]
 +
 +
 +
''Version 2.4, authors: [[Ilya Kuprov]]''

Revision as of 11:32, 10 July 2019

Three-pulse DEER/PELDOR pulse sequence. The sequence uses soft pulses computed with the Fokker-Planck formalism.

Syntax

    echo_stack=deer_3p_soft_deer(spin_system,parameters,H,R,K)

Arguments

     parameters.pulse_frq  - frequencies for the three 
                             pulses, Hz

     parameters.pulse_pwr  - power levels for the three
                             pulses, rad/s

     parameters.pulse_dur  - durations for the three
                             pulses, seconds

     parameters.pulse_phi  - initial phases for the three 
                             pulses, radians

     parameters.pulse_rnk  - Fokker-Planck ranks for the
                             three pulses

     parameters.p1_p3_gap  - time between the first and the
                             third pulses, seconds

     parameters.p2_nsteps  - number of second pulse posi-
                             tions in the interval between
                             the first and the third pulse

     parameters.echo_time  - time to sample around the ex-
                             pected echo position

     parameters.echo_npts  - number of points in the echo
                             discretization

     parameters.rho0       - initial state

     parameters.coil       - detection state

     parameters.method     - soft puse propagation method,
                             'expv' for Krylov propagation,
                             'expm' for exponential propa-
                             gation, 'evolution' for Spin-
                             ach evolution function

   H  - Hamiltonian matrix, received from context function

   R  - relaxation superoperator, received from context function

   K  - kinetics superoperator, received from context function

Outputs

   deer  - DEER echo stack, a matrix with p2_nsteps echoes
           with echo_npts points each

Examples

A number of complete examples for two- and three-electron systems is given in examples/esr_solids folder. A good way to proceed (soft_3_pulse_deer_2e.m example file) is to look at how the pulses affect the system:

Pulse effect.png

then to take a look at the echo and its Fourier transform:

Deer echo.png

and finally at the real and imaginary parts of the DEER trace:

Deer trace.png

Notes

  1. For the method, start with 'expm', change to 'expv' if the calculation runs out of memory, and use 'evolution' as the last resort.
  1. Simulated echoes tend to be sharp and hard to catch because simulation does not have distributions in experimental parameters. Fourier transforming the echo prior to integration is recommended.

See also

deer_3p_hard_echo.m, deer_3p_hard_deer.m, deer_analyt.m, eseem.m, oopeseem.m, deer_3p_soft_diag.m, deer_3p_soft_hole.m


Version 2.4, authors: Ilya Kuprov