The Hückel molecular orbital method is a quantum mechanics approach for calculating the energies of molecular orbitals of π electrons in conjugated hydrocarbon systems, such as ethylene, benzene, and butadiene.

It was proposed by Erich Hückel in 1930, and, subsequently extended and improved but other scientists. It provides the theoretical foundation for Hückel’s rule for the aromaticity of (4n + 2) π electron cyclic, planar systems. These are some of my slides on this topics. In the future, I will add some more explanations. Hückel assumed that for unsaturated organic compounds the π electrons can be treated separately by those involved in the bonds. In fact, a π orbital is antysymmetric for reflection through the plane of the molecule whilst a one is symmetric.

If the assumption is true, in term of energy, we have that . This means that the wavefunction of the molecule is given by the product of the wavefunction describing the and electrons:

We can also assume that the molecular wave function of the whole system can be approximated as a product of 1-electron wave function orbitals

The energy of a sysmt can be evaulated using tha Hamiltonian quantum operator that is defined as

with

- : the kinetic energy of the nucleus.
- : the kinetic energy of the electrons.
- : the proton-electron attraction potential energy.
- : the proton-proton repulsion potential energy.
- : the electron-electron repulsion energy.

In the case of a molecule composed by M atoms, the total hamiltonian can be written as

For a molecular sytem with π electrons, we can further distinguish these electron from the ones as

For not interaction electrons, we can also assume that the total energy of the system can be calculate by a total Hamiltonian (see my blog on the classical mechanics) operator can be expressed as

and the eigenstates of each electron can be calculate by the Schroedingen equations

As further assumption the wavefunction of the electrons can be factorized as follow

Using this approximations the total electronic energy from the electrons can be expressed as

We are now going to consider a simple molecular system

We are now going to consider a simple molecular system to show how to identify the basic atomic orbitals and to construct the Hückel determinant. An allyl molecules has the structural formula and consists of a methylene bridge () attached to a vinyl group (). Its name, derived from the Latin word for garlic, Allium sativum, was given in 1844 by Theodor Wertheim when he was able to isolate an allyl derivative from garlic oil.

We consider only the three orbitals of each carbon since they are antisymmetric for reflection along the molecular plane. The total molecular orbital can be considered as a linear combination of the atomic orbitals.

The total molecular orbital can be considered as a linear combination of the atomic orbitals.

In order to find the energy of the molecular system, we are going to use the result of the so called variational theorem.

The variation theorem states that given a system with a Hamiltonian operator then if is any well-behaved function that satisfies the boundary conditions of the Hamiltonian, then

where is the true value of the lowest energy eigenvalue of This expression is also called **Rayleigh ratio.**

The variational principle allows us to calculate an upper bound for the ground state energy by finding the trial wavefunction for which the integral is minimised. The principle assure us that only if the function is exact in all the other cases .

**Linear variation method**

We are going to use the so-called linear variation method (or Rayleigh-Ritz method) in which the linear variation function is a linear combination of linearly independent atomic orbital functions that satisfy the boundary conditions of the problem (). The coefficients () of the trial wavefunction are parameters to be determined by minimising the variational integral.

By inserting the function of the allyl molecule in the Rayleigh ratio, the numerator become

where is the element of the Hamiltonian matrix.

and the denominator

where is the element of the overlap matrix.

Therefore we can now write

or in the rearranged form as

(1)

Now, we need to find the value of the coefficients that minimize the function using the variational method. This can be obtained using the definition of a mimimum of a function, namely by calculating the equations

for all the functions *i*.

By differentiating both the sides of (1) and using the properties of and and of the coefficients , we finally obtain the so-called secular equation

$\sum_i c_i(H_{ik}-ES_{ik})=0$

for all the *k*.

We therefore have a set of k simultaneous homogeneous equation (secular equations) in k unknowns. In the case of the allyl system the system is the follwing

For a non-trivial solution (i.e. for all i), the determinant of the secular matrix must be equal to zero

Therefore we can set up the values of integral in Hückel determinant as constants and on the bases of the molecular connectivity as

and for the overlap integral

With these semplification the Hückel determinant for the allyl molecule become

Here some more example of Hückel determinat for simple unsaturated molecules

In the second part of this article, we are going to learn how to solve this determinantal equation to find the energy of the molecular orbitals.

**If you have found useful this article, please tell your friends and do not forget to Like it and subscribe for updating!**

**REFERENCES**

J.P. Lowe Quantum Chemistry. 1993, Academic Press.