How LUMO Energy Is Calculated

From structure to orbital energies

Whatever the software, a LUMO calculation follows the same sequence. It starts from a description of the molecule, builds a three-dimensional structure, sets up and solves the electronic structure equations for that structure, and then reads off the orbital energies.

Step 1 and 2: input and geometry

A SMILES string describes connectivity, not shape. To calculate orbitals, a program must first place every atom, including the hydrogens that SMILES leaves implicit, in three dimensions. The usual approach is to generate a reasonable starting conformer from rules and torsion libraries, then relax it with a molecular mechanics force field such as MMFF94 or UFF.

In careful work, the geometry is then reoptimized with the same quantum method used for the orbitals. Bond lengths matter: stretching a C=C bond lowers its π* orbital, so a poor geometry shifts the LUMO. For flexible molecules, different conformers can give different LUMO energies, and a single conformer may not be representative.

Step 3: the electronic structure method

Simple Hückel theory

The oldest method treats only the π electrons of planar conjugated molecules. Each atom contributes one p orbital, overlap between orbitals is ignored, and the energies come out in units of two empirical parameters, α and β. It is ideal for teaching and for comparing conjugated hydrocarbons, but it gives no absolute energies in eV without further assumptions.

Extended Hückel theory

Roald Hoffmann's Extended Hückel theory (1963) includes all valence electrons and the overlap between orbitals. Each valence atomic orbital is a Slater-type function with a parameterized exponent, the diagonal Hamiltonian elements are fixed valence ionization energies, and the off-diagonal elements follow the Wolfsberg-Helmholz approximation:

There is no self-consistent field: the equations are solved once. That makes the method fast and free of convergence problems, and it reproduces orbital shapes and many qualitative trends well. Its absolute energies are tied to the parameters, and LUMO energies are typically much lower than those from other methods. This is the method this site's calculator uses; the Methodology page gives every detail.

Semi-empirical SCF methods

Methods such as AM1, PM3, PM6, PM7 and the GFN-xTB family solve a self-consistent field problem with many integrals approximated or parameterized. They are much faster than Hartree-Fock and are widely used for screening. Their orbital energies depend on the parameterization.

Hartree-Fock

Hartree-Fock theory solves the Roothaan equations self-consistently, with each electron moving in the average field of the others. It has no empirical parameters, but it omits electron correlation. Virtual orbitals, including the LUMO, are computed in the field of all N electrons, which tends to push their energies too high; with diffuse basis functions the LUMO of a neutral molecule is often positive.

Density functional theory

Kohn-Sham density functional theory includes correlation through an approximate exchange-correlation functional. Orbital energies depend strongly on the functional: hybrids such as B3LYP give lower LUMO energies than Hartree-Fock and gaps that are usually smaller, while range-separated and tuned functionals change the picture again. LUMO values from different functionals are not interchangeable.

Basis sets

Except in the empirical methods, the molecular orbitals are built from a basis set of atom-centered functions. Larger basis sets describe the orbitals more flexibly. For the LUMO, diffuse functions (the "+" in 6-31+G* or the "aug-" in aug-cc-pVDZ) have a particularly large effect, because an added electron is loosely held and spreads far from the nuclei. Adding diffuse functions often lowers or raises the LUMO by more than an electronvolt and can change which orbital is lowest.

Charge, spin and environment

The total charge sets the number of electrons and therefore which orbital is the LUMO. The spin multiplicity decides whether electrons are paired: a radical has a singly occupied orbital and separate α and β LUMOs. Most calculations treat the molecule as isolated in the gas phase. Implicit solvent models such as PCM stabilize charged and polar states and shift orbital energies, so a gas-phase LUMO and a solvated one should not be compared directly.

How this site does it

This calculator runs the whole pipeline in your browser. OpenChemLib parses the input and generates a conformer, the MMFF94s force field relaxes it, and an Extended Hückel solver with Hoffmann's Slater-type parameters computes the orbitals. Results report the method variant, the number of basis functions, the geometry source and the engine version, so the number can be reproduced.

References and further reading

1. Hückel, E. Quantentheoretische Beiträge zum Benzolproblem. Z. Phys. 1931, 70, 204-286.
2. Hoffmann, R. An Extended Hückel Theory. I. Hydrocarbons. J. Chem. Phys. 1963, 39, 1397-1412.
3. Wolfsberg, M.; Helmholz, L. The Spectra and Electronic Structure of the Tetrahedral Ions MnO₄⁻, CrO₄²⁻, and ClO₄⁻. J. Chem. Phys. 1952, 20, 837-843.
4. Roothaan, C. C. J. New Developments in Molecular Orbital Theory. Rev. Mod. Phys. 1951, 23, 69-89.
5. Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133-A1138.
6. Halgren, T. A. MMFF VI. MMFF94s Option for Energy Minimization Studies. J. Comput. Chem. 1999, 20, 720-729.
7. Jensen, F. Introduction to Computational Chemistry, 3rd ed. Wiley, 2017.