Spectrum convolution in pepper

A place to report and discuss potential bugs
Post Reply
ovlaguta
Newbie
Posts: 1
Joined: Tue Jun 04, 2024 11:47 pm

Spectrum convolution in pepper

Post by ovlaguta »

Hi all,
It seems there's a convolution problem (or limitation) in pepper (EasySpin version: 6.0.0 (2024-04-01)). Basically, it sets intensity to zero for points below 6x10^-4 relative intensity. The convspec function though works properly. Please see a simple example below.

Sys.S = 1/2;
Sys.g = 2;
Sys.lw = [0 1];

Exp.mwFreq = 9.8;
Exp.Range = [300 400];
Exp.SampleFrame = [0 0 0];
Exp.MolFrame = [0 0 0];
Exp.Harmonic = 0;
Exp.nPoints = 1000;

[m,spc] = pepper(Sys,Exp);
spc = spc/max(spc);
spec = spc*0;
[~, i] = max(spc);
spec(i) = 1;
w = convspec(spec,0.1,1,0,0);
w = w/max(w);
plot(m,spc, m,w)
ylim([0 1e-3])
xlabel('Magnetic field (mT)')
legend('pepper', 'convspec')

Stefan Stoll
EasySpin Creator
Posts: 1106
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: Spectrum convolution in pepper

Post by Stefan Stoll »

This is in purpose. The Lorentzian lineshape is truncated in pepper for perfomance reasons. Is this truncation visible in any practically relevant situation?

Stefan Stoll
EasySpin Creator
Posts: 1106
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: Spectrum convolution in pepper

Post by Stefan Stoll »

It should be fixed in 6.0.1, which is available now.

Post Reply