Selbstenergiefunktionaltheorie

Die Selbstenergiefunktionaltheorie ist ein variationelles Berechnungsverfahren der quantenmechanischen Vielteilchentheorie, die Systeme mit einer großen Zahl von miteinander wechselwirkenden Teilchen beschreibt. Das Verfahren wird vorteilhaft eingesetzt bei stark korrelierten (gekoppelten) Elektronensystemen wie sie bei der Magneten und Supraleitern vorgefunden werden.

Anliegen

Die Rechenzeit für die Berechnung von Mehrteilchensystem steigt stark mit der Anzahl der Teilchen N und das mindestens mit N³. Soll beispielsweise die Genauigkeit durch Verdopplung der Anzahl der betrachteten Teilchen verbessert werden, erhöht sich die Wartezeit vor dem Computer auf das Achtfache. Ein gern genutzter Lösungsansatz fasst Gruppen von Teilchen zusammen, approximiert also mit Clustern unter Vereinfachung der Wechselwirkungen der Teilchen in einem Cluster.

Ein anderer Ansatz besteht in der möglichst genauen Berechnung ähnlicher Vergleichssysteme und "Umrechnung" dieser Systeme in das gegebene System. Diesen Ansatz verfolgt die Selbstenergiefunktionaltheorie (SFT). In der Sprache der SFT heißt das: Die Selbstenergiefunktionaltheorie (SFT) liefert einen Ansatz, mit dem recht exakt die Selbstenergie großer Systeme mithilfe kleiner und vor allem mit relativ geringem Rechenaufwand möglichst exakt lösbarer Referenzsysteme berechnet werden kann. Mit diesem Verfahren entkoppelt der Rechenaufwand weitgehend von der Größe des Systems.

Systembeschreibung

Die SFT, erstmals als Selbstenergiefunktionalansatz (SFA) von M. Potthoff 2003 veröffentlicht, stellt ein Variationsprinzip ohne Approximation der funktionalen Abhängigkeit mit der Selbstenergie $\,\Sigma$ als dynamische Basisvariable dar. Dabei wird das aus der Thermodynamik bekannte großkanonische Potential als Funktional der Selbstenergie ($\,\Omega[\Sigma]$) konstruiert, das bei der exakten Selbstenergie des Systems stationär wird. Für die Test-Selbstenergien werden möglichst exakt lösbare Referenzsysteme $\,H = H_0(t)+H_1(U)$ mit der gleichen Wechselwirkung wie das Originalsystem $H^{\prime} = H_0(t^{\prime})+H_1(U)$ benötigt. Hierbei bezeichnen $\,t$ bzw. $t^{\prime}$ den sog. Hopping-Parameter des nicht-wechselwirkenden Anteils des jeweiligen Hamilton-Operators und $\,U$ die Wechselwirkungsstärke.

Genauere Darstellung

(Die folgende Entwicklung benutzt den Greensfunktionsformalismus der Vielteilchentheorie (siehe Weblinks)).

Das großkanonische Potential $\,\Omega$ ist als Funktional der temperaturabhängigen Greensfunktion (Matsubara-Funktion) G gegeben durch:

$\Omega[G]=\Phi[G]+Tr ln(-G)-Tr((G^{-1}_{0}-G^{-1})G),$

wobei die Spur Tr definiert ist als

$Tr A=\beta^{-1}\sum_\omega e^{i\omega0^{+}} tr A$

und $\,\beta$ die inverse Temperatur ist. Das in dem großkanonischen Potential auftretende Luttinger-Ward-Funktional $\,\Phi[G]$ ist dabei definiert als Summe der Beiträge aller geschlossenen, zusammenhängenden, angezogenen Skelett−Diagramme:

$\Phi=\sum_{n=1}^{\infty}\frac{1}{2n\beta}\sum_{\omega} tr(\Sigma^{(n)} G)$

Da i.A. der exakte funktionale Zusammenhang für $\, \Omega$ unbekannt ist, bietet sich der Übergang zum Funktional der Selbstenergie $\,\Sigma$ an:

$\Sigma[G]\rightarrow G[\Sigma]$

Dabei muss angenommen werden, dass $\,\Sigma[G]$ lokal invertierbar ist, das System sich also nicht an einem kritischen Punkt für einen Phasenübergang befindet. Das großkanonische Potential als Funktional der Selbstenergie ist dann:

$\Omega[\Sigma]=Tr \, ln (-(G^{-1}_0-\Sigma)^{-1})+F[\Sigma]$

wobei $\,F[\Sigma]=\Phi[G[\Sigma]]-Tr(\Sigma G[\Sigma])$ die Legendre-Transformierte des Luttinger-Ward-Potentials $\,\Phi[G]$ ist.

Da das Luttinger-Ward-Funktional $\,\Phi[G]$ formal als Summe der Beiträge aller geschlossenen, zusammenhängenden angezogenen Skelett-Diagramme konstruierbar ist, ist dieses nur von der Wechselwirkung abhängig. Somit ist aber auch die Legendre-Transformierte $\,F[\Sigma]$ nur abhängig von der Wechselwirkung. Haben nun die Referenzsysteme die gleiche Wechselwirkung wie das Originalsystem, so ist $\,F$ universell und lässt sich aus dem großkanonischen Potential eliminieren, indem $\,F[\Sigma]$ durch das Referenzsystem ausgedrückt wird:

$F[\Sigma]=\Omega'[\Sigma]-Tr ln(-(G^{'-1}_0-\Sigma)^{-1})$

Damit ist also das großkanonische Potential des Originalsystems:

$\Omega[\Sigma]=\Omega'[\Sigma]+Tr ln(-(G^{-1}_0-\Sigma)^{-1})-Tr ln(-(G^{'-1}_0-\Sigma)^{-1}))$

Stationäritätsbedingung

Die genaue Selbstenergie $\,\Sigma$ erfüllt die folgende Stationäritätsbedingung:

$\frac{\delta\Omega[\Sigma]}{\delta\Sigma}=0,$

d.h. für die exakte Selbstenergie liegt für das großkanonische Potential ein Extremum oder ein Sattelpunkt vor. Diese Stationäritätsbedingung muss lediglich mithilfe der mit anderen Mitteln möglichst exakt berechneten Selbstenergien der Referenzsysteme getestet werden. Wird eine Selbstenergie gefunden die näherungsweise diese Bedingung erfüllt, so hat man eine Näherungslösung für die Selbstenergie und des großkanonischen Potential des Originalsystems.

Literatur

• S. Bäse, Selbstenergiefunktionaltheorie mit stochastischem Cluster-Solver, Universität Hamburg (2008)
• M. Potthoff, Eur. Phys. J. B 32, 429 (2003)
• M. Potthoff, Phys. Rev. Lett. 91, 206402 (2003)
• M. Potthoff, Eur. Phys. J. B 36, 335 (2003)
• C. Dahnken et al., Phys. Rev. B 70, 245110 (2004)

Weblinks

• engl. WP: Greensche Funktion in der Vielteilchentheorie
• Informationsseite zur SFT der Homepage von M. Potthoff