Molecular Physics, Chemical Bonding & Molecular Spectroscopy
Comprehensive theory of molecular structure and spectroscopy: nature of chemical bonds (ionic, covalent, van der Waals); quantum theory of H2+ and H2 via the LCAO molecular orbital method; Born-Oppenheimer separation of electronic, vibrational, and rotational motions; pure rotational spectroscopy of rigid rotors and bond length determinations; vibrational spectroscopy, Morse potential anharmonicity, and vibration-rotation P/R branches; electronic spectra, Franck-Condon principle, and classical/quantum Raman scattering.
§6.1 Molecular Chemical Bonds & The Molecular Orbital Concept
1. Classification of Molecular Chemical Bonds
Molecules are stable aggregates of two or more atoms bound together by electromagnetic interactions. Chemical bonds are categorized by their physical bonding mechanisms:
- Ionic Bonds: Formed by complete electrostatic electron transfer between atoms of widely differing electronegativities (e.g., $\text{NaCl}, \text{KBr}$). The cohesive energy is governed by Coulomb attraction balanced by Pauli short-range electron core repulsion:
$$U(R) = - \frac{e^2}{4\pi\varepsilon_0 R} + \frac{A}{R^n} \quad (n \approx 8\text{ to }10)$$
- Covalent Bonds: Formed between atoms of similar electronegativities by the quantum mechanical sharing of valence electrons (e.g., $\text{H}_2, \text{O}_2, \text{CH}_4$). The shared electron density accumulates in the internuclear region, electrostatically screening the positive nuclei from mutual repulsion and lowering total quantum energy.
- Van der Waals Bonds: Weak intermolecular bonds ($0.01\text{ to }0.1\,\text{eV}$) arising from fluctuating dipole-induced dipole electrostatic attractions ($U(R) \propto -1/R^6$).
- Hydrogen Bonds: Intermediate dipole-dipole attractions ($0.1\text{ to }0.5\,\text{eV}$) formed between an electropositive hydrogen atom covalently bound to an electronegative atom (N, O, F) and an adjacent lone pair.
2. The Molecular Orbital Concept & LCAO Approximation
In molecular orbital theory, electrons are not confined to individual atoms, but occupy delocalized Molecular Orbitals (MOs) extending across the entire molecule.
The simplest mathematical approximation is the Linear Combination of Atomic Orbitals (LCAO). For a diatomic molecule with nuclei $A$ and $B$, a molecular orbital $\psi_{MO}$ is constructed from atomic orbitals $\phi_A$ and $\phi_B$:
For homonuclear diatomics ($c_A = \pm c_B$), this produces two distinct spatial distributions:
- Bonding Molecular Orbital ($\sigma_g$): Symmetric linear combination $\psi_+ = N_+ (\phi_A + \phi_B)$. Electron density builds up constructively between the nuclei, lowering the electrostatic potential energy and forming a stable chemical bond.
- Antibonding Molecular Orbital ($\sigma_u^*$): Antisymmetric linear combination $\psi_- = N_- (\phi_A - \phi_B)$. A nodal plane where $\psi = 0$ exists midway between the nuclei. Electron density is pushed away from the internuclear region, producing net nuclear repulsion.
§6.2 Quantum Mechanics of the H2+ Ion & The Hydrogen Molecule (H2)
1. The Hydrogen Molecular Ion ($H_2^+$)
The simplest molecular system in nature is the hydrogen molecular ion $H_2^+$, consisting of two positive protons ($A$ and $B$) separated by internuclear distance $R$, and a single electron.
The electronic Hamiltonian in atomic units ($\hbar = m_e = e = 1$) is:
Using the LCAO trial functions formed from hydrogen $1s$ orbitals $\phi_A$ and $\phi_B$:
where the Overlap Integral $S$ is:
The electronic expectation energy $E_{\pm}(R) = \langle \psi_{\pm} | \hat{H} | \psi_{\pm} \rangle$ evaluates to:
where $J$ is the Coulomb Integral (classical electrostatic interaction between the charge cloud around nucleus $A$ and nucleus $B$) and $K$ is the Exchange (Resonance) Integral (representing quantum electron tunneling between the two protons):
Because $K$ is negative, the symmetric state $E_+(R)$ develops a pronounced potential energy minimum at equilibrium bond length $R_e = 1.06\,\text{\AA}$ with a binding dissociation energy $D_e = 2.79\,\text{eV}$.
2. The Neutral Hydrogen Molecule ($H_2$) & Heitler-London Theory
In the neutral hydrogen molecule ($H_2$), there are two electrons. By the Pauli exclusion principle, the total two-electron wavefunction must be antisymmetric under exchange of electrons $1$ and $2$:
- Singlet Ground State ($S = 0$, Antiparallel Spins): The spin state is antisymmetric $\chi_{singlet} = \frac{1}{\sqrt{2}}(\alpha_1 \beta_2 - \beta_1 \alpha_2)$. Consequently, the spatial wavefunction must be symmetric:
$$\psi_S(\vec{r}_1, \vec{r}_2) = \frac{1}{\sqrt{2(1 + S^2)}} [ \phi_A(1) \phi_B(2) + \phi_B(1) \phi_A(2) ]$$This generates a deep potential well with equilibrium bond length $R_e = 0.74\,\text{\AA}$ and strong covalent dissociation energy $D_e = 4.75\,\text{eV}$.
- Triplet Excited State ($S = 1$, Parallel Spins): The spin state is symmetric, requiring an antisymmetric spatial function $\psi_T \propto [\phi_A(1)\phi_B(2) - \phi_B(1)\phi_A(2)]$. Spatial electron density vanishes between the nuclei, producing purely repulsive forces for all $R$.
§6.3 Born-Oppenheimer Approximation & Molecular Energy Hierarchy
1. The Born-Oppenheimer Approximation (1927)
Max Born and J. Robert Oppenheimer recognized that atomic nuclei are vastly more massive than electrons ($M_{nucleus} / m_e \sim 1836\text{ to }10^5$). Consequently, electrons move at velocities hundreds of times faster than nuclei.
On the timescale of electronic orbital motion, the heavy nuclei can be treated as essentially stationary fixed points in space. The full molecular Hamiltonian separates into:
- An Electronic Schrödinger Equation solved at fixed nuclear configurations $R$, generating potential energy curves $V_{el}(R)$.
- A Nuclear Schrödinger Equation governing the vibrational and rotational motions of the nuclei moving on the effective potential energy surface $V_{el}(R)$.
2. Molecular Energy Hierarchy
To excellent approximation, the total internal energy of a molecule is the sum of three independent contributions:
| Motion Type | Energy Order of Magnitude | Spectral Region | Typical Wavelength $\lambda$ |
|---|---|---|---|
| Electronic Transitions | $1\text{ to }10\,\text{eV}$ | Visible / Ultraviolet | $100\text{ to }700\,\text{nm}$ |
| Vibrational Transitions | $0.05\text{ to }0.5\,\text{eV}$ | Infrared (IR) | $2\text{ to }20\,\mu\text{m}$ |
| Rotational Transitions | $10^{-4}\text{ to }10^{-2}\,\text{eV}$ | Microwave / Far-IR | $0.1\text{ to }10\,\text{mm}$ |
Because $\Delta E_{rot} \ll \Delta E_{vib} \ll \Delta E_{el}$, molecular spectra exhibit fine structure: electronic bands contain closely spaced vibrational progressions, which in turn contain ultra-dense rotational lines.
§6.4 Pure Rotational Spectroscopy: Rigid Rotor Model & Bond Length Determinations
1. The Rigid Rotor Model of a Diatomic Molecule
Consider a diatomic molecule consisting of two masses $m_1$ and $m_2$ separated by a fixed equilibrium bond length $r_0$. The classical moment of inertia about the center of mass axis is:
The classical kinetic energy of rotation is $E = \frac{L^2}{2 I}$. In quantum mechanics, orbital angular momentum is quantized: $L^2 = J(J+1) \hbar^2$, where $J = 0, 1, 2, 3, \dots$ is the Rotational Quantum Number.
The quantized rotational energy levels are:
where the Rotational Constant $B$ (expressed in wavenumber units $\text{cm}^{-1}$) is:
2. Selection Rules & Rotational Spectra
For an electric dipole transition to occur in pure rotational spectroscopy:
- Gross Selection Rule: The molecule must possess a permanent electric dipole moment ($\mu_{el} \neq 0$). Homonuclear molecules ($\text{H}_2, \text{N}_2, \text{O}_2$) have zero dipole moment and are completely microwave inactive. Heteronuclear molecules ($\text{CO}, \text{HCl}, \text{NO}$) have permanent dipoles and exhibit intense rotational absorption.
- Specific Selection Rule:
$$\Delta J = \pm 1 \quad (+1 \text{ for absorption, } -1 \text{ for emission})$$
The wavenumber of the transition from level $J$ to $J + 1$ is:
Evaluating for consecutive transitions:
- $J = 0 \to 1$: $\bar{\nu} = 2B$
- $J = 1 \to 2$: $\bar{\nu} = 4B$
- $J = 2 \to 3$: $\bar{\nu} = 6B$
- $J = 3 \to 4$: $\bar{\nu} = 8B$
Fundamental Experimental Signature: The pure rotational absorption spectrum consists of a series of equidistant spectral lines separated by constant spacing $2B$:
By measuring this line separation $\Delta \bar{ u} = 2B$ in the microwave laboratory, one calculates the moment of inertia $I = \frac{h}{8\pi^2 c B}$ and determines the internuclear bond distance $r_0 = \sqrt{I / \mu}$ to four decimal places of precision!
§6.5 Vibrational & Vibration-Rotation Spectra: P-Branch & R-Branch Transitions
1. Harmonic vs. Anharmonic Morse Potential
Near equilibrium separation $r_0$, the molecular potential energy can be approximated as a simple harmonic oscillator with bond force constant $k$:
The quantized vibrational energy levels are:
where $\bar{\nu}_0 = \frac{1}{2\pi c} \sqrt{\frac{k}{\mu}}$. Even in the ground state ($v = 0$), the molecule possesses irreducible zero-point vibrational energy $E_0 = \frac{1}{2} \hbar \omega_0$.
Real chemical bonds dissociate at large separations. A far more realistic model is the Morse Potential:
where $D_e$ is the depth of the potential well. The energy levels of an anharmonic Morse oscillator are:
where $x_e$ is the anharmonicity constant. Anharmonicity relaxes the strict harmonic selection rule $\Delta v = \pm 1$, permitting weaker overtone transitions ($\Delta v = \pm 2, \pm 3$).
2. Vibration-Rotation Spectra & P/R Branch Architecture
Because rotational energy levels are densely packed within each vibrational state, a vibrational transition ($v = 0 \to 1$) is always accompanied by simultaneous rotational transitions ($J \to J'$). The combined energy of a vibration-rotation state is:
For heteronuclear diatomic molecules with zero electronic angular momentum ($\Sigma$ states), the selection rules are:
The resulting spectrum splits into two distinct symmetric branches flanking the missing fundamental vibrational frequency $\bar{ u}_0$:
- The $R$-Branch ($\Delta J = +1$, $J' = J + 1$): Rotational energy increases during vibrational absorption:
$$\bar{\nu}_R(J) = \bar{\nu}_0 + B(J+1)(J+2) - B J(J+1) = \bar{\nu}_0 + 2 B (J + 1) \quad (J = 0, 1, 2, \dots)$$Lines appear at higher frequencies: $\bar{\nu}_0 + 2B, \bar{\nu}_0 + 4B, \bar{\nu}_0 + 6B, \dots$
- The $P$-Branch ($\Delta J = -1$, $J' = J - 1$): Rotational energy decreases:
$$\bar{\nu}_P(J) = \bar{\nu}_0 + B(J-1)J - B J(J+1) = \bar{\nu}_0 - 2 B J \quad (J = 1, 2, 3, \dots)$$Lines appear at lower frequencies: $\bar{\nu}_0 - 2B, \bar{\nu}_0 - 4B, \bar{\nu}_0 - 6B, \dots$
- The Missing $Q$-Branch ($\Delta J = 0$): A line at $\bar{ u} = \bar{ u}_0$ would correspond to $\Delta J = 0$. Because $\Delta J = 0$ is forbidden, there is a distinct central gap of width $4B$ at the fundamental origin!
§6.6 Electronic Spectra, The Franck-Condon Principle & The Raman Effect
1. Electronic Transitions & The Franck-Condon Principle
Electronic transitions involve major rearrangements of electron clouds, shifting the equilibrium internuclear distance from $r_0$ to $r_0'$.
The Franck-Condon Principle states: Because atomic nuclei are vastly heavier than electrons, an electronic transition occurs so rapidly ($\sim 10^{-15}\text{ s}$) that the nuclei do not have time to change their positions or momenta during the transition.
On a potential energy diagram, electronic transitions are represented as strictly vertical lines. The transition probability (intensity) between vibrational state $v$ in the ground electronic state and $v'$ in the excited state is proportional to the Franck-Condon Factor—the square of the vibrational overlap integral:
2. The Raman Effect: Classical Polarizability Model (1928)
Sir C.V. Raman discovered that when a transparent substance is irradiated with intense monochromatic light of frequency $ u_0$, a small fraction ($\sim 10^{-6}$) of the scattered radiation emerges with altered frequencies ($ u_0 \pm u_v$).
Classically, the incident electric field $\vec{E}(t) = \vec{E}_0 \cos(2\pi u_0 t)$ induces an electric dipole moment in the molecule:
where $\alpha$ is the molecular polarizability. As the molecule vibrates with natural frequency $ u_v$, its polarizability fluctuates periodically:
Substituting $\alpha(t)$ into the induced dipole equation:
The oscillating dipole radiates electromagnetic waves at three distinct frequencies:
- Rayleigh Scattering ($\nu_0$): Elastic scattering at the unshifted incident frequency.
- Stokes Lines ($\nu_0 - \nu_v$): Inelastic scattering red-shifted to lower frequency. The incident photon transfers energy to excite a molecular vibration.
- Anti-Stokes Lines ($\nu_0 + \nu_v$): Inelastic scattering blue-shifted to higher frequency. The incident photon absorbs energy from an already-vibrating molecule.
3. Quantum Interpretation & The Rule of Mutual Exclusion
In quantum theory, an incident photon of energy $h u_0$ promotes the molecule to a transient virtual state:
- If it de-excites to a higher vibrational level ($v = 0 \to 1$), the scattered photon carries reduced energy $h(\nu_0 - \nu_v)$ (Stokes).
- If it de-excites from an initial excited level to the ground level ($v = 1 \to 0$), the scattered photon carries increased energy $h(\nu_0 + \nu_v)$ (Anti-Stokes).
Because the thermal population of the excited state $v = 1$ is governed by the Boltzmann distribution $N_1 / N_0 = e^{-h u_v / k_B T} \ll 1$, Stokes lines are always vastly more intense than Anti-Stokes lines:
The Rule of Mutual Exclusion: For molecules with a center of inversion symmetry (such as $\text{CO}_2, \text{C}_2\text{H}_4, \text{N}_2$), vibrations that are Infrared active (change in dipole moment, $\partial\mu/\partial q \neq 0$) are Raman inactive, and vibrations that are Raman active (change in polarizability, $\partial\alpha/\partial q \neq 0$) are Infrared inactive.
The pure rotational absorption spectrum of ${}^{12}\text{C}^{16}\text{O}$ exhibits a series of equidistant absorption lines with a constant wavenumber separation $\Delta\bar{\nu} = 3.842\text{ cm}^{-1}$. (a) Calculate the rotational constant $B$ in $\text{cm}^{-1}$ and in Joules. (b) Determine the moment of inertia $I_0$ of the CO molecule. (c) Given atomic masses $m(^{12}\text{C}) = 1.9926 \times 10^{-26}\text{ kg}$ and $m(^{16}\text{O}) = 2.6567 \times 10^{-26}\text{ kg}$, calculate the reduced mass $\mu$ and the equilibrium internuclear bond distance $r_0$ in angstroms.
Step 1: Calculate Rotational Constant $B$
For a rigid rotor, adjacent spectral line separation is $\Delta \bar{ u} = 2B$:
In energy units:
Step 2: Calculate Moment of Inertia $I$
Step 3: Calculate Reduced Mass $\mu$ and Bond Length $r_0$
Using $I = \mu r_0^2$:
The $\text{C-O}$ triple bond length is determined to five significant figures: $r_0 = 1.1312\,\text{\AA}$.
Step 4: Transition $J = 4 \to 5$
The fundamental infrared absorption band of hydrogen chloride (${}^1\text{H}^{35}\text{Cl}$) has its band origin (missing Q-branch) at wavenumber $\bar{\nu}_0 = 2885.9\text{ cm}^{-1}$. (a) Calculate the fundamental vibrational frequency $\nu_0$ in Hz and the effective chemical bond force constant $k$ in N/m (reduced mass $\mu = 1.6267 \times 10^{-27}\text{ kg}$). (b) Compute the zero-point vibrational energy (ZPE) of the molecule in eV. (c) Predict the fundamental vibrational wavenumber $\bar{\nu}_0'$ for deuterium chloride (${}^2\text{H}^{35}\text{Cl}$) assuming the same force constant.
Step 1: Calculate Reduced Mass $\mu$
Step 2: Calculate Bond Force Constant $k$
From $\bar{ u}_0 = \frac{1}{2\pi c} \sqrt{\frac{k}{\mu}} \implies \omega_0 = 2\pi c \bar{ u}_0$:
The $\text{H-Cl}$ single bond has a force constant of $k \approx 481\,\text{N/m}$.
Step 3: Calculate Wavenumbers of $P$ and $R$ Branches
Formulas: $\bar{ u}_R(J) = \bar{ u}_0 + 2B(J+1)$ and $\bar{ u}_P(J) = \bar{ u}_0 - 2BJ$, with $2B = 2(10.59) = 21.18\,\text{cm}^{-1}$:
- $R(0)$ ($J = 0 \to 1$): $\bar{ u} = 2886.0 + 21.18 = 2907.18\,\text{cm}^{-1}$
- $R(1)$ ($J = 1 \to 2$): $\bar{ u} = 2886.0 + 2(21.18) = 2928.36\,\text{cm}^{-1}$
- $P(1)$ ($J = 1 \to 0$): $\bar{ u} = 2886.0 - 21.18 = 2864.82\,\text{cm}^{-1}$
- $P(2)$ ($J = 2 \to 1$): $\bar{ u} = 2886.0 - 2(21.18) = 2843.64\,\text{cm}^{-1}$
Step 4: Central Gap Separation
The missing central $Q$-branch leaves a double-spacing gap of $4B = 42.36\,\text{cm}^{-1}$.
A gas of nitrogen molecules ($\text{N}_2$) is irradiated with a green Nd:YAG laser beam of wavelength $\lambda_0 = 532.0\text{ nm}$. The vibrational wavenumber of $\text{N}_2$ is $\Delta\bar{\nu} = 2331\text{ cm}^{-1}$. (a) Calculate the wavelengths of the Stokes and anti-Stokes Raman scattered lines in nanometers. (b) At room temperature $T = 300\text{ K}$, calculate the theoretical intensity ratio $I_{\text{anti-Stokes}} / I_{\text{Stokes}}$ using the Boltzmann factor and the $\nu^4$ scattering law. (c) Compute this ratio at $T = 1500\text{ K}$ and explain how Raman spectroscopy is used as an optical thermometer in combustion diagnostics.
Step 1: Calculate Incident Laser Wavenumber $\bar{\nu}_0$
Step 2: Calculate Stokes and Anti-Stokes Lines
Stokes Line (Red-Shifted):
Anti-Stokes Line (Blue-Shifted):
Step 3: Intensity Ratio at $T = 300\,\text{K}$
The intensity ratio is governed by dipole radiation $ u^4$ law and the Boltzmann factor:
At room temperature, the Anti-Stokes line is more than 26,000 times weaker than the Stokes line because virtually all nitrogen molecules reside in the ground vibrational state ($v = 0$).
Step 4: Temperature for 10% Intensity Ratio ($I_{AS}/I_S = 0.10$)
Solved University Examination Problems
Step-by-step mathematical solutions to classic university honors examination questions.