Dear all,
I've been trying to simulate an inversion recovery experiment, but I'm stuck at the very beginning. The pulse sequence does not produce an inverted echo and I'm not sure why, can you please help me?
Here is the script I'm using:
clear, clc
Sys.g = 2; %isotropic g tensor
Sys.Nucs = '1H';
Sys.A = [10 40]*0.5; %MHz
Sys.lwpp = 5; %mT
ExpSpectrum.Field = 324.9; %mT
pepper(Sys,ExpSpectrum);
% inversion
p180.tp = 0.04; %ms
p180.Frequency = 0;
p180.Flip = pi;
p90=p180;
p90.Flip = pi/2;
p90.tp = p180.tp/2;
tau = 0.2;
T = 1;
Exp.Sequence = {p180 T p90 tau p180 tau+p90.tp/2};
Exp.Field = 324.9; %mT
Exp.mwFreq = 9.098;
Exp.DetWindow = [-0.08 +0.08];
Opt.nKnots = 10;
[t,signal] = saffron (Sys,Exp,Opt);
plot (t,signal);
Thank you again,
Edoardo.