Area under Gaussian peak by easy measurements

[fysikjatack]

First I try to find a way of finding area under Gaussian peak by using simple means. Questions at bottom of post.

This concerns peaks with Gaussian distribution, eg. peaks in a chromatogram in some chemical analysis.

The Gaussian distribution is
<br /> f(x) = \frac{1}{\sigma \sqrt{2\pi} } e^{ -\frac{(x-\mu)^2}{2\sigma ^2} }<br />

and the function for a peak in the chromatogram is

<br /> f(x) = \frac{A}{\sigma \sqrt{2\pi} } e^{ -\frac{(x-\mu)^2}{2\sigma ^2} }<br />

where A is the area under the curve. The area is interesting since it represents the number of detected molecules and by using calibration curves one can find the concentration of a previously unknown sample.
The maximum peak height H is easily proven to be found at x =\mu, hence the exponential part is equal to 1;

<br /> H=f(\mu) = \frac{A}{\sigma \sqrt{2\pi} } e^{ -\frac{(\mu-\mu)^2}{2\sigma ^2} }=\frac{A}{\sigma \sqrt{2\pi} } = \frac{0.3989A}{\sigma}<br />

This shows that the area (and hence concentration) is proportional to maximum peak height. But this assumes same /sigma. But I would like to consider that injection and stuff might be different etc and account for peaks having a bit different /sigma.

The full width at half maximum (FWHM) can be shown to be
FWHM=2 sqrt(2 ln(2))\sigma = 2.35\sigma by inserting f(x) = H/2 , find x1 and x2 and then calc the width.
<br /> H = \frac{0.3989A}{\sigma}<br />

<br /> A = \frac{H\sigma}{0.3989}<br />
<br /> FWHM = 2.35\sigma<br />
<br /> \sigma = \frac{FWHM}{2.35}<br />

<br /> A = \frac{H* FWHM}{2.35*0.3989}<br />


Then the area would easily be calculated from a paper with peaks using only ruler and calculator, by measuring the maximum height and FWHM, multiplying them and multiply by a constant.


1.Did I overlook something here, or does it work this way?

It seems very convenient to account for peaks having different /sigma this way, although I doubt it has any large practical significance in the high school lab... Also, this still is under assumption that all peaks are perfectly gaussian in shape.

In school we usually have software to calculate peak areas, or just plot peak maximum height.

2. How does a software compute such an area?


The integral with respect to x of <br /> f(x) = \frac{A}{\sigma \sqrt{2\pi} } e^{ -\frac{(x-\mu)^2}{2\sigma ^2} }<br /> would equal A. Integration yields A=A and does not give the value of the area...

3. If I just sit here with a bunch of peaks on a paper, can I use some method applying integrals to find A?

 

[canucks357]

I just spent some time trying to figure this out. For an equation of the form:

f(x)=N\cdot exp[-\frac{(x-u)^2}{2\sigma^2}]

The area under the curve is given by:

N\cdot\sigma\cdot\sqrt{2\pi}