Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

Technische Zusammenfassung: Sondierung des Grundzustands des antiferromagnetischen Heisenberg-Modells auf dem Kagome-Gitter mittels geometrisch informierter Variationaler Quanten-Eigenwert-Solver

Problemstellung
Die Arbeit befasst sich mit der Herausforderung, den Grundzustand des antiferromagnetischen Heisenberg-Modells auf einem Kagome-Gitter zu finden, ein System, das für geometrische Frustration und potenzielle Quanten-Spin-Flüssigkeitsphasen bekannt ist. Die Bestimmung des Grundzustands solcher $k$-lokaler Hamilton-Operatoren ist ein QMA-hartes Problem, was die klassische Simulation für große Systeme aufgrund des exponentiellen Wachstums des Hilbert-Raums unpraktikabel macht. Während Variationale Quanten-Eigenwert-Solver (VQE) einen hybriden quanten-klassischen Ansatz bieten, der für Noisy Intermediate-Scale Quantum (NISQ)-Geräte geeignet ist, stehen sie vor erheblichen Hürden: der Notwendigkeit tiefer Schaltkreise, der Anfälligkeit für Rauschen (kohärente und inkohärente Fehler), Barren Plateaus, in denen Gradienten verschwinden, und der Schwierigkeit, Ansätze zu entwerfen, die sowohl expressiv als auch hardwareeffizient sind. Speziell zielen die Autoren darauf ab, die Eigenschaften des Grundzustands, wie Spin-Spin-Korrelationen und statische Strukturfaktoren, auf echter Quantenhardware zu charakterisieren, ohne sich auf schwere, problemspezifische Reduktionen oder umfangreiche Fehlerminimierung während der Optimierungsschleife zu verlassen.

Methodik
Die Autoren schlagen einen geometrisch informierten VQE-Ansatz vor, der auf das Kagome-Gitter zugeschnitten ist, wobei der Fokus auf zwei fundamentalen Geometrien liegt: einem einzelnen Dreieck und einem Kagome-Stern (sechs eckenverknüpfte Dreiecke).

1. Geometrisch informierter Ansatz:
   Die zentrale Innovation ist ein maßgeschneiderter, flacher, hardwareeffizienter Ansatz. Im Gegensatz zum Hamiltonian Variational Ansatz (HVA), der für NISQ-Geräte zu tief sein kann, wird dieser Ansatz aus einem Muster parametrisierter Einzelqubit-Rotationen ($R_y$) und Verschränkungsgattern konstruiert, die auf die nativen Hardware-Konnektivitäten beschränkt sind (Vermeidung von SWAP-Gattern).

   • Euklidischer Parameterraum: Ein entscheidendes Merkmal dieses Designs ist, dass der Fubini-Study-Metrik-Tensor analytisch als diagonal und konstant bewiesen wurde. Dies stellt sicher, dass der Parameterraum natürlich euklidisch, singularitätsfrei und frei von Überparametrisierung ist. Folglich ist der Quanten-Natural-Gradient äquivalent zum Standardgradienten, was die Optimierungslandschaft vereinfacht.
   • Hardware-Kompatibilität: Die Schaltungstiefe wird minimiert (eine Verschränkungsschicht), und die Anzahl der Parameter wird durch das abwechselnde Platzieren von $R_y$-Gattern reduziert.

2. Optimierungsstrategie (AQNGD/AGD):
   Die Autoren verwenden einen adaptiven Gradientenabstieg-Optimizer (AGD). Da der Ansatz einen diagonalen und konstanten Fubini-Study-Metrik liefert, reduziert sich der Adaptive Quantum Natural Gradient Descent (AQNGD) zu AGD. Die Lernrate wird unter Verwendung der Armijo-Backtracking-Regel dynamisch angepasst. Dies ermöglicht große initiale Schritte, um die Konvergenz zu beschleunigen, während der Schrittumfang automatisch reduziert wird, wenn sich der Optimierer dem Minimum nähert, um ein Stagnieren in engen Geometrien zu vermeiden. Dies wird mit der Simultaneous Perturbation Stochastic Approximation (SPSA) verglichen.

3. Fehlerminimierung (Post-Processing):
   Die Fehlerminimierung wird nach der VQE-Optimierung angewendet, um die endgültigen Energie-Schätzungen und strukturellen Eigenschaften zu verfeinern, anstatt während der iterativen Trainingsphase.

   • Readout Error Mitigation (REM): Es wird ein partitionierter Response-Matrix-Ansatz verwendet, um SPAM-Fehler (State Preparation and Measurement) zu invertieren, was effizienter skaliert als eine vollständige Matrixinversion für dünnbesetzte Messungen.
   • Zero-Noise Extrapolation (ZNE): Die Autoren wenden ZNE unter Verwendung von quadratischer Extrapolation (und bayesianischer Polynomregression) an, um Ergebnisse auf das Nullrausch-Limit zu extrapolieren. Sie stellen fest, dass ZNE nicht strikt dem Rayleigh-Ritz-Variationsprinzip folgt, hier jedoch eingesetzt wird, um die Genauigkeit zu verbessern.

Wesentliche Ergebnisse
Die Studie wurde auf echter Quantenhardware (IBM Oslo, IBM Kyoto, IBM Torino und IBM Algiers) durchgeführt und mittels State-Vektor-Simulationen validiert.

• Konvergenz und Optimierung:

  • Auf einem 3-Qubit-Dreieckssystem konvergierte der AGD-optimierte VQE glatt zum Grundzustandsenergie ($E \approx -3$) mit hoher Fidelität (99,89 % in der Simulation) unter Verwendung von nur 19 Iterationen. Die finalen Parameter entsprachen eng den theoretischen exakten Werten.
  • Auf einem 12-Qubit-Kagome-Stern zeigte AGD eine überlegene Konvergenz im Vergleich zu SPSA und benötigte weniger Iterationen, um die Grundzustandsenergie ($E \approx -18$) zu erreichen. Die finale Fidelität in der Simulation betrug 96,37 %.
  • Die Schrittgrößen-Adaption via Armijo-Regel zeigte, dass sie effektiv durch die Kostenlandschaft navigiert, indem sie mit großen Schritten beginnt und die Schrittgröße verfeinert, sobald die Konvergenz eintritt.

• Rauschresistenz struktureller Eigenschaften:

  • Eine bedeutende Erkenntnis ist, dass Spin-Spin-Korrelationen und der statische Strukturfaktor ($S(q)$) gegenüber Rauschen resistent sind. Selbst ohne Fehlerminimierung behielten die experimentellen $S(q)$-Maps die qualitative Signatur des Grundzustands bei (z. B. die Unterdrückung bei Null-Impuls, die auf keine ferromagnetische Ordnung hindeutet).
  • Der Pearson-Korrelationskoeffizient zwischen dem verrauschten experimentellen $S(q)$ und dem exakten theoretischen $S(q)$ blieb extrem hoch (> 99,9 %).
  • Experimente in der Simulation zeigten eine „schrittweise“ Abhängigkeit des statischen Strukturfaktors vom Rauschniveau: Die Korrelationsstruktur bleibt bis zu einem kritischen Rauschschwellenwert robust, nach dem sie rapide degradiert, was einer Phasenübergang ähnelt. Dies steht im Gegensatz zur Energie, die monoton mit dem Rauschen abnimmt.

• Wirksamkeit der Fehlerminimierung:

  • Die Anwendung von REM und ZNE verbesserte die Energie-Schätzungen signifikant. Für den Kagome-Stern war die Roh-Energie $\approx -14,9$, während REM allein die Energie auf $\approx -17,3$ verbesserte und ZNE (quadratisch) den Wert näher an den exakten Wert von $-18$ brachte.
  • Die Autoren beobachteten, dass die Kombination von REM und ZNE zu einem Überschießen der Grundzustandsenergie führen kann (Verletzung der Variationsgrenze), ein Verhalten, das konsistent mit der nicht-variationalen Natur von ZNE ist. ZNE allein lieferte oft die beste Genauigkeit unter den getesteten Konfigurationen.

Bedeutung und Ansprüche
Das Paper behauptet, dass ein geometrisch informierter Ansatz-Entwurf, der einen natürlich euklidischen Parameterraum gewährleistet, ein effizientes Training von VQE auf aktueller NISQ-Hardware ermöglicht, ohne dass komplexe, problemspezifische Reduktionen erforderlich sind.

• Robustheit: Die Arbeit zeigt, dass aussagekräftige physikalische Eigenschaften, insbesondere Spin-Spin-Korrelationen und der statische Strukturfaktor, mit ausreichender Genauigkeit aus verrauschter Quantenhardware extrahiert werden können, um den Grundzustand zu charakterisieren, selbst ohne Fehlerminimierung während der Optimierungsschleife.
• Optimierungseffizienz: Die Verwendung von AGD mit adaptiven Schrittgrößen erweist sich für diesen spezifischen hardwareeffizienten Ansatz als effizienter als SPSA und erreicht eine schnellere Konvergenz in Bezug auf die Iterationen.
• Skalierbarkeit: Durch die Vermeidung tiefer Schaltkreise und schwerer Fehlerminimierung während des Trainings bietet der Ansatz einen skalierbaren Pfad zur Untersuchung frustrierter Quantenmagnete. Die Autoren schlagen vor, dass dieser Rahmen auf größere Kagome-Gitter und andere frustrierte Quantensysteme ausgeweitet werden kann.
• Bescheidenheit: Die Autoren räumen ein, dass ihr Verfahren zwar kleine fundamentale Zellen (Dreieck und Stern) erfolgreich charakterisiert, die allgemeine Lösung des unendlichen Kagome-Gitters jedoch weiterhin eine offene Frage in der Literatur bleibt. Sie behaupten nicht, das Problem des unendlichen Gitters gelöst zu haben, sondern haben eine praktikable Methodik zur Untersuchung solcher Systeme auf aktuellen Geräten demonstriert. Sie merken zudem an, dass ZNE die Energie-Schätzungen verbessert, aber das Variationsprinzip nicht bewahrt, und dass die Kombination von REM und ZNE eine sorgfältige Interpretation erfordert.