2.4 Molecular Orbital Theory and Magnetochemistry of Diatomic Molecules

Key Takeaways

  • Molecular orbital theory constructs delocalized wavefunctions via Linear Combination of Atomic Orbitals (LCAO-MO); effective overlap requires comparable orbital energies, compatible spatial symmetry, and spatial proximity.

  • Parity inversion symmetry classifies MOs in centrosymmetric diatomics into gerade (gg, symmetric) and ungerade (uu, antisymmetric); bonding σ\sigma orbitals are σg\sigma_g, antibonding σ∗\sigma^* are σu\sigma_u, bonding π\pi orbitals are πu\pi_u, and antibonding π∗\pi^* are πg\pi_g.

  • Orbital crossover / s−ps-p mixing occurs in Li2,Be2,B2,C2,N2Li_2, Be_2, B_2, C_2, N_2 because their small 2s−2p2s-2p energy gap enables mixing between σ2s\sigma_{2s} and σ2pz\sigma_{2p_z}, pushing σ2pz\sigma_{2p_z} above the degenerate π2p\pi_{2p} orbitals; in O2O_2 and F2F_2, the large energy separation prevents significant mixing, maintaining the standard order.

  • Triplet dioxygen (O2O_2) has bond order 2.0 with two unpaired electrons in degenerate π2p∗\pi_{2p}^* orbitals, directly explaining its experimentally observed paramagnetism (μso=2.83 BM\mu_{so} = 2.83\text{ BM}), a property inexplicable under classical Lewis structures.

  • In heteronuclear diatomics such as COCO, the HOMO (3σ3\sigma) and LUMO (2π∗2\pi^*) are heavily polarized toward carbon, explaining why COCO acts as a carbon-centered σ\sigma-donor and π\pi-acceptor in transition metal complexes.

Last updated: October 2026

Molecular Orbital Theory and Magnetochemistry of Diatomic Molecules

Molecular Orbital (MO) Theory provides an advanced quantum mechanical description of chemical bonding where valence electrons delocalize over the entire molecular framework rather than remaining localized between paired atoms. This section covers the Linear Combination of Atomic Orbitals (LCAO) approximation, orbital symmetry and parity, s−ps-p mixing crossover, paramagnetism, and Frontier Molecular Orbital (FMO) theory tested on the PRC Chemist Licensure Examination.


1. Principles of Linear Combination of Atomic Orbitals (LCAO-MO)

In the LCAO approximation, a molecular orbital ψMO\psi_{MO} is formulated as a linear combination of atomic basis functions ϕi\phi_i:

ψMO=cAϕA+cBϕB\psi_{MO} = c_A \phi_A + c_B \phi_B

where cAc_A and cBc_B are weighting coefficients representing the orbital contribution of each atom. The total electron probability density is given by squaring the wavefunction:

∣ψMO∣2=cA2ϕA2+cB2ϕB2+2cAcBϕAϕB\lvert\psi_{MO}\rvert^2 = c_A^2 \phi_A^2 + c_B^2 \phi_B^2 + 2 c_A c_B \phi_A \phi_B

  • The terms cA2ϕA2c_A^2 \phi_A^2 and cB2ϕB2c_B^2 \phi_B^2 represent electron densities associated with isolated atomic orbitals.
  • The cross-term 2cAcBϕAϕB2 c_A c_B \phi_A \phi_B represents quantum mechanical interference:
    • When wavefunctions combine in-phase (++), constructive interference increases electron probability density in the internuclear region, shielding the nuclei and lowering the potential energy   ⟹  \implies Bonding Molecular Orbital (lower in energy than parent atomic orbitals).
    • When wavefunctions combine out-of-phase (−-), destructive interference produces a nodal plane between the nuclei where ψ=0\psi = 0, exposing the nuclei to mutual electrostatic repulsion   ⟹  \implies Antibonding Molecular Orbital (designated with ∗*, higher in energy than parent atomic orbitals).

The Three Cardinal Rules of Effective Orbital Overlap

  1. Energy Compatibility: Atomic orbitals must have comparable energies (typically ΔE≤12−14 eV\Delta E \le 12-14\text{ eV}). If two atomic orbitals differ greatly in energy, they mix weakly and remain effectively nonbonding.
  2. Symmetry Compatibility: Atomic orbitals must possess identical spatial symmetry with respect to the internuclear axis (zz). The overlap integral SABS_{AB} must be non-zero: SAB=∫ϕAϕB dτ≠0S_{AB} = \int \phi_A \phi_B \, d\tau \neq 0 For example, an ss orbital has zero net overlap with a pxp_x or pyp_y orbital oriented perpendicular to the bond axis because positive and negative phase overlaps cancel identically (S=0S = 0).
  3. Spatial Proximity: Orbitals must be sufficiently close in space for their radial wavefunctions to overlap significantly.

Symmetry Classification and Parity (gg vs uu)

  • σ\sigma (Sigma) Orbitals: Cylindrically symmetric about the internuclear axis. Rotation by any arbitrary angle around the bond axis leaves the wavefunction unchanged. Formed by coaxial overlap (s+s,s+pz,pz+pzs+s, s+p_z, p_z+p_z).
  • π\pi (Pi) Orbitals: Possess a nodal plane containing the internuclear axis. A 180∘180^\circ rotation around the bond axis reverses the sign of the wavefunction. Formed by parallel sideways overlap (px+px,py+pyp_x+p_x, p_y+p_y).
  • Parity / Inversion Symmetry: In centrosymmetric (homonuclear) diatomics, orbitals are classified under the inversion operation i^\hat{i} through the molecular midpoint:
    • Gerade (gg): Symmetric upon inversion: ψ(−x,−y,−z)=+ψ(x,y,z)\psi(-x, -y, -z) = +\psi(x, y, z).
      • σ\sigma bonding orbitals are σg\sigma_g.
      • π∗\pi^* antibonding orbitals are πg∗\pi_g^*.
    • Ungerade (uu): Antisymmetric upon inversion: ψ(−x,−y,−z)=−ψ(x,y,z)\psi(-x, -y, -z) = -\psi(x, y, z).
      • σ∗\sigma^* antibonding orbitals are σu∗\sigma_u^*.
      • π\pi bonding orbitals are πu\pi_u.

2. Homonuclear Diatomics and the s−ps-p Mixing Crossover

Bond Order Formulation

The bond order quantifies the net number of covalent bonds: Bond Order=12(Nb−Na)\text{Bond Order} = \frac{1}{2}(N_b - N_a) where NbN_b is the number of electrons in bonding molecular orbitals and NaN_a is the number of electrons in antibonding molecular orbitals. A positive bond order indicates thermodynamic stability relative to separated atoms; higher bond orders correlate with shorter internuclear bond distances and higher bond dissociation energies.

The Mechanism of s−ps-p Mixing (Orbital Crossover)

In second-period homonuclear diatomics, atomic 2s2s and 2pz2p_z orbitals both form molecular orbitals of Σg+\Sigma_g^+ symmetry ( σ2s and σ2pz \,\sigma_{2s}\,\text{and}\,\sigma_{2p_z}\,, respectively). Quantum mechanics dictates that orbitals sharing the same symmetry interact and repel one another if their energy gap is small:

  1. Without s−ps-p Mixing (O2,F2,Ne2O_2, F_2, Ne_2): In Oxygen and Fluorine, the high effective nuclear charge binds 2s2s electrons very tightly, creating a large energy separation between 2s2s and 2p2p (ΔE2s−2p=16.5 eV\Delta E_{2s-2p} = 16.5\text{ eV} in O, 21.6 eV21.6\text{ eV} in F). Minimal mixing occurs. The standard energy sequence is: σ2s(1σg)<σ2s∗(1σu∗)<σ2pz(2σg)<π2px=π2py(1πu)<π2px∗=π2py∗(1πg∗)<σ2pz∗(2σu∗)\sigma_{2s} (1\sigma_g) < \sigma_{2s}^* (1\sigma_u^*) < \mathbf{\sigma_{2p_z} (2\sigma_g)} < \mathbf{\pi_{2p_x} = \pi_{2p_y} (1\pi_u)} < \pi_{2p_x}^* = \pi_{2p_y}^* (1\pi_g^*) < \sigma_{2p_z}^* (2\sigma_u^*)

  2. With s−ps-p Mixing (Li2,Be2,B2,C2,N2Li_2, Be_2, B_2, C_2, N_2): In lighter elements, lower nuclear charge yields a small 2s−2p2s-2p energy gap (ΔE2s−2p=5.7 eV\Delta E_{2s-2p} = 5.7\text{ eV} in B, 8.8 eV8.8\text{ eV} in C, 12.4 eV12.4\text{ eV} in N). Strong quantum mixing occurs between the 1σg1\sigma_g and 2σg2\sigma_g orbitals: the lower 1σg1\sigma_g is stabilized further down, while the upper 2σg2\sigma_g (σ2pz\sigma_{2p_z}) is pushed upward above the degenerate 1πu1\pi_u orbitals! The modified energy sequence for B2,C2,N2B_2, C_2, N_2 is: σ2s(1σg)<σ2s∗(1σu∗)<π2px=π2py(1πu)<σ2pz(2σg)<π2px∗=π2py∗(1πg∗)<σ2pz∗(2σu∗)\sigma_{2s} (1\sigma_g) < \sigma_{2s}^* (1\sigma_u^*) < \mathbf{\pi_{2p_x} = \pi_{2p_y} (1\pi_u)} < \mathbf{\sigma_{2p_z} (2\sigma_g)} < \pi_{2p_x}^* = \pi_{2p_y}^* (1\pi_g^*) < \sigma_{2p_z}^* (2\sigma_u^*)

Systematic Period 2 Homonuclear Diatomic Series

MoleculeValence e−e^-Ground-State ConfigurationBond OrderMagnetic CharacterExperimental Properties
B2B_26σ2s2σ2s∗2π2px1π2py1\sigma_{2s}^2 \sigma_{2s}^{*2} \pi_{2p_x}^1 \pi_{2p_y}^14−22=1\frac{4-2}{2} = 1Paramagnetic (2 e−2\, e^-)D0=290 kJ/molD_0 = 290\text{ kJ/mol}; provides direct proof of s−ps-p mixing (without mixing, it would be σ2pz2\sigma_{2p_z}^2, diamagnetic).
C2C_28σ2s2σ2s∗2π2p4\sigma_{2s}^2 \sigma_{2s}^{*2} \pi_{2p}^46−22=2\frac{6-2}{2} = 2Diamagnetic (0 e−0\, e^-)D0=599 kJ/molD_0 = 599\text{ kJ/mol}; unusual pure double bond composed of two π\pi bonds with zero net σ\sigma bonding!
N2N_210σ2s2σ2s∗2π2p4σ2pz2\sigma_{2s}^2 \sigma_{2s}^{*2} \pi_{2p}^4 \sigma_{2p_z}^28−22=3\frac{8-2}{2} = 3Diamagnetic (0 e−0\, e^-)D0=945 kJ/mol,r=1.098 A˚D_0 = 945\text{ kJ/mol}, r = 1.098\text{ \AA}; exceptionally inert triple bond.
O2O_212σ2s2σ2s∗2σ2pz2π2p4π2px∗1π2py∗1\sigma_{2s}^2 \sigma_{2s}^{*2} \sigma_{2p_z}^2 \pi_{2p}^4 \pi_{2p_x}^{*1} \pi_{2p_y}^{*1}8−42=2\frac{8-4}{2} = 2Paramagnetic (2 e−2\, e^-)D0=498 kJ/mol,r=1.207 A˚D_0 = 498\text{ kJ/mol}, r = 1.207\text{ \AA}; explains liquid O2O_2 attraction into magnetic poles.
F2F_214σ2s2σ2s∗2σ2pz2π2p4π2p∗4\sigma_{2s}^2 \sigma_{2s}^{*2} \sigma_{2p_z}^2 \pi_{2p}^4 \pi_{2p}^{*4}8−62=1\frac{8-6}{2} = 1Diamagnetic (0 e−0\, e^-)D0=155 kJ/mol,r=1.412 A˚D_0 = 155\text{ kJ/mol}, r = 1.412\text{ \AA}; weak bond due to 33 pairs of repulsive nonbonding lone pairs.
Ne2Ne_216σ2s2σ2s∗2σ2pz2π2p4π2p∗4σ2pz∗2\sigma_{2s}^2 \sigma_{2s}^{*2} \sigma_{2p_z}^2 \pi_{2p}^4 \pi_{2p}^{*4} \sigma_{2p_z}^{*2}8−82=0\frac{8-8}{2} = 0—Unbound van der Waals dimer only.

3. Magnetochemistry and the Triplet Ground State of O2O_2

Paramagnetism vs Diamagnetism

  • Diamagnetism: All electrons in the molecule are paired. When placed in an external magnetic field gradient, paired electrons generate an induced orbital magnetic moment that opposes the applied field, causing weak expulsion.
  • Paramagnetism: The molecule contains one or more unpaired electrons. The intrinsic electron spin produces a permanent magnetic dipole moment that aligns parallel to the external field, drawing the substance strongly into the magnetic field.

Spin-Only Magnetic Moment Formula

The quantitative magnetic moment due to electron spin is given by:

μso=n(n+2) μB\mu_{so} = \sqrt{n(n+2)} \, \mu_B

where nn is the number of unpaired electrons and μB\mu_B is the Bohr Magneton (μB=eℏ2me=9.274×10−24 J/T\mu_B = \frac{e\hbar}{2m_e} = 9.274 \times 10^{-24}\text{ J/T}). Measurements using a Gouy balance or SQUID magnetometer determine μso\mu_{so} directly.

The Triumphant Explanation of O2O_2 Paramagnetism

Under classical Lewis bonding, dioxygen is written as :O¨=O¨::\ddot{O}=\ddot{O}:, predicting an octet-complete, diamagnetic molecule with zero unpaired electrons. However, liquid oxygen is strongly paramagnetic and clings suspended between magnetic pole pieces.

MO theory resolves this discrepancy effortlessly: Dioxygen possesses 12 valence electrons. Filling the molecular orbitals in order of increasing energy: σ2s2  σ2s∗2  σ2pz2  π2p4  (π2px∗1=π2py∗1)\sigma_{2s}^2 \; \sigma_{2s}^{*2} \; \sigma_{2p_z}^2 \; \pi_{2p}^4 \; (\pi_{2p_x}^{*1} = \pi_{2p_y}^{*1}) Ten electrons fill the bonding and core-like antibonding orbitals. The remaining two electrons must enter the degenerate π2p∗\pi_{2p}^* antibonding orbitals. By Hund's rule, they occupy separate degenerate orbitals with parallel spins (S=12+12=1S = \frac{1}{2} + \frac{1}{2} = 1, spin multiplicity 2S+1=32S+1 = 3, forming a triplet ground state (3Σg−^3\Sigma_g^-)): μso(O2)=2(2+2)=8≈2.83 BM\mu_{so}(O_2) = \sqrt{2(2+2)} = \sqrt{8} \approx 2.83\text{ BM}

Dioxygen Redox Series: Bond Order, Length, and Magnetism

Stepwise reduction or oxidation of O2O_2 systematically alters the population of the antibonding π∗\pi^* orbitals:

SpeciesValence e−e^-π∗\pi^* PopulationBond OrderUnpaired e−e^-μso\mu_{so} (BM)Bond Length (A˚\text{\AA})Stretching Frequency (cm−1\text{cm}^{-1})
Dioxygenyl (O2+O_2^+)1112.52.511.731.731.121.1218581858
Dioxygen (O2O_2)1222.02.022.832.831.211.2115561556
Superoxide (O2−O_2^-)1331.51.511.731.731.281.2811081108
Peroxide (O22−O_2^{2-})1441.01.000.000.00 (dia)1.491.49770770

Thermodynamic Deduction: Removing electrons from antibonding π∗\pi^* orbitals (O2→O2+O_2 \rightarrow O_2^+) increases bond order, shortens bond length, and raises the vibrational stretching frequency. Adding electrons into π∗\pi^* (O2→O2−→O22−O_2 \rightarrow O_2^- \rightarrow O_2^{2-}) weakens and lengthens the O−OO-O bond.


4. Heteronuclear Diatomics and Frontier Molecular Orbital (FMO) Theory

In heteronuclear diatomics (A−BA-B), unequal electronegativities cause the atomic orbitals of the more electronegative atom to lie lower in energy. Consequently:

  • Bonding MOs are lower in energy and receive greater mathematical contribution (cB>cAc_B > c_A) from the more electronegative atom.
  • Antibonding MOs are higher in energy and receive greater mathematical contribution (cA∗>cB∗c_A^* > c_B^*) from the less electronegative atom.

Carbon Monoxide (COCO)

Carbon Monoxide has 10 valence electrons (44 from C, 66 from O) and is isoelectronic with N2N_2:

  • Atomic orbital energies: C(2s)=−19.4 eV,C(2p)=−10.7 eVC(2s) = -19.4\text{ eV}, C(2p) = -10.7\text{ eV}; O(2s)=−32.4 eV,O(2p)=−15.9 eVO(2s) = -32.4\text{ eV}, O(2p) = -15.9\text{ eV}.
  • Valence configuration: (1σ)2(2σ)2(1π)4(3σ)2(1\sigma)^2 (2\sigma)^2 (1\pi)^4 (3\sigma)^2
  • HOMO (3σ3\sigma): The Highest Occupied Molecular Orbital is a weakly antibonding / nonbonding σ\sigma orbital. Crucially, the 3σ3\sigma HOMO is polarized heavily toward Carbon (cC≫cOc_C \gg c_O), with a large lone-pair lobe protruding outward from the carbon nucleus into space!
  • LUMO (2π∗2\pi^*): The Lowest Unoccupied Molecular Orbital consists of empty degenerate π∗\pi^* orbitals that are also polarized predominantly on Carbon (cC≫cOc_C \gg c_O).

The Biochemical & Catalytic Consequence: Although Oxygen is far more electronegative, COCO binds transition metals exclusively via Carbon! Carbon monoxide acts as a classical σ\sigma-donor (donating its carbon-centered HOMO 3σ3\sigma electron density into empty metal dd-orbitals) and a powerful π\pi-acceptor (accepting electron density from filled metal dd-orbitals into its empty carbon-centered LUMO 2π∗2\pi^*). This synergistic bonding mechanism (Dewar-Chatt-Duncanson model) explains why COCO binds hemoglobin 200×200\times more tightly than O2O_2, causing lethal carbon monoxide poisoning.

Nitric Oxide (NONO)

Nitric oxide has 11 valence electrons; ground-state configuration: (1σ)2(2σ)2(1π)4(3σ)2(2π∗)1(1\sigma)^2 (2\sigma)^2 (1\pi)^4 (3\sigma)^2 (2\pi^*)^1

  • Bond order = 8−32=2.5\frac{8 - 3}{2} = 2.5.
  • The unpaired electron resides in an antibonding 2π∗2\pi^* orbital polarized primarily on nitrogen.
  • Oxidation to the nitrosonium ion (NO+NO^+) removes this antibonding electron, increasing the bond order to 3.03.0, contracting the bond from 1.15 A˚1.15\text{ \AA} to 1.06 A˚1.06\text{ \AA}, and shifting the vibrational frequency from 1876 cm−11876\text{ cm}^{-1} to 2300 cm−12300\text{ cm}^{-1}.

Hydrogen Fluoride (HFHF)

  • H(1s)H(1s) lies at −13.6 eV-13.6\text{ eV}. F(2s)F(2s) lies at −40.2 eV-40.2\text{ eV} (too low to interact; forms nonbonding 1σ1\sigma).
  • F(2pz)F(2p_z) lies at −18.6 eV-18.6\text{ eV} and shares σ\sigma symmetry with H(1s)H(1s), combining to form a bonding 2σ2\sigma and antibonding 3σ∗3\sigma^* orbital.
  • F(2px)F(2p_x) and F(2py)F(2p_y) have π\pi symmetry and zero overlap with H(1s)H(1s) (S=0S=0), remaining completely nonbonding lone pairs (1π1\pi).
  • Bond order = 1.01.0; the MO description directly mirrors the three nonbonding lone pairs on fluorine in the Lewis structure.
Loading diagram...
Molecular Orbital Energy Level Ordering: Normal vs. s-p Mixed Diatomics
Test Your Knowledge

Which of the following correctly pairs the diatomic oxygen species with their respective bond order and bond length ranking?

A

O2(2-) < O2(-) < O2 < O2(+) in bond order, and O2(+) < O2 < O2(-) < O2(2-) in bond length

B

O2(+) < O2 < O2(-) < O2(2-) in bond order, and O2(+) < O2 < O2(-) < O2(2-) in bond length

C

O2(2-) < O2(-) < O2 < O2(+) in bond order, and O2(2-) < O2(-) < O2 < O2(+) in bond length

D

O2(-) < O2(2-) < O2 < O2(+) in bond order, and O2 < O2(+) < O2(-) < O2(2-) in bond length

Test Your Knowledge

Why does s-p orbital mixing invert the energy ordering of the sigma(2pz) and pi(2p) molecular orbitals in B2, C2, and N2, but not in O2 and F2?

A

In B2, C2, and N2 the small 2s-2p gap lets the sigma_g orbitals mix, pushing sigma(2pz) above pi(2p).

B

s-p mixing occurs in B2, C2, and N2 because these molecules contain unpaired electrons in their ground states, whereas O2 and F2 are diamagnetic.

C

s-p mixing is caused exclusively by spin-orbit coupling that splits the degenerate pi(2p) levels in light diatomic molecules.

D

s-p mixing occurs in O2 and F2 because their higher effective nuclear charge brings the 2s and 2p atomic energy levels closer together.

Test Your Knowledge

Frontier Molecular Orbital (FMO) analysis reveals that Carbon Monoxide (CO) coordinates to transition metal centers through the carbon atom rather than the oxygen atom. What is the fundamental quantum mechanical basis for this behavior?

A

Carbon monoxide binds metals through oxygen because oxygen is more electronegative and possesses a greater partial negative charge.

B

Carbon monoxide coordinates through carbon because both the HOMO (3sigma) and LUMO (2pi*) have their largest orbital coefficients localized on the carbon atom.

C

Carbon monoxide binds metals through oxygen because the 1pi bonding molecular orbitals are concentrated entirely on carbon.

D

Carbon monoxide binds metals through carbon solely because the carbon atom carries an expanded octet of 10 valence electrons.

Sections you finish are checked off in the contents.