Physical Chemistry · Part 8 of 9 · Free

Molecular Quantum Mechanics — formula sheet

Every key expression and definition from Quantum Chemistry, Part 8, on one page. Free to read, no sign-in.

Key expressions

the exact molecular Hamiltonian
Ĥ = −∑_A (1/2M_A)∇²_A − ∑_i ½∇²_i − ∑_i,A Z_A/r_iA + ∑_i
the electronic Hamiltonian at fixed nuclei
Ĥ_el(R) = −∑_i ½∇²_i − ∑_i,A Z_A/r_iA + ∑_i
the electronic Schrödinger equation
Ĥ_el(R) ψ_el(r; R) = E_el(R) ψ_el(r; R)
the potential-energy surface
U(R) = E_el(R) + ∑_A
the nuclear Schrödinger equation
[−∑_A (1/2M_A)∇²_A + U(R)] χ(R) = E_total χ(R)
the Born–Huang expansion (exact)
Ψ(r, R) = ∑_k χ_k(R) ψ_k(r; R)
the exact coupled nuclear equations
[T̂_n + U_j(R) + Λ_jj(R)] χ_j + ∑_k≠j Λ_jk χ_k = E χ_j
the non-adiabatic coupling operators
Λ_jk = −∑_A (1/M_A) ⟨ψ_j|∇_A|ψ_k⟩ · ∇_A − ∑_A (1/2M_A) ⟨ψ_j|∇²_A|ψ_k⟩
the derivative coupling, showing the energy denominator
⟨ψ_j|∇_A|ψ_k⟩ = ⟨ψ_j|(∇_AĤ_el)|ψ_k⟩ / (E_k − E_j)
the H_2^+ electronic Hamiltonian
Ĥ_el = −½∇² − 1/r_A − 1/r_B + 1/R
the LCAO trial function
ψ = c_A 1s_A + c_B 1s_B
the secular equations for a two-orbital problem
(H_AA − ES_AA)c_A + (H_AB − ES_AB)c_B = 0(H_BA − ES_BA)c_A + (H_BB − ES_BB)c_B = 0
the secular determinant
| H_AA − E H_AB − ES || H_AB − ES H_BB − E | = 0
the two LCAO roots of a homonuclear diatomic
E_+ = (α + β)/(1 + S) E_− = (α − β)/(1 − S)
the bonding and antibonding molecular orbitals
σ_g = (1s_A + 1s_B) / √[2(1 + S)] σ_u* = (1s_A − 1s_B) / √[2(1 − S)]
the three two-centre integrals of H_2^+
S(R) = e^−R(1 + R + R²/3)J(R) = (1/R)[1 − (1 + R)e^−2R]K(R) = e^−R(1 + R)
α and β for H_2^+
α = H_AA = E_1s + 1/R − J β = H_AB = E_1sS + S/R − K
the H_2^+ LCAO energies
E_± = E_1s + 1/R − (J ± K)/(1 ± S)
the bonding density, showing the interference term
|σ_g|² = [1s_A² + 1s_B² + 2(1s_A)(1s_B)] / [2(1 + S)]
confocal elliptic coordinates
ξ = (r_A + r_B)/R ∈ [1, ∞) η = (r_A − r_B)/R ∈ [−1, 1] φ ∈ [0, 2π)
the virial theorem at a stationary point of the PES
⟨T⟩ = −E , ⟨V⟩ = 2E , ⟨T⟩/⟨V⟩ = −½
the general LCAO expansion
φ_i = ∑_μ=1^m c_μi χ_μ
the secular problem in matrix form
Hc_i = ε_i S c_i equivalently det|H − εS| = 0
the second-order interaction energy of two orbitals
ΔE^(2) ≈ |H_μν − εS_μν|² / (ε_μ − ε_ν)
the two-state eigenvalues
E_±(Q) = ½(H_11 + H_22) ± √[¼(H_11 − H_22)² + H_12²]
the Landau–Zener transition probability
P_hop = exp[−2πH_12² / (ℏ v |d(H_11 − H_22)/dQ|)]
the two second-row orderings
Li_2 → N_2: σ_g2s < σ_u*2s < π_u2p < σ_g2p < π_g*2p < σ_u*2pO_2, F_2: σ_g2s < σ_u*2s < σ_g2p < π_u2p < π_g*2p < σ_u*2p
Koopmans' theorem
IE_i ≈ −ε_i (Koopmans' theorem — frozen orbitals, no relaxation, no correlation change)
the heteronuclear two-orbital energies
E_± = ½(α_A + α_B) ± √[¼(α_A − α_B)² + β²]
two-configuration CI for H_2
Ψ = c_1|σ_g²⟩ + c_2|σ_u*²⟩
the Coulson–Fischer wavefunction
φ_a = 1s_A + λ1s_B , φ_b = 1s_B + λ1s_A , Ψ = φ_a(1)φ_b(2) + φ_b(1)φ_a(2)
the four sp³ hybrids
h_1 = ½(s + p_x + p_y + p_z)h_2 = ½(s + p_x − p_y − p_z)h_3 = ½(s − p_x + p_y − p_z)h_4 = ½(s − p_x − p_y + p_z)
the hybridisation–angle relation
⟨h_1|h_2⟩ = a² + b² cosθ = 0 ⇒ cosθ = −a²/b² = −1/λ
Bent's rule
Bent's rule: s character accumulates in hybrids directed towards ELECTROPOSITIVE substituents and lone pairs; p character accumulates towards ELECTRONEGATIVE substituents.
the unitary invariance of a single determinant
A determinantal wavefunction is INVARIANT (up to a phase) under any unitary transformation of its occupied orbitals among themselves.
the Hückel secular equation for atom r
(α − E)c_r + β ∑_s bonded to r c_s = 0
the Hückel problem is an adjacency-matrix eigenvalue problem
−x c_r + ∑_s bonded to r c_s = 0 ⇒ A c = x c
the ethene secular determinant
| −x 1 || 1 −x | = x² − 1 = 0 ⇒ x = 1, −1
the π energy of ethene — the universal reference
E_π(ethene) = 2(α + β) = 2α + 2.0000β
the Coulson–Rushbrooke pairing theorem
If x is an eigenvalue of an alternant system, so is −x; and the two eigenvectors differ only by a sign change on every starred atom.
the closed-form solution for a linear polyene
x_j = 2cos(jπ/(n+1)) , c_jr = √(2/(n+1)) sin(jrπ/(n+1)) , j, r = 1 … n
the π energy of butadiene
E_π(butadiene) = 2(α + 1.6180β) + 2(α + 0.6180β) = 4α + 4.4721β
the delocalisation energy of butadiene
DE(butadiene) = (4α + 4.4721β) − 2 × (2α + 2β) = 0.4721β
the π energy of benzene
E_π(benzene) = 2(α + 2β) + 4(α + β) = 6α + 8.0000β
the delocalisation (resonance) energy of benzene
DE(benzene) = (6α + 8.0000β) − 3 × (2α + 2β) = 2.0000β
the π-electron charge density on atom r
q_r = ∑_k n_k c_kr²
the π bond order between atoms r and s
p_rs = ∑_k n_k c_kr c_ks
the free valence of atom r
F_r = √3 − ∑_s p_rs
the bond order – bond length relation (schematic form)
R(p) = R_single − (R_single − R_double) × f(p)
the closed-form solution for a monocyclic polyene
x_k = 2cos(2πk/n) , k = 0, ±1, ±2, … E_k = α + 2βcos(2πk/n)
Hückel's rule, derived from the degeneracy pattern
closed π shell ⇔ 4n + 2 π electrons ⇔ AROMATIC
the Hückel heteroatom parameters
α_X = α + h_Xβ β_CX = k_CXβ
the Möbius (antiperiodic) eigenvalues
x_k = 2cos[(2k + 1)π/n]
the extended Hückel secular problem
Ĥc_i = ε_iSc_i — a GENERALISED eigenvalue problem, because S is no longer the identity

Definitions worth memorising

Born–Oppenheimer approximation: the total molecular wavefunction is written as a product Ψ(r, R) ≈ ψ_el(r; R) χ(R), in which the electronic wavefunction depends on the nuclear coordinates R only parametrically. Equivalently: the nuclear kinetic-energy operator is neglected when the electronic problem is solved, and its eigenvalue E_el(R) then acts as the potential energy for nuclear motion.
Potential-energy surface (PES): the total energy of a molecule — electronic energy plus nuclear repulsion — as a function of nuclear geometry, with the electrons in a specified state (usually the ground state). It has 3N − 6 dimensions for a non-linear molecule of N atoms, 3N − 5 if linear.
LCAO (Linear Combination of Atomic Orbitals): the approximation that a molecular orbital can be written as a weighted sum of atomic orbitals centred on the constituent atoms. It is a basis-set approximation, not a physical claim: the atomic orbitals are simply a convenient, chemically meaningful set of functions in which to expand the unknown MO.
Conical intersection: a point (strictly, a (3N − 8)-dimensional seam) at which two Born–Oppenheimer surfaces of the same spin multiplicity are exactly degenerate, and in whose vicinity the electronic and nuclear motions cannot be separated. Population transfers between the surfaces on a timescale of tens of femtoseconds — faster than a single vibration.
Jahn–Teller theorem: any non-linear molecule in a spatially degenerate electronic state is unstable with respect to some non-totally-symmetric distortion that removes the degeneracy. The theorem states that such a distortion exists; it does not predict which one, or how large.
Bond order: BO = ½(number of electrons in bonding MOs − number in antibonding MOs). It is a bookkeeping device, not an observable, but it correlates strongly and monotonically with bond length (shorter for higher BO) and with dissociation energy (larger for higher BO).
The heteronuclear polarity rule: the bonding MO has the larger coefficient on the more electronegative (lower-energy) atom, and the antibonding MO has the larger coefficient on the less electronegative atom. Electron density in the occupied bonding orbital therefore sits preferentially on the electronegative atom — which is the bond dipole, derived rather than assumed.
Hybrid orbital: a normalised linear combination of atomic orbitals on the same atom, chosen so that it points along a bond direction. Hybridisation is a unitary transformation within one atom's valence space. It changes no energy, creates no new orbital, and describes no physical process.
α (the Coulomb integral): H_rr = ⟨φ_r|Ĥ^eff|φ_r⟩, roughly the energy of an electron in an isolated 2p orbital on carbon. It is negative (a bound electron) and is the zero of the energy scale.
β (the resonance integral): H_rs = ⟨φ_r|Ĥ^eff|φ_s⟩ for bonded r and s. It is negative and is the interaction energy. Both are treated as empirical parameters and are never evaluated; results are quoted in units of β.
Alternant hydrocarbon: a conjugated system whose atoms can be divided into two sets (‘starred’ and ‘unstarred’) such that no two members of the same set are bonded to each other. Equivalently, the molecular graph is bipartite — it contains no odd-membered ring. Benzene, butadiene, naphthalene and allyl are alternant; azulene and fulvene are not.
Delocalisation energy (DE): the difference between the actual Hückel π energy and the π energy of a reference structure built from isolated, localised double bonds. DE = E_π(actual) − E_π(localised reference). Because β is negative, a negative DE coefficient means stabilisation.
Frost–Musulin circle: inscribe the regular n-gon of the ring in a circle of radius 2|β| centred at α, with one vertex pointing down. The vertical position of each vertex is the energy of one MO. It is not a mnemonic but an exact geometric statement of x_k = 2cos(2πk/n).

Where these come from

This sheet is distilled from Quantum Chemistry, Part 8 — 9 sections that derive every one of these results and show you how to use them.

Read Part 8 All formula sheets