Coulson's Directionality Theorem

Coulson's directionality theorem is a fundamental quantum-chemical relationship formulated by Charles Alfred Coulson in 1952 that links the interorbital angle θ_ij between any two orthogonal s–p hybrid atomic orbitals centered on the same nucleus directly to their hybridization mixing parameters λ_i and λ_j. In Linus Pauling's valence bond theory of orbital hybridization, a main-group central atom mixes its valence ns orbital with its three mutually perpendicular np_x, np_y, and np_z orbitals to construct directional hybrid wavefunctions that maximize overlap with surrounding ligand orbitals. Whereas idealized integer hybrids (sp, sp², and sp³) restrict interorbital angles to fixed values of 180°, 120°, and 109.47°, real molecules containing inequivalent ligands or lone pairs—as governed by VSEPR theory and Bent's rule—exhibit continuous, non-integer hybridization indices.

By enforcing quantum-mechanical normalization and mutual orthogonality ⟨ψ_i | ψ_j⟩ = δ_ij across all valence hybrid orbitals on a central atom, Coulson derived the compact analytical condition:

1 + λ_i λ_j cos(θ_ij) = 0

This equation allows chemists to calculate the fractional s character f_s and fractional p character f_p of bonding and lone-pair orbitals directly from experimentally measured bond angles (provided that the hybrid orbitals follow the internuclear axes without bond bending), or conversely to quantify the angular strain ("banana-bond" curvature) in small alicyclic rings such as cyclopropane.

Linear Combination of Atomic Orbitals and the Hybridization Index

Consider two normalized spⁿ hybrid orbitals ψ_i and ψ_j centered on the same atom A. Neglecting d-orbital participation, each hybrid orbital ψ_i is expressed as a linear combination of the spherically symmetric valence s orbital and a single directional p orbital p_i pointing along unit vector u_i = (sin θ_i cos φ_i, sin θ_i sin φ_i, cos θ_i):

ψ_i = c_s,i s + c_p,i p_i = (s + λ_i p_i) / √(1 + λ_i²)

Here, λ_i = c_p,i / c_s,i > 0 is the amplitude mixing coefficient, and its square n_i = λ_i² is the hybridization index appearing in the standard notation sp^n_i (meaning that the orbital possesses n_i parts p character for every 1 part s character). From the normalization condition ⟨ψ_i | ψ_i⟩ = c_s,i² + c_p,i² = 1, the fractional s character f_s,i and fractional p character f_p,i of hybrid ψ_i are given by:

f_s,i = 1 / (1 + λ_i²) = 1 / (1 + n_i), f_p,i = λ_i² / (1 + λ_i²) = n_i / (1 + n_i)

Note that some textbooks define the hybridization index directly as L_i = λ_i² so that ψ_i = (s + √L_i p_i) / √(1 + L_i), in which case Coulson's theorem is written 1 + √(L_i L_j) cos(θ_ij) = 0.

Orthogonality Integral of Directional p Orbitals

Any two distinct valence hybrid orbitals ψ_i and ψ_j on the same atom must be mutually orthogonal, ⟨ψ_i | ψ_j⟩ = 0. Expanding the inner product yields:

⟨ψ_i | ψ_j⟩ = [ ⟨s | s⟩ + λ_i ⟨p_i | s⟩ + λ_j ⟨s | p_j⟩ + λ_i λ_j ⟨p_i | p_j⟩ ] / [ √(1 + λ_i²) √(1 + λ_j²) ] = 0

Because the atomic s orbital is normalized (⟨s | s⟩ = 1) and orthogonal to any p orbital on the same nucleus (⟨s | p_i⟩ = ⟨s | p_j⟩ = 0), the cross terms vanish. Furthermore, decomposing the directional p orbitals into the Cartesian basis (p_i = u_ix p_x + u_iy p_y + u_iz p_z) and using the orthonormal relation ⟨p_α | p_β⟩ = δ_αβ shows that the overlap between two directional p orbitals equals the scalar dot product of their unit direction vectors, which is the cosine of the interorbital angle θ_ij:

⟨p_i | p_j⟩ = u_i · u_j = cos(θ_ij)

Substituting ⟨s | s⟩ = 1 and ⟨p_i | p_j⟩ = cos(θ_ij) into the orthogonality condition yields Coulson's directionality theorem:

1 + λ_i λ_j cos(θ_ij) = 0 ⇒ cos(θ_ij) = −1 / (λ_i λ_j) = −1 / √(n_i n_j)

Because λ_i and λ_j are positive real numbers, −1 / (λ_i λ_j) is always negative, demonstrating that the interorbital angle θ_ij between any two real, orthogonal spⁿ hybrid orbitals must satisfy 90° ≤ θ_ij ≤ 180°. No pair of orthogonal s–p hybrid orbitals can ever enclose an acute angle smaller than 90°.

Formulas for Equivalent Bonding Hybrids

When two or more hybrid orbitals ψ_i and ψ_j are symmetry-equivalent (λ_i = λ_j = λ, and hybridization index n = λ²), Coulson's theorem simplifies to:

cos(θ) = −1 / λ² = −1 / n ⇒ n = −1 / cos(θ)

Substituting n = −1 / cos(θ) into the expressions for fractional s and p character provides direct geometric formulas for equivalent hybrid orbitals:

f_s = cos(θ) / [cos(θ) − 1], f_p = 1 / [1 − cos(θ)]

In addition, because the transformation from the canonical {s, p_x, p_y, p_z} atomic basis to the four localized hybrid orbitals {ψ₁, ψ₂, ψ₃, ψ₄} is unitary, the fractional s and p characters across all four valence hybrids on a central atom must obey the conservation sum rules:

∑_i=1⁴ f_s,i = ∑_i=1⁴ [1 / (1 + λ_i²)] = 1, ∑_i=1⁴ f_p,i = ∑_i=1⁴ [λ_i² / (1 + λ_i²)] = 3

Application to Water (H₂O) and Ammonia (NH₃)

Assuming that the O–H bonding hybrids in water (H₂O) point directly along the internuclear axes at the experimental bond angle θ_HOH = 104.45°, cos(104.45°) = −0.2495 ≈ −0.250. Applying Coulson's formula for two equivalent bonding hybrids gives n_b = −1 / (−0.250) = 4.00, meaning that each O–H bonding orbital on oxygen is sp^4.00 hybridized with f_s,b = 0.200 (20.0% 2s character and 80.0% 2p character).

By the s-character sum rule (2 f_s,b + 2 f_s,LP = 1), the two equivalent lone-pair hybrids on oxygen must share the remaining 60.0% of the oxygen 2s orbital, giving each lone pair f_s,LP = 0.300 (30.0% 2s character and 70.0% 2p character), which corresponds to sp^2.33 hybridization (λ_LP² = 0.70/0.30 = 2.333). Substituting λ_LP² = 2.333 back into Coulson's theorem predicts an inter-lone-pair angle of θ_LP–O–LP = arccos(−1 / 2.333) = 115.4°, in quantitative agreement with Bent's rule (which states that atomic s character concentrates in lone pairs and electropositive bonds) and with the widened lone-pair domains observed in the Electron Localization Function.

Similarly, in ammonia (NH₃, θ_HNH = 107.8°), cos(107.8°) = −0.3057, yielding n_b = 3.27 (sp^3.27 bonding hybrids with f_s,b = 23.4% s character each). The single nitrogen lone pair therefore retains f_s,LP = 1 − 3(0.2342) = 0.2974 (29.7% s character, sp^2.36 hybridization). In heavy p-block hydrides such as H₂S (θ = 92.1°), cos(92.1°) = −0.0366, yielding sp^27.3 bonding orbitals (96.5% p character), confirming Drago's rule that heavy main-group atoms bond via nearly pure np orbitals while retaining the ns² pair in a stereochemically inert s orbital.

Cyclopropane and Interorbital vs. Internuclear Angles

Coulson and William E. Moffitt (1947–1949) applied the directionality theorem to resolve the paradox of cyclopropane (C₃H₆), whose equilateral carbon triangle enforces a geometric C–C–C internuclear angle of 60°. Because Coulson's theorem proves that no two real orthogonal spⁿ hybrid orbitals can form an interorbital angle below 90°, the carbon–carbon bonding hybrids in cyclopropane cannot point along the 60° internuclear axes.

In cyclopropane, the experimental H–C–H bond angle is widened to θ_HCH = 115.1°, which corresponds via cos(115.1°) = −0.4242 to sp^2.36 C–H bonding hybrids (f_s,CH = 0.298 each). By the s-character sum rule, the two equivalent C–C bonding hybrids on each carbon receive the remaining s character: f_s,CC = ½(1 − 2 × 0.298) = 0.202 (20.2% s character), corresponding to sp^3.95 (nearly sp⁴) hybridization. According to Coulson's theorem, two sp^3.95 hybrids enclose an interorbital angle of θ_ij = arccos(−1 / 3.95) = 104.7°. Thus, each C–C hybrid points outward by ½(104.7° − 60.0°) ≈ 22.3° away from the C–C internuclear axis, forming curved bent bonds ("banana bonds") whose outward-bowed ridge of maximum electron density is directly confirmed by X-ray deformation density maps and the (3, −1) bond paths of the Quantum Theory of Atoms in Molecules.

Why can the interorbital angle between two spⁿ hybrid orbitals never be less than 90°?

According to Coulson's directionality theorem, cos(θ_ij) = −1 / (λ_i λ_j). Because λ_i and λ_j are real, non-negative ratios of p to s orbital coefficients, their product λ_i λ_j is strictly positive, making cos(θ_ij) ≤ 0. In the limit where both orbitals are pure p orbitals (λ_i, λ_j → ∞), cos(θ_ij) → 0 and θ_ij → 90°; any admixture of spherical s character makes cos(θ_ij) negative and increases θ_ij into the obtuse range between 90° and 180°.

What is the hybridization of the carbon orbitals in cyclopropane?

Based on the experimental H–C–H angle of 115.1°, the two C–H bonding orbitals on each carbon are sp^2.36 hybridized (29.8% s character each), while the two C–C ring-bonding orbitals are sp^3.95 hybridized (20.2% s character each). Because the sp^3.95 orbitals point at an interorbital angle of 104.7° rather than the 60° geometric ring angle, the C–C bonds are bent outward by ~22° ("banana bonds"), giving cyclopropane alkene-like π reactivity and high C–H stretching frequencies.

How does Coulson's directionality theorem connect to Bent's rule?

Bent's rule states that a central atom directs more s character toward electropositive substituents and lone pairs, and more p character toward electronegative substituents. Through Coulson's formula cos(θ) = f_s / (f_s − 1), increasing the p character (lowering f_s) of two X–A–X bonding orbitals directly decreases the interorbital bond angle θ toward 90°, quantitatively linking substituent electronegativity to molecular bond angles.

1. Coulson, C. A. (1961). Valence (2nd ed.). Oxford University Press. https://doi.org/10.1021/ed039pA607.1
2. Coulson, C. A., & Moffitt, W. E. (1949). "The properties of certain strained hydrocarbons." Philosophical Magazine, 40(300), 1–35. https://doi.org/10.1080/14786444908561208
3. Bent, H. A. (1961). "An appraisal of valence-bond structures and hybridization in compounds of the first-row elements." Chemical Reviews, 61(3), 275–311. https://doi.org/10.1021/cr60211a005
4. Weinhold, F., & Landis, C. R. (2005). Valency and Bonding: A Natural Bond Orbital Donor-Acceptor Perspective. Cambridge University Press. https://doi.org/10.1017/CBO9780511614569