Difference between revisions of "Ngce.m"

From Spinach Documentation Wiki
Jump to: navigation, search
(Sync syntax/arguments/outputs with current Spinach source)
(Update function See also links and function index membership)
 
(3 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
{{DISPLAYTITLE:ngce.m}} __NOTOC__
 
{{DISPLAYTITLE:ngce.m}} __NOTOC__
 
 
Numerical integral route to the Redfield relaxation superoperator.
 
Numerical integral route to the Redfield relaxation superoperator.
  
 
==Syntax==
 
==Syntax==
  
R=ngce(spin_system,H0,H1,dt,tau_est,reg)
+
    R=ngce(spin_system,H0,H1,dt,tau_est,reg)
  
==Arguments==
+
==Parameters==
  
H0 - static laboratory frame Hamiltonian commutation su-
+
  H0 - static laboratory frame Hamiltonian commutation su-
 
       peroperator acting in the background, a matrix
 
       peroperator acting in the background, a matrix
 
   
 
   
Line 29: Line 28:
 
==Outputs==
 
==Outputs==
  
R  - laboratory frame relaxation superoperator
+
  R  - laboratory frame relaxation superoperator
 
 
Note: enough trajectory points must be present to converge
 
      the ensemble averages and Redfield's integral.
 
 
 
Note: the result is returned in the LABORATORY FRAME - eli-
 
      minating non-secular terms is user's responsibility.
 
 
 
ilya.kuprov@weizmann.ac.il
 
jpresteg@uga.edu
 
  
 
==Examples==
 
==Examples==
 
 
See our recent paper with [[Jim Prestegard]]: https://doi.org/10.1016/j.jmr.2020.106891
 
See our recent paper with [[Jim Prestegard]]: https://doi.org/10.1016/j.jmr.2020.106891
  
 
==Notes==
 
==Notes==
 
 
Enough trajectory points must be present to converge the ensemble averages and Redfield's integral.
 
Enough trajectory points must be present to converge the ensemble averages and Redfield's integral.
  
 
==See also==
 
==See also==
 
+
[[relaxation.m]], [[lindbladian.m]], [[magpump.m]], [[adelim.m]], [[blinv.m]], [[blprod.m]], [[corrfun.m]], [[rlx_scalar.m]], [[rlx_split.m]], [[rlx_t1_t2.m]], [[sec2kite.m]], [[spden.m]], [[Kernel_utilities]]
[[relaxation.m]], [[lindbladian.m]], [[magpump.m]]
 
 
 
  
 
''Version 2.6, authors: [[Ilya Kuprov]], [[Jim Prestegard]]''
 
''Version 2.6, authors: [[Ilya Kuprov]], [[Jim Prestegard]]''

Latest revision as of 19:39, 6 June 2026

Numerical integral route to the Redfield relaxation superoperator.

Syntax

    R=ngce(spin_system,H0,H1,dt,tau_est,reg)

Parameters

 H0 - static laboratory frame Hamiltonian commutation su-
      peroperator acting in the background, a matrix

 H1 - stochastic part (zero mean) of the laboratory frame
      Hamiltonian commutation superoperator, a cell array
      of matrices for each point in the MD trajectory.

 dt - time step of the MD trajectory, seconds

 tau_est - correlation tiume estimate for internal safety
           control, seconds

 reg - optional overall relaxation rate, this is added to 
       every eigenvalue of the resulting matrix to prevent
       very small relaxation rates (e.g. singlets) from 
       jumping into positive due to integration accuracy
       limits and then causing problems

Outputs

 R  - laboratory frame relaxation superoperator

Examples

See our recent paper with Jim Prestegard: https://doi.org/10.1016/j.jmr.2020.106891

Notes

Enough trajectory points must be present to converge the ensemble averages and Redfield's integral.

See also

relaxation.m, lindbladian.m, magpump.m, adelim.m, blinv.m, blprod.m, corrfun.m, rlx_scalar.m, rlx_split.m, rlx_t1_t2.m, sec2kite.m, spden.m, Kernel_utilities

Version 2.6, authors: Ilya Kuprov, Jim Prestegard