14.12.03 · genchem-pchem / spectroscopy

Vibrational spectroscopy: harmonic oscillator selection rules, anharmonicity, and IR-active modes

shipped3 tiersLean: nonepending prereqs

Anchor (Master): Herzberg — Infrared and Raman Spectra (1945)

Intuition Beginner

Chemical bonds act like springs. When a bond stretches and compresses, it vibrates at a frequency that depends on two things: how stiff the bond is (the force constant) and how heavy the atoms are (the reduced mass). A stiff, light bond like H-F vibrates fast — its stretching frequency is about 4000 inverse centimetres. A softer, heavier bond like C-Cl vibrates slowly — about 700 inverse centimetres. Infrared light carries exactly the right energies to set bonds vibrating, so when IR light passes through a sample, the frequencies that are absorbed map directly onto the bond-stretching and bond-bending motions present in the molecule.

This is why IR spectroscopy is the chemist's workhorse for identifying functional groups. The C=O stretch always absorbs near 1700 inverse centimetres, whether the carbonyl is in an aldehyde, a ketone, an ester, or a carboxylic acid. The O-H stretch always appears near 3400 inverse centimetres, broadened by hydrogen bonding. A C-H stretch near 3000 inverse centimetres signals the presence of hydrogen atoms bonded to carbon. By reading the positions, shapes, and intensities of the absorption bands, a chemist identifies which bonds and functional groups are present — often within minutes and with a sample of less than a milligram.

Not every vibration absorbs IR light. The selection rule is simple: a vibration is IR-active only if it changes the molecular dipole moment. A homonuclear diatomic like H or N has no dipole moment and stretching the bond does not create one, so these molecules are IR-silent. A heteronuclear diatomic like HCl has a dipole that changes as the bond stretches, so HCl absorbs IR strongly. For polyatomics, the criterion is applied mode by mode: each of the normal modes of a molecule with atoms either changes the dipole moment (IR-active) or does not (IR-silent).

Visual Beginner

A diatomic molecule vibrating along its bond axis traces out a potential-energy curve. The harmonic approximation treats this curve as a parabola, with equally spaced vibrational energy levels like the rungs of a ladder. The spacing between rungs is the vibrational frequency, set by the force constant and reduced mass : .

The harmonic model predicts that only the transition (the fundamental) is allowed. Real molecules are anharmonic — the potential is not a perfect parabola — and this allows weak overtone transitions at roughly , , and so on. The fundamental is always the strongest absorption. The anharmonicity also means the energy levels are not exactly equally spaced: the spacing decreases slightly with increasing , converging to the dissociation limit.

Worked example Beginner

The C-O stretching frequency in carbon monoxide is observed at 2143 inverse centimetres. The reduced mass of CO is kg. Calculate the force constant of the C-O bond and classify the bond order.

Step 1. The vibrational frequency in wavenumber units is . Rearranging for :

N/m.

Step 2. Classify by bond order. Single bonds (C-C, C-O) have force constants around 500 N/m. Double bonds (C=C, C=O) around 1000 N/m. Triple bonds (CC, CO) around 1500-2000 N/m. The value of 1902 N/m confirms that CO has a triple-bond character, consistent with its Lewis structure O: with a formal bond order of 3.

Step 3. The force constant is a measure of bond stiffness. A larger means a stronger, shorter bond that vibrates at higher frequency. This is why triple bonds absorb at higher wavenumbers than double bonds, and double bonds higher than single bonds — the basis of the group-frequency approach to IR spectroscopy.

Check your understanding Beginner

Formal definition Intermediate+

The vibrational motion of a diatomic molecule is modelled by the quantum harmonic oscillator. The Hamiltonian for a diatomic with reduced mass and force constant is

where is the displacement from equilibrium. The energy eigenvalues are

The zero-point energy means the molecule never stops vibrating. In wavenumber units the level spacing is .

Definition (harmonic-oscillator selection rule). The electric-dipole selection rule for the harmonic oscillator is . This follows from the matrix element , which is nonzero only for . The transition is IR-active only if the dipole moment changes during the vibration, i.e., .

Definition (Morse potential and anharmonicity). The Morse potential provides a more realistic model where the potential is asymmetric: steeper than harmonic at short range and flattening toward dissociation at long range. The vibrational energy levels are

where and is the anharmonicity constant. The negative quadratic term causes the level spacing to decrease with increasing , and the number of bound states is finite: .

Definition (IR activity for polyatomic molecules). A molecule with atoms has normal modes ( for linear molecules). Each normal mode is IR-active if the dipole moment changes along that coordinate: . Group theory determines activity from symmetry: a mode is IR-active if it transforms as the same irreducible representation as one of the translation components (, , or ) of the molecular point group.

Counterexamples to common slips

  • Fundamental versus overtone intensity. The selection rule is exact only for the harmonic oscillator. Anharmonicity relaxes this: overtones () are observed but are typically 10 to 100 times weaker than the fundamental. Treating the observed fundamental as "the only band" misses the overtone and combination bands that are diagnostic in gas-phase and matrix-isolation IR spectroscopy.

  • Number of IR-active modes. A molecule with normal modes does not necessarily show IR peaks. Modes may be degenerate (e.g., the two bending modes of CO are degenerate at 667 cm), IR-inactive by symmetry, or too weak to observe. The number of observed IR bands is typically less than .

  • Force constant versus bond energy. A large force constant means a stiff bond, not necessarily a strong one. The C-H force constant (about 500 N/m) is similar to the C-C force constant, but the C-H bond energy (413 kJ/mol) is much larger than the C-C bond energy (348 kJ/mol). Force constants and bond dissociation energies probe different aspects of the potential-energy surface.

Key result Intermediate+

Theorem (normal-mode decomposition and spectroscopic activity). Let a molecule with atoms have a potential energy expanded about the equilibrium geometry as , where are mass-weighted displacement coordinates and is the mass-weighted Hessian. Diagonalising gives eigenvalues and eigenvectors (the normal modes). The translational and rotational eigenvalues are zero; the remaining (or ) positive eigenvalues give the normal-mode frequencies. A mode is IR-active if and only if the transition dipole moment is nonzero, which is equivalent to .

Proof. The molecular potential energy surface near equilibrium is a quadratic form in the nuclear displacement coordinates. Removing the six (or five) zero-frequency translational and rotational degrees of freedom by projecting onto the vibrational subspace leaves a symmetric positive-definite matrix. The spectral theorem guarantees a real orthogonal diagonalisation , where with . The normal-mode coordinates are , and the Hamiltonian separates into independent harmonic oscillators: . The dipole operator expanded to first order in the normal coordinates is . The first term (, the permanent dipole) does not couple different vibrational states. The second term gives transition dipole moments , which are nonzero if and only if . This is the IR-activity criterion.

Bridge. This result connects the purely mathematical diagonalisation of the Hessian matrix to the observable IR spectrum. The normal-mode analysis predicts how many IR-active bands a molecule should have and at what frequencies. The connection to group theory (14.02.01) is that the irreducible representations of the molecular point group determine which normal modes have nonzero , providing a rapid symmetry-based prediction of the IR spectrum without explicit matrix diagonalisation.

Worked example at intermediate level

Determine the number of IR-active normal modes for formaldehyde (HCO, point group ) and predict the symmetry species of each mode.

Formaldehyde has 4 atoms, so normal modes. The reducible representation of the displacements decomposes under as . Subtracting translations ( for ) and rotations ( for ) gives the vibrational modes:

In , the translation components , , transform as , , respectively. All four symmetry species of the vibrations (, , , ) match at least one translation component except . Therefore: five IR-active modes () and one IR-inactive mode (). The mode (a twisting motion) is Raman-active but IR-silent. Experimentally, formaldehyde shows five strong IR absorption bands, consistent with this prediction.

Exercises Intermediate+

The Morse potential, dissociation, and vibrational thermodynamics Master

The harmonic oscillator is the starting point for understanding molecular vibrations, but real bonds are anharmonic. The Morse potential captures the essential physics: the potential is approximately quadratic near equilibrium, but it flattens toward the dissociation energy as the bond is stretched, and it rises steeply for compressed bond lengths. The Schrodinger equation for the Morse potential can be solved exactly, yielding the energy levels

where is the harmonic frequency and is the anharmonicity. The parameter sets the width of the potential well. The Morse model has three experimentally determined parameters (, , and ) compared with two for the harmonic oscillator ( and ), and the extra parameter quantifies the departure from harmonicity.

The Morse potential predicts three observable phenomena absent from the harmonic model. First, the decreasing level spacing produces a convergence of overtone frequencies: the fundamental is at , the first overtone is at , the second overtone is at , and so on. The overtones are not exact integer multiples of the fundamental — each successive overtone is slightly shifted to lower frequency by the cumulative anharmonic correction. Measuring the positions of several overtones determines both and simultaneously.

Second, the dissociation energy can be extracted from the spectroscopic constants. Setting and solving gives in wavenumber units. For H, cm and , giving cm = 4.95 eV, in excellent agreement with the directly measured dissociation energy of 4.48 eV (the small difference arises from the difference between and ). The Birge-Sponer extrapolation — plotting against and extrapolating the straight line to — is the standard graphical method for extracting from vibrational spectra.

Third, the finite number of bound states means that at sufficiently high excitation the molecule dissociates. This connects vibrational spectroscopy to chemical kinetics: photodissociation thresholds correspond to reaching , and the dissociation products carry the excess energy as translational kinetic energy. Predissociation — where the vibrational level lies above a dissociation asymptote but the molecule is temporarily trapped behind a small barrier — produces line broadening in the vibrational spectrum, providing a spectroscopic window onto dissociation dynamics.

Group frequencies and the fingerprint region Master

The normal-mode analysis shows that each molecule has a unique set of vibrational frequencies. In practice, however, certain vibrational frequencies are remarkably transferable across different molecules: the C=O stretch appears near 1700 cm whether the carbonyl is in formaldehyde, acetone, acetic acid, or a peptide bond. This transferability arises because the C=O stretching mode is largely localised to the C=O group — the normal-mode displacement pattern concentrates most of the motion on the carbonyl carbon and oxygen, with only small participation from the rest of the molecule.

This localisation is the basis of group-frequency analysis. The IR spectrum is divided into two regions. The functional-group region (4000-1500 cm) contains stretching vibrations that are diagnostic of specific bonds: O-H and N-H stretches (3200-3600 cm), C-H stretches (2800-3100 cm), CN and CC stretches (2100-2260 cm), C=O stretches (1650-1800 cm), C=C stretches (1500-1680 cm). The fingerprint region (below 1500 cm) contains C-C, C-O, C-N, and C-X single-bond stretches, bending modes, and coupled motions that are specific to the entire molecular skeleton. Two molecules with the same functional groups but different connectivity show the same bands in the functional-group region but different patterns in the fingerprint region.

The positions of group frequencies shift systematically with the electronic environment. A carbonyl adjacent to an electron-withdrawing group is strengthened (the C=O bond has more double-bond character), shifting the stretch to higher wavenumber: acid chlorides (1810 cm), esters (1735 cm), ketones (1715 cm), amides (1680 cm). Conjugation with a C=C double bond weakens the C=O by delocalising the pi electrons, shifting the stretch to lower wavenumber: conjugated ketones absorb near 1685 cm. Hydrogen bonding to the carbonyl oxygen further weakens the bond, producing additional red shifts and peak broadening. These systematic trends, first compiled by Colthup [Colthup1950], form the basis of modern IR spectral interpretation.

Raman spectroscopy: the complementary technique Master

Raman spectroscopy probes molecular vibrations through inelastic light scattering rather than absorption. When monochromatic light of frequency illuminates a sample, most scattered photons have the same frequency (Rayleigh scattering), but a small fraction ( to ) are shifted in frequency by amounts equal to molecular vibrational frequencies. Stokes-shifted photons (frequency ) correspond to vibrational excitation; anti-Stokes photons (frequency ) correspond to vibrational de-excitation. The Stokes lines dominate at room temperature because most molecules are in the vibrational ground state.

The selection rule for Raman activity is different from IR activity: a vibration is Raman-active if the molecular polarisability changes during the vibration (). For centrosymmetric molecules, the mutual-exclusion principle applies: modes that are IR-active are Raman-inactive and vice versa. Combining IR and Raman spectroscopy of a centrosymmetric molecule recovers the complete vibrational spectrum. Even for non-centrosymmetric molecules, the two techniques have different sensitivities: symmetric stretches (which strongly modulate the polarisability) are Raman-strong but often IR-weak, while asymmetric stretches (which strongly modulate the dipole moment) are IR-strong but Raman-weak.

The Raman effect was predicted theoretically by Smekal in 1923 and discovered experimentally by Raman and Krishnan in 1928 [Raman1928], earning Raman the 1930 Nobel Prize in Physics. Modern Raman spectroscopy uses laser excitation and CCD detection to record spectra with high signal-to-noise. Resonance Raman spectroscopy — where the laser frequency matches an electronic transition — enhances Raman scattering by factors of to for modes coupled to the chromophore, providing selective enhancement of specific vibrational modes. Surface-enhanced Raman scattering (SERS) on rough metal surfaces amplifies signals by to , enabling detection of single molecules. Coherent anti-Stokes Raman scattering (CARS) microscopy provides label-free vibrational imaging of biological tissue with subcellular resolution.

Connections Master

  • Rotational spectroscopy 14.12.02. Each vibrational transition carries a rotational fine structure: the P-branch () and R-branch () bracket the vibrational band origin, and the Q-branch () appears when the electronic or vibrational angular momentum permits. The rotational constants and of the ground and excited vibrational states are extracted from the P/R branch line positions, connecting the vibrational and rotational analyses.

  • Electronic spectroscopy 14.12.04. Electronic transitions produce vibronic bands whose spacing equals the vibrational frequency of the excited electronic state. The Franck-Condon principle governs which vibrational levels of the excited state are populated, and the resulting intensity pattern encodes the displacement between ground and excited-state potential-energy surfaces. Understanding vibrational spectroscopy is a prerequisite for interpreting the vibronic structure of UV-Vis spectra.

  • Statistical mechanics 14.07.02. The vibrational partition function per mode determines the vibrational contributions to heat capacity, entropy, and Gibbs free energy. The vibrational frequencies measured by IR spectroscopy feed directly into thermodynamic calculations via the partition function. At low temperature, (all molecules in ); at high temperature, (classical limit).

  • Molecular orbital theory 14.05.02. The force constants of bonds are predicted by MO theory: the bond order determines the electron density in the bonding orbital, which sets the restoring force against displacement. The correlation between bond order and vibrational frequency (single < double < triple) is a direct spectroscopic consequence of the MO picture.

  • Organic functional groups 15.02.01. The IR group-frequency tables used for structure determination map directly onto the functional-group classification of organic chemistry. The carbonyl stretch, the O-H stretch, the C-H stretch — each is a spectroscopic signature of a specific functional group from the organic chemistry taxonomy.

Historical notes Master

The study of molecular vibrations through infrared spectroscopy has its origins in the discovery of infrared radiation itself. William Herschel detected infrared light in 1800 by observing that a thermometer placed beyond the red end of a solar spectrum registered heating. The first IR absorption spectra of molecules were recorded in the 1880s by Abney and Festing, who observed that specific organic functional groups absorbed at characteristic wavelengths. Their work established the qualitative connection between molecular structure and IR absorption that underpins modern spectral interpretation.

The theoretical framework developed over the following decades. The quantum harmonic oscillator was solved by Schrodinger in 1926, providing the energy-level formula and the selection rule . The Morse potential was introduced by Philip Morse in 1929 [Morse1929] as a model for diatomic vibrations that could be solved exactly while capturing anharmonicity. Herzberg's monumental three-volume work Molecular Spectra and Molecular Structure (1939-1966) [Herzberg1945] systematised the entire field of molecular spectroscopy, including the vibrational spectra of hundreds of molecules, and earned him the 1971 Nobel Prize in Chemistry.

The development of IR instrumentation transformed the technique from a specialist research tool into a routine analytical method. The first commercial double-beam IR spectrometer (Perkin-Elmer Model 21, 1944) made IR spectroscopy accessible to industrial and academic chemistry laboratories. The introduction of Fourier-transform IR (FT-IR) spectrometers in the 1970s — replacing the slow mechanical scanning of dispersive instruments with rapid interferometric measurement — improved sensitivity by orders of magnitude (the Jacquinot and Fellgett advantages) and reduced scan times from hours to seconds.

Raman spectroscopy had a slower development trajectory because the effect is inherently weak. The invention of the laser in 1960 provided the intense monochromatic source that Raman spectroscopy required, and by the 1970s laser Raman spectroscopy was routine. The development of notch filters (to reject the intense Rayleigh line), holographic gratings, and CCD detectors in the 1980s and 1990s brought Raman spectroscopy to its current level of sensitivity and convenience.

Colthup's 1950 correlation tables [Colthup1950] — which mapped functional groups to characteristic IR frequency ranges — transformed IR spectroscopy from a research technique into a universal analytical method for organic chemists. The tables have been refined and extended many times since, but the basic framework remains: each functional group has a diagnostic IR signature, and reading the spectrum is reading the molecular structure.

Bibliography Master

@book{Atkins2022,
  author = {Atkins, P. and de Paula, J. and Keeler, J.},
  title = {Physical Chemistry},
  edition = {12},
  publisher = {Oxford University Press},
  year = {2022},
}

@book{Banwell1994,
  author = {Banwell, C. N. and McCash, E. M.},
  title = {Fundamentals of Molecular Spectroscopy},
  edition = {4},
  publisher = {McGraw-Hill},
  year = {1994},
}

@book{Herzberg1945,
  author = {Herzberg, G.},
  title = {Molecular Spectra and Molecular Structure II: Infrared and Raman Spectra of Polyatomic Molecules},
  publisher = {Van Nostrand},
  year = {1945},
}

@book{Pavia2015,
  author = {Pavia, D. L. and Lampman, G. M. and Kriz, G. S. and Vyvyan, J. R.},
  title = {Introduction to Spectroscopy},
  edition = {5},
  publisher = {Cengage},
  year = {2015},
}

@article{Morse1929,
  author = {Morse, P. M.},
  title = {Diatomic Molecules According to the Wave Mechanics II: Vibrational Levels},
  journal = {Physical Review},
  volume = {34},
  year = {1929},
  pages = {57--64},
}

@article{Raman1928,
  author = {Raman, C. V. and Krishnan, K. S.},
  title = {A New Type of Secondary Radiation},
  journal = {Nature},
  volume = {121},
  year = {1928},
  pages = {501--502},
}

@article{Colthup1950,
  author = {Colthup, N. B.},
  title = {Spectra-Structure Correlations in the Infrared Region},
  journal = {Journal of the Optical Society of America},
  volume = {40},
  year = {1950},
  pages = {397--400},
}