Difference between revisions of "Ngce.m"
(Update function See also links and function index membership) |
|||
| (9 intermediate revisions by 2 users not shown) | |||
| Line 1: | Line 1: | ||
| − | {{DISPLAYTITLE:ngce.m}} | + | {{DISPLAYTITLE:ngce.m}} __NOTOC__ |
| − | Numerical | + | Numerical integral route to the Redfield relaxation superoperator. |
==Syntax== | ==Syntax== | ||
| − | R=ngce( | + | 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== | ==Notes== | ||
| − | + | 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]], [[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. | + | ''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
