imported>Anagkai: Assoziative Verweise entfernt: Verweis in Einleitung vorhanden.

2026-02-07T21:14:34Z

Assoziative Verweise entfernt: Verweis in Einleitung vorhanden.

Neue Seite

Unter '''Hartree-Fock-Rechnung''' (beziehungsweise '''Hartree-Fock-Methode''', nach [[Douglas Hartree]]<ref>{{Literatur |Autor=D. R. Hartree |Titel=The Wave Mechanics of an Atom with a Non-Coulomb Central Field. Part I. Theory and Methods |Sammelwerk=Mathematical Proceedings of the Cambridge Philosophical Society |Band=24 |Nummer=1 |Datum=1928-01 |Sprache=en |ISSN=0305-0041 |DOI=10.1017/S0305004100011919 |Seiten=89–110 |Kommentar=Part II: 10.1017/S0305004100011920; Part III: 10.1017/S0305004100015954; Part IV: 10.1017/S0305004100014031 |Online=https://www.cambridge.org/core/product/identifier/S0305004100011919/type/journal_article |Abruf=2024-12-09}}</ref> und [[Wladimir Alexandrowitsch Fock]]<ref>{{Literatur |Autor=V. Fock |Titel=Näherungsmethode zur Lösung des quantenmechanischen Mehrkörperproblems |Sammelwerk=Zeitschrift für Physik |Band=61 |Nummer=1–2 |Datum=1930-01 |ISSN=1434-6001 |DOI=10.1007/BF01340294 |Seiten=126–148 |Online=http://link.springer.com/10.1007/BF01340294 |Abruf=2024-12-09}}</ref>) versteht man eine Methode der [[Quantenmechanik]], in der Systeme mit mehreren gleichartigen [[Teilchen]] in [[Molekularfeldnäherung|Mean-Field-Näherung]] behandelt werden. Sie wird zum Beispiel verwendet in der [[Atomphysik]], der [[Theoretische Chemie|Theoretischen Chemie]] zur Beschreibung von [[Elektron]]en in [[Molekül]]en und der [[Kernphysik]] für Systeme aus [[Proton]]en und [[Neutron]]en.

Sie ermöglicht es, [[Atomorbital|Orbital]]<nowiki />energien und [[Wellenfunktion]]en von quantenmechanischen [[Vielteilchensystem]]en [[Approximation|näherungsweise]] zu berechnen und ist eine so genannte [[Ab initio|Ab-initio]]-Methode, d. h. sie kommt ohne [[empirisch]]e Parameter aus und benötigt nur [[Naturkonstante]]n. Sie ist der Ausgangspunkt für [[Post-Hartree-Fock-Methoden]], welche die Genauigkeit der Berechnungen verbessern.

Die Hartree-Fock-Methode ist die Basis der [[Molekülorbitaltheorie]].

== Funktionsweise ==
Die Hartree-Fock-Methode geht von der zeitunabhängigen [[Schrödinger-Gleichung]] (hier in [[Dirac-Notation]])

:<math>\hat H \, |\psi\rangle = E \, |\psi\rangle</math>

aus, welche die Energie <math>E</math> eines Systems aus der Wellenfunktion <math>|\psi\rangle</math> berechnet, indem die [[Eigenwerte]] des [[Hamilton-Operator]]s <math>\hat H</math> zu dieser Wellenfunktion gesucht werden. Im Hamilton-Operator werden alle Energiebeiträge der Teilchen und Felder im System sowie deren Wechselwirkungen untereinander beschrieben. In vielen praktisch wichtigen Systemen (wie z. B. den [[Elektron]]en in einem [[Molekül]]) sind die Teilchen miteinander [[Quantenverschränkung|korreliert]] und beeinflussen sich gegenseitig. Dadurch kann die Schrödinger-Gleichung für solche Systeme nicht mehr [[Analytische Lösung|exakt]], sondern nur noch [[Näherungslösung|näherungsweise]] gelöst werden.

Die Hartree-Fock-Methode vereinfacht die Wechselwirkungen der Teilchen untereinander so, dass diese nicht mehr jeweils paarweise untereinander wechselwirken, sondern mit einem [[Feld (Physik)|Feld]], das von allen anderen Teilchen im Mittelwert erzeugt wird – dem so genannten ''mean field'' (mittleren Feld). Das Feld hängt zwar immer noch vom Verhalten der einzelnen Teilchen ab, die Lösung kann aber jetzt schrittweise berechnet werden:
* Ein Ausgangszustand wird ausgewählt und daraus das Feld erzeugt.
* Mit diesem wird dann die Schrödingergleichung für jedes einzelne Teilchen gelöst.
* Zusammengenommen ergeben die einzelnen Lösungen dann einen neuen Zustand und ein neues Feld.
Dieser Vorgang wird wiederholt, bis sich aufeinanderfolgende Lösungen nur mehr geringfügig unterscheiden, das Feld also zu Lösungen führt, die das Feld selbst konsistent wieder erzeugen. Daraus leitet sich der Begriff ''self-consistent field'' ab, der für diesen Teil der Hartree-Fock-Methode verwendet wird.

Als Wellenfunktionen für die behandelten Vielteilchensysteme werden bei [[Boson]]en ein symmetrisches Produkt von Einteilchenwellenfunktionen '''(Hartree-Produkt)''' verwendet, bei [[Fermion]]en (wie Elektronen, Protonen und Neutronen) eine antisymmetrische Kombination dieser Produkte (eine sogenannte [[Slater-Determinante]]). Um die Schrödingergleichung zu lösen, werden diese Einteilchenwellenfunktionen so variiert, dass die aus der Gleichung entstehende Energie minimal wird. Aufgrund des [[Rayleigh-Ritz-Prinzip]]s ist diese Energie dann eine obere Grenze für die tatsächliche Energie des Systems. Die dadurch berechnete Wellenfunktion des gesamten Systems ist allerdings nicht notwendigerweise eine Annäherung der tatsächlichen Wellenfunktion.

Bei manchen Molekülen (insbesondere mit [[Radikal (Chemie)|ungepaarten Elektronen]]) wird statt einer einzigen [[Slater-Determinante]] eine symmetrieadaptierte [[Linearkombination]] mehrerer Slater-Determinanten angesetzt, deren [[Koeffizient]]en aber durch die ([[Spin]]-)[[Symmetrie (Physik)|Symmetrie]] des Systems festgelegt sind.

== Hartree-Fock-Gleichung ==
Die Hartree-Fock-Gleichung ist ein [[nichtlinear]]es [[Eigenwertproblem]] mit einem nichtlokalen Integrodifferentialoperator. Sie lautet in [[Dirac-Notation]]
{| cellpadding="5" style="border:2px solid #50C878; background:#ECFCF4; margin-left:2em; text-align:center;"
|
<math> \hat{F}|\phi_m\rangle = \varepsilon_m |\phi_m\rangle</math>
|}
mit dem [[Fock-Operator]] <math> \hat{F}</math>
: <math>
\hat{F} = \hat{h} +\sum\limits_{\gamma}^{N} \left( \left\langle\phi_{\gamma}\right|\hat{w}\left|\phi_{\gamma}\right\rangle-\left|\phi_{\gamma}\right\rangle\left\langle\phi_{\gamma}\right|\hat{w} \right)
</math>
: <math>
\text{mit } \quad \hat{h}= - \frac{\Delta(\mathbf{r})}{2}-\sum\limits_k^{N_k}\frac{Z_k}{\left|{\mathbf{r}-\mathbf{R_k}}\right| }\quad \text{ und } \quad\hat{w}=\frac{1}{\left|{\mathbf{r_1}-\mathbf{r_2}}\right| }
</math>
wobei <math>\hat{h}</math> der Einteilchenanteil des Hamiltonoperators ist und <math>\hat {w}</math> der Anteil der Zweiteilchenwechselwirkung, wie oben erwähnt für den Spezialfall der [[Molekülphysik]] von Elektronen mit [[Coulombwechselwirkung]] untereinander und in [[Atomare Einheiten|atomaren Einheiten]].

Der Index <math> \gamma </math> läuft hierbei über die besetzten elektronischen Zustände, also die mit den <math>N</math> niedrigsten Eigenwerten, wobei <math>N</math> die Zahl der Elektronen angibt. Der Index <math> k </math> läuft über die Atomkerne, wobei <math>N_k</math> die Anzahl der Kerne angibt.
=== Matrixdarstellung ===
Um effiziente Lösungsmethoden verwenden zu können, wird die Gleichung zusätzlich in eine [[Matrizenmechanik|Matrixdarstellung]] übergeführt, indem man <math>|\phi_m\rangle</math> in der [[Basis (Vektorraum)|Basis]] <math>|\varphi_i\rangle</math> darstellt, sodass <math>|\phi_m\rangle = \sum\limits_i^n c_{im}|\varphi_i\rangle </math>. Diese Basis ist typischerweise nicht orthogonal.
:<math>
\sum\limits_i \hat{F}c_{im}|\varphi_i\rangle = \varepsilon_m\sum\limits_i c_{im} |\varphi_i\rangle
</math>
Nach Multiplikation mit <math> \langle \varphi_j|</math> ergibt sich das [[Verallgemeinertes Eigenwertproblem|verallgemeinerte Eigenwertproblem]]
:<math>
\sum\limits_i \langle \varphi_j|\hat{F}|\varphi_i\rangle c_{im}=\varepsilon_m\sum\limits_i \langle \varphi_j|\varphi_i\rangle c_{im}
</math>

{| cellpadding="5" style="border:2px solid #50C878; background:#ECFCF4; margin-left:2em; text-align:center;"
|
<math> F \mathbf{c}_m = \varepsilon_m S \mathbf{c}_m</math>
|}
mit der Fockmatrix <math>F_{ji}=\langle \varphi_j|\hat{F}|\varphi_i\rangle </math>, der Überlappmatrix <math>S_{ji}=\langle \varphi_j|\varphi_i\rangle </math> und den Koeffizientenvektoren <math>\mathbf{c}_{m}</math>. Diese Gleichung ist auch als [[Roothaan-Hall-Gleichungen|Roothaan-Hall-Gleichung]] bekannt. Wird die Basis diagonalisiert (z. B. mit [[Per-Olov Löwdin|Löwdins]] [[Symmetrische Orthogonalisierung|symmetrischer Orthogonalisierung]]), wodurch die Überlappmatrix zu einer Einheitsmatrix wird, vereinfacht sich die Gleichung zu einem einfachen Eigenwertproblem, das von Computern effizient gelöst werden kann.

Als Lösung erhält man <math>n</math> Eigenwerte und Eigenvektoren, wovon man die <math>N</math> niedrigsten Eigenwerte und zugehörigen Eigenvektoren als besetzte Zustände ansieht. Als Basisfunktionen kommen in vielen Fällen [[Linearkombination]]en von [[Gaussian Type Orbitals]] (GTO) oder [[Slater Type Orbitals]] (STO) zum Einsatz. Für Berechnungen an einzelnen Atomen und zweiatomigen oder linearen Molekülen können die Hartree-Fock-Gleichungen auch mit numerischen Verfahren gelöst werden.

== Spin ==
Um die Hartree-Fock-Gleichung zu lösen, muss von den oben verwendeten Spinorbitalen <math>\left|\phi_m\right\rangle</math> noch die Spinwellenfunktion
<math>\left|\chi_m\right\rangle</math> abgespalten werden, sodass <math>\left|\phi_m\right\rangle=\left|\psi_m\right\rangle\left|\chi_m\right\rangle </math>
mit der reinen Ortswellenfunktion <math>\left|\psi_m\right\rangle</math> gilt.
=== Geschlossene-Schalen-Hartree-Fock (RHF) ===
Bei dem Geschlossene-Schalen-Hartree-Fock Ansatz (engl. Restricted Hartree Fock) werden alle Spins als gepaart angenommen,
was natürlich nur bei einer geraden Anzahl von Elektronen möglich ist. Der Grundzustand wird somit als [[Multiplizität|Spin-Singulett]] angenommen.
Für die Wellenfunktionen folgt somit
:<math>
\left|\phi_1\right\rangle = \left|\psi_1\right\rangle \left|\alpha\right\rangle, \; \left|\phi_2\right\rangle = \left|\psi_1\right\rangle \left|\beta\right\rangle,\;\left|\phi_3\right\rangle = \left|\psi_2\right\rangle \left|\alpha\right\rangle, \dots , \left|\phi_N\right\rangle = \left|\psi_{N/2}\right\rangle \left|\beta\right\rangle.
</math>
:<math>
\text{mit } \quad \hat{s}_z \left|\alpha\right\rangle = \frac{1}{2}\left|\alpha\right\rangle \quad \text{ und } \quad \hat{s}_z \left|\beta\right\rangle = -\frac{1}{2}\left|\beta\right\rangle
</math>
Setzt man dies in die Hartree-Fock-Gleichung ein, folgt
:<math>
\hat{F}_{\text{RHF}} \left|\psi_m\right\rangle = \hat{h} \left|\psi_m\right\rangle +\sum\limits_{\gamma}^{N/2}2 \left\langle\psi_\gamma\right|\hat{w}\left|\psi_\gamma\right\rangle\left|\psi_m\right\rangle-\left\langle\psi_\gamma\right|\hat{w}\left|\psi_m\right\rangle\left|\psi_\gamma\right\rangle=\varepsilon_m \left|\psi_m\right\rangle.
</math>
Die Coulombwechselwirkung tritt somit zwischen allen Elektronen auf, die [[Austauschwechselwirkung]] hingegen nur zwischen Elektronen
mit gleichem Spin. Wegen der Symmetrie zwischen Spin up und down ist die HF-Gleichung für beide Spinkonfigurationen gleich,
sodass weiterhin nur eine Eigenwertgleichung gelöst werden muss, wobei nun allerdings nur noch die niedrigsten <math>N/2</math> Eigenwerte
und Eigenvektoren verwendet werden müssen.

=== Offene-Schalen-Hartree-Fock (UHF) ===
Bei dem Offene-Schalen-Hartree-Fock-Ansatz (engl. Unrestricted Hartree Fock) wird im Vergleich zum Geschlossene-Schalen-Ansatz (RHF) die Forderung fallengelassen, dass gleich viele Elektronen im Zustand <math>\left|\alpha\right\rangle</math>, wie im Zustand
<math>\left|\beta\right\rangle </math> sein müssen. Die Spinorbitale werden demnach angesetzt als
:<math>
\left|\phi_1\right\rangle = \left|\psi_1^\alpha\right\rangle \left|\alpha\right\rangle, \; \left|\phi_2\right\rangle = \left|\psi_2^\alpha\right\rangle \left|\alpha\right\rangle , \dots , \left|\phi_{N_\alpha}\right\rangle = \left|\psi_{N_\alpha}^\alpha\right\rangle \left|\alpha\right\rangle,\; \left|\phi_{N_\alpha+1}\right\rangle = \left|\psi_1^\beta\right\rangle \left|\beta\right\rangle , \dots , \left|\phi_{N}\right\rangle = \left|\psi_{N_\beta}^\beta\right\rangle \left|\beta\right\rangle.
</math>
Nach Einsetzen in die ursprüngliche Hartree-Fock-Gleichung ergeben sich zwei verschiedene Gleichungen für <math>\left|\alpha\right\rangle</math>
und <math>\left|\beta\right\rangle</math>.
:<math>
\hat{F}_{\text{UHF}}^\alpha \left|\psi_m^\alpha\right\rangle = \hat{h} \left|\psi_m^\alpha \right\rangle +\sum\limits_{\gamma}^{N_\alpha} \left\langle\psi^\alpha_\gamma\right|\hat{w}\left|\psi^\alpha_\gamma\right\rangle\left|\psi^\alpha_m\right\rangle-\left\langle\psi^\alpha_\gamma\right|\hat{w}\left|\psi^\alpha_m\right\rangle\left|\psi^\alpha_\gamma\right\rangle+\sum\limits_{\gamma}^{N_\beta} \left\langle\psi^\beta_\gamma\right|\hat{w}\left|\psi^\beta_\gamma\right\rangle\left|\psi^\alpha_m\right\rangle=\varepsilon_m \left|\psi^\alpha_m\right\rangle.
</math>
Die Gleichung für <math>\left|\beta\right\rangle</math> folgt aus der Ersetzung <math>\alpha\rightarrow\beta</math> und <math>\beta\rightarrow\alpha</math>.
Hierbei sieht man wieder, dass Elektronen mit gleichem Spin Coulomb- und Austauschwechselwirkung besitzen, Elektronen mit
unterschiedlichem Spin wechselwirken hingegen nur über den Coulombterm. Da die Austauschwechselwirkung die Gesamtenergie
stets verringert, kann somit, im Rahmen von Hartree-Fock, die zweite [[Hundsche Regeln|Hundsche Regel]] erklärt werden.
Diese besagt, dass bei sonstiger Entartung oder Quasientartung die Spins zweier Elektronen möglichst parallel ausgerichtet sind.

== Herleitung für Fermionen ==
Zur Herleitung der Hartree-Fock Gleichungen geht man zunächst von der stationären [[Schrödingergleichung]] aus. Hier wird der Spezialfall eines [[Hamiltonoperator]]s mit Coulombwechselwirkung in der [[Born-Oppenheimer-Näherung]] betrachtet, wie er zum Beispiel für Elektronen in der Molekülphysik auftritt. Das heißt
:<math>
\hat{H}=\sum\limits_{i}^{N}\underbrace{\left(- \frac{\Delta(\mathbf{r}_i)}{2}-
\sum\limits_k^{N_k}\frac{Z_k}{\left|{\mathbf{r_i}-\mathbf{R_k}}\right| } \right)}_{\hat{h}_i}+\frac{1}{2}\sum\limits_{i}^{N}\sum\limits_{j\neq i}^{N}\underbrace{\frac{1}{\left|{\mathbf{r_i}-\mathbf{r_j}}\right| }}_{\hat{w}_{ij}}
</math>
<math>\mathbf{r}</math> bezeichnet hierbei die elektronischen Koordinaten, <math>N</math> die Anzahl der Elektronen, <math>Z_k</math> und <math>\mathbf{R}_k</math> die Ladung und festen
Koordinaten der Kerne. <math>\hat{h}_i</math> ist nun ein Einteilchenoperator und besteht aus der kinetischen Energie und der Wechselwirkung mit allen Kernen des <math>i</math>-ten Elektrons. <math>\hat{w}_{ij}</math> ist hingegen ein Zweiteilchenoperator und stellt die Coulombwechselwirkung des {{nowrap|1=<math>i</math>-ten}} mit dem <math>j</math>-ten Elektron dar. Die stationäre Schrödingergleichung lautet nun
:<math>
\hat{H}\left|\Psi\right\rangle = E \left|\Psi\right\rangle
</math>
Als Näherung für Hartree-Fock schreibt man nun <math>\left|\Psi\right\rangle</math> als [[Slater-Determinante]] von
Einteilchenwellenfunktionen <math>\left|\phi_i\right\rangle</math>. Die Näherung besteht darin, dass man für die exakte Lösung
über alle möglichen Slater-Determinanten summieren müsste, z. B. indem man <math>\phi_N</math> durch <math>\phi_{N+1}</math> ersetzt. Somit gilt
:<math>
\left|\Psi\right\rangle \approx \frac{1}{\sqrt{N!}}
\begin{vmatrix} \phi_1(1) & \phi_2(1) & \dots & \phi_N(1) \\
\phi_1(2) & \phi_2(2) & \dots & \phi_N(2) \\
\vdots & \vdots & \ddots & \vdots \\
\phi_1(N) & \phi_2(N) & \dots & \phi_N(N)
\end{vmatrix}
</math>
und die Energie des Systems lautet
:<math>
E = \left\langle\Psi\right| \hat{H}\left|\Psi\right\rangle
</math>
Dies kann man nun, indem man die [[Orthogonalität]] der <math>\phi_i</math> ausnutzt, zu
:<math>
E[\{ \phi_i\}] =\sum\limits_{\alpha}^{N}\int\phi_\alpha^*(\mathbf{r}_i) \hat{h}_i \phi_\alpha(\mathbf{r}_i) \,\mathrm{d^3}\mathbf{r}_i+
\frac{1}{2}\sum\limits_\alpha^N\sum\limits_{\gamma\neq\alpha}^N\left( \iint \phi_\alpha^*(\mathbf{r}_i) \phi_\gamma^*(\mathbf{r}_j) \hat{w}_{ij} \,\phi_\gamma(\mathbf{r}_j) \phi_\alpha(\mathbf{r}_i) \,d^3\mathbf{r}_i\,d^3\mathbf{r}_j \right.
</math>

:<math>
\quad-\left. \iint \phi_\alpha^*(\mathbf{r}_i) \phi_\gamma^*(\mathbf{r}_j) \hat{w}_{ij} \,\phi_\gamma(\mathbf{r}_i) \phi_\alpha(\mathbf{r}_j) \,d^3\mathbf{r}_i\,d^3\mathbf{r}_j \right)
</math>

:<math>
=\sum\limits_{\alpha}^{N} \left\langle\phi_\alpha\right| \hat{h}\left|\phi_\alpha\right\rangle+
\frac{1}{2}\sum\limits_\alpha^N\sum\limits_{\gamma\neq\alpha}^N\left( \left\langle \phi_\alpha \phi_\gamma\right| \hat{w}\left|\phi_\gamma \phi_\alpha\right\rangle- \left\langle \phi_\alpha \phi_\gamma\right| \hat{w}\left|\phi_\alpha\phi_\gamma \right \rangle \right)
</math>
umformen. Nun wird das [[Rayleigh-Ritz-Prinzip|Ritzsche Variationsprinzip]] verwendet und <math>E</math> als [[Funktional]] nach <math>\phi_i</math> variiert. Um die Orthogonalität der Einteilchenfunktionen zu erhalten, wird allerdings nicht direkt <math>E</math> minimiert, sondern nach der Methode der [[Lagrange-Multiplikator]]en
das Funktional <math>\mathcal{L}</math>
:<math>
\mathcal{L} = E[\{\phi_i \}]-\sum\limits_i^N\sum\limits_{j}^N \varepsilon_{ij}\left(\delta_{ij}- \left\langle \phi_i\right|\left.\phi_j\right\rangle \ \right).
</math>
Man kann nun in die Basis <math>\tilde{\phi}_i</math> wechseln, in der <math>\varepsilon_{ij}</math> diagonal ist, also <math>\varepsilon_{ij}=\varepsilon_i\delta_{ij}</math>.
:<math>
\mathcal{L} = E[\{\tilde{\phi}_i \}]-\sum\limits_i^N \varepsilon_{i}\left(1- \left\langle \tilde{\phi}_i\right|\left.\tilde{\phi}_i\right\rangle \right)
</math>
Die Tilde wird im Weiteren weggelassen. Nun kann <math>\mathcal{L}</math> bezüglich <math>\phi_m</math> minimiert werden.
:<math>
\frac{\partial\mathcal{L}}{\partial \phi_m^*} = \frac{\partial}{\partial\phi_m^*}\left\{\sum\limits_{\alpha}^{N} \left\langle\phi_\alpha\right| \hat{h}\left|\phi_\alpha\right\rangle+
\frac{1}{2} \sum\limits_\alpha^N\sum\limits_{\gamma\neq\alpha}^N\left( \left\langle \phi_\alpha \phi_\gamma\right| \hat{w}\left|\phi_\gamma \phi_\alpha\right\rangle- \left\langle \phi_\alpha \phi_\gamma\right| \hat{w}\left|\phi_\alpha\phi_\gamma \right\rangle\right) -\sum\limits_i^N \varepsilon_{i}\left(1- \left\langle \phi_i\right|\left.\phi_i\right\rangle \ \right) \right\}
</math>
:<math>
= \hat{h} \left|\phi_m\right\rangle + \sum\limits_{\gamma\neq m}^{N} \left( \int \phi_\gamma^*(\mathbf{r}_1)\hat{w}_{12}\phi_\gamma(\mathbf{r}_1) \phi_m(\mathbf{r}_2) \,\mathrm{d^3}\mathbf{r}_1 - \int \phi_\gamma^*(\mathbf{r}_1)\hat{w}_{12}\phi_m(\mathbf{r}_1) \phi_\gamma(\mathbf{r}_2) \,\mathrm{d^3}\mathbf{r}_1 \right) -\varepsilon_m \left|\phi_m\right\rangle
</math>
:<math>
= \hat{h} \left|\phi_m\right\rangle +\sum\limits_{\gamma\neq m}^{N} \left( \underbrace{\left\langle\phi_\gamma\right|\hat{w}\left|\phi_\gamma\right\rangle\left|\phi_m\right\rangle}_{\text{Coulomb-WW.}} - \underbrace{\left\langle\phi_\gamma\right|\hat{w}\left|\phi_m\right\rangle\left|\phi_\gamma\right\rangle}_{\text{Austausch-WW.}} \right) -\varepsilon_m \left|\phi_m\right\rangle
</math>
Da der Summand mit <math>\gamma=m</math> gleich Null ist, kann er hinzugenommen werden, wodurch alle <math>N</math> Gleichungen identisch sind und somit der Index <math>m</math> weggelassen werden kann.
:<math>
= \hat{h} \left|\phi_m\right\rangle +\sum\limits_{\gamma}^{N} \left( \left\langle\phi_\gamma\right|\hat{w}\left|\phi_\gamma\right\rangle\left|\phi_m\right\rangle- \left\langle\phi_\gamma\right|\hat{w}\left|\phi_m\right\rangle\left|\phi_\gamma\right\rangle \right) -\varepsilon_m \left|\phi_m\right\rangle = 0
</math>
Somit folgt
:<math>
\hat{F} \left|\phi_m\right\rangle = \hat{h} \left|\phi_m\right\rangle +\sum\limits_{\gamma}^{N}\left( \left\langle\phi_\gamma\right|\hat{w}\left|\phi_\gamma\right\rangle\left|\phi_m\right\rangle-\left\langle\phi_\gamma\right|\hat{w}\left|\phi_m\right\rangle\left|\phi_\gamma\right\rangle\right)=\varepsilon_m \left|\phi_m\right\rangle,
</math>
die Hartree-Fock-Gleichung mit dem [[Fock-Operator]] <math>\hat{F}</math>. Hierbei besitzen die beiden ersten Terme ein klassisches Analogon. <math>\hat{h}</math> enthält die kinetische Energie und die Coulombwechselwirkung mit den Kernen. Der zweite Term kann als mittleres Coulombpotential aller anderen Elektronen auf das <math>m</math>-te Elektron interpretiert werden. Die instantane Korrelation der Teilchen wird jedoch vernachlässigt. Die Hartree-Fock-Methode ist daher ein [[Molekularfeldnäherung|Mean-Field-Ansatz]]. Der Austauschterm besitzt kein klassisches Analagon. Der Fockoperator für das <math>m</math>-te Elektron enthält die [[Wellenfunktion]]en aller anderer Elektronen, wodurch die Fockgleichungen meist nur mit der Methode der selbstkonsistenten Felder, d. h. iterativ mittels [[Fixpunktiteration]], gelöst werden kann. Zur Konvergenzbeschleunigung kommt hierzu häufig das [[DIIS]]-Verfahren<ref name="DOI10.1016/0009-2614(80)80396-4">{{Literatur |Autor=Péter Pulay |Titel=Convergence acceleration of iterative sequences. the case of scf iteration |Sammelwerk=Chemical Physics Letters |Band=73 |Nummer=2 |Datum=1980-07 |ISSN=0009-2614 |Seiten=393–398 |DOI=10.1016/0009-2614(80)80396-4}}</ref> zum Einsatz.

== Basissätze ==
{{Hauptartikel|Basissatz (Chemie)}}
[[Datei:Gauss sto4.jpeg|mini|Darstellung verschiedener STO-NG-Orbitale und das zugehörige ursprüngliche STO-Orbital]]
Eine direkte numerische Lösung der Hartree-Fock-Gleichung als Differentialgleichung ist bei Atomen und linearen Molekülen möglich. In der Regel werden die Orbitale aber analytisch als [[Linearkombination]]en von Basisfunktionen angesetzt (Basissatz), was wiederum eine Näherung darstellt, die umso besser wird, je größer und intelligenter der Basissatz gewählt wird. Typischerweise bringt jedes Atom im Molekül nun eine vom entsprechenden Basissatz festgelegte Anzahl von Basisfunktionen, die auf ihm zentriert sind, mit. Als grober Ausgangspunkt zur Erstellung solcher Basissätze dienen die analytischen Lösungen des Wasserstoffatoms, welche ein <math>\exp(-\zeta r)</math>-Verhalten für große Kernabstände <math>r</math> zeigen. Ansätze dieses Typs nennt man [[Slater Type Orbitals|Slater Type Orbital]] (STO). Meist haben sie die Form
:<math>\left|\varphi_i^k\right\rangle= \text{Polynom}\cdot \exp\left(-\zeta_i \left|\mathbf{r}-\mathbf{R_k}\right|\right) </math>.
Ein <math>p_z</math>-Orbital besitzt z. B. die Form
:<math>\left|\varphi_i^k\right\rangle= (z-z_k) \cdot \exp\left(-\zeta_i \left|\mathbf{r}-\mathbf{R_k}\right|\right) </math>.
Der große Nachteil der Slater-Type-Orbitale ist jedoch, dass die erforderlichen Matrixelemente <math>\left\langle \varphi_\tau\left|\langle \varphi_\nu\left|\hat{w} \right|\varphi_\eta\rangle\right|\varphi_\gamma\right\rangle</math> nicht im Allgemeinen analytisch berechenbar sind. Deshalb benutzt man fast ausschließlich [[Gaussian Type Orbitals]], d. h. Basisfunktion der Form
:<math> \left|\varphi_i^k\right\rangle= \text{Polynom}\cdot \exp\left(-\zeta_i (\mathbf{r}-\mathbf{R_k})^2\right)</math>.
Hierbei können die Matrixelemente analytisch berechnet werden.<ref>{{Literatur |Autor=Attila Szabo, Neil S. Ostlund |Titel=Modern quantum chemistry: introduction to advanced electronic structure theory |Verlag=Dover Publications |Ort=Mineola NY |Datum=1996 |ISBN=978-0-486-69186-2 |Seiten=}}</ref> Dabei wird u. a. das [[Gaussian Product Theorem]] ausgenutzt, d. h., dass das Produkt zweier Gaußfunktionen wieder eine Gaußfunktion ist. Um die STOs besser zu approximieren, besteht typischerweise eine Basisfunktion aus mehreren Gaußfunktionen mit festen, vom Basissatz festgelegten Parametern („Contraction“). Ein einfacher Basissatz ist z. B. der sog. [[STO-NG-Basissätze|STO-NG]], welcher Slater Type Orbitale mit <math>N</math> Gaußfunktionen annähert. Damit wird die Lösung der Differentialgleichung reduziert auf die analytische Berechnung von Integralen über diese Basisfunktionen und die iterative Lösung des verallgemeinerten Eigenwertproblems mit den Koeffizienten der Basisfunktionen als zu bestimmende Parameter.

Häufig verwendete Basissätze sind die [[Pople-Basen|Pople-]] und die [[Korrelationskonsistente Basen|korrelationskonsistenten Basen]].

== Vor- und Nachteile ==
Die mit der Hartree-Fock-Methode errechnete Energie erreicht nie den exakten Wert, selbst wenn ein unendlich großer Basissatz verwendet werden würde. Bei diesem Grenzfall wird das sogenannte Hartree-Fock-Limit erreicht. Der Grund dafür ist, dass durch die Verwendung des gemittelten Potenzials die [[Elektronenkorrelation]], also die genaue Wechselwirkung der Elektronen untereinander, nicht erfasst wird. Um diesen Makel zu beseitigen, wurden Methoden entwickelt, die in der Lage sind, zumindest einen Teil der Elektronenkorrelation zu erfassen (siehe Artikel [[Korrelierte Rechnungen]]). Von Bedeutung sind insbesondere [[Coupled Cluster|Coupled-Cluster-Methoden]] und die [[Møller-Plesset-Störungstheorie]], die auf der Lösung des Hartree-Fock-Verfahrens aufbauen. Eine andere sehr bedeutende Methode ist die [[Dichtefunktionaltheorie (Quantenphysik)|Dichtefunktionaltheorie]] mit [[Hybridfunktional]]en, bei der der Hartree-Fock-Austausch anteilig in den Austausch-Korrelations-Teil des Dichtefunktionals eingeht.

Die Hartree-Fock-Methode erlaubt aber bei sehr vielen Molekülen eine gute Bestimmung ihrer „groben“ elektronischen Struktur. Daher können z. B. die Molekülorbitale für qualitative Betrachtungen herangezogen werden (z. B. im Falle von [[Grenzorbital]]en). Die Hartree-Fock Methode liefert im Regelfall elektronische Gesamtenergien, die bis auf 0,5 % mit den korrekten elektronischen Energien übereinstimmen (zur Berechnung von Energiedifferenzen, wie z. B. Reaktionsenergien, ist sie aber nur sehr bedingt brauchbar, da diese in der Größenordnung des Fehlers liegen), Dipolmomente, die auf 20 % mit den wirklichen Dipolmomenten übereinstimmen, und sehr genaue Verteilungen der Elektronendichte im Molekül. Aufgrund dieser Eigenschaften werden Hartree-Fock-Rechnungen häufig als Ausgangspunkt für die oben genannten genaueren Rechnungen verwendet.

Ein weiterer Vorteil der Hartree-Fock-Methode ist, dass die erhaltene Energie gemäß dem Variationsprinzip eine obere Schranke für die exakte Grundzustandsenergie darstellt. Durch die Wahl umfangreicherer Basissätze kann die berechnete Wellenfunktion systematisch bis zum sogenannten „Hartree-Fock Limit“ verbessert werden. Eine derartige systematische Betrachtungsweise ist bei [[Dichtefunktionaltheorie (Quantenphysik)|Dichtefunktionalmethoden]] nicht möglich.

== Literatur ==
{{Siehe auch|Theoretische Chemie|Kernphysik}}

=== Originalaufsätze ===
* {{Literatur |Autor=J. C. Slater |Titel=Note on Hartree’s Method |Sammelwerk=Physical Review |Band=35 |Nummer=2 |Datum=1930-01-15 |DOI=10.1103/PhysRev.35.210.2 |Seiten=210–211}}
* {{Literatur |Autor=J. C. Slater |Titel=A Simplification of the Hartree-Fock Method |Sammelwerk=Physical Review |Band=81 |Nummer=3 |Datum=1951-02-01 |Sprache=en |DOI=10.1103/PhysRev.81.385 |Seiten=385–390}}
* {{Literatur |Autor=D. R. Hartree |Titel=The Calculation of Atomic Structures |Verlag=J. Wiley |Ort=New York, New York |Datum=1957 |Reihe=Structure of Matter Series |Online=https://archive.org/details/calculationofato0000hart}}
* {{Literatur |Autor=J. G. Valatin |Titel=Generalized Hartree-Fock Method |Sammelwerk=Physical Review |Band=122 |Nummer=4 |Datum=1961-05-15 |Sprache=en |DOI=10.1103/PhysRev.122.1012 |Seiten=1012–1020}}

== Einzelnachweise ==
<references />

{{Normdaten|TYP=s|GND=4137025-9}}

[[Kategorie:Atomphysik]]
[[Kategorie:Quantenmechanik]]
[[Kategorie:Computerchemie]]
[[Kategorie:Theoretische Chemie]]
[[Kategorie:Kernphysik]]

Hartree-Fock-Methode - Versionsgeschichte

imported>Anagkai: Assoziative Verweise entfernt: Verweis in Einleitung vorhanden.