Hückel-Näherung
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 $ \Psi $ 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 pz-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:
hii = αi
hij = β, falls die beiden Atome benachbart sind (und über Konjugation verknüpft)
hij = 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\quad (i\neq j,{\text{ 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
Der 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 &\beta &0&...&0&\beta \\\beta &\alpha &\beta &0&...&0\\0&\beta &\alpha &\beta &...&0\\...\\\beta &0&...&0&\beta &\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.
- Beispiele für aromatische Verbindungen
_with_MO_diagram.png)
_with_MO_diagram.png)
_with_MO_diagram.png)
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".
- Beispiele für antiaromatische Verbindungen
_with_MO_diagram.png)
_with_MO_diagram.png)
Beispiel
Benzol
Die Hückelmatrix lautet: $ {\begin{pmatrix}\alpha &\beta &0&0&0&\beta \\\beta &\alpha &\beta &0&0&0\\0&\beta &\alpha &\beta &0&0\\0&0&\beta &\alpha &\beta &0\\0&0&0&\beta &\alpha &\beta \\\beta &0&0&0&\beta &\alpha \\\end{pmatrix}} $
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