Subject: Re: [nmr_sparky] exponential decay
From: Mandar T. Naik
Date: Jun 2, 2010

Previous: 722 Next: 725


Hi,
My two cents on an alternate work around. Get the sparky2rate script from
Loria groups website @Yale. Its a pearl script which takes Sparky output and
feeds it to Palmer labs (@Columbia) Curvefit program. Curvefit supports both
functions you want to try. You will need Grace to visualize the output. The
error estimates are calculated differently by Sparky and Curvefit.
The whole exercise is worth but you may also use any commercial software like
Prism or Sigmaplot for the same purpose. I havent dealt with data for long
time so doesnt know any new easier way.
Best regards
-mandar
On 2010 June 3 Thursday 03:43:03 Stefano Ciurli wrote:
Hi Tom,
thanks a lot. I would know how to use eq. b with another fitting
program but I will try to implement it in the python script within
Sparky. However, my question was also a bit more general: I wanted to
know how the community of NMR people deals with this problem, as it
is not clear from the literature and I am dealing with this problem
for the first time.
I hope someone in the list could also respond.
Stefano

On Jun 2, 2010, at 9:29 PM, Tom Goddard wrote:
Hi Stefano,

A modification to the Sparky relax.py file (sparky/python/sparky/
relax.py) could add the parameter C to the fitting. I think I did
this once upon a time but I cant find any code for it. Youd
modify the a_exp_bx_function class in the relax.py code. Maybe an
hours work to figure out change and test if you are familiar with
Python.

Tom

Hello,
when fitting T1 or T2 relaxation data one could use two equations:

a: A*exp(-t/B)

b: A*exp(-t/B) + C

In the first case one assumes that the maximum tau used in the set
of experiments is sufficient to bring to zero the intensity of the
signal, while in the second case one assumes that there may be a
residual intensity. Following my logic, I would prefer the second
case, because the spectrum with the largest tau still has non-zero
signals. However, Sparky uses the first equation. Of course, the
relaxation rates that are calculated in the two cases can be very
different. What is the commen/preferred protocol?
Stefano

--
LA RICERCA C’È E SI VEDE:
5 per mille allUniversità di Bologna - C.F.: 80007010376
http://www.unibo.it/Vademecum5permille.htm

Questa informativa è inserita in automatico dal sistema al fine esclusivo
della realizzazione dei fini istituzionali dell’ente.