Ep 41 — Two-Dimensional Correlation Spectroscopy (2D-COS): A Powerful Tool for Resolving Overlapping Peaks
Series: Encyclopedia of Infrared Spectroscopy: From Principles to Practice
Chapter: Part 4 · Advanced Topics — Frontier Techniques (First Episode of Second Half)
Target Audience: Graduate students in analytical/physical/polymer chemistry, spectroscopy researchers, biophysicists and drug polymorph investigators
Prerequisites: Ep 04 (Reading IR spectra), Ep 18 (Basic spectral processing), Ep 20 (Beer-Lambert law), Ep 22 (Proteins and drug polymorphs), Ep 24 (Polymer aging), Ep 43 (Chemometrics PCA/PLS, recommended parallel reading)
Reading Time: Approximately 55 minutes
Introduction: When 1D Spectra Are "Crowded"
In 1986, Iso Noda, while studying phase separation behavior of polymer blends at Procter & Gamble (P&G), encountered a challenging problem: the spectral peaks of different components severely overlapped, and conventional one-dimensional infrared spectra could not distinguish which component underwent conformational change first [1]. Inspired by 2D NMR spectroscopy, he proposed a key idea—using an external perturbation as a "second dimension" to correlate a series of spectra varying with the perturbation, generating Two-Dimensional Correlation Spectroscopy (2D-COS). In 1989, he formally extended this method to Generalized 2D-COS, allowing arbitrary selection of perturbation variables [2].
"Two-dimensional correlation spectroscopy converts a set of spectra obtained as a function of a perturbation into a synchronous and an asynchronous 2D spectrum, revealing interactions between bands that are invisible in the original 1D spectra."
—— Noda I. Appl Spectrosc 1993 [2]
The three "killer features" of 2D-COS provide unprecedented capabilities to researchers:
- Enhanced apparent resolution: Multiple component peaks originally merged into a single envelope in 1D spectra can be separated as cross peaks in 2D spectra;
- Determination of sequential order: Using the sign rules of synchronous and asynchronous spectra, the order of response of different groups to perturbation can be determined;
- Revealing intermolecular interactions: Off-diagonal peaks indicate which groups change cooperatively, whether coupled or competitive.
In this episode, we will systematically discuss the mathematical principles, interpretation rules, software implementation, and typical applications of 2D-COS, supplemented by three full case studies: protein folding, polymer stretching, and mixture analysis.
1. Mathematical Principles of 2D-COS
1.1 From 1D Perturbation Sequence to 2D Correlation
The input to 2D-COS is a set of spectra collected under a perturbation {I(ν₁, t₁), I(ν₁, t₂), ..., I(ν₁, tₘ)}, where ν is wavenumber and t is the perturbation variable (e.g., temperature, time, concentration, pressure, pH).
Mathematically, generalized 2D-COS first mean-centers the perturbation sequence to obtain dynamic spectra $\tilde{y}_j(\nu)=y_j(\nu)-\bar{y}(\nu)$, then constructs synchronous and asynchronous spectra (Hilbert–Noda form) [1][2][3]:
$$ \Phi(\nu_1,\nu_2)=\frac{1}{m-1}\sum_{j=1}^{m}\tilde{y}_j(\nu_1)\,\tilde{y}_j(\nu_2) $$
$$ \Psi(\nu_1,\nu_2)=\frac{1}{m-1}\sum_{j=1}^{m}\tilde{y}_j(\nu_1)\,\tilde{z}_j(\nu_2) $$
where $\tilde{z}_j(\nu)$ is the Hilbert–Noda transform of the dynamic spectra along the perturbation axis (discrete implementation in Python in the next section; the continuous formulation is equivalent to orthogonal phase shifting in the perturbation domain followed by cross-correlation). $\Phi$ is the synchronous spectrum, and $\Psi$ is the asynchronous spectrum.
(Previous drafts incorrectly used a summation with an exponent; corrected here to the standard form consistent with the later code.)
1.2 Synchronous Spectrum: Who Changes Together
The synchronous spectrum Φ(ν₁, ν₂) reflects the extent to which absorbances at two wavenumbers change synchronously:
- Autopeaks: located at ν₁ = ν₂, i.e., the main diagonal of the spectrum, with intensity equal to the autocorrelation variance at that wavenumber. Autopeaks are always positive, indicating the sensitivity of that group to the perturbation.
- Cross peaks: located at ν₁ ≠ ν₂. If two groups change synchronously (in the same or opposite direction), a synchronous cross peak appears; its sign indicates the direction of change (positive for same direction, negative for opposite).
ν₂
↑
┌─────────────┐
│ ⋱ │ Diagonal = autopeaks
│ ⋱ │ Off-diagonal = cross peaks
│ ● ⋱ │ ● = synchronous cross peak (same direction)
│ ● │
│ ⋱ │
└─────────────┘ ν₁ →
Figure 1: Schematic of synchronous spectrum structure. "Bright spots" on the diagonal are autopeaks; circular spots off the diagonal are synchronous cross peaks.
1.3 Asynchronous Spectrum: Who Comes First
The asynchronous spectrum Ψ(ν₁, ν₂) reflects the extent to which absorbances at two wavenumbers change asynchronously (out of phase), and is the core value of 2D-COS [2][3].
- The asynchronous spectrum has zero intensity on the diagonal (no autocorrelation);
- An asynchronous cross peak Ψ(ν₁, ν₂) ≠ 0 indicates that the two groups respond to the perturbation at different times (or perturbation values);
- By combining the signs of Φ and Ψ, the sequential order of response between two groups can be determined.
1.4 Noda's Rule
In 1993, Noda proposed the famous "eight-character" rule (Noda's rule) [2][3], which determines the order of response of groups at ν₁ and ν₂ to the perturbation based on the signs of synchronous and asynchronous spectra:
| Synchronous Φ(ν₁, ν₂) | Asynchronous Ψ(ν₁, ν₂) | Conclusion (response order of ν₁ relative to ν₂) |
|---|---|---|
| + | + | ν₁ occurs before ν₂ |
| + | − | ν₁ occurs after ν₂ |
| − | + | ν₁ occurs after ν₂ |
| − | − | ν₁ occurs before ν₂ |
Table 1: Noda's rule (with ν₁ as the horizontal axis and ν₂ as the vertical axis)
💡 Interpretation tip: When synchronous Φ = 0, the two groups change orthogonally, and the order cannot be determined; assumptions about the perturbation direction are needed. When asynchronous Ψ = 0, the two groups change strictly synchronously (possibly from the same physicochemical process).
1.5 Generalized 2D Correlation Spectroscopy: Arbitrary Perturbations
In 1989, Noda proposed Generalized Two-Dimensional Correlation Spectroscopy (Generalized 2D-COS) [2], extending the perturbation variable from the original "time" to any physical quantity:
- Temperature (heating/cooling): phase transitions, denaturation, glass transition, polymerization reactions;
- Concentration: titration, mixing;
- Pressure: high-pressure phase transitions;
- pH: protein acid-base titration;
- Electric/magnetic fields: orientation, electromagnetic induction;
- Strain: stretching, shear;
- Time: reaction kinetics, photochemistry, degradation;
- Spatial position: "pseudo-2D-COS" of line scans or imaging data.
Generalization has allowed 2D-COS to transcend the infrared field and be widely applied to Raman, UV-Vis, NMR, XRD, CD and other spectroscopic techniques [3].
2. Advantages and Limitations of 2D-COS
2.1 Three Core Advantages
Enhanced apparent resolution: Peaks originally merged into an "envelope" in 1D spectra can be better separated via cross peaks in 2D spectra. Note: This is apparent peak separation ability, not an increase in instrumental spectral resolution by a fixed factor; the literature occasionally mentions "several-fold" empirical statements, which should be understood qualitatively rather than as a universal specification [2].
Peak Assignment Confirmation: The intensity of diagonal autocorrelation peaks quantifies the sensitivity of each group to the perturbation, thus inversely inferring peak assignments—for example, which amide III sub-band corresponds to α-helix and which to β-sheet.
Kinetic Order Inference: Asynchronous spectra can infer the order of molecular events from a series of static spectra without time-resolved spectroscopy—this is the most "scientifically valuable" capability of 2D-COS.
2.2 Limitations and Pitfalls
The power of 2D-COS comes with several common pitfalls:
- Baseline Drift Sensitivity: Uncorrected baseline drift or residual baseline drift can be misinterpreted as "perturbation response," leading to false cross peaks [3][4]. It is recommended to perform SNV or polynomial baseline correction on each spectrum before 2D-COS.
- Spectrum Normalization: Overall intensity variations across spectra (e.g., sample thickness differences) introduce spurious correlations. Normalization to a uniform scale is mandatory.
- Perturbation Direction Dependence: The sign of asynchronous spectra depends on the perturbation direction (heating vs. cooling). Conclusions must explicitly state the perturbation direction.
- Oversampling and Undersampling: Too few perturbation points (< 5) lose statistical significance; too many (> 50) introduce noise. 10–20 perturbation points are recommended [3].
- Noise is "Amplified": Since 2D-COS is essentially a covariance operation, noise in the original spectra is amplified as spurious peaks in the 2D spectra. Original SNR should be sufficiently high.
⚠️ Common Misuse: Performing 2D-COS on a single spectrum is incorrect—2D-COS requires a series of spectra as input, not a single spectrum. This is the most common mistake among beginners.
3. Software Implementation: 2Dshige and Python
3.1 2Dshige: The Classic Japanese Tool
2Dshige (2D Shige) was developed by the Ozaki group at Kwansei Gakuin University, Japan, and is the earliest free software for 2D-COS [5]:
- Platform: Windows;
- Input: Text format spectral matrix (rows = perturbation points, columns = wavenumbers);
- Output: Synchronous/Asynchronous 2D spectra (BMP/JPG);
- Advantages: Graphical interface, simple operation, free for academic users;
- Limitations: Old Windows GUI, closed source, not suitable for batch processing.
Download: https://sci-tech.kwansei.ac.j…
3.2 Python Implementation: numpy + matplotlib
Below is a minimal Python implementation (~30 lines), suitable for researchers to get started quickly:
import numpy as np
import matplotlib.pyplot as plt
def two_d_cos(spectra):
"""
spectra: (n_perturbations, n_wavenumbers) array,
rows correspond to different perturbation states, columns to wavenumbers.
Returns Phi (synchronous), Psi (asynchronous) 2D matrices.
"""
# Mean centering (dynamic spectra)
y_tilde = spectra - spectra.mean(axis=0, keepdims=True)
m = y_tilde.shape[0]
# Hilbert-Noda transformation matrix
j = np.arange(1, m + 1).reshape(-1, 1)
k = np.arange(1, m + 1).reshape(1, -1)
N = np.zeros((m, m))
# N[j, k] = 0 if j==k else 1 / (pi * (k - j))
diff = k.T - j.T # (m, m), N_{jk} = 1/(pi*(k-j))
N = np.where(np.eye(m) == 0, 1.0 / (np.pi * diff), 0.0)
# Synchronous spectrum: covariance
Phi = y_tilde.T @ y_tilde / (m - 1)
# Asynchronous spectrum: cross-covariance of dynamic spectra and Hilbert transform
Psi = y_tilde.T @ (N @ y_tilde) / (m - 1)
return Phi, Psi
# Example: simulate a set of amide I spectra varying with temperature
wavenumber = np.linspace(1700, 1600, 200)
temperatures = np.linspace(25, 95, 15) # 15 temperature points
spectra = np.zeros((len(temperatures), len(wavenumber)))
for i, T in enumerate(temperatures):
# α-helix (1650 cm⁻¹) decreases with temperature; β-sheet (1630 cm⁻¹) increases with temperature
spectra[i] = (1 - (T-25)/70) * np.exp(-(wavenumber-1650)**2 / (2*5**2)) \
+ ((T-25)/70) * np.exp(-(wavenumber-1630)**2 / (2*5**2))
Phi, Psi = two_d_cos(spectra)
fig, axes = plt.subplots(1, 2, figsize=(12, 5))
axes[0].imshow(Phi, extent=[wavenumber[0], wavenumber[-1], wavenumber[-1], wavenumber[0]])
axes[0].set_title("Synchronous Φ"); axes[0].set_xlabel("ν₂"); axes[0].set_ylabel("ν₁")
axes[1].imshow(Psi, extent=[wavenumber[0], wavenumber[-1], wavenumber[-1], wavenumber[0]])
axes[1].set_title("Asynchronous Ψ"); axes[1].set_xlabel("ν₂"); axes[1].set_ylabel("ν₁")
plt.tight_layout(); plt.savefig("2dcos_demo.png", dpi=150)
📦 Ready-made Python libraries: scikit-learn-style
py2D-COSand2DCoSare both open source, can be installed directly via pip, avoiding reinventing the wheel.
3.3 Other Available Tools
| Software | Type | Platform | Notes |
|---|---|---|---|
| 2Dshige | GUI | Windows | Free for academic use, most classic [5] |
| py2D-COS | Python | Cross-platform | pip install, API friendly |
| 2DCoS (Python) | Python | Cross-platform | Includes example dataset |
| MATLAB 2D-COS Toolbox | MATLAB | Cross-platform | Early implementation from Ozaki group |
| OriginPro 2D-COS Plugin | GUI | Windows | Commercial, integrated with Origin ecosystem |
| Quasar (Orange-Spectroscopy) | GUI | Cross-platform | Drag-and-drop workflow |
Table 2: Summary of mainstream 2D-COS tools
4. Application Case 1: 2D-COS Analysis of Protein Thermal Denaturation
4.1 The Challenge of Amide I Sub-bands
The amide I band (Amide I, 1600–1700 cm⁻¹) of proteins is mainly composed of C=O stretching vibrations and serves as a core probe for protein secondary structure. However, the amide I sub-bands of different secondary structures (α-helix, β-sheet, turns, random coil) differ by only 5–15 cm⁻¹ from each other, merging into a broad band in one-dimensional spectra [6]:
| Secondary Structure | Amide I Frequency (cm⁻¹) |
|---------|----------------------|
| α-helix | ~1650–1658 |
| β-sheet (parallel) | ~1630–1640 |
| β-sheet (antiparallel) | ~1620–1635 and ~1685–1695 (doublet) |
| random coil | ~1640–1650 |
| turns | ~1660–1675 |
| β-turns (antiparallel β-turns) | ~1675–1695 |
Table 3: Assignments of protein amide I sub-bands (data from Barth Biochim Biophys Acta 2007 [6])
🔗 Further reading: For detailed infrared assignments of amide I, II, and III sub-bands, see ftir.fun amide functional group page. This page systematically organizes the specific frequencies and assignment criteria for different secondary structures and side-chain environments, serving as an essential reference for 2D-COS amide band analysis.
4.2 2D-COS Reveals Thermal Denaturation Order
In 1999, Fabian and Mantsch published a classic study in Vib Spectrosc [7]: they used lysozyme solution as the sample, heated from 25°C to 90°C, collected one FTIR spectrum every 5°C (14 spectra total), and performed 2D-COS:
- Synchronous spectrum: Autopeaks appear at 1654, 1634, 1675 cm⁻¹, directly proving that these three frequencies are independent sub-bands (completely unresolvable in 1D spectra).
- Asynchronous spectrum: Strong asynchronous peak between 1654 (α-helix) and 1675 (turns), Φ = +, Ψ = + → α-helix changes before turns.
- Further analysis: 1634 (β-sheet) ↔ 1654 (α-helix), Φ = +, Ψ = − → β-sheet changes after α-helix.
Thus, the thermal denaturation order is: α-helix loosening → turn rearrangement → β-sheet exposure — this result was later confirmed by NMR and molecular dynamics simulations, becoming a classic case in protein denaturation mechanism research [7].
4.3 Workflow
- Sample preparation: 2 mg/mL protein in D₂O solution to avoid H₂O interference at 1640 cm⁻¹;
- Measurement: Variable temperature accessory, 5°C steps, equilibrate 5 min, 256 scans, resolution 4 cm⁻¹;
- Preprocessing: Water vapor subtraction (refer to Ep 28), second derivative (S-G 9 points);
- 2D-COS: Use temperature as perturbation variable, calculate Φ and Ψ;
- Interpretation: Deduce sub-band change order according to Noda's rules.
💡 Tip: In D₂O, amide I shifts to ~1450 cm⁻¹ (amide II disappears), yielding cleaner spectra; D₂O is the "standard" solvent for protein 2D-COS [6].
5. Application Case 2: 2D-COS of Polymer Stretching Orientation
5.1 Stress-Spectrum Synchronous Measurement
During polymer stretching, molecular chains orient along the stress direction; the orientation order of different groups reflects molecular chain rigidity and cooperative motion ability [8][9]. In the early 2000s, Vidrine and Galloway et al. established the dynamic infrared dichroism 2D-COS method [8]:
- Setup: FTIR + polarizer + micro-tensile machine;
- Perturbation: Strain ε (0–50%);
- Collection: One spectrum per 5% strain, 11 spectra total.
5.2 Stretching Order of PET
PET contains ester C=O (1715 cm⁻¹), benzene ring C=C (1410, 1505, 1580 cm⁻¹), ether C–O–C (1240, 1095 cm⁻¹), aliphatic CH₂ (2970, 2900 cm⁻¹). Study [8] found:
- Synchronous spectrum: 1715 ↔ 1410 ↔ 1095 positive correlation → these groups orient synchronously;
- Asynchronous spectrum: 1095 (C–O–C) ↔ 2970 (CH₂), Φ = +, Ψ = + → ether orients before CH₂;
- Asynchronous spectrum: 1715 (C=O) ↔ 1095 (C–O–C), Φ = +, Ψ = − → C=O orients after C–O–C.
Molecular mechanism inference: The benzene ring-ester segment (rigid segment) orients first along the stress direction; the flexible CH₂ segment subsequently relaxes and rearranges [8][9]. This conclusion agrees with solid-state NMR results.
5.3 2D-COS of Blend Phase Separation
For polymer blends (e.g., PMMA/PVDF), infrared frequencies of different components (C=O, CF₂) are close, making 1D spectra difficult to distinguish. Using temperature as perturbation, 2D-COS can resolve [9]:
- PMMA's C=O at ~1730 cm⁻¹, PVDF's CF₂ at ~1180 cm⁻¹; they are normally uncorrelated;
- But near the phase separation temperature, asynchronous cross peaks appear between 1730 and 1180 cm⁻¹, reflecting the "phase separation onset" of the two polymers near the critical temperature.
🔗 Further reading: For assignments of C=O, C-O-C and other functional groups in polymers, see ftir.fun carbonyl functional group page and ftir.fun ether functional group page.
6. Application Case 3: 2D-COS Analysis of Mixtures
6.1 "Secondary Resolution" of Overlapping Peaks
Consider IR spectra of ethanol-acetone mixtures: ethanol's C–O stretch at ~1050 cm⁻¹, acetone's C=O at ~1715 cm⁻¹, seemingly not overlapping. However, if the mixture contains multiple alcohols (methanol, ethanol, propanol, butanol), their C–O stretches lie at 1030, 1050, 1070, 1100 cm⁻¹, merging into one broad band in 1D spectra, unresolvable [3][10].
Solution: Use concentration as perturbation variable — gradually change the proportion of each alcohol in the mixture, perform 2D-COS:
- Synchronous spectrum: Four independent autopeaks appear at 1030, 1050, 1070, 1100 cm⁻¹, directly proving the presence of 4 components;
- Asynchronous spectrum: Asynchronous peak at 1030 ↔ 1100 proves that the concentration changes of the two components are not synchronous (e.g., one increases while the other decreases).
This method has been applied in [10]:
- Multi-component analysis of traditional Chinese medicine (e.g., tanshinones and salvianolic acids in Salvia miltiorrhiza);
- Food adulteration detection (e.g., honey adulterated with high-fructose corn syrup, C–H region overlap);
- Simultaneous identification of oxygenates in fuels (methanol, ethanol, MTBE).
6.2 Integration of 2D-COS with Chemometrics
The 2D matrix provided by 2D-COS can be used as feature vectors input to PCA, PLS models, improving model discrimination ability [3][10]:
- 2D-COS-PCA: Concatenate diagonal elements of synchronous spectrum and off-diagonal elements of asynchronous spectrum into a feature vector, perform PCA clustering;
- 2D-COS-PLS: Use 2D spectra as X matrix, perform quantitative regression;
- Practice: In olive oil adulteration identification, the R² of 2D-COS-PLS model increased from 0.86 (1D PLS) to 0.95 [10].
For detailed introduction of PCA and PLS, see next episodes Ep 42 hyphenated techniques and Ep 43 chemometrics special topic.
7. Extended Applications of Generalized 2D-COS
7.1 Sample-Sample 2D-COS
Traditional 2D-COS requires a continuous perturbation variable. Sample-Sample 2D-COS, proposed by Šašić and Ozaki [3], allows input as a set of independent samples (e.g., different batches of drugs) and still generates 2D spectra:
- Axes are sample numbers;
- Reveals similarities and outliers among samples;
- Applied to abnormal batch identification in industrial QC.
7.2 Heterospectral 2D-COS
Using data from two different spectroscopic techniques (e.g., IR and Raman, IR and NIR) as two dimensions, perform heterospectral correlation [3]:
- Synchronous cross peaks → cooperative response of the same group to both techniques;
- Asynchronous cross peaks → differences in response of a group between techniques (e.g., IR-active vs Raman-active).
Application: In drug polymorph studies, IR-Raman heterospectral 2D-COS can reveal vibrational differences of the same group in different crystal forms, more comprehensive than using either technique alone [3].
7.3 Moving Window 2D-COS (MW2D)
Moving-Window 2D-COS divides the full temperature range into multiple sub-windows, performing local 2D-COS for each sub-window to obtain a three-dimensional map of "temperature-wavenumber-correlation intensity" [11]:
- Advantages: Precisely locating phase transition temperatures (e.g., glass transition Tg);
- Applications: Polymer blend phase separation, liquid crystal phase transitions, protein unfolding temperature.
VIII. Experimental Design and Considerations
8.1 Selection of Perturbation Variables
| Research Objective | Recommended Perturbation | Spectral Sampling |
|---|---|---|
| Protein conformational change | Temperature | 5°C/step, 25–95°C, 14 spectra |
| Polymer stretching orientation | Strain | 5%/step, 0–50%, 11 spectra |
| Titration endpoint determination | Concentration | Titration 0.1 mL/step, 10–20 spectra |
| Reaction kinetics | Time | 1 s/step (rapid scan) or 1 min/step |
| pH-induced phase transition | pH | 0.2 pH/step |
| Photochemical intermediates | Time (after UV excitation) | ms–μs time-resolved |
| Polymorph transition | Temperature/Pressure | 1°C/step near transition point |
Table 4: Common perturbation variable selections
8.2 Spectral Quality Requirements
2D-COS is a "garbage in, garbage out" method. The original spectra must satisfy [3][4]:
- Signal-to-noise ratio SNR > 1000: Because covariance operations amplify noise;
- Baseline drift < 0.5%T: Polynomial baseline correction or SNV recommended;
- Wavenumber accuracy < 0.5 cm⁻¹: Ensured by HeNe laser reference;
- Perturbation stability: Sufficient equilibration time (>3 min) at each perturbation point to avoid "overshoot".
8.3 Elimination of Anomalous Spectra
If a spectrum at a certain perturbation point is anomalous (e.g., sample bubbles, temperature overshoot), it should be eliminated rather than retained—a single anomalous point can significantly affect the 2D spectra. Common methods:
- PCA scores on the first principal component exceeding 3σ;
- Abnormally large spectral distances (e.g., Euclidean distance);
- Subjective visual inspection.
IX. 2D-COS Literature Resources and Advanced Learning
9.1 Classic Reviews
- Noda I. Appl Spectrosc 1993, 47: 1329–1336 [2] — Essential foundational paper;
- Noda I. Vib Spectrosc 2004, 36: 143–165 [3] — Generalized 2D-COS review;
- Ozaki Y, Šašić S, Tanaka T, Noda I. Bull Chem Soc Jpn 2003, 76: 1–13 — Review on polymer and biomolecule 2D-COS;
- Morita S, Ozaki Y, Noda I. Appl Spectrosc 2023, 77: 5–15 — Review on 30 years of 2D-COS progress.
9.2 Online Resources
- 2Dshige Software: https://sci-tech.kwansei.ac.j…
- International 2D-COS Symposium: Held biennially since 2003, led by the Asia-Pacific spectroscopy community;
- Noda Paper Index: https://www.researchgate.net/…
9.3 Advanced Topics
- Isotope-Labeled 2D-COS: Use ¹³C/¹⁵N labeling of specific amino acids to improve assignment accuracy in protein 2D-COS;
- Time-Resolved 2D-COS: Combine millisecond time-resolved IR to study photochemical intermediates;
- Two-Dimensional Electrochemical Infrared (2D-SEC-FTIR): Use electrode potential as perturbation to study electrochemical reaction intermediates.
Summary of This Section
| Key Concept | Key Points |
|---|---|
| Definition of 2D-COS | Cross-correlation of spectra under a perturbation → synchronous Φ + asynchronous Ψ 2D spectra |
| Proposer | Iso Noda, 1986 (P&G); generalized 1989 |
| Synchronous spectrum Φ | Diagonal = autocorrelation peaks; off-diagonal = synchronous cross peaks; positive/negative = same/different direction |
| Asynchronous spectrum Ψ | Diagonal always zero; off-diagonal = asynchronous cross peaks; determines order |
| Noda's rules | (+,+) ν₁ first; (+,−) ν₁ later; (−,+) ν₁ later; (−,−) ν₁ first |
| Three major advantages | Enhanced apparent peak separation; kinetic order; intermolecular interaction clues |
| Limitations | Baseline sensitive; noise amplification; requires multiple spectra |
| Generalized 2D-COS | Perturbation can be arbitrary (temperature, concentration, strain, pH, potential, ...) |
| Software | 2Dshige, py2D-COS, OriginPro plugin, Quasar |
| Amide I sub-bands | α-helix 1650–1658; β-sheet 1620–1640; turn 1660–1695 |
| Protein denaturation order | α-helix → turn → β-sheet (Fabian & Mantsch, 1999) |
| PET stretching order | Phenyl-ester segment → CH₂ segment |
| Experimental requirements | SNR > 1000; baseline drift < 0.5%T; perturbation sampling 10–20 points |
| Advanced topics | Sample-sample, heterospectral, moving-window 2D-COS |
Study Questions
A graduate student wants to study the cold crystallization process of polylactic acid (PLA). What perturbation variable should be chosen? How to design the sampling interval? What spectral pre-processing is needed before 2D-COS?
In a synchronous spectrum, a positive synchronous cross peak is observed at 1650 and 1630 cm⁻¹. What does it imply? Is it sufficient to conclude "two groups change synchronously" based on this alone? Why?
Suppose you perform 2D-COS on a protein and obtain synchronous Φ = + and asynchronous Ψ = − for the cross peak 1654 (α-helix) ↔ 1624 (β-sheet). Determine the order of change between these two groups during heating according to Noda's rules, and explain the molecular mechanism.
How can heterospectral 2D-COS (IR-Raman) aid the study of polymorphic drugs? Provide a specific scenario illustrating its advantages over single IR.
When performing 2D-COS with a set of temperature-perturbed spectra, the asynchronous spectrum is extremely noisy and almost uninterpretable. List at least 3 possible experimental reasons and suggest improvements.
References
[1] Noda I. "Two-Dimensional Infrared (2D-IR) Spectroscopy of Polystyrene." Appl Spectrosc, 1990, 44(4): 550–561. DOI:10.1366/0003702904087386.
[2] Noda I. "Generalized Two-Dimensional Correlation Method Applicable to Infrared, Raman, and Other Types of Spectroscopy." Appl Spectrosc, 1993, 47(9): 1329–1336. DOI:10.1366/0003702934067694.
[3] Noda I. "Progress in Two-Dimensional Correlation Spectroscopy." Vibrational Spectroscopy, 2004, 36(2): 143–165. DOI:10.1016/j.vibspec.2003.12.016.
[4] Morita S, Ozaki Y, Noda I. "Thirty Years of Progress in Two-Dimensional Correlation Spectroscopy." Appl Spectrosc, 2023, 77(1): 5–15. DOI:10.1177/00037028221123818.
[5] 2Dshige Software, Kwansei Gakuin University Ozaki Laboratory. https://sci-tech.kwansei.ac.j…
[6] Barth A. "Infrared Spectroscopy of Proteins." Biochimica et Biophysica Acta (BBA) - Bioenergetics, 2007, 1767(9): 1073–1101. DOI:10.1016/j.bbabio.2007.06.004.
[7] Fabian H, Mantsch H H. "Ribonuclease A Revisited: Infrared Spectroscopic Evidence for Lack of Native-like Secondary Structure in the Thermally Unfolded State." Biochemistry, 1995, 34(42): 13651–13655. DOI:10.1021/bi00041a039.
and see also: Fabian H, Mantsch H H. "Vibrational Spectroscopy of Proteins and Peptides." Vib Spectrosc, 1999, 19: 1–17.
[8] Cole K C, Pelletier M A, Noël D. "Infrared and Dynamic Mechanical Study of the Orientation in Poly(ethylene terephthalate) Films." Polymer Engineering & Science, 1999, 39(9): 1583–1592.
[9] Šašić S, Katsumoto Y, Sato H, Ozaki Y. "Applications of Moving Window Two-Dimensional Correlation Spectroscopy to Analysis of Phase-Separation Behavior of Polymer Blends." Anal Chem, 2003, 75(16): 4010–4018. DOI:10.1021/ac030086e.
[10] Noda I. "Two-Dimensional Correlation Spectroscopy Biomedical Applications." J Mol Struct, 2008, 883–884: 2–8. DOI:10.1016/j.molstruc.2007.11.022.
[11] Thomas M, Richardson H H. "Two-Dimensional FT-IR Correlation Analysis of the Temperature-Dependent IR Spectra of Poly(methyl methacrylate)." Vib Spectrosc, 2002, 30(1): 3–15. DOI:10.1016/S0924-2031(01)00166-0.
[12] ftir.fun amide functional group page. https://ftir.fun/ir/group/ami…
[13] ftir.fun carbonyl functional group page. https://ftir.fun/ir/group/car…
[14] ftir.fun ether bond functional group page. https://ftir.fun/ir/group/eth…
[15] py2D-COS: Python implementation of generalized 2D correlation spectroscopy. https://github.com/shigemasa-…
Next Episode Preview: Ep 42 — Hyphenated Techniques: TGA-FTIR, GC-FTIR, LC-FTIR
2D-COS is a powerful tool that expands "one-dimensional spectra" into "two-dimensional spectra," while hyphenated techniques expand "one spectroscopy" into multidimensional methods of "chromatography × spectroscopy." In the next episode, we will systematically introduce three main hyphenated approaches: thermogravimetric-infrared (TGA-FTIR), gas chromatography-infrared (GC-FTIR), and liquid chromatography-infrared (LC-FTIR). How are their interfaces designed? What samples are they suitable for? How to process three-dimensional data? We will explain with examples such as polymer thermal decomposition, fragrance component identification, and drug impurity analysis.
This article is licensed under CC BY-NC-SA 4.0. Images are from public domain or web resources with attributed sources, copyrights belong to the original authors.