EPR simulation with two spins (radical pair)
Posted: Thu Feb 25, 2016 6:36 am
Dear all,
I'm trying to simulate the EPR spectrum for a radical pair system using garlic. I have tried to match the simulated spectrum as close to the experimental data but the result still show deviations from the experimental one. How can I improve the simulation furthermore, so that it exactly matches the experimental spectrum? The code is as below and the experimental data in ascii as well as raw data (JEOL format) is attached with the query.
Please help.
Thanks & regards,
Abbey
I'm trying to simulate the EPR spectrum for a radical pair system using garlic. I have tried to match the simulated spectrum as close to the experimental data but the result still show deviations from the experimental one. How can I improve the simulation furthermore, so that it exactly matches the experimental spectrum? The code is as below and the experimental data in ascii as well as raw data (JEOL format) is attached with the query.
Code: Select all
clear
SysA.S = 1/2;
SysA.g = 2.012;
SysA.Nucs = '1H,1H,1H,1H,1H';
SysA.A = mt2mhz([15.95,15.95,14.35,14.35,14.1]/10);
SysA.lwpp = 1.33;
SysB.S = 1/2;
SysB.g = 2.012;
SysB.Nucs = '1H,1H,1H,1H,1H';
SysB.A = mt2mhz([15.95,15.95,14.35,14.35,14.1]/10);
SysB.lwpp = 1.33;
%Experimental details
Exp.mwFreq = 9.137;
Exp.Range = [310.569 339.9945];
Exp.nPoints = 64293;
garlic({SysA SysB},Exp);
hold on;
%Loading experimental spectrum.
data=textread('PA_exp_norm B.dat');
x=data(:,1);
y=data(:,2);
plot(x,y,'k');
hold off;
Thanks & regards,
Abbey