Symmetrische Orthogonalisierung

Die Symmetrische Orthogonalisierung ist ein von Per-Olov Löwdin 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, d.h. für je zwei solcher Vektoren ist das Skalarprodukt gleich Null.

Beschreibung

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 S zur Einheitsmatrix 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 s diagonalisiert, und anschließend die Wurzeln der Kehrwerte der Diagonalelemente gezogen (liefert 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