![]() | ||
The Hückel method or Hückel molecular orbital method (HMO), proposed by Erich Hückel in 1930, is a very simple linear combination of atomic orbitals molecular orbitals (LCAO MO) method for the determination of energies of molecular orbitals of pi electrons in conjugated hydrocarbon systems, such as ethene, benzene and butadiene. It is the theoretical basis for the Hückel's rule. It was later extended to conjugated molecules such as pyridine, pyrrole and furan that contain atoms other than carbon, known in this context as heteroatoms. The extended Hückel method developed by Roald Hoffmann is computational and three-dimensional and was used to test the Woodward–Hoffmann rules.
Contents
- Hckel characteristics
- Hckel results
- Mathematics behind the Hckel method
- Hckel solution for ethylene
- Hckel solution for butadiene
- References
It is a very powerful educational tool, and details appear in many chemistry textbooks.
Hückel characteristics
The method has several characteristics:
Hückel results
The results for a few simple molecules are tabulated below:
The theory predicts two energy levels for ethylene with its two pi electrons filling the low-energy HOMO and the high energy LUMO remaining empty. In butadiene the 4 pi electrons occupy 2 low energy MO's, out of a total of 4, and for benzene 6 energy levels are predicted, two of them degenerate.
For linear and cyclic systems (with n atoms), general solutions exist:
The energy levels for cyclic systems can be predicted using the Frost circle mnemonic. A circle centered at α with radius 2β is inscribed with a polygon with one vertex pointing down; the vertices represent energy levels with the appropriate energies. A related mnemonic exists for linear systems.
Many predictions have been experimentally verified:
Mathematics behind the Hückel method
The Hückel method can be derived from the Ritz method, with a few further assumptions concerning the overlap matrix S and the Hamiltonian matrix H.
It is assumed that the overlap matrix S is the identity matrix. This means that overlap between the orbitals is neglected and the orbitals are considered orthogonal. Then the generalised eigenvalue problem of the Ritz method turns into an eigenvalue problem.
The Hamiltonian matrix H = (Hij) is parametrised in the following way:
Hii = α for C atoms and α + hAβ for other atoms A. Hij = β if the two atoms are next to each other and both C, and kAB β for other neighbouring atoms A and B. Hij = 0 in any other case.The orbitals are the eigenvectors, and the energies are the eigenvalues of the Hamiltonian matrix. If the substance is a pure hydrocarbon, the problem can be solved without any knowledge about the parameters. For heteroatom systems, such as pyridine, values of hA and kAB have to be specified.
Hückel solution for ethylene
In the Hückel treatment for ethylene, the molecular orbital
This equation is substituted in the Schrödinger equation:
to give:
This equation is multiplied by
The same equation is multiplied by
This really can be represented as a matrix. After converting this set to matrix notation,
Or more simply as a product of matrices.
where:
All diagonal Hamiltonian integrals
Other assumptions are that the overlap integral between the two atomic orbitals is 0
leading to these two homogeneous equations:
dividing by
Substituting
This is convenient for computation, but it is also convenient as the energy and coefficients can be easily found:
The trivial solution gives both wavefunction coefficients c equal to zero which is not useful so the other (non-trivial) solution is:
which can be solved by expanding its determinant:
Knowing that
The coefficients can be found by using the previous relationship determined:
Only one equation is necessary however:
The second constant can be replaced giving the following wave equation.
After normalization, the coefficient is obtained:
Leaving
The constant β in the energy term is negative; therefore,
Hückel solution for butadiene
In the Hückel treatment for butadiene, the MO
The secular equation is:
which leads to
and: