deerlab.dipolarmodel¶
- dipolarmodel(t, r=None, Pmodel=None, Bmodel=<deerlab.model.Model object>, experiment=None, parametrization='reftimes', npathways=1, spins=2, harmonics=None, excbandwidth=inf, orisel=None, g=[2.00231930436256, 2.00231930436256], gridsize=1000, minamp=0.001, samespins=True, triangles=None, interp=True)[source]¶
Generate a dipolar EPR signal model.
This function generates a model for the dipolar EPR signal arising from the dipolar interactions between the electron spins of the system. The model is defined in terms of the evolution time vector
t
, the spin-spin distancesr
, the distance distribution modelPmodel
, the intermolecular contributionBmodel
, and anexperiment
object containing information about the specific experiment. The number of spins in the systemspins
can be specified, as well as the parametrization strategyparametrization
for the dipolar pathway refocusing times.For multi-spin systems (
spins > 2
), the function allows for the inclusion of three-spin interactions by specifying thetriangles
list of interspin triangles and enabling thesamespins
assumption of spectral permutability. Theminamp
parameter can be used to threshold the amplitude of the three-spin interactions.Additional effects such as orientation selection or limited excitation bandwidth can be included via the
orisel
andexcbandwidth
arguments, respectively.- Parameters:
- tarray_like
Vector of dipolar evolution times, in microseconds.
- rarray_like, optional
Vector of spin-spin distances, in nanometers. If not specified, it defaults to the range 1.5nm-8nm with a resolution of 0.05nm. Has no effect on models with
spins>2
.- Pmodeldeerlab.Model, optional
Model for the distance distribution. If not specified, a non-parametric distance distribution defined over
r
is used. For models withspins>2
a multivariate Gaussian distance distribution is assumed.- Bmodeldeerlab.Model, optional
Model for the intermolecular (background) contribution. If not specified, a background arising from a homogenous 3D distribution of spins is used.
- experimentdeerlab.ExperimentInfo, optional
Experimental information obtained from experiment models (
ex_
). If specified, the boundaries and start values of the dipolar pathways’ refocusing times and amplitudes are refined based on the specific experiment’s delays.- parametrizationstring, optional
Parametrization strategy of the dipolar pathway refocusing times. Must be one of the following:
'reftimes'
- Each refocusing time is represented individually as a parameter.'delays'
- The pulse delays are introduced as parameters from which the refocusing times are computed. Requiresexperiment
to be specified.'shift'
- A time shift is introduced as a parameter to represent the variability of the refocusing times from their theoretical values. Requiresexperiment
to be specified.
The default is
'reftimes'
.- npathwaysinteger scalar, optional
Number of dipolar pathways. If not specified, a single dipolar pathway is used.
- harmonicslist of integers, optional
Harmonic prefactors of the dipolar pathways. Must be a list with
npathways
harmonics, one for each pathway.- spinsinteger scalar, optional
Number of spins in system. If not specified it defaults to two-spin systems. For multi-spins systems indicated by
spins>2
three-spin interactions are accounted for by the model automatically.- triangleslist of lists, optional
Only required when
spins>3
. List of(Nspins)!/(3!*(Nspins-3)!)
triangles formed by the different interspin distances. Each triangle is specified by a list of indices of the interspin distances forming the triangle, e.g. for a three-spin systemtriangles=[[1,2,3]]
where the first, second, and third interspin distances connect to form a triangle.- samespinsboolean, optional
Only when
spins>3
. Enables the assumption of spectral permutability, i.e. the assumption that all spins in the system have the same spectral distribution. If enabled, all pairwise pathways of the different dipolar interactions are assumed to have the same amplitude. If disabled, all the individual pathway amplitudes are parametrized separately. Enabled by default.- minampfloat scalar, optional
Only when
spins>3
. Threshold amplitude for neglecting three-spin interaction pathways. To enhance performance, all dipolar pathways arising from three-spin interactions with an amplitudelam<minamp
are approximated as having zero amplitude. The default threshold is1e-3
.- oriselcallable or
None
, optional Probability distribution of possible orientations of the interspin vector to account for orientation selection. Must be a function taking a value of the angle θ ∈ [0,π/2] between the interspin vector and the external magnetic field and returning the corresponding probability density. If specified as
None
(default), a uniform distribution is used.- gridsizescalar, optional
Number of points on the grid of powder orientations to be used in the
'grid'
kernel calculation method when evaluating models with orientation selection or multi-spin effects. By default set to 1000 points.- excbandwidthscalar, optional
Excitation bandwidth of the pulses in MHz to account for limited excitation bandwidth. If not specified, an infinite excitation bandwidth is used.
- gscalar, 2-element array, optional
Electron g-values of the two spins,
[g1, g2]
. If a single g is specified,[g, g]
is used. If not specified, g = 2.002319… is used for both spins.- interpboolean, optional
Enable dipolar kernel interpolation for computation time reduction. By default enabled.
- Returns:
- Vmodeldeerlab.Model
Dipolar signal model object.