Physical Chemistry · Part 6 of 9 · Free
Approximation Methods — formula sheet
Every key expression and definition from Quantum Chemistry, Part 6, on one page. Free to read, no sign-in.
Key expressions
the variation theorem (variational principle)
W[φ] = ⟨φ|Ĥ|φ⟩ / ⟨φ|φ⟩ ≥ E_0
W[φ] = ⟨φ|Ĥ|φ⟩ / ⟨φ|φ⟩ ≥ E_0
the completeness expansion the proof rests on
φ = ∑_n c_nψ_n , Ĥψ_n = E_nψ_n , E_0 ≤ E_1 ≤ E_2 ≤ …
φ = ∑_n c_nψ_n , Ĥψ_n = E_nψ_n , E_0 ≤ E_1 ≤ E_2 ≤ …
⟨φ|Ĥ|φ⟩ = ∑_n |c_n|² E_n , ⟨φ|φ⟩ = ∑_n |c_n|²
the variation theorem in its transparent form
W − E_0 = [∑_n|c_n|²E_n − E_0∑_n|c_n|²] / ∑_n|c_n|² = ∑_n|c_n|²(E_n − E_0) / ∑_n|c_n|²
W − E_0 = [∑_n|c_n|²E_n − E_0∑_n|c_n|²] / ∑_n|c_n|² = ∑_n|c_n|²(E_n − E_0) / ∑_n|c_n|²
the minimisation conditions
∂W/∂c_i = 0 for every i — the variational conditions
∂W/∂c_i = 0 for every i — the variational conditions
Ĥ = −½∇² − 1/r , E_0 = −½ E_h = −13.6057 eV
⟨φ|φ⟩ = 4π∫_0^∞ e^−2αr² r² dr = (π/2α)^3/2
⟨φ|−1/r|φ⟩ = −4π∫_0^∞ e^−2αr² r dr = −π/α
⟨φ|−½∇²|φ⟩ = (3α/2)(π/2α)^3/2
the variational energy of a Gaussian trial function for H
W(α) = 3α/2 − 2√(2α/π)
W(α) = 3α/2 − 2√(2α/π)
dW/dα = 3/2 − √(2/π)·α^−1/2 = 0 ⇒ √α = (2/3)√(2/π) ⇒ α_opt = 8/9π = 0.28294
the best a single Gaussian can do for the hydrogen atom
W_min = −4/3π = −0.424413 E_h = −11.5489 eV
W_min = −4/3π = −0.424413 E_h = −11.5489 eV
φ(x) = x(L − x)
⟨φ|φ⟩ = ∫_0^L x²(L−x)² dx = L⁵/30
⟨φ|Ĥ|φ⟩ = (ℏ²/2m)∫_0^L (dφ/dx)² dx = (ℏ²/2m)∫_0^L(L−2x)² dx = (ℏ²/2m)(L³/3)
variational energy of the trial function x(L−x)
W = (ℏ²/2m) (L³/3)/(L⁵/30) = 10ℏ²/(2mL²) = 5ℏ²/mL²
W = (ℏ²/2m) (L³/3)/(L⁵/30) = 10ℏ²/(2mL²) = 5ℏ²/mL²
the error in the variational energy is second order in the error in ψ
W = E_0 + ε²(⟨χ|Ĥ|χ⟩ − E_0) + O(ε³)
W = E_0 + ε²(⟨χ|Ĥ|χ⟩ − E_0) + O(ε³)
the scaling form of the variational energy for a Coulomb problem
W(η) = η²⟨T⟩_1 + η⟨V⟩_1
W(η) = η²⟨T⟩_1 + η⟨V⟩_1
the virial theorem, satisfied automatically at the variational optimum
2⟨T⟩ = −⟨V⟩ , E = −⟨T⟩ = ½⟨V⟩
2⟨T⟩ = −⟨V⟩ , E = −⟨T⟩ = ½⟨V⟩
the Hellmann–Feynman theorem
dE/dλ = ⟨ψ_λ|∂Ĥ/∂λ|ψ_λ⟩
dE/dλ = ⟨ψ_λ|∂Ĥ/∂λ|ψ_λ⟩
the two matrices of the linear variation method
H_ij = ⟨f_i|Ĥ|f_j⟩ (matrix elements), S_ij = ⟨f_i|f_j⟩ (overlap integrals)
H_ij = ⟨f_i|Ĥ|f_j⟩ (matrix elements), S_ij = ⟨f_i|f_j⟩ (overlap integrals)
W = (∑_i∑_j c_ic_jH_ij) / (∑_i∑_j c_ic_jS_ij)
W ∑_i∑_j c_ic_jS_ij = ∑_i∑_j c_ic_jH_ij
the secular equations
∑_i c_i(H_ki − W S_ki) = 0 for every k = 1, 2, …, n
∑_i c_i(H_ki − W S_ki) = 0 for every k = 1, 2, …, n
the secular determinant
det | H_ij − W S_ij | = 0
det | H_ij − W S_ij | = 0
| α−W β−WS | | β−WS α−W | = 0 ⇒ (α−W)² = (β−WS)²
the two roots of the symmetric 2×2 secular determinant
W_− = (α + β)/(1 + S) , W_+ = (α − β)/(1 − S)
W_− = (α + β)/(1 + S) , W_+ = (α − β)/(1 − S)
the linear variation method in matrix form
Hc = W Sc
Hc = W Sc
one-centre s-Gaussian integrals
S_ij = (π/p)^3/2 , T_ij = 3α_iα_jπ^3/2/p^5/2 , V_ij = −2πZ/p
S_ij = (π/p)^3/2 , T_ij = 3α_iα_jπ^3/2/p^5/2 , V_ij = −2πZ/p
the perturbation split
Ĥ = Ĥ^(0) + λĤ′
Ĥ = Ĥ^(0) + λĤ′
the Rayleigh–Schrödinger expansions
E_n = E_n^(0) + λE_n^(1) + λ²E_n^(2) + …ψ_n = ψ_n^(0) + λψ_n^(1) + λ²ψ_n^(2) + …
E_n = E_n^(0) + λE_n^(1) + λ²E_n^(2) + …ψ_n = ψ_n^(0) + λψ_n^(1) + λ²ψ_n^(2) + …
the zeroth- and first-order equations
λ^0: Ĥ^(0)ψ_n^(0) = E_n^(0)ψ_n^(0)λ^1: Ĥ^(0)ψ_n^(1) + Ĥ′ψ_n^(0) = E_n^(0)ψ_n^(1) + E_n^(1)ψ_n^(0)
λ^0: Ĥ^(0)ψ_n^(0) = E_n^(0)ψ_n^(0)λ^1: Ĥ^(0)ψ_n^(1) + Ĥ′ψ_n^(0) = E_n^(0)ψ_n^(1) + E_n^(1)ψ_n^(0)
the first-order energy correction
E_n^(1) = ⟨ψ_n^(0)|Ĥ′|ψ_n^(0)⟩ ≡ H′_nn
E_n^(1) = ⟨ψ_n^(0)|Ĥ′|ψ_n^(0)⟩ ≡ H′_nn
the first-order correction to the wavefunction
ψ_n^(1) = ∑_m≠n [H′_mn / (E_n^(0) − E_m^(0))] ψ_m^(0)
ψ_n^(1) = ∑_m≠n [H′_mn / (E_n^(0) − E_m^(0))] ψ_m^(0)
intermediate normalisation
⟨ψ_n^(0)|ψ_n⟩ = 1 ⇒ ⟨ψ_n^(0)|ψ_n^(k)⟩ = 0 for all k ≥ 1
⟨ψ_n^(0)|ψ_n⟩ = 1 ⇒ ⟨ψ_n^(0)|ψ_n^(k)⟩ = 0 for all k ≥ 1
the reduced resolvent of the unperturbed Hamiltonian
R̂_n = ∑_m≠n |ψ_m^(0)⟩⟨ψ_m^(0)| / (E_n^(0) − E_m^(0))
R̂_n = ∑_m≠n |ψ_m^(0)⟩⟨ψ_m^(0)| / (E_n^(0) − E_m^(0))
the Brillouin–Wigner second-order formula (self-consistent in E_n)
E_n = E_n^(0) + H′_nn + ∑_m≠n |H′_mn|² / (E_n − E_m^(0)) + …
E_n = E_n^(0) + H′_nn + ∑_m≠n |H′_mn|² / (E_n − E_m^(0)) + …
the second-order energy correction
E_n^(2) = ∑_m≠n |H′_mn|² / (E_n^(0) − E_m^(0))
E_n^(2) = ∑_m≠n |H′_mn|² / (E_n^(0) − E_m^(0))
the second-order correction to the ground state is never positive
E_0^(2) = ∑_m≠0 |H′_m0|² / (E_0^(0) − E_m^(0)) ≤ 0 always
E_0^(2) = ∑_m≠0 |H′_m0|² / (E_0^(0) − E_m^(0)) ≤ 0 always
the energy through second order — the formula to memorise
E_n ≈ E_n^(0) + H′_nn + ∑_m≠n |H′_mn|² / (E_n^(0) − E_m^(0))
E_n ≈ E_n^(0) + H′_nn + ∑_m≠n |H′_mn|² / (E_n^(0) − E_m^(0))
the Dalgarno–Lewis inhomogeneous equation
(Ĥ^(0) − E_n^(0))ψ_n^(1) = −(Ĥ′ − E_n^(1))ψ_n^(0)
(Ĥ^(0) − E_n^(0))ψ_n^(1) = −(Ĥ′ − E_n^(1))ψ_n^(0)
the Unsöld or average-energy (closure) approximation
E_n^(2) ≈ −(1/ΔE) [⟨n|Ĥ′²|n⟩ − ⟨n|Ĥ′|n⟩²]
E_n^(2) ≈ −(1/ΔE) [⟨n|Ĥ′²|n⟩ − ⟨n|Ĥ′|n⟩²]
second-order energy of a cubic anharmonicity
E_n^(2) = −(c²β⁶/ℏω)(30n² + 30n + 11)
E_n^(2) = −(c²β⁶/ℏω)(30n² + 30n + 11)
the standard spectroscopic vibrational term expansion
G(v) = ν̃_e(v + ½) − ν̃_ex_e(v + ½)² + …
G(v) = ν̃_e(v + ½) − ν̃_ex_e(v + ½)² + …
the secular determinant of degenerate perturbation theory
det | H′_ij − E^(1)δ_ij | = 0
det | H′_ij − E^(1)δ_ij | = 0
the Stark perturbation for a uniform field along z
Ĥ′ = eFz = eFr cosθ
Ĥ′ = eFz = eFr cosθ
the only surviving Stark matrix element in the n = 2 shell
⟨2s|z|2p_z⟩ = −3a_0 ⇒ H′_12 = −3ea_0F
⟨2s|z|2p_z⟩ = −3a_0 ⇒ H′_12 = −3ea_0F
the first-order Stark splitting of the hydrogen n = 2 shell
E^(1) = +3ea_0F, 0, 0, −3ea_0F
E^(1) = +3ea_0F, 0, 0, −3ea_0F
∑_j (H′_ij − E^(1)δ_ij) c_j = 0 for i = 1, …, g
the second-order effective Hamiltonian in a degenerate space
Ĥ^eff_ij = ∑_m∉D H′_imH′_mj / (E^(0)_D − E_m^(0))
Ĥ^eff_ij = ∑_m∉D H′_imH′_mj / (E^(0)_D − E_m^(0))
the time-dependent perturbation problem
Ĥ(t) = Ĥ^(0) + Ĥ′(t)
Ĥ(t) = Ĥ^(0) + Ĥ′(t)
the ansatz of time-dependent perturbation theory
Ψ(t) = ∑_n c_n(t) ψ_n e^−iE_nt/ℏ
Ψ(t) = ∑_n c_n(t) ψ_n e^−iE_nt/ℏ
the first-order transition amplitude
c_f^(1)(t) = (1/iℏ) ∫_0^t H′_fi(t′) e^iω_fit′ dt′ with ω_fi = (E_f − E_i)/ℏ
c_f^(1)(t) = (1/iℏ) ∫_0^t H′_fi(t′) e^iω_fit′ dt′ with ω_fi = (E_f − E_i)/ℏ
the first-order transition probability for a constant perturbation
|c_f(t)|² = (|H′_fi|²/ℏ²) · 4sin²(ω_fit/2)/ω_fi²
|c_f(t)|² = (|H′_fi|²/ℏ²) · 4sin²(ω_fit/2)/ω_fi²
Fermi's golden rule
W_i→f = (2π/ℏ) |H′_fi|² ρ(E_f)
W_i→f = (2π/ℏ) |H′_fi|² ρ(E_f)
the resonance condition from time-dependent perturbation theory
ℏω = E_f − E_i (the Bohr condition, recovered rather than assumed)
ℏω = E_f − E_i (the Bohr condition, recovered rather than assumed)
the electric-dipole interaction with a light wave
Ĥ′(t) = −μ̂·E_0 cosωt
Ĥ′(t) = −μ̂·E_0 cosωt
the transition dipole moment
μ_fi = ⟨ψ_f|μ̂|ψ_i⟩ , rate ∝ |μ_fi|²
μ_fi = ⟨ψ_f|μ̂|ψ_i⟩ , rate ∝ |μ_fi|²
the Einstein B coefficient from first-order time-dependent perturbation theory
B_fi = |μ_fi|² / (6ε_0ℏ²)
B_fi = |μ_fi|² / (6ε_0ℏ²)
the Einstein relation between spontaneous and stimulated emission
A = (8πhν³/c³) B
A = (8πhν³/c³) B
the oscillator strength of a transition
f_fi = (2m_eω_fi/3ℏe²) |μ_fi|²
f_fi = (2m_eω_fi/3ℏe²) |μ_fi|²
lifetime (natural) broadening
δE · τ ≈ ℏ , δν̃ / cm⁻¹ ≈ 5.31 / (τ / ps)
δE · τ ≈ ℏ , δν̃ / cm⁻¹ ≈ 5.31 / (τ / ps)
the helium-atom Hamiltonian (atomic units, Z = 2)
Ĥ = −½∇²_1 − ½∇²_2 − Z/r_1 − Z/r_2 + 1/r_12
Ĥ = −½∇²_1 − ½∇²_2 − Z/r_1 − Z/r_2 + 1/r_12
E^(0) = 2 × (−Z²/2) = −Z² = −4 E_h = −108.85 eV
ψ^(0)(1,2) = (Z³/π) e^−Zr_1 e^−Zr_2
the first-order electron-repulsion energy for a hydrogenic 1s² product
E^(1) = ⟨ψ^(0)|1/r_12|ψ^(0)⟩ = 5Z/8 = 5/4 E_h = 34.01 eV
E^(1) = ⟨ψ^(0)|1/r_12|ψ^(0)⟩ = 5Z/8 = 5/4 E_h = 34.01 eV
helium through first order
E ≈ E^(0) + E^(1) = −Z² + 5Z/8 = −2.75 E_h = −74.83 eV
E ≈ E^(0) + E^(1) = −Z² + 5Z/8 = −2.75 E_h = −74.83 eV
the one-parameter variational trial function for helium
φ(1,2) = (ζ³/π) e^−ζr_1 e^−ζr_2
φ(1,2) = (ζ³/π) e^−ζr_1 e^−ζr_2
the variational energy of helium as a function of the orbital exponent
W(ζ) = ζ² − 2Zζ + 5ζ/8
W(ζ) = ζ² − 2Zζ + 5ζ/8
the optimum effective nuclear charge for helium
ζ_opt = Z − 5/16 = 2 − 0.3125 = 27/16 = 1.68750
ζ_opt = Z − 5/16 = 2 − 0.3125 = 27/16 = 1.68750
the best energy from a single-exponent product function
W_min = −(Z − 5/16)² = −(27/16)² = −2.847656 E_h = −77.489 eV
W_min = −(Z − 5/16)² = −(27/16)² = −2.847656 E_h = −77.489 eV
1/r_12 = ∑_l (r_<^l/r_>^l+1) P_l(cosγ_12)
a Hylleraas-type explicitly correlated trial function
φ = e^−ζ(r_1+r_2) [1 + c_1r_12 + c_2(r_1−r_2)² + …]
φ = e^−ζ(r_1+r_2) [1 + c_1r_12 + c_2(r_1−r_2)² + …]
the local classical momentum
p(x) = √[2m(E − V(x))]
p(x) = √[2m(E − V(x))]
the WKB wavefunction in a classically allowed region
ψ(x) ≈ (C/√p(x)) exp[ ±(i/ℏ) ∫p(x) dx ]
ψ(x) ≈ (C/√p(x)) exp[ ±(i/ℏ) ∫p(x) dx ]
the WKB (Bohr–Sommerfeld) quantisation condition, two soft turning points
∫_x1^x_2 p(x) dx = (n + ½)πℏ , n = 0, 1, 2, …
∫_x1^x_2 p(x) dx = (n + ½)πℏ , n = 0, 1, 2, …
the WKB (Gamow) tunnelling probability
T ≈ exp[ −(2/ℏ) ∫_x1^x_2 √(2m(V(x) − E)) dx ]
T ≈ exp[ −(2/ℏ) ∫_x1^x_2 √(2m(V(x) − E)) dx ]
(S′)² − iℏS″ = 2m(E − V) = p²
the complete non-relativistic molecular Hamiltonian
Ĥ = T̂_n + T̂_e + V̂_ne + V̂_ee + V̂_nn
Ĥ = T̂_n + T̂_e + V̂_ne + V̂_ee + V̂_nn
the Born–Oppenheimer product ansatz
Ψ(r, R) ≈ ψ_el(r; R) χ_nuc(R)
Ψ(r, R) ≈ ψ_el(r; R) χ_nuc(R)
the Born–Oppenheimer expansion parameter
κ = (m_e/M)^1/4 = 0.1528 for a proton
κ = (m_e/M)^1/4 = 0.1528 for a proton
the non-adiabatic coupling terms, which the approximation discards
T̂_n(ψ_elχ) = ψ_elT̂_nχ + [terms in ∇_Rψ_el and ∇²_Rψ_el]
T̂_n(ψ_elχ) = ψ_elT̂_nχ + [terms in ∇_Rψ_el and ∇²_Rψ_el]
Definitions worth memorising
The variation theorem: for any normalisable, well-behaved function φ that satisfies the boundary conditions of the problem, the expectation value of the Hamiltonian is never less than the true ground-state energy: W[φ] = ⟨φ|Ĥ|φ⟩ / ⟨φ|φ⟩ ≥ E_0. Equality holds if and only if φ is itself the ground-state eigenfunction.
The symmetry-restricted variation theorem: if φ belongs to a definite irreducible representation Γ (or a definite eigenvalue of some operator commuting with Ĥ), then W[φ] ≥ E_0^Γ, the lowest exact eigenvalue of that same symmetry. This is what makes variational calculations of excited states possible at all, and it is why quantum-chemistry programs ask you for the state symmetry before they start.
The Hylleraas–Undheim–MacDonald theorem: if the linear variation problem is solved in a basis of n functions, the n roots W_1 ≤ W_2 ≤ … ≤ W_n satisfy W_k ≥ E_k−1 for every k. That is, the k-th root is an upper bound to the k-th exact eigenvalue, not just the lowest. Adding a further basis function pushes every root down or leaves it where it was, and the new roots interlace the old ones.
Trial function (or trial wavefunction): a normalisable function φ that obeys the boundary conditions of the problem and contains one or more adjustable variational parameters. The parameters are chosen so as to make W as small as possible; the resulting minimum is the best estimate the chosen functional form can give.
Linear variation method: take the trial function to be a linear combination of a fixed set of basis functions f_1, f_2, …, f_n, with the expansion coefficients as the variational parameters: φ = ∑_i c_if_i. The basis functions are fixed and known; only the coefficients are varied.
Secular equations: the set of n simultaneous linear homogeneous equations ∑_i(H_ki − W S_ki)c_i = 0 obtained by minimising the linear variational energy. Secular determinant: the determinant of the coefficient matrix, which must vanish for a non-trivial solution to exist.
Unperturbed Hamiltonian Ĥ^(0): an operator whose eigenfunctions ψ_n^(0) and eigenvalues E_n^(0) are known exactly. Perturbation Ĥ′: the remainder, assumed small in the sense that its matrix elements are small compared with the spacings of the unperturbed levels. λ: a dimensionless bookkeeping parameter, run from 0 to 1, whose only job is to let terms of different orders be separated.
First-order energy correction: the expectation value of the perturbation over the unperturbed wavefunction. It requires no new wavefunction, no sum over states and no knowledge of any other level — just one integral.
The second-order ground-state theorem: the second-order correction to the energy of the ground state is always negative or zero, whatever the perturbation is. It is zero only if the perturbation couples the ground state to nothing at all. There is no corresponding rule for excited states.
Degenerate perturbation theory: when the level of interest is g-fold degenerate, form the g×g matrix of the perturbation within that degenerate set, H′_ij = ⟨ψ_i^(0)|Ĥ′|ψ_j^(0)⟩, and diagonalise it. Its eigenvalues are the first-order energy corrections; its eigenvectors are the correct zeroth-order wavefunctions.
Fermi's golden rule: the first-order transition rate from an initial state into a continuum of final states is (2π/ℏ) times the squared coupling matrix element times the density of final states at the energy of the initial state. Energy conservation is not assumed; it emerges from the long-time limit.
Transition dipole moment: the matrix element of the dipole operator between the initial and final states. If it vanishes the transition is forbidden in the electric-dipole approximation; if it is non-zero the transition is allowed, and its intensity is proportional to the square of the moment. Every selection rule in spectroscopy is a statement about when this integral vanishes.
Effective nuclear charge Z_eff and screening constant σ: the charge an electron behaves as though it sees, Z_eff = Z − σ. The variational calculation for helium gives σ = 5/16 = 0.3125 exactly, meaning that one 1s electron screens the nucleus from the other by just under a third of a proton.
Electron correlation: the tendency of electrons to avoid one another beyond what the average (mean-field) repulsion accounts for. Correlation energy: the difference between the exact non-relativistic energy and the Hartree–Fock limit — that is, the best energy obtainable from a single determinant of orbitals. It is always negative.
Validity condition for WKB: the local de Broglie wavelength must change little over one wavelength: |dλ̄/dx| ≪ 1 with λ̄ = ℏ/p, equivalently mℏ|dV/dx| ≪ p³. The condition fails wherever p → 0 — at every classical turning point — and wherever V changes abruptly.
Born–Oppenheimer approximation: because nuclei are far heavier than electrons, the electronic wavefunction can be solved with the nuclei held fixed, with the nuclear positions entering only as parameters. The electronic energies as a function of nuclear geometry constitute the potential-energy surface, on which the nuclei then move as though in an ordinary potential.
Where these come from
This sheet is distilled from Quantum Chemistry, Part 6 — 10 sections that derive every one of these results and show you how to use them.
Read Part 6 All formula sheets