Das Einband-Fermi-Hubbard-Modell gehört zu den grundlegenden Modellen der modernen Vielteilchenphysik. Es beschreibt Fermionen, typischerweise Elektronen oder ultrakalte fermionische Atome, die sich auf den Plätzen eines Gitters bewegen und lokal miteinander wechselwirken. Trotz seiner vergleichsweise einfachen mathematischen Struktur erzeugt das Modell eine außerordentlich vielfältige Physik. Dazu gehören Mott-Isolatoren, Antiferromagnetismus, starke Quantenkorrelationen, ungewöhnliche metallische Zustände und möglicherweise Mechanismen, die mit unkonventioneller Supraleitung zusammenhängen.
Die zentrale Idee ist einfach: Fermionen können zwischen benachbarten Gitterplätzen tunneln, während zwei Fermionen mit unterschiedlichem Spin auf demselben Gitterplatz miteinander wechselwirken. Diese beiden Prozesse konkurrieren miteinander. Das Hopping begünstigt delokalisierte Zustände, während eine starke repulsive Wechselwirkung die Teilchen lokalisiert. Ein großer Teil der Physik des Modells entsteht allein aus diesem Wettbewerb.
Das Einband-Fermi-Hubbard-Modell ist keine vollständige Beschreibung eines realen Festkörpers. Es reduziert die elektronische Struktur bewusst auf ein einziges relevantes Energieband und konzentriert sich auf die wichtigsten lokalen Korrelationsprozesse. Gerade diese Reduktion macht es wertvoll. Sie erlaubt es, fundamentale Mechanismen stark korrelierter Quantensysteme zu untersuchen, ohne jedes mikroskopische Detail eines konkreten Materials berücksichtigen zu müssen.
In der Quantentechnologie besitzt das Modell eine zusätzliche Bedeutung. Ultrakalte Atome in optischen Gittern können den Hubbard-Hamiltonoperator mit hoher Kontrolle realisieren. Gleichzeitig ist das Modell zu einem wichtigen Testproblem für digitale Quantencomputer geworden. Damit verbindet es Festkörperphysik, Atomphysik, Quanteninformation und numerische Vielteilchentheorie.
Die entscheidende Schwierigkeit liegt darin, dass sich viele relevante Zustände des Fermi-Hubbard-Modells nicht als unabhängige Teilchen beschreiben lassen. Die Freiheitsgrade sind quantenmechanisch korreliert, und die Größe des zu behandelnden Hilbertraums wächst exponentiell mit der Anzahl der Gitterplätze. Das Einband-Fermi-Hubbard-Modell ist deshalb zugleich eines der einfachsten und eines der anspruchsvollsten Standardmodelle der Quantenphysik.
Physikalische Grundlagen
Fermionen und Pauli-Prinzip
Fermionen sind Teilchen mit halbzahligem Spin. Elektronen, Protonen, Neutronen sowie Atome wie Lithium-6 und Kalium-40 können unter geeigneten Bedingungen als Fermionen behandelt werden. Ihr wesentliches quantenstatistisches Merkmal ist das Pauli-Ausschlussprinzip: Zwei identische Fermionen können nicht denselben vollständigen Quantenzustand besetzen.
Im Fermi-Hubbard-Modell werden meistens zwei interne Zustände betrachtet, die vereinfachend als Spin-up und Spin-down bezeichnet werden. Ein einzelner Gitterplatz kann deshalb leer sein, ein Fermion mit Spin-up enthalten, ein Fermion mit Spin-down enthalten oder von zwei Fermionen mit entgegengesetztem Spin besetzt sein. Zwei Fermionen mit demselben Spin können aufgrund des Pauli-Prinzips nicht denselben lokalen Einteilchenzustand besetzen.
Fermionen in periodischen Gittern
In einem Kristall bewegen sich Elektronen im periodischen Potential der Atomrümpfe. Die Einteilchenzustände bilden Energiebänder. Für viele Fragestellungen ist nicht die gesamte Bandstruktur relevant. Befindet sich die Fermi-Energie in der Nähe eines einzelnen, relativ schmalen Bandes, kann die Beschreibung auf dieses Band beschränkt werden.
Das Einband-Fermi-Hubbard-Modell macht genau diese Annahme. Alle höheren und tieferen Bänder werden als nicht unmittelbar relevant betrachtet. Übrig bleibt ein effektiver Freiheitsgrad pro Spin und Gitterplatz. Die Komplexität der tatsächlichen elektronischen Wellenfunktionen geht dann in wenige effektive Parameter ein.
Wannier-Funktionen und lokalisierte Zustände
Für die Hubbard-Beschreibung ist eine Darstellung durch Wannier-Funktionen besonders geeignet. Während Bloch-Zustände über den gesamten Kristall ausgedehnt sind, lassen sich Wannier-Funktionen einzelnen Gitterplätzen zuordnen. Sie bilden eine lokalisierte Basis für das betrachtete Energieband.
In dieser Darstellung erhält das Modell seine anschauliche Form. Ein Fermion befindet sich überwiegend an einem bestimmten Gitterplatz, kann aber aufgrund der Überlappung benachbarter Wannier-Funktionen tunneln. Treffen Fermionen mit unterschiedlichem Spin auf demselben Platz zusammen, wirkt eine lokale Wechselwirkungsenergie.
Zweite Quantisierung
Das Modell wird zweckmäßig in zweiter Quantisierung formuliert. Der Operator \(c_{i\sigma}^{\dagger}\) erzeugt ein Fermion mit Spin \(\sigma\) auf dem Gitterplatz \(i\). Der Operator \(c_{i\sigma}\) vernichtet ein entsprechendes Fermion. Der Besetzungszahloperator lautet
\(n_{i\sigma} = c_{i\sigma}^{\dagger}c_{i\sigma}\)
und kann wegen des Pauli-Prinzips für einen gegebenen Spin nur die Werte null oder eins annehmen.
Die fermionischen Erzeugungs- und Vernichtungsoperatoren erfüllen Antikommutatorrelationen:
\(\{c_{i\sigma},c_{j\sigma'}^{\dagger}\} = \delta_{ij}\delta_{\sigma\sigma'}\)
sowie
\(\{c_{i\sigma},c_{j\sigma'}\} = 0\)
Diese Relationen kodieren die Fermi-Statistik und sind für die gesamte Struktur des Modells entscheidend.
Der Hamiltonoperator des Einband-Fermi-Hubbard-Modells
Grundform des Modells
Der übliche Hamiltonoperator des Einband-Fermi-Hubbard-Modells lautet
\(H = -t\sum_{\langle i,j\rangle,\sigma}\left(c_{i\sigma}^{\dagger}c_{j\sigma} + c_{j\sigma}^{\dagger}c_{i\sigma}\right) + U\sum_i n_{i\uparrow}n_{i\downarrow} - \mu\sum_{i,\sigma}n_{i\sigma}\)
Die drei Bestandteile besitzen klar unterscheidbare Aufgaben. Der erste Term beschreibt die Bewegung der Fermionen. Der zweite Term beschreibt die lokale Wechselwirkung. Der dritte Term kontrolliert über das chemische Potential die mittlere Teilchenzahl.
Hopping-Term
Der Parameter \(t\) ist die Hopping- oder Tunnelamplitude. Die Summe \(\langle i,j\rangle\) läuft in der einfachsten Variante über benachbarte Gitterplätze. Der Ausdruck
\(c_{i\sigma}^{\dagger}c_{j\sigma}\)
beschreibt einen Prozess, bei dem ein Fermion vom Platz \(j\) zum Platz \(i\) wechselt. Der hermitesch konjugierte Prozess berücksichtigt die umgekehrte Richtung.
Ein großes \(t\) begünstigt räumlich ausgedehnte Quantenzustände. Die Fermionen gewinnen kinetische Energie, wenn sie sich über viele Gitterplätze verteilen können. Ohne Wechselwirkung entsteht im Translations-invarianten Gitter eine gewöhnliche Bandstruktur.
Lokale Wechselwirkung
Der Wechselwirkungsterm lautet
\(U\sum_i n_{i\uparrow}n_{i\downarrow}\)
Nur wenn sich auf einem Gitterplatz gleichzeitig ein Fermion mit Spin-up und eines mit Spin-down befinden, trägt dieser Term Energie bei. Für \(U>0\) ist die Wechselwirkung repulsiv. Eine Doppelbesetzung kostet dann Energie. Für \(U<0\) ist die Wechselwirkung attraktiv und begünstigt die Paarbildung.
Das Standardproblem stark korrelierter Elektronensysteme wird häufig für \(U>0\) untersucht. Dort konkurriert die energetisch günstige Delokalisierung durch \(t\) mit der Vermeidung von Doppelbesetzung durch \(U\).
Chemisches Potential und Teilchenzahl
Der Term
\(-\mu\sum_{i,\sigma}n_{i\sigma}\)
ermöglicht die Beschreibung des Systems im großkanonischen Ensemble. Durch Variation des chemischen Potentials \(\mu\) kann die mittlere Teilchendichte eingestellt werden.
Besondere Bedeutung besitzt die Halbbesetzung. Bei einem System mit \(N\) Gitterplätzen entspricht sie im üblichen spinhalbigen Hubbard-Modell im Mittel einer Gesamtteilchenzahl
\(N_{\mathrm{F}} = N\)
also durchschnittlich einem Fermion pro Gitterplatz.
Herleitung aus einem kontinuierlichen Vielteilchenmodell
Das Hubbard-Modell kann als effektive Niedrigenergiebeschreibung eines kontinuierlichen Systems verstanden werden. Ausgangspunkt ist ein Feldoperator, der in Wannier-Funktionen des relevanten Bandes entwickelt wird. Werden nur das niedrigste beziehungsweise ein ausgewähltes Band, das Hopping zwischen nahen Gitterplätzen und die dominante lokale Wechselwirkung beibehalten, entsteht der Hubbard-Hamiltonoperator.
Der Hoppingparameter ergibt sich aus Matrixelementen des Einteilchen-Hamiltonoperators zwischen Wannier-Zuständen benachbarter Plätze. Schematisch gilt
\(t = -\int d^3r\, w_i^*(\mathbf{r}) H_0 w_j(\mathbf{r})\)
Die lokale Wechselwirkung wird entsprechend durch den räumlichen Überlapp der Wannier-Funktion mit dem Wechselwirkungspotential bestimmt. Bei einer kurzreichweitigen Kontaktwechselwirkung besitzt sie die Form
\(U = g\int d^3r\, |w_i(\mathbf{r})|^4\)
wobei \(g\) die Stärke der effektiven Zweiteilchenwechselwirkung beschreibt.
Die entscheidenden Modellparameter
Das Verhältnis von Wechselwirkung und Hopping
Für die qualitative Physik ist weniger der absolute Wert von \(U\) oder \(t\) entscheidend als ihr Verhältnis. Die dimensionslose Größe
\(\frac{U}{t}\)
misst, wie stark die Wechselwirkung im Vergleich zur Beweglichkeit der Fermionen ist.
Für
\(\frac{U}{t} \ll 1\)
dominiert die kinetische Energie. Das System lässt sich häufig von einem schwach wechselwirkenden Fermigas ausgehend behandeln. Für
\(\frac{U}{t} \gg 1\)
dominiert die lokale Wechselwirkung. Doppelbesetzungen werden bei repulsivem \(U\) stark unterdrückt, und kollektive Spinprozesse gewinnen an Bedeutung.
Füllung
Die Teilchendichte ist ebenso wichtig wie die Wechselwirkungsstärke. An Halbbesetzung führt eine starke Repulsion zu besonders ausgeprägten Korrelationsphänomenen. Werden Teilchen entfernt oder hinzugefügt, spricht man von Dotierung. Dann können sich Ladungsträger durch einen korrelierten magnetischen Hintergrund bewegen, wodurch zusätzliche konkurrierende Phasen entstehen.
Temperatur
Viele der interessantesten Hubbard-Phänomene sind Niedrigtemperaturphänomene. Die relevante thermische Energieskala ist \(k_B T\). Antiferromagnetische Korrelationen oder supraleitende Paarung können nur beobachtet werden, wenn die Temperatur ausreichend klein gegenüber den charakteristischen Energieskalen des Systems ist.
Im stark wechselwirkenden Bereich wird beispielsweise die magnetische Austauschenergie deutlich kleiner als \(U\). Dadurch ist es wesentlich schwieriger, magnetische Ordnung zu erreichen als lediglich eine Mott-artige Unterdrückung der Doppelbesetzung.
Physikalische Grenzfälle
Nichtwechselwirkendes System
Für
\(U = 0\)
reduziert sich der Hamiltonoperator auf ein Tight-Binding-Modell nichtwechselwirkender Fermionen. Die Zustände können im Impulsraum beschrieben werden. Für ein eindimensionales Gitter mit Gitterkonstante \(a\) und Hopping nur zwischen nächsten Nachbarn ergibt sich beispielsweise die Dispersion
\(\epsilon(k) = -2t\cos(ka)\)
In einem zweidimensionalen quadratischen Gitter lautet sie
\(\epsilon(\mathbf{k}) = -2t\left[\cos(k_xa)+\cos(k_ya)\right]\)
In diesem Grenzfall existiert keine Mott-Physik. Das Verhalten wird durch Bandstruktur, Fermi-Energie und Fermi-Oberfläche bestimmt.
Atomarer Grenzfall
Für
\(t = 0\)
sind die Gitterplätze voneinander entkoppelt. Die Fermionen können nicht tunneln. Die Energie eines Gitterplatzes hängt nur von seiner Besetzung ab. Bei stark repulsivem \(U\) ist Doppelbesetzung energetisch ungünstig.
Dieser Grenzfall zeigt besonders klar, warum eine starke lokale Wechselwirkung zur Lokalisierung führen kann. Das reale Hubbard-Modell liegt jedoch gewöhnlich zwischen dem vollständig delokalisierten und dem vollständig atomaren Grenzfall.
Starke Kopplung
Für große positive Werte von \(U/t\) können reale Hoppingprozesse bei Halbbesetzung stark unterdrückt sein. Virtuelle Prozesse bleiben jedoch möglich. Ein Fermion kann kurzzeitig auf einen bereits besetzten Nachbarplatz tunneln, einen energetisch teuren Zwischenzustand erzeugen und anschließend zurückkehren oder mit dem anderen Fermion effektiv den Platz tauschen.
Diese Prozesse sind die Grundlage des antiferromagnetischen Superaustauschs und zeigen, dass starke Lokalisierung keineswegs bedeutet, dass die Quantendynamik verschwindet.
Attraktive Wechselwirkung
Für
\(U<0\)
wird Doppelbesetzung energetisch begünstigt. Fermionen mit entgegengesetztem Spin können gebundene Paare bilden. Im schwach attraktiven Bereich besteht eine Verbindung zur BCS-Theorie der Paarbildung. Bei starker Anziehung verhalten sich die Paare zunehmend wie kompakte zusammengesetzte Bosonen.
Das attraktive Hubbard-Modell eignet sich deshalb zur Untersuchung des Übergangs zwischen BCS-artiger Paarung und stark gebundenen Paarzuständen.
Mott-Isolator und korrelationsgetriebene Lokalisierung
Warum Bandtheorie allein nicht genügt
Eine der wichtigsten Leistungen des Hubbard-Modells besteht darin, einen isolierenden Zustand zu erklären, der aus einer einfachen Einteilchen-Bandtheorie nicht zwingend folgt. Ein teilweise gefülltes Band sollte nach elementarer Bandtheorie metallisch sein. Eine ausreichend starke Elektron-Elektron-Wechselwirkung kann die Ladungsbewegung jedoch blockieren.
Dieser Zustand wird als Mott-Isolator bezeichnet.
Entstehung des Mott-Zustands
Betrachtet man ein repulsives System nahe Halbbesetzung, befindet sich im stark korrelierten Grenzfall bevorzugt ein Fermion auf jedem Gitterplatz. Für einen Ladungstransport müsste ein Fermion auf einen bereits besetzten Platz wechseln. Dadurch entstünde eine Doppelbesetzung, die ungefähr die Energie \(U\) kostet.
Ist diese Energiestrafe gegenüber der kinetischen Energie groß, wird der Ladungstransport unterdrückt. Das System kann damit elektrisch isolierend werden, obwohl eine nichtwechselwirkende Bandbeschreibung einen metallischen Zustand erwarten ließe.
Hubbard-Bänder
In spektralen Beschreibungen stark korrelierter Systeme kann sich das ursprüngliche Band in ein unteres und ein oberes Hubbard-Band aufspalten. Vereinfacht entspricht das untere Hubbard-Band Zuständen ohne zusätzliche energetisch teure Doppelbesetzung, während das obere Hubbard-Band mit Zuständen verbunden ist, deren Erzeugung die Wechselwirkungsenergie kostet.
Zwischen diesen Bereichen kann eine Ladungslücke entstehen. Die genaue Struktur hängt von Dimension, Gitter, Temperatur, Magnetismus und Näherungsverfahren ab.
Mott-Übergang ist nicht überall identisch
Der Begriff Mott-Übergang muss dimensions- und modellabhängig verwendet werden. Im eindimensionalen repulsiven Hubbard-Modell ist das System bei exakter Halbbesetzung für jedes positive \(U\) im Ladungssektor isolierend. Im zweidimensionalen quadratischen Gitter spielen zusätzlich Fermi-Flächen-Nesting und antiferromagnetische Korrelationen eine zentrale Rolle. In dynamischer Molekularfeldtheorie lässt sich wiederum ein paramagnetischer Mott-Übergang untersuchen, bei dem magnetische Ordnung bewusst ausgeschlossen wird.
Es gibt daher nicht den einen universellen Mott-Übergang, der in jeder Dimension auf identische Weise abläuft.
Magnetismus und Superaustausch
Lokale Momente
Bei starker Repulsion und Halbbesetzung ist Doppelbesetzung selten. Jeder Gitterplatz trägt dann näherungsweise einen Spin. Die Ladungsbewegung ist eingeschränkt, die Spinorientierung bleibt jedoch ein dynamischer Freiheitsgrad.
Damit entsteht aus einem Modell beweglicher Fermionen bei niedriger Energie effektiv ein Spinproblem.
Virtuelles Hopping
Nehmen zwei benachbarte Plätze jeweils ein Fermion auf, kann ein virtueller Hoppingprozess einen kurzlebigen doppelt besetzten Zwischenzustand erzeugen. Dieser Zwischenzustand kostet die Energie \(U\). In zweiter Ordnung der Störungstheorie ergibt sich daraus eine effektive Austauschwechselwirkung der Größenordnung
\(J = \frac{4t^2}{U}\)
für das repulsive Hubbard-Modell mit Hopping zwischen nächsten Nachbarn.
Heisenberg-Grenzfall
Bei Halbbesetzung und
\(U \gg t\)
geht die Niedrigenergiephysik in ein antiferromagnetisches Heisenberg-Modell über. Der effektive Hamiltonoperator besitzt die Form
\(H_{\mathrm{eff}} = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j\)
mit positiver Austauschkonstante \(J\). Benachbarte Spins bevorzugen deshalb eine antiparallele Orientierung.
Antiferromagnetismus
In zwei- und dreidimensionalen Gittern können sich bei ausreichend niedrigen Temperaturen ausgeprägte antiferromagnetische Korrelationen entwickeln. Im idealisierten zweidimensionalen isotropen Heisenberg-Modell verhindert der Mermin-Wagner-Satz bei endlicher Temperatur eine echte langreichweitige kontinuierliche Symmetriebrechung. Am Grundzustand kann dagegen antiferromagnetische Ordnung auftreten.
In dreidimensionalen Systemen ist langreichweitige antiferromagnetische Ordnung auch bei endlicher Temperatur möglich. Für Quantensimulatoren stellt das Erreichen der entsprechend niedrigen Temperaturen eine wesentliche experimentelle Herausforderung dar.
Das dotierte Fermi-Hubbard-Modell
Abweichung von der Halbbesetzung
Besonders schwierig wird das Problem, wenn die Teilchenzahl von der Halbbesetzung abweicht. Werden Fermionen entfernt, entstehen Löcher. Werden zusätzliche Fermionen eingebracht, müssen verstärkt doppelt besetzte Zustände berücksichtigt werden.
Die beweglichen Ladungsträger treffen dabei auf einen korrelierten Spin-Hintergrund. Ihre Bewegung verändert die magnetischen Korrelationen, während diese Korrelationen umgekehrt die Ladungsbewegung beeinflussen. Die Trennung zwischen Ladungs- und Spinphysik wird dadurch wesentlich komplizierter.
Das t-J-Modell als starke-Kopplungs-Näherung
Im Grenzfall großer repulsiver Wechselwirkung kann das dotierte Hubbard-Modell näherungsweise auf das sogenannte t-J-Modell reduziert werden. Seine schematische Form lautet
\(H_{tJ} = -t\sum_{\langle i,j\rangle,\sigma}\left(\tilde{c}_{i\sigma}^{\dagger}\tilde{c}_{j\sigma}+\mathrm{h.c.}\right) + J\sum_{\langle i,j\rangle}\left(\mathbf{S}_i\cdot\mathbf{S}_j-\frac{1}{4}n_in_j\right)\)
Die Operatoren sind dabei auf einen Hilbertraum ohne Doppelbesetzung projiziert. Das Modell zeigt unmittelbar die Konkurrenz zwischen der Bewegungsenergie dotierter Ladungsträger und dem antiferromagnetischen Austausch.
Streifen und Ladungsordnung
Numerische Untersuchungen des zweidimensionalen Hubbard- und t-J-Modells zeigen in verschiedenen Parameterbereichen räumlich modulierte Strukturen. Dazu gehören sogenannte Stripes, bei denen Ladungs- und Spinkorrelationen gekoppelte periodische Muster ausbilden.
Solche Strukturen sind für die Diskussion der Cuprat-Hochtemperatursupraleiter relevant, da auch dort konkurrierende antiferromagnetische, ladungsgeordnete und supraleitende Tendenzen beobachtet werden.
Pseudogap
Ein weiterer wichtiger Bereich ist die Pseudogap-Physik. Dabei wird die Zustandsdichte in bestimmten Energie- und Impulsbereichen bereits oberhalb eines möglichen supraleitenden Übergangs deutlich reduziert, ohne dass eine vollständige gewöhnliche Bandlücke entsteht.
Cluster-Erweiterungen der dynamischen Molekularfeldtheorie und verschiedene numerische Verfahren finden im dotierten zweidimensionalen Hubbard-Modell pseudogapartige Zustände. Die genaue mikroskopische Interpretation bleibt jedoch Teil aktiver Forschung.
Zusammenhang mit Hochtemperatursupraleitung
Das zweidimensionale repulsive Einband-Hubbard-Modell gilt als eines der wichtigsten Minimalmodelle zur Untersuchung der Cuprat-Supraleitung. Besonders interessant ist die Möglichkeit einer Paarung mit \(d_{x^2-y^2}\)-Symmetrie.
Viele numerische Rechnungen liefern Hinweise auf starke d-Wellen-Paarungskorrelationen in relevanten Parameterbereichen. Dennoch ist nicht abschließend geklärt, in welchem Umfang das einfache Einband-Hubbard-Modell sämtliche wesentlichen Eigenschaften realer Cuprate quantitativ reproduziert. Reale Materialien besitzen zusätzliche Orbitale, langreichweitige Wechselwirkungen, Elektron-Phonon-Kopplung und strukturelle Details.
Dimensionsabhängigkeit
Eine Dimension
Das eindimensionale Hubbard-Modell nimmt eine Sonderstellung ein, weil es mit dem Bethe-Ansatz exakt lösbar ist. Seine Physik unterscheidet sich deutlich von gewöhnlichen Fermi-Flüssigkeiten.
Ein bemerkenswertes Phänomen ist die Spin-Ladungs-Trennung. Niedrigenergetische Spin- und Ladungsanregungen können sich mit unterschiedlichen Geschwindigkeiten ausbreiten. Das anschauliche Bild eines Elektrons als stabiles Quasiteilchen mit gleichzeitig fest verbundenem Spin und Ladung verliert damit seine Gültigkeit.
Bei Halbbesetzung führt jede positive repulsive Wechselwirkung zu einer Ladungslücke, während der Spinsektor lückenlos bleiben kann.
Zwei Dimensionen
Das zweidimensionale quadratische Hubbard-Modell ist besonders relevant und zugleich besonders schwer. Es besitzt keine allgemein nutzbare exakte Lösung. An Halbbesetzung dominieren bei niedriger Temperatur antiferromagnetische Korrelationen. Bei Dotierung können Pseudogap, Streifen, ungewöhnliche metallische Zustände und Paarungstendenzen auftreten.
Gerade dieser Bereich ist mit klassischen numerischen Methoden schwierig, weil viele Verfahren entweder durch endliche Systemgrößen, starke Verschränkung oder das fermionische Vorzeichenproblem begrenzt werden.
Drei Dimensionen
In drei Dimensionen kann bei endlicher Temperatur langreichweitige antiferromagnetische Ordnung entstehen. Ultrakalte Atome in dreidimensionalen optischen Gittern ermöglichen eine besonders saubere Realisierung dieses Regimes.
Die höhere Koordinationszahl verändert außerdem die relative Bedeutung räumlicher Fluktuationen. In Grenzfällen sehr hoher Dimension wird die dynamische Molekularfeldtheorie besonders leistungsfähig, weil die Selbstenergie lokal wird.
Numerische und analytische Lösungsverfahren
Warum das Modell rechnerisch schwierig ist
Ein Gitterplatz des spinhalbigen Hubbard-Modells besitzt vier lokale Zustände: leer, Spin-up, Spin-down und doppelt besetzt. Für \(N\) Plätze wächst der vollständige lokale Fock-Raum daher wie
\(4^N\)
Schon für relativ kleine Gitter wird eine direkte Speicherung aller Zustände unmöglich. Symmetrien und feste Teilchenzahlen reduzieren den tatsächlich benötigten Raum, ändern aber nichts am grundlegenden exponentiellen Wachstum.
Exakte Diagonalisierung
Bei der exakten Diagonalisierung wird der Hamiltonoperator für ein endliches Cluster konstruiert und numerisch diagonalisiert. Das Verfahren liefert innerhalb des gewählten Clusters praktisch exakte Eigenwerte und Eigenzustände.
Sein größter Nachteil ist die Systemgröße. Schon einige Dutzend Gitterplätze liegen für eine vollständige Behandlung typischerweise außerhalb realistischer Ressourcen. Trotzdem ist exakte Diagonalisierung sehr wertvoll, um andere numerische Methoden zu überprüfen und kleine Systeme detailliert zu analysieren.
Quanten-Monte-Carlo
Quanten-Monte-Carlo-Verfahren berechnen thermodynamische und Korrelationsgrößen über stochastische Stichproben. Für bestimmte Parameterbereiche des Hubbard-Modells können sie äußerst präzise Ergebnisse liefern.
Ein zentrales Hindernis ist jedoch das fermionische Vorzeichenproblem. Beiträge zur Zustandssumme können positive und negative Gewichte besitzen. Bei ungünstigen Parametern heben sich große Beiträge nahezu gegenseitig auf, wodurch die statistische Unsicherheit exponentiell anwachsen kann.
Am halbbesetzten repulsiven Hubbard-Modell auf einem bipartiten Gitter kann das Vorzeichenproblem aufgrund von Symmetrien verschwinden. Bei Dotierung tritt es im Allgemeinen wieder auf. Genau dort befinden sich jedoch viele der physikalisch interessantesten offenen Fragen.
Dichtematrix-Renormierungsgruppe
Die Dichtematrix-Renormierungsgruppe, kurz DMRG, gehört zu den leistungsfähigsten Methoden für eindimensionale Quantensysteme. Sie nutzt die Tatsache, dass niedrigenergetische Zustände vieler eindimensionaler Systeme trotz des riesigen Hilbertraums durch Matrixproduktzustände effizient beschrieben werden können.
Für zweidimensionale Systeme wird DMRG häufig auf langen Zylindern eingesetzt. Dadurch lassen sich größere Systeme als mit exakter Diagonalisierung untersuchen. Mit wachsender Zylinderbreite steigt der Rechenaufwand jedoch schnell an, da auch die Verschränkung zunimmt.
Tensornetzwerke
Allgemeinere Tensornetzwerkmethoden versuchen, die Struktur der Vielteilchenverschränkung direkt auszunutzen. Matrixproduktzustände sind besonders für eine Dimension geeignet. Für zwei Dimensionen werden unter anderem PEPS-artige Zustände verwendet.
Diese Verfahren sind konzeptionell attraktiv, weil sie nicht auf stochastischen Gewichten beruhen und deshalb kein Vorzeichenproblem im üblichen Monte-Carlo-Sinn besitzen. Ihre numerische Optimierung kann jedoch sehr aufwendig sein.
Dynamische Molekularfeldtheorie
Die dynamische Molekularfeldtheorie, DMFT, bildet das Gitterproblem auf ein selbstkonsistentes Quantenimpurity-Problem ab. Räumliche Korrelationen werden dabei in der einfachsten Form vernachlässigt, lokale zeitabhängige Quantenfluktuationen dagegen vollständig behandelt.
Die zentrale Näherung lautet, dass die Selbstenergie nur von der Frequenz und nicht vom Impuls abhängt:
\(\Sigma(\mathbf{k},\omega) \approx \Sigma(\omega)\)
DMFT eignet sich besonders zur Untersuchung lokaler Korrelationsphänomene, quasiteilchenartiger Zustände und des Mott-Übergangs. In niedrigen Dimensionen sind jedoch räumliche Korrelationen oft wesentlich.
Cluster-DMFT
Cluster-Erweiterungen wie Cellular DMFT oder die Dynamical Cluster Approximation behandeln mehrere Gitterplätze gemeinsam. Dadurch können kurzreichweitige räumliche Korrelationen berücksichtigt werden.
Diese Verfahren spielen eine wichtige Rolle bei Untersuchungen des zweidimensionalen Hubbard-Modells, insbesondere bei Pseudogap-Physik und Paarung. Die erreichbaren Clustergrößen bleiben allerdings begrenzt.
Keine einzelne Methode genügt
Das heutige Verständnis des Hubbard-Modells basiert deshalb auf dem Vergleich verschiedener Methoden. Besonders überzeugend sind Resultate, wenn unabhängige Verfahren im selben Parameterbereich übereinstimmen. Abweichungen zwischen Methoden sind dagegen häufig ein Hinweis darauf, dass endliche Systemgröße, Temperatur, Näherungsstruktur oder numerische Konvergenz noch relevant sind.
Quantensimulation mit ultrakalten Atomen
Optische Gitter
Ultrakalte Atome bieten eine direkte Möglichkeit, den Fermi-Hubbard-Hamiltonoperator experimentell nachzubilden. Mehrere gegenläufige Laserstrahlen erzeugen stehende Lichtwellen. Durch den AC-Stark-Effekt entsteht für die Atome ein periodisches Potential.
In einer Raumrichtung kann ein solches Potential beispielsweise durch
\(V(x) = V_0\sin^2(k_Lx)\)
beschrieben werden. Mehrere orthogonale Laserachsen erzeugen zwei- oder dreidimensionale Gitter.
Bei ausreichend tiefem Gitter sind die Atome in den einzelnen Potentialminima lokalisiert, können aber zwischen benachbarten Minima tunneln. Damit entsteht unmittelbar die für das Hubbard-Modell benötigte Gitterstruktur.
Fermionische Atomsorten
Zu den häufig verwendeten fermionischen Atomen gehören Lithium-6 und Kalium-40. Zwei ausgewählte Hyperfeinzustände übernehmen die Rolle der beiden Spinrichtungen.
Der große Vorteil gegenüber Elektronen in einem Festkörper besteht in der Kontrollierbarkeit. Die Gittertiefe, die Wechselwirkung, die Teilchenzahl, die Geometrie und in vielen Fällen auch lokale Potentiale können gezielt verändert werden.
Kontrolle des Hoppingparameters
Die Tunnelamplitude \(t\) hängt stark von der Tiefe des optischen Gitters ab. Wird das Gitter tiefer, werden die Wannier-Funktionen stärker lokalisiert und ihre Überlappung nimmt ab. Dadurch sinkt die Tunnelamplitude.
Experimentell kann daher das Verhältnis \(U/t\) bereits durch Änderung der Laserintensität verändert werden.
Kontrolle der Wechselwirkung
Die Wechselwirkung zwischen ultrakalten Atomen lässt sich über magnetische Feshbach-Resonanzen einstellen. Dabei verändert ein äußeres Magnetfeld die effektive s-Wellen-Streulänge \(a_s\). Die Kontaktwechselwirkung besitzt näherungsweise die Kopplungsstärke
\(g = \frac{4\pi\hbar^2a_s}{m}\)
mit der Atommasse \(m\). Über \(a_s\) lässt sich damit auch der effektive Hubbard-Parameter \(U\) kontrollieren.
Das ist ein entscheidender Unterschied zu Festkörpern: Im Quantensimulator ist die Wechselwirkung kein unveränderliches Materialmerkmal, sondern ein experimenteller Parameter.
Erzeugung eines Mott-Regimes
Werden fermionische Atome mit starker repulsiver Wechselwirkung bei geeigneter Dichte in ein optisches Gitter geladen, kann ein Bereich mit ungefähr einem Atom pro Gitterplatz entstehen. Doppelbesetzungen werden unterdrückt. Damit wird ein Hubbard-Mott-Regime realisiert.
Die größte Herausforderung besteht anschließend darin, die Temperatur weit genug zu reduzieren, um magnetische Energieskalen zu erreichen. Da
\(J = \frac{4t^2}{U}\)
bei großer Wechselwirkung deutlich kleiner als \(U\) sein kann, erfordert die Beobachtung von Antiferromagnetismus erheblich niedrigere Temperaturen als die Beobachtung bloßer Mott-Lokalisierung.
Quantum-Gas-Mikroskope
Quantum-Gas-Mikroskope haben die experimentellen Möglichkeiten wesentlich erweitert. Mit ihnen können einzelne Atome auf einzelnen Gitterplätzen detektiert werden. Damit wird der Quantenzustand nicht nur über globale Mittelwerte untersucht, sondern mikroskopisch räumlich aufgelöst.
Messbar sind unter anderem lokale Dichten, Doppelbesetzungen, Spin-Korrelationen, Ladungskorrelationen und die Bewegung einzelner Defekte. Durch Wiederholung des Experiments entstehen statistische Momentaufnahmen, aus denen Vielteilchen-Korrelationsfunktionen rekonstruiert werden können.
Direkte Beobachtung magnetischer Korrelationen
Für zwei Gitterplätze können beispielsweise Spin-Korrelationen der Form
\(C_{ij} = \langle S_i^zS_j^z\rangle - \langle S_i^z\rangle\langle S_j^z\rangle\)
bestimmt werden. Antiferromagnetische Korrelationen zeigen sich daran, dass benachbarte Spins bevorzugt entgegengesetzt orientiert sind.
Moderne Experimente haben solche Korrelationen nicht nur an Halbbesetzung, sondern auch in dotierten Systemen sichtbar gemacht. Dadurch können theoretische Vorhersagen des Hubbard-Modells direkt auf der Ebene einzelner Gitterplätze überprüft werden.
Grenzen analoger Quantensimulation
Auch ein Quantensimulator löst das Problem nicht automatisch. Die Temperatur muss ausreichend niedrig sein, die Modellparameter müssen präzise kalibriert werden und Messungen besitzen endliche statistische Unsicherheit. Zusätzlich können Fallenpotentiale, Heizprozesse, Teilchenverluste und nichtideale Gitterparameter die ideale Hubbard-Physik verändern.
Der entscheidende Vorteil liegt dennoch darin, dass das Quantensystem seine Vielteilchendynamik selbst realisiert. Ein klassischer Computer muss den exponentiell großen Zustandsraum nicht vollständig simulieren.
Das Fermi-Hubbard-Modell auf digitalen Quantencomputern
Fermionen müssen auf Qubits abgebildet werden
Digitale Quantencomputer arbeiten gewöhnlich mit Qubits. Ein fermionischer Hamiltonoperator kann deshalb nicht unmittelbar auf der Hardware ausgeführt werden. Die fermionischen Operatoren müssen zunächst in Qubit-Operatoren transformiert werden.
Eine bekannte Methode ist die Jordan-Wigner-Transformation. Für einen fermionischen Modus \(j\) lautet sie schematisch
\(c_j = \left(\prod_{k und \(c_j^{\dagger} = \left(\prod_{k Die Kette aus Pauli-\(Z\)-Operatoren stellt die korrekten fermionischen Vorzeichen sicher. Alternativ kann die Bravyi-Kitaev-Transformation verwendet werden. Sie verteilt Informationen über Besetzung und Parität anders auf die Qubits und kann für bestimmte Algorithmen günstigere Operatorlängen erzeugen. Darüber hinaus existieren speziell auf Gitterfermionen zugeschnittene Kodierungen. Die optimale Wahl hängt von Hardwaretopologie, Fehlerkorrektur, Gattermenge und zu simulierender Observable ab. Ein häufig diskutierter Ansatz für gegenwärtige Quantenprozessoren ist der Variational Quantum Eigensolver. Ein parametrisierter Quantenschaltkreis erzeugt einen Zustand \(|\psi(\boldsymbol{\theta})\rangle\) und minimiert die Energie \(E(\boldsymbol{\theta}) = \langle\psi(\boldsymbol{\theta})|H|\psi(\boldsymbol{\theta})\rangle\) mithilfe einer klassischen Optimierung. Für kleine Hubbard-Systeme ist dies experimentell demonstrierbar. Für große, stark korrelierte Gitter sind jedoch Schaltkreistiefe, Optimierungslandschaft und Messaufwand erhebliche Hindernisse. Quantencomputer können grundsätzlich auch die zeitliche Entwicklung \(|\psi(t)\rangle = e^{-iHt/\hbar}|\psi(0)\rangle\) simulieren. Das ist besonders interessant für Nichtgleichgewichtsphänomene, Quenches, Transport und Spektroskopie. Fehlertolerante Quantencomputer könnten Hamilton-Simulationsalgorithmen mit kontrollierbarer Genauigkeit verwenden. Dafür werden allerdings logische Qubits und umfangreiche Fehlerkorrektur benötigt. Das Hubbard-Modell gilt als aussichtsreiches Ziel für wissenschaftlich relevanten Quantenvorteil, weil klassische Methoden gerade in bestimmten dotierten und dynamischen Regimen erhebliche Schwierigkeiten haben. Ein Quantenvorteil ist jedoch erst dann erreicht, wenn ein Quantenprozessor für eine physikalisch relevante Observable verlässlich bessere Resultate liefert als die besten verfügbaren klassischen Verfahren. Kleine Demonstrationen oder die bloße Implementierung des Hamiltonoperators genügen dafür nicht. In realen Materialien ist Hopping nicht zwingend auf nächste Nachbarn beschränkt. Ein verbreiteter Zusatz ist ein Hoppingparameter \(t'\) zwischen übernächsten Nachbarn. Der kinetische Anteil kann dann beispielsweise durch \(H_t = -t\sum_{\langle i,j\rangle,\sigma}c_{i\sigma}^{\dagger}c_{j\sigma} - t'\sum_{\langle\langle i,j\rangle\rangle,\sigma}c_{i\sigma}^{\dagger}c_{j\sigma} + \mathrm{h.c.}\) erweitert werden. Bereits ein solches zusätzliches Hopping kann Fermi-Oberfläche, magnetische Frustration und Phasendiagramm deutlich verändern. Werden zusätzlich Wechselwirkungen zwischen Fermionen auf unterschiedlichen Gitterplätzen berücksichtigt, entsteht ein erweitertes Hubbard-Modell. Ein typischer Zusatzterm lautet \(V\sum_{\langle i,j\rangle}n_in_j\) Solche Terme können Ladungsordnung und andere konkurrierende Phasen stabilisieren. Besitzen mehrere Orbitale in der relevanten Energieregion wesentlichen Einfluss, reicht das Einbandmodell nicht aus. Mehrbandmodelle enthalten zusätzliche Orbitalindizes sowie intraorbitale und interorbitale Wechselwirkungen, Hund-Kopplungen und gegebenenfalls Paartransferterme. Die Wahl zwischen Einband- und Mehrbandbeschreibung ist daher eine physikalische Modellentscheidung und keine reine Frage mathematischer Bequemlichkeit. Das Bose-Hubbard-Modell besitzt eine ähnliche Gitterstruktur, beschreibt jedoch Bosonen. Aufgrund der unterschiedlichen Quantenstatistik ist seine Physik deutlich anders. Mehrere Bosonen dürfen denselben Einteilchenzustand besetzen, und das Modell zeigt beispielsweise den Übergang zwischen Superfluid und Mott-Isolator. Das Hubbard-Modell wurde entwickelt, um Elektronenkorrelationen in schmalen Energiebändern zu verstehen. Ist die Bandbreite klein, kann die lokale Coulomb-Wechselwirkung mit der kinetischen Energie konkurrieren. Genau dann versagt das Bild nahezu unabhängiger Elektronen besonders deutlich. Relevante Materialklassen umfassen Übergangsmetalloxide, bestimmte organische Leiter, Mott-Materialien und verschiedene magnetische Systeme. Bei Cuprat-Hochtemperatursupraleitern bilden die elektronisch relevanten Kupfer-Sauerstoff-Ebenen annähernd zweidimensionale Strukturen. Aus komplexeren Mehrorbitalmodellen kann unter bestimmten Annahmen eine effektive Einbandbeschreibung abgeleitet werden. Das zweidimensionale Hubbard-Modell reproduziert mehrere qualitative Merkmale, die auch bei Cupraten auftreten: antiferromagnetische Physik nahe Halbbesetzung, starke Korrelationen, Pseudogap-Tendenzen, Ladungsmodulationen und d-Wellen-Paarung. Es wäre jedoch zu weitgehend, daraus eine vollständige mikroskopische Gleichsetzung abzuleiten. Materialabhängige Bandparameter, Sauerstofforbitale, langreichweitige Coulomb-Kräfte, Gitterkopplungen und strukturelle Effekte können relevant sein. Der eigentliche Wert des Einband-Fermi-Hubbard-Modells liegt deshalb nicht darin, jedes reale Material exakt nachzubilden. Seine Stärke besteht darin, klar zu isolieren, welche komplexen Vielteilchenphänomene bereits aus lokal wechselwirkenden Fermionen in einem einzelnen Band entstehen können. Das Fermi-Hubbard-Modell ist eines der wichtigsten Referenzmodelle analoger Quantensimulation. Seine Parameter sind einfach definiert, seine Phasen jedoch nicht trivial. Dadurch lässt sich prüfen, ob ein Quantensimulator tatsächlich bekannte Vielteilchenphysik reproduziert und anschließend in schwer zugängliche Regime vorgedrungen werden kann. Auch digitale Quantencomputer profitieren von diesem Charakter. Das Modell ist komplex genug, um klassische Verfahren herauszufordern, gleichzeitig aber strukturiert genug, um kontrollierte Hamiltonoperatoren, Observablen und bekannte Grenzfälle bereitzustellen. Damit eignet es sich zum Vergleich von Fermionenkodierungen, Variationsansätzen, Zustandspräparation, Hamilton-Simulation, Messstrategien und Fehlerkorrekturverfahren. Stark korrelierte Hubbard-Zustände enthalten räumlich verteilte Quantenkorrelationen, die nicht durch ein einfaches Produkt unabhängiger Einteilchenzustände beschrieben werden können. Ihre kontrollierte Erzeugung und Messung ist deshalb auch aus Sicht der Quanteninformation interessant. Quantum-Gas-Mikroskope ermöglichen beispielsweise die Rekonstruktion höherer Korrelationsfunktionen und die Untersuchung von Verschränkungsindikatoren. Auf digitalen Quantencomputern können verwandte Größen direkt aus Qubitmessungen bestimmt werden. Ein besonders vielversprechender Bereich ist die Dynamik fern vom thermischen Gleichgewicht. Nach einer plötzlichen Parameteränderung kann untersucht werden, wie sich Ladung, Spin, Doppelbesetzung und Korrelationen zeitlich entwickeln. Solche Probleme sind klassisch oftmals noch schwieriger als Grundzustandsberechnungen, weil die quantenmechanische Verschränkung mit der Zeit stark anwachsen kann. Analoge und digitale Quantensimulatoren könnten hier langfristig einen besonders großen wissenschaftlichen Nutzen besitzen. Eine der zentralen offenen Aufgaben ist das quantitative Phasendiagramm des zweidimensionalen repulsiven Hubbard-Modells bei niedriger Temperatur und endlicher Dotierung. Verschiedene numerische Verfahren finden starke Konkurrenz zwischen Antiferromagnetismus, Streifen, Paarung und metallischen beziehungsweise pseudogapartigen Zuständen. Die Unterschiede zwischen einzelnen Ergebnissen sind teilweise gering genug, dass endliche Systemgrößen, Randbedingungen und Näherungen die Zuordnung der stabilsten Phase beeinflussen können. Das Vorzeichenproblem verhindert in vielen wichtigen Parameterbereichen eine einfache Skalierung klassischer Monte-Carlo-Simulationen. Es ist nicht nur ein technisches Implementierungsproblem. Für allgemeine fermionische Quantensysteme kann es eine fundamentale rechnerische Schwierigkeit darstellen. Das erklärt, weshalb gerade das dotierte Hubbard-Modell häufig als möglicher Anwendungsfall für Quantensimulation diskutiert wird. Sowohl numerische als auch experimentelle Verfahren haben Schwierigkeiten mit sehr niedrigen Temperaturen. Klassische Rechnungen werden aufwendiger, während ultrakalte Atomexperimente Entropie effizient aus dem relevanten Gitterbereich entfernen müssen. Gerade supraleitende oder sehr langreichweitige magnetische Korrelationen können jedoch Energieskalen besitzen, die deutlich unter der primären Wechselwirkungsenergie liegen. Die Einfachheit des Modells ist gleichzeitig seine wichtigste Grenze. In realen Materialien können mehrere Orbitale, Spin-Bahn-Kopplung, Gitterschwingungen, nichtlokale Wechselwirkungen und strukturelle Verzerrungen eine wesentliche Rolle spielen. Das Einbandmodell kann daher nur dann quantitativ belastbar sein, wenn die vernachlässigten Freiheitsgrade tatsächlich auf die betrachtete Energieskala integriert werden können, ohne neue dominante Prozesse zu erzeugen. Auch Quantencomputer lösen das Hubbard-Problem nicht automatisch. Zustandspräparation, lange kohärente Zeitentwicklung, präzise Messung und Fehlerkorrektur verursachen erhebliche Ressourcenanforderungen. Der wissenschaftlich entscheidende Meilenstein wird nicht die Simulation eines kleinen Hubbard-Clusters sein, sondern die zuverlässige Berechnung oder Messung eines relevanten Regimes, das mit den besten klassischen Methoden nicht mehr vergleichbar effizient zugänglich ist. Das Einband-Fermi-Hubbard-Modell reduziert die Physik wechselwirkender Fermionen auf zwei fundamentale Prozesse: Bewegung zwischen Gitterplätzen und lokale Wechselwirkung. Der Hamiltonoperator \(H = -t\sum_{\langle i,j\rangle,\sigma}\left(c_{i\sigma}^{\dagger}c_{j\sigma}+c_{j\sigma}^{\dagger}c_{i\sigma}\right)+U\sum_i n_{i\uparrow}n_{i\downarrow}-\mu\sum_{i,\sigma}n_{i\sigma}\) ist kompakt. Seine Vielteilchenphysik ist es nicht. Bei schwacher Wechselwirkung steht die Bewegung der Fermionen im Vordergrund. Bei starker repulsiver Wechselwirkung werden Doppelbesetzungen unterdrückt, Mott-Physik entsteht und virtuelle Tunnelprozesse erzeugen antiferromagnetischen Superaustausch. Bei Dotierung geraten Ladungsbewegung, Magnetismus und Paarung in direkte Konkurrenz. Gerade das zweidimensionale dotierte Modell gehört deshalb zu den anspruchsvollsten Standardproblemen der theoretischen Physik. Es ist eng mit der Frage verbunden, welche minimale mikroskopische Physik für unkonventionelle Supraleitung und andere Phänomene stark korrelierter Elektronensysteme erforderlich ist. Für die Quantentechnologie besitzt das Modell eine zweite Rolle. Ultrakalte Fermionen in optischen Gittern ermöglichen eine kontrollierbare analoge Realisierung. Quantum-Gas-Mikroskope machen einzelne Teilchen und Korrelationen direkt sichtbar. Digitale Quantencomputer verwenden das Hubbard-Modell zugleich als Benchmark für Fermionenkodierung, Vielteilchenalgorithmen und Hamilton-Simulation. Das Einband-Fermi-Hubbard-Modell ist damit kein bloßes Lehrbuchmodell. Es bildet eine gemeinsame Plattform für Festkörperphysik, Quantenmaterie, Atomphysik, Quanteninformation und Quantencomputing. Seine Stärke liegt gerade darin, mit wenigen Parametern eine Physik zu erzeugen, die bis heute nicht in allen relevanten Regimen vollständig verstanden ist. Die wissenschaftliche Literatur zum Einband-Fermi-Hubbard-Modell reicht von den grundlegenden Arbeiten der frühen 1960er-Jahre bis zu modernen numerischen Untersuchungen und Experimenten mit ultrakalten Fermionen. Für eine wissenschaftliche Abhandlung sollten Primärquellen zur Entstehung des Modells mit aktuellen Übersichtsarbeiten, methodischen Vergleichsstudien und experimentellen Arbeiten kombiniert werden. Monographien sind für das Einband-Fermi-Hubbard-Modell vor allem dort wichtig, wo die Abhandlung mathematische Grundlagen, Vielteilchentheorie, Quantenmagnetismus, exakte Lösungen oder Quantensimulation systematisch entwickeln soll. Die folgenden Werke ergänzen die Primärliteratur und eignen sich besonders für Herleitungen und methodische Hintergrundabschnitte. Online-Ressourcen sollten in einer wissenschaftlichen Arbeit vor allem zur Recherche, zur Kontrolle bibliografischer Angaben und zum Auffinden neuer Veröffentlichungen genutzt werden. Für zentrale physikalische Aussagen sollten nach Möglichkeit die Originalartikel, begutachtete Übersichtsarbeiten oder wissenschaftliche Monographien zitiert werden.Bravyi-Kitaev und andere Kodierungen
Variationsalgorithmen
Echtzeitentwicklung
Quantenvorteil
Erweiterungen und verwandte Modelle
Weiterreichendes Hopping
Erweiterte Hubbard-Modelle
Mehrband-Hubbard-Modelle
Bose-Hubbard-Modell
Bedeutung für reale Materialien
Stark korrelierte Elektronensysteme
Cuprate
Minimalmodell statt Materialkopie
Bedeutung für die Quantentechnologie
Referenzproblem für Quantensimulatoren
Benchmark für Quantencomputer
Untersuchung echter Vielteilchenverschränkung
Nichtgleichgewichtsdynamik
Grenzen und offene Forschungsfragen
Das zweidimensionale dotierte Modell
Fermionisches Vorzeichenproblem
Niedrige Temperaturen
Grenzen des Einbandansatzes
Quantenvorteil bleibt nachzuweisen
Fazit
Anhang
Wissenschaftliche Zeitschriften und Artikel
Grundlegende Primärliteratur zum Einband-Fermi-Hubbard-Modell
Übersichts- und Methodenliteratur zu stark korrelierten Fermionen
Quantensimulation und ultrakalte Fermionen
Hintergrundliteratur zur Quantensimulation
Bücher und Monographien
Spezialisierte Arbeiten zum Hubbard-Modell und Quantenmagnetismus
Hintergrundliteratur zur Vielteilchenphysik
Standardwerke zur Quanteninformation
Online-Ressourcen und Datenbanken
Vorlesungsnotizen und Monographie-nahe Ressourcen
Fachjournale und Verlage
Lern- und Forschungsplattformen
