Two-dimensional correlation analysis
Two dimensional correlation analysis is a mathematical technique that is used to study changes in measured signals. As mostly spectroscopic signals are discussed, sometime also two dimensional correlation spectroscopy is used and refers to the same technique.
In 2D correlation analysis, a sample is subjected to an external perturbation while all other parameters of the system are kept at the same value. This perturbation can be a systematic and controlled change in temperature, pressure, pH, chemical composition of the system, or even time after a catalyst was added to a chemical mixture. As a result of the controlled change (the perturbation), the system will undergo variations which are measured by an chemical or physical detection method. The measured signals or spectra will shown systematic variations that are processed with 2D correlation analysis for interpretation.
When one considers spectra that consist of few bands, it is quite obvious to determine which bands are subject to a changing intensity. Such a changing intensity can be caused for example by chemical reactions. However, the interpretation of the measured signal becomes more tricky when spectra are complex and bands are heavily overlapping. Two dimensional correlation analysis allows one to determine at which positions in such a measured signal there is a systematic change in a peak, either continuous rising or drop in intensity. 2D correlation analysis results in two complementary signals, which referred to as the 2D synchronous and 2D asynchronous spectrum. These signals allow amongst others[1][2][3]
- to determine the events that are occurring at the same time (in phase) and those events that are occurring at different times (out of phase)
- to determine the sequence of spectral changes
- to identify various inter- and intramolecular interactions
- band assignments of reacting groups
- to detect correlations between spectra of different techniques, for example near infrared spectroscopy (NIR) and Raman spectroscopy
History
2D correlation analysis originated from 2D NMR spectroscopy. Isao Noda developed perturbation based 2D spectroscopy in the 1980's.[4] This technique required sinusoidal perturbations to the chemical system under investigation. This specific type of the applied perturbation severely limited its possible applications. Following research from several groups, perturbation based 2D spectroscopy was further extended and generalized. Since the development of generalised 2D correlation analysis in 1993 based on fourier transformation of the data, 2D correlation analysis gained widespread use. Alternative techniques that were simpler to calculate, for example the disrelation spectrum, were also developed. Because of its computational efficiency and simplicity, the Hilbert transform is nowadays used for the calculation of the 2D spectra. To date, 2D correlation analysis is not only used for the interpretation of many types of spectroscopic data (including XRF, UV/VIS spectroscopy, fluorescence, infrared, and Raman spectra), although its application is not limited to spectroscopy.
Properties of 2D correlation analysis
2D correlation analysis is frequently used for its main advantage: increasing the spectral resolution by spreading overlapping peaks over two dimensions and as a result simplification of the interpretation of one dimensional spectra that are otherwise visually indistinguishable from each other.[4] Further advantages are its ease of application and the possibility to make the distinction between band shifts and band overlap.[3] Each type of spectral event, band shifting, overlapping bands of which the intensity changes in the opposite direction, band broadening, baseline change, etc. has a particular 2D pattern.
Presence of 2D spectra
2D synchronous and asynchronous spectra are basically 3D-datasets and are generally represented by contour plots. X- and y-axes are identical to the x-axis of the original dataset, whereas the different contours represent the magnitude of correlation between the spectral intensities. The 2D synchronous spectrum is symmetric relative to the main diagonal. The main diagonal thus contains positive peaks. As the peaks at (x,y) in de 2D synchronous spectrum are a measure for the correlation between the intensity changes at x and y in the original data, these main diagonal peaks are also called autopeaks and the main diagonal signal is referred to as autocorrelation signal. The off-diagonal cross-peaks can be either positive or negative. On the other hand the asynchronous spectrum is asymmetric and never has peaks on the main diagonal.
Generally contour plots of 2D spectra are oriented with rising axes from left to right and top to down. Other orientations are possible, but interpretation has to be adapted accordingly.[5]
Calculation of 2D spectra
Suppose the original dataset D contains the n spectra in rows. A 2D synchronous spectrum expresses the similarity between spectral of the data in the original dataset. In generalised 2D correlation spectroscopy this is mathematically expressed as covariance (or correlation).
where:
- Ψ is the 2D synchronous spectrum
- ν1 en ν1 are two spectral channels
- yν is the vector composed of the signal intensities in D in column ν
- n the number of signals in the original dataset
Interpretation
Interpretation of two-dimensional correlation spectra can be considered to consist of several stages.[4]
Detection of peaks of which the intensity changes in the original dataset
As real measurement signals contain a certain level of noise, the derived 2D spectra are influenced and degraded with substantial higher amounts of noise. Hence, interpretation begins with studying the autocorrelation spectrum on the main diagonal of the 2D synchronous spectrum. The intensity of peaks on the autocorrelation spectrum are directly proportional to the relative importance of the intensity change in the original spectra. Hence, if an intense band is present at position x, it is very likely that a true intensity change is occurring and the peak is not due to noise.
Additional techniques help to filter the peaks that can be seen in the 2D synchronous and asynchronous spectra.[6]
Determining the direction of intensity change
It is not always possible to unequivocally determine the direction of intensity change, such as is for example the case for highly overlapping signals next to each other and of which the intensity changes in the opposite direction. This is where the off diagonal peaks in the synchronous 2D spectrum are used for:
- if there is a positive cross-peak at (x, y) in the synchronous 2D spectrum, the intensity of the signals at x and y changes in the same direction
- if there is a negative cross-peak at (x, y) in the synchronous 2D spectrum, the intensity of the signals at x and y changes in the opposite direction
Determining the sequence of events
Most importantly, with the sequential order rules, also referred to as Noda's rules, the sequence of the intensity changes can be determined.[4] By carefully interpreting the signs of the 2D synchronous and asynchronous cross peaks with the following rules, the sequence of spectral events during the experiment can be determined:
- if the intensities of the bands at x and y in the dataset are changing in the same direction, the synchronous 2D cross peak at (x,y) is positive
- if the intensities of the bands at x and y in the dataset are changing in the opposite direction, the synchronous 2D cross peak at (x,y) is negative
- if the change at x mainly precedes the change in the band at y, the asynchronous 2D cross peak at (x,y) is positive
- if the change at x mainly follows the change in the band at y, the asynchronous 2D cross peak at (x,y) is negative
- if the asynchronous 2D cross peak at (x,y) is negative, the interpretation of rule 1 and 2 for the synchronous 2D peak at (x,y) has to be reversed
- where x and y are the positions on the x-xaxis of two bands in the original data that are subject to intensity changes.
See also
- Correlation spectroscopy
- Fluorescence correlation spectroscopy
- Two-dimensional hetero-correlation spectroscopy
- Two-dimensional infrared spectroscopy
References
- ^ Shin-Ichi Morita, Yasuhiro F. Miura, Michio Sugi, and Yukihiro Ozaki (2005). "New correlation indices invariant to band shifts in generalized two-dimensional correlation infrared spectroscopy". Chemical Physics Letters. 402 (251–257): 251–257.
{{cite journal}}
: CS1 maint: multiple names: authors list (link) - ^ Koichi Murayama, Boguslawa Czarnik-Matusewicz, Yuqing Wu, Roumiana Tsenkova, and Yukihiro Ozaki (2000). "Comparison between conventional spectral analysis methods, chemometrics, and two-dimensional correlation spectroscopy in the analysis of near-infrared spectra of protein". Applied Spectroscopy. 54 (7): 978–985.
{{cite journal}}
: CS1 maint: multiple names: authors list (link) - ^ a b Shin-Ichi Morita and Yukihiro Ozaki (2002). "Pattern recognitions of band shifting, overlapping, and broadening using global phase description derived from generalized two-dimensional correlation spectroscopy". Applied Spectroscopy. 56 (4): 502–508.
- ^ a b c d Isao Noda and Yukihiro Ozaki (2004). Two-Dimensional Correlation Spectroscopy - Applications in Vibrational and Optical Spectroscopy. Chichester: John Wiley & Sons Ltd.
- ^ Boguslawa Czarnik-Matusewitcz, Sylwia ilorz, Lorna Ashton, and Ewan W. (2006). "Potential pitfalls concerning visualisation of the 2D results". Journal of molecular structure. 799 (1–3): 253–258. doi:10.1016/j.molstruc.2006.03.064.
{{cite journal}}
: CS1 maint: multiple names: authors list (link) - ^ R. Buchet, Y. Wu, and G. Lachenal, C. Raimbault, and Yukihiro Ozaki (2006). "Selecting two-dimensional cross-correlation functions to enhance interpretation of near-infrared spectra of proteins". Applied Spectroscopy. 55 (2): 155–162.
{{cite journal}}
: CS1 maint: multiple names: authors list (link)