Mims ENDOR with one nucleus: differences between various Opt.EndorMethod

[stestoll]

There is an incorrect extra peak in the 2H Mims ENDOR of a single 1H nucleus when using Opt.EndorMethod=1.

Reported on the forum in this post.

This was introduced in 5.0.21 when Opt.EndorMethod was introduced. 5.2.12 sets the default value for the single-nucleus case to 0 so that the issue doesn't show, but it is still unresolved.

Example script:

clear, clf, clc

Exp.Sequence = 'MimsENDOR';
Exp.Field = 1224;
Exp.tau = 0.200;
Exp.Range = [6 10];

Opt.nKnots = 101;

Sys.Nucs = '2H';
Sys.A = [-1 -1 2];
Sys.Q = 0.18;
Sys.lwEndor = 0.03;

Opt.EndorMethod = 0;
[x,y0] = saffron(Sys,Exp,Opt);
Opt.EndorMethod = 1;
[x,y1] = saffron(Sys,Exp,Opt);
Opt.EndorMethod = 2;
[x,y2] = saffron(Sys,Exp,Opt);
plot(x,y0/sum(y0),x,y1/sum(y1),x,y2/sum(y2));
legend('0','1','2');