Die Dichtefunktionaltheorie, kurz DFT, gehört zu den wichtigsten rechnerischen Methoden der modernen Quantenphysik, Quantenchemie und Materialwissenschaft. Ihr zentrales Ziel besteht darin, die elektronische Struktur von Atomen, Molekülen, Festkörpern, Oberflächen und Nanostrukturen mit quantenmechanischen Methoden zu bestimmen, ohne die vollständige Vielteilchen-Wellenfunktion des Systems explizit berechnen zu müssen.

Das grundlegende Problem ist schnell beschrieben: In einem realen Material wechselwirken viele Elektronen gleichzeitig miteinander und mit den Atomkernen. Die exakte quantenmechanische Beschreibung erfolgt prinzipiell durch die Vielteilchen-Schrödinger-Gleichung. Mit wachsender Zahl der Elektronen steigt deren rechnerische Komplexität jedoch extrem stark an. Bereits für Systeme mit einigen Dutzend Elektronen wird eine direkte Behandlung der vollständigen Vielteilchen-Wellenfunktion praktisch unbeherrschbar.

Die Dichtefunktionaltheorie verfolgt deshalb einen anderen Ansatz. Anstelle einer Wellenfunktion, die von den Koordinaten sämtlicher Elektronen abhängt, verwendet sie die Elektronendichte als zentrale Größe. Diese hängt nur von den drei Raumkoordinaten ab. Dadurch wird die Beschreibung eines quantenmechanischen Vielteilchensystems erheblich kompakter.

Die formale Grundlage liefern die Hohenberg-Kohn-Theoreme. Sie zeigen, dass die Grundzustandseigenschaften eines elektronischen Systems grundsätzlich vollständig durch seine Elektronendichte bestimmt werden können. Die Kohn-Sham-Methode macht diese Aussage rechnerisch nutzbar, indem sie das wechselwirkende Vielteilchensystem auf ein Hilfssystem nicht wechselwirkender Elektronen abbildet.

DFT ist heute ein Standardwerkzeug zur Berechnung von Bindungsenergien, Gleichgewichtsstrukturen, elektronischen Bandstrukturen, magnetischen Eigenschaften, Defektzuständen, Oberflächenreaktionen und zahlreichen weiteren Materialeigenschaften. Für die Quantentechnologie ist sie besonders relevant, weil Eigenschaften von Qubit-Materialien, Defektzentren, Halbleitern, zweidimensionalen Materialien, supraleitenden Systemen und topologischen Materialien häufig zunächst auf DFT-Ebene untersucht werden.

Quantenmechanische Grundlagen

Das elektronische Vielteilchenproblem

Ein Material besteht aus positiv geladenen Atomkernen und negativ geladenen Elektronen. Im nichtrelativistischen Fall kann das gesamte System durch einen Hamiltonoperator beschrieben werden, der die kinetischen Energien sowie sämtliche Coulomb-Wechselwirkungen enthält.

Für die elektronische Strukturrechnung wird meistens zunächst die Born-Oppenheimer-Näherung verwendet. Da Atomkerne wesentlich schwerer als Elektronen sind, können ihre Positionen während der Berechnung des elektronischen Zustands näherungsweise als fest betrachtet werden. Die Bewegung der Kerne wird damit zunächst von der wesentlich schnelleren Elektronenbewegung getrennt.

Der elektronische Hamiltonoperator besitzt in atomaren Einheiten vereinfacht die Form:

\(\hat{H} = -\frac{1}{2}\sum_i \nabla_i^2 - \sum_{i,A}\frac{Z_A}{|\mathbf{r}_i-\mathbf{R}_A|} + \frac{1}{2}\sum_{i\neq j}\frac{1}{|\mathbf{r}_i-\mathbf{r}_j|}\)

Der erste Term beschreibt die kinetische Energie der Elektronen. Der zweite Term enthält die Anziehung zwischen Elektronen und Atomkernen. Der dritte Term beschreibt die Coulomb-Abstoßung zwischen den Elektronen.

Die zugehörige Vielteilchen-Schrödinger-Gleichung lautet:

\(\hat{H}\Psi = E\Psi\)

Die Wellenfunktion hängt von den räumlichen Koordinaten und den Spins aller Elektronen ab. Für ein System mit vielen Elektronen entsteht deshalb ein hochdimensionales Problem. Genau diese Dimensionszunahme ist der wesentliche Grund dafür, dass eine direkte Lösung für realistische Materialien kaum möglich ist.

Die Elektronendichte

Die DFT ersetzt die vollständige Vielteilchen-Wellenfunktion als grundlegende Variable durch die Elektronendichte. Für ein System mit einer Elektronenzahl \(N\) kann die Dichte formal aus der Vielteilchen-Wellenfunktion gewonnen werden:

\(n(\mathbf{r}) = N \sum_{\sigma_1,\ldots,\sigma_N}\int |\Psi(\mathbf{r},\sigma_1,\mathbf{r}_2,\sigma_2,\ldots,\mathbf{r}_N,\sigma_N)|^2 d\mathbf{r}_2 \ldots d\mathbf{r}_N\)

Die Elektronendichte besitzt lediglich drei räumliche Variablen. Unabhängig davon, ob das System zehn, hundert oder tausend Elektronen enthält, bleibt ihre grundlegende Dimension gleich.

Die Gesamtzahl der Elektronen ergibt sich aus der Normierung:

\(\int n(\mathbf{r})\,d\mathbf{r} = N\)

Diese Reduktion der fundamentalen Variable ist der konzeptionelle Kern der Dichtefunktionaltheorie.

Historische Entwicklung

Thomas-Fermi-Modell

Die Idee, elektronische Systeme direkt über ihre Dichte zu beschreiben, entstand bereits vor der modernen DFT. Das Thomas-Fermi-Modell behandelt Elektronen näherungsweise als lokales Fermigas und schreibt die kinetische Energie unmittelbar als Funktional der Elektronendichte.

In vereinfachter Form besitzt die Thomas-Fermi-Kinetik die Struktur:

\(T_{\mathrm{TF}}[n] = C_{\mathrm{TF}}\int n(\mathbf{r})^{5/3}\,d\mathbf{r}\)

Damit wurde erstmals deutlich, dass zumindest bestimmte Eigenschaften elektronischer Systeme direkt aus einer Dichtebeschreibung gewonnen werden können. Das Modell war jedoch zu ungenau, um chemische Bindungen und komplexe elektronische Strukturen zuverlässig zu beschreiben.

Dirac-Austausch

Dirac ergänzte das Dichtemodell um eine lokale Näherung für die Austauschenergie. Die entsprechende Struktur lautet:

\(E_x[n] = -C_x\int n(\mathbf{r})^{4/3}\,d\mathbf{r}\)

Diese Arbeiten waren wichtige Vorläufer der späteren Local Density Approximation. Die moderne DFT erhielt ihre mathematisch belastbare Grundlage jedoch erst durch die Arbeiten von Pierre Hohenberg, Walter Kohn und Lu Jeu Sham.

Die Hohenberg-Kohn-Theoreme

Erstes Hohenberg-Kohn-Theorem

Das erste Hohenberg-Kohn-Theorem besagt, dass die Grundzustands-Elektronendichte eines wechselwirkenden Elektronensystems das äußere Potential bis auf eine additive Konstante eindeutig bestimmt. Damit sind auch der Hamiltonoperator und sämtliche Grundzustandseigenschaften prinzipiell durch die Dichte festgelegt.

Diese Aussage ist fundamental. Sie bedeutet, dass die Elektronendichte nicht nur eine vereinfachte Näherungsgröße ist, sondern grundsätzlich genügend Information enthält, um den Grundzustand eines elektronischen Systems zu bestimmen.

Die Gesamtenergie kann deshalb als Funktional der Dichte geschrieben werden:

\(E[n] = F[n] + \int v_{\mathrm{ext}}(\mathbf{r})n(\mathbf{r})\,d\mathbf{r}\)

Das universelle Funktional enthält die kinetische Energie und die Elektron-Elektron-Wechselwirkung:

\(F[n] = T[n] + V_{ee}[n]\)

Universal bedeutet hier, dass dieses Funktional nicht vom konkreten äußeren Potential abhängt.

Zweites Hohenberg-Kohn-Theorem

Das zweite Hohenberg-Kohn-Theorem formuliert ein Variationsprinzip. Für jede physikalisch zulässige Elektronendichte ist die zugehörige Energie größer oder gleich der exakten Grundzustandsenergie. Die tatsächliche Grundzustandsdichte minimiert das Energiefunktional.

\(E[n] \geq E[n_0] = E_0\)

Damit wird die Lösung des Grundzustandsproblems zu einer Minimierungsaufgabe.

Die Hohenberg-Kohn-Theoreme garantieren jedoch nicht, dass die konkrete mathematische Form des universellen Funktionals bekannt ist. Genau hier liegt die zentrale Schwierigkeit der praktischen DFT. Insbesondere die exakte Behandlung von Austausch und Korrelation bleibt unbekannt.

Die Kohn-Sham-Dichtefunktionaltheorie

Grundidee des Kohn-Sham-Ansatzes

Kohn und Sham entwickelten eine Methode, mit der das Hohenberg-Kohn-Konzept praktisch berechenbar wird. Sie führten ein fiktives System nicht wechselwirkender Elektronen ein, das dieselbe Grundzustandsdichte wie das reale wechselwirkende Elektronensystem besitzt.

Statt eine unbekannte kinetische Energie direkt als Funktional der Dichte approximieren zu müssen, wird ein großer Teil der kinetischen Energie über Einteilchenorbitale berechnet. Diese Kohn-Sham-Orbitale erfüllen effektive Einteilchengleichungen.

Kohn-Sham-Orbitale und Elektronendichte

Die Elektronendichte ergibt sich aus den besetzten Kohn-Sham-Orbitalen:

\(n(\mathbf{r}) = \sum_i f_i |\psi_i(\mathbf{r})|^2\)

Dabei bezeichnet \(f_i\) die Besetzung des jeweiligen Zustands.

Die Kohn-Sham-Gleichungen lauten:

\(\left[-\frac{1}{2}\nabla^2 + v_{\mathrm{eff}}(\mathbf{r})\right]\psi_i(\mathbf{r}) = \varepsilon_i\psi_i(\mathbf{r})\)

Das effektive Potential wird aufgeteilt in:

\(v_{\mathrm{eff}}(\mathbf{r}) = v_{\mathrm{ext}}(\mathbf{r}) + v_H(\mathbf{r}) + v_{xc}(\mathbf{r})\)

Das äußere Potential wird im Wesentlichen durch die Atomkerne erzeugt. Das Hartree-Potential beschreibt die klassische elektrostatische Wechselwirkung mit der gesamten Elektronendichte. Das Austausch-Korrelations-Potential enthält alle verbleibenden quantenmechanischen Vielteilcheneffekte.

Das Kohn-Sham-Energiefunktional

Die Gesamtenergie wird in mehrere Beiträge zerlegt:

\(E[n] = T_s[n] + E_H[n] + E_{xc}[n] + \int v_{\mathrm{ext}}(\mathbf{r})n(\mathbf{r})\,d\mathbf{r}\)

Die nicht wechselwirkende kinetische Energie lautet:

\(T_s = -\frac{1}{2}\sum_i f_i\int \psi_i^*(\mathbf{r})\nabla^2\psi_i(\mathbf{r})\,d\mathbf{r}\)

Die Hartree-Energie ist:

\(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 Term \(E_{xc}[n]\) enthält die Austausch-Korrelations-Energie. Er kompensiert unter anderem die Grenzen der klassischen Hartree-Beschreibung und enthält zudem die Differenz zwischen der tatsächlichen kinetischen Energie des wechselwirkenden Systems und der Kohn-Sham-Kinetik.

Selbstkonsistenz

Das effektive Kohn-Sham-Potential hängt von der Elektronendichte ab. Gleichzeitig wird die Elektronendichte aus den Lösungen der Kohn-Sham-Gleichungen bestimmt. Deshalb müssen die Gleichungen selbstkonsistent gelöst werden.

Zunächst wird eine Anfangsdichte angenommen. Daraus wird das effektive Potential berechnet. Anschließend werden die Kohn-Sham-Gleichungen gelöst und aus den resultierenden Orbitalen eine neue Dichte konstruiert. Dieser Vorgang wird wiederholt, bis sich Eingangs- und Ausgangsdichte innerhalb einer vorgegebenen Toleranz nicht mehr wesentlich unterscheiden.

Formal kann der Zyklus als Abbildung dargestellt werden:

\(n^{(k)}(\mathbf{r}) \rightarrow v_{\mathrm{eff}}^{(k)}(\mathbf{r}) \rightarrow \{\psi_i^{(k)}\} \rightarrow n^{(k+1)}(\mathbf{r})\)

Ist Konvergenz erreicht, spricht man von einer selbstkonsistenten Lösung.

Austausch und Korrelation

Austauschenergie

Elektronen sind Fermionen. Die elektronische Vielteilchen-Wellenfunktion muss deshalb bei der Vertauschung zweier identischer Elektronen antisymmetrisch sein. Diese Eigenschaft führt zur Austauschwechselwirkung.

Austausch ist keine zusätzliche klassische Kraft. Er entsteht unmittelbar aus der fermionischen Symmetrie des quantenmechanischen Zustands und beeinflusst räumliche Verteilungen, Bindungen, Magnetismus und viele weitere Eigenschaften.

Korrelationsenergie

Elektronische Korrelation beschreibt Effekte, die über eine unabhängige Einteilchenbeschreibung hinausgehen. Elektronen bewegen sich nicht unabhängig voneinander. Aufgrund ihrer Coulomb-Abstoßung passen sich ihre Aufenthaltswahrscheinlichkeiten gegenseitig an.

Besonders anspruchsvoll ist die Behandlung stark korrelierter Systeme. Dort können mehrere elektronische Konfigurationen energetisch ähnlich wichtig sein. Beispiele sind bestimmte Übergangsmetalloxide, Mott-Isolatoren sowie Materialien mit lokalisierten d- oder f-Elektronen.

Austausch-Korrelations-Potential

Das Austausch-Korrelations-Potential wird als funktionale Ableitung der Austausch-Korrelations-Energie definiert:

\(v_{xc}(\mathbf{r}) = \frac{\delta E_{xc}[n]}{\delta n(\mathbf{r})}\)

Die Qualität einer DFT-Rechnung hängt daher wesentlich davon ab, wie gut das gewählte Austausch-Korrelations-Funktional die relevante Physik des untersuchten Systems erfasst.

Wichtige Austausch-Korrelations-Funktionale

Local Density Approximation

Die Local Density Approximation, LDA, behandelt die Austausch-Korrelations-Energie an jedem Raumpunkt so, als befände sich dort ein homogenes Elektronengas mit derselben lokalen Dichte.

Die allgemeine Struktur lautet:

\(E_{xc}^{\mathrm{LDA}}[n] = \int n(\mathbf{r})\varepsilon_{xc}(n(\mathbf{r}))\,d\mathbf{r}\)

LDA funktioniert insbesondere dann brauchbar, wenn sich die Elektronendichte räumlich relativ langsam verändert. Bei Festkörpern liefert sie häufig vernünftige strukturelle Eigenschaften. Sie neigt jedoch unter anderem dazu, Bindungen teilweise zu stark zu beschreiben.

Generalized Gradient Approximation

Die Generalized Gradient Approximation berücksichtigt zusätzlich den räumlichen Gradienten der Elektronendichte:

\(E_{xc}^{\mathrm{GGA}} = \int f(n(\mathbf{r}),\nabla n(\mathbf{r}))\,d\mathbf{r}\)

Ein verbreitetes Funktional dieser Klasse ist PBE. GGA-Verfahren stellen in vielen Materialrechnungen einen sinnvollen Kompromiss zwischen Rechenaufwand und Genauigkeit dar.

Ihre Grenzen bleiben jedoch deutlich. Besonders Bandlücken von Halbleitern und Isolatoren werden häufig unterschätzt, und stark lokalisierte elektronische Zustände können unzureichend behandelt werden.

Meta-GGA

Meta-GGA-Funktionale verwenden zusätzliche Informationen wie die kinetische Energiedichte oder höhere Ableitungen der Dichte. Damit erhalten sie mehr Informationen über die lokale elektronische Umgebung.

Ein wichtiges Beispiel ist SCAN. Solche Funktionale können für bestimmte Materialklassen strukturelle und energetische Eigenschaften besser beschreiben als einfache GGA-Verfahren, ohne den vollständigen Rechenaufwand eines Hybridfunktionals zu benötigen.

Hybridfunktionale

Hybridfunktionale kombinieren DFT-Austausch mit einem Anteil exakten Hartree-Fock-Austauschs. Schematisch:

\(E_{xc}^{\mathrm{hybrid}} = aE_x^{\mathrm{HF}} + (1-a)E_x^{\mathrm{DFT}} + E_c^{\mathrm{DFT}}\)

Bekannte Beispiele sind B3LYP, PBE0 und HSE. Hybridfunktionale können insbesondere elektronische Energieniveaus und Bandlücken deutlich verbessern. Dafür steigt der Rechenaufwand erheblich.

Dispersionskorrekturen

Lokale und semilokale Funktionale beschreiben langreichweitige Van-der-Waals-Wechselwirkungen nur unzureichend. Dies ist bei Molekülkristallen, Adsorption, geschichteten Materialien und Van-der-Waals-Heterostrukturen relevant.

Eine häufig verwendete Strategie besteht darin, dem DFT-Energieterm eine Dispersionskorrektur hinzuzufügen. Vereinfacht besitzt ein solcher Beitrag die Form:

\(E_{\mathrm{disp}} = -\sum_{A<B}\frac{C_{6,AB}}{R_{AB}^6}f_{\mathrm{damp}}(R_{AB})\)

Alternativ existieren nichtlokale Funktionale, in denen die Dispersionsphysik unmittelbar in die Dichtebeschreibung integriert wird.

DFT+U

Für stark lokalisierte d- und f-Zustände wird häufig DFT+U eingesetzt. Dabei erhält ein ausgewählter lokalisierter Orbitalraum eine zusätzliche Hubbard-artige Korrektur.

Eine stark vereinfachte Form lautet:

\(E_{\mathrm{DFT}+U} = E_{\mathrm{DFT}} + E_U - E_{\mathrm{dc}}\)

Der Ansatz kann Lokalisierung, magnetische Zustände und Bandlücken verbessern. Die Ergebnisse hängen jedoch teilweise empfindlich vom verwendeten Parameter \(U\) und seiner physikalischen Bestimmung ab.

Numerische Umsetzung

Basisfunktionen

Die Kohn-Sham-Orbitale müssen numerisch dargestellt werden. Dafür stehen unterschiedliche Basissätze zur Verfügung. Bei periodischen Festkörpern werden häufig ebene Wellen verwendet. In molekularen Rechnungen sind Gauß-Basissätze oder lokalisierte Atomorbitale verbreitet.

Eine ebene Welle besitzt beispielsweise die Form:

\(\phi_{\mathbf{G}}(\mathbf{r}) = \frac{1}{\sqrt{\Omega}}e^{i(\mathbf{k}+\mathbf{G})\cdot\mathbf{r}}\)

Die Anzahl der berücksichtigten ebenen Wellen wird durch einen Energie-Cutoff begrenzt:

\(\frac{|\mathbf{k}+\mathbf{G}|^2}{2} \leq E_{\mathrm{cut}}\)

Ein größerer Cutoff verbessert normalerweise die numerische Genauigkeit, erhöht aber gleichzeitig Speicherbedarf und Rechenzeit.

Pseudopotentiale und PAW

Elektronen in inneren Atomschalen sind stark gebunden und beteiligen sich häufig nur indirekt an chemischen Bindungen. Pseudopotentialmethoden ersetzen deshalb den komplizierten Kernbereich durch ein effektives Potential und behandeln hauptsächlich die Valenzelektronen explizit.

Dies reduziert den Rechenaufwand erheblich. Verbreitet sind normerhaltende Pseudopotentiale, ultrasofte Pseudopotentiale und die Projector-Augmented-Wave-Methode.

PAW verbindet viele Vorteile von Pseudopotential- und All-Electron-Verfahren. Im kernnahen Bereich kann die vollständige Wellenfunktion aus einer glatteren Hilfsbeschreibung rekonstruiert werden.

Periodische Randbedingungen

Kristalline Festkörper werden normalerweise über periodische Randbedingungen beschrieben. Anstatt einen makroskopischen Kristall explizit abzubilden, wird eine Einheitszelle periodisch im Raum wiederholt.

Nach dem Bloch-Theorem lassen sich die elektronischen Zustände schreiben als:

\(\psi_{n\mathbf{k}}(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}}u_{n\mathbf{k}}(\mathbf{r})\)

Dabei besitzt \(u_{n\mathbf{k}}(\mathbf{r})\) die Periodizität des Kristallgitters.

Brillouin-Zone und k-Punkte

Elektronische Eigenschaften eines periodischen Festkörpers hängen vom Wellenvektor in der Brillouin-Zone ab. Numerisch wird die Brillouin-Zone durch ein diskretes Netz von k-Punkten abgetastet.

Für viele Kristalle werden Monkhorst-Pack-Gitter verwendet. Die Dichte des k-Punkt-Netzes muss durch Konvergenztests überprüft werden. Metalle benötigen wegen der komplexen Struktur der Zustände nahe der Fermi-Energie häufig dichtere Netze als Isolatoren.

Strukturrelaxation

DFT kann nicht nur elektronische Zustände berechnen, sondern auch die energetisch bevorzugten Atompositionen bestimmen. Die Kraft auf einen Kern ergibt sich grundsätzlich aus der Ableitung der Energie nach seiner Position:

\(\mathbf{F}_A = -\frac{\partial E}{\partial \mathbf{R}_A}\)

Durch iterative Verschiebung der Atomkerne wird eine Struktur gesucht, bei der die Kräfte ausreichend klein sind. Zusätzlich können Gitterparameter und Zellformen optimiert werden.

Eigenschaften aus DFT-Berechnungen

Gesamtenergie und Struktur

Die Gesamtenergie ist eine der wichtigsten DFT-Größen. Unterschiede zwischen Gesamtenergien ermöglichen Aussagen über relative Stabilitäten verschiedener Strukturen, Kristallphasen, Defektkonfigurationen oder Adsorptionsgeometrien.

Aus Energieverläufen können Gleichgewichtsabstände, Bindungsenergien und Reaktionsenergien bestimmt werden. Dabei sind meist Energiedifferenzen zuverlässiger als absolute Gesamtenergien.

Zustandsdichte

Die elektronische Zustandsdichte beschreibt, wie viele elektronische Zustände in einem bestimmten Energiebereich verfügbar sind. Formal kann sie geschrieben werden als:

\(D(E) = \sum_i \delta(E-\varepsilon_i)\)

In Festkörpern wird zusätzlich über die k-Punkte summiert. Die projizierte Zustandsdichte erlaubt es, Beiträge bestimmter Atome oder Orbitale zu analysieren.

Bandstruktur

Die Bandstruktur zeigt die Kohn-Sham-Eigenwerte entlang ausgewählter Richtungen in der Brillouin-Zone:

\(\varepsilon_n = \varepsilon_n(\mathbf{k})\)

Sie liefert Informationen darüber, ob ein Material metallisch, halbleitend oder isolierend ist. Bei Halbleitern lassen sich direkte und indirekte Bandlücken unterscheiden.

Die quantitative Interpretation von Kohn-Sham-Bandlücken verlangt jedoch Vorsicht. Viele Standardfunktionale unterschätzen fundamentale Bandlücken deutlich.

Magnetische Eigenschaften

Spinpolarisierte DFT unterscheidet Elektronendichten verschiedener Spinrichtungen:

\(n(\mathbf{r}) = n_{\uparrow}(\mathbf{r}) + n_{\downarrow}(\mathbf{r})\)

Die Magnetisierungsdichte kann entsprechend definiert werden:

\(m(\mathbf{r}) = n_{\uparrow}(\mathbf{r}) - n_{\downarrow}(\mathbf{r})\)

Damit lassen sich magnetische Momente sowie ferro-, ferri- und antiferromagnetische Konfigurationen vergleichen.

Defekte

Defekte spielen in der Quantentechnologie eine besonders wichtige Rolle. Leerstellen, Fremdatome oder antisite-Defekte können neue elektronische Zustände in einer Bandlücke erzeugen und lokale Spinzentren bilden.

Die Bildungsenergie eines Defekts kann schematisch geschrieben werden als:

\(E_f = E_{\mathrm{def}} - E_{\mathrm{bulk}} - \sum_i n_i\mu_i + q(E_F + E_{\mathrm{VBM}}) + E_{\mathrm{corr}}\)

Bei geladenen Defekten müssen zusätzliche Korrekturen berücksichtigt werden, weil periodische Superzellen künstliche Wechselwirkungen zwischen periodisch wiederholten Defekten erzeugen können.

DFT in der Materialwissenschaft

Materialdesign auf atomarer Ebene

DFT ermöglicht es, Materialeigenschaften ausgehend von der atomaren Zusammensetzung und Struktur vorherzusagen. Dadurch können Kandidaten theoretisch untersucht werden, bevor sie experimentell synthetisiert werden.

Typische Fragestellungen betreffen die Stabilität konkurrierender Kristallstrukturen, elektronische Eigenschaften, magnetische Ordnungen, Defekte, Oberflächen, Adsorptionsprozesse und Phasenübergänge.

Im Hochdurchsatzverfahren werden Tausende oder Millionen struktureller Kandidaten automatisiert berechnet. Solche Datensätze bilden wiederum eine wichtige Grundlage für datengetriebene Materialforschung.

Halbleiter

Bei Halbleitern liefert DFT Informationen über Bandstrukturen, effektive elektronische Zustände, Defekte und Dotierungen. Besonders wichtig ist die Position von Defektzuständen relativ zu Valenz- und Leitungsband.

Für quantitative Bandlücken sind Standard-GGA-Rechnungen jedoch häufig nicht ausreichend. Hybridfunktionale oder Vielteilchenmethoden wie GW können hier notwendig werden.

Oberflächen und Grenzflächen

Viele quantentechnologische Bauelemente werden durch Grenzflächen bestimmt. DFT kann untersuchen, wie sich elektronische Zustände an Materialübergängen verändern, ob Ladungstransfer auftritt und wie stark Atome oder Moleküle an Oberflächen gebunden sind.

Grenzflächen können neue elektronische Zustände erzeugen, Symmetrien brechen oder lokale elektrische Felder verursachen. Diese Effekte können für Qubits und Quantenbauelemente entscheidend sein.

DFT in der Quantentechnologie

Defektzentren als Quantensysteme

Bestimmte Punktdefekte in Festkörpern besitzen lokalisierte elektronische Zustände mit kontrollierbaren Spin-Freiheitsgraden. Sie können als Qubits, Quantensensoren oder Einzelphotonenquellen eingesetzt werden.

DFT wird verwendet, um Defektgeometrien, Ladungszustände, Spinverteilungen, Energieniveaus und lokale Verzerrungen zu untersuchen. Besonders wichtig ist die Frage, ob ein Defekt stabile elektronische Zustände innerhalb der Bandlücke besitzt.

Bekannte Beispiele für die grundlegende Materialklasse sind Farbzentren in Diamant, Siliziumkarbid und weiteren Wide-Bandgap-Materialien. Für eine vollständige Beschreibung optischer Übergänge reicht Grundzustands-DFT allerdings häufig nicht aus.

Materialien für Spin-Qubits

Spin-Qubits reagieren empfindlich auf ihre atomare und elektronische Umgebung. DFT kann Informationen über lokale Magnetfelder, Hyperfeinwechselwirkungen, Spin-Bahn-Kopplung und die elektronische Struktur von Defekten oder Dotieratomen liefern.

Besonders interessant sind Materialien mit geringen Konzentrationen magnetisch aktiver Kerne sowie Systeme, in denen Spin-Zustände räumlich gut lokalisiert werden können.

Supraleitende Quantentechnologie

Supraleitende Qubits werden auf makroskopischer Ebene durch Schaltkreismodelle beschrieben. Auf Materialebene beeinflussen jedoch Grenzflächen, Oxide, Defekte und Verunreinigungen die Kohärenzzeiten.

DFT kann helfen, mögliche mikroskopische Verlustmechanismen zu charakterisieren. Dazu gehören Defektzustände in amorphen oder kristallinen Oxiden, Adsorbate an Oberflächen und elektronische Zustände an Grenzflächen.

Für die eigentliche supraleitende Vielteilchenphysik werden allerdings zusätzliche theoretische Ansätze benötigt.

Zweidimensionale Materialien

Graphen, Übergangsmetalldichalkogenide und andere zweidimensionale Materialien besitzen elektronische Eigenschaften, die sich stark von ihren dreidimensionalen Ausgangsmaterialien unterscheiden können.

DFT ermöglicht die Berechnung ihrer Bandstrukturen, Defektzustände, magnetischen Eigenschaften und Grenzflächen. In Van-der-Waals-Heterostrukturen können unterschiedliche atomar dünne Materialien gezielt kombiniert werden.

Topologische Materialien

Topologische Materialphasen entstehen aus der globalen Struktur elektronischer Zustände. Spin-Bahn-Kopplung spielt dabei häufig eine zentrale Rolle.

DFT liefert die elektronischen Bandstrukturen, aus denen Hinweise auf Bandinversionen und topologisch relevante Zustände gewonnen werden können. Für eine vollständige Klassifikation werden anschließend geeignete topologische Invarianten oder verwandte Größen bestimmt.

Spin, Magnetismus und relativistische Effekte

Spinpolarisierte DFT

Bei magnetischen Systemen muss die Dichtebeschreibung um die Spinstruktur erweitert werden. Die Energie wird dann als Funktional beider Spin-Dichten aufgefasst:

\(E = E[n_{\uparrow},n_{\downarrow}]\)

Durch Vergleich verschiedener magnetischer Anfangskonfigurationen können unterschiedliche magnetische Lösungen identifiziert werden.

Nichtkollinearer Magnetismus

In komplexen magnetischen Systemen müssen Spins nicht parallel oder antiparallel zu einer gemeinsamen Achse ausgerichtet sein. Nichtkollineare DFT beschreibt die lokale Magnetisierung deshalb als Vektorgröße.

Diese Erweiterung ist unter anderem bei Spinspiralen, frustrierten Magneten und Systemen mit starker Spin-Bahn-Kopplung relevant.

Spin-Bahn-Kopplung

Die Spin-Bahn-Kopplung verbindet den Spin eines Elektrons mit seiner orbitalen Bewegung im elektrischen Feld der Atome. Sie wird bei schweren Elementen besonders wichtig.

Auf vereinfachter Ebene kann eine Spin-Bahn-Wechselwirkung geschrieben werden als:

\(\hat{H}_{\mathrm{SOC}} = \lambda\,\mathbf{L}\cdot\mathbf{S}\)

Spin-Bahn-Kopplung beeinflusst Bandaufspaltungen, magnetische Anisotropien und topologische Eigenschaften. Für Spintronik und zahlreiche Quantenmaterialien ist ihre Berücksichtigung daher unverzichtbar.

Grenzen der Dichtefunktionaltheorie

Das unbekannte exakte Funktional

DFT ist formal eine exakte Grundzustandstheorie. Praktische DFT-Rechnungen sind jedoch nur so genau wie die verwendeten Näherungen für Austausch und Korrelation.

Es gibt kein allgemein verfügbares approximatives Funktional, das für sämtliche chemischen, physikalischen und materialwissenschaftlichen Probleme gleichermaßen zuverlässig ist.

Bandlückenproblem

Lokale und semilokale Funktionale unterschätzen die fundamentale Bandlücke vieler Halbleiter und Isolatoren. Ein Grund liegt darin, dass die fundamentale Lücke nicht einfach mit der Differenz zweier gewöhnlicher Kohn-Sham-Eigenwerte identisch ist.

Formal wird die fundamentale Lücke über Ionisationsenergie und Elektronenaffinität definiert:

\(E_g = I - A\)

Sie kann auch über Grundzustandsenergien unterschiedlicher Elektronenzahlen geschrieben werden:

\(E_g = E(N+1) + E(N-1) - 2E(N)\)

Hybridfunktionale und GW-Rechnungen liefern häufig bessere Werte, sind aber rechenintensiver.

Selbstwechselwirkungsfehler

In der Hartree-Energie wechselwirkt die Elektronendichte formal auch mit sich selbst. Im exakten Funktional wird diese Selbstwechselwirkung vollständig kompensiert. Viele Näherungsfunktionale erreichen diese Kompensation jedoch nicht vollständig.

Dadurch können elektronische Zustände künstlich delokalisiert werden. Besonders problematisch ist dies bei lokalisierten Ladungen, Defekten und Übergangsmetallverbindungen.

Stark korrelierte Systeme

Systeme mit starker elektronischer Korrelation können mit konventioneller DFT qualitativ falsch beschrieben werden. Ein bekanntes Beispiel sind Mott-Isolatoren, die aufgrund von Elektron-Elektron-Wechselwirkungen isolierend sind, obwohl eine einfache Einteilchenbeschreibung metallisches Verhalten erwarten lassen kann.

DFT+U, Hybridfunktionale oder Kombinationen von DFT mit Dynamical Mean-Field Theory können solche Systeme teilweise besser behandeln.

Angeregte Zustände

Konventionelle Hohenberg-Kohn-DFT ist eine Grundzustandstheorie. Optische Anregungen, Exzitonen und dynamische elektronische Prozesse können deshalb nicht allgemein direkt aus einer Grundzustandsrechnung abgeleitet werden.

Hierfür werden Erweiterungen wie TDDFT, GW und die Bethe-Salpeter-Gleichung eingesetzt.

Numerische Fehlerquellen

Neben den physikalischen Approximationen existieren rein numerische Fehler. Dazu gehören unzureichende Basisgrößen, zu kleine Energie-Cutoffs, zu grobe k-Punkt-Gitter, zu kleine Superzellen oder ungeeignete Pseudopotentiale.

Eine DFT-Rechnung ist deshalb erst dann belastbar, wenn die relevanten numerischen Parameter systematisch auf Konvergenz geprüft wurden.

Weiterführende Methoden

Zeitabhängige Dichtefunktionaltheorie

Die zeitabhängige Dichtefunktionaltheorie erweitert das Dichtekonzept auf zeitabhängige Systeme. Die Elektronendichte wird damit zu:

\(n = n(\mathbf{r},t)\)

TDDFT wird unter anderem zur Berechnung elektronischer Anregungen, optischer Spektren und zeitabhängiger Reaktionen auf äußere Felder eingesetzt.

GW-Methode

Die GW-Methode gehört zur Vielteilchen-Störungstheorie. Sie beschreibt Quasiteilchenenergien über die elektronische Selbstenergie. Schematisch wird diese als Produkt der Einteilchen-Green-Funktion und der abgeschirmten Coulomb-Wechselwirkung approximiert:

\(\Sigma \approx iGW\)

GW wird häufig auf Basis einer zuvor durchgeführten DFT-Rechnung eingesetzt. Besonders bei Bandlücken und elektronischen Anregungsenergien kann die Methode deutliche Verbesserungen liefern.

Bethe-Salpeter-Gleichung

Die Bethe-Salpeter-Gleichung berücksichtigt die Wechselwirkung angeregter Elektronen mit zurückbleibenden Löchern. Dadurch können gebundene Elektron-Loch-Paare, sogenannte Exzitonen, beschrieben werden.

Für optische Spektren von Halbleitern, Isolatoren und zweidimensionalen Materialien ist dieser Effekt häufig entscheidend.

DFT und DMFT

Die Dynamical Mean-Field Theory behandelt lokale dynamische Korrelationen expliziter als gewöhnliche DFT. In der Kombination DFT+DMFT liefert die DFT die elektronische Grundstruktur, während ein ausgewählter stark korrelierter Orbitalbereich mit einem anspruchsvolleren Vielteilchenmodell behandelt wird.

Diese Methode ist besonders für Übergangsmetallverbindungen und Materialien mit stark korrelierten d- oder f-Elektronen relevant.

DFT und maschinelles Lernen

DFT als Datenquelle

Moderne Materialinformatik verwendet große Mengen berechneter DFT-Daten. Strukturen, Gesamtenergien, Kräfte, Bandlücken und weitere Eigenschaften werden als Trainingsdaten für Machine-Learning-Modelle genutzt.

Der Vorteil liegt in der Skalierbarkeit. Ist ein Modell ausreichend gut trainiert, können bestimmte Eigenschaften wesentlich schneller vorhergesagt werden als mit einer vollständigen DFT-Berechnung.

Machine-Learning-Potentiale

Für atomistische Simulationen werden zunehmend Machine-Learning-Potentiale entwickelt. Sie approximieren die potentielle Energiefläche eines Systems auf Basis quantenmechanischer Referenzrechnungen.

Die Kräfte folgen anschließend aus dem Energiegradienten:

\(\mathbf{F}_i = -\nabla_i E_{\mathrm{ML}}\)

Damit können Molekulardynamik-Simulationen mit wesentlich größeren Systemen und längeren Zeitskalen durchgeführt werden, während eine Genauigkeit nahe der verwendeten Referenzmethode angestrebt wird.

Grenzen datengetriebener Modelle

Machine Learning ersetzt die zugrunde liegende Quantenmechanik nicht automatisch. Ein Modell ist nur innerhalb des chemischen und strukturellen Bereichs zuverlässig, den seine Trainingsdaten ausreichend abdecken.

DFT bleibt deshalb eine zentrale Quelle physikalisch fundierter Referenzdaten.

DFT im Zusammenhang mit Quantencomputern

DFT ist kein Quantenalgorithmus

DFT wird heute überwiegend auf klassischen Hochleistungsrechnern ausgeführt. Der Begriff Quantenmechanik darf daher nicht mit Quantencomputing verwechselt werden. DFT löst ein quantenmechanisches Problem mit klassischen numerischen Verfahren.

DFT als Bestandteil hybrider Workflows

Quantencomputer könnten künftig ausgewählte Teile elektronischer Strukturprobleme übernehmen, die für klassische Methoden besonders schwierig sind. Ein realistischer Ansatz besteht darin, DFT zunächst für die Gesamtstruktur eines Materials einzusetzen und nur einen stark korrelierten Teilbereich mit einem Quantenalgorithmus zu behandeln.

DFT kann beispielsweise Geometrien, Orbitale und effektive elektronische Räume bereitstellen. Daraus lässt sich ein reduzierter Vielteilchen-Hamiltonoperator konstruieren, der anschließend mit einem quantencomputerbasierten Verfahren untersucht wird.

Variational Quantum Eigensolver

Der Variational Quantum Eigensolver verfolgt ein Variationsprinzip. Ein parametrisierter Quantenzustand wird vorbereitet, und die Erwartungsenergie wird minimiert:

\(E(\boldsymbol{\theta}) = \frac{\langle\Psi(\boldsymbol{\theta})|\hat{H}|\Psi(\boldsymbol{\theta})\rangle}{\langle\Psi(\boldsymbol{\theta})|\Psi(\boldsymbol{\theta})\rangle}\)

Die Optimierung der Parameter erfolgt in hybriden Varianten teilweise mit klassischen Algorithmen. Für praktische Materialrechnungen sind jedoch noch erhebliche Fortschritte bei Hardware, Fehlerraten, Schaltkreistiefen und Skalierbarkeit erforderlich.

Langfristige Perspektive

Quantencomputer werden DFT voraussichtlich nicht einfach vollständig ersetzen. Wahrscheinlicher ist eine Arbeitsteilung. DFT bleibt für große Systeme und Materialscreenings effizient, während Quantenalgorithmen gezielt Bereiche behandeln könnten, in denen starke Korrelationen oder komplexe Vielteilcheneffekte dominieren.

Praktischer Ablauf einer DFT-Berechnung

Definition des Systems

Am Anfang steht eine klar formulierte physikalische Fragestellung. Soll eine Kristallstruktur optimiert, eine Bandlücke bestimmt, ein Defekt untersucht oder eine Adsorptionsenergie berechnet werden? Die Fragestellung entscheidet über Modell, Funktional und numerische Genauigkeit.

Auswahl des Funktionals

Für einen ersten Strukturvergleich wird häufig ein GGA-Funktional eingesetzt. Bei Bandlücken können Hybridfunktionale notwendig sein. Bei lokalisierten d- und f-Zuständen kommen beispielsweise DFT+U oder weiterführende Methoden infrage.

Das Funktional sollte daher nicht aus Gewohnheit, sondern anhand des physikalischen Problems gewählt werden.

Konvergenztests

Vor einer umfangreichen Produktionsrechnung müssen die numerischen Parameter geprüft werden. Typische Konvergenztests betreffen Energie-Cutoff, k-Punkt-Dichte, Superzellengröße und elektronische Konvergenzkriterien.

Ist beispielsweise eine Energieänderung kleiner als eine gewünschte Toleranz, kann ein Parameter als ausreichend konvergiert betrachtet werden:

\(|E_{m+1}-E_m| < \epsilon\)

Strukturrelaxation

Die Geometrie wird optimiert, bis Kräfte und gegebenenfalls mechanische Spannungen unter festgelegten Grenzwerten liegen.

Ein mögliches Kraftkriterium lautet:

\(\max_i |\mathbf{F}_i| < F_{\mathrm{tol}}\)

Elektronische Analyse

Auf Basis der relaxierten Struktur können anschließend Zustandsdichte, Bandstruktur, Ladungsdichte, magnetische Momente oder Defektzustände bestimmt werden.

Für reproduzierbare Ergebnisse müssen sämtliche relevanten Rechenparameter dokumentiert werden. Dazu gehören Funktional, Pseudopotentiale, Cutoff, k-Punkt-Netz, Zellgröße, Spinbehandlung und Konvergenzkriterien.

Wichtige DFT-Software

VASP

VASP ist eines der am weitesten verbreiteten Programme für elektronische Strukturrechnungen an Festkörpern, Oberflächen und Materialien. Es verwendet insbesondere ebene Wellen und die PAW-Methode und unterstützt zahlreiche Austausch-Korrelations-Funktionale.

Quantum ESPRESSO

Quantum ESPRESSO ist ein Open-Source-Paket für periodische elektronische Strukturrechnungen. Es bietet Funktionen für Grundzustands-DFT, Phononen, Molekulardynamik und zahlreiche weiterführende Materialberechnungen.

ABINIT

ABINIT ist ebenfalls ein umfangreiches Open-Source-System für Festkörper- und Materialsimulationen. Neben DFT stehen Erweiterungen für Antwortfunktionen, Gitterschwingungen und Vielteilchenmethoden zur Verfügung.

GPAW

GPAW verwendet die PAW-Methode und unterstützt verschiedene numerische Darstellungen, darunter Realraumgitter, ebene Wellen und lokalisierte Orbitale. Die enge Verbindung zum atomistischen Simulationsumfeld ASE erleichtert automatisierte Workflows.

CP2K

CP2K ist besonders für große molekulare und kondensierte Systeme sowie Ab-initio-Molekulardynamik geeignet. Es kombiniert unterschiedliche Basisdarstellungen und ermöglicht effiziente Simulationen komplexer Systeme.

All-Electron-Programme

Programme wie WIEN2k und FHI-aims ermöglichen elektronische Strukturrechnungen ohne klassische Pseudopotentialnäherung. Solche Verfahren sind insbesondere dann wichtig, wenn kernnahe elektronische Eigenschaften oder besonders hohe Genauigkeit benötigt werden.

Materialdatenbanken und Hochdurchsatzrechnungen

Automatisierte DFT

Die klassische DFT-Rechnung untersucht häufig einzelne Materialien detailliert. Moderne Hochdurchsatzverfahren automatisieren dagegen die Berechnung sehr großer Materialmengen.

Ein typischer Workflow umfasst Strukturimport, Relaxation, selbstkonsistente DFT-Rechnung, elektronische Analyse und Speicherung der resultierenden Materialeigenschaften.

Bedeutung großer Datenbestände

Große Materialdatenbanken ermöglichen systematische Vergleiche zwischen Tausenden von Verbindungen. Dadurch können Trends erkannt und Kandidaten für Batterien, Halbleiter, Katalysatoren oder Quantenmaterialien gezielt ausgewählt werden.

Zu den bekannten Ressourcen gehören Materials Project, NOMAD, AFLOW und die Open Quantum Materials Database.

Die Vergleichbarkeit solcher Daten hängt wesentlich davon ab, ob konsistente Funktionale, Pseudopotentiale und numerische Parameter verwendet wurden.

Einordnung und wissenschaftliche Bedeutung

Die große Stärke der Dichtefunktionaltheorie liegt im Verhältnis zwischen physikalischer Aussagekraft und Rechenaufwand. Sie macht quantenmechanische Berechnungen für Systeme möglich, die für explizite Vielteilchenmethoden wesentlich zu groß wären.

DFT ist dabei weder eine einfache empirische Näherung noch eine universell exakte Rechenmethode. Die Grundzustandstheorie besitzt eine exakte formale Grundlage, während praktische Berechnungen von approximativen Austausch-Korrelations-Funktionalen abhängen.

Für viele Struktur- und Energieprobleme liefert DFT hohe Genauigkeit. In anderen Bereichen, insbesondere bei stark korrelierten Elektronensystemen, fundamentalen Bandlücken, angeregten Zuständen und langreichweitigen Wechselwirkungen, müssen geeignete Erweiterungen oder alternative Methoden verwendet werden.

Entscheidend ist deshalb nicht nur die Durchführung einer DFT-Rechnung, sondern die physikalisch begründete Auswahl des Modells. Ein numerisch konvergiertes Ergebnis ist nicht automatisch physikalisch korrekt. Die Grenzen des verwendeten Funktionals müssen Teil der Interpretation sein.

Fazit

Die Dichtefunktionaltheorie hat die elektronische Strukturrechnung grundlegend verändert. Anstelle der vollständigen Vielteilchen-Wellenfunktion stellt sie die wesentlich kompaktere Elektronendichte in den Mittelpunkt. Die Hohenberg-Kohn-Theoreme liefern dafür die formale Grundlage, während die Kohn-Sham-Methode eine praktische numerische Umsetzung ermöglicht.

Der entscheidende Schritt besteht in der Abbildung des wechselwirkenden Elektronensystems auf ein Hilfssystem nicht wechselwirkender Elektronen mit derselben Grundzustandsdichte. Die daraus resultierenden Kohn-Sham-Gleichungen können mit numerischen Methoden für komplexe Materialien gelöst werden.

Die zentrale Näherung steckt im Austausch-Korrelations-Funktional. LDA, GGA, Meta-GGA, Hybridfunktionale, Dispersionskorrekturen und DFT+U stellen unterschiedliche Strategien dar, um relevante quantenmechanische Effekte abzubilden. Keine dieser Methoden ist für sämtliche Systeme optimal.

DFT ermöglicht die Berechnung von Strukturen, Energien, Bandstrukturen, Zustandsdichten, magnetischen Eigenschaften, Defekten, Oberflächen und Grenzflächen. Damit ist sie ein wesentliches Werkzeug der Materialwissenschaft und der modernen Quantentechnologie.

Besonders wichtig ist ihre Rolle bei der Entwicklung neuer Quantenmaterialien. Defektzentren, Spin-Qubits, zweidimensionale Systeme, topologische Materialien und Grenzflächen können auf atomarer Ebene untersucht werden. Gleichzeitig liefert DFT die Ausgangsdaten für weiterführende Verfahren wie GW, Bethe-Salpeter-Rechnungen, DMFT und zunehmend auch Machine-Learning-Modelle.

Auch im Zusammenhang mit Quantencomputern bleibt DFT relevant. Quantenalgorithmen und klassische DFT verfolgen unterschiedliche Ansätze und werden sich langfristig voraussichtlich ergänzen. DFT eignet sich weiterhin für große elektronische Systeme, während Quantencomputer künftig besonders anspruchsvolle korrelierte Teilprobleme übernehmen könnten.

Die Dichtefunktionaltheorie ist damit keine abgeschlossene Einzelmethode, sondern ein zentraler Bestandteil moderner computergestützter Quantenwissenschaft. Ihre Stärke liegt darin, quantenmechanische Materialeigenschaften mit vertretbarem Rechenaufwand zugänglich zu machen und gleichzeitig als Ausgangspunkt für komplexere Vielteilchen- und Quantentechnologieverfahren zu dienen.

Mit freundlichen Grüßen Jörg-Owe Schneppat

Anhang

Wissenschaftliche Zeitschriften und Artikel

Die wissenschaftliche Primärliteratur zur Dichtefunktionaltheorie bildet die Grundlage für ein belastbares Verständnis ihrer mathematischen Fundamente, ihrer praktischen Näherungsverfahren und ihrer Grenzen. Für eine Abhandlung über DFT im Kontext der Quantentechnologie sollten insbesondere die Arbeiten zu den Hohenberg-Kohn-Theoremen, den Kohn-Sham-Gleichungen, Austausch-Korrelations-Funktionalen, Bandlückenproblemen, stark korrelierten Systemen und zeitabhängiger DFT berücksichtigt werden.

Grundlegende Primärliteratur zur Dichtefunktionaltheorie

  • Pierre Hohenberg, Walter Kohn: Inhomogeneous Electron Gas, Physical Review 136, B864–B871, 1964.
    • Diese Arbeit bildet das formale Fundament der modernen Dichtefunktionaltheorie. Die beiden Hohenberg-Kohn-Theoreme zeigen, dass die Grundzustandseigenschaften eines wechselwirkenden Elektronensystems prinzipiell durch dessen Elektronendichte bestimmt werden können. Die Quelle sollte in einer wissenschaftlichen Abhandlung für die theoretische Begründung verwendet werden, warum die Elektronendichte an die Stelle der vollständigen Vielteilchen-Wellenfunktion treten kann.
  • Walter Kohn, Lu Jeu Sham: Self-Consistent Equations Including Exchange and Correlation Effects, Physical Review 140, A1133–A1138, 1965.
    • Die Kohn-Sham-Arbeit macht die abstrakten Hohenberg-Kohn-Theoreme praktisch berechenbar. Sie führt das Hilfssystem nicht wechselwirkender Elektronen und die selbstkonsistenten Kohn-Sham-Gleichungen ein. Diese Primärquelle ist zentral für Abschnitte über Kohn-Sham-Orbitale, effektive Einteilchenpotentiale, Austausch und Korrelation sowie den Self-Consistent-Field-Zyklus moderner DFT-Programme.

Austausch-Korrelations-Funktionale und moderne Näherungsverfahren

  • John P. Perdew, Kieron Burke, Matthias Ernzerhof: Generalized Gradient Approximation Made Simple, Physical Review Letters 77, 3865–3868, 1996.
    • Die Arbeit führt das PBE-Funktional ein, das zu den weltweit meistverwendeten GGA-Funktionalen gehört. Sie ist besonders relevant für die Erklärung, wie über die lokale Elektronendichte hinaus auch deren räumliche Änderung in die Beschreibung von Austausch und Korrelation eingeht. Für materialwissenschaftliche und quantentechnologische DFT-Anwendungen ist PBE ein wichtiger methodischer Referenzpunkt.
  • Jochen Heyd, Gustavo E. Scuseria, Matthias Ernzerhof: Hybrid Functionals Based on a Screened Coulomb Potential, The Journal of Chemical Physics 118, 8207–8215, 2003.
    • Diese Arbeit bildet die Grundlage für die HSE-Familie abgeschirmter Hybridfunktionale. Sie ist für die Abhandlung besonders wichtig, wenn die Grenzen semilokaler Funktionale bei elektronischen Bandstrukturen und Bandlücken diskutiert werden. HSE-Verfahren sind heute in der Festkörper- und Materialphysik eine wichtige Alternative zu einfachen GGA-Funktionalen, wenn eine verbesserte elektronische Struktur benötigt wird.
  • Jianwei Sun, Adrienn Ruzsinszky, John P. Perdew: Strongly Constrained and Appropriately Normed Semilocal Density Functional, Physical Review Letters 115, 036402, 2015.
    • Die Veröffentlichung führt das SCAN-Meta-GGA-Funktional ein. SCAN wurde so konstruiert, dass zahlreiche bekannte exakte Bedingungen für semilokale Funktionale erfüllt werden. Die Arbeit eignet sich als Spezialliteratur für die Entwicklung von LDA über GGA zu Meta-GGA und verdeutlicht, wie zusätzliche Informationen über die elektronische Umgebung genutzt werden können, ohne unmittelbar zu einem Hybridfunktional überzugehen.

Grenzen der DFT und stark korrelierte Elektronensysteme

  • John P. Perdew, Robert G. Parr, Mel Levy, Jose 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 untersucht die DFT für nichtganzzahlige Elektronenzahlen und zeigt die Bedeutung von Ableitungsunstetigkeiten der Energie. Sie ist für eine vertiefte Behandlung des Bandlückenproblems, der fundamentalen Energielücke sowie der Grenzen einfacher Kohn-Sham-Eigenwertinterpretationen von besonderer Bedeutung.
  • 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 häufig verwendete vereinfachte rotationsinvariante DFT+U-Behandlung lokalisierter Elektronenzustände. Sie eignet sich zur Diskussion von Übergangsmetalloxiden, lokalisierten d- und f-Elektronen sowie der Frage, warum konventionelle LDA- oder GGA-Verfahren bei stark korrelierten Materialien scheitern können.

Zeitabhängige Dichtefunktionaltheorie

  • Erich Runge, E. K. U. Gross: Density-Functional Theory for Time-Dependent Systems, Physical Review Letters 52, 997–1000, 1984.
    • Die Runge-Gross-Arbeit bildet das theoretische Fundament der zeitabhängigen Dichtefunktionaltheorie. Sie erweitert die Dichtebeschreibung auf zeitabhängige äußere Potentiale und ist deshalb die zentrale Primärquelle für Abschnitte über TDDFT, elektronische Anregungen, optische Eigenschaften und dynamische Reaktionen quantenmechanischer Systeme.

Bücher und Monographien

Monographien sind für eine DFT-Abhandlung besonders wertvoll, weil sie die Verbindung zwischen der mathematischen Grundlage, den verschiedenen Funktionalklassen und der praktischen numerischen Umsetzung systematisch herstellen. Die folgenden Werke decken unterschiedliche Schwerpunkte ab: von der Vielteilchentheorie und den formalen Grundlagen bis zur praktischen Durchführung von Materialrechnungen.

Standardwerke zu den theoretischen Grundlagen

  • Reiner M. Dreizler, Eberhard K. U. Gross: Density Functional Theory: An Approach to the Quantum Many-Body Problem, Springer, 1990.
    • Dieses Standardwerk behandelt DFT aus der Perspektive des quantenmechanischen Vielteilchenproblems. Es ist besonders geeignet für eine theoretisch anspruchsvolle Darstellung der Hohenberg-Kohn-Theoreme, des Kohn-Sham-Schemas, expliziter Funktionale und relativistischer Erweiterungen. Für die wissenschaftliche Abhandlung dient es vor allem als vertiefende Grundlagenliteratur.
  • Eberhard Engel, Reiner M. Dreizler: Density Functional Theory: An Advanced Course, Springer, 2011.
    • Diese Monographie bietet eine fortgeschrittene und mathematisch fundierte Behandlung moderner DFT. Neben Existenztheoremen und Kohn-Sham-Gleichungen werden Austausch-Korrelations-Funktionale, orbitalabhängige Funktionale, zeitabhängige DFT und relativistische DFT behandelt. Das Werk eignet sich insbesondere zur Absicherung theoretisch anspruchsvoller Abschnitte.
  • Richard M. Martin: Electronic Structure: Basic Theory and Practical Methods, Cambridge University Press, 2004.
    • Martins Werk verbindet elektronische Strukturtheorie mit der praktischen Berechnung von Festkörpern. Es behandelt DFT im größeren Zusammenhang von Bandtheorie, Pseudopotentialen, ebenen Wellen, Brillouin-Zonen und elektronischen Strukturverfahren. Für eine Abhandlung mit Schwerpunkt Quantentechnologie und Materialphysik ist dieses Buch besonders wertvoll.

Praxisorientierte Literatur für DFT-Berechnungen

  • David S. Sholl, Janice A. Steckel: Density Functional Theory: A Practical Introduction, 2. Auflage, Wiley, 2022.
    • Dieses Lehrbuch konzentriert sich auf die praktische Anwendung von DFT, insbesondere auf ebenenwellenbasierte Berechnungen von Festkörpern, Oberflächen und Grenzflächen. Es ist eine geeignete Quelle für Abschnitte zu Konvergenztests, Kristallstrukturen, Oberflächen, Magnetismus und Ab-initio-Simulationen und schlägt eine direkte Brücke zwischen Theorie und tatsächlicher Berechnungspraxis.
  • Wolfram Koch, Max C. Holthausen: A Chemist's Guide to Density Functional Theory, 2. Auflage, Wiley-VCH, 2001.
    • Das Buch erklärt die theoretischen Grundlagen der DFT mit starkem Bezug zu ihrer praktischen Anwendung. Es behandelt Elektronendichte, Hohenberg-Kohn-Theoreme, Kohn-Sham-Theorie, Austausch-Korrelations-Funktionale und die Leistungsfähigkeit verschiedener DFT-Ansätze. Es eignet sich besonders zur verständlichen Einordnung der unterschiedlichen Funktionalklassen und ihrer praktischen Stärken und Schwächen.

Online-Ressourcen und Datenbanken

Online-Ressourcen ergänzen die klassische Fachliteratur durch aktuelle Dokumentationen, reproduzierbare Rechenmethoden und umfangreiche Materialdaten. Für wissenschaftliche Arbeiten sind insbesondere die offiziellen Dokumentationen etablierter DFT-Programme sowie kuratierte Materialdatenbanken nützlich. Sie sollten Primärliteratur nicht ersetzen, können aber zur Überprüfung numerischer Methoden, Rechenparameter und bereits verfügbarer Materialdaten eingesetzt werden.

DFT-Software und technische Dokumentationen

  • Quantum ESPRESSO: Official Documentation, Quantum ESPRESSO Foundation.
    • Quantum ESPRESSO ist ein weit verbreitetes Open-Source-Paket für elektronische Strukturrechnungen mit ebenen Wellen und Pseudopotentialen. Die offizielle Dokumentation ist besonders nützlich für selbstkonsistente Kohn-Sham-Rechnungen, Bandstrukturen, Phononen, Dichtefunktional-Störungstheorie und die praktische Umsetzung periodischer DFT-Berechnungen.
  • VASP: The VASP Manual, Vienna Ab initio Simulation Package.
    • Das offizielle VASP-Wiki bietet detaillierte Informationen zu Kohn-Sham-DFT, PAW, LDA, GGA, Meta-GGA, Hybridfunktionalen, DFT+U, Spin-Bahn-Kopplung, Bandstrukturen und zahlreichen weiterführenden Verfahren. Die Ressource ist besonders hilfreich, um die praktische Umsetzung moderner DFT-Verfahren in der Festkörper- und Quantenmaterialforschung nachzuvollziehen.
  • ABINIT: Official Documentation and Tutorials, ABINIT Project.
    • ABINIT ist ein Open-Source-Paket für elektronische Struktur, DFT, DFT+U, Phononen, GW, Bethe-Salpeter-Verfahren und weitere Methoden. Die umfangreichen Tutorials sind für eine wissenschaftliche Abhandlung besonders nützlich, wenn Zusammenhänge zwischen Grundzustands-DFT und weiterführenden Vielteilchenverfahren erklärt werden sollen.
  • GPAW: Official GPAW Documentation.
    • GPAW ist eine PAW-basierte elektronische Strukturumgebung mit Unterstützung für Realraumgitter, ebene Wellen und lokalisierte Atomorbitale. Die Dokumentation eignet sich zur Vertiefung praktischer DFT-Methoden sowie für weiterführende Themen wie lineare Antwort, elektronische Anregungen und TDDFT.

Materialdatenbanken und Hochdurchsatz-DFT

  • Materials Project: Materials Data and Documentation.
    • Das Materials Project stellt umfangreiche berechnete Materialeigenschaften zur Verfügung und gehört zu den wichtigsten Ressourcen für datengetriebene Materialforschung. Die Datenbank eignet sich zur Recherche von Kristallstrukturen, thermodynamischen Eigenschaften, elektronischen Eigenschaften und bereits berechneten Materialkandidaten. Für eine DFT-Abhandlung ist sie außerdem ein anschauliches Beispiel für die Verbindung von Hochdurchsatz-DFT und moderner Materialentdeckung.
  • NOMAD: Materials Science Data Platform, FAIRmat.
    • NOMAD ist eine offene Forschungsdatenplattform für Materialwissenschaften. Sie verarbeitet Daten unterschiedlicher Simulationsprogramme und ermöglicht deren strukturierte Suche, Analyse und Wiederverwendung. Die Plattform ist besonders relevant für Abschnitte über Reproduzierbarkeit, FAIR-Forschungsdaten, automatisierte Materialsimulation und Machine Learning auf Basis elektronischer Strukturdaten.
  • AFLOW: Automatic FLOW for Materials Discovery.
    • AFLOW verbindet automatisierte elektronische Strukturrechnungen mit einer umfangreichen Materialdatenbank. Die Ressource ist besonders interessant für Hochdurchsatz-DFT, Kristallstrukturen, Bildungsenthalpien, Bandstrukturen, elastische Eigenschaften und automatisiertes Materialscreening. Sie kann in der Abhandlung als Beispiel dafür dienen, wie DFT in großskalige Materialentdeckungsprozesse integriert wird.
  • OQMD: The Open Quantum Materials Database.
    • Die Open Quantum Materials Database stellt DFT-berechnete thermodynamische und strukturelle Eigenschaften einer großen Zahl an Materialien bereit. Sie ist besonders für Vergleiche von Phasenstabilitäten, Bildungsenergien und Materialzusammensetzungen geeignet und demonstriert den praktischen Nutzen standardisierter Hochdurchsatzrechnungen.

Fachjournale und Rechercheplattformen

  • American Physical Society: Physical Review Journals.
    • Die Physical-Review-Zeitschriften enthalten zahlreiche grundlegende und aktuelle Arbeiten zur Dichtefunktionaltheorie, Festkörperphysik, Vielteilchentheorie, Quantenmaterialien und elektronischer Struktur. Für DFT sind insbesondere Physical Review, Physical Review Letters und Physical Review B von hoher Bedeutung. Die Plattform eignet sich vor allem zur gezielten Suche nach Primärliteratur.
  • AIP Publishing: The Journal of Chemical Physics and verwandte Fachpublikationen.
    • AIP Publishing stellt wichtige Fachliteratur zur elektronischen Strukturtheorie, Quantenchemie und Entwicklung von Austausch-Korrelations-Funktionalen bereit. Besonders The Journal of Chemical Physics ist für methodische Arbeiten zur DFT und zu Hybridfunktionalen eine bedeutende Recherchequelle.
  • arXiv: Open-Access Preprint Repository.
    • arXiv ermöglicht den Zugriff auf Preprints aus Physik, Materialwissenschaft, Quantenphysik und angrenzenden Gebieten. Die Plattform ist besonders geeignet, um aktuelle Entwicklungen zu DFT, Quantenmaterialien, elektronischer Struktur, Vielteilchenmethoden und hybriden Quanten-Klassik-Ansätzen zu verfolgen. Bei wissenschaftlichen Zitaten sollte nach Möglichkeit zusätzlich geprüft werden, ob eine begutachtete Journalversion verfügbar ist.
Defekt-Qubits

Defekt-Qubits

...

EU Quantum Flagship

EU Quantum Flagship

...

Anti-Up-Quarks

Anti-Up-Quarks

...

Rotationsgatter

Rotationsgatter

...

Photonische Qubits

Photonische Qubits

...

Repetition Codes

Repetition Codes

...