MO and VB approaches to chemical bonding in diatomics and Huckel theory for conjugated systems

Introduction to Chemical Bonding Theories

Understanding how atoms bond together to form molecules is a cornerstone of chemistry. Two prominent theoretical frameworks have been developed to explain chemical bonding: Valence Bond (VB) theory and Molecular Orbital (MO) theory. While both aim to describe the electronic structure of molecules, they approach the problem from different perspectives. VB theory focuses on the overlap of atomic orbitals to form localized bonds, whereas MO theory considers the combination of atomic orbitals to form delocalized molecular orbitals that span the entire molecule. These theories are particularly powerful when applied to simple systems like diatomic molecules and conjugated systems.

Valence Bond (VB) Theory

Valence Bond theory, developed by Walter Heitler and Fritz London, and later extended by Linus Pauling, describes a chemical bond as the result of the attraction between the nuclei and the overlap of atomic orbitals. When two atoms approach each other, their atomic orbitals containing unpaired electrons can overlap. This overlap creates a region of increased electron density between the two nuclei, leading to a strong attractive force that holds the atoms together. The strength of the bond is directly related to the extent of this overlap.

In VB theory, the electrons involved in bonding are considered to remain localized in the region between the two bonded atoms. For a simple diatomic molecule like H₂, the VB approach considers the overlap of the 1s atomic orbitals of the two hydrogen atoms. Each hydrogen atom has one electron in its 1s orbital. When the two atoms approach, their 1s orbitals overlap, and the two electrons, with opposite spins, become paired in this overlapping region, forming a covalent bond.

VB theory also incorporates the concept of hybridization, where atomic orbitals on the same atom mix to form new hybrid orbitals. These hybrid orbitals have different shapes and orientations, which are more suitable for forming stronger bonds with specific spatial arrangements. For example, in methane (CH₄), the carbon atom undergoes sp³ hybridization to form four equivalent tetrahedral hybrid orbitals, which then overlap with the 1s orbitals of four hydrogen atoms to form four sigma (σ) bonds.

For diatomic molecules, VB theory can accurately predict bond lengths and strengths, especially for homonuclear diatomics (molecules composed of two identical atoms) where the electron distribution is relatively symmetrical. However, it can become more complex to apply to heteronuclear diatomics (molecules with different atoms) due to differences in electronegativity and orbital energies.

Molecular Orbital (MO) Theory

Molecular Orbital theory, pioneered by Robert Mulliken, offers a different perspective. Instead of focusing on localized atomic orbitals, MO theory proposes that atomic orbitals combine to form new molecular orbitals that are delocalized over the entire molecule. These molecular orbitals can be either bonding orbitals (lower in energy, increasing bond stability) or antibonding orbitals (higher in energy, decreasing bond stability).

The formation of molecular orbitals can be understood using the Linear Combination of Atomic Orbitals (LCAO) approximation. According to this approximation, molecular orbitals (ψ_MO) are formed by adding or subtracting the wave functions of the atomic orbitals (ψ_AO):

For a bonding molecular orbital (ψ_bonding), the atomic orbitals combine constructively (in-phase): ψ_bonding = c₁ψ_A + c₂ψ_B

For an antibonding molecular orbital (ψ_antibonding), the atomic orbitals combine destructively (out-of-phase): ψ_antibonding = c₁ψ_A - c₂ψ_B

The coefficients c₁ and c₂ depend on the relative energies and overlap of the atomic orbitals. When atomic orbitals combine, the number of molecular orbitals formed is equal to the number of atomic orbitals that combine.

In a diatomic molecule like H₂, the two 1s atomic orbitals combine to form one bonding molecular orbital (σ_1s) and one antibonding molecular orbital (σ*_1s). The two electrons from the hydrogen atoms fill the lower-energy σ_1s bonding orbital, leading to a stable molecule.

MO theory is particularly useful for explaining the magnetic properties of molecules (paramagnetism and diamagnetism) and predicting the relative stabilities of different molecules. It also provides a more accurate description of bonding in molecules with delocalized electrons, such as conjugated systems.

MO and VB Approaches for Diatomic Molecules

When applied to diatomic molecules, both VB and MO theories provide valuable insights, but they highlight different aspects of bonding.

Valence Bond Approach for Diatomic Molecules

For homonuclear diatomic molecules like H₂, N₂, O₂, and F₂, VB theory effectively describes the formation of sigma (σ) and pi (π) bonds through the overlap of atomic orbitals. For example, in N₂, three bonds are formed: one σ bond from the end-to-end overlap of 2p_z orbitals (assuming the internuclear axis is z), and two π bonds from the side-by-side overlap of 2p_x and 2p_y orbitals. VB theory correctly predicts the triple bond and its strength.

However, VB theory struggles to accurately predict the magnetic properties of O₂. VB theory would suggest that O₂, with a double bond, should be diamagnetic (all electrons paired). In reality, O₂ is paramagnetic, meaning it has unpaired electrons. This is a limitation of the basic VB approach.

Molecular Orbital Approach for Diatomic Molecules

MO theory excels where VB theory falters, particularly with O₂. The MO diagram for O₂ shows that after filling the lower-energy molecular orbitals, the last two electrons occupy two degenerate π*_2p antibonding orbitals, with one electron in each orbital and with parallel spins, according to Hund's rule. This unpaired electron configuration explains the observed paramagnetism of oxygen.

MO theory allows for the calculation of bond order, which is a measure of the number of chemical bonds between two atoms. It is calculated as:

Bond Order = 0.5 * (Number of electrons in bonding MOs - Number of electrons in antibonding MOs)

A higher bond order indicates a stronger and shorter bond. For example, H₂ has a bond order of 1, N₂ has a bond order of 3, and O₂ has a bond order of 2. This aligns well with experimental observations.

For heteronuclear diatomic molecules like CO or HF, MO theory, especially with modifications to account for differing atomic orbital energies and electronegativity, provides a more accurate picture of bond polarity and electronic distribution.

Hückel Molecular Orbital (HMO) Theory for Conjugated Systems

Hückel theory is a simplified version of MO theory specifically developed for delocalized π electron systems in conjugated molecules, such as linear and cyclic polyenes. Conjugated systems are characterized by alternating single and double bonds, where the p-orbitals on adjacent carbon atoms overlap to form a delocalized π system.

The fundamental assumptions of Hückel theory are:

• Only π electrons are considered; σ electrons are assumed to form strong, localized bonds and do not participate in delocalization.
• All carbon atoms in the conjugated system are equivalent in terms of their contribution to the π system.
• The overlap between p-orbitals on non-adjacent atoms is negligible.
• The Coulomb integral (α) for each carbon atom is the same, representing the energy of an electron in an isolated 2p atomic orbital.
• The resonance integral (β) between adjacent carbon atoms is the same, representing the interaction energy between adjacent p-orbitals.
• The resonance integral (β) between non-adjacent carbon atoms is zero.

Using these assumptions, Hückel theory sets up a secular determinant to solve for the energies of the π molecular orbitals and the coefficients of the atomic orbitals in these molecular orbitals. For a conjugated system with 'n' carbon atoms, there will be 'n' π molecular orbitals.

Hückel Theory for Linear Conjugated Systems

Consider a linear conjugated system like butadiene (CH₂=CH-CH=CH₂). It has four carbon atoms, each contributing one 2p orbital to the π system, resulting in four π molecular orbitals. The secular determinant for butadiene is:

In Hückel theory, we often simplify the calculations by setting α = 0 and β = 1 (or any non-zero value). The determinant then simplifies to:

The solutions to this determinant give the energy levels of the π molecular orbitals. For a general linear polyene with 'n' carbon atoms, the energy of the k^th molecular orbital (E_k) is given by:

E_k = α + 2β cos(kπ / (n+1)) , where k = 1, 2, ..., n

For butadiene (n=4):

• k=1: E₁ = α + 2β cos(π/5) ≈ α + 1.618β (bonding)
• k=2: E₂ = α + 2β cos(2π/5) ≈ α + 0.618β (bonding)
• k=3: E₃ = α - 2β cos(2π/5) ≈ α - 0.618β (antibonding)
• k=4: E₄ = α - 2β cos(π/5) ≈ α - 1.618β (antibonding)

Butadiene has 4 π electrons (two from each double bond). These fill the two lowest energy bonding orbitals, E₁ and E₂. The delocalization energy (a measure of the extra stability due to π electron delocalization) can be calculated by comparing the actual π electron energy with the energy if the double bonds were isolated.

Hückel Theory for Cyclic Conjugated Systems (Aromaticity)

Hückel theory is famously used to explain aromaticity in cyclic conjugated systems. For a cyclic system with 'n' carbon atoms, the energy levels are given by:

E_k = α + 2β cos(2kπ / n) , where k = 0, ±1, ±2, ..., ±(n/2-1), n/2 for even n or k = 0, ±1, ±2, ..., ±(n-1)/2 for odd n.

A key outcome of Hückel theory for cyclic systems is the prediction of aromaticity based on the **4n + 2 rule**. A cyclic, planar molecule with a fully conjugated ring system is aromatic if it has (4n + 2) π electrons, where 'n' is a non-negative integer (0, 1, 2, ...).

• n=0: 4(0) + 2 = 2 π electrons (e.g., cyclopropenyl cation)
• n=1: 4(1) + 2 = 6 π electrons (e.g., benzene)
• n=2: 4(2) + 2 = 10 π electrons (e.g., cyclooctatetraene dianion)

Molecules with 4n π electrons are considered antiaromatic, which are typically unstable.

Example: Benzene (C₆H₆)

Benzene is a classic example of a cyclic conjugated system. It has 6 carbon atoms, each contributing one 2p orbital to the π system, forming 6 π molecular orbitals. Applying Hückel theory:

The energy levels are:

• E₁ = α + 2β (bonding, non-degenerate)
• E₂ = α + β (bonding, degenerate)
• E₃ = α + β (bonding, degenerate)
• E₄ = α - β (antibonding, degenerate)
• E₅ = α - β (antibonding, degenerate)
• E₆ = α - 2β (antibonding, non-degenerate)

Benzene has 6 π electrons (one from each double bond). These fill the three lowest energy bonding orbitals: one non-degenerate (α + 2β) and two degenerate (α + β). This results in a highly stable electronic configuration, explaining benzene's aromatic character. The delocalization energy of benzene is significantly higher than what would be expected for three isolated double bonds.

Comparison and Applications

Both MO and VB theories are powerful tools, and their choice often depends on the specific problem. VB theory provides a more intuitive picture of localized bonds and is often used in introductory organic chemistry to explain reaction mechanisms. MO theory, especially in its more advanced forms and Hückel's simplified approach, offers a more quantitatively accurate description of electronic structure, magnetic properties, and spectral transitions, particularly for conjugated and delocalized systems.

Hückel theory, despite its simplifications, provides invaluable qualitative understanding of the electronic properties of conjugated hydrocarbons, predicting stability, reactivity, and spectral properties. It forms the basis for understanding more complex conjugated systems and organic electronic materials.