Tuesday, February 15, 2011

Two Java-based Hückel Theory programs

I here share two examples of online Hückel approximation programs. The first is located here and denoted a) in the following, while the second example, b), may be found here. Please note, that b) may be slightly unstable towards refreshing the site but may still be adequate from a teaching point-of-view, though, as it is readily downloaded as a Java-file. Update: Program a) is available in a newer implementation here (thanks to Rafael R. Pappalardo for pointing this out in the Comments).


The first program is indeed very intuitive in the sense that a molecule is easily sketched on the designated canvas and the calculation carried out by clicking just a single button on the right side of the screen. The individual carbons are marked on the canvas by initially clicking the 'Add' button to the left and subsequently placing the atoms by clicking the mouse at the appropriate locations. The bonds (only one type of bond exists, more on this below) are made by dragging a line between two adjacent atoms, equivalent to GView, Avogadro, Jmol, etc. When the molecule has been build and 'optimized' by clicking the 'Normalize' button in the left column, the calculation is made by clicking the 'Show Data Tbl' tab at the upper right corner. The results of the calculation now appears in a seperate window, divided into:
  • MO number
  • Orbital energies
  • Occupancies of the individual MOs
  • Population tables
  • Bond tables w/ appropriate Hückel bond orders
This program furthermore enables the option of illustrating the MOs. A MO diagram is build 'on-the-fly' as the user builds the bonds between the individual atoms, and by clicking the 'Show Molecule' tab in the lower right corner, the MOs are shown in accordance with the MO diagram. Note here, that the orbital energies are shown in the lower part of the screen when the user shifts between the individual MOs by clicking the 'Up' and 'Down' tabs next to the 'Show Molecule' tab.
As touched upon in the above, this program only offers sp2-hybridized orbitals in accordance with regular Hückel theory. The bonds may be altered, though, in an indirect manner, as the alpha- and beta-parameters are modifiable.


As this program has shown to be rather unstable towards running directly from the browser, I would recommend downloading it by clicking the download link below the applet. The file is just a regular .jar file but advantages with respect to the browser applet in the fact that it is faster and more reliable. Furthermore, I will not comment on the Lewis part of the program, only the Hückel part.
In this program, the molecule is build following the same guidelines as for program a). Press 'Build' and draw on the canvas. As bonds are generated automatically in this program, the MO diagram to the right is dynamic as is the total energy, shown in the left column. By pressing the 'Results' button to the left, a seperate window appears, including:

  • The Hückel Hamiltonian (in contrast to program a)). Note that the form of the matrix is in accordance with the numbering within the molecule and may thus not resemble the expected.
  • The orbital energies in terms of the alpha- and beta-parameters.
  • The electron density matrix.
  • The atomic charges as well as the total charge.
As in program a), the degree of hybridization may be altered but in a fancier manner. By clicking the box 'Parameters' in the left side column, and subsequently clicking on one of the individual carbon atoms, the program lets the user change the atom type and hence the final Hamiltonian, etc.
UPDATE: As pointed out by Markus in the comments, you can also display MOs in HuLiS by clicking on specific MO's on the energy graph. Thanks to Markus for this comment.


- Benzene: Use this link (see below) to illustrate the benzene MOs as a supplement to the MOs from program a).
Edit. As Jan mentions in the comments, this section of chemtube3d also serves as a great tool of visualization for the MOs of benzene.

- Naphthalene vs. Azulene. Calculate the differences in the HOMO-LUMO gap. Why does these differences exist and why is azulene dark blue, while naphthalene is colorless?

- Why is it impossible to design and describe 1-methylnaphthalene within the ordinary Hückel approximation? What is the 'Hückel equivalent' to 1-methylnaphthalene?

- Topic of discussion: If we really were to attempt a description of 1-methylnaphthalene, how could the program be manipulated (Hint: Think in terms of alternative ways of describing hybridization).


Jan Jensen said...

Very cool. For a more 3D representation of MOs see this subpage of ChemTube3D, left-hand menu.

That page also has some vibrational subpages. I suggest adding this to the KemiKS DGU page.

Janus J. Eriksen said...

@Jan. Thanks for the link. I have edited the original post!

Anonymous said...

Thanks for these programs. I just have tried them. Now I can solve my problems.
I like the both (SHMO and HuLiS)

Just a notice : In HuLiS, you can display MOs, too, directly clicking on each MO on the energy graph !


Jan Jensen said...

Hi Markus,

Glad to heat you found the post useful, and thanks for the HuLiS tip.

Anonymous said...

Hello Jan, thanks to your answer.

Maybe, you could correct your post, because you wrote : 'As opposed to program b), this program enables the option of illustrating the MOs. '

I hope the both will add 3D representations...

I put your blog in my favorites. Nice job !


Jan Jensen said...

I should hasten to add that Janus Eriksen wrote this post. I have forwarded you comments to him. Glad to hear you like the blog!

Janus J. Eriksen said...


I have made corrections in accordance to your comments in the descriptions of both programs. Thanks for the comments and the nice words.


Rafael R. Pappalardo said...

There is a newer implementation of SHMO at http://wchem.cup.uni-muenchen.de/SHMO/index.html by the same authors.

Rafael R. Pappalardo

Janus J. Eriksen said...

@Rafael. Thanks for the hint. I have updated the Introduction and acknowledged you in the text.