Compendium

Physical Chemistry — Thermodynamics, Kinetics, Quantum, Statistical Mech

33 min read

Physical Chemistry — Thermodynamics, Kinetics, Quantum, Statistical Mechanics

Physical chemistry sits at the intersection of physics and chemistry — the quantitative theory of why molecules behave the way they do. It explains why ice melts at 0 °C (32 °F), why catalysts speed reactions without being consumed, why a hydrogen atom’s emission spectrum has only certain frequencies, and why a polymer melt flows like honey while its solid form fractures. Four pillars support the field: classical thermodynamics (energy, entropy, free energy), chemical kinetics (rates and mechanisms), statistical mechanics (bridging microscopic states to macroscopic observables), and quantum chemistry (the electronic-structure theory underlying all of chemistry).

This note covers each pillar with the load-bearing equations, then surveys the spectroscopic toolkit — IR, Raman, UV-Vis, NMR, EPR, mass spectrometry — that lets us measure what theory predicts. Sections on polymer physics and surface chemistry close the loop into materials, where physical chemistry meets chemical-engineering and materials-overview.

1. Classical thermodynamics — the four laws

Thermodynamics describes how energy redistributes in systems of many particles, without needing to know what those particles are. Its four laws were assembled empirically over the 19th century and remain among the most rigorously tested statements in physical science.

Zeroth law — thermal equilibrium and temperature

If system A is in thermal equilibrium with system C, and B is also in equilibrium with C, then A and B are in equilibrium with each other. This transitivity is what defines temperature as a state variable: a thermometer (C) reads the same value when placed in contact with two systems already in equilibrium. The Kelvin scale (K) is the SI thermodynamic temperature, with 0 K = -273.15 °C (-459.67 °F) being absolute zero. The 2019 SI redefinition fixes the Boltzmann constant exactly at k_B = 1.380649 × 10^-23 J/K, which in turn defines the kelvin.

First law — energy conservation

ΔU = q + w, where U is internal energy (J), q is heat added to the system, and w is work done on the system. Both q and w are path-dependent (process variables), while U is a state function (depends only on initial and final states). For an ideal gas undergoing a reversible expansion, w_rev = -∫P dV, and for an isothermal ideal-gas process w = -nRT ln(V_f/V_i).

Heat capacities: C_V = (∂U/∂T)_V (constant volume), C_P = (∂H/∂T)_P (constant pressure). For an ideal gas, C_P - C_V = nR (Mayer’s relation). Monatomic ideal gas: C_V = (3/2)R = 12.47 J/(mol·K), C_P = (5/2)R = 20.79 J/(mol·K), γ = C_P/C_V = 5/3 ≈ 1.667. Diatomic gas at room temperature: C_V ≈ (5/2)R (two rotational degrees of freedom active), γ ≈ 1.4.

Second law — entropy and the arrow of time

dS_universe = dS_system + dS_surroundings ≥ 0 for any real process; equality holds only at reversible equilibrium. The Clausius statement: heat does not flow spontaneously from cold to hot. The Kelvin-Planck statement: no cyclic process converts heat entirely to work. Both are equivalent to dS_univ ≥ 0.

Carnot efficiency for a heat engine operating between T_h and T_c: η_Carnot = 1 - T_c/T_h. A coal-fired plant with T_h ≈ 833 K (560 °C / 1040 °F) and T_c ≈ 300 K (27 °C / 80 °F) has a Carnot limit of about 64%; real plants achieve 33-45% due to irreversibilities. Combined-cycle gas turbines with T_h up to 1700 K can hit 64% real efficiency (GE 9HA.02, Siemens SGT-9000HL). See thermodynamics-heat-transfer.

Entropy as disorder is a useful but imprecise picture; entropy is more rigorously the logarithm of the number of accessible microstates (Boltzmann’s S = k_B ln W on his tombstone in Vienna). For an ideal gas, Sackur-Tetrode gives S = nR[ln(V/N · (mk_BT/2πℏ²)^(3/2)) + 5/2], the most-celebrated closed-form entropy expression.

Third law — absolute entropy

S → 0 as T → 0 for a perfect crystal (Walther Nernst 1906; Nobel 1920). Empirically, this means a perfect crystal at absolute zero has only one microstate (W = 1, ln 1 = 0). Real materials have residual entropy from defects, isotope mixing, and frozen-in disorder — ice has S_residual ≈ 3.4 J/(mol·K) from the geometric frustration of hydrogen bonds (Pauling 1935).

2. State functions and thermodynamic potentials

Four state functions cover most situations:

  • U internal energy — natural variables S, V; dU = TdS - PdV (closed system, no work other than PV)
  • H = U + PV enthalpy — natural variables S, P; dH = TdS + VdP; ΔH = q_P at constant pressure
  • A = U - TS Helmholtz free energy — natural variables T, V; dA = -SdT - PdV; minimum at equilibrium for constant T, V
  • G = H - TS Gibbs free energy — natural variables T, P; dG = -SdT + VdP; minimum at equilibrium for constant T, P (most relevant to chemistry under atmospheric conditions)

For open systems with mass exchange:

dG = -SdT + VdP + Σ μ_i dn_i

where μ_i = (∂G/∂n_i){T,P,n{j≠i}} is the chemical potential of species i. Chemical potential is the central concept in solution thermodynamics, phase equilibria, and chemical reactions — a species flows from regions of high μ to low μ, just as heat flows from high T to low T.

Maxwell relations arise from equality of mixed partial derivatives. From dG = -SdT + VdP, (∂S/∂P)_T = -(∂V/∂T)_P. Useful for converting hard-to-measure quantities (entropy) into easier ones (thermal expansion coefficient).

Reaction thermodynamics

ΔG°_rxn = Σ ν_i ΔG°_f,i (products) - Σ ν_i ΔG°_f,i (reactants), using standard Gibbs energies of formation at 1 bar, typically 298.15 K. The standard state for solutes is 1 mol/kg activity; for gases, 1 bar fugacity.

ΔG = ΔG° + RT ln Q (Q is the reaction quotient)

At equilibrium ΔG = 0, so ΔG° = -RT ln K. This connects free-energy tables (thermodynamic) to equilibrium constants (measurable). For a reaction with ΔG° = -10 kJ/mol at 298 K, K = exp(10000/(8.314·298)) ≈ 57.

3. Phase equilibria

Clausius-Clapeyron equation

For coexistence between two phases (liquid-vapor, solid-vapor, solid-liquid):

dP/dT = ΔH_trans / (T ΔV_trans)

For liquid-vapor with ideal gas vapor and negligible liquid volume:

ln(P_2/P_1) = -(ΔH_vap/R)(1/T_2 - 1/T_1)

Water has ΔH_vap = 40.7 kJ/mol at 100 °C (212 °F); the equation correctly predicts vapor pressure curves over modest temperature ranges. Calibrated forms include the Antoine equation: log_10 P = A - B/(T + C), with empirical A, B, C for each substance (NIST WebBook).

Gibbs phase rule

F = C - P + 2, where F is the number of degrees of freedom (intensive variables freely chosen), C is the number of components, P is the number of phases coexisting. For pure water (C=1) at a triple point (P=3), F = 0 — temperature and pressure are uniquely determined (273.16 K, 611.657 Pa for the H2O triple point, used as a former kelvin reference).

Notable phase diagrams

  • Water: triple point 273.16 K / 611.657 Pa; critical point 647.096 K (374 °C / 705 °F), 22.064 MPa (3200 psi); ice-Ih is less dense than liquid water (anomaly); supercritical water above critical point is a strong solvent for non-polar species, exploited in supercritical water oxidation (SCWO) for waste destruction.
  • CO₂: triple point 216.59 K / 5.18 bar; critical point 304.13 K (31 °C / 88 °F) and 7.39 MPa (1071 psi); supercritical CO₂ is widely used in decaffeination (Café HAG, Maxwell House), dry cleaning, and as a process fluid in heat pumps and Allam-cycle power plants (NET Power’s 8 Rivers technology).
  • Sulfur: two solid allotropes (rhombic and monoclinic), liquid and gas; classic illustration of polymorphism.

See phase-diagrams and cryogenics-low-temperature.

4. Solutions and mixtures

Ideal mixing

For an ideal mixture: ΔH_mix = 0, ΔV_mix = 0, but ΔS_mix = -R Σ x_i ln x_i (always positive), so ΔG_mix = RT Σ x_i ln x_i (always negative). Entropy of mixing is what drives miscibility in the absence of enthalpic barriers.

Raoult’s law and Henry’s law

For an ideal liquid mixture in equilibrium with its vapor: P_i = x_i · P_i*, where P_i* is the pure-component vapor pressure (Raoult). For dilute solutes: P_i = K_H · x_i (Henry’s law constant K_H), describing gas dissolution in liquids (CO₂ in soda, O₂ in blood). At 298 K, K_H for O₂ in water is ≈ 4.4 × 10⁴ bar (mole-fraction basis).

Real mixtures

Activity a_i = γ_i · x_i, where γ_i is the activity coefficient (γ = 1 for ideal). Models:

  • Margules — empirical polynomial fit to ln γ vs composition; two-parameter form widely used
  • van Laar — derived assuming regular solution
  • Wilson (Grant Wilson 1964) — handles strongly polar mixtures, accurate but cannot represent liquid-liquid splits
  • NRTL (Non-Random Two-Liquid; Renon & Prausnitz 1968) — handles liquid-liquid equilibria
  • UNIQUAC and UNIFAC — group-contribution methods (Abrams & Prausnitz 1975; Fredenslund 1975); UNIFAC predicts γ from molecular structure with no experimental data, used in Aspen Plus, COSMOtherm, ChemCAD

Activity-coefficient models drive process simulators for distillation, extraction, and absorption design — covered in separations-distillation-extraction.

Colligative properties

Depending on solute concentration but not solute identity:

  • Vapor pressure depression: ΔP = -x_solute · P_solvent*
  • Boiling point elevation: ΔT_b = K_b · m (water K_b = 0.512 K·kg/mol)
  • Freezing point depression: ΔT_f = -K_f · m (water K_f = 1.86 K·kg/mol); salt on roads, antifreeze in radiators
  • Osmotic pressure: π = M · R · T (van’t Hoff; van’t Hoff 1901 first Nobel chemistry prize); biological membranes, reverse osmosis desalination (water-treatment)

5. Acid-base, buffers, solubility

pH = -log_10[H⁺·activity] (strictly), -log_10[H⁺] in dilute aqueous solution. K_w = 10⁻¹⁴ at 25 °C, so [H⁺][OH⁻] = K_w, and pH + pOH = 14 in pure water.

K_a = [H⁺][A⁻]/[HA]; pK_a = -log K_a. Strong acids (HCl, HNO₃, H₂SO₄ first ionization, HClO₄) have pK_a < 0 — essentially fully dissociated. Weak acids span pK_a 2-12: acetic acid 4.76, hydrofluoric 3.17, ammonium 9.25, phenol 9.95, water 15.7.

Henderson-Hasselbalch for buffer pH:

pH = pK_a + log_10([A⁻]/[HA])

Buffer capacity is highest within ±1 unit of pK_a. Common biological buffers: phosphate (pK_a 7.2, used in PBS), HEPES (pK_a 7.5, no metal binding, life sciences standard), TRIS (pK_a 8.06, but T-dependent), bicarbonate (blood, pK_a 6.35).

Titration curves for monoprotic, polyprotic, and complex equilibria — sharp inflection at equivalence point, half-equivalence point gives pK_a directly. Endpoint indicators (phenolphthalein pH 8.3-10.0, methyl orange 3.1-4.4) or pH-meter detection.

Solubility product: for AgCl(s) ⇌ Ag⁺ + Cl⁻, K_sp = [Ag⁺][Cl⁻] = 1.8 × 10⁻¹⁰ at 25 °C. Common-ion effect, pH-dependent solubility of hydroxides and sulfides, complex equilibria (Ag(NH₃)₂⁺ formation increases Ag salt solubility) — quantitatively predicted by simultaneous mass-balance and charge-balance equations.

6. Chemical kinetics — rate laws and mechanisms

Rate laws

For a reaction aA + bB → products: rate = -(1/a)d[A]/dt = -(1/b)d[B]/dt = k[A]^m[B]^n. The orders m and n are empirical and need not match stoichiometric coefficients. Overall order = m + n.

Integrated rate laws:

  • Zero order: [A] = [A]₀ - kt; t_{1/2} = [A]₀/(2k); k has units of M/s
  • First order: ln[A] = ln[A]₀ - kt; t_{1/2} = ln 2 / k = 0.693/k; k has units of 1/s; independent of [A]₀
  • Second order (in A): 1/[A] = 1/[A]₀ + kt; t_{1/2} = 1/(k[A]₀); k has units of 1/(M·s)
  • n-th order (n ≠ 1): 1/(n-1) · ([A]^(1-n) - [A]₀^(1-n)) = kt

Arrhenius and Eyring

k = A · exp(-E_a / RT) — Svante Arrhenius 1889; Nobel 1903 (the third Nobel in chemistry ever). Plotting ln k vs 1/T yields slope -E_a/R, intercept ln A. The pre-exponential factor A reflects collision frequency and orientation.

Eyring (Henry Eyring, Princeton, 1935) transition-state theory:

k = (k_B T / h) · exp(-ΔG‡/RT) = (k_B T / h) · exp(ΔS‡/R) · exp(-ΔH‡/RT)

Activation parameters ΔS‡ and ΔH‡ are extracted from a plot of ln(k/T) vs 1/T. Negative ΔS‡ implies an associative transition state (bimolecular); positive ΔS‡ implies dissociative.

Collision theory (Trautz, Lewis): k = Z · ρ · exp(-E_a/RT), where Z is the collision frequency and ρ is the steric factor (often 10⁻⁵ to 1).

Reaction mechanisms

A multi-step mechanism is a series of elementary reactions. The rate-determining step (RDS) is the slowest; the overall rate law is derived from it, possibly with pre-equilibrium or steady-state assumptions.

  • Pre-equilibrium: fast equilibrium followed by slow step; pre-equilibrium gives concentrations of intermediates
  • Steady-state approximation (Bodenstein): d[I]/dt ≈ 0 for short-lived intermediates; rearrange to eliminate [I]
  • Lindemann-Hinshelwood unimolecular: A + M ⇌ A*+ M (M = inert third body); A* → products; explains pressure dependence of unimolecular reactions (RRKM theory extends this with statistical microcanonical rate constants)

Enzyme kinetics — Michaelis-Menten

v = V_max [S] / (K_M + [S])

where V_max = k_cat · [E]total, K_M = (k-1 + k_cat)/k_1 ≈ K_d at high k_cat. K_M is the substrate concentration giving v = V_max/2.

Lineweaver-Burk double reciprocal: 1/v = (K_M/V_max)·(1/[S]) + 1/V_max — historically used to extract K_M and V_max from a linear plot (intercept 1/V_max, slope K_M/V_max). Modern practice: nonlinear regression of the Michaelis-Menten equation directly (more accurate, doesn’t over-weight low-[S] data). Diagnostic specificity constant k_cat/K_M — for “perfect” enzymes (catalase, triose phosphate isomerase, fumarase, acetylcholinesterase) approaches the diffusion limit ~10⁸-10⁹ M⁻¹s⁻¹.

Heterogeneous catalysis

  • Langmuir-Hinshelwood: both reactants adsorb on surface, react, products desorb — typical for CO oxidation on Pt, Fischer-Tropsch on Co/Fe
  • Eley-Rideal: one reactant adsorbs, the other strikes from gas phase — H + D-on-Cu (less common in industry)
  • Mars-van Krevelen: lattice oxygen of metal oxide is the active reagent; oxide is reduced, then reoxidized — selective oxidation of olefins on V₂O₅, MoO₃; SCR DeNOx on V₂O₅/WO₃/TiO₂

Sabatier principle (Paul Sabatier; Nobel 1912): the best catalyst binds the intermediate just right — not too weakly (no activation) and not too strongly (poisoning). Plotting activity vs binding energy traces out a volcano plot — the maximum identifies the optimum. Modern computational catalysis (Nørskov, Stanford/DTU) uses DFT-computed adsorption energies as descriptors to screen catalysts; this approach found platinum-group alternatives for fuel cells and developed activity-selectivity relationships for the oxygen evolution reaction. See electrochemistry for ORR/OER details.

7. Statistical mechanics — bridging micro and macro

Statistical mechanics derives thermodynamics from molecular degrees of freedom. The central object is the partition function:

Z = Σ_i exp(-ε_i / k_B T) (over states i)

For independent particles: Z_total = Z_single^N / N! (Boltzmann statistics, dilute classical limit).

Macroscopic observables from Z:

  • ⟨E⟩ = k_B T² (∂ ln Z / ∂T)_V
  • A = -k_B T ln Z (Helmholtz free energy)
  • S = k_B ln Z + ⟨E⟩/T = -k_B Σ p_i ln p_i (Gibbs entropy)
  • P = k_B T (∂ ln Z / ∂V)_T

The Boltzmann distribution p_i = exp(-ε_i/k_B T) / Z is the most-probable distribution of particles over energy levels at thermal equilibrium — the foundation of every statistical-mechanics calculation.

Ensembles

  • Microcanonical (NVE): closed isolated system; fixed energy
  • Canonical (NVT): closed thermostatted system; fixed temperature; partition function Z(N,V,T)
  • Grand canonical (μVT): open system; particle exchange; partition function Ξ(μ,V,T) = Σ exp(βμN) Z(N,V,T)

Most chemistry calculations live in the canonical ensemble. Molecular dynamics simulations use Nosé-Hoover or Langevin thermostats to sample NVT.

Equipartition

Each quadratic degree of freedom contributes (1/2)k_B T per molecule to the energy. A diatomic gas has 3 translational + 2 rotational + 2 vibrational (potential + kinetic), giving 7·(1/2)k_B T = (7/2)k_B T per molecule — but vibrational modes are typically frozen out at room temperature because ℏω >> k_B T, so the effective value is (5/2)k_B T, matching the C_V = (5/2)R observation.

Partition functions for ideal gas

Translational: Z_trans = V/Λ³, where Λ = h/√(2π m k_B T) is the thermal de Broglie wavelength. For Ar at 298 K, Λ ≈ 16 pm (much smaller than interatomic spacing — classical regime).

Rotational (linear molecule): Z_rot = T/(σ Θ_rot), with Θ_rot = ℏ²/(2 I k_B). For N₂, Θ_rot ≈ 2.88 K, σ = 2 (symmetry number).

Vibrational (single mode): Z_vib = 1/(1 - exp(-Θ_vib/T)), with Θ_vib = ℏω/k_B. For N₂, Θ_vib ≈ 3374 K — vibrational mode frozen out at 298 K.

Putting it together, the Sackur-Tetrode equation for an ideal monatomic gas:

S = N k_B [ln(V/N · (mk_B T/2πℏ²)^(3/2)) + 5/2]

Predicted entropies match experimental values to ~0.1 J/(mol·K) for noble gases — a stunning success of statistical mechanics.

Maxwell-Boltzmann speed distribution

f(v) = 4π (m/2πk_B T)^(3/2) v² exp(-mv²/2k_B T)

Most probable speed v_mp = √(2k_B T/m); mean speed ⟨v⟩ = √(8k_B T/πm); root-mean-square v_rms = √(3k_B T/m). For N₂ at 298 K, v_rms ≈ 515 m/s (1690 ft/s — about Mach 1.5).

8. Quantum chemistry — postulates and foundations

Chemistry is a quantum phenomenon. Bonds, spectra, magnetism, color, reactivity — all require quantum mechanics for a quantitative account.

Postulates

  1. The state of a system is described by a wavefunction Ψ(r,t) satisfying iℏ ∂Ψ/∂t = Ĥ Ψ.
  2. Observables correspond to Hermitian operators (eigenvalues are real); position x̂, momentum p̂ = -iℏ∂/∂x, energy Ĥ, angular momentum L̂.
  3. Measurement of an observable yields one of the operator’s eigenvalues, with probability |⟨ψ_n|Ψ⟩|².
  4. After measurement, the state collapses to the corresponding eigenstate.

Commutators [Â, B̂] = ÂB̂ - B̂Â. Position-momentum: [x̂, p̂] = iℏ → Heisenberg uncertainty Δx Δp ≥ ℏ/2.

Solved one-particle models

  • Particle in a box (1D, length L): E_n = n²h²/(8mL²); separation Δε = (2n+1)h²/(8mL²). Quantum dots (CdSe/InP, Quantum Dot Corporation acquired by Thermo Fisher; Nanosys; Nanoco) exploit this — bandgap tuned by particle size (Brus 1980s; Bawendi, Brus, Ekimov shared 2023 Nobel chemistry).
  • Harmonic oscillator: E_v = ℏω(v + 1/2), v = 0,1,2,…; ω = √(k/μ); used for vibrational modes; selection rule Δv = ±1 for IR.
  • Rigid rotor: E_J = J(J+1)ℏ²/(2I); rotational spectroscopy at microwave frequencies; selection rule ΔJ = ±1.
  • Hydrogen atom: E_n = -13.6 eV / n² (Rydberg); quantum numbers n, l, m_l, m_s; orbitals 1s, 2s/2p, 3s/3p/3d, 4s/4p/4d/4f. Schrödinger 1926 solved this analytically — the first triumph of quantum chemistry.

Born-Oppenheimer approximation

Electrons move much faster than nuclei (mass ratio ≥ 1836 for H; >50000 for heavy atoms). Separate the electronic Schrödinger equation (nuclei fixed) from nuclear motion (on the resulting potential energy surface). Underlies almost all of computational chemistry — fails near conical intersections (photochemistry, internal conversion).

Multi-electron atoms and molecules

Pauli exclusion + electron correlation make this hard. The independent-electron approximation:

  • Hartree-Fock (HF, SCF) — antisymmetric Slater determinant of one-electron orbitals; self-consistent iteration; mean-field. Captures ~99% of total energy but misses correlation energy (~1%, but this is what makes chemistry).
  • Post-HF methods systematically improve on HF:
    • Configuration Interaction (CI) — linear combination of excited determinants; CISD adds singles & doubles; truncated CI is non-size-extensive
    • Møller-Plesset perturbation (MP2, MP3, MP4) — MP2 is the workhorse for ground states of organic molecules
    • Coupled Cluster — CCSD, CCSD(T); the latter is the “gold standard” of quantum chemistry, often within 1 kcal/mol of experiment when paired with large basis sets
  • Multireference methods (CASSCF, CASPT2, NEVPT2) — needed for biradicals, transition metals, excited states with strong static correlation

Basis sets

The wavefunction is expanded in atom-centered Gaussian functions (Boys 1950). Common sets:

  • Minimal: STO-3G (3 primitives per orbital — qualitative only)
  • Pople split-valence: 6-31G, 6-31G(d) = 6-31G*, 6-31G(d,p) = 6-31G**, 6-311G**, 6-311+G(2d,2p) — workhorses
  • Dunning correlation-consistent: cc-pVDZ, cc-pVTZ, cc-pVQZ, cc-pV5Z; augmented (aug-cc-pVnZ) for anions, polarizabilities, excited states; extrapolate to CBS (complete basis-set) limit
  • Effective core potentials (ECP, pseudopotentials): LANL2DZ, def2-TZVP, SBKJC — replace core electrons with potential, treat valence explicitly; essential for transition metals and heavy elements

Density functional theory (DFT)

Hohenberg-Kohn theorems (1964): the ground-state electron density determines the external potential (and hence everything). Kohn-Sham (1965): replace the interacting many-electron problem with non-interacting electrons in an effective potential. Walter Kohn shared the 1998 Nobel chemistry prize with John Pople.

The exchange-correlation functional is the unknown; approximations climb Jacob’s ladder (Perdew):

  1. LDA (Local Density Approximation): SVWN, depends on ρ only; underestimates bond lengths
  2. GGA (Generalized Gradient Approximation): PBE (Perdew-Burke-Ernzerhof 1996 — most-cited paper in physics), BLYP, BP86; depend on ρ and ∇ρ
  3. meta-GGA: TPSS, SCAN (Sun-Ruzsinszky-Perdew 2015); depend additionally on kinetic energy density τ
  4. Hybrid GGA: B3LYP (Becke 3-parameter Lee-Yang-Parr — the most-used functional in chemistry; despite known flaws, dominates organic chemistry papers since 1993), PBE0, M06-2X (Truhlar); mix a fraction of exact (Hartree-Fock) exchange
  5. Hybrid meta-GGA: M06, ωB97X-V — range-separated for long-range Coulomb

Dispersion corrections are essential for non-covalent interactions (London dispersion, π-stacking, host-guest). Grimme D3, D3(BJ), D4 add a pair-additive C₆/R⁶ + C₈/R⁸ term post-SCF — used in DFT-D3, B3LYP-D3(BJ), ωB97X-D, etc.

Computational chemistry software

  • Gaussian (Gaussian Inc., John Pople; commercial; Gaussian 16 / Gaussian 23) — the most-cited quantum chemistry program; ubiquitous in academia and industry
  • Q-Chem (Marty Head-Gordon, UC Berkeley) — strong DFT, TDDFT, fragment methods; Q-Chem 6
  • ORCA (Frank Neese, MPI Mülheim; free for academia) — outstanding for EPR, multireference, ROHF/QROHF; ORCA 5/6
  • NWChem (PNNL; free, open) — massively parallel; plane-wave + Gaussian
  • MOLPRO (Werner, Knowles; commercial) — best-in-class CCSD(T) and MRCI
  • GAMESS-US (Mark Gordon, Iowa State; free) — legacy but maintained
  • Turbomole (commercial) — efficient RI-DFT and RI-MP2
  • Psi4 (Sherrill, Georgia Tech; open) — Python-scriptable, modern

Periodic / solid-state (plane-wave or LCAO):

  • VASP (Vienna; commercial; Kresse, Hafner) — dominant in materials, PAW pseudopotentials
  • Quantum ESPRESSO (open, Giannozzi et al.) — plane-wave DFT
  • CP2K (open) — Gaussian + plane-wave hybrid, fast for biomolecular and condensed-phase MD
  • CASTEP (commercial, BIOVIA) — plane-wave, used in industry
  • CRYSTAL (Torino) — Gaussian LCAO for periodic systems

Surveyed broadly in computational-materials-science.

9. Spectroscopy — measuring electronic, vibrational, rotational, magnetic states

Spectroscopy probes energy-level differences by photon absorption/emission. Fermi’s golden rule gives transition rate ∝ |⟨f|Ĥ_int|i⟩|² ρ(E_f).

Beer-Lambert law: A = ε l c, where A = log₁₀(I₀/I), ε is molar absorptivity (M⁻¹cm⁻¹), l is path length, c is concentration. Underlies UV-Vis quantitation, NIR moisture/protein in agriculture and food, IR analysis.

Line shapes: Lorentzian (lifetime-limited, natural broadening; FWHM = 1/πτ); Gaussian (inhomogeneous, Doppler in gases at low pressure); Voigt (convolution of both). Pressure broadening dominates IR at atmospheric pressure.

Vibrational spectroscopy

Infrared (IR) absorption: a mode is IR-active if its vibration changes the dipole moment (dμ/dq ≠ 0). Polar bonds (C=O at 1700-1750 cm⁻¹, O-H at 3200-3600 cm⁻¹, N-H at 3300-3500 cm⁻¹, C-H at 2850-3100 cm⁻¹) dominate. Mid-IR 4000-400 cm⁻¹ (2.5-25 μm) is the workhorse range.

FTIR (Fourier-transform IR) — Michelson interferometer with a moving mirror creates an interferogram; FFT gives the spectrum. Vastly faster, higher-SNR than dispersive. Standard since 1970s; Bruker Vertex 70/Tensor II, Thermo Scientific Nicolet iS50, Agilent Cary 630 / 660. Sampling: transmission, ATR (attenuated total reflection — diamond, ZnSe; Bruker Alpha II), DRIFTS, PAS.

Raman spectroscopy: mode is Raman-active if the polarizability changes (dα/dq ≠ 0). Stokes scattering (red-shifted) and anti-Stokes (blue-shifted, weaker). Complementary to IR — homonuclear bonds (C-C, S-S) are IR-silent but Raman-active. Mutual exclusion for centrosymmetric molecules: a mode cannot be both IR and Raman active.

Confocal Raman: Renishaw inVia Qontor, Horiba LabRAM HR Evolution, WITec alpha300. SERS (surface-enhanced Raman; rough Ag/Au) — Fleischmann 1974, ~10⁶-10¹⁰ enhancement, single-molecule sensitivity. TERS adds AFM tip for nanoscale. CARS (coherent anti-Stokes) and stimulated Raman scattering (SRS) — fast bio-imaging without labels.

Selection rules depend on molecular symmetry (point group analysis). For a normal mode to be IR active, it must transform as x, y, or z; Raman active, it must transform as a quadratic (xx, xy, etc.).

Group frequencies (functional-group fingerprints):

  • C=O 1650-1850 cm⁻¹ (acid > ester > aldehyde ≈ ketone > amide); conjugation lowers
  • O-H 3200-3650 cm⁻¹ broad (H-bonded)
  • N-H 3300-3500 cm⁻¹; two bands for primary amines
  • C≡C 2100-2260 cm⁻¹ weak; C≡N 2200-2260 cm⁻¹
  • C-H stretch 2850-3100 (alkane), 3000-3100 (alkene/aromatic), 3300 (alkyne)
  • Aromatic C=C 1580, 1600 cm⁻¹

Fermi resonance: accidental degeneracy between a fundamental and an overtone of similar symmetry causes mixing and splitting — classic example CO₂ ν₁ (1340 cm⁻¹) and 2ν₂ (1340 cm⁻¹) give Fermi doublet at 1285/1388 cm⁻¹.

Rotational spectroscopy

Microwave region (1 GHz - 1 THz, wavenumbers 0.03 - 30 cm⁻¹). For a linear rotor: ν = 2B(J+1), where B = ℏ/(4π c I) is the rotational constant (cm⁻¹). Centrifugal distortion: ν = 2B(J+1) - 4D(J+1)³. Used in atmospheric and astrochemistry — Sgr B2, IRC+10216 spectra revealed ~250 interstellar molecules. Lab instruments: chirped-pulse FT microwave spectroscopy (Brooks Pate, UVA) — Hz resolution.

Electronic spectroscopy

UV-Vis (200-800 nm; 50000-12500 cm⁻¹): electronic transitions. Common: π→πin aromatics (250-280 nm for benzene), n→π in ketones (270-290 nm, weak), d-d in transition metal complexes (visible, weak — color of Cu²⁺ pentaaquo at 800 nm), charge-transfer (intense, e.g. KMnO₄ purple at 525 nm). Instruments: Agilent Cary 60/3500, Shimadzu UV-2700, Thermo Evolution 220.

Franck-Condon principle: electronic transitions are vertical on a potential-energy surface (nuclei fixed during electron jump). Vibrational structure in absorption/emission spectra reflects overlap between vibrational wavefunctions of ground and excited electronic states.

Fluorescence (S₁ → S₀): typical lifetime ns. Stokes shift (red-shifted emission relative to absorption) — reorganization in the excited state. Jablonski diagram organizes radiative (fluorescence, phosphorescence) and non-radiative (internal conversion, intersystem crossing) processes.

Phosphorescence (T₁ → S₀, spin-forbidden): lifetime μs-s; quenched by O₂. Used in OLED triplet harvesting (Forrest, Thompson; UDC PHOLED, used in Samsung/Apple displays).

FRET (Förster resonance energy transfer): R⁻⁶ distance dependence, R₀ typically 1-10 nm; molecular ruler for protein conformation, single-molecule biophysics.

TDDFT (time-dependent DFT) is the workhorse for excited-state calculations on medium-sized molecules; CASPT2/NEVPT2 for the demanding cases.

NMR — nuclear magnetic resonance

Nuclei with non-zero spin (I ≠ 0) precess in a magnetic field at the Larmor frequency ω = γ B₀. ¹H (γ = 42.58 MHz/T, I = 1/2, abundance 99.99%), ¹³C (10.71 MHz/T, I = 1/2, abundance 1.1%), ¹⁵N (-4.32 MHz/T, I = 1/2, abundance 0.37%), ³¹P (17.25 MHz/T, I = 1/2, abundance 100%), ¹⁹F (40.05 MHz/T, I = 1/2, abundance 100%) are the most-used.

Chemical shift δ (ppm) = (ν_sample - ν_ref)/ν_spectrometer × 10⁶, referenced to TMS (tetramethylsilane) at 0 ppm for ¹H/¹³C. Typical ranges: alkyl H 0-2 ppm, vinyl H 5-7 ppm, aromatic H 6-9 ppm, aldehyde H 9-10 ppm; aliphatic C 0-50 ppm, sp²-C 100-150 ppm, carbonyl C 160-220 ppm.

J coupling (scalar): through-bond, Hz, splits peaks into multiplets; n+1 rule for equivalent neighbors. Vicinal ³J_HH 0-15 Hz, Karplus relation predicts torsion angle.

Relaxation: T₁ (longitudinal, spin-lattice; ~100 ms - seconds for ¹H, longer for ¹³C), T₂ (transverse, spin-spin; line width = 1/(πT₂)). Field gradients, viscosity, paramagnetic centers all affect T₁/T₂.

2D NMR revolutionized structure determination (Richard Ernst, Nobel 1991; Kurt Wüthrich, Nobel 2002 for protein NMR):

  • COSY (correlation spectroscopy): ¹H-¹H J-coupling correlations
  • NOESY (nuclear Overhauser): through-space (< 5 Å) ¹H-¹H proximities; gives 3D structure
  • HSQC (heteronuclear single-quantum coherence): ¹H-¹³C or ¹H-¹⁵N single-bond
  • HMBC (heteronuclear multiple-bond): 2-3 bond ¹H-¹³C correlations; backbone connectivity
  • TOCSY, HOESY, triple-resonance HNCA/HNCACB for protein backbone assignment

High-field magnets: 400-600 MHz benchtop standard (cryogen-free); 700-900 MHz routine; 1.0 GHz, 1.1 GHz, 1.2 GHz (28.2 T) at Bruker Avance Neo and JEOL ECZ — superconducting joints, Bi-2212 HTS inserts; for biomolecular NMR. NMR data analysis: MestReNova (Mestrelab), TopSpin (Bruker), Sparky, NMRPipe.

EPR — electron paramagnetic resonance

Same physics as NMR but for unpaired electrons. g-factor (free electron g_e = 2.0023; deviations probe spin-orbit coupling), hyperfine coupling (A, MHz or G) to nuclei (¹H, ¹⁴N, ⁵⁵Mn, ⁶³,⁶⁵Cu) reveals electronic structure of radicals, transition-metal complexes, defects in semiconductors. ENDOR (electron-nuclear double resonance) for resolved hyperfine; DEER/PELDOR for distance measurements (1-8 nm) in biomolecules. X-band (9.5 GHz), Q-band (35 GHz), W-band (94 GHz). Instruments: Bruker EMX, ELEXSYS-II E580 pulsed.

Mass spectrometry

Separates ions by mass-to-charge ratio (m/z). Three stages: ion source, mass analyzer, detector.

Ion sources:

  • EI (electron ionization, 70 eV) — fragmentation, gas-phase volatile small molecules; classical GC-MS
  • CI (chemical ionization) — softer than EI
  • ESI (electrospray; John Fenn, Nobel 2002) — biomolecules in solution; multiple charging gives large m/z range; standard for proteomics
  • MALDI (matrix-assisted laser desorption ionization; Koichi Tanaka, Nobel 2002) — peptides, polymers; singly-charged usually; MALDI-TOF imaging mass spec
  • APCI — atmospheric pressure chemical ionization, intermediate; common in LC-MS for less polar analytes
  • DESI, DART — ambient ionization, no sample prep

Mass analyzers:

  • Quadrupole (Q): RF + DC voltages select m/z; unit mass resolution; cheap, widely-used (Agilent 6470 triple-quad)
  • Time-of-flight (TOF): pulsed acceleration, drift tube; high resolution (50000-80000); Bruker timsTOF Pro, Sciex X500B, Waters Vion IMS
  • Ion trap (3D Paul trap, 2D linear trap): tandem MS in time
  • FT-ICR (Fourier-transform ion cyclotron resonance; Alan Marshall, NHMFL): highest resolution (>10⁶), highest mass accuracy (sub-ppm); 7-21 T magnets; Bruker SolariX
  • Orbitrap (Alexander Makarov, Thermo Fisher; introduced 2005): high resolution (240000-500000), benchtop; Thermo Q Exactive, Orbitrap Eclipse, Orbitrap Astral (2023) — dominant in proteomics and metabolomics

Software: Thermo Xcalibur, Bruker DataAnalysis, MZmine, Skyline (open), MaxQuant for proteomics. Top-down intact-protein MS (FT-ICR/Orbitrap, ETD fragmentation) vs bottom-up tryptic peptide MS.

10. Polymer physics

Flory-Huggins theory: ΔG_mix / k_B T = (φ₁/N₁) ln φ₁ + (φ₂/N₂) ln φ₂ + χ φ₁ φ₂; the χ parameter measures monomer-monomer interaction strength. χ_critical = (1/2)(1/√N₁ + 1/√N₂)² — polymers mix poorly compared to small molecules. Paul Flory, Nobel 1974.

Coil dimensions:

  • Mean-square end-to-end distance: ⟨R²⟩ = N a² (ideal chain), N^(6/5) a² (good solvent — swollen), N^(2/3) a² (poor solvent — collapsed globule)
  • Radius of gyration R_g; theta solvent (χ = 1/2) recovers ideal chain behavior; cyclohexane is a theta solvent for polystyrene at 34.5 °C
  • Persistence length ℓ_p: ds-DNA 50 nm (~150 bp), polystyrene 0.9 nm, F-actin 17 μm

Polymer dynamics:

  • Rouse model: unentangled chains, viscosity η ∝ N (linear)
  • Zimm model: includes hydrodynamic interactions in dilute solution
  • Reptation (de Gennes 1971; Nobel 1991): entangled chains in melt slither along a “tube” defined by entanglement constraints; viscosity η ∝ N³·⁴ experimentally; relaxation time τ_d ∝ N³. M_e (entanglement molecular weight) is material-specific: PE ~1.25 kg/mol, PS ~13 kg/mol, PDMS ~10 kg/mol.

Polymer physics flows into organic-chemistry-foundations for polymerization mechanisms and polymers for processing and properties.

11. Surface chemistry

Adsorption isotherms

  • Langmuir (1916; Nobel 1932): monolayer, equivalent sites, no lateral interaction: θ = K·P / (1 + K·P)
  • Freundlich: empirical, heterogeneous surfaces: q = K_F · P^(1/n) (1/n typically 0.1-1)
  • BET (Brunauer-Emmett-Teller 1938): multilayer adsorption, the standard for surface area measurement via N₂ adsorption at 77 K; isotherm linearized to extract monolayer capacity → surface area (m²/g). BET areas typical: activated carbon 500-2000, zeolites 300-900, MOFs (NU-1501-Al 7310 m²/g — record), silica gel 200-800.

Instruments: Micromeritics ASAP 2020/2460, Quantachrome (Anton Paar) Autosorb, BELSORP. IUPAC classification — Type I (microporous, zeolites), II (non-porous), III (weak), IV (mesoporous w/ hysteresis), V (weak mesoporous), VI (stepwise).

Wetting and capillarity

Young’s equation for contact angle θ: γ_sv = γ_sl + γ_lv cos θ. θ < 90° wetting (hydrophilic for water); > 90° non-wetting (hydrophobic). Lotus effect (Wilhelm Barthlott, Bonn) — superhydrophobic surfaces (θ > 150°) from hierarchical roughness + low surface energy; commercial: Lotusan paint (Sto), Nikon Crystal Coat.

Capillary rise: h = 2 γ cos θ / (ρ g r) — water in a 1 mm glass capillary rises ~15 mm. Drives groundwater movement, plant xylem transport (with Münch pressure flow), oil recovery, microfluidics.

Surface tension & interfacial phenomena

Surfactants self-assemble at interfaces (CMC, critical micelle concentration); HLB scale 1-20 (water vs oil affinity). Industrial applications: Henkel Persil (anionic + nonionic blend), Pfizer formulations, Unilever consumer products. Langmuir-Blodgett films (single-layer assemblies), self-assembled monolayers (alkanethiols on Au — Whitesides, Harvard).

12. Transport phenomena in physical chemistry

Transport — of mass, momentum, energy, and charge — connects molecular dynamics to macroscopic measurement.

Diffusion

Fick’s first law: J = -D ∇c. Fick’s second law: ∂c/∂t = D ∇²c. Stokes-Einstein relation for spherical particle in viscous medium:

D = k_B T / (6 π η r)

A 1 nm protein in water at 25 °C (η = 0.89 mPa·s): D ≈ 2 × 10⁻¹⁰ m²/s. Glucose in water D ≈ 6.7 × 10⁻¹⁰ m²/s. Self-diffusion of H₂O ≈ 2.3 × 10⁻⁹ m²/s. Diffusion through air ~10⁵× faster than through water for same species. Knudsen diffusion dominates in nanopores when mean free path > pore diameter — important for catalyst pellets and porous-membrane separations.

Viscosity and rheology

Newtonian fluids: τ = η (dv/dy). Water at 25 °C: η = 0.89 mPa·s. Honey: 2-10 Pa·s. Glycerol: 1.4 Pa·s. Polymer melts and solutions are typically non-Newtonian — shear thinning (most polymer melts, paint, ketchup, blood at low shear), shear thickening (cornstarch suspensions, “Oobleck”), Bingham plastic (toothpaste, drilling mud). Power-law model τ = K (du/dy)^n.

The Einstein viscosity equation: η = η_solvent (1 + 2.5 φ) for dilute hard-sphere suspensions, volume fraction φ. Useful starting point for nano-fluids, colloids.

Heat conduction

Fourier’s law: q = -k ∇T. Thermal conductivity at 25 °C: Ag 429, Cu 401, Al 237, Si 149, Fe 80, stainless 304 16, glass 1.05, water 0.6, air 0.026 (all W/(m·K)). Phonon transport in solids; molecular kinetic theory for gases (k_gas ≈ (1/3) c̄ ℓ ρ c_V); Wiedemann-Franz law for metals (k/σ = LT, Lorenz number L ≈ 2.44 × 10⁻⁸ V²/K²).

Onsager reciprocal relations

For coupled transport (heat and mass, charge and mass, etc.), Onsager (Nobel 1968) showed the cross-coefficient matrix is symmetric near equilibrium: L_ij = L_ji. Underlies thermoelectric effects (Seebeck, Peltier), thermo-diffusion (Soret), electrokinetic phenomena (streaming potential, electroosmosis).

13. Photochemistry and photophysics

Photochemistry: absorption of light drives a chemical change. Photophysics: relaxation without chemistry (fluorescence, phosphorescence, IC, ISC, energy transfer).

Selection rules and excited-state dynamics

Allowed electronic transitions: ΔS = 0 (spin), Δl = ±1 (orbital), Δm_l = 0, ±1. Forbidden transitions become weakly allowed via spin-orbit coupling (heavy-atom effect — explains why Ir(ppy)₃ phosphoresces brightly), vibronic coupling, symmetry breaking.

Excited-state lifetimes:

  • Internal conversion (IC, same multiplicity): ps - ns
  • Fluorescence S₁ → S₀: ns (typical 1-10 ns)
  • Intersystem crossing (ISC, S → T): ns - μs
  • Phosphorescence T₁ → S₀: μs - s (long because spin-forbidden)
  • TADF (thermally activated delayed fluorescence; Adachi 2012): exploits small ΔE(S₁-T₁); used in OLED displays without precious metals — Cynora, Kyulux, Idemitsu Kosan

Photocatalysis

  • TiO₂ (anatase, rutile, P25 mixed) under UV: organic pollutant degradation, antimicrobial coatings, self-cleaning glass (Pilkington Activ); discovered Honda-Fujishima 1972 for water splitting on TiO₂ under UV
  • Visible-light organic photoredox: Ru(bpy)₃²⁺, Ir(ppy)₃, organic acridinium and 4CzIPN — David MacMillan (Nobel 2021 for organocatalysis, also pioneered metallaphotoredox), Tehshik Yoon, Corey Stephenson
  • CO₂ reduction: photoelectrochemical → solar fuels; key open problem of energy chemistry

Solar cells (photovoltaic) physics

  • Crystalline silicon: Shockley-Queisser limit ~33% for single-junction (1.34 eV); real c-Si modules 21-24%; LONGi HPBC 26.81% record cell 2024; TOPCon and HJT replacing PERC
  • Multi-junction III-V: NREL/Spectrolab tandems exceed 47% under concentration
  • Perovskite: methylammonium lead iodide (CH₃NH₃PbI₃) and mixed-cation analogs reached 26.7% (UNIST 2024); stability remains the engineering challenge; Oxford PV silicon-perovskite tandem 28.6%

Surveyed in photovoltaic-systems.

14. Coordination chemistry as physical chemistry

Crystal-field and ligand-field theory predict d-orbital splitting, magnetism, color, and spectroscopy of transition-metal complexes.

  • Octahedral field: t₂g (xy, xz, yz) lower energy + e_g (z², x²-y²) higher; Δ_o (10 Dq) splitting
  • Tetrahedral: e (lower) + t₂ (upper); Δ_t = (4/9) Δ_o
  • Spectrochemical series: I⁻ < Br⁻ < Cl⁻ < F⁻ < OH⁻ < H₂O < NH₃ < en < CN⁻ < CO (weak → strong field)
  • High-spin / low-spin: based on Δ vs P (pairing energy); Fe²⁺(H₂O)₆ high-spin (paramagnetic), Fe²⁺(CN)₆⁴⁻ low-spin (diamagnetic)
  • Tanabe-Sugano diagrams for d-d transitions in different field strengths

Color: hexaaquacopper(II) blue (d-d absorption around 800 nm transmits blue); KMnO₄ purple (intense LMCT — ligand-to-metal charge transfer, ε ≈ 2400 at 525 nm); ferrocene yellow-orange; Ru(bpy)₃²⁺ orange (MLCT). UV-Vis spectroscopy is a routine diagnostic of metal-complex electronic structure.

See inorganic-chemistry for in-depth coordination chemistry, organometallics, and catalysis.

15. Molecular simulation methods

Beyond the static electronic-structure calculations of section 8, dynamic simulation is the second pillar of computational physical chemistry.

Molecular dynamics (MD)

Integrate Newton’s equations F = -∇V for ~10³-10⁹ atoms with classical force fields (additive pair potentials + bonded terms). Timestep 1-2 fs (constrained bonds via SHAKE / LINCS); routine timescales ns-μs on GPUs, ms on dedicated hardware (Anton 3 at D. E. Shaw Research).

Force fields:

  • Biomolecular: AMBER (ff14SB, ff19SB; Case et al.), CHARMM (c36m; Karplus, MacKerell — Martin Karplus shared Nobel 2013 with Michael Levitt and Arieh Warshel for multi-scale models), OPLS-AA (Jorgensen), GROMOS (van Gunsteren)
  • Materials: ReaxFF (van Duin — bond-breaking), COMPASS, UFF, PCFF, EAM/MEAM (metals), Tersoff/Brenner/AIREBO (carbon), LJ models (rare gases, generic studies)
  • Water models: TIP3P (simplest, fast), TIP4P/2005 (best balance), SPC/E, OPC, polarizable AMOEBA, MB-pol (many-body water)
  • Machine-learned potentials: ANI (Roitberg/Isayev), SchNet, NequIP/Allegro (Boris Kozinsky MIT/Harvard), MACE (Csányi Cambridge), GAP, DeePMD — DFT accuracy at near-classical cost; emerging for catalysis, electrolytes, materials discovery

Software: GROMACS (Hess, Lindahl, KTH; free), AMBER (commercial+free), NAMD (Schulten/Tajkhorshid UIUC; free), LAMMPS (Plimpton Sandia; free; materials-focused), OpenMM (Pande/Eastman Stanford; GPU-native, used by Folding@home), Desmond (Schrödinger LLC), ASE (Atomic Simulation Environment, Python).

Enhanced sampling

Plain MD samples a few ns-μs around the starting state; biological events (protein folding, ligand binding, conformational changes) take ms-s. Enhanced sampling methods bias the sampling to overcome free-energy barriers:

  • Umbrella sampling + WHAM analysis — free energy along chosen reaction coordinate
  • Metadynamics (Laio & Parrinello 2002) — history-dependent biasing potential fills wells
  • Replica-exchange MD (REMD/T-REMD) — parallel copies at different temperatures swap configurations
  • Steered MD, accelerated MD, well-tempered metadynamics, conformational flooding

Monte Carlo

Markov-chain Monte Carlo sampling of the Boltzmann distribution; Metropolis-Hastings acceptance. Useful for lattice models (Ising — magnetic, lattice gas), polymer configurations (Rosenbluth-style), and grand canonical adsorption studies (GCMC; predicts MOF/zeolite gas uptake).

Multiscale models

QM/MM (Warshel-Levitt 1976; Nobel 2013): quantum-mechanical region (active site, ~30-100 atoms with DFT or HF/MP2) embedded in classical environment (rest of protein, solvent). Standard for enzyme reaction mechanism studies; software ChemShell, ONIOM (Gaussian), Q-Chem QM/MM, Amber QM/MM.

Coarse-grained: Martini (Marrink Groningen — 4 heavy atoms ≈ 1 bead), SDK CG water, dissipative particle dynamics (DPD), implicit-solvent Gō-models for protein folding.

16. Worked numerical anchors

Brief illustrative numbers students and practitioners can carry in their heads.

  • k_B T at 298 K: 4.114 × 10⁻²¹ J = 25.7 meV; convenient for comparing to bond energies (eV) or photon energies (UV-Vis 1.5-6 eV, IR 0.05-0.5 eV).
  • RT at 298 K: 2.479 kJ/mol = 0.5924 kcal/mol; “kT per mole” benchmark for chemistry energies.
  • PV = nRT: 1 mol gas at 1 bar, 273 K → 22.71 L (molar volume); at STP (1 bar, 25 °C) → 24.79 L. (Avoid stale 22.4 L value — uses outdated 101.325 kPa STP.)
  • Henderson-Hasselbalch sanity check: pH = pK_a when [A⁻] = [HA]; 1 unit above pK_a means 10:1 conjugate base.
  • Diffusion length L = √(D t): protein in cytoplasm (D ≈ 10⁻¹¹ m²/s) covers 4 μm in 1 s; ion in water (D ≈ 10⁻⁹ m²/s) covers 30 μm in 1 s; gas in air (D ≈ 10⁻⁵ m²/s) covers 3 mm in 1 s.
  • Wavelength to energy: λ (nm) × E (eV) = 1240; so 500 nm ≈ 2.48 eV; 1240 nm = 1.0 eV (silicon band-gap edge near 1100 nm).
  • NMR field: 1 T → 42.58 MHz for ¹H. A 600 MHz spectrometer has B₀ = 14.09 T.
  • Beer’s law sanity check: for ε = 10,000 M⁻¹cm⁻¹ (typical aromatic π-π*) and a 1 cm path, [analyte] = 10⁻⁴ M gives A = 1 → I/I₀ = 0.1.

17. Connections forward

Physical chemistry underlies every applied chemistry domain:

Adjacent

Graph View

Backlinks