Setup

ColorBlindSafeColorScheme=Table[RGBColor[i],{i,{"#000000","#4477AA","#EE6677","#228833","#CCBB44","#66CCEE","#AA3377","#BBBBBB"}}];
SetOptions[Plot,PlotStyle->Table[Directive[Thick,n],{n,ColorBlindSafeColorScheme}],FrameStyle->Directive[Black,AbsoluteThickness[1.5],FontSize->30,FontWeight->"Light",FontFamily->"Helvetica"],ImageSize->600];
SetOptions[ListPlot,PlotStyle->Table[Directive[Thick,n],{n,ColorBlindSafeColorScheme}],FrameStyle->Directive[Black,AbsoluteThickness[1.5],FontSize->30,FontWeight->"Light",FontFamily->"Helvetica"],ImageSize->600];SetOptions[ListLinePlot,PlotStyle->Table[Directive[Thick,n],{n,ColorBlindSafeColorScheme}],AxesStyle->Directive[Black,AbsoluteThickness[1.5],FontSize->30,FontWeight->"Light",FontFamily->"Helvetica"],FrameStyle->Directive[Black,AbsoluteThickness[1.5],FontSize->30,FontWeight->"Light",FontFamily->"Helvetica"],ImageSize->600]

Needs["MaTeX`"]

SetOptions[MaTeX,FontSize->36]

SD

Needs["SpinDynamica`"]

SetSpinSystem[2];

SetBasis[SingletTripletBasis[]]

spin system

H = 2π (J opI[1].opI[2]+Δ (((opI[1,"z"]-opI[2,"z"]))/2));

J=255.94;
Δ=17.8;

INADEQUATE plot (Fig. 1)

func2[t_,θ_]=(((Sin[2 π t] (1-2 ((Cos[θST])^2) ((Sin[π t Sec[θST]])^2))-Cos[2 π t] Cos[θST] Sin[2 π t Sec[θST]]))^2)/.θST->π/180 θ

[WARNING: IMAGE SIZE IS VERY LARGE. RASTERIZED FOR UPLOADABILITY]

Rasterize[Show[{DensityPlot[Evaluate[func2[t,θ]],{θ,0,90},{t,0,2},PlotRange->{{-0.5,90+0.5},{-0.02,2+0.02}},ColorFunctionScaling->False,ColorFunction->(Blend[{White,Black},#]&),PlotLegends->Placed[BarLegend[{(Blend[{White,Black},#]&),{0,1}},None,Method->{Frame->True,FrameStyle->Directive[AbsoluteThickness[1.5]],TicksStyle->Directive[Black,AbsoluteThickness[1.5]],FrameTicks->{0,0.2},ImagePadding -> None},LabelStyle->{Black,FontSize->36,FontFamily->"Helvetica",FontWeight->"Light"},LegendLayout->"Column",LegendMarkerSize->{30,512},LegendLabel->],{{1,0.562}}],PlotRangePadding->0,FrameLabel->{{,None},{,}},FrameTicks->{{{{0.`,"0.0"},{0.5`,"0.5"},{1.`,"1.0"},{1.5`,"1.5"},{2.`,"2.0"}},True},{{{0,0},{20,20},{40,40},{60,60},{80,80}},{{5.710593137499643`,"0.1"},{26.56505117707799`,"0.5"},{45,"1"},{(180 ArcTan[2])/π,"2"},{(180 ArcTan[10])/π,"10"}}}},PlotPoints->{300,300},ImageSize->800,MaxRecursion->4,Mesh->None,FrameStyle->Directive[Black,AbsoluteThickness[1.5],FontSize->36,FontWeight->"Light",FontFamily->"Helvetica"]],ListPlot[{{87,1/4}},PlotMarkers->★,PlotRange->{{0,90},{0,2}}]}]]

Simulated trajectories (Fig.5)

SLIC[ω1_,τSLIC_]={ω1 opI["x"],τSLIC}

κPulsePol=(((2(1+Sqrt[2])))/(3 π))//N

PulsePol[τR_,nR_]:=If[EvenQ[nR]==True,PulseSequence[{
Repeat[{{Pulse[{π/2,3π/4}],{None,τR/2},Pulse[{π,π/4}],{None,τR/2},Pulse[{π/2,3π/4}],Pulse[{π/2,-3π/4}],{None,τR/2},Pulse[{π,3π/4}],{None,τR/2},Pulse[{π/2,-3π/4}]}},nR/2]
}],PulseSequence[{
Repeat[{{Pulse[{π/2,3π/4}],{None,τR/2},Pulse[{π,π/4}],{None,τR/2},Pulse[{π/2,3π/4}],Pulse[{π/2,-3π/4}],{None,τR/2},Pulse[{π,3π/4}],{None,τR/2},Pulse[{π/2,-3π/4}]}},Floor[nR/2]],Pulse[{π/2,3π/4}],{None,τR/2},Pulse[{π,π/4}],{None,τR/2},Pulse[{π/2,3π/4}]
}]]

ppoltab1=Table[{nR 3/(4 255.94),TransformationAmplitude[opI["x"]->-2opI[1,"x"].opI[2,"z"]-2opI[1,"z"].opI[2,"x"],{PulsePol[3/(4 255.94),nR]},BackgroundGenerator->H]},{nR,0,144,1}];

slictab1=Table[{tslic,TransformationAmplitude[opI["z"]->-2opI[1,"x"].opI[2,"y"]-2opI[1,"y"].opI[2,"x"],{{2π 255.94 opI["x"],tslic}},BackgroundGenerator->H]},{tslic,0+1/(4 255.94),108/255.94+1/(4 255.94),1/(4 255.94)}];

Clear[custmplotmark]
custmplotmark[name_String,Fcol_,Opact_:White,size_:6]:=Show[ResourceFunction["PolygonMarker"][name,Offset@size,{FaceForm[Opact],Dynamic@EdgeForm[{CurrentValue["Color"],AbsoluteThickness[1.618],Opacity[1]}]}],ResourceFunction["PolygonMarker"][name,Offset[Fcol size],EdgeForm[None]]];

ResourceFunction["PlotGrid"][{{Show[Plot[{((Sin[κPulsePol ωΔ t/4])^2)Cos[2(π/4) J (4/3)t]/.J->255.94/.ωΔ->2π 17.8,((Sin[κPulsePol ωΔ t/4])^2)/.ωΔ->2π 17.8},{t,0,0.42},PlotStyle->{Directive[Gray],Directive[Blue,Thickness[0.01]]},PlotRange->{{0,0.42},{1.1,-1.1}},PlotPoints->200],ListPlot[ppoltab1,Joined->False,PlotMarkers->{custmplotmark["Circle",0,Opacity[0.95,White],5]}],FrameLabel->{{,None},{"T / s",""}},FrameTicks->{{{{-1,"-1.0"},{0,"0.0"},{1,"1.0"}},Automatic},{{0,0.1,0.2,0.3,0.4},Table[{3/ 255.94 m,m},{m,0,30,10}]}}]},{Show[Plot[{Sin[2π J t]((Sin[(1/(Sqrt[2])) ωΔ t/4])^2)/.ωΔ->2π 17.8/.J->255.94,((Sin[(1/(Sqrt[2])) ωΔ t/4])^2)/.ωΔ->2π 17.8},{t,0,0.42},PlotStyle->{Directive[Gray],Directive[Blue,Thickness[0.01]]},PlotRange->{{0,0.42},{1.1,-1.1}},PlotPoints->200,FrameLabel->{{,None},{"T / s",""}},FrameTicks->{{{{-1,"-1.0"},{0,"0.0"},{1,"1.0"}},True},{{0,0.1,0.2,0.3,0.4},Table[{1/ 255.94 n,n},{n,{0,20,40,60,80,100}}]}}],ListPlot[Transpose[{Downsample[slictab1[[All,1]],1],Downsample[slictab1[[All,2]],1]}],Joined->False,PlotMarkers->{custmplotmark["Circle",0,Opacity[0.95,White],5]},PlotRange->{{0,0.42},{1.1,-1.1}},FrameLabel->{{,None},{"T / s",""}},FrameTicks->{{{{-1,"-1.0"},{0,"0.0"},{1,"1.0"}},True},{{0,0.1,0.2,0.3,0.4},Table[{1/ 255.94 n,n},{n,{0,20,40,60,80,100}}]}}]]}},Spacings->{0,200},ImageSize->600,AspectRatio->5/3,"ShowFrameLabels"->Full,PlotLabels->Placed[Row->"a",{Left,Above}],LabelStyle->{Black,35,FontFamily->"Helvetica"}]

experimental data vs simulations (Fig. 6)

BB1[{β_,ϕ_}]={Pulse[{β/2,ϕ}],Pulse[{π,ArcCos[-β/(4π)]+ϕ}],Pulse[{2π,3ArcCos[-β/(4π)]+ϕ}],Pulse[{π,ArcCos[-β/(4π)]+ϕ}],Pulse[{β/2,ϕ}]}

PulsePolWithBB1[τR_,nR_,ωnut_]:=PulseSequence[{Repeat[{{Pulse[{π/2,3π/4}],{None,τR/2-((3π)/( ωnut))},Sequence@@BB1[{π,π/4}],{None,τR/2-((3π)/(ωnut))},Pulse[{π/2,3π/4}],Pulse[{π/2,-3π/4}],{None,τR/2-((3π)/(ωnut))},Sequence@@BB1[{π,3π/4}],{None,τR/2-((3π)/(ωnut))},Pulse[{π/2,-3π/4}]}},nR/2]},NutationFrequency->ωnut]

pulsepolbuildups="Dataset :
 C:/Users/Mohamed/Downloads/20251024_19F2/862/pdata/1
AREA fit :
 I[t]=I[0]+P*exp(-t/T1)

91 points for Integral 1,  Integral Region from -102.847 to -103.018 ppm
Results     Comp. 1    

I[0]  =   -6.045e+03
P     =    6.045e+03
T1    =    1216.040s
SD    =    2.540e-01

    tau    ppm     integral    intensity

    2.000m  -102.929   4.9491e+07     eliminated
    3.000m  -102.929   6.3594e+07     eliminated
    4.000m  -102.929   9.2086e+05     eliminated
    5.000m  -102.929   2.1805e+08     eliminated
    6.000m  -102.929   6.7514e+08     eliminated
    7.000m  -102.929   1.1333e+09     eliminated
    8.000m  -102.929   2.7473e+09   3.7225e+07
    9.000m  -102.929   4.3201e+09   5.7007e+07
   10.000m  -102.929   6.7509e+09     eliminated
   11.000m  -102.929   9.3942e+09     eliminated
   12.000m  -102.930   1.3036e+10   1.6748e+08
   13.000m  -102.930   1.6658e+10   2.1615e+08
   14.000m  -102.930   2.0601e+10   2.6752e+08
   15.000m  -102.929   2.4324e+10   3.1546e+08
   16.000m  -102.929   2.7951e+10   3.6164e+08
   17.000m  -102.929    3.099e+10    4.012e+08
   18.000m  -102.930   3.3357e+10   4.3027e+08
   19.000m  -102.929   3.4818e+10   4.4869e+08
   20.000m  -102.929   3.5326e+10   4.5536e+08
   21.000m  -102.930   3.4813e+10   4.4995e+08
   22.000m  -102.929   3.3491e+10   4.3319e+08
   23.000m  -102.929   3.1219e+10   4.0441e+08
   24.000m  -102.929    2.841e+10   3.6738e+08
   25.000m  -102.929   2.5023e+10   3.2046e+08
   26.000m  -102.929   2.1321e+10   2.7279e+08
   27.000m  -102.929   1.7455e+10   2.2457e+08
   28.000m  -102.929   1.3976e+10    1.802e+08
   29.000m  -102.930   1.0611e+10   1.3746e+08
   30.000m  -102.929   7.4952e+09   9.6134e+07
   31.000m  -102.929   5.1724e+09   6.5506e+07
   32.000m  -102.929   3.1853e+09   4.0219e+07
   33.000m  -102.930   1.8205e+09   2.2151e+07
   34.000m  -102.929   1.0731e+09     eliminated
   35.000m  -102.929    6.101e+08     eliminated
   36.000m  -102.929   1.5068e+08     eliminated
   37.000m  -102.929   1.3872e+08     eliminated
   38.000m  -102.929    7.972e+07     eliminated
   39.000m  -102.929   4.0073e+07     eliminated
   40.000m  -102.929   4.8668e+07     eliminated
   41.000m  -102.929   4.1002e+07     eliminated
   42.000m  -102.929   2.4416e+08     eliminated
   43.000m  -102.929    8.794e+07     eliminated
   44.000m  -102.929   1.7344e+08     eliminated
   45.000m  -102.929   2.2682e+08     eliminated
   46.000m  -102.929   5.0733e+08   6.5118e+06
   47.000m  -102.929   9.6536e+08     eliminated
   48.000m  -102.929   1.8415e+09     eliminated
   49.000m  -102.929   2.6989e+09     eliminated
   50.000m  -102.929   3.8336e+09     eliminated
   51.000m  -102.929   5.0437e+09     eliminated
   52.000m  -102.929   6.3809e+09     eliminated
   53.000m  -102.929    7.822e+09     eliminated
   54.000m  -102.929   9.0024e+09     eliminated
   55.000m  -102.929    1.003e+10     eliminated
   56.000m  -102.929    1.084e+10     eliminated
   57.000m  -102.929   1.1287e+10     eliminated
   58.000m  -102.929    1.141e+10     eliminated
   59.000m  -102.929   1.1142e+10     eliminated
   60.000m  -102.929   1.0632e+10     eliminated
   61.000m  -102.929   9.7796e+09     eliminated
   62.000m  -102.929   8.7075e+09     eliminated
   63.000m  -102.929   7.4788e+09     eliminated
   64.000m  -102.929    6.241e+09     eliminated
   65.000m  -102.929   4.9141e+09     eliminated
   66.000m  -102.929   3.7547e+09     eliminated
   67.000m  -102.929    2.761e+09     eliminated
   68.000m  -102.929   1.8413e+09     eliminated
   69.000m  -102.929   1.1632e+09     eliminated
   70.000m  -102.929   6.4914e+08     eliminated
   71.000m  -102.929   3.1959e+08     eliminated
   72.000m  -102.929   1.3859e+08     eliminated
   73.000m  -102.929   6.9635e+07     eliminated
   74.000m  -102.929   1.5693e+08  -1.9159e+06
   75.000m  -102.930   6.4566e+07  -1.5951e+06
   76.000m  -102.929   1.8146e+08     eliminated
   77.000m  -102.930   1.3331e+08   -2.326e+06
   78.000m  -102.929   1.0839e+08     eliminated
   79.000m  -102.929   4.2838e+07     eliminated
   80.000m  -102.930  -5.8038e+07  -3.0425e+06
   81.000m  -102.929  -5.0892e+07  -2.2342e+06
   82.000m  -102.929   8.8075e+07     eliminated
   83.000m  -102.929   7.0647e+07     eliminated
   84.000m  -102.929    2.704e+08     eliminated
   85.000m  -102.929   5.3103e+08     eliminated
   86.000m  -102.929   7.7175e+08     eliminated
   87.000m  -102.929   1.2331e+09     eliminated
   88.000m  -102.929   1.7056e+09     eliminated
   89.000m  -102.929   2.0989e+09     eliminated
   90.000m  -102.929   2.5918e+09     eliminated
   91.000m  -102.929   3.0251e+09     eliminated
   92.000m  -102.929   3.3664e+09     eliminated
";

slicbuildups="Dataset :
 C:/Users/Mohamed/Downloads/20251024_19F2/867/pdata/1
AREA fit :
 I[t]=I[0]+P*exp(-t/T1)

128 points for Integral 1,  Integral Region from -102.847 to -103.018 ppm
Results     Comp. 1    

I[0]  =   -1.621e-01
P     =    8.599e-01
T1    =      41.076m
SD    =    2.326e-01

    tau    ppm     integral    intensity

    9.000m  -102.930    9.943e+06     eliminated
   10.000m  -102.930   2.4489e+08     eliminated
   11.000m  -102.930   -1.916e+08     eliminated
   12.000m  -102.930   9.0253e+08     eliminated
   13.000m  -102.930   9.1347e+08     eliminated
   14.000m  -102.930   3.2358e+09     eliminated
   15.000m  -102.930   4.9214e+09     eliminated
   16.000m  -102.930   7.1957e+09   9.8521e+07
   17.000m  -102.930    1.119e+10   1.5532e+08
   18.000m  -102.930   1.4162e+10     eliminated
   19.000m  -102.930   1.8274e+10   2.6207e+08
   20.000m  -102.930    2.191e+10   3.1227e+08
   21.000m  -102.930   2.5397e+10   3.6218e+08
   22.000m  -102.930   2.9026e+10     4.11e+08
   23.000m  -102.930   3.1646e+10    4.447e+08
   24.000m  -102.930   3.4136e+10   4.7405e+08
   25.000m  -102.930   3.6056e+10   4.9733e+08
   26.000m  -102.930   3.7449e+10   5.0563e+08
   27.000m  -102.930   3.7775e+10   5.0182e+08
   28.000m  -102.930   3.6157e+10   4.7204e+08
   29.000m  -102.930   3.4133e+10   4.4084e+08
   30.000m  -102.930   3.1157e+10   4.0339e+08
   31.000m  -102.930   2.7869e+10   3.6355e+08
   32.000m  -102.930   2.3637e+10   3.1277e+08
   33.000m  -102.930    1.969e+10   2.6624e+08
   34.000m  -102.930   1.6169e+10   2.1926e+08
   35.000m  -102.930   1.2551e+10   1.7482e+08
   36.000m  -102.930   9.8315e+09   1.4021e+08
   37.000m  -102.930   7.9958e+09   1.0997e+08
   38.000m  -102.930   6.5083e+09   8.6801e+07
   39.000m  -102.930   5.5225e+09   6.8889e+07
   40.000m  -102.931   4.5145e+09     5.29e+07
   41.000m  -102.931   3.5363e+09     4.02e+07
   42.000m  -102.930   3.5427e+09     eliminated
   43.000m  -102.930   2.7772e+09     eliminated
   44.000m  -102.930   2.2161e+09     eliminated
   45.000m  -102.930   9.3243e+08     eliminated
   46.000m  -102.930  -4.6848e+08     eliminated
   47.000m  -102.930  -1.1986e+09     eliminated
   48.000m  -102.930  -1.8223e+09     eliminated
   49.000m  -102.930  -2.3868e+09  -2.3809e+07
   50.000m  -102.930  -3.4173e+09     eliminated
   51.000m  -102.930  -4.1695e+09  -4.7759e+07
   52.000m  -102.930  -4.9888e+09     eliminated
   53.000m  -102.930  -6.1315e+09  -7.7172e+07
   54.000m  -102.930  -6.7436e+09  -8.5806e+07
   55.000m  -102.930  -7.1741e+09  -9.2202e+07
   56.000m  -102.930  -7.6902e+09  -9.8594e+07
   57.000m  -102.930  -7.6408e+09     eliminated
   58.000m  -102.931  -6.8346e+09  -8.9842e+07
   59.000m  -102.930  -7.2202e+09  -9.7148e+07
   60.000m  -102.931  -7.3126e+09  -9.6543e+07
   61.000m  -102.930  -6.5283e+09  -8.7584e+07
   62.000m  -102.930   -6.089e+09  -8.0767e+07
   63.000m  -102.930  -6.0076e+09     eliminated
   64.000m  -102.930  -5.6147e+09  -7.6003e+07
   65.000m  -102.930  -5.4205e+09     eliminated
   66.000m  -102.930  -5.4976e+09   -7.136e+07
   67.000m  -102.930  -5.4733e+09  -6.9982e+07
   68.000m  -102.930  -5.3851e+09  -6.8828e+07
   69.000m  -102.930  -5.3878e+09  -6.8509e+07
   70.000m  -102.930  -5.4886e+09  -7.1186e+07
   71.000m  -102.931   -5.485e+09  -7.2257e+07
   72.000m  -102.930  -5.6055e+09     eliminated
   73.000m  -102.930  -5.2739e+09     eliminated
   74.000m  -102.931  -5.1098e+09  -6.9894e+07
   75.000m  -102.930  -5.3757e+09     eliminated
   76.000m  -102.930  -5.0745e+09  -7.1107e+07
   77.000m  -102.930  -4.7248e+09     eliminated
   78.000m  -102.930  -4.7807e+09    -6.77e+07
   79.000m  -102.930  -4.4111e+09  -6.3784e+07
   80.000m  -102.930  -4.5512e+09     eliminated
   81.000m  -102.931  -4.1863e+09   -5.961e+07
   82.000m  -102.931  -4.0736e+09  -5.8261e+07
   83.000m  -102.930  -3.9598e+09     eliminated
   84.000m  -102.930  -3.7287e+09     eliminated
   85.000m  -102.930  -3.7171e+09     eliminated
   86.000m  -102.931  -3.5266e+09  -4.9646e+07
   87.000m  -102.930  -3.4338e+09     eliminated
   88.000m  -102.930  -3.4218e+09  -4.8945e+07
   89.000m  -102.930  -3.2252e+09  -4.8111e+07
   90.000m  -102.930  -3.2514e+09     eliminated
   91.000m  -102.930  -3.5212e+09     eliminated
   92.000m  -102.930  -3.3903e+09  -4.9613e+07
   93.000m  -102.931  -3.2189e+09  -4.8306e+07
   94.000m  -102.930  -3.2823e+09  -4.4842e+07
   95.000m  -102.930  -3.2275e+09  -4.7153e+07
   96.000m  -102.930  -3.1901e+09     eliminated
   97.000m  -102.930  -3.0495e+09     eliminated
   98.000m  -102.930   -3.133e+09     eliminated
   99.000m  -102.930  -2.8453e+09     eliminated
  100.000m  -102.930  -2.9124e+09     eliminated
  101.000m  -102.930  -2.8027e+09     eliminated
  102.000m  -102.930  -2.5563e+09     eliminated
  103.000m  -102.930  -2.6408e+09     eliminated
  104.000m  -102.930  -2.5166e+09  -3.3169e+07
  105.000m  -102.930  -2.3181e+09  -3.1187e+07
  106.000m  -102.930  -2.2693e+09  -2.9684e+07
  107.000m  -102.930   -2.292e+09     eliminated
  108.000m  -102.930  -2.1185e+09   -3.001e+07
  109.000m  -102.930  -2.0771e+09     eliminated
  110.000m  -102.930  -1.9878e+09     eliminated
  111.000m  -102.930  -1.8224e+09  -2.4667e+07
  112.000m  -102.930  -2.0274e+09  -2.6135e+07
  113.000m  -102.930  -1.8456e+09     eliminated
  114.000m  -102.930  -1.9048e+09  -2.5799e+07
  115.000m  -102.931  -2.0579e+09  -2.7266e+07
  116.000m  -102.930  -1.8003e+09     eliminated
  117.000m  -102.930  -1.7979e+09     eliminated
  118.000m  -102.930  -1.7828e+09     eliminated
  119.000m  -102.930  -1.6494e+09  -2.2382e+07
  120.000m  -102.930  -1.5629e+09  -2.3487e+07
  121.000m  -102.930  -1.7428e+09     eliminated
  122.000m  -102.930   -1.467e+09     eliminated
  123.000m  -102.930  -1.4154e+09       -2e+07
  124.000m  -102.931  -1.3679e+09  -1.7651e+07
  125.000m  -102.930    -1.26e+09     eliminated
  126.000m  -102.930  -1.2532e+09     eliminated
  127.000m  -102.930  -1.1016e+09   -1.624e+07
  128.000m  -102.930  -1.0438e+09     eliminated
  129.000m  -102.930  -9.1216e+08     eliminated
  130.000m  -102.930  -9.3399e+08  -1.3836e+07
  131.000m  -102.930  -8.2902e+08     eliminated
  132.000m  -102.930  -8.7654e+08  -1.1493e+07
  133.000m  -102.930  -8.7599e+08  -1.1844e+07
  134.000m  -102.930   -8.181e+08  -1.2165e+07
  135.000m  -102.930  -8.2186e+08     eliminated
  136.000m  -102.931  -7.0804e+08  -1.3745e+07";

cslicbuildups="Dataset :
 C:/Users/Mohamed/Downloads/20251024_19F2/868/pdata/1
AREA fit :
 I[t]=I[0]+P*exp(-t/T1)

152 points for Integral 1,  Integral Region from -102.847 to -103.018 ppm
Results     Comp. 1    

I[0]  =   -7.175e-02
P     =    5.516e-01
T1    =      41.673m
SD    =    2.168e-01

    tau    ppm     integral    intensity

    2.000m  -102.930   4.9875e+07     eliminated
    3.000m  -102.929   4.4636e+07  -2.2362e+06
    4.000m  -102.930   8.0469e+07     eliminated
    5.000m  -102.930    4.652e+07     eliminated
    6.000m  -102.930   4.4028e+08     eliminated
    7.000m  -102.930   1.1008e+09     eliminated
    8.000m  -102.930   2.0815e+09   2.5873e+07
    9.000m  -102.930    3.431e+09     eliminated
   10.000m  -102.930   5.1485e+09   6.8391e+07
   11.000m  -102.930   7.6072e+09    1.012e+08
   12.000m  -102.930     1.06e+10   1.3854e+08
   13.000m  -102.930   1.3711e+10   1.8141e+08
   14.000m  -102.930   1.7207e+10   2.2806e+08
   15.000m  -102.930   2.1016e+10   2.7685e+08
   16.000m  -102.930   2.4799e+10   3.2781e+08
   17.000m  -102.930   2.8211e+10   3.7322e+08
   18.000m  -102.930   3.0863e+10    4.124e+08
   19.000m  -102.930   3.3214e+10   4.4495e+08
   20.000m  -102.930   3.4651e+10   4.6283e+08
   21.000m  -102.930   3.5163e+10   4.7546e+08
   22.000m  -102.930   3.4751e+10   4.6894e+08
   23.000m  -102.930   3.3242e+10   4.5109e+08
   24.000m  -102.930   3.1119e+10   4.2419e+08
   25.000m  -102.930   2.8298e+10   3.8603e+08
   26.000m  -102.930   2.4879e+10   3.4186e+08
   27.000m  -102.930   2.1073e+10     2.92e+08
   28.000m  -102.930   1.7085e+10   2.3847e+08
   29.000m  -102.930   1.2978e+10     eliminated
   30.000m  -102.930   9.1815e+09    1.322e+08
   31.000m  -102.930   6.0567e+09     eliminated
   32.000m  -102.930   4.1845e+09   6.3806e+07
   33.000m  -102.930   1.9432e+09     eliminated
   34.000m  -102.930   3.3289e+08     eliminated
   35.000m  -102.930  -6.0351e+08     eliminated
   36.000m  -102.930  -1.3127e+09     eliminated
   37.000m  -102.930  -1.7176e+09  -2.0885e+07
   38.000m  -102.930  -1.7021e+09  -2.1171e+07
   39.000m  -102.930  -1.2134e+09     eliminated
   40.000m  -102.930  -6.3221e+08  -7.6211e+06
   41.000m  -102.930  -1.4208e+08     eliminated
   42.000m  -102.930   5.3008e+08     eliminated
   43.000m  -102.929   2.9122e+08   5.3932e+06
   44.000m  -102.930   4.9287e+08     eliminated
   45.000m  -102.930   2.5717e+08     eliminated
   46.000m  -102.930  -1.4569e+08     eliminated
   47.000m  -102.930   4.0344e+08     eliminated
   48.000m  -102.930  -7.2283e+07     eliminated
   49.000m  -102.930  -5.0373e+08     eliminated
   50.000m  -102.929  -9.8085e+08  -1.5083e+07
   51.000m  -102.930  -1.1798e+09     eliminated
   52.000m  -102.930  -1.7752e+09   -2.588e+07
   53.000m  -102.930  -1.9421e+09  -2.8657e+07
   54.000m  -102.929  -2.3376e+09  -3.4762e+07
   55.000m  -102.930  -2.9566e+09     eliminated
   56.000m  -102.930  -3.3365e+09  -4.8734e+07
   57.000m  -102.930  -4.0548e+09  -5.9681e+07
   58.000m  -102.930  -4.3971e+09  -6.4855e+07
   59.000m  -102.929  -4.8903e+09  -7.2274e+07
   60.000m  -102.930  -4.8918e+09  -7.2235e+07
   61.000m  -102.930  -5.6283e+09  -8.2896e+07
   62.000m  -102.930   -6.001e+09  -8.9228e+07
   63.000m  -102.930  -6.0251e+09  -9.0299e+07
   64.000m  -102.930  -6.0083e+09  -9.0329e+07
   65.000m  -102.930  -5.9497e+09  -9.0009e+07
   66.000m  -102.930  -5.6575e+09  -8.6323e+07
   67.000m  -102.930  -5.1395e+09  -7.8959e+07
   68.000m  -102.930  -4.4553e+09  -6.8584e+07
   69.000m  -102.930  -3.8863e+09    -6.03e+07
   70.000m  -102.930  -3.2241e+09  -5.0751e+07
   71.000m  -102.930  -2.6172e+09  -4.1562e+07
   72.000m  -102.930  -1.8913e+09     eliminated
   73.000m  -102.929  -1.1893e+09  -1.9239e+07
   74.000m  -102.931  -8.3374e+08  -1.3694e+07
   75.000m  -102.930  -4.1716e+08     eliminated
   76.000m  -102.930  -1.2121e+08     eliminated
   77.000m  -102.930  -2.8492e+07     eliminated
   78.000m  -102.930  -1.8628e+08     eliminated
   79.000m  -102.930  -4.6888e+08  -3.0317e+06
   80.000m  -102.930   -4.909e+08     eliminated
   81.000m  -102.930  -7.4109e+08     eliminated
   82.000m  -102.930   -1.055e+09     eliminated
   83.000m  -102.930    -1.34e+09     eliminated
   84.000m  -102.930  -1.6526e+09   -1.813e+07
   85.000m  -102.930  -1.6702e+09   -1.949e+07
   86.000m  -102.929  -1.8575e+09  -2.1646e+07
   87.000m  -102.930  -1.9191e+09  -2.2642e+07
   88.000m  -102.930  -1.9324e+09  -2.3502e+07
   89.000m  -102.930  -1.8117e+09  -2.2251e+07
   90.000m  -102.930  -1.6157e+09  -1.9843e+07
   91.000m  -102.930  -1.5864e+09  -2.0016e+07
   92.000m  -102.930  -1.2679e+09  -1.5771e+07
   93.000m  -102.930  -1.1602e+09    -1.53e+07
   94.000m  -102.930  -8.8089e+08     eliminated
   95.000m  -102.930  -7.7341e+08  -9.9282e+06
   96.000m  -102.930   -5.845e+08  -7.7121e+06
   97.000m  -102.930  -2.8125e+08     eliminated
   98.000m  -102.930   -5.734e+07     eliminated
   99.000m  -102.930   8.8731e+07     eliminated
  100.000m  -102.930   2.2062e+08     eliminated
  101.000m  -102.930   3.9923e+08     eliminated
  102.000m  -102.930   6.4884e+08   1.0227e+07
  103.000m  -102.929   7.2869e+08   1.1137e+07
  104.000m  -102.930   7.1186e+08   1.1542e+07
  105.000m  -102.930   8.7938e+08     eliminated
  106.000m  -102.930   6.3424e+08   1.0064e+07
  107.000m  -102.930   5.6114e+08   9.2756e+06
  108.000m  -102.930   4.0482e+08   7.2406e+06
  109.000m  -102.930   2.6833e+08     eliminated
  110.000m  -102.930   1.1243e+08     eliminated
  111.000m  -102.930  -1.6522e+08     eliminated
  112.000m  -102.930  -3.2346e+08     eliminated
  113.000m  -102.931  -5.2643e+08  -7.8216e+06
  114.000m  -102.929   -8.369e+08  -1.2591e+07
  115.000m  -102.929  -9.1013e+08  -1.4483e+07
  116.000m  -102.930  -9.8441e+08  -1.5584e+07
  117.000m  -102.929  -1.0132e+09  -1.6514e+07
  118.000m  -102.929  -8.7345e+08    -1.46e+07
  119.000m  -102.930  -7.7383e+08  -1.3382e+07
  120.000m  -102.930  -5.0957e+08  -1.0336e+07
  121.000m  -102.930  -3.0658e+08  -6.8386e+06
  122.000m  -102.930  -1.0488e+08     eliminated
  123.000m  -102.930    4.176e+07     eliminated
  124.000m  -102.930   2.7125e+08     eliminated
  125.000m  -102.929   4.1905e+08    5.048e+06
  126.000m  -102.930   4.9524e+08     eliminated
  127.000m  -102.930     4.74e+08     eliminated
  128.000m  -102.929   5.5593e+08   8.1128e+06
  129.000m  -102.930   5.4432e+08   8.3087e+06
  130.000m  -102.930   4.6542e+08     eliminated
  131.000m  -102.930   3.3042e+08     eliminated
  132.000m  -102.930   1.9531e+08     eliminated
  133.000m  -102.930     4.41e+07     eliminated
  134.000m  -102.930  -5.2454e+07     eliminated
  135.000m  -102.930  -2.7328e+07     eliminated
  136.000m  -102.929  -8.7858e+07  -2.1463e+06
  137.000m  -102.930  -6.9246e+07     eliminated
  138.000m  -102.930  -1.3521e+08  -1.9875e+06
  139.000m  -102.929   -6.728e+07  -1.3155e+06
  140.000m  -102.930  -5.5859e+07     eliminated
  141.000m  -102.930   -1.977e+06     eliminated
  142.000m  -102.930   4.6275e+07     eliminated
  143.000m  -102.930   1.5744e+08   2.6183e+06
  144.000m  -102.930   1.7268e+08   2.9511e+06
  145.000m  -102.929   1.5363e+08   3.2818e+06
  146.000m  -102.930   1.1603e+08     eliminated
  147.000m  -102.930   1.2484e+08   3.3408e+06
  148.000m  -102.929   9.1881e+07   2.1693e+06
  149.000m  -102.930   1.2026e+07   1.5626e+06
  150.000m  -102.930   6.2462e+07     eliminated
  151.000m  -102.930    8.138e+06     eliminated
  152.000m  -102.930  -1.1122e+08     eliminated
  153.000m  -102.930  -6.5983e+07     eliminated


";

inadequatebuildups="Dataset :
 C:/Users/Mohamed/Downloads/20251024_19F2/866/pdata/1
AREA fit :
 I[t]=I[0]+P*exp(-t/T1)

50 points for Integral 1,  Integral Region from -102.847 to -103.018 ppm
Results     Comp. 1    

I[0]  =   -1.201e+04
P     =    1.201e+04
T1    =    1060.542s
SD    =    2.543e-01

    tau    ppm     integral    intensity

    3.000m  -102.931    3.804e+09     eliminated
    4.000m  -102.931   1.6923e+09     eliminated
    5.000m  -102.931  -2.8106e+09     eliminated
    6.000m  -102.931   5.5263e+09     eliminated
    7.000m  -102.931   9.3476e+09   1.1181e+08
    8.000m  -102.931    1.131e+10     eliminated
    9.000m  -102.931    1.197e+10     eliminated
   10.000m  -102.931   1.4245e+10     eliminated
   11.000m  -102.931   1.7544e+10   2.3473e+08
   12.000m  -102.931   1.8599e+10   2.4968e+08
   13.000m  -102.931   2.1873e+10   2.8364e+08
   14.000m  -102.931   2.2204e+10     eliminated
   15.000m  -102.931   2.2674e+10     eliminated
   16.000m  -102.931   2.3364e+10     eliminated
   17.000m  -102.931   2.7196e+10   3.6389e+08
   18.000m  -102.931   2.4462e+10   3.3324e+08
   19.000m  -102.931    2.228e+10    3.049e+08
   20.000m  -102.931   2.3431e+10   3.2657e+08
   21.000m  -102.931   2.1485e+10   3.0685e+08
   22.000m  -102.931   2.1178e+10   3.0791e+08
   23.000m  -102.931   1.9396e+10     eliminated
   24.000m  -102.931    1.803e+10   2.6327e+08
   25.000m  -102.931   1.8388e+10     eliminated
   26.000m  -102.931   1.6585e+10   2.4349e+08
   27.000m  -102.931   1.5564e+10   2.3842e+08
   28.000m  -102.931   1.4725e+10   2.2017e+08
   29.000m  -102.931   1.3096e+10   2.0237e+08
   30.000m  -102.931   1.1275e+10   1.7668e+08
   31.000m  -102.931   1.1792e+10   1.7696e+08
   32.000m  -102.931   9.8814e+09   1.5287e+08
   33.000m  -102.931   9.1861e+09   1.4297e+08
   34.000m  -102.931   8.5482e+09   1.3534e+08
   35.000m  -102.931   7.5793e+09   1.2054e+08
   36.000m  -102.931   6.7855e+09   1.0937e+08
   37.000m  -102.931   4.5003e+09   8.0739e+07
   38.000m  -102.931   5.2372e+09   8.7348e+07
   39.000m  -102.932   4.5377e+09   7.8899e+07
   40.000m  -102.931   4.1864e+09   6.5015e+07
   41.000m  -102.931   3.6797e+09   6.5741e+07
   42.000m  -102.931   3.1479e+09   5.6422e+07
   43.000m  -102.931   2.8733e+09     eliminated
   44.000m  -102.931   2.6714e+09   4.7126e+07
   45.000m  -102.932   2.0444e+09   3.9082e+07
   46.000m  -102.931    1.942e+09     eliminated
   47.000m  -102.930   1.6444e+09   3.3317e+07
   48.000m  -102.931     1.32e+09   2.6385e+07
   49.000m  -102.931   1.5791e+09     eliminated
   50.000m  -102.931   1.0183e+09     eliminated
   51.000m  -102.932   1.0046e+09   2.5529e+07
   52.000m  -102.932   8.3945e+08   1.8906e+07";

TopspinIntegralImporter[m_]:= ToExpression/@ StringReplace[Select[(Drop[StringSplit[m],First@@Position[StringSplit[m,Whitespace],"intensity"]]//Partition[#,4]&),Length[#]==4 &][[All,3]],"e"->"*^"]

IntegralImporterNormalized[m_]:=TopspinIntegralImporter[m]/Max@TopspinIntegralImporter[m]

IntegralImporterNormalized[m_,normalizationfactor_]:=IntegralImporterNormalized[m] * normalizationfactor

PPolData=Transpose[{Table[n 0.002935 ,{n,0,180,2}],IntegralImporterNormalized[pulsepolbuildups,0.5647]}];

SLICData = Transpose[{Range[0,0.5,0.5/(128-1)],IntegralImporterNormalized[slicbuildups,0.1496]}];

INADEQUATEData = Transpose[{Range[0.01,0.99,0.02],IntegralImporterNormalized[inadequatebuildups,0.0543]}];

CSLICdata=Transpose[{Table[0.0039156n ,{n,0,151,1}],IntegralImporterNormalized[cslicbuildups,0.5621]}];

inadsims=Table[{t,TransformationAmplitude[opI["x"]->opI["x"],{{None,t/2},Pulse[{π,0}],{None,t/2},Pulse[{π/2,π/2}],CoherenceOrderFiltrationSuperoperator[{-2,2}],Pulse[{π/2,π/2}],{None,t/2},Pulse[{π,0}],{None,t/2}},BackgroundGenerator->2π 255.94 opI[1].opI[2]+2π 17.8SingleTransitionOperator[{1,3},"x"]]},{t,0,1,0.05}];

ppolsims=Table[{nR 0.002935,TransformationAmplitude[opI["x"]->opI["x"],{PulsePolWithBB1[0.002935,nR,2π 25000],Pulse[{π/2,0}],CoherenceOrderFiltrationSuperoperator[{-2,2}],Pulse[{π/2,π}],PulsePol5[0.002935,nR,2π 25000]},BackgroundGenerator->(2π 255.94 opI[1].opI[2]+2π 17.8 SingleTransitionOperator[{1,3},"x"])]},{nR,0,168,4}];

slicsims=TransformationAmplitudeTable[opI["z"]->-opI["z"],{SLIC[ 2π 255.94,t],Pulse[{π/2,0}],CoherenceOrderFiltrationSuperoperator[{-2,2}],Pulse[{π/2,0}],SLIC[ 2π 255.94,t]},{t,0,128/255.94,1/(255.94)},BackgroundGenerator->2π 255.94 opI[1].opI[2]+2π 17.8 SingleTransitionOperator[{1,3},"x"],Parallel->False];

{d10,d11,drest}={0.003868/2,0.000039596/2,2*(10^(-6))};

cSLICexp={{None,drest},{2π 255.94 opI["x"],d10},{None,drest},{-2π 25000 opI["x"],d11},{None,drest},{-2π 25000 opI["x"],d11},{None,drest},{2π 255.94 opI["x"],d10}};

cSLICsims=Table[{n 0.0039156,TransformationAmplitude[opI["z"]->opI["z"],{Repeat[cSLICexp,n],Pulse[{π/2,0}],CoherenceOrderFiltrationSuperoperator[{2,-2}],Pulse[{π/2,0}],Repeat[cSLICexp,n]},BackgroundGenerator->2π 255.94 opI[1].opI[2]+2π 17.8 SingleTransitionOperator[{1,3},"x"]]},{n,0,128,2}];

func[T_]=(((Cos[θST] Sin[Ω τ1]Cos[ωJ τ1]-(((Sin[θST])^2)+((Cos[θST])^2) Cos[Ω τ1])Sin[ωJ τ1]))^2) /.Ω->Sqrt[((ωΔ)^2)+((ωJ)^2)]/. θST->ArcTan[ωΔ/ωJ]/.ωΔ ->2π 17.8/.ωJ->2π 255.94/.τ1->T/2

{SLICmax,PPolmax,CSLICmax,INADmax}=Flatten[Take[#,Ordering[Last/@#,-1]]&/@{SLICData,PPolData,CSLICdata,INADEQUATEData},1]

Flatten[Take[#,Ordering[Last/@#,1]]&/@{CSLICdata},1]

Clear[custmplotmark]
custmplotmark[name_String,Fcol_,Opact_:White,size_:6]:=Show[ResourceFunction["PolygonMarker"][name,Offset@size,{FaceForm[Opact],Dynamic@EdgeForm[{CurrentValue["Color"],AbsoluteThickness[1.618],Opacity[1]}]}],ResourceFunction["PolygonMarker"][name,Offset[Fcol size],EdgeForm[None]]];

gridlineindicator[max_]:={Directive[Dashed,Thickness[0.004]],Line[{{First@max,-1},max}]}

axes={Directive[Black,Thickness[0.0015]],Line[{{0,0},{0.5,0}}]}

PlotMarkers->{custmplotmark["Circle",0,Black,6]}

ListPlot[INADEQUATEData,PlotStyle->Directive[Black,Thickness[0.006]],Joined->True,PlotRange->{{0,0.5},{-0.115,1.1}},Epilog->{gridlineindicator[INADmax],axes}]

ResourceFunction["PlotGrid"][{{Show[ListPlot[INADEQUATEData,PlotStyle->Directive[Black,Thickness[0.006]],Joined->True,PlotRange->{{0,0.5},{-0.12,1.1}},Epilog->{gridlineindicator[INADmax],axes},MaxPlotPoints->200],Plot[func[T] ,{T,0,0.5},PlotPoints->40,PlotStyle->{Directive[Blue,Thickness[0.006]]}],ListPlot[inadsims,Joined->False,PlotStyle->{Directive[Black,Thickness[0.02]]},PlotMarkers->{custmplotmark["Circle",0,White,6]},PlotRange->{{0,0.5},{-0.115,1.1}}],FrameLabel->{{"Normalized integral",None},{"Time (s)",None}},FrameTicks->{{{{0,"0.0"},{0.2,0.2},{0.4,0.4},{0.6,0.6},{0.8,0.8},{1.0,"1.0"}},Automatic},{{{0,"0.0"},{0.2,0.2},{0.4,0.4}},Automatic}}]},{Show[ListPlot[PPolData,Axes->True,Joined->True,PlotStyle->{Directive[Black,Thickness[0.006]]},AxesStyle->Black,Epilog->{gridlineindicator[PPolmax],axes},MaxPlotPoints->200,PlotRange->{{0,0.5},{-0.12,1.1}},FrameLabel->{{"Normalized integral",None},{"Time (s)",None}},FrameTicks->{{{{0,"0.0"},{0.2,0.2},{0.4,0.4},{0.6,0.6},{0.8,0.8},{1.0,"1.0"}},True},{{{0,"0.0"},{0.2,0.2},{0.4,0.4}},Automatic}}],Plot[((Sin[(1/4)ωΔ κ t])^4)/.ωΔ ->2π 17.8/.κ->κPulsePol,{t,0,0.5},PlotStyle->{Directive[Blue,Thickness[0.006]]}],ListPlot[Abs@ppolsims,Joined->False,PlotStyle->{Directive[Black,Thickness[0.02]]},PlotMarkers->{custmplotmark["Circle",0,White,6]},PlotRange->{{0,0.5},All}]]},{Show[ListPlot[SLICData,Axes->True,Joined->True,AxesStyle->Black,PlotStyle->{Directive[Black,Thickness[0.006]]},Epilog->{gridlineindicator[SLICmax],axes},PlotRange->{{0,0.5},{-0.12,1.1}}],Plot[((Sin[(1/4)ωΔ κ t])^4)/.ωΔ ->2π 17.8/.κ->1/Sqrt[2],{t,0,0.5},PlotStyle->{Directive[Blue,Thickness[0.006]]}],ListPlot[Abs@slicsims,Joined->False,PlotMarkers->{custmplotmark["Circle",0,White,6]},PlotRange->{{0,0.5},All}],FrameLabel->{{"Normalized integral",None},{"Time (s)",None}},FrameTicks->{{{{0,"0.0"},{0.2,0.2},{0.4,0.4},{0.6,0.6},{0.8,0.8},{1.0,"1.0"}},True},{{{0,"0.0"},{0.2,0.2},{0.4,0.4}},Automatic}}]},{Show[ListPlot[CSLICdata,Axes->True,Joined->True,AxesStyle->Black,PlotStyle->{Directive[Black,Thickness[0.006]]},PlotRange->{{0,0.5},{-0.12,1.1}},Epilog->{gridlineindicator[CSLICmax],axes},FrameLabel->{{"Normalized integral",None},{"T / s",None}},FrameTicks->{{{{0,"0.0"},{0.2,0.2},{0.4,0.4},{0.6,0.6},{0.8,0.8},{1.0,"1.0"}},True},{Automatic,Automatic}}],Plot[((Sin[(1/4)ωΔ κ t])^4)/.ωΔ ->2π 17.8/.κ->1/Sqrt[2],{t,0,0.5},PlotStyle->{Directive[Blue,Thickness[0.006]]}],ListPlot[Abs@cSLICsims,Joined->False,PlotMarkers->{custmplotmark["Circle",0,White,6]},PlotRange->{{0,0.5},{-0.12,1.1}}]]}},Spacings->{0,0},ImageSize->600,AspectRatio->2.5,"ShowFrameLabels"->Full,PlotRangePadding->False,PlotLabels->Placed[Row->"a",{0.055,0.926}],LabelStyle->{Black,35,FontFamily->"Helvetica"}]

sims for rf errors (Figure 8)

nocompsims=Table[{ϵ,TransformationAmplitude[opI["z"]->-2(opI[1,"x"].opI[2,"y"]+opI[1,"y"].opI[2,"x"]),{{2π 255(1+ϵ) opI["x"],20/255.94 +1/(4 255.94)}},BackgroundGenerator->H]},{ϵ,-3,1,0.0125}];

Solve[k (1/(255.94)+1/(25000)) == 1/(255.94),k]//N

supercompsims=Table[{ϵRF,TransformationAmplitude[opI["z"]->-2opI[1,"x"].opI[2,"y"]-2opI[1,"y"].opI[2,"x"],{Repeat[{{(1+ϵRF) 2π 255.94 opI["x"],0.9898661463402273`/(2 255.94)},{-(1+ϵRF) 2π 25000 opI["x"],0.9898661463402273`/( 25000)},{(1+ϵRF) 2π 255.94 opI["x"],0.9898661463402273`/(2 255.94)}},20],Pulse[{π/2,0}]},BackgroundGenerator->H]},{ϵRF,-1,1,0.1}];

RFerrorExpression["supercompensation",ωΔ_,{ϵ_,t_}]=((4(1+ϵ) ((Sin[((Sqrt[2])/(8π))ωΔ ((Sin[π ϵ])/ϵ)t (Sqrt[1+((ϵ^2)/(((2+ϵ))^2))])])^2))/((ϵ^2)+(((2+ϵ))^2)));

RFerrorExpression["nocompensation",κ_,{ωJ_,ωΔ_},{ϵrf_,t_}]=(1/2) (Sin[t (1+ϵrf) ωJ]-Cos[(1/2) t(Sqrt[ ((ϵrf)^2) ((ωJ)^2)+(κ^2) ((ωΔ)^2)])] Sin[(1/2) t (2+ϵrf) ωJ]-((ϵrf ωJ Cos[(1/2) t (2+ϵrf) ωJ] Sin[(1/2) t(Sqrt[ ((ϵrf)^2) ((ωJ)^2)+(κ^2) ((ωΔ)^2)])])/(Sqrt[ ((ϵrf)^2) ((ωJ)^2)+(κ^2) ((ωΔ)^2)])));

Show[{Plot[{RFerrorExpression["supercompensation",ωΔ,{ϵ,τseq}]/.ωΔ->2π 17.8/.τseq->20/255.94,RFerrorExpression["nocompensation",κΔ,{ωJ,ωΔ},{ϵ,t}]/.κΔ->1/Sqrt[2]/.ωΔ->2π 17.8/.ωJ->2π 255.94/.t->20/255.94+1/(4 255.94)},{ϵ,-1,1},PlotLegends->{,}],ListPlot[{supercompsims,nocompsims},Joined->False,PlotMarkers->{"OpenMarkers",8}],PlotRange->{{-1,1},All}},FrameLabel->{,}]

Resonance offset errors (Fig. 9)

nocompdetuningsims=Table[{ϵ,TransformationAmplitude[opI["z"]->-2(opI[1,"x"].opI[2,"y"]+opI[1,"y"].opI[2,"x"]),{{2π 255.94 opI["x"],20/255.94+1/(4 255.94)}},BackgroundGenerator->(H+2π J ϵ opI["z"])]},{ϵ,-1,1,0.005}];

supercompdetuningsims=Table[{ϵ,TransformationAmplitude[opI["z"]->-2opI[1,"x"].opI[2,"y"]-2opI[1,"y"].opI[2,"x"],{Repeat[{{(1) 2π 255.94 opI["x"],0.9898661463402273`/(2 255.94)},{-(1) 2π 25000 opI["x"],0.9898661463402273`/( 25000)},{(1) 2π 255.94 opI["x"],0.9898661463402273`/(2 255.94)}},20],Pulse[{π/2,0}]},BackgroundGenerator->(H+2π J ϵ opI["z"])]},{ϵ,-1,1,0.005}];

MaTeX["\\Delta \\omega /\\omega_{J}"]

ListPlot[{supercompdetuningsims,nocompdetuningsims},Joined->True,PlotStyle->{Directive[Black,Thickness[0.005]],Directive[ColorBlindSafeColorScheme[[2]],Thickness[0.005]]},FrameLabel->{"",""},PlotLegends->{,}]