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

From Spinach Documentation Wiki
Jump to: navigation, search
(Sync syntax/arguments/outputs with current Spinach source)
(Restore pre-bulk-edit wiki content)
Line 1: Line 1:
 
{{DISPLAYTITLE:deer_4p_soft_deer.m}} __NOTOC__
 
{{DISPLAYTITLE:deer_4p_soft_deer.m}} __NOTOC__
−
 
 
Four-pulse DEER/PELDOR pulse sequence. The sequence uses soft pulses computed with the Fokker-Planck formalism.
 
Four-pulse DEER/PELDOR pulse sequence. The sequence uses soft pulses computed with the Fokker-Planck formalism.
  
 
==Syntax==
 
==Syntax==
  
−
echo_stack=deer_4p_soft_deer(spin_system,parameters,H,R,K)
+
    echo_stack=deer_4p_soft_deer(spin_system,parameters,H,R,K)
  
 
==Arguments==
 
==Arguments==
−
 
+
−
parameters.pulse_frq  - frequencies for the four
+
    parameters.pulse_frq  - frequencies for the four  
 
                             pulses, Hz
 
                             pulses, Hz
−
 
+
 
     parameters.pulse_pwr  - power levels for the four
 
     parameters.pulse_pwr  - power levels for the four
 
                             pulses, rad/s
 
                             pulses, rad/s
−
 
+
 
     parameters.pulse_dur  - durations for the four
 
     parameters.pulse_dur  - durations for the four
 
                             pulses, seconds
 
                             pulses, seconds
−
 
+
−
     parameters.pulse_phi  - initial phases for the four
+
     parameters.pulse_phi  - initial phases for the four  
 
                             pulses, radians
 
                             pulses, radians
−
 
+
 
     parameters.pulse_rnk  - Fokker-Planck ranks for the
 
     parameters.pulse_rnk  - Fokker-Planck ranks for the
 
                             four pulses
 
                             four pulses
−
 
+
−
     parameters.p1_p2_gap  - time between the end of the
+
     parameters.p1_p2_gap  - time between the end of the  
 
                             first and the start of the
 
                             first and the start of the
 
                             second pulse, seconds
 
                             second pulse, seconds
−
 
+
−
     parameters.p2_p4_gap  - time between the end of the
+
     parameters.p2_p4_gap  - time between the end of the  
 
                             second the start of the third
 
                             second the start of the third
 
                             pulse, seconds
 
                             pulse, seconds
−
 
+
 
     parameters.p3_nsteps  - number of third pulse posi-
 
     parameters.p3_nsteps  - number of third pulse posi-
 
                             tions in the interval between
 
                             tions in the interval between
 
                             the first echo and the fourth
 
                             the first echo and the fourth
 
                             pulse
 
                             pulse
−
 
+
 
     parameters.echo_time  - time to sample around the ex-
 
     parameters.echo_time  - time to sample around the ex-
 
                             pected second echo position
 
                             pected second echo position
−
 
+
 
     parameters.echo_npts  - number of points in the second
 
     parameters.echo_npts  - number of points in the second
 
                             echo discretization
 
                             echo discretization
−
 
+
 
     parameters.rho0      - initial state
 
     parameters.rho0      - initial state
−
 
+
 
     parameters.coil      - detection state
 
     parameters.coil      - detection state
−
 
+
−
    parameters.spins      - irradiated spins, normally {'E'}
 
−
 
 
 
     parameters.method    - soft puse propagation method,
 
     parameters.method    - soft puse propagation method,
 
                             'expv' for Krylov propagation,
 
                             'expv' for Krylov propagation,
Line 54: Line 51:
 
                             gation, 'evolution' for Spin-
 
                             gation, 'evolution' for Spin-
 
                             ach evolution function
 
                             ach evolution function
−
 
+
 
   H  - Hamiltonian matrix, received from context function
 
   H  - Hamiltonian matrix, received from context function
−
 
+
 
   R  - relaxation superoperator, received from context function
 
   R  - relaxation superoperator, received from context function
−
 
+
 
   K  - kinetics superoperator, received from context function
 
   K  - kinetics superoperator, received from context function
  
 
==Outputs==
 
==Outputs==
  
−
echo_stack  - DEER echo stack, a matrix with p3_nsteps echoes
+
    echo_stack  - DEER echo stack, a matrix with p3_nsteps echoes
 
                   with echo_npts points each
 
                   with echo_npts points each
−
 
−
Note: for the method, start with 'expm', change to 'expv' if the
 
−
      calculation runs out of memory, and use 'evolution' as the
 
−
      last resort.
 
−
 
−
Note: simulated echoes tend to be sharp and hard to catch becau-
 
−
      se simulation does not have distributions in experimental
 
−
      parameters. Fourier transforming the echo prior to integ-
 
−
      ration is recommended.
 
−
 
−
Note: the time in the DEER trace refers to the third pulse inser-
 
−
      tion point, after end of the second pulse.
 
−
 
−
ilya.kuprov@weizmann.ac.il
 
  
 
==Examples==
 
==Examples==
−
 
 
A number of complete examples for two- and three-electron systems are given in examples/esr_solids folder. A good way to proceed (soft_4_pulse_deer_2e.m example file) is to look at how the pulses affect the system:
 
A number of complete examples for two- and three-electron systems are given in examples/esr_solids folder. A good way to proceed (soft_4_pulse_deer_2e.m example file) is to look at how the pulses affect the system:
  
Line 95: Line 77:
  
 
==Notes==
 
==Notes==
−
 
 
# For the method, start with 'expm', change to 'expv' if the calculation runs out of memory, and use 'evolution' as the last resort.
 
# 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.
 
# 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.
Line 101: Line 82:
  
 
==See also==
 
==See also==
−
 
 
[[deer_3p_hard_echo.m]], [[deer_3p_hard_deer.m]], [[deer_analyt.m]], [[eseem.m]], [[oopeseem.m]], [[deer_4p_soft_hole.m]], [[deer_4p_soft_diag.m]], [[deer_3p_soft_deer.m]], [[deer_3p_soft_hole.m]], [[deer_3p_soft_diag.m]]
 
[[deer_3p_hard_echo.m]], [[deer_3p_hard_deer.m]], [[deer_analyt.m]], [[eseem.m]], [[oopeseem.m]], [[deer_4p_soft_hole.m]], [[deer_4p_soft_diag.m]], [[deer_3p_soft_deer.m]], [[deer_3p_soft_hole.m]], [[deer_3p_soft_diag.m]]
  
  
 
''Version 2.4, authors: [[Ilya Kuprov]]''
 
''Version 2.4, authors: [[Ilya Kuprov]]''

Revision as of 15:48, 5 April 2026

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

Syntax

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

Arguments

   parameters.pulse_frq  - frequencies for the four 
                           pulses, Hz

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

   parameters.pulse_dur  - durations for the four
                           pulses, seconds

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

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

   parameters.p1_p2_gap  - time between the end of the 
                           first and the start of the
                            second pulse, seconds

   parameters.p2_p4_gap  - time between the end of the 
                           second the start of the third
                           pulse, seconds

   parameters.p3_nsteps  - number of third pulse posi-
                           tions in the interval between
                           the first echo and the fourth
                           pulse

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

   parameters.echo_npts  - number of points in the second
                           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

   echo_stack  - DEER echo stack, a matrix with p3_nsteps echoes
                 with echo_npts points each

Examples

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

Pulse effect 4p.png

then to inspect the echo stack as a function of the pump pulse position:

Deer echo stack 4p.png

and finally to inspect the principal components of the echo stack. The most prominent component is the DEER trace, the other components come from the finite size of the spherical integration grid.

Deer echo 4p.png Deer trace 4p.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.
  2. 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.
  3. The time in the DEER trace refers to the second pulse insertion point, after end of first pulse.

See also

deer_3p_hard_echo.m, deer_3p_hard_deer.m, deer_analyt.m, eseem.m, oopeseem.m, deer_4p_soft_hole.m, deer_4p_soft_diag.m, deer_3p_soft_deer.m, deer_3p_soft_hole.m, deer_3p_soft_diag.m


Version 2.4, authors: Ilya Kuprov