Unit 4: Hybridization, Shape, and Hückel Molecular Orbital Theory

1. Quantum Mechanical Approach to Hybridization

Hybridization is a quantum mechanical concept where atomic orbitals on a single atom mix to form new, equivalent hybrid orbitals. This model explains molecular geometries and bond angles that cannot be understood through pure atomic orbitals (s, p, d, etc.) alone.

Mathematical Framework and Bond Angles

The general form of a hybrid orbital constructed from s and p orbitals is represented as:

Ψ = c_sψ_s + c_pψ_p

For any spⁿ hybrid system, the mixing coefficients must satisfy normalization and orthogonality conditions. The angle θ between any two equivalent hybrid orbitals is mathematically related to the hybridization index n by the formula:

cos θ = -1 / n

sp³ Hybridization

In sp³ hybridization, one 2s orbital and three 2p orbitals (2p_x, 2p_y, and 2p_z) mix to form four equivalent hybrid orbitals directed toward the corners of a regular tetrahedron.

The four normalized wavefunctions (Ψ₁ to Ψ₄) are given by:

Ψ₁ = (1/2) (s + p_x + p_y + p_z)
Ψ₂ = (1/2) (s - p_x - p_y + p_z)
Ψ₃ = (1/2) (s - p_x + p_y - p_z)
Ψ₄ = (1/2) (s + p_x - p_y - p_z)

Here, the hybridization index n = 3. Substituting this into the angle formula:

cos θ = -1 / 3 ≈ -0.3333
θ = 109.5°

This explains the perfect tetrahedral geometry of molecules like methane (CH₄).

sp² Hybridization

In sp² hybridization, one 2s orbital mixes with two 2p orbitals (typically 2p_x and 2p_y), leaving the 2p_z orbital unhybridized. This forms three equivalent hybrid orbitals in a single plane.

The wavefunctions are mathematically represented as:

Ψ₁ = (1/√3) s + (√2/√3) p_x
Ψ₂ = (1/√3) s - (1/6) p_x + (1/√2) p_y
Ψ₃ = (1/√3) s - (1/6) p_x - (1/√2) p_y

Here, the hybridization index n = 2. Substituting into the angle formula:

cos θ = -1 / 2 = -0.5
θ = 120°

This results in a planar trigonal arrangement, characteristic of molecules like ethylene (C₂H₄).

sp Hybridization

In sp hybridization, one 2s orbital mixes with one 2p orbital (typically 2p_x), leaving the 2p_y and 2p_z orbitals unhybridized to form π bonds. This forms two equivalent hybrid orbitals oriented in opposite directions.

The wavefunctions are:

Ψ₁ = (1/√2) (s + p_x)
Ψ₂ = (1/2) (s - p_x)

Here, the hybridization index n = 1. Substituting into the angle formula:

cos θ = -1 / 1 = -1.0
θ = 180°

This results in a linear geometry, characteristic of molecules like acetylene (C₂H₂).

2. The π-Electron Approximation and Hückel MO Theory

The π-Electron Approximation

For planar conjugated molecules, the total wavefunction can be separated into two independent parts: σ-electrons and π-electrons. This separation is known as the π-electron approximation. The σ-bonds form a rigid, symmetric framework that lies in the molecular plane, while the π-system is formed by the lateral overlap of 2p_z atomic orbitals perpendicular to this plane. Under this approximation, the properties of the molecule's chemical reactivity and UV-visible spectra are assumed to be dominated entirely by the π-electrons moving in the average potential of the core σ-skeleton.

The Hückel Molecular Orbital (HMO) Approximation

Erich Hückel simplified the quantum mechanical treatment of π-electrons by introducing several key approximations to solve the Schrödinger equation for conjugated systems:

• LCAO Method: Each molecular orbital Ψ_i is expressed as a Linear Combination of Atomic Orbitals (φ_r) of the carbon 2p_z orbitals:
  Ψ_i = ∑ c_ir φ_r
• Coulomb Integral (α): The energy of an electron in an isolated carbon 2p_z orbital is constant for all carbon atoms:
  H_rr = ∫ φ_r H φ_r dτ = α
• Resonance (Overlap) Integral (β): The interaction energy between two adjacent atomic orbitals. It is assumed to be constant for bonded atoms and zero for non-bonded atoms:
  H_rs = β (if r and s are directly bonded)
  H_rs = 0 (if r and s are not bonded)
• Overlap Integral (S_rs): The overlap between different atomic orbitals is neglected to simplify calculations:
  S_rs = 1 (if r = s)
  S_rs = 0 (if r ≠ s)

The secular equations are derived from the variation principle, resulting in the secular determinant:

| H_rs - E S_rs | = 0

3. Simple Hückel Treatment of Specific Systems

Ethylene System (and Note on Ethane)

Note: In physical chemistry syllabi, the term "ethane" in the context of Hückel theory is a common historical typo for the ethylene (ethene) system. Because ethane is saturated and has no π-electrons, the Hückel treatment applies to the unsaturated counterpart, ethylene, representing an isolated two-center carbon-carbon π-bond.

For ethylene (C₁=C₂), the secular determinant is:

| α - E     β |
|    β     α - E | = 0

Dividing by β and setting x = (α - E) / β, we get:

| x   1 |
| 1   x | = 0

Solving this determinant gives:

x² - 1 = 0   →   x = ±1

Since β is a negative quantity, the energy levels are:

E₁ = α + β (Bonding MO, lower energy)
E₂ = α - β (Antibonding MO, higher energy)

The total π-electron energy of ethylene (holding 2 electrons in E₁) is:

E_π = 2(α + β) = 2α + 2β

Allyl System

The allyl system (C₁-C₂-C₃) consists of three carbon atoms. The secular determinant is:

| x   1   0 |
| 1   x   1 | = 0
| 0   1   x |

Expanding the determinant gives the polynomial:

x³ - 2x = 0   →   x(x² - 2) = 0

The roots are x = -√2, 0, and +√2. This yields three energy levels:

E₁ = α + √2 β (Bonding MO)
E₂ = α (Non-bonding MO)
E₃ = α - √2 β (Antibonding MO)

Depending on the species, these energy levels are occupied by different numbers of π-electrons:

Butadiene System (1,3-Butadiene)

The 1,3-butadiene system is a conjugated chain of four carbon atoms (C₁-C₂-C₃-C₄). The secular determinant is:

| x   1   0   0 |
| 1   x   1   0 | = 0
| 0   1   x   1 |
| 0   0   1   x |

Solving this determinant yields the characteristic equation:

x⁴ - 3x² + 1 = 0

The roots of this quadratic-in-x² equation are:

x = ±1.618   and   x = ±0.618

The resulting energy levels in order of increasing energy are:

E₁ = α + 1.618 β
E₂ = α + 0.618 β
E₃ = α - 0.618 β
E₄ = α - 1.618 β

The four π-electrons of butadiene occupy the two lowest molecular orbitals (E₁ and E₂):

E_π = 2(α + 1.618 β) + 2(α + 0.618 β) = 4α + 4.472 β

We can calculate the Delocalization Energy (DE) by comparing this to two isolated ethylene molecules:

E_isolated = 2 × (2α + 2β) = 4α + 4β
DE = E_π - E_isolated = (4α + 4.472 β) - (4α + 4β) = 0.472 β

Since β is negative, this difference represents a significant stabilization of the molecule due to conjugation.

4. Hückel's Rule of Aromaticity and Cyclic Systems

Hückel's Rule of Aromaticity

Hückel's rule provides a quantum mechanical criterion for predicting whether a planar, cyclic, conjugated monocyclic system will exhibit aromatic stabilization:

A planar, cyclic, completely conjugated system containing (4n + 2) π-electrons (where n = 0, 1, 2, 3...) exhibits exceptional stability and is classified as aromatic. Systems with 4n π-electrons exhibit instability and are classified as antiaromatic.

Delocalization Energy of Cyclic Systems

Let us compare the delocalization energy of a 4n cyclic system (Cyclobutadiene) with a (4n + 2) cyclic system (Benzene) to observe this effect quantitatively.

Cyclobutadiene (4n system, n = 1)

Cyclobutadiene has four carbon atoms in a closed ring. Its secular determinant is:

| x   1   0   1 |
| 1   x   1   0 | = 0
| 0   1   x   1 |
| 1   0   1   x |

Solving this determinant yields the roots:

x = -2, 0, 0, +2

The energy levels are:

E₁ = α + 2β
E₂, E₃ = α (Degenerate non-bonding levels)
E₄ = α - 2β

For 4 π-electrons, the total energy is:

E_π = 2(α + 2β) + 2(α) = 4α + 4β

Comparing this to two isolated ethylene molecules (4α + 4β):

DE = E_π - E_isolated = (4α + 4β) - (4α + 4β) = 0

Because the delocalization energy is zero, cyclobutadiene does not enjoy thermodynamic stabilization and is highly unstable (antiaromatic).

Benzene (4n + 2 system, n = 1)

Benzene is a six-membered carbon ring. Solving its secular determinant yields the roots:

x = -2, -1, -1, +1, +1, +2

The energy levels are:

E₁ = α + 2β
E₂, E₃ = α + β (Degenerate bonding levels)
E₄, E₅ = α - β (Degenerate antibonding levels)
E₆ = α - 2β

With 6 π-electrons filling the bonding levels, the total π-energy is:

E_π = 2(α + 2β) + 4(α + β) = 6α + 8β

Comparing this to three isolated ethylene molecules (6α + 6β):

DE = E_π - E_isolated = (6α + 8β) - (6α + 6β) = 2β

Since β ≈ -75 kJ/mol, this yields a immense delocalization stabilization energy (around 150 kJ/mol), explaining benzene's high chemical stability and classic aromatic character.