Configuration Interaction

Configuration Interaction (CI) bezeichnet eine Methode zur Lösung der Schrödinger-Gleichung (bzw. ihrer relativistischen Verallgemeinerungen), die besonders in der Quantenchemie verwendet wird. Die Vielteilchen-Wellenfunktion wird dabei in eine Basis aus Slater-Determinanten entwickelt, wodurch die Schrödinger-Gleichung auf ein Matrix-Eigenwertproblem reduziert wird. Die (teilweise) Diagonalisierung dieser Matrix liefert dann die Eigenzustände des quantenmechanischen Systems.

Basisentwicklung, Slater-Determinanten

Die Schrödingergleichung

$ {\hat {H}}|\Psi \rangle \ =\ E|\Psi \rangle $

stellt eine Operatorengleichung für abstrakte Vektoren in einem Hilbertraum dar. Zu deren Lösung wählt man eine bestimmte Darstellung der Wellenfunktion. Eine Einteilchenwellenfunktion stellt man z.B. dar durch Entwicklung in eine Basis Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): | \phi_k \rangle der Größe Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): N_b eines Einteilchenhilbertraums Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathcal H ,

$ |\Psi \rangle \ =\ \sum _{k}c_{k}\,|\phi _{k}\rangle $

N-Teilchenwellenfunktionen sind Funktionen aus Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathcal H_N \ = \ \mathcal H \times \mathcal H \times \cdots \times \mathcal H , dem N-fachen kartesischen Produkt eines Einteilchen-Hilbertraums. Eine Basis von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathcal H_N ist durch alle möglichen Produkte der Einteilchenbasis gegeben, so dass die Wellenfunktion wie folgt entwickelt werden kann:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): | \Psi \rangle \ = \ \sum_{k_1 \cdots k_N} c_{k_1 \cdots k_N} \, | \phi_{k_1} \cdots \phi_{k_N} )

wobei man die Basisvektoren

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): | \phi_{k_1} \cdots \phi_{k_N} ) \ = \ | \phi_{k_1} \rangle \cdots | \phi_{k_N} \rangle

Hartree-Produkte nennt.

Aufgrund des Pauliprinzips muss die elektronische Wellenfunktion antisymmetrisch gegenüber Vertauschung zweier Teilchenkoordinaten sein, d.h. Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): | \Psi \rangle lebt nur in dem Unterraum Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathcal H_N^- der antisymmetrischen Funktionen. Die Hartree-Produkte erfüllen diese Forderung nicht, weswegen auch die Wellenfunktion nicht antisymmetrisch sein muss. Um die Antisymmetrisierung zu gewährleisten, kann man die Wellenfunktion auf Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathcal H_N^- projizieren. Weitaus häufiger jedoch projiziert man bereits vorher die Basisvektoren auf Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathcal H_N^- , wodurch man aus den Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): N_b^{\, N} Hartree-Produkten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \binom{2N_b}{N} Slater-Determinanten erhält,

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): | \phi_{k_1} \cdots \phi_{k_N} \rangle \ = \ \frac{1}{\sqrt{N!}} \, \sum_{\sigma \in {\mathcal S}_N} \text{sign}(\sigma) | \phi_{\sigma(k_1)} \rangle \cdots | \phi_{\sigma(k_N)} \rangle

wobei die Summe über alle möglichen Permutationen geht. Durch die Slater-Determinanten erhält man eine geeignete Basis zur Entwicklung der Wellenfunktion,

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): | \Psi \rangle \ = \ \sum_{k_1 < k_2 < \cdots < k_N} c_{k_1 \cdots k_N} \, | \phi_{k_1} \cdots \phi_{k_N} \rangle

Slater-Determinanten sind Eigenfunktionen des projizierten Spins $ {\hat {S}}_{z} $, jedoch im Allgemeinen keine Eigenfunktionen des Gesamtspins Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \hat S^2 . In der Praxis wählt man deshalb häufig auch Configuration State Functions (CSF) als Basisfunktionen. Eine CSF lässt sich als eine Linearkombinationen von einigen wenigen Slater-Determinanten angeben. Ihr Vorteil liegt darin, dass die Wellenfunktion automatisch Eigenfunktion des Spins ist, und dass man weniger CSFs als Determinanten zur Entwicklung braucht. Es sollte jedoch erwähnt werden, dass die zur Zeit erfolgreichsten CI Codes mit Slater-Determinanten arbeiten.

Full Configuration Interaction

Die Configuration Interaction Methode erhält man nun sehr einfach. Man setzt die Entwicklung der Wellenfunktion in die Schrödingergleichung ein,

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \sum_{k_1 <\cdots <k_N} \hat H c_{k_1 \cdots k_N} | \phi_{k_1} \cdots \phi_{k_N} \rangle \ = \ E \sum_{k_1 <\cdots <k_N} c_{k_1 \cdots k_N} | \phi_{k_1} \cdots \phi_{k_N} \rangle

und multipliziert sie mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \langle \phi_{j_1} \cdots \phi_{j_N} | . Wegen der Orthonormalität der Slater-Determinante (folgt aus der orthonormalen Einteilchenbasis) erhält man

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \sum_{k_1 <\cdots <k_N} c_{k_1 \cdots k_N} \langle \phi_{j_1} \cdots \phi_{j_N} | \hat H | \phi_{k_1} \cdots \phi_{k_N} \rangle \ = \ E c_{j_1 \cdots j_N}

und damit ein Matrix-Eigenwertproblem,

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf H \mathbf c \ = \ E \mathbf c

In der Quantenchemie ist der Hamiltonian häufig gegeben durch

$ {\hat {H}}=\sum _{pq}h_{pq}\sum _{\sigma }{\hat {a}}_{p\sigma }^{\dagger }{\hat {a}}_{q\sigma }+{\frac {1}{2}}\sum _{pqrs}g_{pqrs}\sum _{\sigma \tau }{\hat {a}}_{p\sigma }^{\dagger }{\hat {a}}_{r\tau }^{\dagger }{\hat {a}}_{s\tau }{\hat {a}}_{q\sigma }\,, $

d.h. als Summe aus Einteilchentermen (kinetische + potentielle Energie) sowie der Zweiteilchen-Coulomb-Wechselwirkung. Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \sigma und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \tau bezeichnen die Spinvariablen.

Um das Eigenwertproblem zu bestimmen, müssen Matrixelemente der Form

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \langle \phi_{j_1} \cdots \phi_{j_N} | \hat H | \phi_{k_1} \cdots \phi_{k_N} \rangle

berechnet werden. Die Auswertung dieser Matrixelemente geschieht mit den Slater-Condon-Regeln.

Eigenschaften

Die Methode ist im Prinzip exakt, die einzige Näherung besteht in der Wahl einer endlich großen Einteilchenbasis. Eine große Einschränkung ist allerdings durch die Skalierung der Hamiltonmatrix gegeben. Für eine gewählte Anzahl an Teilchen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): N und Anzahl an Basisfunktionen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): N_b hat die Matrix die Dimension Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \binom {2N_b} N . Durch Ausnutzung von Symmetrien, z.B. $ [{\hat {H}},{\hat {S}}_{z}]=[{\hat {H}},{\hat {S}}^{2}]=...=0 $ kann diese Zahl zwar reduziert werden, die exponentielle Skalierung bleibt aber bestehen.

In der Praxis verwendet man deswegen iterative Methoden zur Lösung des Eigenwertproblems (z. B. Arpack), oder andere Minimierungsmethoden (z.B. Formen des Newton-Verfahrens), mit denen man nur einige wenige Eigenfunktionen erhält, typischerweise den Grundzustand.

In vielen Fällen wird dabei die Hamiltonmatrix nicht explizit gebildet, sondern nur ihre Wirkung auf den Koeffizientenvektor berechnet, eine Variante die man "Direct CI" nennt.

Die CI-Methode ist size consistent, d.h. die Energie zweier Untersysteme ist immer gleich der Energie des Gesamtsystems.

Einbettung in die Quantenchemie

Verwandte Methoden sind:

• Coupled Cluster (CC),
• self-consistent field method (SCF), siehe Hartree-Fock-Methode,
• Møller-Plesset-Störungsrechnung (MP) sowie
• Multiconfiguration Self-Consistent-Field Algorithmen (MCSCF).