Das erweiterte Hubbard-Modell gehört zu den zentralen theoretischen Werkzeugen für die Untersuchung stark korrelierter Quantensysteme. Es beschreibt Teilchen auf einem Gitter, deren Verhalten nicht allein durch ihre Bewegungsenergie bestimmt wird, sondern wesentlich durch ihre gegenseitige Wechselwirkung. Während das ursprüngliche Hubbard-Modell vor allem die Konkurrenz zwischen Teilchenbewegung und lokaler Abstoßung erfasst, berücksichtigt die erweiterte Variante zusätzliche Wechselwirkungen zwischen verschiedenen Gitterplätzen.
Diese Erweiterung ist physikalisch wichtig. Elektronen in einem realen Material wechselwirken nicht nur dann miteinander, wenn sie sich am selben Ort befinden. Die Coulomb-Wechselwirkung besitzt grundsätzlich eine räumliche Reichweite. In vielen Festkörpern, künstlichen Quantensystemen und Quantensimulatoren können daher Wechselwirkungen zwischen benachbarten oder weiter entfernten Gitterplätzen nicht vernachlässigt werden.
Das erweiterte Hubbard-Modell ermöglicht dadurch die Beschreibung einer deutlich größeren Vielfalt kollektiver Quantenzustände. Dazu gehören Mott-Isolatoren, Ladungsdichtewellen, Spindichtewellen, magnetische Phasen, Bindungsordnungen und verschiedene Formen der Supraleitung. Besonders interessant ist die Konkurrenz dieser Zustände. Schon kleine Änderungen einzelner Modellparameter können das System von einer Phase in eine andere überführen.
Für die Quantentechnologie besitzt das Modell eine doppelte Bedeutung. Einerseits dient es als theoretisches Modell für Quantenmaterialien und stark korrelierte Materie. Andererseits ist es ein wichtiges Zielproblem für analoge Quantensimulatoren und digitale Quantencomputer. Gerade weil seine Lösung mit klassischen Computern schnell sehr aufwendig wird, eignet es sich dazu, die Leistungsfähigkeit zukünftiger Quantentechnologien zu untersuchen.
Ausgangspunkt: Das klassische Hubbard-Modell
Grundidee
Das klassische Hubbard-Modell beschreibt typischerweise Fermionen, etwa Elektronen, auf diskreten Gitterplätzen. Die Teilchen können zwischen benachbarten Plätzen wechseln. Gleichzeitig wirkt eine lokale Wechselwirkung, wenn sich zwei Teilchen mit entgegengesetztem Spin am selben Gitterplatz befinden.
Der Standard-Hamiltonoperator 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}\)
Der erste Term beschreibt die Bewegung der Teilchen. Der Parameter \(t\) ist die Hopping-Amplitude. Je größer sein Betrag ist, desto stärker neigen die Teilchen zur Delokalisierung über das Gitter.
Der zweite Term beschreibt die lokale Wechselwirkung. Der Parameter \(U\) gibt die Energieänderung an, wenn sich zwei Fermionen mit unterschiedlichem Spin am gleichen Gitterplatz befinden. Bei \(U>0\) ist diese Doppelbesetzung energetisch ungünstig. Bei \(U<0\) ist sie dagegen begünstigt.
Die Operatoren \(c_{i\sigma}^{\dagger}\) und \(c_{i\sigma}\) erzeugen beziehungsweise vernichten ein Fermion mit Spin \(\sigma\) auf dem Gitterplatz \(i\). Der Besetzungsoperator ist
\(n_{i\sigma}=c_{i\sigma}^{\dagger}c_{i\sigma}\)
und die gesamte Besetzung eines Gitterplatzes kann als
\(n_i=n_{i\uparrow}+n_{i\downarrow}\)
geschrieben werden.
Konkurrenz von Bewegung und Wechselwirkung
Die zentrale Physik des Hubbard-Modells entsteht durch die Konkurrenz zwischen \(t\) und \(U\). Eine große Hopping-Amplitude fördert bewegliche, delokalisierte Teilchen. Eine starke lokale Abstoßung verhindert dagegen Doppelbesetzungen und fördert lokalisierte Zustände.
Ein wichtiges Maß für die Stärke der Korrelationen ist deshalb das Verhältnis
\(\frac{U}{t}\)
Bei kleinen Werten kann das System häufig näherungsweise als schwach korreliertes Metall behandelt werden. Bei großen positiven Werten dominieren lokale Korrelationen. Insbesondere bei geeigneter Füllung kann dadurch ein Mott-Isolator entstehen: Ein Zustand, der nach einer einfachen Bandtheorie leitfähig erscheinen müsste, aufgrund starker Teilchenwechselwirkungen jedoch isolierend wird.
Grenzen des ursprünglichen Modells
Das klassische Hubbard-Modell reduziert die Elektron-Elektron-Wechselwirkung im Wesentlichen auf einen lokalen Term. Für viele grundlegende Fragestellungen ist diese Vereinfachung äußerst wertvoll. Sie kann jedoch zu grob werden, wenn Ladungskorrelationen über mehrere Gitterplätze eine bedeutende Rolle spielen.
Gerade in niederdimensionalen Materialien, organischen Leitern, molekularen Festkörpern, künstlichen Gittern und Systemen mit relativ schwacher Abschirmung können Wechselwirkungen zwischen benachbarten Teilchen erheblich sein. An diesem Punkt beginnt das eigentliche erweiterte Hubbard-Modell.
Definition des erweiterten Hubbard-Modells
Einführung einer nichtlokalen Wechselwirkung
Die einfachste und wichtigste Erweiterung besteht darin, zusätzlich zur lokalen Wechselwirkung \(U\) eine Dichte-Dichte-Wechselwirkung \(V\) zwischen benachbarten Gitterplätzen einzuführen.
Eine häufig verwendete Form lautet:
\(H = -t\sum_{\langle i,j\rangle,\sigma}\left(c_{i\sigma}^{\dagger}c_{j\sigma}+\mathrm{h.c.}\right)+U\sum_i n_{i\uparrow}n_{i\downarrow}+V\sum_{\langle i,j\rangle}n_i n_j\)
Der Zusatzterm
\(V\sum_{\langle i,j\rangle}n_i n_j\)
bestraft bei positivem \(V\) Konfigurationen, in denen benachbarte Gitterplätze stark besetzt sind. Dadurch entsteht ein neuer Wettbewerb. Die lokale Wechselwirkung \(U\) bestimmt, wie teuer eine Doppelbesetzung eines einzelnen Platzes ist. Die nichtlokale Wechselwirkung \(V\) bestimmt dagegen, wie günstig oder ungünstig bestimmte räumliche Ladungsverteilungen sind.
Warum die Erweiterung qualitativ neue Physik erzeugt
Ohne den Parameter \(V\) liegt der Schwerpunkt des Modells auf lokaler Korrelation und Magnetismus. Mit einer ausreichend starken nichtlokalen Abstoßung kann es dagegen energetisch günstiger werden, Ladung räumlich zu ordnen. Ein Teil der Gitterplätze wird dann bevorzugt stärker, der andere schwächer besetzt.
Damit kann eine Ladungsdichtewelle entstehen. Sie konkurriert unter anderem mit magnetischer Ordnung und mit der durch \(U\) geförderten gleichmäßigeren Einzelbesetzung der Gitterplätze.
Das erweiterte Modell besitzt deshalb kein einziges universelles Verhalten. Seine Physik hängt von Füllung, Dimensionalität, Gittergeometrie, Temperatur und den relativen Größen von \(t\), \(U\) und \(V\) ab.
Allgemeinere Form
Eine weiter verallgemeinerte Wechselwirkung kann als
\(H_V=\frac{1}{2}\sum_{i\neq j}V_{ij}n_i n_j\)
geschrieben werden. Dabei beschreibt \(V_{ij}\) die effektive Wechselwirkung zwischen den Gitterplätzen \(i\) und \(j\). Sie kann von deren Abstand und von den Eigenschaften des jeweiligen Systems abhängen.
Das erweiterte Hubbard-Modell ist deshalb weniger ein einzelner Hamiltonoperator als eine Familie verwandter Vielteilchenmodelle.
Physikalische Bedeutung der Modellparameter
Hopping-Amplitude
Der Parameter \(t\) setzt eine wesentliche Energieskala. Er bestimmt, wie leicht ein Teilchen seinen Gitterplatz verlassen und einen benachbarten Platz besetzen kann. Großes \(t\) begünstigt Delokalisierung und damit metallisches Verhalten.
Neben dem nächsten-Nachbar-Hopping können weitere Bewegungsprozesse auftreten. Ein nächstnächster Nachbar kann beispielsweise durch einen Parameter \(t'\) beschrieben werden. Solche Terme verändern die Bandstruktur und können magnetische oder supraleitende Tendenzen erheblich beeinflussen.
Lokale Wechselwirkung
\(U\) beschreibt die lokale Korrelationsenergie. Bei starker positiver Wechselwirkung versuchen die Teilchen, Doppelbesetzungen zu vermeiden. Im Bereich
\(U\gg t\)
ist die Bewegung stark durch die Wechselwirkung eingeschränkt. Dennoch können virtuelle Hopping-Prozesse stattfinden. Daraus kann ein effektiver magnetischer Austausch entstehen, dessen charakteristische Energieskala im einfachen Grenzfall ungefähr durch
\(J\approx\frac{4t^2}{U}\)
gegeben ist.
Nichtlokale Wechselwirkung
Der Parameter \(V\) verändert vor allem die räumliche Organisation der Ladung. Bei positiver Wechselwirkung werden stark besetzte Nachbarplätze energetisch ungünstig. Das kann regelmäßige Ladungsmuster erzeugen.
Die entscheidende Frage lautet daher nicht nur, ob \(U\) groß oder klein ist, sondern wie sich lokale und nichtlokale Wechselwirkungen zueinander verhalten. Das Verhältnis
\(\frac{V}{U}\)
kann bestimmen, ob eher magnetische oder ladungsgeordnete Phasen entstehen.
Chemisches Potential und Füllung
Die Teilchenzahl kann über ein chemisches Potential kontrolliert werden. Der entsprechende Term lautet beispielsweise
\(H_{\mu}=-\mu\sum_i n_i\)
Besonders wichtig ist die Halbfüllung. Bei einem Einbandmodell mit zwei möglichen Spinrichtungen entspricht sie im Mittel einem Fermion pro Gitterplatz:
\(\langle n_i\rangle=1\)
Abweichungen von dieser Situation werden als Dotierung beschrieben. Dotierung kann bestehende Ordnungszustände schwächen, neue metallische Phasen erzeugen oder Paarungstendenzen verändern.
Zusätzliche Terme im erweiterten Hubbard-Modell
Weiterreichende Coulomb-Wechselwirkungen
Die Beschränkung auf direkte Nachbarn ist ebenfalls eine Näherung. In realen Systemen kann die effektive Coulomb-Wechselwirkung mehrere Gitterabstände umfassen. Dann werden Terme der Form
\(\sum_{i,j}V_{ij}n_i n_j\)
relevant. Ihre Bedeutung hängt wesentlich von der elektronischen Abschirmung ab.
Weiterreichendes Hopping
Teilchen können auch zu nächstnächsten oder weiter entfernten Gitterplätzen tunneln. Ein entsprechender zusätzlicher Term kann beispielsweise geschrieben werden als
\(H_{t'}=-t'\sum_{\langle\langle i,j\rangle\rangle,\sigma}\left(c_{i\sigma}^{\dagger}c_{j\sigma}+\mathrm{h.c.}\right)\)
Solche Prozesse verändern nicht nur die Beweglichkeit, sondern auch die Geometrie der möglichen quantenmechanischen Pfade. Dadurch können Frustration und konkurrierende Ordnungen entstehen.
Austausch- und Pair-Hopping-Terme
Komplexere Modelle berücksichtigen einen direkten Spin-Austausch zwischen verschiedenen Gitterplätzen. Ein typischer Ausdruck ist
\(H_J=J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j\)
Daneben können Paar-Hopping-Prozesse auftreten, bei denen zwei korrelierte Teilchen gemeinsam zwischen Gitterplätzen übertragen werden.
Solche Erweiterungen sind vor allem dann relevant, wenn ein einfaches Dichte-Dichte-Modell die beobachtete Magnetik oder Paarbildung nicht ausreichend beschreibt.
Spin-Bahn-Kopplung
In Materialien mit starker Spin-Bahn-Kopplung sind Spin und räumliche Bewegung nicht unabhängig. Erweiterte Hubbard-Modelle können daher spinabhängige Hopping-Terme enthalten. Damit lassen sich unter anderem magnetische Anisotropien und bestimmte topologische Vielteilchenzustände untersuchen.
Mehrere Orbitale
Viele reale Materialien besitzen mehrere relevante Orbitale pro Gitterplatz. Dann reicht ein Einband-Hubbard-Modell nicht aus. Multi-Orbital-Modelle enthalten zusätzliche intra- und interorbitale Wechselwirkungen sowie Hund-Kopplungen.
Die Grundidee bleibt jedoch dieselbe: Bewegungsprozesse konkurrieren mit lokalen und nichtlokalen Wechselwirkungen.
Konkurrenz zwischen verschiedenen Quantenphasen
Metallische Zustände
Wenn die kinetische Energie gegenüber den Wechselwirkungen dominiert oder das System ausreichend von einer kommensurablen Füllung entfernt ist, können bewegliche Ladungsträger erhalten bleiben. Das System zeigt dann metallisches Verhalten.
Ein Metall in einem stark korrelierten System ist jedoch nicht automatisch mit einem einfachen freien Elektronengas gleichzusetzen. Quasiteilchen können stark renormiert sein, und ihre effektive Masse, Lebensdauer und spektrale Struktur werden durch Vielteilcheneffekte bestimmt.
Mott-Isolator
Der Mott-Isolator ist eines der bekanntesten Ergebnisse der Hubbard-Physik. Bei geeigneter Füllung verhindert eine starke lokale Abstoßung den freien Ladungstransport.
Die Blockierung entsteht nicht aus einer vollständig gefüllten Bandstruktur, sondern aus elektronischer Korrelation. Gerade deshalb besitzt die Mott-Physik große Bedeutung für die Theorie stark korrelierter Quantenmaterialien.
Spindichtewellen
Bei starker lokaler Abstoßung kann die Ladung relativ gleichmäßig verteilt bleiben, während sich die Spins räumlich ordnen. Ein typischer Fall ist antiferromagnetische Ordnung, bei der benachbarte Spins bevorzugt entgegengesetzt ausgerichtet sind.
Eine entsprechende Spin-Korrelationsfunktion kann etwa durch
\(C_S(i,j)=\langle S_i^z S_j^z\rangle\)
charakterisiert werden.
Ladungsdichtewellen
Eine starke nichtlokale Wechselwirkung kann die Bildung einer periodischen Ladungsstruktur fördern. In einem einfachen eindimensionalen Bild können stärker und schwächer besetzte Gitterplätze alternieren.
Ein möglicher Ladungsordnungsparameter besitzt beispielsweise die Form
\(m_{\mathrm{CDW}}=\frac{1}{L}\sum_j(-1)^j\langle n_j\rangle\)
mit der Systemgröße \(L\). Ein endlicher Wert im thermodynamischen Grenzfall zeigt eine entsprechende Ladungsordnung an.
Bindungsordnung
Besonders interessant ist, dass zwischen magnetisch dominierter und ladungsgeordneter Phase weitere Zustände auftreten können. In eindimensionalen erweiterten Hubbard-Modellen ist die Bond-Order-Wave-Phase ein wichtiges Beispiel.
Dabei liegt die wesentliche Ordnung nicht primär in der lokalen Ladung oder im lokalen Spin, sondern in alternierenden Bindungsstärken zwischen benachbarten Plätzen. Die Erwartungswerte der kinetischen Bindungsoperatoren bilden ein periodisches Muster.
Supraleitende Phasen
Je nach Parameterbereich und Modellvariante können auch supraleitende Korrelationen verstärkt werden. Besonders außerhalb der Halbfüllung kann die Konkurrenz zwischen Spin-, Ladungs- und Paarungskorrelationen komplex werden.
Nichtlokale Wechselwirkungen können Paarbildung sowohl unterdrücken als auch indirekt bestimmte unkonventionelle Paarungssymmetrien begünstigen. Deshalb ist nicht allein das Vorzeichen eines Parameters entscheidend, sondern die vollständige Vielteilchendynamik.
Phasendiagramm des erweiterten Hubbard-Modells
Der Parameterraum
Phasendiagramme ordnen verschiedene Grundzustände bestimmten Bereichen des Parameterraums zu. Für das einfache erweiterte Hubbard-Modell ist insbesondere die Ebene aus \(U\) und \(V\) von Interesse.
Bei positiver lokaler und nichtlokaler Abstoßung konkurrieren typischerweise Spinordnung und Ladungsordnung. Große Werte von \(U\) bestrafen Doppelbesetzungen. Große Werte von \(V\) bestrafen dagegen gleichzeitig besetzte Nachbarplätze.
Starkkopplungsbild
In einem vereinfachten halbgefüllten eindimensionalen Starkkopplungsbild kann man die Energien idealisierter Zustände vergleichen. Ein Zustand mit ungefähr einem Teilchen pro Platz profitiert von der Vermeidung lokaler Doppelbesetzungen. Eine stark ladungsgeordnete Konfiguration kann dagegen abwechselnd stärker und schwächer besetzte Plätze besitzen.
Die häufig diskutierte Konkurrenz liegt dadurch in der Nähe einer Relation der Größenordnung
\(U\approx 2V\)
Diese Beziehung ist jedoch keine universelle exakte Phasengrenze. Endliches Hopping, Quantenfluktuationen, Dimension und zusätzliche Modellterme verschieben und verändern das tatsächliche Phasendiagramm.
Zwischenphasen
Gerade in eindimensionalen Systemen kann zwischen einer Spin-Dichtewellen- und einer Ladungsdichtewellenphase eine Bindungsordnungsphase auftreten. Dies zeigt, dass ein einfacher Vergleich klassischer Energien nicht ausreicht.
Quantenfluktuationen erzeugen kollektive Zustände, die im rein statischen Grenzfall nicht offensichtlich sind. Das ist einer der Gründe, warum numerisch exakte oder hochpräzise Vielteilchenmethoden für das erweiterte Hubbard-Modell notwendig sind.
Eindimensionale Systeme
Quantenfluktuationen und Dimensionalität
Eindimensionale Quantensysteme besitzen eine besonders ausgeprägte Vielteilchenphysik. Quantenfluktuationen sind stärker als in höheren Dimensionen, und viele Vorstellungen aus der gewöhnlichen Fermi-Flüssigkeitstheorie verlieren ihre Gültigkeit.
Das eindimensionale erweiterte Hubbard-Modell ist deshalb kein bloßes vereinfachtes Lehrmodell. Es stellt einen eigenen physikalischen Grenzfall dar, der zugleich mit modernen Quantensimulatoren experimentell untersucht werden kann.
Luttinger-Flüssigkeit
Metallische eindimensionale Fermionensysteme werden häufig durch die Luttinger-Flüssigkeit beschrieben. Ihre elementaren Anregungen entsprechen nicht einfach gewöhnlichen Elektronenquasiteilchen.
Eine charakteristische Eigenschaft ist die mögliche Trennung von Spin- und Ladungsfreiheitsgraden. Störungen können sich dadurch mit unterschiedlichen Geschwindigkeiten durch das System ausbreiten.
Halbfüllung
Bei Halbfüllung treten besonders starke Wechselwirkungseffekte auf. Lokale und nichtlokale Abstoßung konkurrieren direkt miteinander. Abhängig vom Parameterbereich entstehen Spin-Dichtewellen, Ladungsdichtewellen und Bindungsordnung.
Die genaue Bestimmung der Phasengrenzen erfordert empfindliche numerische Größen wie Energieableitungen, Korrelationslängen, Ordnungsparameter und Anregungslücken.
Zwei- und höherdimensionale Systeme
Zweidimensionale Gitter
Zweidimensionale Hubbard-Systeme sind für die Theorie korrelierter Materialien besonders relevant. Ein prominentes Beispiel ist das quadratische Gitter. Hier konkurrieren magnetische Ordnung, Ladungsordnung, metallische Zustände und mögliche supraleitende Phasen.
Bereits die Lösung des gewöhnlichen zweidimensionalen Hubbard-Modells ist numerisch äußerst anspruchsvoll. Zusätzliche nichtlokale Wechselwirkungen vergrößern den relevanten Parameterraum weiter.
Gittergeometrie und Frustration
Auf einem Dreiecksgitter können nicht alle antiferromagnetischen Bindungen gleichzeitig optimal erfüllt werden. Dadurch entsteht geometrische Frustration.
Auf Waben-, Kagome-, Dreiecks- und anderen Gittern können deshalb qualitativ andere Grundzustände auftreten als auf einem eindimensionalen oder quadratischen Gitter. Das erweiterte Hubbard-Modell verbindet elektronische Korrelationen damit unmittelbar mit der Geometrie des zugrunde liegenden Quantensystems.
Dreidimensionale Systeme
In drei Dimensionen werden langreichweitige Ordnungszustände häufig stabiler gegenüber Quantenfluktuationen. Gleichzeitig steigt der numerische Aufwand erheblich.
Für reale Festkörper ist eine dreidimensionale Behandlung oftmals notwendig. Dennoch bleiben ein- und zweidimensionale Modelle wichtig, weil zahlreiche Materialien stark anisotrope elektronische Strukturen besitzen.
Stark korrelierte Elektronen und Quantenmaterialien
Jenseits unabhängiger Elektronen
Viele Standardmodelle der Festkörperphysik basieren zunächst auf weitgehend unabhängigen Elektronen, die sich in einem effektiven Potential bewegen. Dieser Ansatz erklärt Bandstrukturen und zahlreiche Eigenschaften gewöhnlicher Metalle und Halbleiter sehr erfolgreich.
Bei stark korrelierten Materialien reicht diese Näherung jedoch nicht aus. Der Zustand eines Elektrons hängt wesentlich davon ab, wo sich andere Elektronen befinden und welchen Spin sie besitzen.
Mott-Physik
Die Mott-Physik demonstriert besonders deutlich, dass ein Vielteilchensystem nicht aus der Summe unabhängiger Einteilchenzustände verstanden werden kann. Elektronische Abstoßung kann den Transport blockieren und gleichzeitig lokale magnetische Momente erzeugen.
Das erweiterte Hubbard-Modell ergänzt dieses Bild um räumliche Ladungskorrelationen. Dadurch können Mott-artige und ladungsgeordnete Zustände unmittelbar miteinander konkurrieren.
Relevante Materialklassen
Hubbard-artige Modelle werden unter anderem zur theoretischen Beschreibung von Übergangsmetallverbindungen, organischen Leitern, molekularen Kristallen, niedrigdimensionalen Materialien und anderen stark korrelierten Systemen eingesetzt.
Die Parameter eines effektiven Hubbard-Modells sind dabei keine universellen Naturkonstanten. Sie müssen für das jeweilige Material aus mikroskopischen Rechnungen, experimentellen Daten oder einer Kombination beider Ansätze bestimmt werden.
Modell und reales Material
Ein erweitertes Hubbard-Modell ist kein vollständiger atomistischer Ersatz für ein Material. Es reduziert das reale System bewusst auf die Freiheitsgrade, die für eine bestimmte Fragestellung als wesentlich angesehen werden.
Gerade diese Reduktion macht das Modell wissenschaftlich wertvoll: Sie erlaubt es, grundlegende Mechanismen voneinander zu trennen und gezielt zu untersuchen.
Supraleitung im erweiterten Hubbard-Modell
Paarbildung
Supraleitung erfordert korrelierte Ladungsträgerpaare und eine makroskopisch kohärente Phase. Im konventionellen BCS-Bild wird die effektive Anziehung zwischen Elektronen hauptsächlich durch Gitterschwingungen vermittelt.
In stark korrelierten Systemen können dagegen elektronische Wechselwirkungen selbst wesentlich zur Paarungsphysik beitragen. Das Hubbard-Modell und seine Erweiterungen dienen deshalb als Minimalmodelle für die Untersuchung unkonventioneller Supraleitung.
Paaroperatoren
Ein lokaler Singulett-Paaroperator kann beispielsweise geschrieben werden als
\(\Delta_i=c_{i\downarrow}c_{i\uparrow}\)
Eine zugehörige Paar-Korrelationsfunktion ist
\(P(i,j)=\langle\Delta_i^{\dagger}\Delta_j\rangle\)
Für nichtlokale Paarung werden Operatoren auf verschiedenen Gitterplätzen kombiniert. Ihre räumliche Struktur unterscheidet unter anderem verschiedene Paarungssymmetrien.
Einfluss von nichtlokaler Abstoßung
Eine starke positive Wechselwirkung \(V\) kann Paarbildung auf benachbarten Gitterplätzen direkt erschweren. Gleichzeitig verändert sie Spin- und Ladungsfluktuationen. Diese Fluktuationen können wiederum die effektive Paarungswechselwirkung beeinflussen.
Deshalb kann nicht aus einem einzelnen Parameter unmittelbar auf das Vorhandensein oder Fehlen von Supraleitung geschlossen werden. Entscheidend sind Gitter, Füllung, Wechselwirkungsstärke und Paarungskanal.
Konkurrenz mit Ladungsordnung
Supraleitung und Ladungsordnung konkurrieren häufig um dieselben elektronischen Freiheitsgrade. In manchen Parameterregimen kann die Nähe zu einer Ladungsordnungsinstabilität Paarungskorrelationen verändern oder verstärken. In anderen Bereichen wird Supraleitung durch eine stabile Ladungsdichtewelle unterdrückt.
Diese Konkurrenz ist einer der Gründe, warum das erweiterte Hubbard-Modell für die Theorie komplexer Quantenmaterialien interessant bleibt.
Mathematische und numerische Lösungsverfahren
Das Vielteilchenproblem
Die größte Schwierigkeit des Hubbard-Modells liegt nicht in der Länge seines Hamiltonoperators, sondern in der Größe seines Hilbertraums. Für ein Gitter mit \(L\) Plätzen besitzt jeder Platz im fermionischen Einbandmodell vier lokale Zustände: leer, Spin oben, Spin unten oder doppelt besetzt.
Ohne Einschränkung wächst die Zahl der Basiszustände daher wie
\(4^L\)
Diese exponentielle Skalierung macht eine direkte Behandlung größerer Systeme unmöglich.
Exakte Diagonalisierung
Bei kleinen Gittern kann der Hamiltonoperator als Matrix dargestellt und numerisch diagonalisiert werden. Exakte Diagonalisierung liefert sehr präzise Ergebnisse für endliche Systeme und erlaubt den direkten Zugang zu Grundzuständen und angeregten Zuständen.
Der Nachteil ist die exponentielle Skalierung des Speicher- und Rechenaufwands. Symmetrien und Teilchenzahlerhaltung reduzieren den relevanten Hilbertraum, lösen das grundlegende Skalierungsproblem aber nicht.
Mean-Field- und Hartree-Fock-Verfahren
Mean-Field-Methoden ersetzen komplizierte Wechselwirkungsterme teilweise durch effektive Einteilchenfelder. Beispielsweise können Produkte von Operatoren näherungsweise durch Erwartungswerte entkoppelt werden.
Diese Methoden sind rechnerisch vergleichsweise günstig und können mögliche Ordnungszustände identifizieren. Starke Quantenfluktuationen und hochgradig verschränkte Zustände werden jedoch nur unzureichend erfasst.
Dichtematrix-Renormierungsgruppe
Die Dichtematrix-Renormierungsgruppe, kurz DMRG, ist für viele eindimensionale stark korrelierte Systeme außerordentlich leistungsfähig. Moderne DMRG-Verfahren lassen sich als Optimierung von Matrix-Product-States verstehen.
Die Effizienz beruht darauf, dass niedrigenergetische Zustände eindimensionaler lokaler Hamiltonoperatoren häufig eine vergleichsweise begrenzte Verschränkungsstruktur besitzen.
Für breite zweidimensionale Systeme steigt der benötigte Rechenaufwand jedoch stark an.
Quanten-Monte-Carlo
Quanten-Monte-Carlo-Verfahren berechnen Vielteilchengrößen statistisch. Für geeignete Parameterregime können sie große Systeme mit hoher Genauigkeit untersuchen.
Bei Fermionen tritt jedoch häufig das Vorzeichenproblem auf. Beiträge zur statistischen Summe können unterschiedliche Vorzeichen besitzen und sich stark gegenseitig aufheben. Dadurch wächst der statistische Aufwand in problematischen Fällen exponentiell.
Das Vorzeichenproblem ist einer der wesentlichen Gründe, weshalb bestimmte fermionische Vielteilchenprobleme klassisch besonders schwierig bleiben.
Dynamical Mean-Field Theory
Die Dynamical Mean-Field Theory, DMFT, bildet ein Gitterproblem auf ein selbstkonsistentes Quanten-Impurity-Problem ab. Lokale zeitabhängige Korrelationen werden dabei nichtperturbativ behandelt.
Standard-DMFT besitzt jedoch eine lokale Selbstenergie. Für das erweiterte Hubbard-Modell sind nichtlokale Wechselwirkungen besonders wichtig. Deshalb wurden Verfahren wie Extended Dynamical Mean-Field Theory und Cluster-Erweiterungen entwickelt.
Funktionale Renormierungsgruppe
Die funktionale Renormierungsgruppe verfolgt, wie effektive Wechselwirkungen beim schrittweisen Eliminieren energetischer Freiheitsgrade verändert werden. Dadurch können konkurrierende Instabilitäten etwa in magnetischen, ladungsgeordneten oder supraleitenden Kanälen analysiert werden.
Besonders in schwacher bis mittlerer Kopplung liefert sie wertvolle Informationen über die führenden Ordnungsmechanismen.
Tensor-Netzwerke
Tensor-Netzwerk-Methoden komprimieren den Vielteilchenzustand, indem sie seine Verschränkungsstruktur gezielt ausnutzen. Matrix-Product-States sind die wichtigste eindimensionale Variante. Für höhere Dimensionen existieren komplexere Tensor-Netzwerk-Ansätze.
Ihr Vorteil liegt darin, quantenmechanische Korrelationen explizit in der Struktur der Wellenfunktion abzubilden. Der numerische Aufwand kann in stark verschränkten Systemen allerdings rasch steigen.
Quantensimulation des erweiterten Hubbard-Modells
Warum Quantensimulatoren?
Richard Feynmans grundlegende Motivation für Quantensimulation war einfach: Ein kontrolliertes Quantensystem sollte ein anderes Quantensystem effizienter darstellen können als ein klassischer Computer.
Hubbard-Modelle eignen sich hervorragend für diesen Ansatz. Sie besitzen überschaubare Hamiltonoperatoren, erzeugen aber hochkomplexe Vielteilchenzustände.
Ultrakalte Atome in optischen Gittern
Optische Gitter entstehen durch interferierende Laserfelder. Die resultierende periodische Potentiallandschaft kann für ultrakalte Atome eine künstliche Kristallstruktur erzeugen.
Die Gittertiefe beeinflusst die Beweglichkeit der Atome und damit die effektive Hopping-Amplitude. Wechselwirkungen zwischen den Atomen können beispielsweise über Feshbach-Resonanzen kontrolliert werden.
Dadurch entsteht eine bemerkenswerte Situation: Parameter, die in einem realen Kristall weitgehend durch die Materialchemie festgelegt sind, werden im Quantensimulator experimentell einstellbar.
Nichtlokale Wechselwirkungen
Für ein erweitertes Hubbard-Modell werden zusätzlich kontrollierbare nichtlokale Wechselwirkungen benötigt. Eine Möglichkeit bieten dipolare Atome oder polare Moleküle. Ihre Wechselwirkung besitzt eine natürliche räumliche Reichweite.
Auch Rydberg-Anregungen können starke langreichweitige Wechselwirkungen erzeugen. Weitere Ansätze verwenden durch optische Resonatoren vermittelte Kopplungen.
Damit lassen sich Hamiltonoperatoren realisieren, die über das einfache lokale Hubbard-Modell hinausgehen.
Quantengasmikroskopie
Quantengasmikroskope ermöglichen die Beobachtung einzelner Gitterplätze. Dadurch lassen sich lokale Besetzungen und daraus abgeleitete Korrelationsfunktionen direkt untersuchen.
Statt nur makroskopische Mittelwerte zu messen, kann experimentell verfolgt werden, wie sich Teilchen, Löcher und Spins über ein künstliches Quantengitter verteilen.
Diese Einzelplatzauflösung hat die Hubbard-Physik erheblich verändert. Viele theoretische Korrelationsgrößen sind damit zu experimentell zugänglichen Observablen geworden.
Umsetzung auf digitalen Quantencomputern
Fermionen auf Qubits
Ein digitaler Quantencomputer arbeitet typischerweise mit Qubits. Fermionische Erzeugungs- und Vernichtungsoperatoren müssen deshalb auf Qubit-Operatoren abgebildet werden.
Eine bekannte Methode ist die Jordan-Wigner-Transformation. Vereinfacht kann ein fermionischer Operator durch eine Kombination aus Pauli-Operatoren und einer Zeichenkette dargestellt werden.
Für einen Modus \(j\) besitzt die Transformation beispielsweise die Struktur
\(c_j=\frac{1}{2}\left(X_j+iY_j\right)\prod_{k Die Kette von \(Z\)-Operatoren stellt die korrekte fermionische Antikommutation sicher. Die Bravyi-Kitaev-Transformation verteilt Informationen über Besetzung und Parität anders auf die Qubits. Je nach Problem und Hardware kann sie günstigere Operatorlängen beziehungsweise Schaltungsstrukturen ermöglichen. Welche Abbildung optimal ist, hängt von Systemgröße, Konnektivität, Algorithmus und Fehlerkorrekturarchitektur ab. Die quantenmechanische Zeitentwicklung wird durch \(|\psi(t)\rangle=e^{-iHt}|\psi(0)\rangle\) bestimmt. Besteht der Hamiltonoperator aus nicht miteinander kommutierenden Teilen, kann die Evolution näherungsweise zerlegt werden. Für \(H=A+B\) lautet eine einfache Trotter-Näherung \(e^{-i(A+B)t}\approx\left(e^{-iA\Delta t}e^{-iB\Delta t}\right)^N\) mit \(\Delta t=\frac{t}{N}\) Je kleiner der Zeitschritt, desto geringer der Trotterfehler, gleichzeitig steigt jedoch die Zahl der benötigten Quantengatter. Der Variational Quantum Eigensolver versucht, die Grundzustandsenergie mithilfe einer parametrierten Quantenschaltung zu bestimmen. Ein Quantenzustand \(|\psi(\boldsymbol{\theta})\rangle\) wird auf dem Quantenprozessor erzeugt. Anschließend wird die Energie \(E(\boldsymbol{\theta})=\langle\psi(\boldsymbol{\theta})|H|\psi(\boldsymbol{\theta})\rangle\) gemessen. Ein klassischer Optimierer verändert die Parameter, um die Energie zu minimieren. Das Verfahren ist konzeptionell attraktiv für Hubbard-Modelle. Seine praktische Leistungsfähigkeit hängt jedoch stark von Ansatz, Messaufwand, Optimierungslandschaft und Hardwarefehlern ab. Auf einem fehlertoleranten Quantencomputer könnte Quantum Phase Estimation wesentlich präzisere Energieinformationen liefern. Voraussetzung ist eine ausreichend gute Approximation eines relevanten Eigenzustands sowie eine präzise kontrollierte Zeitentwicklung. Der Ressourcenbedarf ist beträchtlich. Deshalb gilt die Methode eher als Perspektive für zukünftige fehlertolerante Systeme als für heutige verrauschte Quantenprozessoren. Ein guter Benchmark für Quantenhardware sollte physikalisch relevant sein, schwierig genug bleiben und zugleich über einen klar definierten Hamiltonoperator verfügen. Das Hubbard-Modell erfüllt diese Bedingungen. Seine erweiterte Variante erhöht den Schwierigkeitsgrad, weil zusätzliche nichtlokale Wechselwirkungen und konkurrierende Phasen berücksichtigt werden müssen. Eine Behauptung über einen quantenmechanischen Rechenvorteil ist nur sinnvoll, wenn das Ergebnis mit den jeweils besten klassischen Verfahren verglichen wird. Für kleine Systeme sind exakte Methoden möglich. In eindimensionalen Fällen ist DMRG häufig extrem leistungsfähig. Andere Parameterregime können mit Quanten-Monte-Carlo- oder Tensor-Netzwerk-Methoden untersucht werden. Ein Quantenvorteil muss deshalb problem- und parameterabhängig bewertet werden. Die bloße Ausführung eines Hubbard-Hamiltonoperators auf Qubits ist noch kein Beleg für einen praktischen Vorteil. Hubbard-artige Probleme können außerdem genutzt werden, um Quantengatter, Qubit-Konnektivität, Fehlerunterdrückung und Zustandpräparation zu testen. Bekannte Ergebnisse kleiner Systeme dienen dabei als Referenz. Ultrakalte fermionische Atome gehören zu den direktesten Plattformen für Hubbard-Quantensimulation. Sie ermöglichen eine saubere Kontrolle von Gitterpotential, Teilchenzahl und Wechselwirkungsstärke. Da die Atome elektrisch neutral sind, unterscheidet sich ihre mikroskopische Physik von Elektronen im Festkörper. Entscheidend ist jedoch, dass der effektive Hamiltonoperator dieselbe Hubbard-Struktur annehmen kann. Für erweiterte Hubbard-Modelle sind Systeme mit dipolaren Wechselwirkungen besonders interessant. Die Wechselwirkung nimmt mit dem Abstand ab, endet aber nicht nach dem nächsten Gitterplatz. Dadurch können nichtlokale Kopplungen gezielt realisiert und in ihrer Stärke verändert werden. Rydberg-Atome besitzen stark angeregte elektronische Zustände und können sehr starke Wechselwirkungen über vergleichsweise große Abstände zeigen. Rydberg-Plattformen werden deshalb für programmierbare Quantensimulation eingesetzt. Nicht jede Rydberg-Implementierung bildet direkt ein fermionisches erweitertes Hubbard-Modell ab. Sie kann jedoch verwandte Gitter-Hamiltonoperatoren und wechselwirkende Vielteilchenphysik realisieren. Halbleiter-Quantenpunkte erlauben die kontrollierte Lokalisierung einzelner Elektronen. Kopplungen zwischen benachbarten Quantenpunkten können Hubbard-artige Hamiltonoperatoren erzeugen. Solche Systeme besitzen den besonderen Reiz, dass Ladung, Spin und Tunnelkopplungen innerhalb einer Festkörperplattform technisch kontrolliert werden können. Supraleitende Qubits simulieren Hubbard-Systeme nicht notwendigerweise durch identische mikroskopische Teilchen. Stattdessen wird der gewünschte Hamiltonoperator durch programmierbare Qubit-Wechselwirkungen nachgebildet. Damit gehören sie vor allem zur Kategorie digitaler oder digital-analoger Quantensimulatoren. Die lokale Teilchendichte ist eine der grundlegendsten Observablen: \(\langle n_i\rangle\) Sie zeigt unmittelbar, ob die Ladung gleichmäßig oder räumlich moduliert verteilt ist. Die Wahrscheinlichkeit einer Doppelbesetzung wird durch \(D_i=\langle n_{i\uparrow}n_{i\downarrow}\rangle\) beschrieben. Bei starker positiver lokaler Abstoßung wird sie typischerweise reduziert. Die doppelte Besetzung ist deshalb eine wichtige diagnostische Größe für lokale Korrelationen und Mott-Physik. Eine verbundene Dichtekorrelationsfunktion kann geschrieben werden als \(C_n(i,j)=\langle n_i n_j\rangle-\langle n_i\rangle\langle n_j\rangle\) Sie zeigt, wie stark Dichtefluktuationen an unterschiedlichen Orten miteinander korreliert sind. Spinkorrelationen können beispielsweise über \(C_S(i,j)=\langle \mathbf{S}_i\cdot\mathbf{S}_j\rangle\) untersucht werden. Das Vorzeichen und die räumliche Struktur liefern Informationen über ferro- oder antiferromagnetische Tendenzen. Korrelationsfunktionen lassen sich in den Impulsraum transformieren. Ein Dichtestrukturfaktor kann beispielsweise als \(S_c(q)=\frac{1}{L}\sum_{j,k}e^{iq(j-k)}\langle(n_j-\bar{n})(n_k-\bar{n})\rangle\) definiert werden. Ausgeprägte Maxima bei bestimmten Wellenvektoren liefern Hinweise auf periodische Ordnungsstrukturen. Die Einteilchen-Spektralfunktion beschreibt, bei welchen Energien und Impulsen elektronische Anregungen auftreten. Sie stellt damit eine wichtige Verbindung zwischen Vielteilchentheorie und spektroskopischen Experimenten dar. In korrelierten Systemen können sich Spektren stark von einfachen Bandstrukturen unterscheiden. Wechselwirkungen erzeugen Renormierungen, zusätzliche Strukturen und Energielücken. Ein Quantensystem muss nicht im thermischen Gleichgewicht untersucht werden. Wird ein Modellparameter plötzlich verändert, spricht man häufig von einem Quantum Quench. Ein System kann beispielsweise zunächst im Grundzustand eines Hamiltonoperators \(H_0\) präpariert werden. Ab dem Zeitpunkt \(t=0\) entwickelt es sich unter einem anderen Hamiltonoperator \(H_1\): \(|\psi(t)\rangle=e^{-iH_1t}|\psi(0)\rangle\) Damit lässt sich untersuchen, wie schnell sich Korrelationen ausbreiten und ob sich das System langfristig einem thermischen Zustand annähert. Lokale Änderungen wirken sich nicht augenblicklich auf ein beliebig weit entferntes lokales Gebiet aus. Bei kurzreichweitigen Gittermodellen existieren effektive Grenzen für die Geschwindigkeit, mit der sich Information und Korrelationen ausbreiten. Quantensimulatoren erlauben es, diese Dynamik räumlich und zeitlich direkt zu verfolgen. Wird ein System periodisch angeregt, kann seine Dynamik mit Floquet-Theorie beschrieben werden. Ein zeitperiodischer Hamiltonoperator erfüllt \(H(t+T)=H(t)\) mit der Periodendauer \(T\). Durch geeignetes periodisches Driving können effektive Hopping-Amplituden oder Wechselwirkungen verändert werden. Damit entsteht die Möglichkeit des Hamiltonian Engineering: Eigenschaften eines Quantensystems werden dynamisch erzeugt, statt ausschließlich durch seine statischen Materialparameter bestimmt zu sein. Das erweiterte Hubbard-Modell verbindet Grundlagenphysik direkt mit Quantentechnologie. Analoge Quantensimulatoren versuchen, seine Dynamik physisch nachzubilden. Digitale Quantencomputer zerlegen denselben Hamiltonoperator in Quantengatter. Damit dient ein einziges theoretisches Modell als gemeinsame Sprache für Festkörperphysik, Atomphysik und Quanteninformatik. Ein wesentliches langfristiges Ziel besteht darin, komplexe Materialeigenschaften besser aus mikroskopischen Wechselwirkungen abzuleiten. Besonders relevant sind Materialien, in denen Magnetismus, Ladungsordnung und Supraleitung eng miteinander konkurrieren. Quantencomputer könnten klassische Methoden dabei ergänzen, wenn künftig Systemgrößen und Genauigkeiten erreicht werden, die klassisch nicht mehr effizient zugänglich sind. In naher Zukunft sind hybride Verfahren besonders relevant. Dabei übernimmt ein Quantenprozessor bestimmte Zustandpräparationen oder Messungen, während klassische Rechner Optimierung, Datenanalyse oder Teile der Vielteilchenrechnung durchführen. Hubbard-Modelle eignen sich hervorragend zur Erforschung solcher Arbeitsaufteilungen, weil zahlreiche klassische Referenzmethoden existieren. Eine besonders weitreichende Entwicklung besteht darin, nicht nur reale Materialien nachzubilden, sondern Hamiltonoperatoren gezielt zu konstruieren. In optischen Gittern, Rydberg-Arrays oder anderen Plattformen können Parameterbereiche realisiert werden, die in natürlichen Festkörpern schwer oder überhaupt nicht zugänglich sind. Das erweiterte Hubbard-Modell wird damit von einem Beschreibungsmodell zu einem Konstruktionsprinzip künstlicher Quantenmaterie. Auch eine erweiterte Hubbard-Beschreibung bleibt eine Vereinfachung. Ein reales Material enthält Atomkerne, viele elektronische Orbitale, Phononen, langreichweitige Coulomb-Wechselwirkungen, Defekte, Unordnung und weitere Freiheitsgrade. Ein Hubbard-Modell behält nur diejenigen Freiheitsgrade, die für die untersuchte Niedrigenergiephysik als entscheidend betrachtet werden. Die Aussagekraft des Modells hängt davon ab, wie zuverlässig die effektiven Parameter bestimmt werden. Werte für \(t\), \(U\), \(V\) und weitere Kopplungen können beispielsweise aus elektronischen Strukturrechnungen oder experimenteller Anpassung gewonnen werden. Die Trennung zwischen abgeschirmten und nicht abgeschirmten Wechselwirkungen ist dabei nicht trivial. Ein Einbandmodell kann versagen, wenn mehrere Orbitale in einem ähnlichen Energiebereich liegen. Ebenso kann die Vernachlässigung von Phononen problematisch sein, wenn Elektron-Gitter-Kopplungen eine zentrale Rolle spielen. In solchen Fällen muss das Modell erweitert oder durch andere effektive Hamiltonoperatoren ersetzt werden. Das Modell ist leicht aufzuschreiben, aber schwer exakt zu lösen. Hohe Dimensionen, Frustration, lange Wechselwirkungsreichweiten und Fermionen können den klassischen Rechenaufwand massiv erhöhen. Auch Quantencomputer lösen dieses Problem nicht automatisch. Heutige Geräte besitzen begrenzte Qubitzahlen, endliche Kohärenzzeiten, Gatterfehler und erheblichen Messaufwand. Ein entscheidender Fortschritt wird aus der Skalierung kontrollierter Quantensimulatoren auf größere Gitter erwartet. Je größer das System, desto besser können thermodynamische Eigenschaften, kritische Phänomene und langreichweitige Korrelationen untersucht werden. Entscheidend ist jedoch nicht nur die Teilchenzahl. Präparation, Homogenität, Temperatur und Messgenauigkeit müssen gleichzeitig verbessert werden. Viele besonders interessante Hubbard-Phasen entstehen erst bei sehr niedrigen effektiven Temperaturen. Das Erreichen dieser Temperaturbereiche gehört daher zu den wesentlichen experimentellen Herausforderungen bei ultrakalten fermionischen Atomen. Fortschritte bei Kühlung und Entropieumverteilung könnten den direkten Zugang zu komplexeren magnetischen und möglicherweise paarungsdominierten Zuständen erheblich erweitern. Dipolare Atome, Moleküle und Rydberg-Systeme eröffnen die Möglichkeit, nicht nur ein einzelnes \(V\), sondern ganze räumliche Wechselwirkungsprofile zu kontrollieren. Damit könnten Hubbard-Modelle mit \(V_{ij}=V(|\mathbf{r}_i-\mathbf{r}_j|)\) systematisch untersucht werden. Das erweitert die zugängliche Physik von lokaler Korrelation zu kontrollierter langreichweitiger Vielteilchenwechselwirkung. Ein wichtiger Forschungsweg verbindet detaillierte elektronische Strukturrechnungen mit effektiven Vielteilchenmodellen. Zunächst werden relevante Orbitale und Wechselwirkungen aus einer realistischen Materialbeschreibung gewonnen. Anschließend wird das resultierende korrelierte Modell mit spezialisierten Vielteilchenverfahren untersucht. Quantentechnologien könnten langfristig die besonders schwierigen Schritte innerhalb dieser Prozesskette übernehmen. Der entscheidende Technologiesprung wäre ein fehlertoleranter Quantencomputer mit ausreichend vielen logischen Qubits. Dann könnten tiefe Schaltungen, präzise Zeitentwicklungen und Algorithmen wie Quantum Phase Estimation wesentlich realistischer eingesetzt werden. Hubbard-artige Hamiltonoperatoren gelten als natürliche Kandidaten für solche Anwendungen, weil sie sowohl wissenschaftlich relevant als auch klassisch schwierig sind. Die zukünftige Hubbard-Forschung wird sich nicht ausschließlich auf bekannte Gleichgewichtsphasen konzentrieren. Programmierbare Quantensysteme ermöglichen zeitabhängige Wechselwirkungen, künstliche Eichfelder, periodisches Driving und kontrollierte Zustandspräparation. Damit können Zustände untersucht werden, die in einem gewöhnlichen thermischen Material nicht ohne Weiteres existieren. Das erweitert die Aufgabe von der Erklärung vorhandener Quantenmaterie zur gezielten Erzeugung neuer Quantenmaterie. Das erweiterte Hubbard-Modell gehört zu den wichtigsten Minimalmodellen stark korrelierter Quantenmaterie. Seine Stärke liegt nicht in einer möglichst vollständigen Beschreibung jedes Details eines realen Materials, sondern in der präzisen Erfassung zentraler konkurrierender Mechanismen. Das gewöhnliche Hubbard-Modell enthält bereits zwei entscheidende Prozesse: die Bewegung von Teilchen mit der Energieskala \(t\) und die lokale Wechselwirkung \(U\). Das erweiterte Hubbard-Modell ergänzt insbesondere nichtlokale Wechselwirkungen wie \(V\). Dadurch wird die räumliche Verteilung der Ladung zu einem eigenständigen Bestandteil der Vielteilchenphysik. Aus dieser scheinbar kleinen Erweiterung entsteht ein reiches Phasendiagramm. Magnetische Zustände, Mott-Isolatoren, Ladungsdichtewellen, Bindungsordnungen, metallische Phasen und supraleitende Korrelationen können miteinander konkurrieren. Gittergeometrie, Dimensionalität, Füllung und zusätzliche Wechselwirkungsterme entscheiden darüber, welcher Zustand energetisch bevorzugt wird. Das Modell zeigt damit ein Grundprinzip der Quantenphysik vieler Teilchen: Komplexes kollektives Verhalten benötigt keinen komplizierten mikroskopischen Hamiltonoperator. Schon wenige lokale Regeln können im Zusammenspiel vieler Quantenteilchen völlig unterschiedliche makroskopische Zustände hervorbringen. Für klassische Computer wird die exakte Beschreibung solcher Systeme mit zunehmender Größe schnell schwierig. Verfahren wie DMRG, Quanten-Monte-Carlo, DMFT, funktionale Renormierungsgruppe und Tensor-Netzwerke lösen wichtige Teilbereiche, besitzen jedoch jeweils eigene Grenzen. Genau deshalb hat das erweiterte Hubbard-Modell eine besondere Bedeutung für die Quantentechnologie. Ultrakalte Atome, polare Moleküle, Rydberg-Systeme, Quantenpunkte und programmierbare Quantenprozessoren eröffnen neue Möglichkeiten, Hubbard-artige Vielteilchensysteme experimentell oder digital nachzubilden. Der entscheidende wissenschaftliche Wert liegt dabei nicht allein darin, einen bekannten Hamiltonoperator auf neuer Hardware auszuführen. Ziel ist es, jene Parameterbereiche, Systemgrößen und Dynamiken zugänglich zu machen, bei denen klassische Verfahren an praktische Grenzen stoßen. Das erweiterte Hubbard-Modell steht somit an einer Schnittstelle zwischen theoretischer Festkörperphysik, Quantenmaterialforschung und Quanteninformation. Es ist zugleich Grundlagenmodell, numerische Herausforderung, experimentelles Designprinzip und Benchmark für zukünftige Quantentechnologien. Die wissenschaftliche Literatur zum erweiterten Hubbard-Modell reicht von den grundlegenden Arbeiten zur lokalen Elektronenkorrelation bis zu modernen Untersuchungen nichtlokaler Wechselwirkungen, konkurrierender Quantenphasen und quantentechnologischer Simulationen. Für eine wissenschaftliche Abhandlung ist es sinnvoll, zwischen historischer Primärliteratur, Arbeiten zum Phasendiagramm, methodischer Literatur und experimentell beziehungsweise quantentechnologisch orientierten Veröffentlichungen zu unterscheiden. Monographien sind für das erweiterte Hubbard-Modell vor allem dort wichtig, wo die Abhandlung über einzelne Forschungsergebnisse hinaus die theoretischen Grundlagen erklären muss. Besonders relevant sind Werke zur Hubbard-Physik, eindimensionalen Vielteilchensystemen, Quantenmagnetismus, Vielteilchentheorie, Quantenphasenübergängen sowie zur numerischen und quantenalgorithmischen Behandlung wechselwirkender Fermionen. Online-Ressourcen erfüllen im Zusammenhang mit dem erweiterten Hubbard-Modell unterschiedliche Funktionen. Preprint-Server ermöglichen den Zugriff auf aktuelle Forschung, Fachverlage stellen begutachtete Journalversionen bereit, Literaturdatenbanken unterstützen die systematische Recherche und wissenschaftliche Softwareplattformen helfen bei der numerischen Reproduktion beziehungsweise Weiterentwicklung von Modellrechnungen.Bravyi-Kitaev-Transformation
Digitale Zeitentwicklung
Variational Quantum Eigensolver
Quantum Phase Estimation
Das erweiterte Hubbard-Modell als Benchmark
Ein anspruchsvolles Referenzproblem
Klassischer Vergleich
Hardwarevalidierung
Experimentelle Plattformen
Ultrakalte Fermionen
Polare Moleküle und dipolare Atome
Rydberg-Systeme
Quantenpunkt-Arrays
Supraleitende Quantenschaltungen
Messbare Größen und Observablen
Teilchendichte
Doppelte Besetzung
Ladungskorrelationen
Spinkorrelationen
Strukturfaktoren
Spektralfunktion
Nichtgleichgewichtsphysik und Quantendynamik
Quenches
Korrelationsausbreitung
Periodisches Driving
Bedeutung für aktuelle Quantentechnologien
Quantensimulation
Erforschung von Quantenmaterialien
Hybride Algorithmen
Programmierbare Quantenmaterie
Grenzen des erweiterten Hubbard-Modells
Reduziertes Modell
Parameterbestimmung
Fehlende Freiheitsgrade
Numerische Grenzen
Forschungsperspektiven
Größere Quantensimulatoren
Kühlere fermionische Systeme
Langreichweitige Wechselwirkungen
Verbindung zu Ab-initio-Verfahren
Fehlertolerantes Quantencomputing
Nichtgleichgewicht und neue Phasen
Fazit
Anhang
Wissenschaftliche Zeitschriften und Artikel
Grundlegende Primärliteratur zum Hubbard-Modell
Spezialisierte Arbeiten zum Phasendiagramm des erweiterten Hubbard-Modells
Aktuelle Spezialliteratur zu nichtlokalen Wechselwirkungen und Supraleitung
Hintergrundliteratur zu numerischen Vielteilchenmethoden
Quantensimulation und experimentelle Hubbard-Physik
Bücher und Monographien
Monographien zum Hubbard-Modell und zu eindimensionalen Systemen
Hintergrundliteratur zu stark korrelierten Elektronen
Hintergrundliteratur zu Quantenphasen und kritischem Verhalten
Standardwerke zur Quanteninformation
Online-Ressourcen und Datenbanken
Vorlesungsnotizen und Monographie-nahe Ressourcen
Fachjournale und Verlage
Lern- und Forschungsplattformen
