Ein Projected Entangled Pair Operator, kurz PEPO, ist eine Tensor-Netzwerk-Darstellung für Operatoren in quantenmechanischen Vielteilchensystemen. PEPOs sind insbesondere für zwei- und höherdimensionale Gittersysteme von Bedeutung. Sie übertragen die Grundidee der Projected Entangled Pair States, kurz PEPS, von Quantenzuständen auf Operatoren.

Während ein PEPS die Koeffizienten eines Vielteilchenzustands durch ein Netzwerk lokaler Tensoren beschreibt, repräsentiert ein PEPO die Matrixelemente eines Vielteilchenoperators. Jeder lokale Tensor besitzt dabei physikalische Ein- und Ausgangsindizes sowie virtuelle Indizes, über die er mit benachbarten Tensoren verbunden ist. Durch die Kontraktion aller virtuellen Indizes entsteht der vollständige Operator.

Ein allgemeiner Operator auf einem Gitter kann formal als

\(O = \sum_{\{i_r,j_r\}} C_{\{i_r,j_r\}} |\{i_r\}\rangle \langle\{j_r\}|\)

geschrieben werden. Die Aufgabe eines PEPO besteht darin, den Koeffiziententensor \(C_{\{i_r,j_r\}}\) nicht explizit zu speichern, sondern ihn als Kontraktion vieler kleiner lokaler Tensoren darzustellen.

Diese Zerlegung ist entscheidend, weil der vollständige Operatorraum eines Vielteilchensystems exponentiell mit der Zahl der Freiheitsgrade wächst. PEPOs versuchen, die tatsächlich vorhandene lokale Struktur, Korrelation und Operatorverschränkung auszunutzen und dadurch eine kompaktere Beschreibung zu ermöglichen.

Bedeutung für die Quantentechnologie

PEPOs spielen eine Rolle bei der klassischen Simulation von Quantenmaterie, thermischen Zuständen, Zeitentwicklungsoperatoren, offenen Quantensystemen und zweidimensionalen Quantengittern. Sie sind besonders interessant, wenn ein Problem zwar auf einem exponentiell großen Hilbertraum definiert ist, der relevante Operator aber eine begrenzte lokale oder korrelative Struktur besitzt.

Damit gehören PEPOs zu einer größeren Familie von Tensor-Netzwerk-Methoden, die versuchen, nicht die gesamte mathematisch mögliche Komplexität eines Quantensystems abzubilden, sondern dessen physikalisch relevante Struktur.

Mathematische und physikalische Grundlagen

Tensoren und Tensor-Netzwerke

Ein Tensor ist eine mehrdimensionale Verallgemeinerung von Skalar, Vektor und Matrix. Ein Skalar besitzt keinen freien Index, ein Vektor einen Index, eine Matrix zwei Indizes und ein allgemeiner Tensor entsprechend mehr.

Ein lokaler Tensor eines zweidimensionalen PEPO kann beispielsweise sechs Indizes besitzen: zwei physikalische und vier virtuelle. Bei einem quadratischen Gitter entsprechen die vier virtuellen Indizes typischerweise den Verbindungen nach links, rechts, oben und unten.

Tensor-Netzwerke stellen solche Tensoren grafisch als Knoten dar. Indizes erscheinen als Linien. Werden zwei Tensoren über einen gemeinsamen Index verbunden, wird über diesen Index summiert. Dieser Vorgang heißt Tensor-Kontraktion.

Für zwei Tensoren \(A_{i\alpha}\) und \(B_{\alpha j}\) lautet eine einfache Kontraktion beispielsweise

\(C_{ij} = \sum_{\alpha} A_{i\alpha} B_{\alpha j}\)

Ein PEPO erweitert dieses Prinzip auf ein gesamtes mehrdimensionales Gitter.

Hilberträume von Vielteilchensystemen

Die zentrale Schwierigkeit der Quanten-Vielteilchenphysik ist das exponentielle Wachstum des Hilbertraums. Besitzt jeder von N Gitterplätzen eine lokale Hilbertraumdimension d, dann hat der vollständige Hilbertraum die Dimension

\(\dim(\mathcal{H}) = d^N\)

Für ein System aus Spin-ein-halb-Freiheitsgraden gilt beispielsweise \(d=2\). Schon bei einer moderaten Anzahl von Gitterplätzen wird eine vollständige Speicherung eines allgemeinen Zustandsvektors auf klassischen Rechnern praktisch unmöglich.

Für Operatoren verschärft sich das Problem. Ein allgemeiner linearer Operator auf einem Hilbertraum der Dimension \(d^N\) kann als Matrix mit

\(d^N \times d^N = d^{2N}\)

Matrixelementen aufgefasst werden. Eine naive Darstellung skaliert damit noch ungünstiger als die Speicherung eines Zustandsvektors.

Tensor-Netzwerke setzen genau an diesem Punkt an. Statt alle Koeffizienten unabhängig voneinander zu speichern, zerlegen sie die globale Struktur in lokale Tensoren und virtuelle Verbindungen.

Verschränkung und Korrelation

Der Erfolg von Tensor-Netzwerken hängt eng mit der Struktur quantenmechanischer Verschränkung zusammen. Physikalisch relevante Vielteilchenzustände sind häufig nicht beliebige Vektoren ihres Hilbertraums. Insbesondere Grundzustände lokaler Hamiltonoperatoren weisen vielfach eine deutlich strukturiertere Verschränkung auf als ein zufälliger Zustand.

In vielen kurzreichweitig wechselwirkenden Systemen folgt die Verschränkungsentropie näherungsweise einem Flächengesetz. Die relevante Verschränkung wird dann stärker durch die Größe der Grenze eines Teilgebiets als durch dessen Volumen bestimmt.

PEPS wurden gerade mit Blick auf solche höherdimensionalen Verschränkungsstrukturen entwickelt. PEPOs übertragen eine verwandte Idee auf Operatoren. Hier ist nicht nur Zustandsverschränkung relevant, sondern auch die Frage, wie komplex die korrelative Struktur eines Operators über verschiedene Raumregionen hinweg ist.

Lokale und globale Operatoren

Ein lokaler Operator wirkt nur auf einen einzelnen oder wenige Gitterplätze. Ein Vielteilchen-Hamiltonoperator kann beispielsweise als Summe lokaler Terme geschrieben werden:

\(H = \sum_{\langle r,s\rangle} h_{r,s} + \sum_r h_r\)

Dabei beschreibt \(h_r\) einen lokalen Term und \(h_{r,s}\) eine Wechselwirkung zwischen Gitterplätzen.

Obwohl der gesamte Hamiltonoperator auf dem vollständigen Hilbertraum wirkt, besitzt er damit eine stark strukturierte Form. Genau diese Lokalität ermöglicht häufig eine kompakte Tensor-Netzwerk-Darstellung.

Von MPS und MPO zu PEPS und PEPO

Matrix Product States

Matrix Product States, MPS, sind Tensor-Netzwerke für eindimensionale Quantensysteme. Jeder Gitterplatz wird durch einen Tensor mit einem physikalischen Index und virtuellen Indizes zu seinen Nachbarn dargestellt.

Für eine offene eindimensionale Kette kann ein Zustand schematisch als

\(|\psi\rangle = \sum_{i_1,\ldots,i_N} A^{[1]i_1} A^{[2]i_2} \cdots A^{[N]i_N} |i_1,\ldots,i_N\rangle\)

geschrieben werden. Die Matrixmultiplikation steht hier für die Kontraktion der virtuellen Indizes.

Die Dimension dieser virtuellen Indizes wird Bond-Dimension genannt. Sie bestimmt maßgeblich, wie komplex die durch das Netzwerk darstellbaren Korrelationen und Verschränkungsstrukturen sein können.

Matrix Product Operators

Ein Matrix Product Operator, MPO, ist das Operatorgegenstück zum MPS. An jedem Gitterplatz besitzt der Tensor zwei physikalische Indizes: einen für den Eingang und einen für den Ausgang des lokalen Operators.

Ein MPO eignet sich unter anderem zur Darstellung lokaler Hamiltonoperatoren, Dichtematrizen und Zeitentwicklungsoperatoren in eindimensionalen Systemen.

Das Grundprinzip ist unmittelbar: Ein MPS kodiert Amplituden eines Zustands, ein MPO kodiert Matrixelemente eines Operators.

Projected Entangled Pair States

PEPS verallgemeinern MPS auf zwei oder mehr räumliche Dimensionen. Auf einem quadratischen zweidimensionalen Gitter besitzt ein innerer PEPS-Tensor typischerweise einen physikalischen und vier virtuelle Indizes.

Dadurch entsteht kein eindimensionales Kettennetzwerk mehr, sondern ein Netzwerk mit Schleifen. Diese Änderung ist physikalisch wertvoll, numerisch aber entscheidend: Die exakte Kontraktion eines zweidimensionalen PEPS ist im Allgemeinen wesentlich schwieriger als die Kontraktion eines MPS.

Der Übergang zum PEPO

Ein PEPO entsteht konzeptionell aus demselben Schritt, der vom MPS zum MPO führt. An die Stelle eines einzigen physikalischen Index eines PEPS treten physikalische Ein- und Ausgangsindizes.

Ein PEPO-Tensor auf einem quadratischen Gitter kann daher schematisch als

\(W^{ij}_{l r u d}\)

geschrieben werden. Dabei sind \(i\) und \(j\) physikalische Indizes, während \(l\), \(r\), \(u\) und \(d\) virtuelle Verbindungen zu benachbarten Tensoren kennzeichnen.

Damit ist ein PEPO im Kern die höherdimensionale Verallgemeinerung eines MPO.

Definition und Struktur eines Projected Entangled Pair Operator

Lokale Tensorstruktur

Betrachten wir ein zweidimensionales Gitter mit lokalen PEPO-Tensoren. Jeder Tensor enthält Informationen darüber, wie ein lokaler Eingabezustand in einen lokalen Ausgabezustand überführt wird und wie diese Operation mit der Umgebung korreliert ist.

Ein allgemeiner PEPO kann symbolisch als

\(O = \sum_{\{i_r,j_r\}} \mathrm{tTr}\left(\prod_r W^{[r]i_rj_r}\right) |\{i_r\}\rangle\langle\{j_r\}|\)

geschrieben werden. Die Bezeichnung \(\mathrm{tTr}\) steht hier für die vollständige Kontraktion der virtuellen Tensorindizes entsprechend der Gittergeometrie.

Entscheidend ist, dass der riesige globale Koeffiziententensor nicht unmittelbar gespeichert wird. Er entsteht erst mathematisch aus dem Netzwerk lokaler Tensoren.

Physikalische Indizes

Da ein PEPO einen Operator und keinen Zustandsvektor beschreibt, benötigt jeder lokale Tensor grundsätzlich einen Eingangs- und einen Ausgangsindex.

Bei einer lokalen Hilbertraumdimension d besitzt der lokale Operatorraum damit formal die Dimension

\(d^2\)

pro Gitterplatz. Man kann einen PEPO deshalb auch als Tensor-Netzwerk in einem lokal vergrößerten Raum betrachten, bei dem die beiden physikalischen Indizes zu einem gemeinsamen Operatorindex zusammengefasst werden.

Virtuelle Indizes

Die virtuellen Indizes haben keine direkte Interpretation als lokal messbare Freiheitsgrade. Sie dienen dazu, Korrelationen und globale Strukturen im Netzwerk zu kodieren.

Auf einem quadratischen Gitter hat ein Tensor im Inneren typischerweise vier virtuelle Nachbarn. Randtensoren besitzen entsprechend weniger Verbindungen.

Je größer die Bond-Dimension dieser virtuellen Indizes ist, desto komplexere Operatorstrukturen können dargestellt werden. Gleichzeitig wachsen Speicherbedarf und Rechenaufwand.

Bond-Dimension

Die Bond-Dimension ist einer der wichtigsten Parameter eines PEPO. Ein Produktoperator ohne nichttriviale räumliche Korrelation kann mit sehr kleiner Bond-Dimension dargestellt werden. Komplexe Vielteilchenoperatoren benötigen im Allgemeinen größere Werte.

Die Bond-Dimension ist jedoch nicht einfach ein Maß für Systemgröße. Ein großes lokales System kann unter günstigen physikalischen Bedingungen einen strukturierten Operator mit moderater Bond-Dimension besitzen. Umgekehrt kann ein kleineres, aber stark korreliertes Problem eine hohe Bond-Dimension erfordern.

Der praktische Kern vieler PEPO-Verfahren besteht deshalb darin, die Bond-Dimension groß genug für eine brauchbare Genauigkeit, aber klein genug für eine beherrschbare Berechnung zu wählen.

Randbedingungen

PEPOs können für offene und periodische Randbedingungen definiert werden. Offene Randbedingungen vereinfachen viele numerische Verfahren, weil Randtensoren weniger virtuelle Verbindungen besitzen.

Periodische Randbedingungen bilden hingegen Systeme ohne offene Kanten besser ab, erhöhen aber häufig den Kontraktionsaufwand.

Auch unendliche Gitter können mit translationsinvarianten Tensor-Netzwerken beschrieben werden. Dabei wird nicht für jeden Gitterplatz ein eigener Tensor gespeichert. Stattdessen wiederholt sich eine endliche Elementarzelle periodisch.

Konstruktion von PEPOs

Produktoperatoren

Der einfachste Fall ist ein Produkt lokaler Operatoren:

\(O = O_1 \otimes O_2 \otimes \cdots \otimes O_N\)

Hier besteht keine zusätzliche virtuelle Korrelation zwischen den lokalen Operatoren. Eine Tensor-Netzwerk-Darstellung kann deshalb mit minimaler Bond-Dimension auskommen.

Interessanter werden Operatoren, die Summen verschiedener Produktstrukturen enthalten. Die virtuellen Indizes übernehmen dann die Aufgabe, festzulegen, welche lokalen Bestandteile zusammen einen globalen Term ergeben.

Lokale Hamiltonoperatoren

Ein typischer Hamiltonoperator eines Gittersystems hat die Form

\(H = \sum_r h_r + \sum_{\langle r,s\rangle} h_{r,s}\)

Obwohl die Zahl der Terme mit der Systemgröße wächst, ist jeder einzelne Term lokal. Dadurch kann der Gesamtoperator oft wesentlich effizienter dargestellt werden als eine allgemeine Matrix des vollständigen Hilbertraums.

Die virtuelle Struktur des PEPO übernimmt dabei gewissermaßen die Buchhaltung darüber, wo ein Wechselwirkungsterm beginnt, über welche Gitterplätze er verläuft und wo er endet.

Langreichweitige Wechselwirkungen

Für langreichweitige Wechselwirkungen wird die Konstruktion anspruchsvoller. Ein Hamiltonoperator kann beispielsweise Terme der Form

\(H = \sum_{r\neq s} J_{rs} A_r B_s\)

enthalten. Wenn die Kopplung \(J_{rs}\) nur langsam mit der Entfernung abfällt, müssen Informationen über größere räumliche Distanzen im Tensor-Netzwerk transportiert werden.

Eine effiziente PEPO-Darstellung hängt dann davon ab, ob die Kopplungsstruktur selbst komprimiert oder durch geeignete Summen einfacherer Funktionen approximiert werden kann.

Faktorisierung und Tensorzerlegung

Ist ein Operator zunächst in einer anderen Darstellung bekannt, können Tensorzerlegungen eingesetzt werden, um eine PEPO-Struktur zu erzeugen. Ein grundlegendes Werkzeug ist die Singularwertzerlegung.

Für eine Matrix M gilt

\(M = U S V^\dagger\)

Die Singularwerte in S geben an, welche Komponenten einer Zerlegung besonders wichtig sind. Kleine Singularwerte können unter kontrolliertem Genauigkeitsverlust verworfen werden. Entsprechende Ideen treten bei der Zerlegung und Kompression höherdimensionaler Tensoren immer wieder auf.

Exakte und approximative PEPOs

Nicht jeder relevante Operator besitzt bei vorgegebener kleiner Bond-Dimension eine exakte PEPO-Darstellung. In vielen praktischen Anwendungen wird deshalb eine Approximation gesucht.

Das Ziel besteht dann darin, einen Operator \(\widetilde{O}\) zu finden, der den ursprünglichen Operator O in den für die konkrete Anwendung wichtigen Eigenschaften hinreichend gut wiedergibt.

Entscheidend ist dabei die Fehlerdefinition. Eine Approximation kann hinsichtlich einer globalen Norm mäßig sein und trotzdem lokale Erwartungswerte sehr genau reproduzieren. Umgekehrt garantiert eine gute Beschreibung einzelner Observablen nicht automatisch eine global genaue Darstellung aller Matrixelemente.

PEPOs für Dichteoperatoren und gemischte Quantenzustände

Der Dichteoperator

Ein reiner Zustand wird durch einen Zustandsvektor beschrieben. Ein gemischter Zustand benötigt dagegen einen Dichteoperator. Allgemein gilt

\(\rho = \sum_k p_k |\psi_k\rangle\langle\psi_k|\)

mit

\(p_k \geq 0\)

und

\(\sum_k p_k = 1\)

Ein physikalischer Dichteoperator erfüllt insbesondere

\(\rho^\dagger = \rho\)

sowie

\(\rho \geq 0\)

und

\(\mathrm{Tr}(\rho)=1\)

Da ein Dichteoperator selbst ein Vielteilchenoperator ist, kann er grundsätzlich als PEPO dargestellt werden.

Lokale Observablen

Ist ein Dichteoperator als PEPO gegeben, lässt sich der Erwartungswert einer Observable A über

\(\langle A\rangle = \mathrm{Tr}(\rho A)\)

bestimmen.

Formal ist diese Gleichung einfach. Numerisch muss jedoch das entsprechende zweidimensionale Tensor-Netzwerk kontrahiert werden. Genau hier liegt eine der zentralen Schwierigkeiten der Methode.

Reduzierte Dichteoperatoren

Für lokale Eigenschaften wird häufig nicht der vollständige Dichteoperator benötigt. Interessiert nur ein Teilgebiet A, kann über die übrigen Freiheitsgrade B partiell gespurt werden:

\(\rho_A = \mathrm{Tr}_B(\rho)\)

Die lokale Struktur eines PEPO kann solche Berechnungen unterstützen, weil viele physikalische Fragestellungen ohnehin nur kleine Teilsysteme oder lokale Korrelationsfunktionen betreffen.

Positivität als Herausforderung

Ein beliebiges Tensor-Netzwerk mit zwei physikalischen Indizes pro Gitterplatz stellt nicht automatisch einen positiven Dichteoperator dar. Dies ist ein wichtiger Unterschied zwischen einer formalen Operatorapproximation und einer garantiert physikalischen Dichtematrix.

Ein Ansatz besteht darin, Positivität bereits durch die Parametrisierung sicherzustellen. Beispielsweise kann ein Dichteoperator über eine Faktorisierung

\(\rho = X X^\dagger\)

konstruiert werden.

Eine solche Struktur garantiert Positivität, kann aber den numerischen Aufwand erhöhen. Daher muss zwischen mathematisch garantierten Eigenschaften und praktischer Effizienz abgewogen werden.

Thermische Zustände und statistische Quantenmechanik

Gibbs-Zustände

Ein besonders wichtiger Einsatzbereich von PEPOs sind thermische Quantenzustände. Für einen Hamiltonoperator H und die inverse Temperatur \(\beta\) lautet der Gibbs-Zustand

\(\rho_\beta = \frac{e^{-\beta H}}{Z}\)

mit der Zustandssumme

\(Z = \mathrm{Tr}(e^{-\beta H})\)

und

\(\beta = \frac{1}{k_B T}\)

Der nicht normierte thermische Operator

\(e^{-\beta H}\)

ist selbst ein Vielteilchenoperator und kann daher durch einen PEPO approximiert werden.

Imaginärzeitentwicklung

Eine verbreitete Idee besteht darin, bei unendlich hoher Temperatur zu beginnen. Dort ist der Dichteoperator proportional zur Identität. Anschließend wird die inverse Temperatur schrittweise erhöht.

Formal kann geschrieben werden:

\(e^{-\beta H} = \left(e^{-\Delta\beta H}\right)^M\)

mit

\(\beta = M\Delta\beta\)

Der Operator für einen kleinen Schritt kann mithilfe lokaler Zerlegungen approximiert und sukzessive auf das Tensor-Netzwerk angewendet werden.

Suzuki-Trotter-Zerlegung

Kann der Hamiltonoperator in mehrere Teile zerlegt werden,

\(H = H_A + H_B\)

deren lokale Terme innerhalb der jeweiligen Gruppe geeignet gewählt sind, lässt sich für einen kleinen Schritt beispielsweise die Approximation

\(e^{-\Delta\beta(H_A+H_B)} \approx e^{-\Delta\beta H_A} e^{-\Delta\beta H_B}\)

verwenden.

Eine symmetrische Variante höherer Genauigkeit lautet

\(e^{-\Delta\beta(H_A+H_B)} \approx e^{-\frac{\Delta\beta}{2}H_A} e^{-\Delta\beta H_B} e^{-\frac{\Delta\beta}{2}H_A}\)

Für endliche Schrittweiten entsteht ein Trotter-Fehler. Kleinere Schritte reduzieren diesen Fehler, erhöhen jedoch die Zahl der notwendigen Operationen.

Wachstum der Bond-Dimension

Bei der wiederholten Anwendung lokaler Operatoren wächst typischerweise die Bond-Dimension des Netzwerks. Ohne Gegenmaßnahmen würde die Rechnung schnell unbeherrschbar.

Deshalb folgt auf die Operatoranwendung gewöhnlich eine Kompression. Die Bond-Dimension wird auf einen vorgegebenen Wert reduziert, wobei versucht wird, die physikalisch wichtigsten Informationen zu erhalten.

Die Genauigkeit einer thermischen PEPO-Rechnung wird damit im Wesentlichen durch mehrere Faktoren bestimmt: Trotter-Schrittweite, Bond-Dimension, Kontraktionsgenauigkeit und Kompressionsverfahren.

Thermodynamische Observablen

Ist der thermische Zustand hinreichend genau bestimmt, können daraus thermodynamische Größen berechnet werden. Die mittlere Energie lautet beispielsweise

\(E = \mathrm{Tr}(\rho_\beta H)\)

und die Wärmekapazität folgt aus der Temperaturabhängigkeit der Energie.

Auch Korrelationsfunktionen und Ordnungsparameter lassen sich untersuchen. Damit ermöglichen PEPO-basierte Verfahren den Zugang zu endlichen Temperaturen in zweidimensionalen Quantengittern, einem Bereich, in dem viele andere numerische Methoden erhebliche Schwierigkeiten haben.

Anwendung von PEPOs auf PEPS

Operatorwirkung auf einen Zustand

Eine besonders natürliche Situation entsteht, wenn ein Zustand als PEPS und ein Operator als PEPO vorliegt. Formal gilt

\(|\phi\rangle = O|\psi\rangle\)

Lokale physikalische Indizes des PEPO werden dabei mit den entsprechenden physikalischen Indizes des PEPS kontrahiert.

Aus jedem Paar aus PEPS- und PEPO-Tensor entsteht ein neuer lokaler Tensor. Dessen virtuelle Dimension ist im Allgemeinen größer als zuvor.

Wachstum der virtuellen Dimension

Besitze das ursprüngliche PEPS eine Bond-Dimension D und der angewendete PEPO eine Bond-Dimension \(\chi\), dann kann die direkte Operatoranwendung virtuelle Dimensionen erzeugen, die sich schematisch wie

\(D' = D\chi\)

verhalten.

Bei wiederholten Anwendungen wächst die Darstellung daher sehr schnell. Praktische Algorithmen müssen das resultierende Netzwerk regelmäßig komprimieren.

Erwartungswerte

Für einen Zustand \(|\psi\rangle\) und einen Operator O lautet der Erwartungswert

\(\langle O\rangle = \frac{\langle\psi|O|\psi\rangle}{\langle\psi|\psi\rangle}\)

Wenn sowohl Zustand als auch Operator als Tensor-Netzwerke vorliegen, wird dieser Ausdruck selbst zu einem Tensor-Netzwerk. Dabei werden Bra-PEPS, PEPO und Ket-PEPS übereinandergelegt.

Die physikalischen Indizes werden lokal kontrahiert. Übrig bleibt ein Netzwerk aus virtuellen Indizes, dessen Kontraktion den gesuchten Erwartungswert ergibt.

Korrelationsfunktionen

Zweipunktkorrelationen besitzen die Form

\(C_{AB}(r,s) = \langle A_r B_s\rangle - \langle A_r\rangle\langle B_s\rangle\)

Sie geben Auskunft darüber, wie stark zwei räumlich getrennte Bereiche eines Quantensystems miteinander korreliert sind.

In Tensor-Netzwerken können solche Observablen lokal eingesetzt werden. Die eigentliche Herausforderung besteht wiederum darin, den Einfluss des umgebenden unendlichen oder sehr großen Gitters ausreichend genau zu approximieren.

Kontraktion von PEPO-Netzwerken

Warum die Kontraktion schwierig ist

Bei einem eindimensionalen MPS oder MPO können Tensoren schrittweise entlang der Kette kontrahiert werden. Die Netzwerkgeometrie ist vergleichsweise einfach.

Ein zweidimensionales PEPO enthält dagegen geschlossene Schleifen. Wird eine Zeile oder Spalte exakt in die nächste kontrahiert, wächst die effektive Dimension der entstehenden Randstruktur fortlaufend.

Eine exakte Kontraktion eines allgemeinen großflächigen zweidimensionalen Tensor-Netzwerks ist deshalb im Regelfall nicht effizient möglich. In der Praxis werden approximative Verfahren eingesetzt.

Boundary-MPS

Eine verbreitete Strategie interpretiert den Rand eines zweidimensionalen Netzwerks als eindimensionales Tensor-Netzwerk. Eine Zeile des PEPO wird dann schrittweise auf diesen Randzustand angewendet.

Nach jeder Anwendung würde die Bond-Dimension des Rand-MPS anwachsen. Daher wird es anschließend komprimiert.

Das Verfahren reduziert die zweidimensionale Kontraktion auf eine Folge eindimensionaler Operationen. Seine Genauigkeit hängt maßgeblich von der für den Rand zugelassenen Bond-Dimension ab.

Corner Transfer Matrix

Ein weiteres wichtiges Verfahren basiert auf Corner-Transfer-Matrizen. Die Umgebung eines lokalen Bereichs wird dabei durch effektive Eck- und Randtensoren beschrieben.

Statt das gesamte unendliche Netzwerk explizit zu kontrahieren, konstruiert man iterativ eine effektive Umgebung. Nach Erweiterung der Umgebung werden die wachsenden Tensoren wieder trunciert.

Dieses Vorgehen ist besonders für translationsinvariante PEPS- und PEPO-Berechnungen auf unendlichen Gittern wichtig.

Tensor-Renormierungsverfahren

Tensor-Renormierungsverfahren verfolgen eine andere Perspektive. Das Netzwerk wird schrittweise grobkörniger gemacht. Gruppen benachbarter Tensoren werden zusammengefasst und anschließend durch Tensoren geringerer effektiver Dimension approximiert.

Nach mehreren Schritten repräsentiert ein einzelner effektiver Tensor einen immer größeren ursprünglichen Gitterbereich.

Die Herausforderung besteht darin, bei der Dimensionsreduktion genau die Informationen zu bewahren, die für die großskalige Physik relevant sind.

Umgebungsgenauigkeit

Bei jeder Trunkierung stellt sich die Frage, welche Komponenten verworfen werden dürfen. Eine rein lokale Entscheidung ist schnell, berücksichtigt aber nicht vollständig, wie wichtig ein virtueller Freiheitsgrad für den Rest des Netzwerks ist.

Genauere Verfahren beziehen deshalb die effektive Umgebung des zu optimierenden Tensors ein. Dies erhöht die Rechenkosten, kann aber insbesondere bei stark korrelierten Systemen deutlich bessere Ergebnisse liefern.

Konvergenzkontrolle

Eine einzelne Rechnung bei einer festen Bond-Dimension liefert noch keinen zuverlässigen Nachweis numerischer Genauigkeit. Typischerweise werden Rechnungen mit steigenden Bond-Dimensionen wiederholt.

Konvergieren relevante Observablen gegen stabile Werte, steigt das Vertrauen in das Ergebnis. Besondere Vorsicht ist in der Nähe kritischer Punkte erforderlich, weil dort lange Korrelationslängen größere virtuelle Dimensionen verlangen können.

Numerische Optimierung und Approximation

Trunkierung

Die Trunkierung ist für praktische PEPO-Algorithmen unverzichtbar. Nach Operatorprodukten, Zeitentwicklung oder thermischer Evolution enthält ein Netzwerk häufig mehr virtuelle Freiheitsgrade, als dauerhaft gespeichert werden können.

Der Kern der Trunkierung besteht darin, die virtuelle Dimension von einem größeren Wert auf eine Zielgröße zu reduzieren.

Schematisch kann dies als Abbildung

\(\chi_{\mathrm{groß}} \rightarrow \chi_{\mathrm{ziel}}\)

verstanden werden.

Die entscheidende Frage lautet nicht, ob Informationen verloren gehen, sondern welche Informationen mit möglichst geringem Einfluss auf die gesuchten Observablen verworfen werden.

Lokale und globale Optimierung

Lokale Verfahren optimieren einzelne Tensoren oder Bindungen bei festgehaltener Umgebung. Sie sind vergleichsweise schnell, können aber die globale Struktur nur näherungsweise berücksichtigen.

Aufwendigere Verfahren beziehen eine effektive Umgebung ein und minimieren einen globaler definierten Fehler. Ein mögliches Ziel ist beispielsweise

\(\min_{\widetilde{O}} ||O-\widetilde{O}||^2\)

unter der Bedingung, dass \(\widetilde{O}\) eine vorgegebene maximale Bond-Dimension besitzt.

Fehlerquellen

PEPO-Rechnungen können mehrere voneinander unabhängige Fehlerquellen enthalten. Dazu gehören eine endliche Bond-Dimension, approximative Netzwerkkontraktion, Trotter-Fehler, endliche Systemgrößen und numerische Rundungsfehler.

Für belastbare Ergebnisse müssen diese Effekte soweit möglich getrennt untersucht werden.

Speicher- und Rechenaufwand

Schon ein einzelner lokaler PEPO-Tensor kann bei vier virtuellen Indizes eine Anzahl von Elementen besitzen, die näherungsweise wie eine hohe Potenz der Bond-Dimension wächst. Bei gleichen virtuellen Dimensionen kann die reine Tensorgröße schematisch proportional zu

\(d^2 \chi^4\)

sein.

Kontraktionen mehrerer solcher Tensoren können noch deutlich höhere Potenzen von \(\chi\) erzeugen. Deshalb ist eine scheinbar moderate Vergrößerung der Bond-Dimension numerisch keineswegs harmlos.

PEPOs in der Simulation von Quantensystemen

Zweidimensionale Spinmodelle

Zweidimensionale Spinmodelle gehören zu den naheliegendsten Einsatzbereichen von PEPOs. Ein Beispiel ist das transversale Ising-Modell mit einem Hamiltonoperator der Form

\(H = -J\sum_{\langle r,s\rangle}\sigma_r^z\sigma_s^z - h\sum_r\sigma_r^x\)

Der Hamiltonoperator besteht aus lokalen und nächsten-nachbarschaftlichen Termen und eignet sich deshalb grundsätzlich für eine kompakte Tensor-Netzwerk-Repräsentation.

Von Interesse sind Grundzustände, thermische Zustände, Korrelationsfunktionen und das Verhalten in der Umgebung von Phasenübergängen.

Heisenberg-artige Modelle

Ein weiterer grundlegender Fall ist der Heisenberg-Hamiltonoperator

\(H = J\sum_{\langle r,s\rangle}\left(S_r^xS_s^x + S_r^yS_s^y + S_r^zS_s^z\right)\)

Solche Modelle beschreiben wechselwirkende magnetische Freiheitsgrade und spielen eine zentrale Rolle bei der Untersuchung von Quantenmagnetismus.

Der Operator selbst besitzt eine lokale Struktur. Die Schwierigkeit liegt daher weniger in der formalen PEPO-Darstellung des Hamiltonoperators als in der genauen Simulation des daraus resultierenden stark korrelierten Quantenzustands.

Fermionische Systeme

Fermionische Modelle stellen zusätzliche Anforderungen, weil Fermionoperatoren antikommutieren. Für Erzeugungs- und Vernichtungsoperatoren gelten Beziehungen wie

\(\{c_i,c_j^\dagger\} = \delta_{ij}\)

und

\(\{c_i,c_j\}=0\)

Die Vorzeichenstruktur muss in einem fermionischen Tensor-Netzwerk korrekt berücksichtigt werden. Dafür existieren spezielle fermionische PEPS- und PEPO-Formulierungen.

Solche Methoden sind für Modelle wechselwirkender Elektronen und andere fermionische Gitterprobleme relevant.

Reale Zeitentwicklung

Neben imaginärer Zeit kann auch die reale Quantendynamik betrachtet werden. Der Zeitentwicklungsoperator lautet

\(U(t) = e^{-iHt}\)

Ein anfänglicher Zustand entwickelt sich gemäß

\(|\psi(t)\rangle = e^{-iHt}|\psi(0)\rangle\)

Für kurze Zeiten kann der Evolutionsoperator häufig durch ein Tensor-Netzwerk mit begrenzter Bond-Dimension approximiert werden.

Das grundlegende Problem besteht darin, dass sich Korrelationen und Verschränkung während der Dynamik ausbreiten können. Dadurch steigt die benötigte Bond-Dimension mit der Zeit häufig deutlich an.

Quenches

Bei einem Quantum Quench wird ein Parameter des Hamiltonoperators plötzlich verändert. Das System befindet sich unmittelbar danach nicht im Gleichgewicht des neuen Hamiltonoperators und beginnt sich zeitlich zu entwickeln.

Tensor-Netzwerke ermöglichen die Untersuchung dieser Nichtgleichgewichtsdynamik, solange die entstehende Komplexität numerisch kontrollierbar bleibt.

PEPOs können dabei unter anderem zur Darstellung von Evolutionsoperatoren oder von Observablen in der Heisenberg-Darstellung eingesetzt werden.

PEPOs und offene Quantensysteme

Wechselwirkung mit der Umgebung

Reale Quantensysteme sind selten vollständig isoliert. Wechselwirkungen mit einer Umgebung führen zu Dekohärenz, Dissipation und Rauschen.

Der Zustand eines offenen Systems wird durch eine Dichtematrix beschrieben. Damit liegt bereits eine natürliche Verbindung zu Operator-Tensor-Netzwerken vor.

Lindblad-Gleichung

Eine verbreitete Beschreibung Markovscher offener Quantensysteme ist die Lindblad-Gleichung:

\(\frac{d\rho}{dt} = -i[H,\rho] + \sum_\mu \left(L_\mu\rho L_\mu^\dagger - \frac{1}{2}L_\mu^\dagger L_\mu\rho - \frac{1}{2}\rho L_\mu^\dagger L_\mu\right)\)

Hier beschreibt H die kohärente Dynamik und \(L_\mu\) die dissipativen Kanäle.

Wird die Dichtematrix als PEPO dargestellt, kann ihre zeitliche Entwicklung prinzipiell direkt innerhalb eines Operator-Tensor-Netzwerks behandelt werden.

Stationäre Zustände

Von besonderem Interesse sind stationäre Zustände \(\rho_{\mathrm{ss}}\), für die

\(\frac{d\rho_{\mathrm{ss}}}{dt}=0\)

gilt.

Solche Zustände können Nichtgleichgewichtsphänomene darstellen, die in geschlossenen Systemen kein direktes Gegenstück besitzen. PEPO-Methoden bieten eine Möglichkeit, entsprechende Dichteoperatoren in großen zweidimensionalen Gittersystemen näherungsweise zu behandeln.

Physikalische Konsistenz

Bei offenen Systemen ist die physikalische Konsistenz der Dichtematrix besonders wichtig. Näherungsfehler dürfen nicht unkontrolliert zu negativen Eigenwerten oder einer falschen Normierung führen.

Die Bedingung

\(\mathrm{Tr}(\rho)=1\)

kann leicht überprüft und korrigiert werden. Positivität ist dagegen wesentlich anspruchsvoller. Deshalb werden bei bestimmten Algorithmen Darstellungen bevorzugt, in denen Positivität konstruktiv eingebaut ist.

PEPO und verwandte Tensor-Netzwerk-Strukturen

PEPO und MPO

MPO und PEPO beruhen auf derselben Grundidee. Der Unterschied liegt vor allem in der Netzwerkgeometrie.

Ein MPO beschreibt typischerweise eine eindimensionale Kette. Ein PEPO erweitert das Konzept auf zweidimensionale oder allgemein höherdimensionale Gitter.

Diese Erweiterung verändert die numerische Situation grundlegend. Die zusätzliche räumliche Dimension führt zu Schleifen im Netzwerk und macht die Kontraktion deutlich schwieriger.

PEPO und PEPS

PEPS und PEPO besitzen nahezu dieselbe virtuelle Netzwerkgeometrie. Der wesentliche Unterschied liegt in den physikalischen Indizes.

Ein PEPS-Tensor besitzt pro lokalem Freiheitsgrad einen physikalischen Index. Ein PEPO-Tensor besitzt typischerweise zwei: einen Eingangs- und einen Ausgangsindex.

Dadurch kann ein PEPO auf einen PEPS angewendet werden, ohne die räumliche Grundstruktur des Gitters zu verändern.

Matrix Product Density Operators

Matrix Product Density Operators sind MPO-Darstellungen von Dichtematrizen in einer Dimension. Ein PEPO kann als entsprechende Verallgemeinerung auf höhere Dimensionen betrachtet werden.

Die physikalischen Anforderungen an einen Dichteoperator bleiben unabhängig von der Netzwerkgeometrie gleich. Hermitizität, Positivität und Normierung müssen weiterhin berücksichtigt werden.

Tree Tensor Networks

Tree Tensor Networks besitzen eine baumartige Struktur ohne geschlossene Schleifen. Dadurch lassen sich viele Kontraktionen effizienter durchführen.

Die Geometrie entspricht jedoch nicht direkt einem zweidimensionalen lokalen Gitter. PEPOs bilden die lokale Nachbarschaft eines solchen Gitters natürlicher ab, bezahlen diese Nähe zur physikalischen Geometrie aber mit schwierigeren Kontraktionen.

MERA

Die Multiscale Entanglement Renormalization Ansatz, MERA, organisiert Tensoren hierarchisch über verschiedene Längenskalen. Sie wurde insbesondere entwickelt, um Verschränkungsstrukturen kritischer Systeme effizient abzubilden.

PEPOs verfolgen dagegen in ihrer Standardform eine lokale Gitterstruktur. Sie sind keine explizite Hierarchie verschiedener Skalen, sondern ordnen lokale Operator-Tensoren entsprechend der räumlichen Struktur des betrachteten Systems an.

Bedeutung für die Quantentechnologie

Klassische Simulation von Quantenmaterie

Ein wesentlicher Nutzen von PEPOs liegt in der klassischen Simulation komplexer Quantensysteme. Moderne Quantentechnologie benötigt nicht nur Quantenhardware, sondern ebenso verlässliche klassische Referenzmethoden.

Tensor-Netzwerke können für bestimmte Klassen von Vielteilchenproblemen Ergebnisse liefern, die zur Entwicklung, Prüfung und Interpretation von Quantensimulatoren beitragen.

Thermische Quantensimulation

Viele reale Anwendungen interessieren sich nicht ausschließlich für den Grundzustand bei absoluter Nulltemperatur. Materialien, Quantensimulatoren und experimentelle Plattformen arbeiten bei endlichen Temperaturen.

PEPOs bieten einen direkten Rahmen für Dichteoperatoren und thermische Operatoren. Dadurch sind sie für die Untersuchung temperaturabhängiger Quanteneffekte besonders relevant.

Verifikation von Quantensimulatoren

Ein Quantensimulator soll ein komplexes Vielteilchensystem kontrolliert nachbilden. Solange sich bestimmte Parameterbereiche noch klassisch mit Tensor-Netzwerken berechnen lassen, können PEPO- und PEPS-Ergebnisse als Vergleich dienen.

Weichen experimentelle und numerische Daten voneinander ab, kann untersucht werden, ob die Ursache in experimentellem Rauschen, Modellannahmen oder einer unzureichenden klassischen Approximation liegt.

Quantenkanäle

Nicht nur Zustände und Hamiltonoperatoren, sondern auch Quantenkanäle können Tensor-Netzwerk-Strukturen besitzen. Ein allgemeiner Quantenkanal E wirkt auf eine Dichtematrix gemäß

\(\rho' = E(\rho)\)

Lokale oder räumlich strukturierte Kanäle können mithilfe von Operator-Tensor-Netzwerken beschrieben werden. Dies ist für Rauschmodelle, dissipative Dynamik und die Untersuchung räumlich verteilter Quantenprozesse von Bedeutung.

Verbindung zu Quantencomputern

PEPOs ersetzen keinen universellen Quantencomputer. Ihre Stärke liegt vielmehr darin, bestimmte strukturierte Quantensysteme auf klassischen Rechnern näherungsweise handhabbar zu machen.

Gerade diese Grenze ist technologisch interessant. Solange eine Quantenschaltung, ein Quantenzustand oder eine Observable eine ausreichend begrenzte Tensor-Netzwerk-Komplexität besitzt, kann eine klassische Simulation möglich bleiben. Mit wachsender Verschränkung oder Operator-Komplexität steigt hingegen der klassische Aufwand.

Tensor-Netzwerke sind deshalb auch ein Werkzeug, um zu verstehen, welche Quantenprozesse tatsächlich schwer klassisch simulierbar sind.

Vorteile und Grenzen von PEPOs

Kompakte Darstellung strukturierter Operatoren

Der wichtigste Vorteil eines PEPO besteht darin, dass ein globaler Vielteilchenoperator nicht als vollständige exponentiell große Matrix gespeichert werden muss.

Bei lokaler Struktur und moderater Bond-Dimension kann die Zahl der gespeicherten Parameter drastisch geringer sein.

Dieser Vorteil ist jedoch strukturell und nicht universell. Ein völlig allgemeiner Operator auf einem großen Hilbertraum lässt sich nicht allein dadurch effizient darstellen, dass er PEPO genannt wird. Die Effizienz entsteht nur, wenn eine ausreichend kompakte Tensor-Netzwerk-Struktur existiert.

Natürliche Beschreibung zweidimensionaler Geometrien

PEPOs bilden die lokale Nachbarschaft zweidimensionaler Gitter direkt ab. Ein Tensor kann unmittelbar mit Tensoren der räumlich angrenzenden Gitterplätze verbunden werden.

Diese geometrische Natürlichkeit ist ein wesentlicher Unterschied zu eindimensionalen Darstellungen, in denen ein zweidimensionales Gitter zunächst in eine lange Kette umgeordnet werden müsste. Eine solche Umordnung kann ursprünglich lokale Wechselwirkungen effektiv langreichweitig erscheinen lassen.

Flexible Operatorrepräsentation

PEPOs können Hamiltonoperatoren, Dichtematrizen, thermische Operatoren, Evolutionsoperatoren und andere strukturierte Vielteilchenoperatoren darstellen.

Damit bilden sie keine Methode für nur eine einzige physikalische Aufgabe, sondern eine allgemeine Repräsentationsklasse.

Schwierige Kontraktion

Der größte Nachteil ist die Kontraktion. Selbst wenn ein PEPO kompakt gespeichert werden kann, bedeutet dies nicht automatisch, dass alle gewünschten Größen effizient berechnet werden können.

Ein Tensor-Netzwerk kann lokal aus relativ kleinen Tensoren bestehen und trotzdem eine global schwierige Kontraktionsaufgabe darstellen.

Dieser Unterschied zwischen kompakter Repräsentation und effizienter Auswertung ist für das Verständnis von PEPOs zentral.

Wachstum der Bond-Dimension

Produkte von Operatoren, Zeitentwicklung und thermische Evolution vergrößern typischerweise die Bond-Dimension.

Wird nicht komprimiert, wächst der Rechenaufwand schnell über jede praktische Grenze. Wird stark komprimiert, entstehen Approximationsfehler.

Die Qualität eines PEPO-Verfahrens hängt daher wesentlich davon ab, wie gut die relevante Information bei begrenzter Bond-Dimension erhalten werden kann.

Kritische Systeme

In der Nähe eines kontinuierlichen Phasenübergangs können Korrelationslängen sehr groß werden. Eine endliche Bond-Dimension begrenzt jedoch die Komplexität der darstellbaren Korrelationsstruktur.

Deshalb erfordern kritische Systeme oft deutlich größere Bond-Dimensionen und sorgfältige Skalierungsanalysen.

Ein scheinbar konvergiertes Ergebnis bei kleiner Bond-Dimension kann sonst Eigenschaften des numerischen Ansatzes widerspiegeln statt die tatsächliche kritische Physik.

Langzeitdynamik

Auch die reale Zeitentwicklung stößt an fundamentale praktische Grenzen. Nach einem Quench können sich Verschränkung und Operator-Komplexität schnell ausbreiten.

Die für eine genaue Darstellung notwendige Bond-Dimension steigt dann mit der Zeit. PEPO-basierte Zeitentwicklung ist deshalb häufig für kurze oder mittlere Zeiten besonders leistungsfähig, während sehr lange Zeiten wesentlich schwieriger werden.

Aktuelle und weiterführende Forschungsrichtungen

Verbesserte Kontraktionsalgorithmen

Ein großer Teil der Weiterentwicklung von PEPS- und PEPO-Methoden konzentriert sich auf effizientere Kontraktionen. Jede Verbesserung der Umgebungsapproximation kann entweder die Genauigkeit bei gleichem Rechenaufwand steigern oder größere Bond-Dimensionen zugänglich machen.

Von besonderem Interesse sind Algorithmen, die unnötige Zwischen-Dimensionen vermeiden und Tensor-Kontraktionsreihenfolgen automatisch optimieren.

Symmetrien

Viele Quantenmodelle besitzen Erhaltungsgrößen oder Symmetrien. Wird eine solche Struktur explizit in den Tensoren kodiert, entstehen Blockstrukturen.

Statt einen vollständig dichten Tensor zu speichern, müssen dann nur die symmetry-kompatiblen Blöcke berücksichtigt werden.

Für eine Erhaltungsgröße Q kann beispielsweise eine Bedingung der Form

\([H,Q]=0\)

gelten.

Die Nutzung entsprechender Symmetrien kann den Speicherbedarf reduzieren, Rechenoperationen einsparen und sicherstellen, dass die numerische Entwicklung im gewünschten Symmetriesektor bleibt.

Automatische Differenzierung

Tensor-Netzwerke können auch als parametrisierte Rechenmodelle betrachtet werden. Eine Zielfunktion hängt von vielen lokalen Tensorparametern ab. Automatische Differenzierung ermöglicht es, Gradienten dieser Zielfunktion systematisch zu bestimmen.

Für eine Zielfunktion L und Tensorparameter \(\theta\) wird beispielsweise

\(\frac{\partial L}{\partial\theta}\)

benötigt.

Dies eröffnet Möglichkeiten für gradientenbasierte Optimierung, variationale Tensor-Netzwerk-Algorithmen und Verbindungen zu Methoden des maschinellen Lernens.

Hochleistungsrechnen

Viele Tensoroperationen bestehen aus großen linearen Algebraoperationen und lassen sich gut auf moderne Hardware abbilden. GPUs und verteilte Rechensysteme können daher den praktischen Anwendungsbereich von PEPO-Verfahren erweitern.

Allerdings beseitigt zusätzliche Rechenleistung nicht das grundlegende Skalierungsproblem. Wenn der Aufwand mit einer hohen Potenz der Bond-Dimension wächst, kann auch leistungsfähige Hardware nur einen begrenzten Bereich zusätzlich erschließen.

Lokal genaue Operatorapproximationen

Eine wichtige konzeptionelle Richtung besteht darin, nicht zwingend den vollständigen globalen Operator mit hoher Genauigkeit zu approximieren. Für viele physikalische Fragestellungen genügt es, lokale Observablen korrekt wiederzugeben.

Gesucht wird dann ein PEPO \(\widetilde{O}\), für das relevante lokale Erwartungswerte des ursprünglichen Operators O mit kontrolliertem Fehler reproduziert werden.

Diese Perspektive ist besonders bedeutsam, weil Quantentechnologie in der Praxis häufig lokale Messgrößen, Korrelationsfunktionen und begrenzte räumliche Ausschnitte untersucht, nicht sämtliche Einträge einer exponentiell großen Operator-Matrix.

Tomographie und Rekonstruktion

Eine weitere Forschungsrichtung betrifft die Rekonstruktion strukturierter Quantenzustände und Dichteoperatoren aus Messdaten. Wenn bekannt oder plausibel ist, dass ein Zustand eine PEPS- oder PEPO-Struktur mit begrenzter Bond-Dimension besitzt, kann diese Struktur als Vorwissen genutzt werden.

Das Ziel besteht darin, nicht einen völlig allgemeinen exponentiell großen Dichteoperator zu rekonstruieren, sondern direkt innerhalb einer strukturierten Tensor-Netzwerk-Klasse zu arbeiten.

Damit berührt die PEPO-Theorie auch Fragen der Quantentomographie, der Datenkompression und der Charakterisierung großer Quantensysteme.

Physikalische Interpretation eines PEPO

Mehr als eine Datenkompression

Es wäre zu kurz gegriffen, einen PEPO lediglich als mathematische Kompression einer großen Matrix zu betrachten. Die Netzwerkstruktur spiegelt häufig die lokale Organisation der zugrunde liegenden Physik wider.

Die physikalischen Indizes repräsentieren lokale Ein- und Ausgänge. Die virtuellen Indizes vermitteln die Abhängigkeiten zwischen verschiedenen Regionen. Die Bond-Dimension beschreibt, wie umfangreich diese vermittelten Informationen sein können.

Damit verbindet die PEPO-Darstellung lokale Physik und globale Korrelation in einer einzigen mathematischen Struktur.

Virtuelle Freiheitsgrade

Die virtuellen Indizes sind keine zusätzlichen realen Teilchen. Sie sind Hilfsvariablen der mathematischen Darstellung.

Trotzdem tragen sie physikalisch relevante Information. Wenn eine Bond-Dimension erhöht werden muss, um ein Problem genauer zu beschreiben, zeigt dies, dass eine komplexere nichtlokale Struktur zwischen verschiedenen Gitterregionen berücksichtigt werden muss.

Operatorverschränkung

Ein Operator kann mathematisch selbst als Vektor in einem vergrößerten Raum aufgefasst werden. Dazu werden Eingangs- und Ausgangsindizes geeignet zusammengefasst.

Für einen Operator

\(O = \sum_{ij} O_{ij}|i\rangle\langle j|\)

kann formal ein zugehöriger Vektor

\(|O\rangle = \sum_{ij} O_{ij}|i\rangle\otimes|j\rangle\)

betrachtet werden.

Diese Sichtweise ermöglicht es, Konzepte der Zustandsverschränkung auf Operatoren zu übertragen. Eine geringe Operatorverschränkung begünstigt kompakte Tensor-Netzwerk-Darstellungen. Stark wachsende Operatorverschränkung erklärt dagegen, warum bestimmte Zeitentwicklungen zunehmend schwer klassisch simulierbar werden.

Lokalität als Ressource

Lokalität ist eine der entscheidenden Ressourcen hinter der praktischen Nutzbarkeit von PEPOs. Ein lokaler Hamiltonoperator enthält nicht beliebige Kopplungen zwischen sämtlichen Freiheitsgraden. Informationen breiten sich über die Wechselwirkungsstruktur des Systems aus.

Diese Einschränkung erlaubt es, globale Quantenprobleme über lokale Tensoren aufzubauen. Die Methode ist daher besonders leistungsfähig, wenn Gittergeometrie, Wechselwirkungsreichweite und Korrelationen zueinander passen.

Praktischer Ablauf einer PEPO-Simulation

Definition des Modells

Am Anfang steht die physikalische Modellierung. Festgelegt werden Gittergeometrie, lokale Hilbertraumdimension, Hamiltonoperator, Randbedingungen und gegebenenfalls Temperatur oder Dynamik.

Bereits an dieser Stelle sollte geprüft werden, welche Symmetrien und lokalen Strukturen genutzt werden können.

Konstruktion des Operators

Der relevante Operator wird in eine PEPO-Struktur übersetzt. Bei einem Hamiltonoperator erfolgt dies typischerweise aus lokalen Wechselwirkungstermen. Bei einem thermischen Zustand wird häufig zunächst ein einfacher Operator initialisiert und anschließend imaginärzeitentwickelt.

Festlegung der Bond-Dimension

Eine kleine Bond-Dimension reduziert Rechenzeit und Speicherbedarf, begrenzt aber die darstellbare Struktur. Deshalb ist die Bond-Dimension ein Konvergenzparameter und sollte nicht als beliebige feste Konstante behandelt werden.

Wahl des Kontraktionsverfahrens

Anschließend muss entschieden werden, wie die Umgebung des Netzwerks approximiert wird. Boundary-MPS, Corner-Transfer-Verfahren und verschiedene Renormierungsansätze besitzen unterschiedliche Stärken.

Die geeignete Wahl hängt von Gittergeometrie, Translationsinvarianz, Bond-Dimension und gewünschter Observable ab.

Iteration und Kompression

Bei Zeitentwicklung oder thermischen Rechnungen wechseln sich Operatoranwendung und Kompression ab. Nach jedem Schritt muss verhindert werden, dass die Bond-Dimension unkontrolliert wächst.

Diese Kompression ist keine nebensächliche technische Operation. Sie bestimmt wesentlich die Qualität der gesamten Simulation.

Messung von Observablen

Nach ausreichender Konvergenz werden physikalische Größen bestimmt. Dazu gehören lokale Magnetisierungen, Energien, Korrelationsfunktionen, Dichteprofile und thermodynamische Größen.

Für jede Observable entsteht eine Tensor-Netzwerk-Kontraktion, die mit der zuvor konstruierten Umgebung ausgewertet wird.

Konvergenzprüfung

Die Berechnung sollte anschließend für verschiedene Bond-Dimensionen, Umgebungsdimensionen und gegebenenfalls Trotter-Schrittweiten wiederholt werden.

Erst wenn relevante Observablen unter diesen Änderungen stabil bleiben, kann von numerischer Konvergenz gesprochen werden.

Einordnung gegenüber vollständiger exakter Berechnung

Exakte Diagonalisierung

Bei kleinen Systemen kann ein Hamiltonoperator vollständig als Matrix aufgebaut und diagonalisiert werden. Dies liefert sehr genaue Referenzergebnisse.

Die Methode scheitert jedoch schnell am exponentiellen Wachstum des Hilbertraums.

PEPOs verfolgen einen grundsätzlich anderen Ansatz. Sie verzichten auf die vollständige Speicherung aller Freiheitsgrade und nutzen stattdessen die Struktur des Operators.

Was PEPOs nicht leisten

Ein PEPO macht ein beliebiges Quantenproblem nicht automatisch einfach. Es gibt keine Garantie, dass jeder physikalisch interessante Operator mit kleiner Bond-Dimension genau dargestellt werden kann.

Ebenso wenig ist garantiert, dass eine vorhandene kompakte PEPO-Darstellung effizient exakt kontrahiert werden kann.

Diese beiden Punkte müssen klar getrennt werden: Darstellungskomplexität und Kontraktionskomplexität sind unterschiedliche Probleme.

Wo PEPOs besonders sinnvoll sind

Besonders geeignet sind PEPOs, wenn mehrere günstige Eigenschaften zusammentreffen: lokale Wechselwirkungen, begrenzte Operatorverschränkung, zweidimensionale Gitterstruktur, relevante lokale Observablen und eine numerisch beherrschbare Bond-Dimension.

Unter diesen Bedingungen können sie Größenordnungen von Systemen zugänglich machen, die einer direkten Matrixdarstellung vollständig entzogen sind.

Fazit

Der Projected Entangled Pair Operator ist die natürliche Operatorerweiterung der Projected Entangled Pair States. Während ein PEPS einen quantenmechanischen Vielteilchenzustand über ein Netzwerk lokaler Tensoren beschreibt, repräsentiert ein PEPO einen Vielteilchenoperator durch lokale Tensoren mit physikalischen Ein- und Ausgangsindizes sowie virtuellen Verbindungen.

Der entscheidende Vorteil liegt in der strukturierten Darstellung. Ein allgemeiner Operator eines Systems aus N lokalen Freiheitsgraden besitzt eine exponentiell große Zahl möglicher Matrixelemente. Lokale Hamiltonoperatoren, thermische Zustände und viele andere physikalisch relevante Operatoren besitzen jedoch eine wesentlich stärkere innere Struktur. PEPOs versuchen, genau diese Struktur auszunutzen.

Die Bond-Dimension ist dabei der zentrale Kontrollparameter. Sie begrenzt die Informationsmenge, die über virtuelle Verbindungen zwischen verschiedenen Regionen des Systems transportiert werden kann. Kleine Bond-Dimensionen ermöglichen effiziente Rechnungen, können aber komplexe Korrelationen unzureichend darstellen. Große Bond-Dimensionen steigern die Genauigkeit, lassen Speicherbedarf und Rechenaufwand jedoch stark anwachsen.

Die größte Herausforderung ist nicht allein die Darstellung, sondern die Kontraktion des zweidimensionalen Tensor-Netzwerks. Anders als eindimensionale Matrix-Product-Strukturen enthalten PEPOs Schleifen. Exakte Kontraktionen werden deshalb für große Systeme schnell unpraktikabel. Boundary-MPS, Corner-Transfer-Verfahren und Tensor-Renormierung ersetzen die exakte Umgebung durch kontrollierte Approximationen.

Besonders wichtig sind PEPOs für Dichteoperatoren und thermische Quantenzustände. Der Gibbs-Operator

\(e^{-\beta H}\)

ist selbst ein Vielteilchenoperator und kann über imaginäre Zeitentwicklung schrittweise aufgebaut werden. Dadurch lassen sich zweidimensionale Quantensysteme bei endlicher Temperatur untersuchen. Weitere Anwendungen reichen von Hamiltonoperatoren und realer Zeitentwicklung über offene Quantensysteme bis zu Quantenkanälen und lokalen Operatorapproximationen.

Für die Quantentechnologie sind PEPOs deshalb vor allem als klassische Struktur- und Simulationswerkzeuge relevant. Sie ermöglichen Referenzrechnungen für Quantensimulatoren, Untersuchungen stark korrelierter Quantenmaterie und Analysen darüber, wann ein Quantenprozess noch klassisch durch Tensor-Netzwerke beschrieben werden kann.

PEPOs umgehen die exponentielle Komplexität der Quantenmechanik nicht grundsätzlich. Ihre Stärke besteht vielmehr darin, jene Vielteilchenoperatoren zu identifizieren und effizient darzustellen, deren physikalische Struktur erheblich kleiner ist als der vollständige mathematische Operatorraum. Genau darin liegt ihre Bedeutung: Sie ersetzen die vollständige Beschreibung eines exponentiell großen Systems durch eine lokal organisierte Darstellung seiner tatsächlich relevanten Quantenstruktur.

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

Anhang

Wissenschaftliche Zeitschriften und Artikel

Die wissenschaftliche Literatur zu Projected Entangled Pair Operators (PEPOs) ist eng mit der Entwicklung von Projected Entangled Pair States (PEPS), Matrix Product Operators (MPOs), Tensor-Netzwerk-Kontraktionen und Verfahren für thermische sowie offene Quantensysteme verbunden. Für eine belastbare wissenschaftliche Abhandlung sollten deshalb sowohl die grundlegenden Arbeiten zur PEPS-Struktur als auch Publikationen berücksichtigt werden, in denen PEPOs ausdrücklich als Operatoransatz für thermische Zustände, Zeitentwicklung oder dissipative Vielteilchensysteme eingesetzt werden.

Grundlegende Primärliteratur zu PEPS und PEPO

  • F. Verstraete und J. I. Cirac: Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions, arXiv Preprint, 2004.
    • Diese Arbeit gehört zu den grundlegenden Publikationen zur Einführung von Projected Entangled Pair States als Verallgemeinerung eindimensionaler Matrix Product States auf zwei und mehr Dimensionen. Sie ist für eine PEPO-Abhandlung besonders wichtig, weil die Netzwerkgeometrie und die virtuellen Bindungen eines PEPO unmittelbar auf der PEPS-Struktur aufbauen. Die Arbeit eignet sich als Primärquelle für die historische und konzeptionelle Entstehung der höherdimensionalen Tensor-Netzwerk-Darstellung.
  • F. Verstraete, M. M. Wolf, D. Pérez-García und J. I. Cirac: Criticality, the Area Law, and the Computational Power of Projected Entangled Pair States, Physical Review Letters 96, 220601, 2006.
    • Die Veröffentlichung untersucht grundlegende Eigenschaften von PEPS, darunter Flächengesetze, kritische Zustände und die rechnerische Komplexität. Für PEPOs ist sie vor allem als theoretischer Hintergrund relevant: Die Leistungsfähigkeit und die Grenzen der Operatornetzwerke hängen eng mit den Verschränkungs- und Korrelationsstrukturen der zugrunde liegenden PEPS-Geometrie zusammen. Die Quelle eignet sich für Abschnitte über Bond-Dimension, Verschränkung und die prinzipiellen Grenzen klassischer Tensor-Netzwerk-Simulationen.
  • F. Verstraete, V. Murg und J. I. Cirac: Matrix product states, projected entangled pair states, and variational renormalization group methods for quantum spin systems, Advances in Physics 57, 143–224, 2008.
    • Diese umfangreiche Übersichtsarbeit verbindet Matrix Product States, PEPS, Renormierungsverfahren und die Simulation stark korrelierter Quantensysteme. Sie ist eine wichtige Grundlagenquelle, um den Übergang von eindimensionalen MPS/MPO-Strukturen zu höherdimensionalen PEPS/PEPO-Strukturen systematisch herzuleiten. Besonders geeignet ist sie für die Erklärung virtueller Bindungen, Bond-Dimensionen, Tensor-Kontraktionen und variationaler Verfahren.
  • J. Jordan, R. Orús, G. Vidal, F. Verstraete und J. I. Cirac: Classical Simulation of Infinite-Size Quantum Lattice Systems in Two Spatial Dimensions, Physical Review Letters 101, 250602, 2008.
    • Die Arbeit entwickelt einen Tensor-Netzwerk-Ansatz für translationsinvariante zweidimensionale Quantensysteme direkt im thermodynamischen Limes. Sie ist für PEPOs methodisch relevant, weil unendliche PEPO- beziehungsweise iPEPO-Verfahren auf derselben Grundidee wiederholter lokaler Tensoren und approximierter Umgebungen beruhen. Die Publikation kann zur Einordnung unendlicher Tensor-Netzwerke und zweidimensionaler Zeitentwicklungsverfahren genutzt werden.

Übersichtsarbeiten zu Tensor-Netzwerken und PEPS

  • Román Orús: A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States, Annals of Physics 349, 117–158, 2014.
    • Diese didaktisch ausgerichtete Übersichtsarbeit ist besonders geeignet, um Tensor-Netzwerk-Notation, Kontraktionen, MPS und PEPS verständlich einzuführen. Für eine PEPO-Abhandlung bietet sie die mathematischen und numerischen Grundlagen, die notwendig sind, bevor der zusätzliche physikalische Ausgangsindex eines Operatornetzwerks eingeführt wird. Sie eignet sich hervorragend als Grundlagenliteratur für Tensoren, virtuelle Indizes, Bond-Dimensionen und zweidimensionale Kontraktionsverfahren.
  • Román Orús: Tensor networks for complex quantum systems, Nature Reviews Physics 1, 538–550, 2019.
    • Die Arbeit bietet einen kompakten Überblick über moderne Tensor-Netzwerk-Methoden, Symmetrien, Fermionen, zweidimensionale Vielteilchensysteme und weitere Anwendungen. Für eine wissenschaftliche PEPO-Abhandlung eignet sie sich als Überblicksquelle, um PEPOs in das größere Feld der Tensor-Netzwerk-Methoden und der klassischen Simulation komplexer Quantensysteme einzuordnen.
  • J. Ignacio Cirac, David Pérez-García, Norbert Schuch und Frank Verstraete: Matrix Product States and Projected Entangled Pair States: Concepts, Symmetries, and Theorems, Reviews of Modern Physics 93, 045003, 2021.
    • Diese umfassende Review behandelt die mathematische Struktur von MPS und PEPS, Symmetrien, Verschränkung, Parent-Hamiltonians und zentrale theoretische Resultate. Obwohl der Schwerpunkt auf Zustandsnetzwerken liegt, liefert sie eine besonders solide Grundlage für die theoretische Einordnung von PEPOs. Sie sollte verwendet werden, wenn Begriffe wie virtuelle Symmetrie, Tensorstruktur, lokale Projektionen und strukturelle Eigenschaften von PEPS präzise begründet werden sollen.

Spezialisierte Arbeiten zu PEPO-Kontraktion und Operatornetzwerken

  • Román Orús: Exploring corner transfer matrices and corner tensors for the classical simulation of quantum lattice systems, Physical Review B 85, 205117, 2012.
    • Diese Veröffentlichung ist für das PEPO-Thema besonders wertvoll, weil sie Projected Entangled Pair Operators ausdrücklich in Verbindung mit Corner-Transfer-Matrizen, Zeitentwicklungsoperatoren und zweidimensionalen Tensor-Netzwerk-Verfahren diskutiert. Die Arbeit eignet sich als Spezialliteratur für die Kapitel über Netzwerkkontraktion, effektive Umgebungen, thermodynamischen Limes und die numerische Auswertung zweidimensionaler Operatornetzwerke.
  • Álvaro M. Alhambra und J. Ignacio Cirac: Locally Accurate Tensor Networks for Thermal States and Time Evolution, PRX Quantum 2, 040331, 2021.
    • Diese Arbeit gehört zu den direktesten theoretischen Quellen für PEPOs. Die Autoren konstruieren Projected Entangled Pair Operators zur Approximation thermischer Erwartungswerte und der Heisenberg-Zeitentwicklung lokaler Observablen. Besonders relevant ist das Resultat, dass für lokal genaue Approximationen Bond-Dimensionen betrachtet werden können, deren erforderliche Größe von Temperatur beziehungsweise Zeit, nicht jedoch unmittelbar von der Gesamtsystemgröße abhängt. Die Publikation ist zentrale Primärliteratur für Abschnitte über lokale Genauigkeit, thermische Zustände und Zeitentwicklung.

Thermische Zustände und imaginäre Zeitentwicklung

  • Piotr Czarnik, Lukasz Cincio und Jacek Dziarmaga: Projected Entangled Pair States at Finite Temperature: Imaginary Time Evolution with Ancillas, Physical Review B 86, 245101, 2012.
    • Die Arbeit behandelt thermische Zustände zweidimensionaler Quantensysteme durch imaginäre Zeitentwicklung eines PEPS mit Hilfsfreiheitsgraden. Obwohl hier eine Purifikationsdarstellung verwendet wird, ist die Publikation für PEPOs unmittelbar relevant, weil sie einen alternativen Tensor-Netzwerk-Zugang zur Darstellung thermischer Dichteoperatoren zeigt. Sie eignet sich besonders für einen Vergleich zwischen direkter PEPO-Darstellung einer Dichtematrix und Purifikationsverfahren.
  • András Molnár, Norbert Schuch, Frank Verstraete und J. Ignacio Cirac: Approximating Gibbs states of local Hamiltonians efficiently with projected entangled pair states, Physical Review B 91, 045138, 2015.
    • Diese theoretische Arbeit analysiert systematisch, unter welchen Bedingungen Gibbs-Zustände lokaler Hamiltonoperatoren effizient durch Projected-Entangled-Pair-Strukturen approximiert werden können. Sie ist wichtig für die mathematische Begründung thermischer Tensor-Netzwerk-Ansätze und liefert Hintergrund für die Frage, wie Bond-Dimension, Temperatur, Systemgröße und Approximationsfehler zusammenhängen.
  • Piotr Czarnik und Jacek Dziarmaga: Projected Entangled Pair States at Finite Temperature: Iterative Self-Consistent Bond Renormalization for Exact Imaginary Time Evolution, Physical Review B 92, 035120, 2015.
    • Diese Arbeit konzentriert sich auf die kontrollierte Reduktion wachsender Bond-Dimensionen während der imaginären Zeitentwicklung. Sie ist für PEPO-Rechnungen relevant, weil auch bei der Multiplikation und Evolution von Operatornetzwerken virtuelle Dimensionen anwachsen und anschließend komprimiert werden müssen. Die Publikation eignet sich daher als Spezialquelle zu Bond-Renormierung, Umgebungsinformation und numerischer Fehlerkontrolle.
  • Bram Vanhecke, Laurens Vanderstraeten und Frank Verstraete: Symmetric cluster expansions with tensor networks, Physical Review A 103, L020402, 2021.
    • Die Arbeit entwickelt Cluster-Expansionen für Exponentialoperatoren lokaler Hamiltonians. Das Verfahren ist eine wichtige Alternative zu reinen Suzuki-Trotter-Zerlegungen und erhält räumliche sowie interne Symmetrien auf natürliche Weise. Für PEPOs ist diese Publikation insbesondere für die Konstruktion von Zeitentwicklungs- und thermischen Operatoren sowie für höhere Ordnungen und komplexere Wechselwirkungsstrukturen relevant.

Offene Quantensysteme und Infinite PEPO

  • C. Mc Keever und M. H. Szymańska: Stable iPEPO Tensor-Network Algorithm for Dynamics of Two-Dimensional Open Quantum Lattice Models, Physical Review X 11, 021035, 2021.
    • Diese Publikation ist eine zentrale Spezialquelle für Infinite Projected Entangled Pair Operators. Der iPEPO-Ansatz wird direkt für die Dynamik und stationäre Zustände offener zweidimensionaler Quantengitter im thermodynamischen Limes eingesetzt. Die Arbeit ist besonders geeignet für Abschnitte über Dichtematrizen, Lindblad-Dynamik, dissipative Systeme, Trunkierung und PEPO-Methoden für offene Quantensysteme.

Aktuelle Spezialliteratur zu fermionischen PEPOs

  • Sander De Meyer, Atsushi Ueda, Yuchi He, Nick Bultinck und Jutho Haegeman: Lowering the temperature of two-dimensional fermionic tensor networks with cluster expansions, Physical Review B 114, 175113, 2026.
    • Diese aktuelle Arbeit erweitert Cluster-Expansionen auf zweidimensionale fermionische Systeme und konstruiert ausdrücklich PEPO-Approximationen von Gibbs-Zuständen. Zusätzlich werden unterschiedliche Trunkierungsverfahren beim Multiplizieren mehrerer PEPO-Schichten untersucht. Die Publikation eignet sich als aktuelle Spezialliteratur zu fermionischen PEPOs, endlichen Temperaturen, Gibbs-Zuständen und der praktischen Kompression mehrschichtiger Operatornetzwerke.

Bücher und Monographien

Monographien zu PEPOs im engeren Sinn sind wesentlich seltener als wissenschaftliche Artikel. Für eine fundierte Abhandlung sind deshalb Werke besonders wertvoll, die Tensor-Netzwerk-Kontraktionen, Quanteninformation, Vielteilchenphysik, statistische Mechanik und offene Quantensysteme systematisch behandeln. Die folgenden Titel bilden die theoretische Grundlage, auf der PEPO-Verfahren mathematisch und physikalisch verstanden werden können.

Tensor-Netzwerke und Kontraktionsmethoden

  • Shi-Ju Ran, Emanuele Tirrito, Cheng Peng, Xi Chen, Luca Tagliacozzo, Gang Su und Maciej Lewenstein: Tensor Network Contractions: Methods and Applications to Quantum Many-Body Systems, Lecture Notes in Physics 964, Springer, 2020.
    • Dieses Open-Access-Fachbuch ist für eine PEPO-Abhandlung besonders empfehlenswert. Es behandelt Tensor-Netzwerke von den mathematischen Grundlagen bis zu zweidimensionalen Kontraktionsalgorithmen und höherdimensionalen Quantengittern. Die Kapitel zu zweidimensionalen Tensor-Netzwerken sind unmittelbar relevant für PEPO-Kontraktionen, Umgebungsapproximationen, Renormierung und die numerische Behandlung großer Gittersysteme.

Standardwerke zur Quanteninformation

  • Michael A. Nielsen und Isaac L. Chuang: Quantum Computation and Quantum Information, 10th Anniversary Edition, Cambridge University Press, 2010.
    • Das Werk liefert die Standardgrundlagen zu Hilberträumen, Dichtematrizen, Quantenoperationen, Quantenkanälen, Verschränkung und Quanteninformation. Diese Konzepte sind für das Verständnis eines PEPO als Darstellung eines Vielteilchenoperators oder Dichteoperators unverzichtbar. Das Buch sollte insbesondere für grundlegende Definitionen und quanteninformationstheoretische Begriffe herangezogen werden.

Hintergrundliteratur zur Vielteilchenphysik

  • Henrik Bruus und Karsten Flensberg: Many-Body Quantum Theory in Condensed Matter Physics: An Introduction, Oxford University Press, 2004.
    • Das Buch vermittelt die physikalischen Grundlagen quantenmechanischer Vielteilchensysteme. Für PEPOs ist es insbesondere als Hintergrundliteratur zu lokalen Freiheitsgraden, Wechselwirkungen, Vielteilchenoperatoren und zeitabhängiger Quantentheorie geeignet. Es hilft, Tensor-Netzwerk-Methoden nicht nur als numerisches Verfahren, sondern im Kontext realer Vielteilchenmodelle zu verstehen.
  • Subir Sachdev: Quantum Phase Transitions, 2. Auflage, Cambridge University Press, 2011.
    • Das Standardwerk behandelt Quantenphasenübergänge, kritische Systeme und das Verhalten bei endlichen Temperaturen. Für eine PEPO-Abhandlung ist es vor allem dann relevant, wenn Tensor-Netzwerke zur Untersuchung von Quantenphasen, Korrelationslängen und kritischem Verhalten eingesetzt werden. Es bietet die physikalische Grundlage zur Interpretation numerischer PEPS- und PEPO-Ergebnisse.

Statistische Mechanik und Corner-Transfer-Methoden

  • Rodney J. Baxter: Exactly Solved Models in Statistical Mechanics, Academic Press, 1982.
    • Baxters Monographie ist eine klassische Quelle zu zweidimensionalen Gittermodellen und Corner-Transfer-Matrizen. Diese Methoden bilden einen wichtigen historischen und mathematischen Hintergrund moderner CTM- und CTMRG-Verfahren, die bei der approximativen Kontraktion von PEPS- und PEPO-Netzwerken eingesetzt werden. Das Buch eignet sich daher besonders für die theoretische Einordnung von Corner-Transfer-Verfahren.

Offene Quantensysteme

  • Heinz-Peter Breuer und Francesco Petruccione: The Theory of Open Quantum Systems, Oxford University Press, 2002.
    • Das Werk ist eine Standardreferenz für offene Quantensysteme, Mastergleichungen, dissipative Dynamik und Dichteoperatoren. Es ist besonders relevant für die wissenschaftliche Einordnung von iPEPO-Verfahren, die Dichtematrizen zweidimensionaler offener Quantengitter direkt entwickeln. Für Kapitel zu Lindblad-Dynamik, Dissipation und stationären Zuständen bietet dieses Buch die notwendige physikalische Grundlage.

Online-Ressourcen und Datenbanken

Für PEPOs ist eine laufende Literaturrecherche besonders wichtig, da viele methodische Weiterentwicklungen zunächst als Preprints erscheinen und später in Fachzeitschriften veröffentlicht werden. Online-Datenbanken sollten deshalb nicht als Ersatz für Primärliteratur, sondern als Werkzeug zur Identifikation, Überprüfung und Nachverfolgung relevanter Arbeiten verwendet werden.

Vorlesungsnotizen und Monographie-nahe Ressourcen

  • Jacob Biamonte: Lectures on Quantum Tensor Networks, arXiv, 2019.
    • Die umfangreichen Vorlesungsnotizen behandeln Tensor-Netzwerke als gemeinsame Sprache für Quantenzustände, Operatoren und Quantenprozesse. Sie sind besonders hilfreich, um die diagrammatische Darstellung von Operator-Tensor-Netzwerken zu verstehen und die Verbindung zwischen Quanteninformation und Tensor-Netzwerken herzustellen. Sie sollten vor allem als ergänzende Lern- und Hintergrundressource verwendet werden.
  • TensorNetwork.org: Review Articles and Learning Resources, laufend aktualisierte Fachressource.
    • TensorNetwork.org bündelt Übersichtsarbeiten, Lernmaterialien, Tensor-Netzwerk-Notation, Softwarehinweise und weiterführende Literatur zu MPS, MPO, PEPS und verwandten Methoden. Die Plattform eignet sich besonders zur Orientierung und zum Auffinden einschlägiger Review-Artikel. Für wissenschaftliche Aussagen sollten anschließend die dort referenzierten Originalpublikationen zitiert werden.

Fachjournale und Verlage

  • American Physical Society: Physical Review Journals.
    • Die APS veröffentlicht einen großen Teil der grundlegenden und spezialisierten Literatur zu PEPS, PEPO, Quanteninformation, Vielteilchenphysik und offenen Quantensystemen. Besonders relevant sind Physical Review Letters, Physical Review B, Physical Review X, PRX Quantum und Reviews of Modern Physics.
  • Springer Nature: SpringerLink.
    • SpringerLink bietet Fachbücher, Lecture Notes und wissenschaftliche Artikel zu Tensor-Netzwerken, mathematischer Physik, statistischer Mechanik und Quanteninformation. Besonders nützlich ist die Plattform für Monographien und systematische Darstellungen von Tensor-Kontraktionsverfahren.
  • Elsevier: ScienceDirect.
    • ScienceDirect enthält unter anderem Annals of Physics und weitere Zeitschriften, in denen grundlegende Übersichtsarbeiten zu Tensor-Netzwerken, MPS, PEPS und numerischer Vielteilchenphysik erschienen sind. Die Plattform eignet sich zur Recherche begutachteter Journalversionen und bibliographischer Angaben.
  • Nature Portfolio: Nature Reviews Physics und verwandte Fachzeitschriften.
    • Nature Portfolio veröffentlicht Übersichts- und Forschungsartikel zu Quantentechnologie, Vielteilchenphysik und modernen Tensor-Netzwerk-Anwendungen. Für PEPOs sind insbesondere Review-Artikel hilfreich, die den Stand des Feldes und Verbindungen zu Quantencomputing, Quantensimulation und stark korrelierten Systemen einordnen.

Wissenschaftliche Literaturdatenbanken

  • arXiv: Open-Access-Preprint-Repository.
    • arXiv ist für Tensor-Netzwerk-Forschung besonders wichtig, weil neue Arbeiten häufig bereits vor der endgültigen Journalveröffentlichung verfügbar sind. Für PEPO-Recherchen sind insbesondere die Bereiche Quantum Physics und Strongly Correlated Electrons relevant. Bei wissenschaftlichen Zitaten sollte nach Möglichkeit geprüft werden, ob inzwischen eine begutachtete Journalversion vorliegt.
  • INSPIRE: Wissenschaftliche Literaturdatenbank.
    • INSPIRE ist besonders nützlich für die Suche nach Veröffentlichungen an den Schnittstellen von theoretischer Physik, Quanteninformation, mathematischer Physik und Tensor-Netzwerken. Die Datenbank stellt häufig DOI, arXiv-Version und Journalangaben gemeinsam dar und erleichtert dadurch die bibliographische Überprüfung.
  • Google Scholar: Wissenschaftliche Such- und Zitationsdatenbank.
    • Google Scholar eignet sich zur Suche nach Arbeiten, Folgepublikationen und Zitationsketten. Für PEPO-Recherchen ist dies besonders hilfreich, um von grundlegenden PEPS- und PEPO-Arbeiten zu neueren Algorithmen, Anwendungen und Erweiterungen zu gelangen. Bibliographische Angaben sollten anschließend anhand der Originalzeitschrift oder des DOI überprüft werden.
  • Web of Science: Zitations- und Literaturdatenbank.
    • Web of Science ermöglicht strukturierte Zitationsanalysen und eignet sich zur Nachverfolgung wissenschaftlicher Entwicklungslinien von grundlegenden PEPS-Arbeiten bis zu modernen PEPO-Anwendungen. Die Datenbank ist insbesondere für eine systematische Literaturrecherche und die Überprüfung häufig zitierter Arbeiten sinnvoll.
  • Scopus: Abstract- und Zitationsdatenbank.
    • Scopus bietet eine breite Abdeckung wissenschaftlicher Veröffentlichungen und kann ergänzend verwendet werden, um Zitationsbeziehungen, Autorenprofile und thematisch verwandte Arbeiten zu Tensor-Netzwerken und Quantentechnologie zu identifizieren.

Lern- und Forschungsplattformen

  • TensorNetwork.org: Tensor Network Algorithms, Theory and Software.
    • Die Plattform bietet strukturierte Einführungstexte zu Tensor-Diagrammen, Symmetrien, MPS, MPO, PEPS und Tensor-Renormierungsverfahren. Sie ist insbesondere hilfreich, um die Verbindung zwischen MPO und PEPO sowie grundlegende Kontraktionskonzepte nachzuvollziehen. Für die wissenschaftliche Abhandlung dient sie primär als Lern- und Recherchehilfe, nicht als Ersatz für Primärquellen.
  • ITensor: Tensor-Network-Software und Dokumentation.
    • ITensor ist eine etablierte Softwareumgebung für Tensor-Netzwerk-Berechnungen. Die Dokumentation behandelt Tensor-Kontraktionen, Singularwertzerlegungen, Quantenzahlsymmetrien, MPS und MPO. Auch wenn viele bereitgestellte Algorithmen auf eindimensionale MPS/MPO-Probleme ausgerichtet sind, ist die Dokumentation wertvoll, um die numerischen Grundoperationen zu verstehen, aus denen auch höherdimensionale PEPO-Verfahren aufgebaut werden.
European Quantum Technology Network (EWTN)

European Quantum Technology Network (EWTN)

...

D-Mesonen

D-Mesonen

...

Einzel-Spin-Qubits

Einzel-Spin-Qubits

...

Tetraquarks

Tetraquarks

...

Quantum Information Processing (QIP)

Quantum Information Processing (QIP)

...

Bold Diagrammatic Monte Carlo (BDMC)

Bold Diagrammatic Monte Carlo (BDMC)

...