§2 The Idea, Step by Step
From a jiggling spring to a molecular barcode
Picture every chemical bond as two balls joined by a tiny spring. Flick it and it jiggles — quickly if the spring is stiff and the balls are light, slowly if the spring is floppy and the balls are heavy. Molecules jiggle like this all the time, and infrared (IR) light is the "flick." Shine a spread of IR colours through a sample and each bond soaks up exactly the colour that matches its own natural jiggle rate. The colours that go missing form a barcode that tells you which bonds — and therefore which molecule — are present.
To put numbers on it we need two ideas. The force constant $k$ measures how stiff the bond is (a stiffer bond springs back harder), and the reduced mass $\mu = \dfrac{m_1 m_2}{m_1+m_2}$ is the single "effective mass" of the two atoms bouncing against each other. The jiggle shows up in the spectrum at a wavenumber $\tilde{\nu}$, measured in cm⁻¹ — literally how many wave crests fit in a centimetre, so a bigger number means higher-energy light. The simplest rules of thumb: a stiffer bond pushes $\tilde{\nu}$ up, and heavier atoms drag it down. As a worked feel for the numbers, a C–H bond ($k\approx480$ N/m, with a feather-light hydrogen) sits near 3000 cm⁻¹, while a C–Cl bond — heavier chlorine, floppier single bond — drops all the way to about 700 cm⁻¹.
The precise law comes from treating the bond as a quantum harmonic oscillator: $\tilde{\nu} = \dfrac{1}{2\pi c}\sqrt{\dfrac{k}{\mu}}$. Two consequences fall straight out. First, bond order sets $k$ (single < double < triple), so the carbonyl family climbs as the bond tightens: C–C near 1000, C=O near 1715, C≡N near 2200 cm⁻¹. Second, a vibration is only visible in IR if it changes the molecule's dipole moment, $\partial\mu/\partial Q \neq 0$ — which is why symmetric O₂ and N₂ are invisible to IR while a lopsided, polar C=O bond absorbs strongly. The control sliders map directly onto the formula: $k$ is bond stiffness, and $m_1$ and $m_2$ together fix $\mu$.
Try this in the sim above
(1) Drag k upward and watch the green marker on the spectrum slide toward higher wavenumber — stiffness wins. (2) Hold $k$ fixed and raise m₂ from 12 (carbon) toward 32 (sulfur): $\tilde{\nu}$ falls as $\mu$ grows. (3) Load the C=O preset, then the C–C preset, and watch the peak leap from the ~1000 cm⁻¹ region up to the carbonyl region — same kinds of atoms, just a stiffer double bond.
§3 Equation Derivation
Harmonic Oscillator — Vibrational Wavenumber
\[\tilde{\nu} = \frac{1}{2\pi c}\sqrt{\frac{k}{\mu}}\qquad \mu = \frac{m_1 m_2}{m_1+m_2}\]
| Symbol | Meaning | Unit |
|---|
| \(\tilde{\nu}\) | Wavenumber (cm⁻¹) | cm⁻¹ |
| \(k\) | Force constant (bond stiffness) | N m⁻¹ |
| \(\mu\) | Reduced mass = m₁m₂/(m₁+m₂) | kg |
| \(c\) | Speed of light = 3×10¹⁰ cm/s | cm s⁻¹ |
| \(E_v\) | Vibrational energy = hν(v+½) | J |
| \(v\) | Vibrational quantum number | 0,1,2,… |
| \(D_e\) | Dissociation energy | J mol⁻¹ |
Step-by-Step Derivation
Step 1 — Hooke's Law restoring force
For a diatomic bond stretched by x from equilibrium: F = −kx. Potential energy: V(x) = ½kx²
Step 2 — Reduce two-body to one-body problem
Instead of two masses m₁, m₂ oscillating, use reduced mass \(\mu = m_1m_2/(m_1+m_2)\) moving in V(x). This converts the two-body problem to a single harmonic oscillator.
Step 3 — Classical angular frequency
\(\omega = \sqrt{k/\mu}\) rad/s → frequency \(\nu = \omega/2\pi = \frac{1}{2\pi}\sqrt{k/\mu}\) Hz
Step 4 — Convert to wavenumber (IR convention)
\(\tilde{\nu} = \nu/c = \frac{1}{2\pi c}\sqrt{k/\mu}\) cm⁻¹. Using μ in kg and c = 2.998×10¹⁰ cm/s.
Step 5 — Quantum energy levels (Schrödinger for QHO)
\(E_v = h\nu(v+\tfrac{1}{2}) = hc\tilde{\nu}(v+\tfrac{1}{2})\), v = 0,1,2,…
Zero-point energy: \(E_0 = \tfrac{1}{2}hc\tilde{\nu} > 0\) even at v=0
Step 6 — Selection rule and IR activity
Δv = ±1 (harmonic oscillator). For IR activity: vibration must cause a change in dipole moment (∂μ/∂x ≠ 0). Symmetric stretches of homonuclear diatomics (N₂, O₂) are IR-inactive. CO₂ symmetric stretch is IR-inactive; asymmetric stretch is IR-active.
Worked Example — C=O stretch wavenumber
k(C=O) ≈ 1190 N/m, m_C = 12 amu, m_O = 16 amu
\(\mu = \dfrac{12\times16}{12+16}\times1.661\times10^{-27} = 1.139\times10^{-26}\) kg
\(\tilde{\nu} = \dfrac{1}{2\pi\times2.998\times10^{10}}\sqrt{\dfrac{1190}{1.139\times10^{-26}}} = 1715\) cm⁻¹ ✓
(The carbonyl double bond has k ≈ 1190 N/m. A stiffer value like ~1860 N/m belongs to the C≡O triple bond of carbon monoxide, which absorbs near 2143 cm⁻¹.)
Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., Chapter 12: "Vibrational Spectroscopy of Diatomic Molecules" | Skoog et al. — Fundamentals of Analytical Chemistry, 9th Ed., Chapter 16
§4 Frequently Asked Questions
🔬 SimulationWhat does each simulation mode show?▼
The Bond Vibration mode shows two atoms connected by a spring, oscillating at the frequency ν̃ = (1/2πc)√(k/μ), with the dipole moment arrow changing as the bond length changes. Stretching Modes shows symmetric and asymmetric stretch for a triatomic (like H₂O). Bending Modes shows scissoring, rocking, wagging, and twisting. IR Spectrum displays a realistic IR spectrum with characteristic group frequencies marked. Harmonic Potential shows the V = ½kx² well with quantized energy levels.
Key: Frequency depends on k (bond stiffness) and μ (masses) — stiffer bonds and lighter atoms vibrate at higher wavenumbers.
🌍 Real LifeWhere is IR spectroscopy used in real life?▼
IR spectroscopy is the primary tool for identifying functional groups in unknown organic compounds in pharmaceutical synthesis, polymer characterization, and forensic chemistry. In environmental science, atmospheric CO₂ is monitored by its IR absorption at 2349 cm⁻¹ — the basis of most climate sensors and greenhouse gas detectors. Breathalyzers use IR to detect ethanol in breath (C–H stretch at ~2950 cm⁻¹). Food quality control uses FTIR to detect adulteration (e.g., melamine in milk protein). In Bangladesh, BRTC pharmaceutical companies use FTIR as standard QC for API identification.
Key: IR is the molecular fingerprint — the region below 1500 cm⁻¹ is uniquely specific to each molecule and used for definitive identification.
🧪 ConceptualWhy must a vibration change the dipole moment to be IR-active?▼
IR absorption occurs when the oscillating electric field of the IR photon couples to an oscillating dipole moment in the molecule. For this coupling (resonance) to occur, the vibration must cause the dipole moment to change: ∂μ/∂x ≠ 0. If the dipole doesn't change during vibration, the molecule cannot "grab" the photon's energy. This is why N₂ and O₂ (symmetric, no permanent dipole) are IR-transparent — their stretches don't change the dipole. CO₂'s symmetric stretch doesn't change dipole (IR-inactive), but its bending and asymmetric stretch do (IR-active). This selectivity makes CO₂ a greenhouse gas but N₂ and O₂ transparent to IR.
Key: IR activity requires ∂μ/∂Q ≠ 0 — the bond must be polar AND its polarity must change during vibration.
🧮 MathematicalHow do you predict the wavenumber of a bond from k and masses?▼
Use ν̃ = (1/2πc)√(k/μ). Step 1: calculate μ = m₁m₂/(m₁+m₂) in kg (1 amu = 1.661×10⁻²⁷ kg). Step 2: compute √(k/μ) in rad/s. Step 3: divide by 2π×2.998×10¹⁰ to get cm⁻¹. Example: O–H bond, k ≈ 780 N/m, μ = (1×16)/(17)×1.661×10⁻²⁷ = 1.560×10⁻²⁷ kg. √(780/1.56×10⁻²⁷) = 7.06×10¹⁴ rad/s. ν̃ = 7.06×10¹⁴/(2π×3×10¹⁰) = 3750 cm⁻¹ ≈ observed O–H 3400–3650 cm⁻¹ (the harmonic estimate runs a few percent high) ✓
Key: The formula predicts wavenumbers to within ~10% — good enough to understand trends, but experimental values are needed for identification.
💡 Non-ObviousWhy does H/D substitution (deuteration) shift IR peaks?▼
Replacing H (1 amu) with D (2 amu) in a C–D or O–D bond changes μ but not k (same bond, same electrons). For C–H: μ = (1×12)/13 = 0.923 amu → ν̃(CH) ∝ √(1/0.923). For C–D: μ = (2×12)/14 = 1.714 amu → ν̃(CD) ∝ √(1/1.714). The ratio ν̃(CD)/ν̃(CH) = √(0.923/1.714) = √(0.538) = 0.734. So C–D appears at about 73% of C–H wavenumber: if C–H is at 3000 cm⁻¹, C–D appears at ~2200 cm⁻¹. This isotope effect is used to identify H-containing bonds and to study reaction mechanisms.
Key: Deuteration shifts X–H peaks by factor ~√(μ_H/μ_D) ≈ 0.73 — this is a powerful diagnostic tool in mechanistic chemistry.
🌌 DeepWhy is the real bond anharmonic, and what are the consequences?▼
The harmonic potential V=½kx² predicts equal energy spacing for all vibrational levels, but real bonds are described by the Morse potential V = De(1−e^{−ax})², which is asymmetric — it rises steeply for compression but allows bond breaking at large extension. The Morse potential gives unequally spaced levels (ΔE decreases as v increases) and allows overtone transitions (Δv = 2, 3) at roughly 2ν̃, 3ν̃ — not observable in the harmonic approximation. Anharmonicity also causes Fermi resonance (coupling of fundamentals near same energy) and combination bands, which are important for interpreting complex IR spectra.
Key: Real bonds are anharmonic — this is why overtones exist, why bonds break at finite energy, and why IR peaks are not perfectly predicted by the harmonic formula.
Reference: Skoog et al. — Fundamentals of Analytical Chemistry, 9th Ed., Chapter 16: "An Introduction to Infrared Spectrometry" | LibreTexts Chemistry — IR Spectroscopy https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Supplemental_Modules/Spectroscopy/Vibrational_Spectroscopy | Chemguide — Infrared Spectroscopy https://chemguide.co.uk/analysis/ir.html
§5 Common Misconceptions
❌ Misconception: "Every bond in a molecule gives a separate peak in the IR spectrum."
✅ Correction: Only IR-active vibrations give peaks (requires ∂μ/∂Q ≠ 0). In a molecule with N atoms, there are 3N−6 (nonlinear) or 3N−5 (linear) vibrational modes. Many modes are overlapping, degenerate, or so weak as to be invisible. Symmetric bonds (C–C in cyclohexane) give very weak peaks, while polar bonds (C=O, O–H, N–H) give very strong, characteristic peaks. Additionally, the fingerprint region (500–1500 cm⁻¹) contains many overlapping C–C, C–N, C–O stretches that are difficult to assign individually.
📖 Reference: Skoog et al. — Fundamentals of Analytical Chemistry, 9th Ed., §16B: "Theory of Infrared Spectroscopy"
❌ Misconception: "The C=O peak always appears at exactly 1715 cm⁻¹ regardless of the molecule."
✅ Correction: The C=O stretch appears in the range 1650–1850 cm⁻¹ depending on the molecular environment. Ketone C=O: ~1715 cm⁻¹; aldehyde: ~1720–1740 cm⁻¹; ester: ~1735–1750 cm⁻¹; amide: ~1650–1680 cm⁻¹ (lower due to resonance); acid chloride: ~1800 cm⁻¹; carboxylic acid: ~1700–1725 cm⁻¹. Conjugation (with C=C or aromatic rings) lowers the wavenumber by ~30–50 cm⁻¹ because the C=O bond order decreases. This sensitivity to environment is what makes IR useful for distinguishing functional groups.
📖 Reference: Clayden, Greeves & Warren — Organic Chemistry, 2nd Ed., Chapter 3: "Determining Organic Structures"
❌ Misconception: "CO₂ and N₂ don't absorb IR — so all symmetric molecules are IR transparent."
✅ Correction: N₂ and O₂ (homonuclear diatomics) have no dipole moment change during stretching — truly IR-inactive for their stretch. CO₂ (linear, symmetric) has its symmetric stretch IR-inactive, BUT its bending mode (doubly degenerate at ~667 cm⁻¹) and asymmetric stretch (~2349 cm⁻¹) ARE IR-active because they change the dipole. This is why CO₂ is a greenhouse gas despite being "symmetric." A molecule need not be entirely asymmetric to have IR-active modes — it only needs some vibrational modes that change the dipole.
📖 Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., §12D.1: "The vibrational selection rule"
❌ Misconception: "Heavier atoms always vibrate at lower wavenumber — so C–Cl is always lower than C–H."
✅ Correction: Generally true (μ is larger for heavier atoms, so ν̃ is lower), but force constant k also matters. C–Cl (single bond, k ≈ 300 N/m) appears at ~700 cm⁻¹. But C=O (double bond, k ≈ 1860 N/m, μ ≈ 11.4×10⁻²⁷ kg) appears at 1715 cm⁻¹ — higher than C–Cl despite heavier oxygen, because the larger k dominates. C≡N (triple bond, k ≈ 2400 N/m) appears near 2200 cm⁻¹. Bond order (which determines k) often outweighs mass effects.
📖 Reference: Skoog et al. — Fundamentals of Analytical Chemistry, 9th Ed., Table 16.2: "Force constants for various bonds"
❌ Misconception: "The fingerprint region (500–1500 cm⁻¹) is useless because it is too complex to interpret."
✅ Correction: The fingerprint region is actually the most useful part for definitive molecular identification. While individual peaks are hard to assign, the entire pattern is unique to each molecule — like a fingerprint. Comparing the fingerprint region of an unknown to a reference spectrum in databases (SDBS, NIST) provides unambiguous identification even for isomers that may have identical functional group regions. The fingerprint region is also where skeletal C–C stretches, C–O–C vibrations, and ring deformations appear.
📖 Reference: Clayden, Greeves & Warren — Organic Chemistry, 2nd Ed., Chapter 3: "The fingerprint region"
Section 4 reference: Taber, K.S. — Chemical Misconceptions (RSC, 2002) | Henary, M. & Russell, A.A. — J. Chem. Educ. 2007, 84, 1442 (IR misconceptions in undergraduate courses)