Energy-Level Diagrams and Molecular Orbitals for Conjugated Polyenes
Hückel molecular orbital (HMO) energy diagrams can be constructed using linear combinations of -orbitals for conjugated polyenes. Quantum chemistry packages have been used to model the molecules and molecule-orbitals via RHF/6-31G(d) or U-B3LYP/6-31G(d).
This Demonstration shows characteristics for a set of seven conjugated polyenes: the HMO energy diagram, a 3D view of the probability density (with boundary surface, phase-coloring included), and a ball-and-stick model for each example.
The simple Hückel method (SMO) is a rudimentary method for energies and orbitals of -electron systems.
Possible sets of corresponding LCAO molecular orbitals for cyclobutadiene and benzene are
To visualize the different linear combinations of wavefunctions for -electron systems of conjugated polyenes, a shaded version of the orbital is used . For example, the coefficients in the linear combinations for cyclobutadiene are all equal by size and differ only by sign. In benzene they can be ordered by size into three groups. In an HMO energy diagram for benzene, a combination of three different relative sizes of is used. The last column in the energy diagram shows no symbol or the symbols ↑ and ⥮ to indicate that the orbital is not occupied, singly occupied, or doubly occupied.
ChemSpider  is the basis for the coordinates in the input files of the nonradical structures. The ChemSpider IDs for each structure are listed in the following table.
An input file for the allyl- and 1,4-pentadienyl-radical is created in Avogadro beginning with a sketch .
The non-radicals' optimizations are run in GAMESS-US with an RHF/6-31G(d) model via GamesQ .
For both radicals, an optimization with unrestricted Hartree–Fock is necessary. The model of choice is U-B3LYP/6-31G(d).
Here is a table for the final energy from the log-files of the optimizations for each molecule.
fin. energy [Hartree]
From each log-file the last frame is exported as a mol-file via MacMolPlt . The Demonstration "Displaying Molecules with Multiple Bonds" by Bianca Eifert makes it possible to show the double-bonds of molecules in Mathematica . Through slightly modifying her code, the modeled molecules can be imported as a mol-file with double-bonds on.
The corresponding molecule-orbitals for the molecules are created in Gabedit .
The results are copied to this Demonstration into the initalization code and named MOOrbitals and MOMolecules.
 M. Hanwell, D. Curtis, D. Lonie, T. Vandermeersch, E. Zurek, and G. Hutchison, "Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform," Journal of Chemoinformatics, 4(17), 2004 pp. 1758–2946. doi:10.1186/1758-2946-4-17.
 M. Schmidt, K. Baldridge, J. Boatz, S. Elbert, M. Gordon, J. Jensen, S. Koseki, N. Matsunaga, K. Nguyen, S. Su, T. Windus, M. Dupuis, and J. Montgomery, "General Atomic and Molecular Electronic Structure System," Journal of Computational Chemistry, 14, 1993 pp. 1347–1363. doi:10.1002/jcc.540141112.
 B. Bode and M. Gordon, "MacMolPlt," Journal of Molecular Graphics and Modelling,16, 1998 pp. 133–138.