HMO-Modell

Die von Erich Hückel vorgeschlagene Hückel-Näherung, auch: Hückel-Molekülorbital-Modell/Methode (HMO-Modell/Methode), ist eine Methode der semiempirischen Quantenchemie, die mit wenig Rechenaufwand erstaunlich gute Resultate liefert. Sie erlaubt es, Molekülorbitale in konjugierten Systemen zu approximieren. Die beiden wichtigsten Folgerungen sind die Hückel-Regel und die Woodward-Hoffmann-Regeln.

Anwendung

Die Näherung leitet sich vom Ritz'schen Verfahren ab. Dabei wird formal das Eigenwertproblem $\mathbb H _m \Psi = E \Psi$ reduziert auf $H \vec{x}=E \vec{x}$. Im ersten Fall ist Ψ eine Eigenfunktion des Molekül-Hamiltonoperators $\mathbb H _m$. Um sie genau zu bestimmen müsste bei einem n-atomigen Molekül eine 6n-dimensionale partielle Differentialgleichung gelöst werden, was analytisch nicht möglich ist. Im zweiten Fall ist $\vec{x}$ der n-dimensionale Vektor der die Koeffizienten zur Linearkombination enthält.

Nachdem $\vec{x}$ berechnet wurde, kann das Molekülorbital als entsprechende Linearkombination der einzelnen p_z-Orbitale angegeben werden. Der Wert E gibt die Energie des Orbitals wieder.

Im Ritz'schen Verfahren wird nach der Hartree-Fock Methode iterativ eine möglichst gute Lösung gefunden. Dafür müssen in jedem Schritt eine große Zahl von Integralen gelöst werden, was hohen Rechenaufwand bedeutet. Die Vereinfachung der Hückelnäherung ist, dass alle Integrale parametrisiert werden.

Die n konjugierten Atome im Molekül werden durchnummeriert. Die Matrix H = (h_ij) ist eine n x n Matrix. Man setzt:

h_ii = α_i

h_ij = β, falls die beiden Atome benachbart sind (und über Konjugation verknüpft)

h_ij = 0 sonst

α_i ist dabei das Coulombintegral des Atoms i im Molekül (α_i<0)
$\alpha_i=\int_{V}~\Psi_i^* \mathbb H _m \Psi_i~dV$
(* ist dabei die komplexe Konjugation, V bezeichnet das ganze Volumen)

β das Resonanzintegral zwischen zwei Atomen i,j (wird für alle Atompaare als gleich angenommen) (β<0)
$\beta=\int_{V}~\Psi_j^* \mathbb H _m \Psi_i~dV (i \neq j, beliebig)$

Der Sinn der Näherung ist, dass diese beiden Integrale nicht berechnet werden. Sie können zum Beispiel anhand von empirischen Daten geschätzt werden. Für zwei Atome von derselben Art werden die Coulombintegrale gleichgesetzt. Besonders einfach kann man daher reine Kohlenwasserstoffe berechnen. Es bleiben nur zwei Konstante α und β übrig. Die Eigenvektoren sind unabhängig von deren Wert.

Ableitung der Hückel-Regel

Die Frost-Kreis bietet eine einfache Möglichkeit, die Stabilität cyclisch-konjugierter ebener Moleküle abzuschätzen. Sie basiert auf der überraschenden Tatsache, dass das Energieniveauschema eines cyclisch-konjugierten ebenen n-atomaren Moleküls als regelmäßiges n-Eck dargestellt werden kann.

Die zu einem n-Ring gehörige n x n Matrix wird gebildet, indem man in die Hauptdiagonale α schreibt, links und rechts daneben β und in die links untere und rechts obere Ecke β:

$\begin{pmatrix} \alpha &amp;amp;amp; \beta &amp;amp;amp; 0 &amp;amp;amp; ... &amp;amp;amp; 0 &amp;amp;amp; \beta \\ \beta &amp;amp;amp; \alpha &amp;amp;amp; \beta &amp;amp;amp; 0 &amp;amp;amp; ... &amp;amp;amp; 0 \\ 0 &amp;amp;amp; \beta &amp;amp;amp; \alpha &amp;amp;amp; \beta &amp;amp;amp; ... &amp;amp;amp; 0 \\ ... \\ \beta &amp;amp;amp; 0 &amp;amp;amp; ...&amp;amp;amp; 0 &amp;amp;amp; \beta &amp;amp;amp; \alpha \end{pmatrix}$

Die Eigenwerte dieser Matrix ergeben sich zu $\alpha + 2 \beta cos(2 \pi k / n); 0 \leq k \leq n-1$. Trägt man nun die Energie entlang der y-Achse auf, beachtet, dass β negativ ist und hält entsprechende Abstände in x-Richtung ein, so erhält man ein auf der Spitze stehendes n-Eck.

Dieses n-Eck kann nun mit einer ungeraden Zahl von Elektronenpaaren (4n+2 π-Elektronen) unter hohem Energiegewinn besetzt werden. Diese Moleküle werden als "Aromaten" bezeichnet. Bei einer geraden Zahl von Elektronenpaaren (4n π-Elektronen) wären die beiden obersten Niveaus halbbesetzt, das Molekül paramagnetisch und der gesamte Energiegewinn gering. Diese Verbindungen sind instabil und heißen "Antiaromaten".

Beispiel

Benzol

Die Hückelmatrix lautet: Die Eigenwerte ergeben sich zu: (α + 2 β, α + β, α + β, α - β, α - β, α - 2 β) Die Eigenvektoren zu den jeweiligen Eigenwerten: ((1, 1, 1, 1, 1, 1), (1, 0, -1, -1, 0, 1), (-1, -1, 0, 1, 1, 0), (-1, 1, 0, -1, 1, 0), (-1, 0, 1, -1, 0, 1), (-1, 1, -1, 1, -1, 1)) Ein Schnitt (parallel zu xy-Ebene in 75 pm Entfernung) durch die entsprechenden Orbitale ist rechts gezeichnet. Die Koordinaten der Vektoren bedeuten jeweils das Vorzeichen der pz-Wellenfunktionen. Diese sind angefangen von rechts gegen den Uhrzeigersinn durchgezählt. Gemäß einem auf der Spitze stehenden 6-eck gibt es 3 stabilisierte und drei destabilisierte Orbitale. 6 Elektronen können also unter hohem Energiegewinn die unteren drei Orbitale besetzen. Die Gesamtbindungsenergie ist 2*2β+4*β=8β. Dieser Wert ist deutlich höher als 6β für drei isolierte π-Bindungen.

Quellen

• Peter W. Atkins: Physikalische Chemie. Weinheim [u.a.] : VCH-Verl.-Ges.
• Dietmar Dorninger: Mathematische Grundlagen für Chemiker. Eisenstadt : Prugg