Die Kohn-Sham-Dichtefunktionaltheorie, meist als Kohn-Sham-DFT bezeichnet, gehört zu den wichtigsten theoretischen und numerischen Verfahren zur Beschreibung elektronischer Strukturen in Atomen, Molekülen, Festkörpern, Oberflächen und nanoskaligen Systemen. Ihr besonderer Wert liegt darin, dass sie das quantenmechanische Vielteilchenproblem auf eine mathematisch und rechnerisch handhabbare Form reduziert. Statt die vollständige Vielteilchen-Wellenfunktion aller Elektronen direkt zu bestimmen, wird die räumliche Elektronendichte als zentrale Größe verwendet.
Damit verbindet die Kohn-Sham-DFT zwei Ziele, die in der Quantenphysik häufig im Konflikt stehen: eine möglichst realistische Beschreibung der elektronischen Struktur und einen vertretbaren Rechenaufwand. Sie bildet deshalb heute eine wesentliche Grundlage der computergestützten Materialforschung, Quantenchemie, Festkörperphysik und Quantentechnologie.
Bedeutung der elektronischen Struktur
Nahezu sämtliche chemischen und viele physikalische Eigenschaften eines Materials werden durch die Elektronen bestimmt. Bindungsenergien, elektrische Leitfähigkeit, magnetische Eigenschaften, optische Übergänge, Reaktivität, Gitterstabilität und viele weitere Größen hängen von der elektronischen Struktur ab.
Wer die Elektronenverteilung eines Systems kennt, besitzt daher Zugang zu einer großen Zahl messbarer Eigenschaften. Das grundlegende Problem besteht jedoch darin, dass Elektronen miteinander wechselwirken. Jedes Elektron beeinflusst alle anderen Elektronen über die Coulomb-Wechselwirkung. Dadurch entsteht ein Vielteilchenproblem, dessen direkte Lösung bereits bei relativ kleinen Systemen außerordentlich aufwendig wird.
Grundidee der Kohn-Sham-DFT
Die Kohn-Sham-Methode ersetzt das reale System wechselwirkender Elektronen durch ein mathematisches Referenzsystem nicht wechselwirkender Teilchen. Dieses Referenzsystem wird so konstruiert, dass es dieselbe Grundzustandsdichte besitzt wie das reale System.
Die komplizierten Wechselwirkungen werden dabei nicht einfach ignoriert. Sie werden vielmehr in ein effektives Potential und insbesondere in das sogenannte Austausch-Korrelationsfunktional eingebaut. Genau dieser Schritt macht die Kohn-Sham-DFT zugleich leistungsfähig und näherungsabhängig.
Das quantenmechanische Vielteilchenproblem
Die Vielteilchen-Schrödinger-Gleichung
Die elektronische Struktur eines quantenmechanischen Systems wird grundsätzlich durch die Schrödinger-Gleichung beschrieben. Für einen stationären Zustand lautet sie:
\(\hat{H}\Psi = E\Psi\)
Dabei bezeichnet \(\hat{H}\) den Hamiltonoperator, \(\Psi\) die Vielteilchen-Wellenfunktion und \(E\) die Gesamtenergie des Systems.
Für ein System aus Elektronen und Atomkernen enthält der Hamiltonoperator mehrere Beiträge. Dazu gehören die kinetische Energie der Elektronen, die kinetische Energie der Kerne sowie die Coulomb-Wechselwirkungen zwischen Elektronen und Kernen, zwischen Elektronen untereinander und zwischen den Atomkernen.
Born-Oppenheimer-Näherung
Da Atomkerne wesentlich schwerer als Elektronen sind, bewegen sie sich auf deutlich langsameren Zeitskalen. In vielen elektronischen Strukturrechnungen können die Kerne daher zunächst als fest angenommen werden. Diese sogenannte Born-Oppenheimer-Näherung trennt die elektronische Bewegung von der Kernbewegung.
Der elektronische Hamiltonoperator kann dann vereinfacht als
\(\hat{H} = \hat{T} + \hat{V}_{ext} + \hat{V}_{ee}\)
geschrieben werden. Hier bezeichnet \(\hat{T}\) die kinetische Energie der Elektronen, \(\hat{V}_{ext}\) die Wechselwirkung mit dem äußeren Potential der Kerne und \(\hat{V}_{ee}\) die Elektron-Elektron-Wechselwirkung.
Komplexität der Vielteilchen-Wellenfunktion
Für ein System mit \(N\) Elektronen hängt die Wellenfunktion von den räumlichen Koordinaten sämtlicher Elektronen sowie deren Spin ab. Bereits ohne explizite Spinvariablen besitzt sie die Form
\(\Psi(\mathbf{r}_1,\mathbf{r}_2,\ldots,\mathbf{r}_N)\)
und ist damit eine Funktion in einem \(3N\)-dimensionalen Konfigurationsraum.
Mit zunehmender Elektronenzahl wächst der numerische Aufwand deshalb extrem schnell. Genau hier setzt die Dichtefunktionaltheorie an. Anstelle einer hochdimensionalen Vielteilchen-Wellenfunktion verwendet sie die Elektronendichte
\(n(\mathbf{r})\)
die unabhängig von der Zahl der Elektronen nur von drei räumlichen Koordinaten abhängt.
Grundlagen der Dichtefunktionaltheorie
Die Elektronendichte
Die Elektronendichte beschreibt, wie die elektronische Ladung räumlich verteilt ist. Sie kann aus der Vielteilchen-Wellenfunktion gewonnen werden und erfüllt die Normierungsbedingung
\(\int n(\mathbf{r})\,d\mathbf{r} = N\)
wobei \(N\) die Gesamtzahl der Elektronen bezeichnet.
Die Elektronendichte enthält wesentlich weniger Variablen als die Vielteilchen-Wellenfunktion. Die entscheidende Frage lautet daher, ob sie dennoch genügend Informationen enthält, um die grundlegenden Eigenschaften des quantenmechanischen Systems zu bestimmen. Die Hohenberg-Kohn-Theoreme beantworten diese Frage für den Grundzustand grundsätzlich mit Ja.
Energie als Funktional der Dichte
In der Dichtefunktionaltheorie wird die Gesamtenergie als Funktional der Elektronendichte aufgefasst:
\(E = E[n]\)
Ein Funktional ordnet einer vollständigen Funktion, hier also der Elektronendichte \(n(\mathbf{r})\), einen Zahlenwert zu.
Das grundlegende Energiefunktional lässt sich formal schreiben als
\(E[n] = F[n] + \int v_{ext}(\mathbf{r})n(\mathbf{r})\,d\mathbf{r}\)
Dabei enthält \(F[n]\) die kinetische Energie und die Elektron-Elektron-Wechselwirkung. Das äußere Potential \(v_{ext}(\mathbf{r})\) wird insbesondere durch die Atomkerne bestimmt.
Die Hohenberg-Kohn-Theoreme
Erstes Hohenberg-Kohn-Theorem
Das erste Hohenberg-Kohn-Theorem besagt, dass die Grundzustandsdichte das äußere Potential bis auf eine additive Konstante eindeutig festlegt. Daraus folgt, dass die Grundzustandsdichte grundsätzlich alle Informationen enthält, die zur Bestimmung der Eigenschaften des Grundzustands notwendig sind.
Wenn \(n_0(\mathbf{r})\) bekannt ist, kann prinzipiell das äußere Potential bestimmt werden. Damit sind auch der Hamiltonoperator, die Grundzustands-Wellenfunktion und alle aus ihr ableitbaren Observablen festgelegt.
Zweites Hohenberg-Kohn-Theorem
Das zweite Hohenberg-Kohn-Theorem formuliert ein Variationsprinzip. Für jede physikalisch zulässige Elektronendichte \(n(\mathbf{r})\) gilt
\(E[n] \geq E[n_0]\)
wobei \(n_0(\mathbf{r})\) die exakte Grundzustandsdichte bezeichnet.
Die tatsächliche Elektronendichte kann somit durch Minimierung des Energiefunktionals bestimmt werden. Das Problem besteht allerdings darin, dass die exakte Form des universellen Funktionals \(F[n]\) nicht bekannt ist.
Vom Existenztheorem zum praktischen Verfahren
Die Hohenberg-Kohn-Theoreme beweisen die grundsätzliche Existenz einer dichtebasierten Beschreibung, liefern jedoch keine direkte Methode zur Konstruktion des exakten Funktionals. Der entscheidende praktische Fortschritt erfolgte deshalb durch den Ansatz von Walter Kohn und Lu Jeu Sham.
Der Kohn-Sham-Ansatz
Ein nicht wechselwirkendes Referenzsystem
Kohn und Sham führten ein Hilfssystem aus nicht wechselwirkenden Elektronen ein. Dieses künstliche System wird so gewählt, dass seine Grundzustandsdichte genau mit der Grundzustandsdichte des realen wechselwirkenden Systems übereinstimmt.
Für die nicht wechselwirkenden Elektronen können Einteilchenorbitale verwendet werden. Dadurch lässt sich ein großer Teil der kinetischen Energie explizit berechnen.
Zerlegung des Energiefunktionals
Das Kohn-Sham-Energiefunktional wird in mehrere Beiträge zerlegt:
\(E[n] = T_s[n] + \int v_{ext}(\mathbf{r})n(\mathbf{r})\,d\mathbf{r} + E_H[n] + E_{xc}[n]\)
Der Ausdruck \(T_s[n]\) ist die kinetische Energie des nicht wechselwirkenden Kohn-Sham-Referenzsystems. Der zweite Term beschreibt die Wechselwirkung der Elektronen mit dem äußeren Potential.
\(E_H[n]\) ist die klassische elektrostatische Elektron-Elektron-Energie, die Hartree-Energie. Sie lautet
\(E_H[n] = \frac{1}{2}\int\int\frac{n(\mathbf{r})n(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\,d\mathbf{r}\,d\mathbf{r}'\)
Der verbleibende Term \(E_{xc}[n]\) ist die Austausch-Korrelationsenergie.
Bedeutung der Austausch-Korrelationsenergie
Die Austausch-Korrelationsenergie enthält alle komplizierten quantenmechanischen Effekte, die nicht bereits durch die übrigen Terme berücksichtigt werden. Dazu zählen insbesondere der quantenmechanische Austausch, elektronische Korrelationen sowie die Differenz zwischen der tatsächlichen kinetischen Energie des wechselwirkenden Systems und der kinetischen Energie des nicht wechselwirkenden Referenzsystems.
Formal kann man schreiben
\(E_{xc}[n] = T[n] - T_s[n] + E_{ee}[n] - E_H[n]\)
Gerade weil die exakte Form von \(E_{xc}[n]\) unbekannt ist, muss sie in praktischen Berechnungen approximiert werden.
Die Kohn-Sham-Gleichungen
Herleitung durch Variation
Die Kohn-Sham-Gleichungen entstehen durch Minimierung des Energiefunktionals unter der Bedingung, dass die Kohn-Sham-Orbitale orthonormal bleiben. Das Ergebnis ist ein Satz gekoppelter Einteilchengleichungen:
\(\left[-\frac{\hbar^2}{2m_e}\nabla^2 + v_{eff}(\mathbf{r})\right]\phi_i(\mathbf{r}) = \epsilon_i\phi_i(\mathbf{r})\)
Die Größe \(\phi_i(\mathbf{r})\) bezeichnet ein Kohn-Sham-Orbital und \(\epsilon_i\) den zugehörigen Kohn-Sham-Eigenwert.
Effektives Potential
Das effektive Potential setzt sich aus drei wesentlichen Beiträgen zusammen:
\(v_{eff}(\mathbf{r}) = v_{ext}(\mathbf{r}) + v_H(\mathbf{r}) + v_{xc}(\mathbf{r})\)
Das Hartree-Potential ist
\(v_H(\mathbf{r}) = \int\frac{n(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\,d\mathbf{r}'\)
und das Austausch-Korrelationspotential ergibt sich als funktionale Ableitung der Austausch-Korrelationsenergie:
\(v_{xc}(\mathbf{r}) = \frac{\delta E_{xc}[n]}{\delta n(\mathbf{r})}\)
Elektronendichte aus den Orbitalen
Nach Lösung der Kohn-Sham-Gleichungen wird die Elektronendichte aus den besetzten Orbitalen berechnet:
\(n(\mathbf{r}) = \sum_i f_i|\phi_i(\mathbf{r})|^2\)
Die Faktoren \(f_i\) geben die Besetzung der einzelnen Zustände an.
Hier entsteht die charakteristische Selbstkonsistenz der Methode: Die Kohn-Sham-Gleichungen benötigen das effektive Potential, dieses hängt von der Elektronendichte ab, und die Elektronendichte wird wiederum erst aus den Lösungen der Kohn-Sham-Gleichungen bestimmt.
Bedeutung der Kohn-Sham-Eigenwerte
Die Kohn-Sham-Eigenwerte sind zunächst mathematische Größen des nicht wechselwirkenden Referenzsystems. Sie dürfen nicht allgemein mit experimentell messbaren Anregungsenergien gleichgesetzt werden.
Eine wichtige Ausnahme betrifft unter idealen Bedingungen den höchsten besetzten Zustand des exakten Funktionals. Für die Ionisierungsenergie besteht die Beziehung
\(\epsilon_{HOMO} = -I\)
wobei \(I\) die Ionisierungsenergie bezeichnet.
Für Bandlücken und Anregungsspektren sind einfache Differenzen von Kohn-Sham-Eigenwerten dagegen im Allgemeinen nicht exakt.
Der selbstkonsistente Lösungsprozess
Self-Consistent Field-Verfahren
Die Kohn-Sham-Gleichungen sind nichtlinear, weil ihr Potential von der gesuchten Elektronendichte abhängt. In der Praxis werden sie deshalb iterativ gelöst.
Wahl einer Anfangsdichte
Zunächst wird eine geschätzte Elektronendichte \(n^{(0)}(\mathbf{r})\) vorgegeben. Sie kann beispielsweise aus überlagerten atomaren Dichten konstruiert werden.
Aufbau des effektiven Potentials
Aus der aktuellen Elektronendichte werden Hartree- und Austausch-Korrelationspotential berechnet. Zusammen mit dem äußeren Potential entsteht daraus \(v_{eff}(\mathbf{r})\).
Lösung der Kohn-Sham-Gleichungen
Für das aktuelle Potential werden die Kohn-Sham-Orbitale und ihre Eigenwerte bestimmt.
Berechnung einer neuen Dichte
Aus den besetzten Orbitalen entsteht eine neue Elektronendichte \(n_{out}(\mathbf{r})\).
Mischung und Wiederholung
Um numerische Instabilitäten zu vermeiden, wird die neue Dichte häufig mit der bisherigen Dichte gemischt:
\(n^{(k+1)} = (1-\alpha)n^{(k)} + \alpha n_{out}^{(k)}\)
Der Mischungsparameter \(\alpha\) beeinflusst die Stabilität und Geschwindigkeit der Konvergenz.
Konvergenz
Die Iteration wird fortgesetzt, bis sich Energie, Dichte oder Potential nur noch um einen vorgegebenen kleinen Betrag ändern. Ein mögliches Energiekriterium lautet
\(|E^{(k+1)}-E^{(k)}| < \epsilon\)
Die konkrete Wahl des Konvergenzkriteriums hängt von der später gewünschten Genauigkeit ab. Für Geometrieoptimierungen müssen beispielsweise auch die Kräfte auf die Atomkerne ausreichend genau konvergiert sein.
Austausch- und Korrelationsfunktionale
Warum Näherungen erforderlich sind
Der theoretische Rahmen der Kohn-Sham-DFT ist exakt, wenn das exakte Austausch-Korrelationsfunktional bekannt wäre. In der Praxis existiert jedoch keine allgemeine geschlossene Formel für \(E_{xc}[n]\). Deshalb entscheidet die Wahl des Funktionals wesentlich über die Qualität einer Berechnung.
Local Density Approximation
Die Local Density Approximation, kurz LDA, betrachtet die Elektronendichte lokal so, als gehöre sie zu einem homogenen Elektronengas. Die Austausch-Korrelationsenergie wird dabei näherungsweise als
\(E_{xc}^{LDA}[n] = \int n(\mathbf{r})\epsilon_{xc}(n(\mathbf{r}))\,d\mathbf{r}\)
geschrieben.
LDA ist besonders bei Systemen mit relativ langsam variierender Elektronendichte erfolgreich. Bei Festkörpern liefert sie häufig brauchbare Strukturen und Bindungseigenschaften. Typischerweise neigt sie jedoch dazu, Bindungen etwas zu stark vorherzusagen und Gleichgewichtsabstände etwas zu unterschätzen.
Generalized Gradient Approximation
Die Generalized Gradient Approximation berücksichtigt zusätzlich den räumlichen Gradienten der Elektronendichte:
\(E_{xc}^{GGA}[n] = \int f(n(\mathbf{r}),\nabla n(\mathbf{r}))\,d\mathbf{r}\)
Ein verbreitetes GGA-Funktional ist PBE. GGA-Verfahren liefern bei vielen Molekülen, Oberflächen und Festkörpern verbesserte Strukturparameter und Energien gegenüber einfachen LDA-Berechnungen.
Sie lösen jedoch nicht alle grundsätzlichen Probleme der Dichtefunktionaltheorie. Insbesondere Bandlücken, stark lokalisierte Elektronen und langreichweitige Dispersionswechselwirkungen können weiterhin unzureichend beschrieben werden.
Meta-GGA
Meta-GGA-Funktionale verwenden neben der Dichte und ihrem Gradienten zusätzliche lokale Größen. Häufig wird die kinetische Energiedichte
\(\tau(\mathbf{r}) = \frac{1}{2}\sum_i |\nabla\phi_i(\mathbf{r})|^2\)
einbezogen.
Dadurch können verschiedene elektronische Umgebungen besser unterschieden werden. Meta-GGA-Verfahren erreichen in vielen Fällen eine höhere Genauigkeit, ohne zwingend den hohen Rechenaufwand klassischer Hybridfunktionale zu verursachen.
Hybridfunktionale
Hybridfunktionale mischen einen Anteil exakten Hartree-Fock-Austauschs mit einem DFT-Austausch-Korrelationsfunktional. Schematisch kann dies als
\(E_{xc}^{hybrid} = aE_x^{HF} + (1-a)E_x^{DFT} + E_c^{DFT}\)
geschrieben werden.
Der Parameter \(a\) bestimmt den Anteil des exakten Austauschs. Bekannte Vertreter sind B3LYP, PBE0 und HSE.
Hybridfunktionale können insbesondere molekulare Energien und Bandstrukturen verbessern. Ihr Nachteil liegt im höheren numerischen Aufwand, da die nichtlokale Austauschwechselwirkung explizit ausgewertet werden muss.
Dispersionswechselwirkungen
Langreichweitige Van-der-Waals-Kräfte werden von einfachen lokalen oder semilokalen Funktionalen nicht zuverlässig erfasst. Das ist insbesondere bei Molekülkristallen, Adsorption, geschichteten Materialien und biologischen Systemen problematisch.
Zur Behandlung dieser Effekte werden zusätzliche Dispersionskorrekturen oder nichtlokale Funktionale verwendet. Die Wahl eines solchen Verfahrens muss an die physikalischen Eigenschaften des untersuchten Systems angepasst werden.
Numerische Umsetzung der Kohn-Sham-DFT
Darstellung der Orbitale
Um Kohn-Sham-Gleichungen auf einem Computer zu lösen, müssen die Orbitale in einer endlichen mathematischen Darstellung beschrieben werden. Häufig erfolgt dies durch Entwicklung in einer Basis:
\(\phi_i(\mathbf{r}) = \sum_{\mu} c_{\mu i}\chi_{\mu}(\mathbf{r})\)
Dabei sind \(\chi_{\mu}(\mathbf{r})\) Basisfunktionen und \(c_{\mu i}\) die zu bestimmenden Koeffizienten.
Ebene Wellen
Für periodische Festkörper werden besonders häufig ebene Wellen verwendet. Eine ebene Welle besitzt die Form
\(\exp(i(\mathbf{k}+\mathbf{G})\cdot\mathbf{r})\)
wobei \(\mathbf{k}\) ein Wellenvektor innerhalb der ersten Brillouin-Zone und \(\mathbf{G}\) ein reziproker Gittervektor ist.
Die Zahl der verwendeten ebenen Wellen wird über einen Energie-Cutoff begrenzt. Nur Funktionen mit
\(\frac{\hbar^2|\mathbf{k}+\mathbf{G}|^2}{2m_e} \leq E_{cut}\)
werden berücksichtigt.
Ein höherer Cutoff verbessert im Allgemeinen die numerische Genauigkeit, erhöht aber den Speicherbedarf und die Rechenzeit.
Lokalisierte Basisfunktionen
In molekularen Rechnungen werden häufig lokalisierte atomare Orbitale oder Gaußfunktionen verwendet. Solche Basisfunktionen können elektronische Zustände mit vergleichsweise wenigen Freiheitsgraden darstellen.
Die Wahl des Basissatzes beeinflusst allerdings ebenfalls das Ergebnis. Ein zu kleiner Basissatz kann systematische Fehler verursachen.
Pseudopotentiale
Die stark gebundenen inneren Elektronen eines Atoms verändern sich bei chemischen Bindungen meist nur wenig. Deshalb werden sie in vielen Rechnungen nicht explizit behandelt. Stattdessen wird ihre Wirkung zusammen mit dem Kern durch ein effektives Pseudopotential beschrieben.
Dadurch müssen nur die chemisch relevanten Valenzelektronen explizit berechnet werden. Der Rechenaufwand kann besonders bei schweren Elementen erheblich reduziert werden.
Projector-Augmented-Wave-Methode
Die Projector-Augmented-Wave-Methode verbindet Eigenschaften von Pseudopotentialverfahren mit einer Rekonstruktion der vollständigen Wellenfunktion in der Nähe der Atomkerne. Sie wird deshalb häufig für präzise Festkörper- und Materialrechnungen eingesetzt.
Periodische Systeme und k-Punkte
Für einen periodischen Kristall gilt nach dem Bloch-Theorem
\(\phi_{n\mathbf{k}}(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}}u_{n\mathbf{k}}(\mathbf{r})\)
wobei \(u_{n\mathbf{k}}(\mathbf{r})\) dieselbe Periodizität wie das Kristallgitter besitzt.
Integrale über die Brillouin-Zone werden numerisch durch eine endliche Menge von k-Punkten angenähert. Für Metalle sind häufig dichtere k-Punkt-Gitter erforderlich als für Isolatoren mit großer Bandlücke.
Spin, Magnetismus und relativistische Effekte
Spinpolarisierte Dichtefunktionaltheorie
Bei magnetischen Systemen müssen Elektronen unterschiedlicher Spinorientierung getrennt behandelt werden. Die Gesamtdichte wird dann als
\(n(\mathbf{r}) = n_{\uparrow}(\mathbf{r}) + n_{\downarrow}(\mathbf{r})\)
geschrieben.
Die lokale Magnetisierung kann entsprechend durch
\(m(\mathbf{r}) = n_{\uparrow}(\mathbf{r}) - n_{\downarrow}(\mathbf{r})\)
charakterisiert werden.
Damit können ferromagnetische und antiferromagnetische Zustände sowie lokale magnetische Momente untersucht werden.
Nichtkollinearer Magnetismus
In komplexen magnetischen Materialien genügt eine feste globale Quantisierungsachse nicht immer. Bei nichtkollinearem Magnetismus kann sich die lokale Richtung der Magnetisierung räumlich verändern. Die elektronische Struktur muss dann durch Spinoren beschrieben werden.
Spin-Bahn-Kopplung
Die Spin-Bahn-Kopplung verbindet den Spin eines Elektrons mit seiner orbitalen Bewegung. Sie wird insbesondere bei schweren Elementen und in Materialien mit starken relativistischen Effekten wichtig.
Eine vereinfachte Darstellung eines Spin-Bahn-Terms lautet
\(\hat{H}_{SO} \propto \mathbf{L}\cdot\mathbf{S}\)
mit Bahndrehimpuls \(\mathbf{L}\) und Spin \(\mathbf{S}\).
Spin-Bahn-Kopplung ist entscheidend für zahlreiche moderne Forschungsgebiete, darunter topologische Materialien, Spintronik und bestimmte Festkörper-Qubit-Konzepte.
Berechenbare physikalische Eigenschaften
Gesamtenergie und Struktur
Eine der wichtigsten Anwendungen der Kohn-Sham-DFT ist die Berechnung relativer Gesamtenergien. Unterschiedliche atomare Strukturen können miteinander verglichen werden, um stabile oder metastabile Konfigurationen zu bestimmen.
Die auf einen Kern wirkende Kraft ergibt sich grundsätzlich aus der Ableitung der Gesamtenergie:
\(\mathbf{F}_I = -\frac{\partial E}{\partial\mathbf{R}_I}\)
Durch iterative Verschiebung der Atompositionen können Gleichgewichtsgeometrien und Kristallstrukturen berechnet werden.
Zustandsdichte
Die elektronische Zustandsdichte gibt an, wie viele elektronische Zustände innerhalb eines bestimmten Energieintervalls vorhanden sind. Formal kann sie als
\(D(E) = \sum_i \delta(E-\epsilon_i)\)
geschrieben werden.
In realen numerischen Berechnungen wird die Deltafunktion meist durch eine endliche Verbreiterungsfunktion ersetzt.
Bandstruktur
In Kristallen hängen die elektronischen Energien vom Wellenvektor \(\mathbf{k}\) ab. Die Funktionen
\(\epsilon_n(\mathbf{k})\)
bilden die elektronische Bandstruktur.
Aus ihr lässt sich erkennen, ob ein ideales Material innerhalb der verwendeten Näherung metallisch, halbleitend oder isolierend erscheint. Auch effektive Massen, Talstrukturen und bestimmte Transportmerkmale können daraus abgeleitet werden.
Phononen und Gitterschwingungen
Durch Berechnung von Kräften und Kraftkonstanten können Gitterschwingungen bestimmt werden. Im harmonischen Grenzfall führt dies auf ein Eigenwertproblem der Form
\(\sum_{J\beta}D_{I\alpha,J\beta}(\mathbf{q})e_{J\beta}^{(\nu)}(\mathbf{q}) = \omega_{\nu}^2(\mathbf{q})e_{I\alpha}^{(\nu)}(\mathbf{q})\)
Dabei bezeichnet \(D\) die dynamische Matrix und \(\omega_{\nu}(\mathbf{q})\) die Frequenz eines Phononenmodus.
Phononenrechnungen liefern Informationen über dynamische Stabilität, Wärmetransport und Elektron-Phonon-Wechselwirkungen.
Kohn-Sham-DFT in der Material- und Quantentechnologie
Materialdesign auf atomarer Ebene
Kohn-Sham-DFT ermöglicht die Untersuchung eines Materials noch bevor es experimentell hergestellt wurde. Atomare Strukturen können virtuell konstruiert, optimiert und auf ihre elektronischen Eigenschaften untersucht werden.
Dadurch lassen sich Kandidaten für Halbleiter, magnetische Materialien, Elektroden, Katalysatoren und quantentechnologisch relevante Festkörper systematisch vergleichen.
Defekte und Dotierungen
Punktdefekte spielen in Halbleitern und Quantenmaterialien eine zentrale Rolle. Ein Defekt kann zusätzliche elektronische Zustände erzeugen, Ladung binden oder magnetische Momente verursachen.
Die Bildungsenergie eines Defekts wird häufig schematisch als
\(E_f = E_{def} - E_{bulk} - \sum_i n_i\mu_i + q(E_F + E_{VBM}) + E_{corr}\)
geschrieben.
Dabei ist \(E_{def}\) die Energie des defekthaltigen Systems, \(E_{bulk}\) die Energie des idealen Kristalls, \(\mu_i\) sind chemische Potentiale, \(q\) ist der Ladungszustand und \(E_{corr}\) bezeichnet mögliche Korrekturterme.
Solche Rechnungen sind unter anderem für Farbzentren, Einzelphotonenquellen und Festkörper-Qubits relevant.
Zweidimensionale Materialien
Graphen, Übergangsmetalldichalkogenide und andere zweidimensionale Materialien besitzen elektronische Eigenschaften, die stark von Schichtzahl, Defekten, mechanischer Spannung und Substratwechselwirkungen abhängen.
Kohn-Sham-DFT ermöglicht die systematische Untersuchung dieser Faktoren. Durch Variation der atomaren Struktur lassen sich beispielsweise Veränderungen von Bandlücken, magnetischen Eigenschaften und Spin-Bahn-Effekten berechnen.
Materialien für Quantentechnologien
Für Quantencomputer, Quantensensoren und Quantenkommunikation werden Materialien benötigt, deren mikroskopische Eigenschaften äußerst präzise kontrolliert werden müssen. DFT kann hierbei Informationen über Defektzustände, Ladungszustände, lokale Geometrien, magnetische Wechselwirkungen und Spin-Bahn-Kopplung liefern.
Die Methode simuliert jedoch nicht automatisch die vollständige kohärente Dynamik eines Qubits. Sie liefert vielmehr die elektronische und energetische Grundlage, auf der weiterführende Modelle aufgebaut werden können.
Grenzen der Kohn-Sham-DFT
Abhängigkeit vom Funktional
Die größte praktische Unsicherheit der Kohn-Sham-DFT liegt in der unbekannten exakten Form des Austausch-Korrelationsfunktionals. Unterschiedliche Näherungen können für dasselbe System unterschiedliche Resultate liefern.
Es gibt deshalb kein einzelnes Standardfunktional, das für sämtliche Materialien und physikalischen Fragestellungen gleichermaßen optimal ist.
Selbstwechselwirkungsfehler
Im Hartree-Term wechselwirkt die Elektronendichte formal auch mit sich selbst. Im exakten Funktional wird diese unphysikalische Selbstwechselwirkung vollständig durch den Austauschanteil kompensiert. Näherungsfunktionale erreichen diese Aufhebung jedoch oft nur unvollständig.
Das kann zu einer zu starken Delokalisierung elektronischer Zustände und zu Fehlern bei Ladungstransferprozessen führen.
Bandlückenproblem
Semilokale DFT-Funktionale unterschätzen bei vielen Halbleitern und Isolatoren die fundamentale Bandlücke. Die fundamentale Lücke ist definiert als
\(E_g = I - A\)
wobei \(I\) die Ionisierungsenergie und \(A\) die Elektronenaffinität bezeichnet.
Die Differenz zwischen niedrigstem unbesetzten und höchstem besetzten Kohn-Sham-Eigenwert
\(\epsilon_{LUMO} - \epsilon_{HOMO}\)
ist im Allgemeinen nicht identisch mit dieser fundamentalen Lücke.
Für das exakte Funktional tritt zusätzlich eine sogenannte Ableitungsdiskontinuität auf. Formal lässt sich schreiben
\(E_g = \epsilon_{LUMO} - \epsilon_{HOMO} + \Delta_{xc}\)
wobei \(\Delta_{xc}\) die Austausch-Korrelations-Diskontinuität bezeichnet.
Stark korrelierte Systeme
Bei Materialien mit lokalisierten d- oder f-Elektronen kann die unabhängige Kohn-Sham-Teilchenbeschreibung qualitativ unzureichend werden. Beispiele sind zahlreiche Übergangsmetalloxide, Mott-Isolatoren und Materialien mit schweren Fermionen.
In solchen Fällen reichen gewöhnliche LDA- oder GGA-Verfahren häufig nicht aus. Zusätzliche Korrekturen oder weiterführende Vielteilchenmethoden werden benötigt.
Angeregte Zustände
Die ursprüngliche Hohenberg-Kohn- und Kohn-Sham-Theorie ist eine Grundzustandstheorie. Optische Anregungen, Exzitonen und zeitabhängige Prozesse sind deshalb nicht unmittelbar Gegenstand der statischen Grundzustands-DFT.
Für diese Fragestellungen werden beispielsweise zeitabhängige DFT oder Vielteilchenmethoden eingesetzt.
Erweiterungen und weiterführende Verfahren
DFT+U
DFT+U ergänzt ein gewöhnliches DFT-Funktional um einen zusätzlichen Term für stark lokalisierte Elektronenzustände. In vereinfachter Form kann eine Korrektur als
\(E_{DFT+U} = E_{DFT} + E_U\)
geschrieben werden.
Eine häufig verwendete vereinfachte Form lautet
\(E_U = \frac{U_{eff}}{2}\sum_m n_m(1-n_m)\)
mit einer effektiven Wechselwirkungsstärke \(U_{eff}\).
DFT+U kann die Lokalisierung von d- und f-Elektronen verbessern, führt jedoch einen zusätzlichen Parameter ein, dessen Wahl begründet werden muss.
GW-Näherung
Die GW-Methode ist eine Vielteilchenmethode zur Berechnung von Quasiteilchenenergien. Ihre Selbstenergie wird näherungsweise als Produkt der Einteilchen-Green-Funktion \(G\) und der abgeschirmten Coulomb-Wechselwirkung \(W\) behandelt:
\(\Sigma \approx iGW\)
GW-Rechnungen können deutlich bessere Bandlücken liefern als einfache semilokale Kohn-Sham-DFT. Sie sind allerdings wesentlich rechenintensiver.
Bethe-Salpeter-Gleichung
Für optische Spektren und Exzitonen muss die Wechselwirkung zwischen angeregtem Elektron und zurückbleibendem Loch berücksichtigt werden. Die Bethe-Salpeter-Gleichung baut typischerweise auf Quasiteilchenenergien auf und beschreibt solche gebundenen Elektron-Loch-Zustände.
Zeitabhängige Dichtefunktionaltheorie
Die zeitabhängige Dichtefunktionaltheorie erweitert das Dichtekonzept auf zeitabhängige Prozesse. Die zeitabhängigen Kohn-Sham-Gleichungen besitzen die Form
\(i\hbar\frac{\partial}{\partial t}\phi_i(\mathbf{r},t) = \left[-\frac{\hbar^2}{2m_e}\nabla^2 + v_{eff}(\mathbf{r},t)\right]\phi_i(\mathbf{r},t)\)
TDDFT wird insbesondere zur Untersuchung optischer Anregungen, elektronischer Dynamik und zeitabhängiger Felder eingesetzt.
DFT und Dynamical Mean-Field Theory
Die Kombination aus DFT und Dynamical Mean-Field Theory wird für stark korrelierte Materialien eingesetzt. DFT beschreibt dabei die elektronische Struktur des Materials, während DMFT lokale zeitabhängige Vielteilchenkorrelationen expliziter behandelt.
Diese Kombination ist vor allem bei Systemen interessant, in denen gewöhnliche Bandtheorien allein nicht ausreichen.
Machine Learning und neue Funktionale
Maschinelles Lernen wird zunehmend verwendet, um Energiefunktionale, Potentiale und elektronische Strukturgrößen aus umfangreichen Referenzdaten zu approximieren. Ziel ist es, Genauigkeiten auf dem Niveau aufwendiger Methoden bei deutlich geringerem Rechenaufwand zu erreichen.
Ein maschinell gelerntes Funktional kann schematisch als
\(E_{xc}^{ML} = F_{\theta}[n]\)
geschrieben werden, wobei \(F_{\theta}\) ein durch Parameter \(\theta\) beschriebenes Modell ist.
Eine zentrale Herausforderung bleibt dabei die physikalische Übertragbarkeit auf Systeme außerhalb der Trainingsdaten.
Rechenaufwand und Skalierung
Numerische Kosten
Die klassische Kohn-Sham-DFT ist erheblich effizienter als viele explizite Vielteilchenverfahren. Dennoch steigt ihr Rechenaufwand mit zunehmender Systemgröße stark an.
Bei traditionellen Implementierungen kann die Lösung des Eigenwertproblems näherungsweise mit
\(O(N^3)\)
skalieren, wobei \(N\) eine Größe proportional zur Zahl der Elektronen oder Basisfunktionen bezeichnet.
Diese Skalierung ist keine unvermeidbare Eigenschaft sämtlicher DFT-Algorithmen, beschreibt aber den Aufwand vieler Standardverfahren.
Parallelisierung
Moderne DFT-Programme verteilen Berechnungen auf viele Prozessorkerne und Rechenknoten. Parallelisiert werden können beispielsweise k-Punkte, Bänder, schnelle Fourier-Transformationen und bestimmte Matrixoperationen.
Dadurch lassen sich Systeme mit Hunderten oder in spezialisierten Verfahren auch Tausenden Atomen untersuchen. Für sehr große Systeme werden zusätzlich linear skalierende Methoden entwickelt.
Genauigkeit und Kosten
Die Auswahl numerischer Parameter ist immer ein Kompromiss. Ein größerer Basissatz, ein höherer Energie-Cutoff, ein dichteres k-Punkt-Gitter oder ein aufwendigeres Funktional kann die Genauigkeit verbessern, erhöht jedoch die Rechenzeit.
Seriöse DFT-Berechnungen erfordern deshalb Konvergenztests. Ein Ergebnis ist nur dann belastbar, wenn gezeigt werden kann, dass es gegenüber den relevanten numerischen Parametern ausreichend stabil ist.
Bedeutung der Kohn-Sham-DFT für die moderne Quantenforschung
Verbindung zwischen Theorie und Experiment
Die Kohn-Sham-DFT besitzt eine besondere Stellung, weil sie atomare Modelle direkt mit experimentell relevanten Materialeigenschaften verbindet. Berechnete Strukturen, Bindungsenergien, Schwingungen, Zustandsdichten und Bandstrukturen können mit Messdaten verglichen werden.
Umgekehrt können theoretische Rechnungen experimentelle Beobachtungen erklären, indem sie mikroskopische Mechanismen sichtbar machen, die experimentell nur indirekt zugänglich sind.
High-Throughput-Berechnungen
Da Kohn-Sham-DFT im Vergleich zu hochgenauen Vielteilchenmethoden relativ effizient ist, können große Mengen verschiedener Materialien automatisch untersucht werden.
Bei einem High-Throughput-Ansatz werden für Tausende oder sogar mehr hypothetische und bekannte Kristallstrukturen standardisierte Berechnungen durchgeführt. Anschließend können Kandidaten anhand vorgegebener physikalischer Eigenschaften gefiltert werden.
Damit hat sich DFT von einer Methode zur Analyse einzelner Materialien zu einem Werkzeug für systematische Materialsuchen entwickelt.
DFT und künstliche Intelligenz
DFT-Daten bilden außerdem eine wichtige Grundlage für maschinelles Lernen in der Materialforschung. Modelle können mit berechneten Energien, Kräften und elektronischen Eigenschaften trainiert werden und anschließend neue Strukturen wesentlich schneller vorhersagen.
Dadurch entsteht eine Arbeitsteilung: DFT liefert physikalisch fundierte Referenzdaten, während Machine-Learning-Modelle große Konfigurationsräume effizient durchsuchen können.
Rolle in der Quantentechnologie
Die Entwicklung quantentechnologischer Bauelemente verlangt präzise Kenntnisse über atomare Defekte, elektronische Zustände, magnetische Wechselwirkungen und relativistische Effekte. Genau auf dieser mikroskopischen Ebene liefert die Kohn-Sham-DFT besonders wertvolle Informationen.
Sie ist deshalb kein Quantencomputer-Algorithmus und auch keine vollständige Theorie kohärenter Qubit-Dynamik. Ihre Rolle besteht vielmehr darin, die materielle Grundlage vieler Quantentechnologien quantitativ zu untersuchen.
Fazit
Die Kohn-Sham-Dichtefunktionaltheorie löst eines der zentralen praktischen Probleme der Quantenmechanik: die elektronische Struktur komplexer Vielteilchensysteme mit vertretbarem numerischem Aufwand zu bestimmen.
Ihr entscheidender Schritt besteht darin, die vollständige Vielteilchen-Wellenfunktion durch die wesentlich einfachere Elektronendichte als fundamentale Variable zu ersetzen. Die Hohenberg-Kohn-Theoreme zeigen, dass die Grundzustandsdichte prinzipiell sämtliche Informationen über den Grundzustand enthält.
Der Kohn-Sham-Ansatz macht diese Idee praktisch nutzbar. Das reale System wechselwirkender Elektronen wird durch ein nicht wechselwirkendes Referenzsystem mit identischer Grundzustandsdichte ersetzt. Die Elektronen bewegen sich dabei in einem effektiven Potential. Die resultierenden Kohn-Sham-Gleichungen besitzen die Form von Einteilchengleichungen und können numerisch effizient gelöst werden.
Die zentrale Unsicherheit liegt im Austausch-Korrelationsfunktional. Für dieses Funktional existiert keine allgemeine exakte geschlossene Form. LDA, GGA, Meta-GGA, Hybridfunktionale und weitere Ansätze stellen unterschiedliche Näherungsstufen dar. Ihre Qualität hängt von der untersuchten Stoffklasse und der gewünschten physikalischen Größe ab.
Damit ist Kohn-Sham-DFT weder eine einfache Näherungsmethode noch eine universell exakte Lösung. Der theoretische Rahmen ist grundsätzlich exakt, die praktische Umsetzung hängt jedoch entscheidend von näherungsweisen Austausch-Korrelationsfunktionalen und numerischen Parametern ab.
Gerade das Verhältnis von Aussagekraft und Rechenaufwand erklärt ihren außergewöhnlichen Erfolg. Mit Kohn-Sham-DFT lassen sich Kristallstrukturen, molekulare Geometrien, Bindungsenergien, Defekte, Magnetismus, Bandstrukturen, Phononen und zahlreiche weitere Eigenschaften quantitativ untersuchen.
Für stark korrelierte Materialien, hochpräzise Quasiteilchenenergien oder komplexe angeregte Zustände reichen gewöhnliche Kohn-Sham-Verfahren nicht immer aus. Methoden wie DFT+U, Hybrid-DFT, GW, TDDFT oder DFT+DMFT erweitern den Ansatz für spezielle Fragestellungen.
In der modernen Material- und Quantenforschung bleibt die Kohn-Sham-DFT damit eines der zentralen Werkzeuge zur Verbindung fundamentaler Quantenmechanik mit realen Materialien. Sie ermöglicht nicht nur die Erklärung bekannter Stoffeigenschaften, sondern auch die computergestützte Suche nach neuen Materialien und bildet eine wichtige Grundlage für moderne Hochdurchsatzrechnungen und datengetriebene Materialentwicklung.
Mit freundlichen Grüßen
Anhang
Wissenschaftliche Zeitschriften und Artikel
Die wissenschaftliche Primärliteratur zur Kohn-Sham-Dichtefunktionaltheorie umfasst sowohl die fundamentalen Arbeiten zur Dichtefunktionaltheorie als auch Veröffentlichungen zu Austausch-Korrelationsfunktionalen, Bandlücken, elektronischer Korrelation und Erweiterungen der ursprünglichen Kohn-Sham-Methode. Für eine wissenschaftliche Abhandlung sollten die Originalarbeiten von Hohenberg, Kohn und Sham den theoretischen Ausgangspunkt bilden. Neuere Arbeiten dienen dazu, die praktischen Näherungen, ihre Leistungsfähigkeit und ihre bekannten Grenzen fachlich einzuordnen.
Grundlegende Primärliteratur zur Dichtefunktionaltheorie
- P. Hohenberg, W. Kohn: Inhomogeneous Electron Gas, Physical Review 136, B864–B871, 1964.
- Diese Arbeit bildet das theoretische Fundament der modernen Grundzustands-Dichtefunktionaltheorie. Die beiden Hohenberg-Kohn-Theoreme zeigen, dass das externe Potential und damit die Grundzustandseigenschaften eines wechselwirkenden Elektronensystems grundsätzlich durch seine Grundzustandsdichte bestimmt werden. Für eine Abhandlung zur Kohn-Sham-DFT ist diese Veröffentlichung die wichtigste Primärquelle für die theoretische Rechtfertigung des Dichteansatzes.
- W. Kohn, L. J. Sham: Self-Consistent Equations Including Exchange and Correlation Effects, Physical Review 140, A1133–A1138, 1965.
- Die Originalarbeit von Kohn und Sham führt die heute zentrale praktische Form der Dichtefunktionaltheorie ein. Das wechselwirkende Elektronensystem wird auf ein Hilfssystem nicht wechselwirkender Elektronen mit derselben Grundzustandsdichte abgebildet. Die Quelle sollte insbesondere für die Einführung der Kohn-Sham-Orbitale, des effektiven Potentials, des selbstkonsistenten Lösungsverfahrens und des Austausch-Korrelationsanteils herangezogen werden.
Grundlegende Arbeiten zu Austausch-Korrelationsfunktionalen
- J. P. Perdew, K. Burke, M. Ernzerhof: Generalized Gradient Approximation Made Simple, Physical Review Letters 77, 3865–3868, 1996.
- Diese Veröffentlichung führt das PBE-Funktional ein, eine der einflussreichsten Generalized-Gradient-Approximationen der DFT. Die Arbeit eignet sich zur Erklärung des Übergangs von lokalen Näherungen zu gradientenabhängigen Funktionalen und ist besonders relevant für Festkörper-, Oberflächen- und Materialberechnungen.
- J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, K. Burke: Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces, Physical Review Letters 100, 136406, 2008.
- Die Arbeit führt PBEsol als speziell für dicht gepackte Festkörper und Oberflächen verbesserte Variante von PBE ein. Sie ist für eine quantentechnologisch ausgerichtete Abhandlung besonders nützlich, weil sie verdeutlicht, dass die Wahl des Austausch-Korrelationsfunktionals von der untersuchten Stoffklasse und der gewünschten physikalischen Größe abhängt.
- URL: https://journals.aps.org/...
- arXiv: https://arxiv.org/...
- DOI: https://doi.org/...
- Die Arbeit führt PBEsol als speziell für dicht gepackte Festkörper und Oberflächen verbesserte Variante von PBE ein. Sie ist für eine quantentechnologisch ausgerichtete Abhandlung besonders nützlich, weil sie verdeutlicht, dass die Wahl des Austausch-Korrelationsfunktionals von der untersuchten Stoffklasse und der gewünschten physikalischen Größe abhängt.
- J. Sun, A. Ruzsinszky, J. P. Perdew: Strongly Constrained and Appropriately Normed Semilocal Density Functional, Physical Review Letters 115, 036402, 2015.
- Mit dieser Arbeit wurde das SCAN-Funktional als moderne Meta-GGA eingeführt. Die Veröffentlichung eignet sich, um die nächste Entwicklungsstufe oberhalb einfacher GGA-Funktionale zu erläutern und zu zeigen, wie zusätzliche Informationen über die elektronische Umgebung in semilokale Funktionale eingebaut werden können.
- J. Tao, J. P. Perdew, V. N. Staroverov, G. E. Scuseria: Climbing the Density Functional Ladder: Nonempirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids, Physical Review Letters 91, 146401, 2003.
- Die Arbeit ist eine wichtige Quelle für das Verständnis von Meta-GGA-Funktionalen und der sogenannten Jacob's Ladder der Dichtefunktionalnäherungen. Sie eignet sich insbesondere dazu, verschiedene Funktionalklassen systematisch hinsichtlich ihrer verwendeten Informationen und ihrer theoretischen Komplexität einzuordnen.
Hybridfunktionale und exakter Austausch
- A. D. Becke: Density-functional thermochemistry. III. The role of exact exchange, The Journal of Chemical Physics 98, 5648–5652, 1993.
- Diese vielzitierte Arbeit ist grundlegend für das Verständnis moderner Hybridfunktionale. Sie zeigt, wie ein Anteil exakten Hartree-Fock-Austauschs mit DFT-Komponenten kombiniert werden kann. Die Quelle ist insbesondere für Abschnitte zu Hybrid-DFT, zur Verbesserung von Energien und zu den Grenzen rein semilokaler Funktionale relevant.
- DOI: https://doi.org/...
- Diese vielzitierte Arbeit ist grundlegend für das Verständnis moderner Hybridfunktionale. Sie zeigt, wie ein Anteil exakten Hartree-Fock-Austauschs mit DFT-Komponenten kombiniert werden kann. Die Quelle ist insbesondere für Abschnitte zu Hybrid-DFT, zur Verbesserung von Energien und zu den Grenzen rein semilokaler Funktionale relevant.
- J. P. Perdew, M. Ernzerhof, K. Burke: Rationale for Mixing Exact Exchange with Density Functional Approximations, The Journal of Chemical Physics 105, 9982–9985, 1996.
- Die Veröffentlichung liefert eine theoretische Begründung für die Mischung von exaktem Austausch mit Dichtefunktionalnäherungen. Sie ist besonders geeignet, um die physikalische Motivation von Hybridfunktionalen zu vertiefen und den Übergang von empirischen Mischungen zu systematischer begründeten Ansätzen darzustellen.
- DOI: https://doi.org/...
- Die Veröffentlichung liefert eine theoretische Begründung für die Mischung von exaktem Austausch mit Dichtefunktionalnäherungen. Sie ist besonders geeignet, um die physikalische Motivation von Hybridfunktionalen zu vertiefen und den Übergang von empirischen Mischungen zu systematischer begründeten Ansätzen darzustellen.
Grenzen der Kohn-Sham-DFT und Bandlückenproblem
- J. P. Perdew, R. G. Parr, M. Levy, J. L. Balduz Jr.: Density-Functional Theory for Fractional Particle Number: Derivative Discontinuities of the Energy, Physical Review Letters 49, 1691–1694, 1982.
- Diese grundlegende Arbeit behandelt die Abhängigkeit der Gesamtenergie von nichtganzzahligen Elektronenzahlen und die daraus entstehende Ableitungsdiskontinuität. Sie ist eine zentrale Quelle zur fachlich korrekten Erklärung, warum die fundamentale Bandlücke nicht allgemein mit der Differenz zwischen unbesetztem und besetztem Kohn-Sham-Zustand gleichgesetzt werden darf.
Spezialisierte Arbeiten zu stark korrelierten Systemen
- S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, A. P. Sutton: Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study, Physical Review B 57, 1505–1509, 1998.
- Diese Arbeit ist eine zentrale Referenz für die heute weit verbreitete vereinfachte DFT+U-Formulierung. Sie zeigt anhand von Nickeloxid, weshalb gewöhnliche lokale DFT-Näherungen bei lokalisierten d-Elektronen versagen können und wie eine zusätzliche lokale Wechselwirkungskorrektur die Beschreibung verbessert. Für den Abschnitt über stark korrelierte Materialien und DFT+U ist diese Quelle besonders relevant.
Historische und fachliche Einordnung der Kohn-Sham-DFT
- A. D. Becke: Perspective: Fifty years of density-functional theory in chemical physics, The Journal of Chemical Physics 140, 18A301, 2014.
- Diese Übersichtsarbeit ordnet fünf Jahrzehnte Kohn-Sham-DFT historisch und fachlich ein. Sie behandelt die Entwicklung des Kohn-Sham-Konzepts, Fortschritte bei Austausch-Korrelationsfunktionalen sowie zentrale Herausforderungen der Methode. Sie eignet sich besonders für Einleitung, wissenschaftliche Einordnung und Diskussion der langfristigen Entwicklung der DFT.
- DOI: https://doi.org/...
- Diese Übersichtsarbeit ordnet fünf Jahrzehnte Kohn-Sham-DFT historisch und fachlich ein. Sie behandelt die Entwicklung des Kohn-Sham-Konzepts, Fortschritte bei Austausch-Korrelationsfunktionalen sowie zentrale Herausforderungen der Methode. Sie eignet sich besonders für Einleitung, wissenschaftliche Einordnung und Diskussion der langfristigen Entwicklung der DFT.
Bücher und Monographien
Monographien sind für die Kohn-Sham-DFT besonders wichtig, weil die mathematische Herleitung, die Interpretation der Funktionale und die numerische Umsetzung in einzelnen Fachartikeln meist nur ausschnittsweise behandelt werden. Die folgenden Werke ergänzen die Primärliteratur durch systematische Darstellungen von Grundlagentheorie, elektronischer Strukturrechnung und praktischer Anwendung.
Grundlagenwerke zur Dichtefunktionaltheorie
- R. G. Parr, W. Yang: Density-Functional Theory of Atoms and Molecules, Oxford University Press, Ausgabe 1995.
- Das Werk bietet eine mathematisch fundierte Darstellung der Dichtefunktionaltheorie mit besonderem Schwerpunkt auf Atomen und Molekülen. Es eignet sich zur Vertiefung der Hohenberg-Kohn-Theoreme, der Variationsprinzipien, chemischer Potentiale und grundlegender Dichtefunktionalkonzepte. Für theoretische Definitionen und konzeptionelle Hintergründe ist es eine wichtige Referenz.
- E. Engel, R. M. Dreizler: Density Functional Theory: An Advanced Course, Springer, 2011.
- Diese Monographie bietet eine besonders gründliche theoretische Behandlung der DFT. Sie umfasst Existenztheoreme, effektive Einteilchengleichungen, Austausch-Korrelationsfunktionale, orbitalabhängige Funktionale sowie zeitabhängige und relativistische DFT. Das Werk eignet sich für vertiefende theoretische Argumentationen und für Abschnitte, die über eine rein anwendungsorientierte Darstellung hinausgehen.
Standardwerke zur elektronischen Struktur von Materialien
- R. M. Martin: Electronic Structure: Basic Theory and Practical Methods, Cambridge University Press, 2004.
- Martin verbindet die Theorie der elektronischen Struktur mit ihrer praktischen numerischen Umsetzung. Das Buch behandelt unter anderem Elektronengasmodelle, Dichtefunktionaltheorie, periodische Festkörper, ebene Wellen, Pseudopotentiale und numerische Lösungsverfahren. Für eine Abhandlung mit Schwerpunkt auf Quantentechnologie und Materialwissenschaft ist es besonders wertvoll, da es die Brücke zwischen theoretischer Kohn-Sham-DFT und realen Festkörperberechnungen schlägt.
- DOI: https://doi.org/...
- Martin verbindet die Theorie der elektronischen Struktur mit ihrer praktischen numerischen Umsetzung. Das Buch behandelt unter anderem Elektronengasmodelle, Dichtefunktionaltheorie, periodische Festkörper, ebene Wellen, Pseudopotentiale und numerische Lösungsverfahren. Für eine Abhandlung mit Schwerpunkt auf Quantentechnologie und Materialwissenschaft ist es besonders wertvoll, da es die Brücke zwischen theoretischer Kohn-Sham-DFT und realen Festkörperberechnungen schlägt.
Praxisorientierte Literatur für DFT-Berechnungen
- D. S. Sholl, J. A. Steckel: Density Functional Theory: A Practical Introduction, 2. Auflage, Wiley, 2022.
- Dieses Lehrbuch konzentriert sich auf die praktische Durchführung und Bewertung von DFT-Berechnungen. Behandelt werden unter anderem ebene Wellen, Festkörper, Oberflächen, magnetische Systeme, Konvergenztests und die Interpretation berechneter Eigenschaften. Es eignet sich hervorragend zur Absicherung des methodischen und numerischen Teils einer wissenschaftlichen Abhandlung.
- W. Koch, M. C. Holthausen: A Chemist's Guide to Density Functional Theory, 2. Auflage, Wiley-VCH, 2001.
- Das Werk bietet eine anwendungsorientierte Darstellung der DFT mit ausführlicher Behandlung von Austausch-Korrelationsfunktionalen, Selbstwechselwirkung, Hybridfunktionalen und praktischer Methodenwahl. Es eignet sich besonders zur verständlichen Erklärung unterschiedlicher Funktionalklassen und ihrer jeweiligen Stärken und Schwächen.
Online-Ressourcen und Datenbanken
Online-Ressourcen ergänzen die theoretische Literatur insbesondere bei praktischen Kohn-Sham-DFT-Rechnungen, beim Vergleich berechneter Materialeigenschaften und bei der Auswahl von Austausch-Korrelationsfunktionalen. Für wissenschaftliche Arbeiten sollten solche Plattformen primär als Recherche-, Vergleichs- und Dokumentationswerkzeuge verwendet werden. Theoretische Kernaussagen sollten nach Möglichkeit auf Primärliteratur oder etablierte Monographien zurückgeführt werden.
Materialdatenbanken und Forschungsplattformen
- Materials Project: Datenbank für berechnete Materialeigenschaften und elektronische Strukturen.
- Das Materials Project stellt umfangreiche computergestützte Materialdaten bereit, darunter Kristallstrukturen, thermodynamische Größen, elektronische Eigenschaften, Bandstrukturen und Zustandsdichten. Die Plattform eignet sich zum Vergleich eigener DFT-Ergebnisse mit vorhandenen Berechnungen und zur Identifikation interessanter Materialien für quantentechnologische Anwendungen.
- NOMAD: Forschungsdatenplattform für Materialwissenschaft und elektronische Strukturrechnungen.
- NOMAD sammelt, strukturiert und erschließt Forschungsdaten aus der computergestützten Materialwissenschaft. Die Plattform ist besonders nützlich für die Suche nach vorhandenen elektronischen Strukturrechnungen, für Datenvergleich und für Fragen der Nachvollziehbarkeit und Wiederverwendbarkeit wissenschaftlicher Simulationsdaten.
- OQMD: Open Quantum Materials Database.
- Die Open Quantum Materials Database enthält eine große Zahl DFT-berechneter thermodynamischer und struktureller Materialeigenschaften. Sie eignet sich für Materialscreening, Stabilitätsvergleiche und die Untersuchung chemischer Zusammensetzungsräume. Für eine Abhandlung kann sie als Beispiel für den Einsatz der Kohn-Sham-DFT in modernen High-Throughput-Materialstudien dienen.
- AFLOW: Automatic FLOW for Materials Discovery.
- AFLOW kombiniert eine umfangreiche Materialdatenbank mit automatisierten Verfahren für elektronische Struktur- und Materialberechnungen. Die Plattform verdeutlicht die Bedeutung standardisierter Kohn-Sham-DFT-Verfahren für High-Throughput-Screening, Phasenstabilität, Bandstrukturen und datengetriebene Materialforschung.
Dokumentationen für praktische Kohn-Sham-DFT-Rechnungen
- Quantum ESPRESSO: Offizielle Dokumentation für elektronische Struktur- und Materialberechnungen.
- Quantum ESPRESSO ist eine etablierte Open-Source-Programmsammlung für DFT-Berechnungen mit ebenen Wellen und Pseudopotentialen. Die Dokumentation ist eine wertvolle technische Referenz für Self-Consistent-Field-Rechnungen, k-Punkt-Abtastungen, Phononen, Bandstrukturen, Pseudopotentiale und zahlreiche weiterführende Methoden. Sie eignet sich für die methodische Konkretisierung praktischer DFT-Verfahren.
- VASP Wiki: Theoretische und praktische Dokumentation zur elektronischen Strukturrechnung.
- Die VASP-Dokumentation behandelt zentrale Aspekte moderner Festkörper-DFT, darunter Kohn-Sham-Gleichungen, Self-Consistent-Field-Verfahren, Austausch-Korrelationsfunktionale, DFT+U, Hybridfunktionale, Spin-Bahn-Kopplung, Bandstrukturen und elektronische Zustandsdichten. Sie ist besonders nützlich, um den Zusammenhang zwischen theoretischen Methoden und ihrer numerischen Implementierung nachvollziehbar darzustellen.
- URL: https://vasp.at/...
- Die VASP-Dokumentation behandelt zentrale Aspekte moderner Festkörper-DFT, darunter Kohn-Sham-Gleichungen, Self-Consistent-Field-Verfahren, Austausch-Korrelationsfunktionale, DFT+U, Hybridfunktionale, Spin-Bahn-Kopplung, Bandstrukturen und elektronische Zustandsdichten. Sie ist besonders nützlich, um den Zusammenhang zwischen theoretischen Methoden und ihrer numerischen Implementierung nachvollziehbar darzustellen.
- Libxc: Library of Exchange-Correlation Functionals for Density-Functional Theory.
- Libxc stellt eine umfangreiche, softwareunabhängige Sammlung von Austausch-, Korrelations- und kinetischen Energiefunktionalen zur Verfügung. Die Ressource ist besonders hilfreich, um Funktionale eindeutig zu identifizieren, ihre Klassifikation zu prüfen und die große Vielfalt moderner DFT-Näherungen zu überblicken. Für Abschnitte über LDA, GGA, Meta-GGA und weitere Funktionalklassen ist sie eine wertvolle Recherchehilfe.
Fachjournale und Verlage
- American Physical Society: Physical Review Journals.
- Die Physical-Review-Zeitschriften enthalten einen wesentlichen Teil der historischen und aktuellen Primärliteratur zu Dichtefunktionaltheorie, elektronischer Struktur, Festkörperphysik und Quantenmaterialien. Für eine wissenschaftliche Abhandlung sollte diese Plattform insbesondere zur Prüfung von Originalarbeiten, bibliografischen Angaben und Folgepublikationen genutzt werden.
- AIP Publishing: The Journal of Chemical Physics.
- The Journal of Chemical Physics gehört zu den zentralen Fachzeitschriften für elektronische Strukturtheorie, Quantenchemie und Entwicklung von Dichtefunktionalmethoden. Die Plattform ist insbesondere für Arbeiten zu Austausch-Korrelationsfunktionalen, Hybridmethoden und methodischen Weiterentwicklungen der Kohn-Sham-DFT relevant.
