Spectral functions on a quantum computer through system-environment interaction

Dieser Beitrag stellt einen effizienten Quantenalgorithmus vor, der System-Umgebungs-Wechselwirkungen modelliert, um Spektralfunktionen mit einer $O(N)$-Reduktion des Stichprobenaufwands im Vergleich zu Standardtechniken zu messen, und demonstriert seine Wirksamkeit an einem 27-Platz-System unter Verwendung von 54 Qubits auf einem Quantencomputer mit Ionenfalle von Quantinuum.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, die „musikalischen Noten" (Energieniveaus) zu verstehen, die ein komplexes Material spielen kann. In der realen Welt verwenden Wissenschaftler eine High-Tech-Kamera namens ARPES (Winkelaufgelöste Photoemissionsspektroskopie), um ein Bild dieser Noten zu machen. Dazu schießen sie Licht auf das Material, schlagen Elektronen heraus und messen dann, wie schnell und in welche Richtung diese Elektronen fliegen.

Das Problem ist, dass die Simulation dieses Prozesses auf einem Computer unglaublich schwierig ist. Es ist wie der Versuch, den Klang einer Symphonie vorherzusagen, indem man jedes einzelne Instrument nacheinander in absoluter Stille anhört und dann versucht, das ganze Lied zu erraten. Auf einem Quantencomputer war die alte Methode wie die Aufforderung an einen Musiker, eine Note zu spielen, zu stoppen, zurückzusetzen, die nächste Note zu spielen, zu stoppen und erneut zurückzusetzen. Wenn Sie 1.000 Instrumente (oder „Orte" im Material) haben, müssen Sie diesen Vorgang 1.000 Mal wiederholen, nur um ein vollständiges Bild zu erhalten. Das dauert ewig und verschwendet eine enorme Menge an Zeit.

Die neue Idee: Eine „gefälschte" Umgebung

Die Autoren dieses Papers haben einen cleveren Abkürzungsweg gefunden. Anstatt den Computer zu bitten, die Noten einzeln zu berechnen, beschlossen sie, das tatsächliche Experiment direkt auf dem Quantencomputer zu simulieren.

Stellen Sie es sich so vor:

• Das System: Dies ist das Material, das Sie untersuchen möchten (das Orchester).
• Die Umgebung: Dies ist die „Kamera" oder das „Vakuum", das die Elektronen auffängt (das Publikum).

In ihrer neuen Methode verbinden sie das „Orchester" mit einem „gefälschten Publikum" (einer Umgebung) innerhalb des Computers. Sie lassen das Orchester für kurze Zeit mit diesem Publikum interagieren. Anstatt das Orchester direkt zu messen, schauen sie einfach auf das Publikum, um zu sehen, wer eine Note aufgefangen hat.

Da das Publikum gleichzeitig mit dem gesamten Orchester verbunden ist, verrät eine einzige Messung ihnen die „Noten" für das gesamte Orchester gleichzeitig.

Der große Gewinn: Geschwindigkeit und Effizienz

Das Paper behauptet, dies sei ein Wendepunkt für eine bestimmte Art von Quantencomputer, den Ionenfallen-Computer (der gefangene Atome als Qubits verwendet).

• Der alte Weg: Um ein klares Bild zu erhalten, müssten Sie möglicherweise 1.000 Fotos (Messungen) machen, weil die Kamera langsam und unscharf ist.
• Der neue Weg: Sie benötigen nur ein Foto.

Die Autoren sagen, dies spare eine enorme Menge an Zeit. Wenn die alte Methode 100 Stunden dauerte, könnte diese neue Methode nur 1 Stunde dauern. Sie nennen dies eine O(N)-Verbesserung, was bedeutet, dass wenn Sie die Größe des untersuchten Materials verdoppeln, die alte Methode doppelt so langsam wird, diese neue Methode jedoch genauso schnell bleibt.

Der Haken: Sie benötigen mehr „Qubits"

Es gibt einen Kompromiss. Um diesen Trick auszuführen, müssen Sie die Anzahl der „Qubits" (die Grundeinheiten des Quantencomputers) verdoppeln, da Sie sowohl das Material als auch die gefälschte Umgebung simulieren müssen. Es ist wie der Bedarf an einem größeren Raum, um sowohl die Band als auch das Publikum unterzubringen. Die Autoren argumentieren jedoch, dass für diese spezifischen Computer die Zeitersparnis bei den Messungen viel wichtiger ist als ein paar zusätzliche Qubits.

Der „magische" Trick: Die Fermionische Fourier-Transformation

Um die „gefälschte Umgebung" funktionsfähig zu machen, muss der Computer einen komplexen mathematischen Tanz namens Fermionische Fourier-Transformation (FFT) aufführen. Stellen Sie sich vor, Sie mischen ein Kartenspiel so, dass alle Herzen zusammen, alle Pik zusammen usw. liegen, aber dies auf eine Weise tun, die den seltsamen Regeln von Quantenteilchen (Fermionen) gerecht wird.

Die Autoren haben nicht einfach einen Standard-Mix verwendet; sie erfanden eine effizientere Art, diese spezifischen Quantenkarten zu mischen, insbesondere für ein Setup, bei dem die Anzahl der Karten keine Potenz von 2 ist (wie 27 Karten). Sie testeten diesen Mix auf einer echten Maschine (Quantinuums H2) und bewiesen, dass er funktioniert.

Der Realwelt-Test

Das Team schrieb nicht nur Theorie; sie führten das Experiment auf einem echten Quantencomputer mit 54 Qubits durch (27 für das Material, 27 für die Umgebung). Sie maßen erfolgreich die „Spektralfunktion" (die musikalischen Noten) einer 27-Orte-Kette von Teilchen.

Obwohl der echte Computer etwas „Rauschen" hat (wie statisches Rauschen im Radio), waren die Ergebnisse klar genug, um die Hauptmerkmale des Materials zu erkennen. Das „Rauschen" machte das Signal etwas schwächer, verzerrte aber nicht die Form der Noten, was bedeutet, dass die gesuchte Physik genau blieb.

Zusammenfassung

Kurz gesagt führt dieses Paper eine neue Methode ein, um zu simulieren, wie Materialien mit Licht interagieren. Indem sie das gesamte Experiment (System + Umgebung) gleichzeitig simulieren, anstatt Teile davon separat zu berechnen, können sie die Antwort N-mal schneller erhalten (wobei N die Größe des Systems ist). Dies macht es auf heutigen Quantencomputern, insbesondere der Ionenfallen-Art, viel praktikabler, große, komplexe Materialien zu untersuchen.

Technische Zusammenfassung

Problemstellung
Spektralfunktionen, die mittels winkelauflösender Photoemissionsspektroskopie (ARPES) gemessen werden, sind entscheidend für die Aufklärung der Bandstruktur von Materialien. Theoretisch entsprechen diese Funktionen dynamischen Einpunktfunktionen in einem Gleichgewichtszustand. Obwohl ihre Berechnung klassisch herausfordernd ist, stehen ihrer Messung auf Quantencomputern erhebliche Hürden entgegen. Standardansätze, die die Berechnung dynamischer Korrelationen $\langle c^\dagger_j(t)c_0(0)\rangle$ und die Durchführung von Fourier-Transformationen beinhalten, leiden unter einem hohen Stichproben-Overhead. Da die Operatoren $c_j + c^\dagger_j$ für verschiedene Gitterplätze $j$ nicht kommutieren, kann pro Messschuss (Shot) nur eine Korrelation gemessen werden. Um die Spektralfunktion für alle $N$ Impulse wiederherzustellen, ist dies mit einem Stichproben-Overhead von $O(N)$ verbunden. Dies ist insbesondere für Ionenfallen-Quantencomputer problematisch, die zwar eine hohe Gate-Fidelity und Konnektivität aufweisen, aber relativ langsame Shot-Zeiten haben, wodurch die $O(N)$-Anzahl an Schüssen zu einem Engpass für Anwendungen im Nutzenmaßstab wird.

Methodik
Die Autoren schlagen eine Methode vor, um Spektralfunktionen $A(k, \omega)$ direkt zu messen, indem sie die physikalische Wechselwirkung zwischen dem System und der Umgebung modellieren, wie sie in ARPES-Experimenten auftritt. Der Ansatz umfasst:

1. System-Umgebung-Kopplung: Das System von Interesse (Hamiltonoperator $H_{sys}$) wird über einen Wechselwirkungsterm $H_{int}$ mit einer Umgebung von $N$ Gitterplätzen (Hamiltonoperator $H_{env}$) gekoppelt. Der totale Hamiltonoperator lautet $H = H_{sys} + \epsilon H_{int} + \omega H_{env}$.
2. Zustandspräparation und Evolution: Das System wird in einem Eigenzustand $|E\rangle$ von $H_{sys}$ initialisiert, während die Umgebung in einem Vakuumzustand $|0\rangle$ (für $A_+$) oder einem vollständig gefüllten Zustand $|all\rangle$ (für $A_-$) initialisiert wird. Das Gesamtsystem entwickelt sich unter $H$ für eine feste Zeit $t$.
3. Messung: Anstatt dynamische Korrelationen zu messen, misst das Protokoll die Impulsbesetzungszahl $n(k)$ in der Umgebung. Dies erfordert die Durchführung einer fermionischen Fourier-Transformation (FFT) auf den Umgebung-Qubits am Ende der Evolution, gefolgt von der Messung lokaler $Z$-Pauli-Matrizen.
4. Theoretisches Ergebnis: Im Grenzfall schwacher Kopplung ($\epsilon \to 0$) und langer Zeit ($t \to \infty$) ist der Erwartungswert $\langle n(k) \rangle$ proportional zur Spektralfunktion $A(k, \omega)$ (gefaltet mit einem Kernel, der für $t \to \infty$ gegen eine Delta-Funktion konvergiert).
5. Implementierung der fermionischen FFT: Um beliebige Systemgrößen (insbesondere $N=27$) zu handhaben, entwickeln die Autoren eine effiziente Gate-Zerlegung für eine generische radix-$n$ fermionische FFT. Sie nutzen eine rekursive Struktur, die „Interleave"-Operationen beinhaltet, um nicht-adjazente fermionische Modi für die Transformation benachbart zu machen, optimiert durch Graph-Reduktionsverfahren für Architekturen mit All-zu-All-Konnektivität.

Hauptbeiträge

• Effizientes Sampling: Der primäre Beitrag ist eine Reduktion des Stichproben-Overheads. Durch die Messung der Impulsbesetzung der Umgebung erhält die Methode Spektraldaten für alle $N$ Impulse $k$ in einem einzigen Schuss und reduziert die erforderliche Anzahl an Schüssen um einen Faktor von $O(N)$ im Vergleich zu Standardansätzen für dynamische Korrelationen.
• Laufzeitverbesserung: Diese Reduktion der Schüsse führt zu einer $O(N)$-Verbesserung der gesamten Laufzeit, was für Plattformen wie Ionenfallen-Quantencomputer entscheidend ist, bei denen die Shot-Zeit der limitierende Faktor ist.
• Robustheit gegenüber Trotter-Fehlern: Die Autoren argumentieren, dass ihr Ansatz robuster gegenüber Evolutionsfehlern (Trotterisierung) ist als Standardmethoden. Während Standardmethoden aufgrund von Fehlern unphysikalische negative Werte für Spektralfunktionen liefern können, misst die vorgeschlagene Methode eine Besetzungszahl, die selbst bei unvollkommener Evolution eine positive Größe bleibt.
• Hardware-Benchmarking: Die Autoren implementierten eine radix-3 fermionische FFT auf 27 Qubits und verglichen verschiedene Kompilierungsstrategien (lokale CZs, CX-Leitern, Graph-Reduktion) auf den Ionenfallen-Systemen H2-1 und H2-2 von Quantinuum. Sie identifizierten die Kompilierung mit $\log_2(N)$-Skalierung als diejenige, die die niedrigste Energie und In-Fidelity liefert.
• Demonstration des vollständigen Algorithmus: Der vollständige Algorithmus wurde auf einem 54-Qubit-System (27 Systemplätze + 27 Umgebungsplätze) auf der Quantinuum H2-2-Hardware demonstriert und berechnete erfolgreich die Spektralfunktion eines eindimensionalen freien-Fermionen-Hamiltonoperators.

Ergebnisse

• Hardware-Leistung: Auf der H2-2-Hardware reproduzierte die Methode erfolgreich die Form des Spektralbands für ein 27-Platz-System. Obwohl Hardware-Rauschen eine allgemeine Abschwächung des Signals verursachte (ein „Rauschboden" nahe 0,5), blieben die physikalischen Merkmale, wie Bandpositionen und -breiten, unverzerrt.
• Fehleranalyse: Vergleiche mit exakten Simulationen zeigten, dass die vorgeschlagene Methode mit deutlich weniger Trotter-Schritten als der Standardansatz für dynamische Korrelationen die gleiche Präzision erreicht, insbesondere für kleine Kopplungen $\epsilon$.
• Kostenanalyse: Für das spezifische 27-Platz-Experiment verursachte die Methode einen Overhead an Zwei-Qubit-Gates (aufgrund der Umgebung und der FFT), sparte jedoch einen Faktor von $\sim 27$ an der Anzahl der Schüsse. Die Autoren schätzen, dass ohne diese Methode die Gesamtkosten der Laufzeit bei gleicher Präzision etwa 11-mal höher gewesen wären.
• Effekte endlicher Kopplung: Theoretische Analysen und Simulationen freier Fermionen zeigen, dass eine endliche Kopplung $\epsilon$ zu einer Bandverbreiterung und „Geisterbändern" führt. Für $\epsilon t \lesssim \pi$ sind diese Effekte jedoch vernachlässigbar, und die Methode approximiert die Spektralfunktion getreu.

Bedeutung
Die Arbeit behauptet, dass dieser Ansatz die Skalierbarkeit der Berechnung von Spektralfunktionen auf Quantencomputern grundlegend verbessert. Indem ein Qubit- und Gate-Overhead gegen eine massive Reduktion der Sampling-Anforderungen getauscht wird, adressiert die Methode einen kritischen Engpass für Ionenfallen-Quantencomputer. Die Autoren gehen davon aus, dass mit zunehmender Systemgröße in Richtung des Nutzenmaßstabs (hunderte oder tausende Plätze) die $O(N)$-Reduktion der Schüsse zu einem entscheidenden Faktor für die praktische Machbarkeit der Simulation von Materialeigenschaften auf Quantenhardware werden wird. Die Arbeit liefert zudem eine konkrete, effiziente Implementierung fermionischer Fourier-Transformationen für Systemgrößen, die keine Zweierpotenzen sind, eine notwendige Komponente für die Simulation realistischer Materialmodelle.