Curve decomposition Some spectroscopic curves can be approximated by the sum of a set of component curves. For example, when
Beer's law :A_\lambda=\sum_k \epsilon_{k,\lambda} c_k applies, the total absorbance,
A, at wavelength λ, is a
linear combination of the absorbance due to the individual components,
k, at
concentration,
ck. ε is an
extinction coefficient. In such cases the curve of experimental data may be decomposed into sum of component curves in a process of
curve fitting. This process is also widely called deconvolution.
Curve deconvolution and
curve fitting are completely different mathematical procedures. • Parameters of the line shape are unknown. The intensity of each component is a function of at least 3 parameters, position, height and half-width. In addition one or both of the line shape function and baseline function may not be known with certainty. When two or more parameters of a fitting curve are not known the method of non-linear
least squares must be used. The reliability of curve fitting in this case is dependent on the separation between the components, their shape functions and relative heights, and the
signal-to-noise ratio in the data. When Gaussian-shaped curves are used for the decomposition of set of
Nsol spectra into
Npks curves, the p_0 and w parameters are common to all
Nsol spectra. This allows to calculated the heights of each Gaussian curve in each spectrum (
Nsol·
Npks parameters) by a (fast) linear least squares fitting procedure, while the p_0 and
w parameters (2·
Npks parameters) can be obtained with a non-linear least-square fitting on the data from all spectra simultaneously, thus reducing dramatically the correlation between optimized parameters.
Derivative spectroscopy Spectroscopic curves can be subjected to
numerical differentiation. When the data points in a curve are equidistant from each other the
Savitzky–Golay convolution method may be used. The best convolution function to use depends primarily on the signal-to-noise ratio of the data. The first derivative (slope, \frac{dy}{dx}) of all single line shapes is zero at the position of maximum height. This is also true of the third derivative; odd derivatives can be used to locate the position of a peak maximum. The second derivatives, \frac{d^2y}{dx^2}, of both Gaussian and Lorentzian functions have a reduced half-width. This can be used to apparently improve
spectral resolution. The diagram shows the second derivative of the black curve in the diagram above it. Whereas the smaller component produces a shoulder in the spectrum, it appears as a separate peak in the 2nd. derivative. Fourth derivatives, \frac{d^4y}{dx^4}, can also be used, when the signal-to-noise ratio in the spectrum is sufficiently high.
Deconvolution Deconvolution can be used to apparently improve
spectral resolution. In the case of NMR spectra, the process is relatively straight forward, because the line shapes are Lorentzian, and the convolution of a Lorentzian with another Lorentzian is also Lorentzian. The
Fourier transform of a Lorentzian is an exponential. In the co-domain (time) of the spectroscopic domain (frequency) convolution becomes multiplication. Therefore, a convolution of the sum of two Lorentzians becomes a multiplication of two exponentials in the co-domain. Since, in FT-NMR, the measurements are made in the time domain division of the data by an exponential is equivalent to deconvolution in the frequency domain. A suitable choice of exponential results in a reduction of the half-width of a line in the frequency domain. This technique has been rendered all but obsolete by advances in NMR technology. A similar process has been applied for resolution enhancement of other types of spectra, with the disadvantage that the spectrum must be first Fourier transformed and then transformed back after the deconvoluting function has been applied in the spectrum's co-domain. == See also ==