PhysSandbox
Classical MechanicsWaves & SoundElectricity & MagnetismOptics & LightGravity & OrbitsLabs
🌙Astronomy & The Sky🌡️Thermodynamics🌍Biophysics, Fluids & Geoscience📐Math Visualization🔧Engineering🧪Chemistry

Related simulators

Continue with similar topics in this category — or all 65 in Chemistry.

View category →
NewSchool

Frost Circle & Aromaticity (4n+2)

Launch Simulator

Inscribed n-gon on Hückel π energies; π count vs 4k+2 / 4k; Aufbau filling.

NewSchool

Collision Theory (2D Particles)

Launch Simulator

Hard-disk gas: elastic hits vs an activation speed threshold; T and Eₐ vs exp(−Eₐ/RT).

NewSchool

Radioactive Decay & Chain

Launch Simulator

N(t) from half-life T₁/₂; optional parent → daughter → stable Bateman ODEs (RK4, normalized).

School

Electron Configuration

Launch Simulator

Fill orbitals visually with Aufbau principle animation.

NewUniversity / research

WKB / Bohr–Sommerfeld Quantization

Launch Simulator

Semiclassical bound-state energies for arbitrary V(x): harmonic, quartic, Morse, double-well, asymmetric and square wells. Bisection on the action ∫√(2m(E−V))dx = (n+½)πℏ gives the WKB ladder; compare with the exact harmonic ladder ℏω(n+½).

NewSchool

Nernst Equation

Launch Simulator

E = E° − (RT/nF) ln Q: sliders for E°, n, T, and reaction quotient.

PhysSandbox

Interactive physics, chemistry, and engineering simulators for students, teachers, and curious minds.

Physics

  • Classical Mechanics
  • Waves & Sound
  • Electricity & Magnetism

Science

  • Optics & Light
  • Gravity & Orbits
  • Astronomy & The Sky

More

  • Thermodynamics
  • Biophysics, Fluids & Geoscience
  • Math Visualization
  • Engineering
  • Chemistry

© 2026 PhysSandbox. Free interactive science simulators.

PrivacyTermsContact
Home/Chemistry/Hückel π-MO (Butadiene & Benzene)

Hückel π-MO (Butadiene & Benzene)

Secular matrix H = αI + βA; eigen-energies and LCAO maps on the π skeleton.

π system

0
-1

α is set to 0 as reference; energies are E = α + β μ.

Measured values

ε = (E−α)/β 2.0000·E − α -2.0000

Adjacency matrix A (H = αI + βA)

010001
101000
010100
001010
000101
100010

π MO energies (β < 0)

MO map (LCAO coefficients on carbons)

0.41C10.41C20.41C30.41C40.41C50.41C6

About this model

The Hückel molecular orbital model approximates π conjugation by a tight-binding Hamiltonian on one p_z-like orbital per conjugated carbon: H = α I + β A, where A is the graph adjacency matrix (β for bonded neighbors, zero otherwise). Diagonalizing H yields MO energies E_k = α + β μ_k with μ_k eigenvalues of A, and eigenvector components are LCAO coefficients on each carbon. Butadiene uses a four-site path graph; benzene uses a six-membered cycle, producing degenerate π levels. The page displays A, sorted MO energies for β < 0, and a schematic amplitude map for the selected MO. Electron count, Coulomb differences between atoms, σ–π separation, overlap non-orthogonality, and electron correlation are omitted — this is the textbook π-only introduction.

Who it's for: Organic chemistry alongside MO theory; complements homonuclear diatomic and crystal-field pages.

Key terms

  • Hückel theory
  • secular equation
  • LCAO
  • resonance integral
  • conjugation
  • degenerate orbitals

How it works

Simple Hückel π-electron model: H = α I + β A with A the carbon adjacency matrix (one per p_z basis). Eigenvalues give ε_k = (E_k − α)/β = μ_k (μ are eigenvalues of A). With β < 0, larger μ means lower π energy E = α + β μ. Eigenvectors are LCAO coefficients on each carbon for that MO. Butadiene is a path of four π centers; benzene is a six-cycle. Degeneracies appear for the ring. This is a one-electron, orthogonal-π cartoon — no σ framework, no electron–electron, no non-orthogonality.

Key equations

Secular: det(H − E I) = 0 , H_{ij} = α δ_{ij} + β A_{ij}
E_k = α + β μ_k (μ_k eigenvalue of A)

Frequently asked questions

Why are eigenvalues of A shown as μ = (E−α)/β?
With uniform α on carbons, subtracting α and dividing by β diagonalizes the same eigenvectors as A; μ is dimensionless and β sets the energy scale.
Does this predict benzene bond lengths?
Not quantitatively. Bond-order hints from π occupation can be discussed, but geometry needs going beyond a minimal Hückel parameterization.
Why can math eigenvalues differ slightly from lecture chalkboard numbers?
Numerical diagonalization uses tolerances; degenerate pairs may appear as close splits. The qualitative level pattern and degeneracies match the analytic Hückel solution.