Higher Berry curvature, second Chern numbers and magnetoelectric coupling in crystalline insulators

Technische Zusammenfassung: Höhere Berry-Krümmung, zweite Chern-Zahlen und magnetoelektrische Kopplung in kristallinen Isolatoren

Problemstellung
Topologische Invarianten sind essenziell für die Klassifizierung von Quanten-Vielteilchen-Grundzuständen. Während die Berry-Krümmung (eine 2-Form) und die assoziierte Chern-Zahl erfolgreich Familien von nuldimensionalen Grundzuständen (z. B. Bloch-Zustände in der Brillouin-Zone) charakterisieren, deuten jüngste theoretische Entwicklungen auf eine Verallgemeinerung für räumlich ausgedehnte Systeme hin. In diesem Rahmen werden Familien von kurzreichweitig verschränkten Zuständen in $d$ räumlichen Dimensionen durch Differentialformen höherer Ordnung (z. B. eine 3-Form für 1D-Systeme) charakterisiert. Speziell für eindimensionale Systeme, die durch einen dreidimensionalen Raum parametrisiert sind, ist die relevante Invariante die Dixmier–Douady–Kapustin–Spodyneiko (DDKS)-Zahl $\nu$, die aus einer höheren Berry-Krümmung abgeleitet wird.

Es besteht eine signifikante Lücke in der expliziten, manifest quantisierten Berechnung dieser höheren Invarianten für generische wechselwirkende oder Free-Fermion-Systeme. Während analytische Ausdrücke für spezifische Free-Fermion-Modelle (wie das vierdimensionale Dirac-Modell) existieren, beruhen numerische Methoden zur Berechnung der zweiten Chern-Zahl ($ch_2$) in vier Dimensionen typischerweise auf Gitterfeldtheorie oder Realraum-Magnetfeldantworten. Diese Standardmethoden lassen oft eine „manifeste Quantisierung“ vermissen, was bedeutet, dass sie erst im Kontinuumslimit zu ganzzahligen Werten konvergieren oder eine extensive Verfeinerung des Gitters erfordern, anstatt konstruktionsbedingt ganze Zahlen zurückzugeben. Die physikalische Interpretation der DDKS-Zahl in drei räumlichen Dimensionen, insbesondere ihre Beziehung zur orbitalen magnetoelektrischen Kopplung und dem Chern-Simons-Axion-Winkel ($\theta_{CS}$), bleibt eine offene Frage.

Methodik
Die Autoren adressieren diese Probleme, indem sie das Gittermodell eines vierdimensionalen (4D) Chern-Isolators als eine Familie von translationsinvarianten eindimensionalen (1D) Ketten umformulieren, die durch eine dreidimensionale Brillouin-Zone ($T^3$) parametrisiert sind.

1. Modellumformulierung: Ausgehend vom 4D-Dirac-Hamiltonian führen die Autoren eine partielle Fourier-Transformation auf eine Impulskomponente ($k_4$) durch, was zu einer Familie von 1D-Hamiltonianen $\{H_\lambda\}$ führt, die von Parametern $\lambda = (\lambda_1, \lambda_2, \lambda_3)$ abhängen. Diese Hamiltonianen beschreiben 1D-Ketten mit zwei Orbitalen pro Einheitszelle.
2. Numerisches Framework (iMPS): Die Grundzustände dieser 1D-Ketten werden durch unendliche Matrix-Produkt-Zustände (iMPS) approximiert. Die Autoren nutzen den Algorithmus, der in Ref. [13] vorgeschlagen wurde, welcher die Fukui-Hatsugai-Suzuki-Methode zur Berechnung höherer Berry-Krümmungen generalisiert.
3. Diskrete höhere Berry-Krümmung:
   • Die Parameter-Raum $T^3$ wird in ein Gitter diskretisiert und in Tetraeder trianguliert.
   • Für jeden Knoten wird ein iMPS-Grundzustand via iDMRG ermittelt.
   • Übergangsfunktionen zwischen benachbarten iMPS-Tensoren werden unter Verwendung der dominanten Eigenvektoren gemischter Transfermatrizen berechnet.
   • Eine gauge-invariante Phase $\phi(\Delta_2)$ wird für trianguläre Flächen definiert, gewichtet durch das Schmidt-Spektrum, um die variierende Bindungsdimension zu handhaben.
   • Die höhere Berry-Krümmung $F(\Delta_3)$ für ein Tetraeder ist als orientierte Summe dieser Phasen über dessen Rand definiert.
   • Die DDKS-Zahl $\nu$ wird durch Summation von $F(\Delta_3)$ über die gesamte Triangulation von $T^3$ erhalten, was eine manifest quantisierte ganze Zahl liefert.
4. Vergleichende Analyse: Die Autoren berechnen die zweite Chern-Zahl $ch_2$ für dasselbe Modell unter Verwendung eines Standard-Gitterfeldtheorie-Ansatzes (Ref. [34]) und vergleichen die resultierenden Phasendiagramme.
5. Magnetoelektrische Kopplung: Die Studie untersucht die Beziehung zwischen der DDKS-Zahl und dem Chern-Simons-Axion-Winkel $\theta_{CS}$, der die orbitale magnetoelektrische Kopplung in 3D steuert. Sie vergleichen die Variation von $\theta_{CS}$ (abgeleitet aus der zweiten Chern-Form) mit der Variation einer „höheren Berry-Phase“ $\theta_{HBP}$ (abgeleitet aus der iMPS-Konstruktion) über einen Pumping-Zyklus.

Wesentliche Beiträge und Ergebnisse

• Manifest quantisierte Berechnung von $ch_2$: Die Autoren zeigen, dass die DDKs-Zahl $\nu$, berechnet über den iMPS-basierten Algorithmus der höheren Berry-Krümmung, ein topologisches Phasendiagramm für den 4D-Chern-Isolator erzeugt, das exakt kongruent zum Phasendiagramm der zweiten Chern-Zahl $ch_2$ ist. Diese Übereinstimmung gilt über verschiedene Massenregime ($m$) hinweg, einschließlich Regionen, in denen Standard-Gitterfeldtheorie-Berechnungen von $ch_2$ unter erheblichen endlicher-Gitter-Fehlern leiden. Die iMPS-Methode liefert exakte ganze Zahlen selbst auf groben Gittern (z. B. $4^3$), während die Gitterfeldmethode wesentlich feinere Gitter ($10^4$ Punkte) benötigt, um zu konvergieren.
• Theoretische Verbindung via Künneth-Formel: Das Paper etabliert eine rigorose Verbindung zwischen der 4D-zweiten Chern-Zahl und der 3D-DDKS-Zahl mittels der Künneth-Formel. Durch die Zerlegung des 4D-Torus $T^4$ als $S^1 \times T^3$ zeigen die Autoren, dass die zweite Chern-Klasse auf $T^4$ dem Produkt des Erzeugers von $H^1(S^1, \mathbb{Z})$ und der Klasse in $H^3(T^3, \mathbb{Z})$ entspricht, welche durch $\nu$ gemessen wird. Folglich gilt für diese spezifische Zerlegung $ch_2 = \nu$.
• Magnetoelektrische Kopplung im 4D-Dirac-Modell: Für das 4D-Dirac-Modell finden die Autoren, dass die Änderung des Chern-Simons-Axion-Winkels $\Delta\theta_{CS}$ entlang eines Pumping-Zyklus präzise durch die Änderung der höheren Berry-Phase $\Delta\theta_{HBP}$ reproduziert wird. Dies deutet auf eine direkte Korrespondenz zwischen der höheren Berry-Phase und der magnetoelektrischen Antwort in Systemen hin, in denen die ersten Chern-Zahlen der BZ-Schnitte verschwinden.
• Gegenbeispiel im Hopf-Isolator: Die Autoren testen diese Korrespondenz am Hopf-Isolator, einem Zwei-Band-Modell mit einer nicht-trivialen Hopf-Invariante. In diesem Fall ist, obwohl der Chern-Simons-Winkel durch die Hopf-Invariante quantisiert ist, die berechnete höhere Berry-Phase $\theta_{HBP}$ weder quantisiert noch gleich $\theta_{CS}$. Die Autoren führen diese Diskrepanz auf die „delikate“ Natur der Hopf-Invariante zurück, die streng genommen nur für Zwei-Band-Modelle im Impulsraum wohldefiniert ist – eine Bedingung, die in der Realraum-iMPS-Beschreibung des thermodynamischen Limits möglicherweise nicht erhalten bleibt.

Bedeutung und Ansprüche
Das Paper beansprucht, die erste explizite, manifest quantisierte numerische Berechnung der zweiten Chern-Zahl für ein 4D-Gittermodell unter Verwendung höherer Berry-Krümmung und iMPS zu liefern. Dies validiert den theoretischen Rahmen, wonach höhere Berry-Krümmungen als Sonden für höhergradige topologische Invarianten in Familien räumlich ausgedehnter Zustände dienen können.

Die Autoren ordnen ihre Ergebnisse in den Kontext der Quantengeometrie und der topologischen Eigenschaften ein und merken an, dass während die Korrespondenz zwischen höheren Berry-Phasen und magnetoelektrischer Kopplung für das 4D-Dirac-Modell gilt, sie nicht universell ist, wie der Fall des Hopf-Isolators zeigt. Sie betonen, dass ihre Methode besonders für wechselwirkende Systeme geeignet ist, in denen die Single-Particle-Bloch-Beschreibungen versagen, wenngleich sie sich derzeit auf Free-Fermion-Benchmarks konzentrieren. Die Arbeit hebt das Potenzial synthetischer Dimensionen und quasi-periodischer Anregung hervor, um diese höherdimensionalen topologischen Phänomene in experimentellen Plattformen zu realisieren und zu untersuchen.

Das Paper schließt bescheiden mit dem Hinweis, dass die Beziehung zwischen höherer Berry-Krümmung und physikalischer Antwort in drei Dimensionen eine offene Frage bleibt, insbesondere hinsichtlich der Universalität des Zusammenhangs mit $\theta_{CS}$. Es legt nahe, dass höhere Berry-Phasen geometrische Informationen erfassen könnten, die sich von Standard-Chern-Simons-Antworten unterscheiden, was weitere Untersuchungen zu ihrer Rolle in der nichtlinearen Transportphysik und der Grundzustandsphysik rechtfertigt.