Die Local Density Approximation, kurz LDA, ist eine der grundlegenden Näherungen innerhalb der Dichtefunktionaltheorie. Sie dient dazu, die Austausch- und Korrelationsenergie eines quantenmechanischen Vielteilchensystems aus der lokalen Elektronendichte zu bestimmen. Anstatt die vollständige Vielteilchen-Wellenfunktion aller Elektronen direkt zu berechnen, wird die räumliche Verteilung der Elektronendichte verwendet. Dadurch lässt sich ein ansonsten extrem aufwendiges quantenmechanisches Problem auf eine rechnerisch wesentlich handhabbarere Form reduzieren.
Die zentrale Annahme der LDA ist einfach: An jedem Punkt im Raum wird die reale, möglicherweise stark inhomogene Elektronendichte so behandelt, als gehöre sie lokal zu einem homogenen Elektronengas mit derselben Dichte. Für jeden Punkt kann dadurch eine bekannte Austausch-Korrelationsenergie angesetzt werden. Anschließend werden die lokalen Beiträge über den gesamten Raum integriert.
Diese Annahme ist insbesondere dann sinnvoll, wenn sich die Elektronendichte räumlich nur langsam verändert. Genau deshalb funktioniert die LDA bei vielen metallischen Festkörpern und kompakten Materialien überraschend gut. Obwohl sie konzeptionell einfach ist, bildete sie über Jahrzehnte eine zentrale Grundlage der rechnergestützten Festkörperphysik und Materialwissenschaft.
Bedeutung für die Quantenphysik und Quantentechnologie
Moderne Quantentechnologien beruhen auf Materialien, deren elektronische Eigenschaften auf quantenmechanischen Effekten basieren. Dazu gehören Halbleiter, supraleitende Materialien, Defektzentren in Kristallen, magnetische Systeme und nanostrukturierte Festkörper. Um solche Systeme gezielt entwickeln zu können, muss ihre elektronische Struktur verstanden werden.
Die LDA ermöglicht unter anderem die Berechnung von Elektronendichten, Gesamtenergien, Gleichgewichtsgeometrien, Kristallstrukturen und elektronischen Zuständen. Sie ist damit kein Verfahren, das ausschließlich aus historischem Interesse betrachtet wird, sondern weiterhin Teil vieler elektronischer Strukturberechnungen und Vergleichsstudien.
Ihre Stärke liegt nicht darin, jedes Materialproblem exakt zu lösen. Entscheidend ist vielmehr, dass sie ein klar definiertes, effizientes und gut verstandenes Näherungsverfahren bereitstellt. An ihr lässt sich besonders deutlich erkennen, wie die Dichtefunktionaltheorie aus einem formal exakten theoretischen Rahmen zu einem praktisch einsetzbaren Berechnungsverfahren wird.
Ziel der Abhandlung
Diese Abhandlung erläutert die physikalischen Grundlagen der LDA, ihre mathematische Formulierung, ihre Einbindung in die Kohn-Sham-Dichtefunktionaltheorie sowie ihre praktischen Anwendungen. Ebenso werden typische Fehlerquellen, Grenzen und Weiterentwicklungen behandelt. Ziel ist es, die LDA nicht isoliert als Formel zu betrachten, sondern als Bestandteil einer umfassenden Methode zur Beschreibung elektronischer Vielteilchensysteme.
Quantenmechanische Grundlagen
Das quantenmechanische Vielteilchenproblem
Ein Material besteht aus Atomkernen und Elektronen, die über elektromagnetische Kräfte miteinander wechselwirken. Bereits bei wenigen Elektronen wird eine direkte quantenmechanische Beschreibung kompliziert. Für ein System mit vielen Elektronen müsste grundsätzlich eine Vielteilchen-Wellenfunktion bestimmt werden, die von den Koordinaten aller beteiligten Elektronen abhängt.
Für ein System aus \(N\) Elektronen besitzt die Wellenfunktion eine Struktur der Form
\(\Psi(\mathbf r_1,\mathbf r_2,\ldots,\mathbf r_N)\)
wobei die Spins der Elektronen zusätzlich berücksichtigt werden müssen. Die Anzahl der Freiheitsgrade nimmt mit wachsender Elektronenzahl stark zu. Dadurch steigt der rechnerische Aufwand so schnell an, dass eine direkte Lösung der Vielteilchen-Schrödinger-Gleichung für realistische Festkörper meist nicht praktikabel ist.
Die zeitunabhängige Schrödinger-Gleichung lautet allgemein
\(\hat H \Psi = E \Psi\)
Der Hamilton-Operator enthält dabei die kinetische Energie der Elektronen, das Potential der Atomkerne und insbesondere die Coulomb-Wechselwirkung zwischen den Elektronen. Gerade diese Elektron-Elektron-Wechselwirkung verhindert in der Regel eine Zerlegung des Problems in völlig unabhängige Einteilchenprobleme.
Elektronendichte statt Vielteilchen-Wellenfunktion
Die Dichtefunktionaltheorie verfolgt einen grundlegend anderen Ansatz. Anstelle der vollständigen Wellenfunktion wird die Elektronendichte
\(n(\mathbf r)\)
als zentrale Größe verwendet. Sie gibt an, wie viele Elektronen sich im Mittel in einem kleinen Raumbereich um den Punkt \(\mathbf r\) befinden.
Während die Vielteilchen-Wellenfunktion eines Systems mit vielen Elektronen von sehr vielen Variablen abhängt, ist die Elektronendichte lediglich eine Funktion der drei Raumkoordinaten. Diese Reduktion ist der entscheidende praktische Vorteil der Dichtefunktionaltheorie.
Die Gesamtzahl der Elektronen ergibt sich aus
\(N = \int n(\mathbf r)\,d^3r\)
Die Elektronendichte enthält dennoch genügend Information, um im Grundzustand prinzipiell die relevanten Eigenschaften des Systems festzulegen. Dieser Zusammenhang bildet die theoretische Grundlage für die gesamte Dichtefunktionaltheorie.
Energie eines elektronischen Systems
Die Gesamtenergie eines elektronischen Vielteilchensystems setzt sich aus mehreren Beiträgen zusammen. Dazu gehören die kinetische Energie der Elektronen, die Wechselwirkung mit einem externen Potential und die Elektron-Elektron-Wechselwirkung.
Das externe Potential stammt bei Atomen, Molekülen und Festkörpern hauptsächlich von den Atomkernen. Zusätzlich kann ein System äußeren elektrischen oder magnetischen Feldern ausgesetzt sein.
Die klassische Coulomb-Wechselwirkung beschreibt die elektrostatische Abstoßung zwischen den Elektronen. Darüber hinaus treten rein quantenmechanische Effekte auf. Dazu gehören der Austausch, der aus der Fermionennatur der Elektronen folgt, sowie Korrelationen zwischen den Bewegungen verschiedener Elektronen.
Gerade Austausch und Korrelation machen eine genaue Behandlung des Vielteilchenproblems schwierig. In der Dichtefunktionaltheorie werden diese komplexen Effekte in einem Austausch-Korrelationsfunktional zusammengefasst.
Dichtefunktionaltheorie als Grundlage der LDA
Grundidee der Dichtefunktionaltheorie
Die Dichtefunktionaltheorie beschreibt die Energie eines elektronischen Systems als Funktional der Elektronendichte. Ein Funktional ordnet einer Funktion einen Zahlenwert zu. Im vorliegenden Fall wird einer bestimmten Elektronendichte eine Gesamtenergie zugeordnet.
Formal lässt sich dies schreiben als
\(E = E[n]\)
Der entscheidende Vorteil besteht darin, dass nicht mehr die vollständige Vielteilchen-Wellenfunktion bestimmt werden muss. Stattdessen wird nach der Elektronendichte gesucht, für die die Energie des Systems minimal ist.
Erstes Hohenberg-Kohn-Theorem
Das erste Hohenberg-Kohn-Theorem besagt, dass die Grundzustands-Elektronendichte das externe Potential eines elektronischen Systems bis auf eine additive Konstante eindeutig bestimmt. Damit legt die Grundzustandsdichte auch den Hamilton-Operator und folglich sämtliche Grundzustandseigenschaften fest.
Diese Aussage ist von grundlegender Bedeutung. Sie zeigt, dass die Elektronendichte nicht lediglich eine vereinfachte Hilfsgröße ist. Im Grundzustand enthält sie prinzipiell die notwendige Information zur vollständigen Charakterisierung des Systems.
Zweites Hohenberg-Kohn-Theorem
Das zweite Hohenberg-Kohn-Theorem liefert ein Variationsprinzip. Für jede zulässige Elektronendichte ist die daraus berechnete Energie größer oder gleich der exakten Grundzustandsenergie.
Formal lässt sich dies als
\(E[n] \geq E_0\)
schreiben, wobei \(E_0\) die exakte Grundzustandsenergie bezeichnet. Das Energieminimum wird bei der exakten Grundzustandsdichte erreicht.
Damit wird das quantenmechanische Grundzustandsproblem zu einem Optimierungsproblem: Gesucht ist die Elektronendichte, welche das Energiefunktional minimiert.
Das unbekannte Funktional
Die Hohenberg-Kohn-Theoreme beweisen die Existenz eines universellen Dichtefunktionals, geben aber nicht dessen vollständige mathematische Form an. Genau hier liegt das zentrale praktische Problem der Dichtefunktionaltheorie.
Insbesondere die exakte Form der Austausch-Korrelationsenergie ist im Allgemeinen unbekannt. Deshalb werden Näherungen benötigt. Die Local Density Approximation stellt eine der einfachsten und wichtigsten dieser Näherungen dar.
Das Kohn-Sham-Verfahren
Konzept der Kohn-Sham-Dichtefunktionaltheorie
Walter Kohn und Lu Jeu Sham entwickelten einen Ansatz, der die Dichtefunktionaltheorie praktisch berechenbar macht. Sie führten ein fiktives System nicht wechselwirkender Elektronen ein, dessen Grundzustandsdichte mit der Dichte des realen wechselwirkenden Systems übereinstimmt.
Die Elektronen dieses Hilfssystems bewegen sich in einem effektiven Potential. Dadurch kann das komplizierte Vielteilchenproblem auf eine Gruppe von Einteilchengleichungen zurückgeführt werden.
Aufteilung der Gesamtenergie
Die Kohn-Sham-Gesamtenergie wird üblicherweise in mehrere Beiträge zerlegt:
\(E[n] = T_s[n] + E_{\mathrm{ext}}[n] + E_H[n] + E_{xc}[n]\)
Dabei bezeichnet \(T_s[n]\) die kinetische Energie des nicht wechselwirkenden Kohn-Sham-Systems. Der Term \(E_{\mathrm{ext}}[n]\) beschreibt die Wechselwirkung mit dem externen Potential.
Die klassische elektrostatische Elektron-Elektron-Wechselwirkung wird durch die Hartree-Energie beschrieben:
\(E_H[n] = \frac{1}{2}\int\int \frac{n(\mathbf r)n(\mathbf r')}{|\mathbf r-\mathbf r'|}\,d^3r\,d^3r'\)
Alle Effekte, die durch diese Terme nicht vollständig erfasst werden, werden im Austausch-Korrelationsfunktional \(E_{xc}[n]\) zusammengefasst.
Austausch-Korrelationsfunktional
Das Austausch-Korrelationsfunktional enthält zwei physikalisch unterschiedliche Beiträge. Der Austausch entsteht durch die antisymmetrische Natur der Wellenfunktion identischer Fermionen. Die Korrelation beschreibt darüber hinausgehende Wechselwirkungen zwischen den Elektronen.
Formal kann geschrieben werden:
\(E_{xc}[n] = E_x[n] + E_c[n]\)
Die Schwierigkeit besteht darin, dass die exakte Form dieser beiden Funktionale für allgemeine inhomogene Systeme nicht bekannt ist. Die LDA setzt genau an dieser Stelle an.
Kohn-Sham-Gleichungen
Durch Variation der Gesamtenergie erhält man die Kohn-Sham-Gleichungen:
\(\left[-\frac{1}{2}\nabla^2 + v_{\mathrm{eff}}(\mathbf r)\right]\phi_i(\mathbf r) = \epsilon_i\phi_i(\mathbf r)\)
Das effektive Potential lautet
\(v_{\mathrm{eff}}(\mathbf r) = v_{\mathrm{ext}}(\mathbf r) + v_H(\mathbf r) + v_{xc}(\mathbf r)\)
Die Elektronendichte wird anschließend aus den besetzten Kohn-Sham-Orbitalen bestimmt:
\(n(\mathbf r) = \sum_i^{\mathrm{besetzt}} |\phi_i(\mathbf r)|^2\)
Das Austausch-Korrelationspotential ergibt sich aus der funktionalen Ableitung
\(v_{xc}(\mathbf r) = \frac{\delta E_{xc}[n]}{\delta n(\mathbf r)}\)
Für praktische Berechnungen muss deshalb eine konkrete Näherung für \(E_{xc}[n]\) gewählt werden.
Grundprinzip der Local Density Approximation
Die zentrale Annahme
Die LDA betrachtet jeden Punkt eines inhomogenen elektronischen Systems so, als befände er sich in einem unendlich ausgedehnten homogenen Elektronengas mit derselben lokalen Elektronendichte.
Für einen Punkt \(\mathbf r\) wird somit angenommen, dass die Austausch-Korrelationsenergie pro Elektron näherungsweise jener eines homogenen Elektronengases mit der Dichte \(n(\mathbf r)\) entspricht.
Diese Idee ist lokal. Informationen darüber, wie stark sich die Dichte in unmittelbarer Umgebung verändert, gehen in der einfachsten Form der LDA nicht ausdrücklich ein.
Homogenes Elektronengas
Das homogene Elektronengas ist ein theoretisches Modellsystem, in dem sich Elektronen in einem konstanten positiven Hintergrund bewegen. Dieser Hintergrund sorgt für elektrische Neutralität, ohne die Elektronen durch diskrete Atomkerne lokal zu binden.
Die Elektronendichte ist räumlich konstant:
\(n(\mathbf r) = n_0\)
Das Modell wird häufig auch als Jellium-Modell bezeichnet. Obwohl reale Materialien nicht vollständig homogen sind, eignet sich dieses System hervorragend als Referenz für die Beschreibung lokaler Austausch- und Korrelationseffekte.
Mathematische Form der LDA
Das Austausch-Korrelationsfunktional der LDA lautet
\(E_{xc}^{\mathrm{LDA}}[n] = \int n(\mathbf r)\epsilon_{xc}(n(\mathbf r))\,d^3r\)
Die Größe \(\epsilon_{xc}(n)\) bezeichnet die Austausch-Korrelationsenergie pro Elektron eines homogenen Elektronengases mit der Dichte \(n\).
Das gesamte reale System wird somit gedanklich in unendlich kleine Raumbereiche zerlegt. Jedem Raumbereich wird entsprechend seiner lokalen Elektronendichte eine Austausch-Korrelationsenergie zugeordnet. Die Summe aller lokalen Beiträge ergibt die gesamte angenäherte Austausch-Korrelationsenergie.
Physikalische Interpretation
Die LDA geht davon aus, dass Austausch- und Korrelationseffekte hauptsächlich durch die unmittelbare lokale Umgebung eines Elektrons bestimmt werden. Räumliche Änderungen der Elektronendichte werden nur indirekt berücksichtigt, nämlich dadurch, dass an unterschiedlichen Punkten unterschiedliche lokale Dichten eingesetzt werden.
Diese Annahme funktioniert besonders gut, wenn die charakteristische Länge, über die sich die Elektronendichte verändert, deutlich größer ist als typische elektronische Abstände. In stark inhomogenen Regionen kann diese Voraussetzung dagegen verletzt sein.
Austauschenergie innerhalb der LDA
Ursprung der Austauschwechselwirkung
Elektronen sind Fermionen und unterliegen dem Pauli-Prinzip. Die vollständige elektronische Wellenfunktion muss bei Vertauschung zweier identischer Elektronen antisymmetrisch sein. Dieser mathematische Zwang führt dazu, dass Elektronen mit gleichem Spin räumlich eine verringerte Wahrscheinlichkeit besitzen, sich am selben Ort aufzuhalten.
Dieser Effekt wird als Austausch bezeichnet. Es handelt sich nicht um eine zusätzliche klassische Kraft, sondern um eine direkte Konsequenz der Quantenstatistik.
Dirac-Austauschfunktional
Für ein homogenes Elektronengas kann die Austauschenergie analytisch bestimmt werden. Daraus ergibt sich das LDA-Austauschfunktional:
\(E_x^{\mathrm{LDA}}[n] = -\frac{3}{4}\left(\frac{3}{\pi}\right)^{1/3}\int n(\mathbf r)^{4/3}\,d^3r\)
Dieser Ausdruck zeigt unmittelbar, dass die Austauschenergie nur von der lokalen Elektronendichte abhängt. Dichtegradienten oder weiter entfernte Bereiche des Systems treten nicht explizit auf.
Bedeutung für die Gesamtenergie
Der Austauschbeitrag ist negativ und senkt die Gesamtenergie. Er korrigiert einen Teil der klassischen Hartree-Beschreibung, in der Elektronen lediglich als kontinuierliche Ladungsverteilungen behandelt werden.
Ohne den Austausch wäre eine physikalisch realistische Beschreibung elektronischer Systeme nicht möglich. Die LDA liefert für diesen Beitrag eine mathematisch besonders kompakte Näherung.
Korrelationsenergie innerhalb der LDA
Elektronenkorrelation
Auch Elektronen mit unterschiedlichem Spin beeinflussen ihre Bewegungen gegenseitig. Aufgrund der Coulomb-Abstoßung versuchen Elektronen generell, ungünstige räumliche Konfigurationen zu vermeiden. Diese dynamische Anpassung der Elektronenbewegung wird als Korrelation bezeichnet.
Die Korrelationsenergie ist im Vergleich zum Austausch häufig kleiner, kann jedoch für quantitative Vorhersagen entscheidend sein.
Numerische Bestimmung
Anders als der Austausch besitzt die Korrelationsenergie des homogenen Elektronengases keine einfache geschlossene exakte Formel für alle Dichten. Deshalb wurden hochpräzise numerische Vielteilchenmethoden eingesetzt.
Besonders wichtig wurden Quantum-Monte-Carlo-Berechnungen des homogenen Elektronengases. Sie lieferten sehr genaue Referenzwerte der Korrelationsenergie für unterschiedliche Elektronendichten.
Parametrisierungen
Damit diese numerischen Daten effizient in DFT-Programmen genutzt werden können, werden sie durch analytische Funktionen parametrisiert. Bekannte Beispiele sind die Perdew-Zunger-Parametrisierung und die Vosko-Wilk-Nusair-Parametrisierung.
Diese Ansätze beschreiben die Korrelationsenergie pro Elektron als Funktion geeigneter Dichteparameter. Häufig wird dafür der Wigner-Seitz-Radius verwendet:
\(\frac{4}{3}\pi r_s^3 n = 1\)
Daraus folgt
\(r_s = \left(\frac{3}{4\pi n}\right)^{1/3}\)
Der Parameter \(r_s\) beschreibt im Wesentlichen den mittleren Abstand zwischen Elektronen in Einheiten einer geeigneten Längenskala.
Gesamtes Austausch-Korrelationsfunktional
Die Austausch- und Korrelationsanteile werden kombiniert:
\(\epsilon_{xc}(n) = \epsilon_x(n) + \epsilon_c(n)\)
Damit erhält man ein vollständiges lokales Austausch-Korrelationsfunktional, das direkt in den Kohn-Sham-Gleichungen eingesetzt werden kann.
Selbstkonsistente Berechnung mit der LDA
Ausgangsdichte
Eine praktische DFT-Berechnung beginnt mit einer geschätzten Elektronendichte. Diese kann beispielsweise aus überlagerten atomaren Elektronendichten konstruiert werden.
Aus dieser Startdichte werden das Hartree-Potential und das Austausch-Korrelationspotential berechnet. Zusammen mit dem externen Potential entsteht das effektive Kohn-Sham-Potential.
Lösung der Kohn-Sham-Gleichungen
Mit diesem Potential werden die Kohn-Sham-Gleichungen gelöst. Daraus erhält man neue Orbitale und Eigenwerte. Die besetzten Orbitale liefern eine aktualisierte Elektronendichte.
Symbolisch lässt sich der Ablauf schreiben als
\(n^{(k)}(\mathbf r) \rightarrow v_{\mathrm{eff}}^{(k)}(\mathbf r) \rightarrow \{\phi_i^{(k)}\} \rightarrow n^{(k+1)}(\mathbf r)\)
Selbstkonsistenz
Die neue Dichte wird mit der vorherigen Dichte verglichen. Stimmen beide noch nicht ausreichend überein, beginnt der Zyklus erneut. Dieser Prozess wird so lange wiederholt, bis eine festgelegte Konvergenzbedingung erfüllt ist.
Eine einfache Bedingung kann beispielsweise lauten
\(|E^{(k+1)}-E^{(k)}| < \delta\)
wobei \(\delta\) eine vorher festgelegte Toleranz ist.
Das Ergebnis wird als selbstkonsistente Lösung bezeichnet. Erst diese Lösung besitzt innerhalb der gewählten Näherung eine physikalisch konsistente Elektronendichte.
Berechenbare Eigenschaften
Aus einer konvergierten LDA-Berechnung können zahlreiche Größen abgeleitet werden. Dazu gehören Gesamtenergie, Elektronendichte, Kräfte auf Atomkerne, Gleichgewichtsgeometrien und elektronische Zustandsdichten.
Durch Variation der Atompositionen lassen sich energetisch bevorzugte Kristall- oder Molekülstrukturen bestimmen. Ebenso können Volumen, Gitterkonstanten und elastische Eigenschaften berechnet werden.
Anwendungsbereiche der LDA
Atome und Moleküle
Die LDA kann auf Atome und Moleküle angewendet werden. Allerdings weisen solche Systeme häufig stärker variierende Elektronendichten auf als kompakte Festkörper. Dadurch treten bestimmte systematische Fehler deutlicher hervor.
Insbesondere Bindungsenergien können überschätzt werden. Dennoch liefert die LDA häufig qualitativ sinnvolle Elektronendichten und stellt einen nützlichen Vergleichspunkt für komplexere Funktionale dar.
Festkörper
Eine ihrer größten Stärken zeigt die LDA bei Festkörpern mit vergleichsweise gleichmäßiger Elektronenverteilung. Bei vielen einfachen Metallen und dicht gepackten Materialien können strukturelle Eigenschaften erstaunlich genau berechnet werden.
Typische Anwendungen umfassen die Bestimmung von Gitterkonstanten, Volumina, Kohäsionsenergien, Kompressionsmodulen und Phasenstabilitäten.
Die Gesamtenergie kann für verschiedene Gitterparameter berechnet werden. Das Energieminimum liefert das theoretische Gleichgewichtsvolumen.
Halbleiter
Auch für Halbleiter wurde und wird die LDA intensiv eingesetzt. Sie liefert häufig vernünftige Kristallstrukturen und qualitative Bandstrukturen.
Ein bekanntes Problem ist jedoch die systematische Unterschätzung fundamentaler Bandlücken. Eine berechnete Bandstruktur kann daher qualitativ richtig aussehen, während die Größe der Bandlücke deutlich vom experimentellen Wert abweicht.
Metalle
Metalle gehören zu den Systemen, bei denen die physikalische Grundidee der LDA besonders plausibel ist. Durch delokalisierte Leitungselektronen kann die Elektronendichte über größere Bereiche relativ gleichmäßig sein.
Das reale System ähnelt damit lokal stärker dem homogenen Elektronengas, auf dem die LDA basiert. Entsprechend kann die Näherung bei einfachen Metallen teilweise sehr gute Ergebnisse liefern.
Oberflächen und Grenzflächen
An Oberflächen nimmt die Elektronendichte vom Materialinneren zum Vakuum stark ab. Damit liegt gerade dort eine Situation vor, bei der die Annahme einer langsam variierenden Dichte problematisch werden kann.
Trotzdem wurde die LDA umfangreich zur Untersuchung von Oberflächenenergien, Adsorption und Grenzflächen verwendet. Die Ergebnisse müssen jedoch abhängig vom untersuchten System sorgfältig bewertet werden.
Nanostrukturen
Nanopartikel, dünne Schichten und Quantendrähte besitzen einen besonders hohen Anteil an Oberflächen und Grenzflächen. Dadurch können lokale Dichteänderungen eine größere Rolle spielen als in makroskopischen Festkörpern.
Die LDA kann dennoch als erste Näherung oder Referenzberechnung verwendet werden. Für quantitative Aussagen werden häufig weiterentwickelte Funktionale oder zusätzliche Korrekturverfahren eingesetzt.
Bedeutung der LDA für die Quantentechnologie
Materialentwicklung
Quantentechnologische Bauelemente benötigen Materialien mit genau definierten elektronischen Eigenschaften. DFT-Berechnungen ermöglichen es, Materialkandidaten theoretisch zu untersuchen, bevor sie aufwendig synthetisiert werden.
Die LDA kann dabei grundlegende Informationen über Kristallstrukturen, elektronische Zustände und energetische Stabilitäten liefern. Auch wenn spätere Berechnungen häufig komplexere Funktionale verwenden, kann die LDA einen effizienten Ausgangspunkt darstellen.
Quantenmaterialien
Quantenmaterialien zeigen makroskopische Eigenschaften, die wesentlich durch quantenmechanische Vielteilcheneffekte bestimmt werden. Beispiele sind topologische Materialien, magnetische Systeme und bestimmte supraleitende Verbindungen.
Die LDA kann elektronische Bandstrukturen und Grundzustandseigenschaften solcher Systeme berechnen. Bei stark korrelierten Materialien reicht sie jedoch häufig nicht aus. In diesen Fällen bildet sie oft den Ausgangspunkt für weiterführende Methoden.
Halbleiterbasierte Quantensysteme
Viele Qubit-Konzepte basieren auf Halbleitern, Nanostrukturen oder Defektzuständen. Die elektronische Struktur des umgebenden Materials beeinflusst dabei das Verhalten der relevanten Quantenzustände.
LDA-Berechnungen können Informationen über Kristallstrukturen, lokale Ladungsverteilungen und Defektzustände bereitstellen. Für präzise Anregungsenergien oder stark lokalisierte Zustände müssen jedoch häufig genauere Verfahren eingesetzt werden.
Supraleitende Quantentechnologien
Supraleitende Schaltkreise gehören zu den wichtigsten Plattformen für Quantencomputer. Die mikroskopische Beschreibung der verwendeten Materialien beginnt häufig mit Berechnungen ihrer elektronischen Struktur.
Die LDA kann hier Grundzustandseigenschaften metallischer und intermetallischer Systeme beschreiben. Die eigentliche Supraleitung erfordert allerdings zusätzliche theoretische Konzepte, da das konventionelle LDA-Funktional nicht automatisch eine vollständige Beschreibung supraleitender Paarbildung liefert.
Defektzentren und Festkörper-Qubits
Defekte in Kristallen können diskrete elektronische Zustände erzeugen, die sich als Qubits nutzen lassen. Ihre Eigenschaften hängen von der lokalen chemischen Umgebung, der Ladungsverteilung und der elektronischen Struktur ab.
Hier stößt die LDA häufig an Grenzen, weil Defektzustände stark lokalisiert sein können. Trotzdem liefert sie eine grundlegende Beschreibung des Wirtskristalls und kann für erste Strukturrelaxationen oder Vergleichsrechnungen eingesetzt werden.
Stärken der Local Density Approximation
Einfache mathematische Struktur
Die LDA hängt ausschließlich von der lokalen Elektronendichte ab. Dadurch besitzt das Funktional eine vergleichsweise einfache mathematische Form und lässt sich effizient numerisch auswerten.
Moderate Rechenkosten
Im Vergleich zu komplexeren Austausch-Korrelationsfunktionalen ist der zusätzliche Rechenaufwand gering. Besonders bei großen Festkörpersystemen oder umfangreichen Parameterstudien kann dies relevant sein.
Gute strukturelle Ergebnisse
Für viele kompakte Festkörper liefert die LDA zuverlässige qualitative und teilweise sehr gute quantitative Ergebnisse. Häufig werden Gitterstrukturen, Bindungslängen und Schwingungseigenschaften zufriedenstellend beschrieben.
Gut untersuchte Fehlerstruktur
Da die LDA seit Jahrzehnten intensiv verwendet wird, sind ihre typischen Fehler gut bekannt. Dies ist ein praktischer Vorteil. Ein berechnetes Ergebnis kann besser eingeordnet werden, wenn bekannt ist, in welche Richtung ein Verfahren bei bestimmten Materialklassen systematisch abweicht.
Referenz für weiterentwickelte Funktionale
Die LDA bildet einen natürlichen Ausgangspunkt zum Verständnis moderner Austausch-Korrelationsfunktionale. Viele Verbesserungen können als Erweiterungen der lokalen Idee verstanden werden.
Grenzen und typische Fehler der LDA
Stark inhomogene Elektronendichten
Die wichtigste konzeptionelle Einschränkung folgt direkt aus der Grundannahme der LDA. Sie verwendet lediglich den lokalen Wert der Elektronendichte. Wenn sich die Dichte sehr schnell im Raum ändert, enthält dieser Einzelwert nicht genügend Information über die Umgebung.
Solche Situationen treten beispielsweise an Oberflächen, in Molekülen, bei schwachen Bindungen und in stark lokalisierten elektronischen Zuständen auf.
Überbindung
Die LDA zeigt häufig eine Tendenz zur Überbindung. Chemische Bindungen können stärker erscheinen als sie experimentell sind. Dadurch können Bindungslängen und Gitterkonstanten zu klein ausfallen.
Diese Tendenz ist nicht universell, aber charakteristisch genug, dass sie bei der Interpretation von LDA-Ergebnissen berücksichtigt werden sollte.
Bandlückenproblem
Eine besonders bekannte Schwäche betrifft Bandlücken von Halbleitern und Isolatoren. Kohn-Sham-Eigenwertdifferenzen sind nicht allgemein mit fundamentalen Anregungsenergien identisch.
Zusätzlich besitzen einfache lokale und semilokale Funktionale nicht die vollständige Ableitungsdiskontinuität des exakten Austausch-Korrelationsfunktionals.
Die fundamentale Bandlücke kann formal als
\(E_g = I - A\)
geschrieben werden, wobei \(I\) die Ionisierungsenergie und \(A\) die Elektronenaffinität bezeichnet.
Die einfache Differenz zwischen niedrigstem unbesetztem und höchstem besetztem Kohn-Sham-Zustand kann hiervon abweichen. Deshalb sind LDA-Bandlücken häufig deutlich zu klein.
Stark korrelierte Systeme
In Materialien mit stark lokalisierten \(d\)- oder \(f\)-Elektronen können lokale Coulomb-Wechselwirkungen besonders wichtig sein. Beispiele sind zahlreiche Übergangsmetalloxide und Verbindungen seltener Erden.
In solchen Fällen kann die LDA Elektronen zu stark delokalisieren. Dies kann zu einer falschen Beschreibung metallischer oder isolierender Zustände, magnetischer Eigenschaften und lokaler Momente führen.
Van-der-Waals-Wechselwirkungen
Langreichweitige Dispersionskräfte entstehen durch korrelierte Ladungsfluktuationen über größere Entfernungen. Ein rein lokales Funktional kann diese nicht systematisch korrekt erfassen.
Die LDA kann in manchen Fällen scheinbar brauchbare Bindungsenergien liefern, doch dies kann aus einer zufälligen Fehlerkompensation entstehen und stellt keine zuverlässige Beschreibung von Van-der-Waals-Wechselwirkungen dar.
Selbstwechselwirkungsfehler
In der Hartree-Energie wechselwirkt die Elektronendichte mathematisch auch mit sich selbst. Im exakten Austausch-Korrelationsfunktional wird diese unphysikalische Selbstwechselwirkung vollständig kompensiert.
Bei der LDA ist diese Kompensation unvollständig. Dies kann insbesondere bei lokalisierten Zuständen erhebliche Fehler verursachen.
Angeregte Zustände
Die gewöhnliche Grundzustands-DFT ist primär auf Grundzustandseigenschaften ausgerichtet. Kohn-Sham-Orbitale und Eigenwerte werden zwar häufig zur Interpretation elektronischer Strukturen verwendet, ihre direkte Gleichsetzung mit realen Anregungsenergien ist jedoch grundsätzlich eingeschränkt.
Für optische Spektren, elektronische Anregungen oder quasiteilchenartige Zustände werden deshalb häufig weiterführende Verfahren eingesetzt.
LDA im Vergleich zu weiterentwickelten Funktionalen
Generalized Gradient Approximation
Die Generalized Gradient Approximation, kurz GGA, erweitert die rein lokale Beschreibung. Zusätzlich zur Elektronendichte wird deren räumlicher Gradient berücksichtigt.
Das Funktional besitzt dann grundsätzlich die Form
\(E_{xc}^{\mathrm{GGA}}[n] = \int f(n(\mathbf r),\nabla n(\mathbf r))\,d^3r\)
Dadurch kann das Funktional unterscheiden, ob eine bestimmte Elektronendichte Teil eines nahezu homogenen Bereichs ist oder sich in einer Region mit starkem Dichtegradienten befindet.
Bekannte GGA-Funktionale werden heute in vielen Bereichen standardmäßig eingesetzt.
Meta-GGA
Meta-GGA-Funktionale verwenden zusätzliche lokale oder semilokale Größen. Dazu können etwa höhere Ableitungen der Dichte oder die kinetische Energiedichte gehören.
Die zusätzliche Information ermöglicht eine differenziertere Beschreibung verschiedener Bindungssituationen, ohne unmittelbar den Rechenaufwand nichtlokaler Hybridfunktionale zu erreichen.
Hybridfunktionale
Hybridfunktionale kombinieren DFT-Austausch mit einem Anteil exakten Hartree-Fock-Austauschs. Allgemein kann eine solche Konstruktion schematisch geschrieben werden als
\(E_{xc}^{\mathrm{hybrid}} = aE_x^{\mathrm{HF}} + (1-a)E_x^{\mathrm{DFT}} + E_c^{\mathrm{DFT}}\)
wobei \(a\) den Anteil des exakten Austauschs beschreibt.
Hybridfunktionale können unter anderem bei Molekülen, Halbleitern und Defektzuständen deutliche Verbesserungen liefern. Der Rechenaufwand ist jedoch erheblich höher als bei der LDA.
DFT+U
Für stark lokalisierte Elektronenzustände wird häufig eine zusätzliche Hubbard-Korrektur verwendet. Ein vereinfachter zusätzlicher Energiebeitrag kann beispielsweise in der Form
\(E_{\mathrm{DFT}+U} = E_{\mathrm{DFT}} + E_U\)
geschrieben werden.
Die DFT+U-Methode soll die übermäßige Delokalisierung bestimmter Orbitale reduzieren. Sie wird häufig bei Übergangsmetalloxiden und Materialien mit lokalisierten \(d\)- oder \(f\)-Elektronen verwendet.
Van-der-Waals-korrigierte Funktionale
Um dispersive Wechselwirkungen zu erfassen, wurden Funktionale und Korrekturverfahren entwickelt, die explizit nichtlokale Beiträge oder zusätzliche Dispersionskorrekturen enthalten.
Diese Erweiterungen sind für Molekülkristalle, Schichtmaterialien, Adsorption und viele biologische oder chemische Systeme besonders wichtig.
Einordnung der LDA
Die LDA kann als unterste Stufe einer Hierarchie von Dichtefunktionalen verstanden werden. Mit zunehmender methodischer Komplexität werden zusätzliche Informationen über das elektronische System berücksichtigt.
Mehr Komplexität bedeutet jedoch nicht automatisch bessere Ergebnisse für jede Eigenschaft. Das geeignete Funktional muss anhand des konkreten physikalischen Problems ausgewählt werden.
Praktische Implementierung
Basisdarstellungen
Zur numerischen Lösung der Kohn-Sham-Gleichungen müssen die Orbitale in einer geeigneten Basis dargestellt werden. In Festkörperrechnungen werden häufig ebene Wellen verwendet.
Eine Kohn-Sham-Wellenfunktion kann schematisch als Entwicklung
\(\phi_i(\mathbf r) = \sum_{\mathbf G} c_{i,\mathbf G}e^{i(\mathbf k+\mathbf G)\mathbf r}\)
dargestellt werden.
Alternativ können lokalisierte atomähnliche Funktionen oder Realraumgitter verwendet werden.
Pseudopotentiale
Innerhalb eines Atoms verändern sich kernnahe Wellenfunktionen sehr schnell. Eine direkte Darstellung aller Elektronen wäre insbesondere bei schweren Elementen rechnerisch aufwendig.
Pseudopotentiale ersetzen daher häufig den Einfluss der inneren Elektronen und des Atomkerns durch ein effektives Potential. Dadurch müssen hauptsächlich die chemisch relevanten Valenzelektronen explizit berechnet werden.
Typische Software
Die LDA ist in nahezu allen bedeutenden DFT-Programmen implementiert. Zu den bekannten Programmpaketen gehören Quantum ESPRESSO, VASP, ABINIT, GPAW, CASTEP und CP2K.
Diese Programme unterscheiden sich unter anderem hinsichtlich ihrer Basisdarstellungen, Pseudopotentialverfahren und numerischen Algorithmen. Das zugrunde liegende LDA-Konzept bleibt jedoch identisch.
Numerische Konvergenz
Eine physikalisch sinnvolle Rechnung muss numerisch konvergiert sein. Bei ebenen Wellen ist beispielsweise der Energie-Cutoff entscheidend.
Es werden nur Basisfunktionen berücksichtigt, deren kinetische Energie eine gewählte Grenze nicht überschreitet:
\(\frac{|\mathbf k+\mathbf G|^2}{2} \leq E_{\mathrm{cut}}\)
Auch die Anzahl der verwendeten \(k\)-Punkte im Brillouin-Raum beeinflusst das Ergebnis. Insbesondere bei Metallen kann eine zu grobe Abtastung erhebliche Fehler verursachen.
Reproduzierbarkeit
Die Bezeichnung LDA allein reicht nicht aus, um eine Berechnung vollständig zu dokumentieren. Angegeben werden sollten die genaue Parametrisierung des Austausch-Korrelationsfunktionals, verwendete Pseudopotentiale, Basisparameter, \(k\)-Punkt-Netz und Konvergenzkriterien.
Nur bei ausreichender Dokumentation können Ergebnisse reproduziert und sinnvoll mit anderen Rechnungen verglichen werden.
LDA und moderne quantenmechanische Simulationsmethoden
LDA als Referenzmethode
Auch wenn moderne Berechnungen häufig komplexere Funktionale verwenden, bleibt die LDA eine wichtige Referenz. Sie ist mathematisch einfach, effizient und umfangreich getestet.
Bei der Entwicklung neuer Methoden ist es daher üblich, Ergebnisse zunächst mit LDA oder anderen etablierten Standardfunktionalen zu vergleichen.
GW-Verfahren
Für präzisere quasiteilchenartige Anregungsenergien können GW-Verfahren eingesetzt werden. Dabei wird die Selbstenergie der Elektronen mit Hilfe der Einteilchen-Greenschen Funktion und einer abgeschirmten Coulomb-Wechselwirkung berechnet.
Schematisch wird die Selbstenergie häufig als
\(\Sigma = iGW\)
geschrieben.
Eine LDA-Berechnung kann dabei als Ausgangspunkt dienen, auf dessen Orbitalen und Energien die weiterführende Vielteilchenrechnung aufbaut.
Dynamical Mean-Field Theory
Für stark korrelierte Materialien kann die Dichtefunktionaltheorie mit der Dynamical Mean-Field Theory kombiniert werden. Verfahren wie LDA+DMFT verbinden eine realistische elektronische Bandstruktur mit einer expliziteren Behandlung lokaler dynamischer Korrelationen.
Solche Methoden werden unter anderem bei Übergangsmetallverbindungen und Materialien mit stark korrelierten Elektronenzuständen eingesetzt.
Machine Learning und Materialdaten
DFT-Berechnungen erzeugen große Mengen physikalisch strukturierter Daten. Diese Daten können zum Training maschineller Lernverfahren genutzt werden.
Maschinelle Lernmodelle können beispielsweise Energien, Kräfte oder Materialeigenschaften aus atomaren Strukturen vorhersagen. Die Qualität solcher Modelle hängt wesentlich von der Qualität und Konsistenz der Trainingsdaten ab.
Da die LDA effizient ist, eignet sie sich grundsätzlich auch zur Erzeugung großer homogener Datensätze. Für Anwendungen mit hohen Genauigkeitsanforderungen werden jedoch häufig modernere Funktionale gewählt.
High-Throughput-Berechnungen
Bei High-Throughput-Verfahren werden tausende oder sogar deutlich mehr Materialstrukturen automatisiert berechnet. Ziel ist es, Trends und geeignete Materialkandidaten zu identifizieren.
In solchen Szenarien ist das Verhältnis zwischen Rechenaufwand und Genauigkeit entscheidend. Ein einfaches Funktional kann dann trotz bekannter Einschränkungen wertvoll sein, sofern seine Fehler für die untersuchte Eigenschaft ausreichend verstanden sind.
Fachliche Einordnung: Wann ist die LDA geeignet?
Geeignete Systeme
Die LDA ist besonders geeignet, wenn die Elektronendichte räumlich relativ gleichmäßig ist. Dies gilt näherungsweise für viele einfache Metalle und dicht gepackte Festkörper.
Auch für strukturelle Optimierungen kann die Methode sinnvoll sein, sofern die bekannte Tendenz zur Überbindung berücksichtigt wird.
Problematische Systeme
Vorsicht ist erforderlich bei stark lokalisierten Elektronen, stark korrelierten Materialien, schwach gebundenen Molekülen, ausgeprägten Van-der-Waals-Systemen und elektronischen Anregungen.
Auch bei Oberflächen und Defekten sollte geprüft werden, ob die lokale Näherung für die gewünschte physikalische Größe ausreichend ist.
Funktionalwahl als physikalische Entscheidung
Es gibt kein Austausch-Korrelationsfunktional, das für jede Materialklasse und jede beobachtbare Größe automatisch die beste Wahl darstellt.
Die Auswahl sollte anhand der relevanten Physik erfolgen. Für strukturelle Eigenschaften eines einfachen Metalls kann die LDA zweckmäßig sein. Für eine genaue Halbleiter-Bandlücke kann dagegen ein Hybridfunktional oder ein GW-Verfahren geeigneter sein.
Für stark lokalisierte Zustände können DFT+U, Hybridfunktionale oder Vielteilchenmethoden erforderlich sein.
Vergleich mit Experimenten und Benchmarks
Berechnete Werte sollten nicht allein deshalb als korrekt angesehen werden, weil eine numerische Lösung konvergiert ist. Numerische Konvergenz bedeutet lediglich, dass die gewählten Gleichungen innerhalb der angegebenen Einstellungen stabil gelöst wurden.
Die physikalische Genauigkeit hängt zusätzlich vom Funktional, der Basis, den verwendeten Potentialen und dem Modell des Systems ab.
Wo möglich sollten Ergebnisse deshalb mit experimentellen Daten oder höherwertigen theoretischen Verfahren verglichen werden.
Fazit
Die Local Density Approximation ist eine der fundamentalsten Näherungen der Dichtefunktionaltheorie. Ihre zentrale Idee besteht darin, die Austausch- und Korrelationsenergie eines inhomogenen elektronischen Systems lokal durch die bekannten Eigenschaften eines homogenen Elektronengases gleicher Dichte zu approximieren.
Mathematisch führt dies zu einem besonders einfachen Funktional:
\(E_{xc}^{\mathrm{LDA}}[n] = \int n(\mathbf r)\epsilon_{xc}(n(\mathbf r))\,d^3r\)
Diese kompakte Formel verdeckt, wie wirkungsvoll die zugrunde liegende Idee ist. Ein hochkomplexes quantenmechanisches Vielteilchenproblem wird auf eine Berechnung zurückgeführt, deren zentrale Variable lediglich die räumliche Elektronendichte ist.
Die LDA besitzt klare Stärken. Sie ist rechnerisch effizient, mathematisch gut verstanden und liefert insbesondere für Systeme mit langsam variierender Elektronendichte häufig brauchbare oder sogar sehr gute Ergebnisse. Viele grundlegende Entwicklungen der rechnergestützten Festkörperphysik wären ohne diese Näherung kaum in derselben Form möglich gewesen.
Gleichzeitig sind ihre Grenzen eindeutig. Die rein lokale Behandlung kann starke Dichteinhomogenitäten, langreichweitige Korrelationen, dispersive Wechselwirkungen und stark lokalisierte elektronische Zustände nicht vollständig erfassen. Bandlücken von Halbleitern und Isolatoren werden häufig unterschätzt, stark korrelierte Materialien können qualitativ falsch beschrieben werden und Van-der-Waals-Wechselwirkungen erfordern in der Regel zusätzliche theoretische Ansätze.
Gerade diese klare Kombination aus Leistungsfähigkeit und bekannten Grenzen macht die LDA wissenschaftlich wertvoll. Sie ist weder als universell exakte Methode noch als überholte historische Näherung zu betrachten. Vielmehr stellt sie einen elementaren Baustein der modernen elektronischen Strukturtheorie dar.
Für die Quantentechnologie liegt ihre Bedeutung vor allem in der Materialsimulation. Sie ermöglicht die Untersuchung elektronischer Grundzustände, Kristallstrukturen, Ladungsverteilungen und energetischer Stabilitäten. In komplexeren Fällen bildet sie einen Ausgangspunkt für GGA, Hybridfunktionale, DFT+U, GW oder dynamische Vielteilchenverfahren.
Das Verständnis der Local Density Approximation vermittelt damit weit mehr als die Kenntnis eines einzelnen Funktionals. Es zeigt exemplarisch, wie physikalische Modellbildung funktioniert: Ein exakt unlösbares Problem wird auf einen kontrollierten, nachvollziehbaren und praktisch nutzbaren Näherungsansatz reduziert. Genau darin liegt die bleibende Bedeutung der LDA für Quantenphysik, Materialwissenschaft und Quantentechnologie.
Mit freundlichen Grüßen
Anhang
Wissenschaftliche Zeitschriften und Artikel
Die wissenschaftliche Grundlage der Local Density Approximation liegt in der Entwicklung der Dichtefunktionaltheorie, der Theorie des homogenen Elektronengases und der Konstruktion praktisch nutzbarer Austausch-Korrelationsfunktionale. Die folgenden Arbeiten bilden eine belastbare Quellenbasis für eine wissenschaftliche Abhandlung zur LDA. Dabei wird zwischen fundamentaler Primärliteratur, Arbeiten zum Elektronengas und weiterführender Literatur zu Grenzen und Verbesserungen der lokalen Näherung unterschieden.
Grundlegende Primärliteratur zur Dichtefunktionaltheorie und LDA
- P. Hohenberg, W. Kohn: Inhomogeneous Electron Gas, Physical Review, 1964.
- Fundamentale Primärliteratur der Dichtefunktionaltheorie. Die Arbeit enthält die später nach Hohenberg und Kohn benannten Theoreme und zeigt, dass die Grundzustandsenergie eines wechselwirkenden Elektronensystems grundsätzlich als Funktional der Elektronendichte formuliert werden kann. Für eine Abhandlung zur LDA ist diese Veröffentlichung die theoretische Ausgangsbasis für die Ablösung der Vielteilchen-Wellenfunktion durch die Elektronendichte als zentrale Variable.
- W. Kohn, L. J. Sham: Self-Consistent Equations Including Exchange and Correlation Effects, Physical Review, 1965.
- Diese Primärarbeit führt den Kohn-Sham-Formalismus ein und ist für das Verständnis der LDA unverzichtbar. Sie zeigt, wie das wechselwirkende Vielteilchenproblem auf selbstkonsistent zu lösende Einteilchengleichungen zurückgeführt werden kann. Besonders relevant ist die Einführung eines Austausch-Korrelationspotentials, in das die Local Density Approximation unmittelbar eingeht.
- P. A. M. Dirac: Note on Exchange Phenomena in the Thomas Atom, Mathematical Proceedings of the Cambridge Philosophical Society, 1930.
- Historische und theoretische Grundlagenliteratur zur lokalen Beschreibung der Austauschenergie. Diracs Behandlung des Austauschs in einem Elektronengas liefert die Grundlage für den bis heute in der LDA verwendeten lokalen Austauschterm. Die Arbeit eignet sich insbesondere für den Abschnitt über die Herkunft und physikalische Interpretation des LDA-Austauschfunktionals.
Homogenes Elektronengas und Korrelationsenergie
- D. M. Ceperley, B. J. Alder: Ground State of the Electron Gas by a Stochastic Method, Physical Review Letters, 1980.
- Zentrale Spezial- und Primärliteratur für die Korrelationskomponente moderner LDA-Funktionale. Ceperley und Alder berechneten die Grundzustandsenergie des homogenen Elektronengases mit stochastischen Vielteilchenmethoden. Die daraus gewonnenen Daten wurden zur Referenz für mehrere später verbreitete LDA- und LSDA-Parametrisierungen. Die Arbeit ist besonders wichtig, wenn erklärt werden soll, woher die numerischen Korrelationsenergien der LDA stammen.
- S. H. Vosko, L. Wilk, M. Nusair: Accurate Spin-Dependent Electron Liquid Correlation Energies for Local Spin Density Calculations: A Critical Analysis, Canadian Journal of Physics, 1980.
- Eine der klassischen Arbeiten zur praktischen Parametrisierung der Korrelationsenergie des spinpolarisierten homogenen Elektronengases. Das sogenannte VWN-Funktional ist für Local-Spin-Density-Berechnungen von großer Bedeutung. Die Quelle eignet sich zur Erläuterung des Übergangs von numerischen Elektronengasdaten zu analytisch auswertbaren Funktionalen.
- J. P. Perdew, A. Zunger: Self-Interaction Correction to Density-Functional Approximations for Many-Electron Systems, Physical Review B, 1981.
- Diese Arbeit ist sowohl für die praktische LDA als auch für die Analyse ihrer Grenzen von großer Bedeutung. Sie enthält eine verbreitete Parametrisierung der Ceperley-Alder-Daten und behandelt zugleich den Selbstwechselwirkungsfehler approximativer Dichtefunktionale. Die Quelle kann deshalb sowohl für die Beschreibung der Perdew-Zunger-LDA als auch für den Abschnitt über typische Fehler lokaler Funktionale genutzt werden.
- J. P. Perdew, Y. Wang: Accurate and Simple Analytic Representation of the Electron-Gas Correlation Energy, Physical Review B, 1992.
- Spezialisierte Literatur zur analytischen Darstellung der Korrelationsenergie des homogenen Elektronengases. Die Arbeit entwickelt eine präzise Parametrisierung in Abhängigkeit von Elektronendichte und Spinpolarisation und vergleicht sie mit früheren Ansätzen wie VWN und Perdew-Zunger. Sie eignet sich besonders für eine vertiefte Darstellung verschiedener LDA- beziehungsweise LSDA-Korrelationsfunktionale.
Einordnung, Grenzen und Weiterentwicklungen
- R. O. Jones, O. Gunnarsson: The Density Functional Formalism, Its Applications and Prospects, Reviews of Modern Physics, 1989.
- Umfassende Übersichtsarbeit zur Dichtefunktionaltheorie und ihrer praktischen Anwendung. Die Autoren behandeln ausführlich die Leistungsfähigkeit der Local Density Approximation und ordnen ihre überraschend guten Ergebnisse für zahlreiche Festkörpersysteme in den theoretischen Gesamtkontext ein. Die Arbeit eignet sich als Hintergrund- und Überblicksquelle für Bewertung, Anwendung und historische Einordnung der LDA.
- J. P. Perdew, K. Burke, M. Ernzerhof: Generalized Gradient Approximation Made Simple, Physical Review Letters, 1996.
- Diese Primärarbeit ist für die Abgrenzung zwischen LDA und Generalized Gradient Approximation besonders wichtig. Das PBE-Funktional berücksichtigt zusätzlich zur lokalen Dichte deren räumliche Variation und gehört zu den einflussreichsten Weiterentwicklungen semilokaler Dichtefunktionale. Die Quelle kann genutzt werden, um präzise zu erklären, welche zusätzliche Information eine GGA gegenüber der LDA verwendet.
Bücher und Monographien
Bücher und Monographien sind besonders geeignet, um die LDA nicht nur als einzelne Näherungsformel, sondern innerhalb der gesamten elektronischen Strukturtheorie zu verstehen. Die folgenden Werke decken sowohl die mathematischen Grundlagen der Dichtefunktionaltheorie als auch Austausch-Korrelationsfunktionale, numerische Verfahren und praktische Festkörperberechnungen ab.
Grundlagenwerke zur Dichtefunktionaltheorie
- Reiner M. Dreizler, Eberhard K. U. Gross: Density Functional Theory – An Approach to the Quantum Many-Body Problem, Springer, 1990.
- Eine grundlegende Monographie zur mathematischen und physikalischen Struktur der Dichtefunktionaltheorie. Das Werk behandelt den Kohn-Sham-Formalismus, Austausch- und Energiefunktionale sowie ausdrücklich die Local Density Approximation und ihre Erweiterungen. Für eine wissenschaftliche Abhandlung eignet es sich besonders als Grundlagenquelle für Herleitungen und theoretische Zusammenhänge.
- Eberhard Engel, Reiner M. Dreizler: Density Functional Theory – An Advanced Course, Springer, 2011.
- Fortgeschrittenes Standardwerk zur Dichtefunktionaltheorie. Besonders relevant sind die ausführlichen Kapitel zu Existenztheoremen, effektiven Einteilchengleichungen und Austausch-Korrelationsfunktionalen. Das Buch eignet sich, wenn die LDA in einen modernen theoretischen Rahmen eingeordnet und mit orbitalabhängigen, zeitabhängigen oder weiterentwickelten Funktionalen verglichen werden soll.
Elektronische Struktur und Festkörperberechnungen
- Richard M. Martin: Electronic Structure – Basic Theory and Practical Methods, 2. Auflage, Cambridge University Press, 2020.
- Umfassendes Standardwerk zur elektronischen Struktur von Materialien. Das Buch verbindet Dichtefunktionaltheorie mit dem homogenen Elektronengas, Austausch-Korrelationsfunktionalen, Bandstrukturen, Pseudopotentialen, ebenen Wellen und numerischen Berechnungsverfahren. Für eine Abhandlung mit Bezug zur Quantentechnologie ist es besonders wertvoll, weil es die LDA unmittelbar mit realen Festkörper- und Materialsimulationen verbindet.
Quantenchemischer Zugang zur LDA
- Wolfram Koch, Max C. Holthausen: A Chemist's Guide to Density Functional Theory, 2. Auflage, Wiley-VCH, 2001.
- Praxisorientiertes Grundlagenwerk mit ausführlicher Darstellung von Elektronendichte, Hohenberg-Kohn-Theoremen, Kohn-Sham-Methode und approximativen Austausch-Korrelationsfunktionalen. Besonders nützlich ist die Behandlung der Local Density und Local Spin Density Approximation sowie ihrer Abgrenzung gegenüber GGA- und Hybridfunktionalen. Das Werk eignet sich für die verständliche physikalische Interpretation der LDA und für den Vergleich unterschiedlicher Funktionalklassen.
Online-Ressourcen und Datenbanken
Online-Ressourcen sind insbesondere für die praktische Arbeit mit LDA-Funktionalen, die Überprüfung konkreter Implementierungen und die Recherche berechneter Materialdaten wichtig. Bevorzugt werden hier offizielle Dokumentationen, etablierte wissenschaftliche Funktionalbibliotheken und Forschungsdatenplattformen verwendet. Sie ersetzen die Primärliteratur nicht, ergänzen sie jedoch um technische und reproduzierbare Informationen.
Funktionalbibliotheken und Referenzdaten
- Libxc: Library of Exchange-Correlation Functionals for Density-Functional Theory.
- Libxc ist eine wissenschaftliche Bibliothek für Austausch-Korrelationsfunktionale und wird von zahlreichen elektronischen Strukturprogrammen verwendet. Die Funktionsübersicht enthält verschiedene LDA-, GGA-, Meta-GGA- und Hybridfunktionale einschließlich Literaturreferenzen und eindeutiger Funktionskennungen. Für eine Abhandlung eignet sich Libxc besonders zur Überprüfung, welche konkrete Variante eines LDA-Funktionals in einer Rechnung gemeint ist.
DFT-Software und technische Dokumentationen
- Quantum ESPRESSO: Offizielle Dokumentation der Exchange-Correlation-Funktionale.
- Quantum ESPRESSO gehört zu den etablierten Open-Source-Programmpaketen für elektronische Strukturberechnungen mit ebenen Wellen und Dichtefunktionaltheorie. Die Dokumentation zeigt, wie Austausch-Korrelationsfunktionale einschließlich LDA-Parametrisierungen ausgewählt und in praktischen Rechnungen eingesetzt werden. Sie ist vor allem als technische Recherchehilfe für Implementierung, Reproduzierbarkeit und konkrete Berechnungsparameter geeignet.
- GPAW: Exchange-Correlation Functionals Documentation.
- Die GPAW-Dokumentation beschreibt die praktische Verwendung verschiedener Austausch-Korrelationsfunktionale und weist die Local Density Approximation ausdrücklich als implementierte Funktionalklasse aus. Die Ressource ist geeignet, um die theoretische Darstellung der LDA mit einer konkreten modernen DFT-Implementierung zu verbinden und Unterschiede zwischen LDA, GGA, Meta-GGA und Hybridfunktionalen nachzuvollziehen.
- ABINIT: Exchange-Correlation Functionals Documentation.
- ABINIT ist ein etabliertes Open-Source-Paket für First-Principles-Berechnungen von Materialien. Die Dokumentation zu Austausch und Korrelation enthält Hinweise zu verfügbaren Funktionalen sowie zu LDA- und GGA-Berechnungen. Sie eignet sich als technische Ergänzung für Abschnitte über die numerische Implementierung und den praktischen Vergleich unterschiedlicher DFT-Näherungen.
Materialdatenbanken und Forschungsplattformen
- NOMAD: Materials Science Data Platform and Repository.
- NOMAD ist eine Forschungsdatenplattform für Materialwissenschaften und besitzt historische Wurzeln als Repository für DFT-Berechnungen. Die Plattform speichert strukturierte Informationen zu Berechnungsmethoden, Materialien und numerischen Ergebnissen. Für eine wissenschaftliche Abhandlung eignet sie sich als Recherchehilfe für veröffentlichte elektronische Strukturrechnungen und zur Untersuchung, welche DFT-Methoden und Funktionale bei konkreten Materialien verwendet wurden.
- Materials Project: Materials Methodology and Computational Materials Database.
- Das Materials Project stellt umfangreiche berechnete Materialdaten sowie eine Dokumentation der zugrunde liegenden elektronischen Strukturmethoden bereit. Die Plattform ist keine Primärquelle für die theoretische Herleitung der LDA, eignet sich aber hervorragend als Recherchehilfe für den Zusammenhang zwischen DFT, Bandstrukturen, Kristallstrukturen, thermodynamischer Stabilität und quantentechnologisch relevanten Materialeigenschaften.
Fachjournale und Verlage
- American Physical Society: Physical Review Journals.
- Die Zeitschriften der American Physical Society enthalten einen großen Teil der grundlegenden und weiterführenden Literatur zur Dichtefunktionaltheorie, zum homogenen Elektronengas, zur LDA und zu modernen Austausch-Korrelationsfunktionalen. Physical Review, Physical Review Letters, Physical Review B und Reviews of Modern Physics sind für die vertiefende Recherche besonders relevant.
- Springer Nature: Fachliteratur zu Dichtefunktionaltheorie und elektronischer Struktur.
- Springer Nature stellt zahlreiche Monographien und Fachbeiträge zur Dichtefunktionaltheorie, Vielteilchentheorie, Festkörperphysik und Computational Materials Science bereit. Die Plattform ist insbesondere für die Recherche vertiefender theoretischer Literatur und wissenschaftlicher Monographien geeignet.
- Cambridge University Press: Electronic Structure and Condensed Matter Physics.
- Cambridge University Press bietet wissenschaftliche Monographien und Fachliteratur zur elektronischen Struktur, Festkörperphysik und Vielteilchentheorie. Für die LDA sind insbesondere Werke relevant, die das homogene Elektronengas, Kohn-Sham-DFT, Austausch-Korrelationsfunktionale und numerische Materialberechnung gemeinsam behandeln.
Antiproton
Ein Antiproton ist das Antiteilchen des Protons, also ein Teilchen mit derselben Masse wie das Proton, jedoch entgegengesetzter elektrischer Ladung....
