Chemical Bond

A chemical bond is a stabilizing physical attraction between two or more atoms, ions, or atomic cores that lowers the total Gibbs or electronic potential energy of the multi-atom assembly relative to separated neutral atoms, binding them into discrete molecules, polyatomic ions, or extended crystalline lattices. First conceptualized through Gilbert N. Lewis's electron-pair sharing model (1916) and Irving Langmuir's octet rule (1919), and given its rigorous quantum-mechanical foundation in 1927 by Walter Heitler and Fritz London (for H₂) alongside Friedrich Hund and Robert S. Mulliken (molecular orbital theory), chemical bonding arises from the interplay of Coulombic electrostatic interactions, the quantum kinetic energy of delocalized electrons, and the Pauli exclusion principle acting on the electronic wave function.

Within physical chemistry, chemical thermodynamics, and electronic structure theory—including photoelectron spectroscopy, the electron localization function (ELF), and Richard Bader's quantum theory of atoms in molecules (QTAIM)—primary chemical bonds span a continuous triad of electronic localization regimes (covalent, ionic, and metallic bonding, with bond dissociation enthalpies D₀ of 150 to 1,072 kJ mol⁻¹), complemented by weaker non-covalent intermolecular interactions (hydrogen bonds, dipole–dipole forces, and London dispersion forces of 1 to 40 kJ mol⁻¹) parameterized in molecular dynamics force fields.

The Morse Potential Curve (r_e, D_e) and Zero-Point Vibrational Energy

Under the Born–Oppenheimer approximation, clamping two nuclei A and B at internuclear distance R and solving the electronic Schrödinger equation yields the adiabatic potential energy curve V(R), modeled accurately across bound and dissociative regimes by Philip M. Morse's (1929) anharmonic potential:

where r_e is the equilibrium bond length (the bottom of the potential well where net force −*dV*/dR = 0), D_e is the spectroscopic well depth measured from the potential minimum to separated atoms, and a = (k_e / 2*D*_e)^1/2 controls the width of the well given harmonic force constant k_e. Because quantum nuclear motion cannot have zero momentum by the Heisenberg uncertainty principle, the lowest vibrational eigenstate (v = 0) sits above D_e by the zero-point vibrational energy E_ZPE = (1/2)ℏω₀, defining the experimental 0 K bond dissociation enthalpy: D₀ = D_e − (1/2)ℏω₀ (for H₂, r_e = 0.7414 Å, D_e = 458.1 kJ mol⁻¹, and D₀ = 432.1 kJ mol⁻¹).

Quantum Kinetic Energy Lowering and the Molecular Virial Theorem

A widespread misconception treats covalent bonding purely as electrostatic attraction between shared electrons and two positively charged nuclei. Detailed quantum-chemical energy decomposition (by Klaus Ruedenberg and John C. Slater) using the molecular virial theorem—which at the equilibrium bond length r_e (dE/dR = 0) enforces 2⟨*T*⟩ + ⟨*V*⟩ = 0 and Δ*E* = −Δ⟨*T*⟩ = (1/2)Δ⟨*V*⟩—reveals a deeper two-stage mechanism:

1. Delocalization Lowering of Quantum Kinetic Energy: As two H atoms approach from infinity, allowing each electron's wave function ψ to delocalize across both nuclei expands the effective spatial box length Δ*x*, reducing the curvature ∇²ψ and lowering the electronic quantum kinetic energy ⟨*T*⟩ (an interference effect driven by the cross term 2*S* ψ_Aψ_B).
2. Orbital Contraction Toward the Nuclei: In the final stage near r_e, the bonding orbital contracts radially around the two nuclei; this contraction raises electronic quantum kinetic energy Δ⟨*T*⟩ > 0 while lowering electrostatic electron–nucleus potential energy Δ⟨*V*⟩ = −2Δ⟨*T*⟩ < 0 by twice as much, satisfying the virial theorem at the energy minimum.

Heitler–London Spin Exchange, Pauling Hybridization (sp, sp², sp³), and Resonance

Valence Bond (VB) theory—initiated by Walter Heitler and Fritz London (1927) and expanded by John C. Slater and Linus Pauling (1931)—constructs the multi-electron wave function from localized atomic orbitals undergoing pairwise spin-singlet coupling (quantum entanglement of opposite electron spins, α(1)β(2) − β(1)α(2)):

Whereas the antisymmetric spatial triplet state Ψ₋ has a nodal plane midway between the nuclei and is purely repulsive, the symmetric spatial singlet state Ψ₊ concentrates probability density in the internuclear region via the quantum exchange integral J, stabilizing H₂. To explain polyatomic directional geometry (such as the 109.47° tetrahedral angles of methane, CH₄, the 120° trigonal planar geometry of ethene, C₂H₄, and the 180° linear geometry of ethyne, C₂H₂), Pauling introduced orbital hybridization (sp³, sp², sp linear combinations of valence 2*s* and 2*p* orbitals) and resonance superposition Ψ = c_covΨ_A–B + c_ion(Ψ_A+B− + Ψ_A−B+).

Linear Combination of Atomic Orbitals (LCAO-MO) and Formal Bond Order

Molecular Orbital (MO) theory—developed by Friedrich Hund, Robert S. Mulliken, and John Lennard-Jones—dispenses with localized electron pairs and builds canonical molecular orbitals ψ_i delocalized over the entire molecule as a Linear Combination of Atomic Orbitals (LCAO):

where S_AB = ⟨ϕ_A | ϕ_B⟩ is the overlap integral:

• Bonding Molecular Orbitals (σ, π): Constructive in-phase addition (ϕ_A + ϕ_B) lacks an internuclear nodal plane, lowering orbital energy below the parent atomic orbitals. Cylindrically symmetric head-on overlap forms sigma (σ) bonds, while side-by-side parallel p- or d-orbital overlap with a nodal plane containing the internuclear axis forms pi (π) bonds.
• Antibonding Molecular Orbitals (σ, π): Destructive out-of-phase subtraction (ϕ_A − ϕ_B) creates a nodal plane perpendicular to the bond axis between the nuclei, raising energy more strongly than the bonding orbital lowers it (ΔE^* > ΔE due to 1 − S_AB in the denominator).

Filling canonical MOs according to the Aufbau principle, Pauli exclusion, and Hund's rule defines the MO bond order:

MO theory effortlessly explains why He₂ (2 σ_1s electrons and 2 σ^*_1s electrons, bond order = 0) does not form a covalent bond, why N₂ (bond order = 3, D₀ = 945 kJ mol⁻¹, r_e = 1.098 Å) is exceptionally strong, and why ground-state molecular oxygen (O₂) (bond order = 2, with two unpaired electrons occupying degenerate π_2px and π_2py antibonding orbitals in a ³Σ_g⁻ triplet ground state) is paramagnetic—a milestone prediction validated by photoelectron spectroscopy.

Pauling Electronegativity (Δχ), Born–Haber Cycles, and Madelung Energy

When two bonded atoms A and B belong to different chemical elements, unequal nuclear effective charge (Z_eff) skews the shared electron cloud toward the more electronegative atom, generating a polar covalent bond with electric dipole moment μ = q_eff r_e. In 1932, Linus Pauling noticed that heteronuclear bonds A–B are nearly always stronger than the arithmetic or geometric mean of homonuclear A–A and B–B bonds by an ionic resonance stabilization Δ*D*_AB, defining the Pauling electronegativity scale (χ) (in eV^1/2, where F = 3.98, O = 3.44, Cl = 3.16, C = 2.55, H = 2.20, Cs = 0.79):

When Δχ exceeds ~1.7 to 2.0 (as between an low-ionization-energy alkali/alkaline-earth metal and a high-electron-affinity halogen or chalcogen), valence electron transfer forms a 3D crystalline ionic lattice stabilized by long-range Coulomb summation quantified by the Born–Landé / Kapustinskii lattice energy (U_L):

where M is the geometric Madelung constant (M = 1.74756 for rocksalt NaCl; 1.76267 for CsCl) and n ≈ 8–12 is the Born Pauli-repulsion exponent, determined experimentally via thermodynamic Born–Haber cycles.

Bader's QTAIM Bond Critical Points (3, −1) and ELF Basins

To classify chemical bonds rigorously from experimental X-ray charge densities or ab initio wave functions without arbitrary orbital choices, Richard F. W. Bader's Quantum Theory of Atoms in Molecules (QTAIM) analyzes the topology of the total electron density scalar field ρ(r). A chemical bond corresponds to a ridge of maximum electron density—a bond path—connecting two nuclear attractors through a saddle point where ∇ρ(r_c) = 0, termed a (3, −1) bond critical point (BCP). Evaluating the Laplacian of the electron density ∇²ρ(r_c) at the BCP separates shared-shell covalent bonds (where charge is concentrated in the internuclear region, ∇²ρ(r_c) < 0 and high disynaptic electron localization function basin population V(A,B)) from closed-shell ionic, metallic, and hydrogen bonds (where Pauli exclusion depletes charge at the boundary, ∇²ρ(r_c) > 0).

Does breaking a chemical bond release energy, or does forming a chemical bond release energy?

Breaking a stable chemical bond is always endothermic (absorbs energy, Δ*H* > 0), whereas forming a chemical bond is always exothermic (releases energy, Δ*H* < 0). Because bonded atoms sit in a potential energy well (D₀) below the energy of separated atoms, energy must be supplied to pull them apart against their attractive forces. Exothermic chemical and biochemical reactions—such as the combustion of methane (CH₄ + 2 O₂ → CO₂ + 2 H₂O) or the hydrolysis of ATP—release net free energy only because the new bonds formed in the products (and their hydration shells) are stronger and lower in energy than the bonds broken in the reactants.

Why does Molecular Orbital (MO) theory succeed where Lewis structures fail in explaining liquid oxygen (O₂)?

A standard Lewis dot structure for O₂ (:Ö=Ö:) pairs all 12 valence electrons into one σ bond, one π bond, and four lone pairs, predicting that O₂ should be diamagnetic (weakly repelled by a magnetic field). In reality, pale-blue liquid O₂ is strongly paramagnetic and hangs suspended between the poles of a magnet. Molecular Orbital theory explains this directly: in the O₂ MO diagram, the highest occupied molecular orbitals are a pair of degenerate π_2p antibonding orbitals holding two electrons; by Hund's rule of maximum multiplicity, these two electrons occupy separate π orbitals with parallel spins (S = 1 triplet state), giving O₂ a permanent magnetic dipole moment while maintaining a net bond order of (8 − 4)/2 = 2.

What is the difference between Valence Bond (VB) theory and Molecular Orbital (MO) theory at high levels of quantum chemistry?

At their simplest single-determinant/single-configuration textbook levels, VB theory (Heitler–London) overemphasizes electron correlation (excluding ionic terms A⁺B⁻ so H₂ dissociates correctly to neutral H + H) whereas MO theory (Hartree–Fock LCAO) overemphasizes delocalization (including 50% ionic character so restricted HF H₂ dissociates incorrectly at R → ∞). However, when VB theory is extended by adding ionic resonance structures and when MO theory is extended by adding Configuration Interaction (CI) (mixing the doubly excited (σ^*)² determinant into the ground (σ)² state), both theories converge to the exact same multi-configurational wave function and identical physical predictions.

1. Lewis, G. N. (1916). "The atom and the molecule." Journal of the American Chemical Society, 38(4), 762–785. https://doi.org/10.1021/ja02261a002
2. Heitler, W., & London, F. (1927). "Wechselwirkung neutraler Atome und homöopolare Bindung nach der Quantenmechanik." Zeitschrift für Physik, 44(6–7), 455–472. https://doi.org/10.1007/BF01397394
3. Pauling, L. (1931). "The nature of the chemical bond. Application of results obtained from the quantum mechanics and from a theory of paramagnetic susceptibility to the structure of molecules." Journal of the American Chemical Society, 53(4), 1367–1400. https://doi.org/10.1021/ja01355a027
4. Ruedenberg, K. (1962). "The physical nature of the chemical bond." Reviews of Modern Physics, 34(2), 326–376. https://doi.org/10.1103/RevModPhys.34.326