Spotlight #39 – Physical Chemistry & Reaction Dynamics: Transition States, Marcus Theory, Spectroscopy and Femtochemistry

Physical chemistry occupies the overlap between physics and chemistry: it applies the rigorous quantitative framework of physics to understand the structure, dynamics, and transformations of chemical systems. The field spans from equilibrium thermodynamics and electrochemistry through molecular quantum mechanics to ultrafast laser spectroscopy. This spotlight covers six central topics, each with the mathematical underpinning that makes physical chemistry genuinely predictive rather than merely descriptive.

1. Transition State Theory

Classical collision theory gives A-factor estimates but lacks the geometry of the potential-energy surface (PES). Henry Eyring’s transition state theory (TST, 1935) treats the saddle-point configuration (the “transition state” or activated complex) as being in quasi-thermodynamic equilibrium with reactants, then calculates the flux over the barrier using statistical mechanics.

Arrhenius rate law:
  k(T) = A ⋅ e^(−E_a/RT)
  A = pre-exponential (frequency) factor [s−¹ or M−¹s−¹]
  E_a = activation energy [kJ/mol]

Eyring (transition state) equation:
  k(T) = (k_B T / h) ⋅ (Κ)† ⋅ exp(−ΔG‡/RT)
  k_B T/h ≈ 6.25 × 10^12 s−¹ at 298 K
  ΔG‡ = ΔH‡ − TΔS‡   (free energy of activation)
  κ = transmission coefficient (≤1; accounts for recrossing)

Eyring plot:
  ln(k/T) = −ΔH‡/R ⋅ (1/T) + ln(k_B/h) + ΔS‡/R
  Slope = −ΔH‡/R,  intercept gives ΔS‡

Hammond postulate:
  Exothermic: TS is reactant-like → ΔE‡ ≈ small
  Endothermic: TS is product-like → ΔE‡ ≈ ΔH_rxn
  Evans-Polanyi: ΔE‡ = αΔH_rxn + β (linear free-energy relation)

Quantum tunnelling corrections:
  Wigner correction: κ_tun ≈ 1 + (1/24)(hν‡/k_BT)²
  Bell correction (wider parabola): more elaborate; significant for proton/H transfer
  DFT/path-integral methods: used for enzyme reactions and H-tunnelling
          

2. Marcus Electron Transfer Theory

Rudolph Marcus (Nobel Prize 1992) developed a quantitative theory of electron transfer (ET) reactions in solution that predicted a remarkable “inverted region”: above a threshold driving force, increasing the thermodynamic driving force decreases the rate. This counterintuitive prediction, confirmed in the 1980s, is central to photosynthesis, solar cells, and charge transport in organic electronics.

Outer-sphere electron transfer:
  D + A → D•+ + A•−   (no bond breaking, solvent reorganisation)

Marcus rate constant:
  k_ET = (2π/ℏ) |H_DA|² (1/√(4πλk_BT)) exp(−(ΔG° + λ)²/(4λk_BT))

  H_DA = electronic coupling matrix element (donor-acceptor overlap)
  λ   = total reorganisation energy = λ_inner + λ_outer
  ΔG° = standard free energy change of ET reaction

Reorganisation energies:
  λ_inner (intramolecular, bond length change):
    λ_i = (1/2)(f_D f_A / (f_D+f_A)) ⋅ (Δq_eq)²
  λ_outer (solvent polarisation):
    λ_o = (e²/2)(1/r_D + 1/r_A − 1/r_DA)(1/ε_op − 1/ε_s)

Three regimes (|H_DA| small → non-adiabatic):
  Normal:    |ΔG°| < λ  → barrier = (ΔG°+λ)²/4λ > 0
  Activationless: ΔG° = −λ → maximum rate (zero barrier)
  Inverted:  |ΔG°| > λ  → rate DECREASES as |ΔG°| increases further
    Confirmed by Miller, Calcaterra, Closs (1984) for intramolecular ET

Photosynthetic ET in Photosystem II reaction centre:
  P680* → pheophytin → QA → QB  (sequential ET within ~200 ps)
  Inverted region prevents back-ET at each step: quantum efficiency ~99%

Semiconductor solar cells:
  Grätzel cell (DSSC): dye excited state injects e− into TiO&sub2; CB
  Marcus λ engineered for fast forward ET, slow back ET via ΔG° in inverted region
          

3. Chemical Thermodynamics

The equilibrium state of a chemical system is determined by minimising the Gibbs free energy G at constant temperature and pressure. The chemical potential μ_i of each species is the gateway between thermodynamic tables and equilibrium constants, Le Chatelier’s principle, and phase diagrams.

Gibbs free energy:
  G = H − TS = U + pV − TS
  dG = −S dT + V dp + Σ_i μ_i dn_i
  Equilibrium at constant T,p: dG = 0

Chemical potential:
  μ_i = (∂G/∂n_i)_{T,p,n_j}
  Ideal mixture: μ_i = μ_i°(T) + RT ln x_i   (mole fraction x_i)
  Real solution: μ_i = μ_i°(T) + RT ln a_i   (activity a_i = γ_i x_i)
  Fugacity (real gases): μ = μ° + RT ln(f/p°)

Equilibrium constant:
  K(T) = exp(−Δ_r G°/RT) = ∏_i a_i^(ν_i)
  Van't Hoff: d ln K / dT = Δ_r H°/(RT²)

Clapeyron equation (phase boundary slope):
  dP/dT = ΔH/(TΔV)
  Solid-liquid: large dP/dT (small ΔV); 1-3 MPa/K for most
  Clausius-Clapeyron (vapour-liquid, ideal gas approximation):
    d ln P/dT = Δ_vap H / RT²
    P(T) = P_0 exp(−Δ_vap H/R ⋅ (1/T − 1/T_0))

Gibbs phase rule:
  F = C − P + 2    (F=degrees of freedom, C=components, P=phases)
  Binary two-phase (e.g. liquid+vapour): F=1 (isothermal ↔ isobaric equivalence)
  Triple point: F=0 (unique T,P for pure substance)
          

4. Spectroscopy

Spectroscopy is the primary experimental tool of physical chemistry. Beer-Lambert quantifies UV-visible absorption; rotational and vibrational spectroscopy (IR and Raman) probe molecular structure; NMR resolves chemical environments and coupling networks. Each technique connects to quantum mechanics via selection rules and energy level expressions.

Beer-Lambert law:
  A = εlc = −log(I/I_0)
  ε = molar absorption coefficient [M−¹cm−¹],  l = path length [cm]
  Linear in concentration for ideal dilute solutions

Rigid-rotor rotational energy:
  E_J = hcBJ(J+1),  B = h/(8π²Ic)  [cm−¹]
  I = μr_e²,  μ = reduced mass
  Selection rule: ΔJ = ±1 (microwave active → requires permanent dipole)
  Rotational spacing: Δν = 2B(J+1)

Harmonic oscillator vibrational energy:
  E_v = hν_e(v + 1/2),  v = 0,1,2,...
  ν_e = (1/2πc)√(k/μ)  [cm−¹],  k = force constant
  IR selection rule: Δv = ±1, requires dμ/dq ≠ 0 (dipole change)
  Raman selection rule: requires dα/dq ≠ 0 (polarisability change)
  Mutual exclusion rule: IR active ↔ Raman inactive for centrosymmetric molecules

Anharmonicity:
  E_v = hν_e(v+1/2) − hν_eχ_e(v+1/2)²
  D_e = hν_e/(4χ_e)  (dissociation energy from Morse potential)
  Overtones: Δv = ±2,±3 allowed (weak)

Pulsed FT-NMR:
  Larmor frequency: ω_L = γ B_0  (proton γ/2π = 267.5 MHz/T)
  Chemical shift δ (ppm): reports electron shielding environment
  Spin-spin coupling J (Hz): through-bond orbital overlap
  FID → Fourier transform → chemical shift spectrum
  2D COSY (coupled spins), NOESY (spatial proximity <5Å): structure determination
          

5. Quantum Chemistry

Quantum chemistry solves the electronic Schrödinger equation to obtain molecular energies, structures, and properties. The Born-Oppenheimer approximation separates electronic and nuclear motion; molecular orbital (MO) theory builds wavefunctions from atomic basis functions; and density functional theory (DFT) reformulates the problem in terms of electron density rather than many-electron wavefunctions.

Born-Oppenheimer approximation:
  Ψ_total ≈ Ψ_electronic(r;R) ċ Ψ_nuclear(R)
  Electronic Schrödinger: H_el Ψ_el = E_el(R) Ψ_el   (parametric in nuclear coords R)
  PES = E_el(R) = potential energy surface on which nuclei move

LCAO-MO theory (molecule = linear combination of atomic orbitals):
  ψ = Σ_μ c_μ φ_μ
  Secular equation: det(H − ES) = 0
  H_μν = ⟨φ_μ|H|φ_ν⟩,  S_μν = ⟨φ_μ|φ_ν⟩

Hückel method (π-electrons only, S_μν = δ_μν):
  H_μμ = α,  H_μν{ adjacent } = β (<0)
  Benzene C_6: E_k = α + 2βcos(2πk/6), k=0,1,2,3,4,5
  Delocalization energy = 2|β| per π bond vs ethylene reference

Hartree-Fock (HF):
  Approximate Ψ as single Slater determinant
  Fock operator: F(φ_i) = hφ_i + Σ_j(2J_j − K_j)φ_i
  Correlation energy: E_corr = E_exact − E_HF  (missing in HF; ~1 eV/bond)
  MP2, CCSD, CCSD(T): systematic improvements

Density Functional Theory (DFT):
  Hohenberg-Kohn theorems (1964): E uniquely determined by ρ(r)
  Kohn-Sham equations: auxiliary non-interacting system
    [−ℏ²∇²/(2m) + V_ext + V_H + V_xc] ψ_i = ε_i ψ_i
  ρ(r) = Σ_i|ψ_i(r)|²
  Exchange-correlation functional E_xc[ρ]: approximated (LDA, GGA, hybrid B3LYP)
  Computational cost: O(N³) vs O(N^5) for MP2; workhorse of computational chemistry
          

6. Femtochemistry

In 1999 Ahmed Zewail received the Nobel Prize in Chemistry for developing femtosecond spectroscopy — the use of 10⁻¹⁵-second laser pulses to “film” chemical reactions in real time. Femtochemistry revealed that transition states are not just mathematical constructs but observable transient structures with lifetimes of 100–500 fs.

Pump-probe experiment:
  Pump pulse: excites molecule to excited electronic state at t=0
  Probe pulse (delayed by Δt): detects product or intermediate via LIF or absorption
  Time resolution limited by pulse duration Δt ~ 50–200 fs

Franck-Condon principle:
  Electronic transition: nuclei do not move during electronic excitation (fast)
  FC factor: S_vv' = |⟨χ_v'|χ_v⟩|²  (overlap of nuclear wavefunctions)
  Large FC factor for transitions preserving nuclear geometry

Wavepacket dynamics on excited PES:
  Ψ(R,t) = Σ_v c_v χ_v(R) exp(−iE_vt/ℏ)  (coherent superposition)
  Wavepacket oscillates with period T_vib ~ 200–1000 fs
  Coherent oscillations visible as probe signal modulations vs Δt

Conical intersections:
  Degeneracy point between two electronic PES (S1 and S0 for ethylene)
  Ultra-fast non-radiative relaxation (10−100 fs) via conical intersection
  Key to photostability of DNA bases (thymine decay: ~1 ps via S2→CI→S0)

XFEL (X-ray Free Electron Laser) – 2010s-present:
  LCLS (SLAC), European XFEL, SACLA: sub-100 fs X-ray pulses
  Serial crystallography: molecular movies of protein dynamics
  Solvated electron dynamics (LiF crystal), nanocrystal diffraction
  Δt_XRD ≈ 25 fs → structural dynamics of bond lengths during chemistry
          

Try These Simulations