Huckel Molecular Orbital (HMO) Calculations for Butadiene

The Butadiene System ($C_4H_6$)

1,3-Butadiene represents a conjugated linear framework containing four carbon $2p$ atomic orbitals in a continuous chain. The determinantal form of the secular equation characterising this open system is given by:

$$\begin{vmatrix} x & 1 & 0 & 0 \\ 1 & x & 1 & 0 \\ 0 & 1 & x & 1 \\ 0 & 0 & 1 & x \end{vmatrix} = 0$$

Expanding this $4 \times 4$ matrix yields a fourth-degree characteristic polynomial equation:

$$x^4 - 3x^2 + 1 = 0 \tag{1}$$

Solving this polynomial equation reveals four distinct roots for the system:

$$x = \pm 1.618, \quad \pm 0.618$$

These roots map directly onto four quantized $\pi$ energy levels:

$$\begin{aligned} E_1 &= \alpha + 1.618\beta \quad \text{(Strongly Bonding MO)} \\ E_2 &= \alpha + 0.618\beta \quad \text{(Weakly Bonding MO)} \\ E_3 &= \alpha - 0.618\beta \quad \text{(Weakly Antibonding MO)} \\ E_4 &= \alpha - 1.618\beta \quad \text{(Strongly Antibonding MO)} \end{aligned}$$

Four $\pi$-electrons occupy these orbitals in pairs within the ground-state configuration ($\psi_1^2 \psi_2^2$). The total net $\pi$-electron energy ($E_\pi$) for butadiene is calculated as:

$$E_\pi = 2(\alpha + 1.618\beta) + 2(\alpha + 0.618\beta) = 4\alpha + 4.472\beta$$

Since the baseline energy of four isolated electrons residing across two separate non-interacting ethylene units evaluates to $2 \times (2\alpha + 2\beta) = 4\alpha + 4\beta$, the net Delocalization Energy (DE) or resonance stabilization energy gained by conjugation in butadiene is computed as:

$$\text{DE} = (4\alpha + 4.472\beta) - (4\alpha + 4\beta) = 0.472\beta$$

Molecular Orbital Wave Functions

By substituting the roots back into the linear combinations, we calculate the precise numerical coefficients for each of the four normalized $\pi$ molecular orbitals. Arranged in order of increasing energy, they are expressed as:

$$\begin{aligned} \psi_1 &= 0.3755(\phi_1 + \phi_4) + 0.6070(\phi_2 + \phi_3) \\ \psi_2 &= 0.6070(\phi_1 - \phi_4) + 0.3755(\phi_2 - \phi_3) \\ \psi_3 &= 0.6070(\phi_1 + \phi_4) - 0.3755(\phi_2 - \phi_3) \\ \psi_4 &= 0.3755(\phi_1 - \phi_4) - 0.6070(\phi_2 - \phi_3) \end{aligned}$$

Nodal Symmetry and Orbital Planes

The $\pi$-molecular orbitals always possess the horizontal molecular skeleton plane as a primary nodal plane. In addition to this, they contain vertical nodal planes passing perpendicular to the molecular plane between the carbon atoms. As predicted by quantum boundary conditions, the number of vertical nodal planes increases sequentially with the energy of the system:

• $\psi_1$ (Lowest Energy): 0 vertical nodes. Entirely bonding across all centers.
• $\psi_2$ (HOMO): 1 vertical node located right at the central C₂–C₃ bond.
• $\psi_3$ (LUMO): 2 vertical nodes located symmetrically between C₁–C₂ and C₃–C₄.
• $\psi_4$ (Highest Energy): 3 vertical nodes alternating between every single adjacent carbon pair.

Calculations of Electron Density, Bond Order, and Free Valence

Evaluating the ground-state electron configuration ($\psi_1^2 \psi_2^2$) yields the following highly informative structural values for the conjugated molecule:

$$\begin{aligned} \text{Electron Densities:} \quad & q_1 = q_2 = q_3 = q_4 = 1.0 \\ \text{\(\pi\)-Bond Orders:} \quad & p_{12} = p_{34} = 0.912; \quad p_{23} = 0.436 \\ \text{Free Valence Indices:} \quad & F_1 = F_4 = 0.820; \quad F_2 = F_3 = 0.384 \end{aligned}$$

Chemical Interpretation and Reactivity

These values reveal a set of crucial insights regarding the behavior of conjugated dienes:

• Bond Length Discrepancy: Although the overall $\pi$-charge distribution ($q_i$) is completely uniform across all four carbon centers in the ground state, the calculated central $\pi$-bond order ($p_{23} = 0.436$) is drastically lower than the outer double bonds ($p_{12} = 0.912$). This confirms that the middle link is weaker, longer, and maintains significant single-bond character, highlighting that assuming a uniform resonance integral ($\beta$) for all bonds oversimplifies the system.
• Free Valence and Diels-Alder Additions: The free valence index displays that the central carbon atoms have significantly less partial unshared valency than the terminal ends ($F_{\text{terminal}} = 0.820$ vs. $F_{\text{central}} = 0.384$). Consequently, the terminal atoms 1 and 4 are far more susceptible to radical attacks and chemical interactions. This perfectly rationalizes why 1,3-butadiene undergoes highly selective 1,4-addition reactions during electrocyclic processes like the Diels-Alder reaction.

Related Topics
Huckel Molecular Orbital (HMO) Calculations for Ethylene
Huckel Molecular Orbital (HMO) Calculations for Allyl System