Study molecular orbitals, their formation from atomic orbitals, and their importance in explaining chemical bonding and molecular structure.
Molecular orbitals (MOs) embody the quantum mechanical description of electrons within molecules, extending the concept of atomic orbitals beyond isolated atoms to systems composed of multiple nuclei. Introduced by Robert S. Mulliken in 1932 as one-electron orbital wave functions, molecular orbitals represent mathematical functions that describe electron locations and their wave-like behavior across an entire molecule [1]. Unlike atomic orbitals confined to single atoms, molecular orbitals result from a constructive or destructive combination of atomic orbitals when atoms chemically bond.
At the core of molecular orbital theory lies the Schrödinger equation, which governs the quantum state of electrons. Exact solutions for multi-electron molecules remain computationally intractable; thus, molecular orbitals serve as approximate solutions to this equation in the electric field produced by all nuclei and averaged electron interactions within a molecule. These approximations often utilize methods like Hartree–Fock or self-consistent field (SCF), which treat electrons independently but account for average repulsions between them [1].
The fundamental approach to constructing molecular orbitals is through the linear combination of atomic orbitals (LCAO). This method expresses each molecular orbital as a weighted sum or difference of atomic wavefunctions from constituent atoms:
\[
\Psi = c_{a}\psi_{a} + c_{b}\psi_{b}
\]
for bonding combinations and
\[
\Psi^{*} = c_{a}\psi_{a} - c_{b}\psi_{b}
\]
for antibonding combinations, where \(c_a\), \(c_b\) are coefficients representing contributions from each atomic orbital \(\psi_a\), \(\psi_b\), respectively [1]. This formalism captures how electron probability densities overlap and interfere across different atoms.
The LCAO approximation was pioneered by Sir John Lennard-Jones in 1929, building upon the introduction of molecular orbital theory by Friedrich Hund and Robert S. Mulliken in 1927 and 1928. Lennard-Jones's work demonstrated how to derive the electronic structure of fluorine and oxygen molecules, laying the groundwork for modern quantum chemistry models [1].
Molecular orbitals fall into three categories based on their energy relative to the original atomic orbitals:
- Bonding orbitals possess lower energy than their parent atomic orbitals due to constructive interference increasing electron density between nuclei. This stabilizes the molecule by promoting chemical bonds.
- Antibonding orbitals have higher energy because destructive interference creates nodes—regions with zero electron density—between nuclei, destabilizing bonds.
- Nonbonding orbitals retain energy levels similar to their constituent atomic orbitals and do not contribute significantly to bonding interactions [1].
Electrons fill these molecular orbitals according to principles analogous to atomic electronic configurations but considering spin pairing per the Pauli exclusion principle. Occupation of bonding MOs strengthens bonds; occupation of antibonding MOs weakens or prevents bonding.
Molecular symmetry plays a crucial role in determining allowed combinations of atomic orbitals into molecular orbitals. Only atomic orbitals with compatible symmetries can interact effectively. If a molecule possesses symmetry elements such as reflections or rotations, its nondegenerate molecular orbitals will be symmetric (unchanged sign upon symmetry operation \(S\), \(S\psi=+\psi\)) or antisymmetric (sign reversal upon \(S\), \(S\psi=-\psi\)) with respect to those operations [1].
Planar molecules illustrate this clearly: molecular orbitals are classified as sigma (σ), symmetric under reflection through the plane, or pi (π), antisymmetric under that reflection. Degenerate molecular orbital energies correspond to bases for irreducible representations of the molecule’s symmetry group—a formalism rooted in group theory that predicts possible orbital interactions and degeneracies.
Canonical molecular orbitals provide delocalized descriptions spanning entire molecules; they are eigenfunctions of the Fock operator in Hartree–Fock theory and possess definite energies linked directly to physical observables. Alternatively, localized molecular orbitals emerge from unitary transformations on canonical sets designed to resemble classical Lewis structures more closely by concentrating electron density into bonds or lone pairs. While these localized forms aid chemical intuition by mimicking localized bonding patterns, their energies lose direct physical meaning due to this transformation’s mathematical nature [1][3].
Molecular orbital diagrams graphically represent relative energies and occupancy of MOs formed from combining atomic orbitals during bond formation. These diagrams track how electrons distribute among bonding and antibonding levels as atoms combine.
For diatomic molecules involving elements with atomic numbers less than or equal to 7 (e.g., B₂, C₂, N₂), the ordering of valence molecular orbital energies follows:
\[
\sigma (1s) < \sigma^{*}(1s) < \sigma(2s) < \sigma^{*}(2s) < \pi(2p_y)= \pi(2p_z)< \sigma(2p_x)< \pi^{*}(2p_y)= \pi^{*}(2p_z)< \sigma^{*}(2p_x)
\]
Here s–p mixing causes π(2p) MOs to lie lower in energy than σ(2p). In contrast, for elements with atomic numbers greater than or equal to 8 (e.g., O₂, F₂, Ne₂), the ordering changes due to negligible s–p mixing:
\[
\sigma (1s)< \sigma^{*}(1s)< \sigma(2s)< \sigma^{*}(2s)< \sigma(2p_x)< \pi(2p_y)= \pi(2p_z)< \pi^{*}(2p_y)= \pi^{*}(2p_z)< \sigma^{*}(2p_x)
\]
In this regime, σ(2p_x) lies below π(2p_y/z). These variations impact predicted bond orders and magnetic properties derived from MO occupancy patterns [4].
Bond order provides an approximate measure correlating with bond strength and length. It is computed as half the difference between electrons occupying bonding versus antibonding MOs:
\[
\text{Bond Order} = \frac{N_b - N_a}{2}
\]
where \(N_b\) is number of electrons in bonding MOs and \(N_a\) those in antibonding MOs. Higher bond orders imply stronger bonds with shorter lengths.
Examples include:
- H₂ has a bond order indicating stable covalent bonding.
- He₂ does not form a stable molecule due to equal filling of bonding and antibonding MOs resulting in zero net bond order.
- O₂ exhibits paramagnetism because unpaired electrons occupy degenerate π* antibonding MOs.
Magnetic properties also correlate directly with MO occupancy: fully paired electrons yield diamagnetism; presence of unpaired electrons yields paramagnetism detectable experimentally [4].
Modern computational chemistry relies heavily on calculating canonical molecular orbitals as foundational building blocks for electronic structure methods. These calculations underpin density functional theory (DFT), post-Hartree–Fock correlation treatments, and excited-state simulations.
However, limitations arise since canonical MOs represent single-particle approximations ignoring explicit electron correlation beyond mean-field effects. More accurate approaches incorporate configuration interaction where wavefunctions extend beyond simple MO products but lose straightforward orbital interpretations altogether [1].
---
Molecular orbital theory remains indispensable for understanding chemical bonding at a fundamental level through quantum mechanics. Its rigorous framework connects microscopic electron behavior expressed by wavefunctions with macroscopic chemical phenomena including bond formation energies, magnetic properties, reactivity trends, and spectroscopic signatures.
Generating summary…