Difference between revisions of "Ngce.m"
(Link contributor names to author pages) |
(Sync syntax/arguments/outputs with current Spinach source) |
||
| 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) | |
==Arguments== | ==Arguments== | ||
| − | + | H0 - static laboratory frame Hamiltonian commutation su- | |
peroperator acting in the background, a matrix | peroperator acting in the background, a matrix | ||
| Line 28: | Line 29: | ||
==Outputs== | ==Outputs== | ||
| − | + | 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]] | [[relaxation.m]], [[lindbladian.m]], [[magpump.m]] | ||
''Version 2.6, authors: [[Ilya Kuprov]], [[Jim Prestegard]]'' | ''Version 2.6, authors: [[Ilya Kuprov]], [[Jim Prestegard]]'' | ||
Revision as of 15:06, 5 April 2026
Numerical integral route to the Redfield relaxation superoperator.
Syntax
R=ngce(spin_system,H0,H1,dt,tau_est,reg)
Arguments
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
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
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
Version 2.6, authors: Ilya Kuprov, Jim Prestegard
