Two-dimensional correlation analysis

Last updated

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.

Contents

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 a chemical or physical detection method. The measured signals or spectra will show 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]

  1. 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)
  2. to determine the sequence of spectral changes
  3. to identify various inter- and intramolecular interactions
  4. band assignments of reacting groups
  5. 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 1980s. [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 done by several groups of scientists, perturbation based 2D spectroscopy could be developed to a more extended and generalized broader base. Since the development of generalized 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 simultaneously. 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 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

Demo dataset consisting of signals at specific intervals (1 out of 3 signals on a total of 15 signals is shown for clarity), peaks at 10 and 20 are rising in intensity whereas the peaks at 30 and 40 have a decreasing intensity 2DCorrelationDemoDataset.svg
Demo dataset consisting of signals at specific intervals (1 out of 3 signals on a total of 15 signals is shown for clarity), peaks at 10 and 20 are rising in intensity whereas the peaks at 30 and 40 have a decreasing intensity

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. See also the figure with the original dataset on the right and the corresponding 2D spectrum in the figure below.

Presence of 2D spectra

Schematic presence of a 2D correlation spectrum with peak positions represented by dots. Region A is the main diagonal containing autopeaks, off-diagonal regions B contain cross-peaks. 2DCorrelationSpectrumPresence.svg
Schematic presence of a 2D correlation spectrum with peak positions represented by dots. Region A is the main diagonal containing autopeaks, off-diagonal regions B contain cross-peaks.

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 the 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. The signals of the original dataset are generally preprocessed. The original spectra are compared to a reference spectrum. By subtracting a reference spectrum, often the average spectrum of the dataset, so called dynamic spectra are calculated which form the corresponding dynamic dataset E. The presence and interpretation may be dependent on the choice of reference spectrum. The equations below are valid for equally spaced measurements of the perturbation.

Calculation of the synchronous spectrum

A 2D synchronous spectrum expresses the similarity between spectral of the data in the original dataset. In generalized 2D correlation spectroscopy this is mathematically expressed as covariance (or correlation). [6]

where:

Calculation of the asynchronous spectrum

Orthogonal spectra to the dynamic dataset E are obtained with the Hilbert-transform:

where:

The values of N, Nj, k are determined as follows:

where:

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

Autocorrelation signal on the main diagonal of the synchronous 2D spectrum of the figure below (arbitrary axis units) 2DCorrelationSpectrumAC.svg
Autocorrelation signal on the main diagonal of the synchronous 2D spectrum of the figure below (arbitrary axis units)

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. In the 2D synchronous main diagonal signal on the right 4 peaks are visible at 10, 20, 30, and 40 (see also the 4 corresponding positive autopeaks in the 2D synchronous spectrum on the right). This indicates that in the original dataset 4 peaks of changing intensity are present. 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. [7]

Determining the direction of intensity change

Example of a two-dimensional correlation spectrum. Open circles in this simplified view represent positive peaks, while discs represent negative peaks 2DCorrelationSpectrum.svg
Example of a two-dimensional correlation spectrum. Open circles in this simplified view represent positive peaks, while discs represent negative peaks

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:

  1. 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
  2. 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

As can be seen in the 2D synchronous spectrum on the right, the intensity changes of the peaks at 10 and 30 are related and the intensity of the peak at 10 and 30 changes in the opposite direction (negative cross-peak at (10,30)). The same is true for the peaks at 20 and 40.

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:

  1. 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
  2. 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
  3. if the change at x mainly precedes the change in the band at y, the asynchronous 2D cross peak at (x,y) is positive
  4. if the change at x mainly follows the change in the band at y, the asynchronous 2D cross peak at (x,y) is negative
  5. if the synchronous 2D cross peak at (x,y) is negative, the interpretation of rule 3 and 4 for the asynchronous 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.

Following the rules above. It can be derived that the changes at 10 and 30 occur simultaneously and the changes in intensity at 20 and 40 occur simultaneously as well. Because of the positive asynchronous cross-peak at (10, 20), the changes at 10 and 30 (predominantly) occur before the intensity changes at 20 and 40.

In some cases the Noda rules cannot be so readily implied, predominately when spectral features are not caused by simple intensity variations. This may occur when band shifts occur, or when a very erratic intensity variation is present in a given frequency range.

See also

Related Research Articles

<span class="mw-page-title-main">Infrared spectroscopy</span> Interaction of infrared radiation with matter

Infrared spectroscopy is the measurement of the interaction of infrared radiation with matter by absorption, emission, or reflection. It is used to study and identify chemical substances or functional groups in solid, liquid, or gaseous forms. It can be used to characterize new materials or identify and verify known and unknown samples. The method or technique of infrared spectroscopy is conducted with an instrument called an infrared spectrometer which produces an infrared spectrum. An IR spectrum can be visualized in a graph of infrared light absorbance on the vertical axis vs. frequency, wavenumber or wavelength on the horizontal axis. Typical units of wavenumber used in IR spectra are reciprocal centimeters, with the symbol cm−1. Units of IR wavelength are commonly given in micrometers, symbol μm, which are related to the wavenumber in a reciprocal way. A common laboratory instrument that uses this technique is a Fourier transform infrared (FTIR) spectrometer. Two-dimensional IR is also possible as discussed below.

<span class="mw-page-title-main">Spectroscopy</span> Study involving matter and electromagnetic radiation

Spectroscopy is the field of study that measures and interprets the electromagnetic spectra that result from the interaction between electromagnetic radiation and matter as a function of the wavelength or frequency of the radiation. Matter waves and acoustic waves can also be considered forms of radiative energy, and recently gravitational waves have been associated with a spectral signature in the context of the Laser Interferometer Gravitational-Wave Observatory (LIGO)

Fourier-transform spectroscopy is a measurement technique whereby spectra are collected based on measurements of the coherence of a radiative source, using time-domain or space-domain measurements of the radiation, electromagnetic or not. It can be applied to a variety of types of spectroscopy including optical spectroscopy, infrared spectroscopy, nuclear magnetic resonance (NMR) and magnetic resonance spectroscopic imaging (MRSI), mass spectrometry and electron spin resonance spectroscopy.

<span class="mw-page-title-main">X-ray photoelectron spectroscopy</span> Spectroscopic technique

X-ray photoelectron spectroscopy (XPS) is a surface-sensitive quantitative spectroscopic technique based on the photoelectric effect that can identify the elements that exist within a material or are covering its surface, as well as their chemical state, and the overall electronic structure and density of the electronic states in the material. XPS is a powerful measurement technique because it not only shows what elements are present, but also what other elements they are bonded to. The technique can be used in line profiling of the elemental composition across the surface, or in depth profiling when paired with ion-beam etching. It is often applied to study chemical processes in the materials in their as-received state or after cleavage, scraping, exposure to heat, reactive gasses or solutions, ultraviolet light, or during ion implantation.

Rotational–vibrational spectroscopy is a branch of molecular spectroscopy concerned with infrared and Raman spectra of molecules in the gas phase. Transitions involving changes in both vibrational and rotational states can be abbreviated as rovibrational transitions. When such transitions emit or absorb photons, the frequency is proportional to the difference in energy levels and can be detected by certain kinds of spectroscopy. Since changes in rotational energy levels are typically much smaller than changes in vibrational energy levels, changes in rotational state are said to give fine structure to the vibrational spectrum. For a given vibrational transition, the same theoretical treatment as for pure rotational spectroscopy gives the rotational quantum numbers, energy levels, and selection rules. In linear and spherical top molecules, rotational lines are found as simple progressions at both higher and lower frequencies relative to the pure vibration frequency. In symmetric top molecules the transitions are classified as parallel when the dipole moment change is parallel to the principal axis of rotation, and perpendicular when the change is perpendicular to that axis. The ro-vibrational spectrum of the asymmetric rotor water is important because of the presence of water vapor in the atmosphere.

In physics and physical chemistry, time-resolved spectroscopy is the study of dynamic processes in materials or chemical compounds by means of spectroscopic techniques. Most often, processes are studied after the illumination of a material occurs, but in principle, the technique can be applied to any process that leads to a change in properties of a material. With the help of pulsed lasers, it is possible to study processes that occur on time scales as short as 10−16 seconds. All time-resolved spectra are suitable to be analyzed using the two-dimensional correlation method for a correlation map between the peaks.

Resonance Raman spectroscopy is a Raman spectroscopy technique in which the incident photon energy is close in energy to an electronic transition of a compound or material under examination. The frequency coincidence can lead to greatly enhanced intensity of the Raman scattering, which facilitates the study of chemical compounds present at low concentrations.

<span class="mw-page-title-main">Nuclear magnetic resonance spectroscopy</span> Laboratory technique

Nuclear magnetic resonance spectroscopy, most commonly known as NMR spectroscopy or magnetic resonance spectroscopy (MRS), is a spectroscopic technique to observe local magnetic fields around atomic nuclei. The sample is placed in a magnetic field and the NMR signal is produced by excitation of the nuclei sample with radio waves into nuclear magnetic resonance, which is detected with sensitive radio receivers. The intramolecular magnetic field around an atom in a molecule changes the resonance frequency, thus giving access to details of the electronic structure of a molecule and its individual functional groups. As the fields are unique or highly characteristic to individual compounds, in modern organic chemistry practice, NMR spectroscopy is the definitive method to identify monomolecular organic compounds.

<span class="mw-page-title-main">Electron paramagnetic resonance</span> Technique to study materials that have unpaired electrons

Electron paramagnetic resonance (EPR) or electron spin resonance (ESR) spectroscopy is a method for studying materials that have unpaired electrons. The basic concepts of EPR are analogous to those of nuclear magnetic resonance (NMR), but the spins excited are those of the electrons instead of the atomic nuclei. EPR spectroscopy is particularly useful for studying metal complexes and organic radicals. EPR was first observed in Kazan State University by Soviet physicist Yevgeny Zavoisky in 1944, and was developed independently at the same time by Brebis Bleaney at the University of Oxford.

The heteronuclear single quantum coherence or heteronuclear single quantum correlation experiment, normally abbreviated as HSQC, is used frequently in NMR spectroscopy of organic molecules and is of particular significance in the field of protein NMR. The experiment was first described by Geoffrey Bodenhausen and D. J. Ruben in 1980. The resulting spectrum is two-dimensional (2D) with one axis for proton (1H) and the other for a heteronucleus, which is usually 13C or 15N. The spectrum contains a peak for each unique proton attached to the heteronucleus being considered. The 2D HSQC can also be combined with other experiments in higher-dimensional NMR experiments, such as NOESY-HSQC or TOCSY-HSQC.

Two-dimensional nuclear magnetic resonance spectroscopy is a set of nuclear magnetic resonance spectroscopy (NMR) methods which give data plotted in a space defined by two frequency axes rather than one. Types of 2D NMR include correlation spectroscopy (COSY), J-spectroscopy, exchange spectroscopy (EXSY), and nuclear Overhauser effect spectroscopy (NOESY). Two-dimensional NMR spectra provide more information about a molecule than one-dimensional NMR spectra and are especially useful in determining the structure of a molecule, particularly for molecules that are too complicated to work with using one-dimensional NMR.

<span class="mw-page-title-main">Angle-resolved photoemission spectroscopy</span> Experimental technique to determine the distribution of electrons in solids

Angle-resolved photoemission spectroscopy (ARPES) is an experimental technique used in condensed matter physics to probe the allowed energies and momenta of the electrons in a material, usually a crystalline solid. It is based on the photoelectric effect, in which an incoming photon of sufficient energy ejects an electron from the surface of a material. By directly measuring the kinetic energy and emission angle distributions of the emitted photoelectrons, the technique can map the electronic band structure and Fermi surfaces. ARPES is best suited for the study of one- or two-dimensional materials. It has been used by physicists to investigate high-temperature superconductors, graphene, topological materials, quantum well states, and materials exhibiting charge density waves.

In statistical signal processing, the goal of spectral density estimation (SDE) or simply spectral estimation is to estimate the spectral density of a signal from a sequence of time samples of the signal. Intuitively speaking, the spectral density characterizes the frequency content of the signal. One purpose of estimating the spectral density is to detect any periodicities in the data, by observing peaks at the frequencies corresponding to these periodicities.

In spectroscopy, spectral flux density is the quantity that describes the rate at which energy is transferred by electromagnetic radiation through a real or virtual surface, per unit surface area and per unit wavelength. It is a radiometric rather than a photometric measure. In SI units it is measured in W m−3, although it can be more practical to use W m−2 nm−1 or W m−2 μm−1, and respectively by W·m−2·Hz−1, Jansky or solar flux units. The terms irradiance, radiant exitance, radiant emittance, and radiosity are closely related to spectral flux density.

<span class="mw-page-title-main">Two-dimensional infrared spectroscopy</span>

Two-dimensional infrared spectroscopy is a nonlinear infrared spectroscopy technique that has the ability to correlate vibrational modes in condensed-phase systems. This technique provides information beyond linear infrared spectra, by spreading the vibrational information along multiple axes, yielding a frequency correlation spectrum. A frequency correlation spectrum can offer structural information such as vibrational mode coupling, anharmonicities, along with chemical dynamics such as energy transfer rates and molecular dynamics with femtosecond time resolution. 2DIR experiments have only become possible with the development of ultrafast lasers and the ability to generate femtosecond infrared pulses.

<span class="mw-page-title-main">Nuclear magnetic resonance</span> Spectroscopic technique based on change of nuclear spin state

Nuclear magnetic resonance (NMR) is a physical phenomenon in which nuclei in a strong constant magnetic field are perturbed by a weak oscillating magnetic field and respond by producing an electromagnetic signal with a frequency characteristic of the magnetic field at the nucleus. This process occurs near resonance, when the oscillation frequency matches the intrinsic frequency of the nuclei, which depends on the strength of the static magnetic field, the chemical environment, and the magnetic properties of the isotope involved; in practical applications with static magnetic fields up to ca. 20 tesla, the frequency is similar to VHF and UHF television broadcasts (60–1000 MHz). NMR results from specific magnetic properties of certain atomic nuclei. Nuclear magnetic resonance spectroscopy is widely used to determine the structure of organic molecules in solution and study molecular physics and crystals as well as non-crystalline materials. NMR is also routinely used in advanced medical imaging techniques, such as in magnetic resonance imaging (MRI).

<span class="mw-page-title-main">Kolmogorov–Zurbenko filter</span>

Within statistics, the Kolmogorov–Zurbenko (KZ) filter was first proposed by A. N. Kolmogorov and formally defined by Zurbenko. It is a series of iterations of a moving average filter of length m, where m is a positive, odd integer. The KZ filter belongs to the class of low-pass filters. The KZ filter has two parameters, the length m of the moving average window and the number of iterations k of the moving average itself. It also can be considered as a special window function designed to eliminate spectral leakage.

Vibronic spectroscopy is a branch of molecular spectroscopy concerned with vibronic transitions: the simultaneous changes in electronic and vibrational energy levels of a molecule due to the absorption or emission of a photon of the appropriate energy. In the gas phase, vibronic transitions are accompanied by changes in rotational energy also.

Spectral line shape describes the form of a feature, observed in spectroscopy, corresponding to an energy change in an atom, molecule or ion. This shape is also referred to as the spectral line profile. Ideal line shapes include Lorentzian, Gaussian and Voigt functions, whose parameters are the line position, maximum height and half-width. Actual line shapes are determined principally by Doppler, collision and proximity broadening. For each system the half-width of the shape function varies with temperature, pressure and phase. A knowledge of shape function is needed for spectroscopic curve fitting and deconvolution.

Two-dimensional electronic spectroscopy (2DES) is an ultrafast laser spectroscopy technique that allows the study of ultrafast phenomena inside systems in condensed phase. The term electronic refers to the fact that the optical frequencies in the visible spectral range are used to excite electronic energy states of the system; however such a technique is also used in the IR optical range and in this case the method is called two-dimensional infrared spectroscopy (2DIR). This technique records the signal which is emitted from a system after an interaction with a sequence of 3 laser pulses. Such pulses usually have a time duration of few hundred femtosecond and this high time resolution allows capturing of dynamics inside the system that evolves with the same time scale. The main result of this technique is a two-dimensional absorption spectrum that shows the correlation between excitation and detection frequencies.

References

  1. ↑ Shin-Ichi Morita; Yasuhiro F. Miura; Michio Sugi & 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. Bibcode:2005CPL...402..251M. doi:10.1016/j.cplett.2004.12.038.
  2. ↑ Koichi Murayama; Boguslawa Czarnik-Matusewicz; Yuqing Wu; Roumiana Tsenkova & 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. Bibcode:2000ApSpe..54..978M. doi:10.1366/0003702001950715. S2CID   95843070.
  3. 1 2 Shin-Ichi Morita & 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. Bibcode:2002ApSpe..56..502M. doi:10.1366/0003702021954953. S2CID   95679157.
  4. 1 2 3 4 Isao Noda & Yukihiro Ozaki (2004). Two-Dimensional Correlation Spectroscopy - Applications in Vibrational and Optical Spectroscopy. John Wiley & Sons Ltd. ISBN   978-0-471-62391-5.
  5. ↑ Boguslawa Czarnik-Matusewicz; Sylwia Pilorz; Lorna Ashton & Ewan W. Blanch (2006). "Potential pitfalls concerning visualization of the 2D results". Journal of Molecular Structure. 799 (1–3): 253–258. Bibcode:2006JMoSt.799..253C. doi:10.1016/j.molstruc.2006.03.064.
  6. ↑ Noda, Isao (1993). "Generalized Two-Dimensional Correlation Method Applicable to Infrared, Raman, and Other Types of Spectroscopy". Applied Spectroscopy. 47 (9): 1329–1336. Bibcode:1993ApSpe..47.1329N. doi:10.1366/0003702934067694. S2CID   94722664.
  7. ↑ R. Buchet, Y. Wu; G. Lachenal; C. Raimbault & Yukihiro Ozaki (2006). "Selecting two-dimensional cross-correlation functions to enhance interpretation of near-infrared spectra of proteins". Applied Spectroscopy. 55 (2): 155–162. Bibcode:2001ApSpe..55..155B. doi:10.1366/0003702011951452. S2CID   95827191.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.