Huckel molecular orbital theory#

Conjugated systems#

  • Molecules with extensive π bonding systems, such as benzene, are not described very well by the valence bond theory because the π electrons are delocalized over the whole molecule. σ and π bonds are demonstrated below for ethylene (C2H4 with sp2 carbons):

  • Note that we have chosen z-axis along the internuclear axis. Because both σ and π bonding occurs between the two carbon atoms, we say that this is a double bond. Note that the hybrid orbitals here also explain the geometry. For triple bonds, one σ and two π bonds are formed.

Huckel MO theory#

  • Huckel molecular orbital theory assumes that the π electrons, which are responsible for the special properties of conjugated and aromatic hydrocarbons, do not interact with one another and the total wavefunction is just a product of the one-electron molecular orbitals. The π molecular orbital of the two carbons in C2H4 can be written approximately as:

ψ=c1ϕ1+c2ϕ2
  • where ϕ1 and ϕ2 are the 2py atomic orbitals for carbon 1 and 2, respectively. By using the variational principle gives the following secular determinant:

|H11ES11H12ES12H21ES21H22ES22|=0 with Hij=ϕiHϕjdτ and Sij=ϕiϕjdτ

In H”uckel theory, the secular equation is simplified by assuming:

  1. All the overlap integrals Sij are set to zero unless i=j, when Sii=1.

  2. All diagonal matrix elements Hii are set to a constant denoted by α.

  3. The resonance integrals Hij (ij) are set to zero except for those on the neighboring atoms, which are set equal to a constant (β). Note that the indices here also identify atoms because the atomic orbitals are centered on atoms.

|αEββαE|=0
  • In Huckel theory, the Coulomb integral α and the resonance integral β are regarded as empirical parameters. They can be obtained, for example, from experimental data. Thus, in the Huckel theory it is not necessary to specify the Hamiltonian operator! Expansion of the determinant in Eq. (???) leads to a quadratic equation for E. The solutions are found to be E=α±β. In general, it can be shown that β<0, which implies that the lowest orbital energy is E1=α+β. There are two π electrons and therefore the total energy is Etot=2E1=2α+2β. Do not confuse α and β here with electron spin.

Solving for Huckel MOs#

The wavefunctions (i.e. the coefficients c1 and c2 can be obtained by substituting the two values of E into the original linear equations:

c1(αE)+c2β=0
c1β+c2(αE)=0

For the lowest energy orbital (E1=α+β), we get (including normalization):

ψ1=12(ϕ1+ϕ2)(i.e. c1=c2=12)

and for the highest energy orbital (E2=αβ) (including normalization):

ψ2=12(ϕ1ϕ2)(i.e. c1=12,c2=12)

These orbitals resemble the H2+ LCAO MOs discussed previously. This also gives us an estimate for one of the excited states where one electron is promoted from the bonding to the antibonding orbital. The excitation energy is found to be 2|β|, which allows for, for example, estimation of β from UV/VIS absorption spectroscopy.

HOMO orbital

= The highest occupied molecular orbital.

LUMO orbital

= The lowest unoccupied molecular orbital.

Example: 1,3-butadiene#

  • Lets alculate the π electronic energy for 1,3-butadiene (CH2=CHCH=CH2) by using the Huckel theory. First we have to write the secular determinant using the rules given earlier. In order to do this, it is convenient to number the carbon atoms in the molecule:

CH21=CH2CH3=CH24

In this case there are two scenarios that should be considered:

  1. A localized solution where the π electrons are shared either with atoms 1 and 2 or 3 and 4. This would imply that the β parameter should not be written between nuclei 2 and 3.

  2. A delocalized solution where the π electrons are delocalized over all four carbons. This would imply that the β parameters should be written between nuclei 2 and 3.

Here it turns out that scenario 2) gives a lower energy solution and we will study that in more detail. In general, however, both cases should be considered. The energy difference between 1) and 2) is called the {resonance stabilization energy}. The secular determinant is:

1234|αEβ00βαEβ00βαEβ00βαE|=0

To simplify notation, we divide each row by β and denote x=(αE)/β:

|x1001x1001x1001x|=0

Expansion of this determinant gives x43x2+1=0. There are four solutions x=±0.618 and x=±1.618. Thus there are four possible orbital energy levels:

E1=α+1.618β(lowest energy)E2=α+0.618βE3=α0.618βE4=α1.618β(highest energy)

There are four π electrons, which occupy the two lowest energy orbitals. This gives the total π electronic energy for the molecule:

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

and the lowest excitation energy is 1.236|β|.

The four Huckel MO wavefunctions are (calculations not shown):

ψ1=0.372ϕ1+0.602ϕ2+0.602ϕ3+0.372ϕ4ψ2=0.602ϕ1+0.372ϕ20.372ϕ30.602ϕ4ψ3=0.602ϕ10.372ϕ20.372ϕ3+0.602ϕ4ψ4=0.372ϕ10.602ϕ2+0.602ϕ30.372ϕ4

Example: Benzene Molecule#

  • Let’s now apply the H”uckel method for benzene molecule. The secular determinant for benzene is (electrons delocalized):

|αEβ000ββαEβ0000βαEβ0000βαEβ0000βαEββ000βαE|=0

The solutions are (where the six π electrons should be placed):

E1=α+2β  (lowest energy)E2=E3=α+βE4=E5=αβE6=α2β  (highest energy)