Symmetrische Orthogonalisierung

Die Symmetrische Orthogonalisierung ist ein von Per-Olov Löwdin (1916–2000) entwickeltes, in der Quantenchemie häufig eingesetztes Orthogonalisierungsverfahren. Als solches dient es dazu, aus einem gegebenen nichtorthogonalen Satz von Vektoren einen orthogonalen Satz zu erzeugen, bei dem für je zwei verschiedene Vektoren das Skalarprodukt gleich Null ist.

Beschreibung

Gegeben sei eine Basis für einen Untervektorraum eines reellen oder komplexen endlichdimensionalen Vektorraums mit Skalarprodukt ( $\mathbb{C}^n$ oder $\mathbb {R} ^{n}$ ). Es sei die Matrix, deren Spaltenvektoren die Basisvektoren von sind.

Man bilde die Gram-Matrix $A^{{\dagger }}A$ . Die Gram-Matrix ist quadratisch, symmetrisch und positiv definit (da die Zeilen von linear unabhängig sind und das Skalarprodukt positiv definit ist) und kann somit unitär diagonalisiert werden. Dabei ist eine unitäre Matrix und eine Diagonalmatrix.

$A^{{\dagger }}A=U^{{\dagger }}DU$

und man kann die Matrix $H:=U^{{\dagger }}D^{{-{\frac {1}{2}}}}U$ bilden. Anschließend bildet man die Matrix ${\tilde {A}}:=AH$ . Die Spaltenvektoren von ${\tilde {A}}$ bilden ein Orthonormalsystem, da:

${\tilde {A}}^{{\dagger }}{\tilde {A}}=(AH)^{{\dagger }}(AH)=(H^{{\dagger }}A^{{\dagger }})(AH)=(HA^{{\dagger }})(AH)=HH^{{-2}}H=Id$

Die Spalten von ${\tilde {A}}$ bilden also die gesuchte Orthonormalbasis von .

Anwendung in der Quantenchemie

In der Quantenchemie führt die approximative, d.h. näherungsweise Lösung der elektronischen Schrödingergleichung auf generalisierte Matrix-Eigenwertprobleme der Form

${\mathbf {FC=SC\epsilon }}$ ,

mit der Fock-Matrix $\mathbf {F}$ , der Koeffizientenmatrix ${\mathbf {C}}$ , welche die LCAO-Koeffizienten der Molekülorbitale enthält und der Diagonalmatrix der Orbitalenergien ${\mathbf {\epsilon }}$ .
Um dieses Eigenwertproblem zu lösen, wird die Gleichung so transformiert, dass die sogenannte Überlappungsmatrix ${\mathbf S}$ zur Einheitsmatrix ${\mathbf E}$ wird. Damit wäre das generalisierte Eigenwertproblem auf ein gewöhnliches Eigenwertproblem

${\mathbf {F'C'=C'\epsilon }}$

reduziert. Dazu werden die Überlappungsmatrix $\mathbf{S}$ mittels einer unitären Transformation zur Matrix ${\mathbf s}$ diagonalisiert, und anschließend die Wurzeln der Kehrwerte der Diagonalelemente gezogen (liefert ${\mathbf s}^{{-1/2}}$ ). Danach wird die Matrix mittels der Rücktransformation wieder "entdiagonalisiert". Mit der so erhaltenen Matrix ${\mathbf {X}}$ und dem Zusammenhang ${\mathbf {C=XC'}}$ kann nun die ursprüngliche Gleichung wie folgt modifiziert werden:

${\mathbf {FXC'=SXC'\epsilon }}$ .

Durch Multiplikation von der linken Seite mit der adjungierten Matrix ${\mathbf {X^{{\dagger }}}}$ erhält man:

${\mathbf {X^{{\dagger }}FXC'=X^{{\dagger }}SXC'\epsilon }}$ .

${\mathbf {X^{{\dagger }}SX}}$ ist aber gerade wieder die Einheitsmatrix, und wir definieren

${\mathbf {F'=X^{{\dagger }}FX}}$ .

Damit erhalten wir das Endergebnis:

${\mathbf {F'C'=C'\epsilon }}$ .

Literatur

A. Szabo, N. S. Ostlund: Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. McGraw-Hill, 1989, ISBN 0-07-062739-8

Basierend auf einem Artikel in:

Wikipedia.de