Nuclear scattering phase shifts with factorized geometric-time RODEO on a quantum processor

Technische Zusammenfassung: Kernstreuungsphasenverschiebungen mit faktorisierter geometrischer Zeit-RODEO auf einem Quantenprozessor

Problemstellung
Die Rekonstruktion von Kernstreuungsphasenverschiebungen aus quantencomputierten Spektren stellt eine erhebliche Herausforderung dar. Während Quanteneigenwertlöser diskrete Energieniveaus in einem beschränkten System (einer „Falle“) schätzen können, erfordert die Verbindung dieser diskreten Zustände mit Kontinuumsobservablen (Freiraumstreuung) eine klassische Nachbearbeitung, spezifisch die Modifizierte Effektive Reichweitenexpansion (MERE). Die Leistung eines Quantenschaltkreises kann in diesem Kontext nicht allein anhand des quadratischen Mittelwertfehlers (RMS-Fehler) einzelner Energieniveaus beurteilt werden. Spektrale Fehler können während der anschließenden Energieinterpolation und der Null-Einschluss-Extrapolation verstärkt, unterdrückt oder umorganisiert werden. Daher muss ein ressourcenreduzierter Quantenschaltkreis daran gemessen werden, ob er die für die endgültige Streurekonstruktion notwendigen Informationen bewahrt, anstatt nur an der spektralen Fidelität.

Methodik
Die Autoren rekonstruieren elastische s-Wellen-Neutron-Proton (np)-Phasenverschiebungen unter Verwendung von Daten eines IBM Aachen Quantenprozessors. Die Studie vergleicht zwei Implementierungen des RODEO-Algorithmus (Repetitive Operation for Discrete Energy Optimization):

1. Dynamic Direct R10: Eine Baseline-Implementierung unter Verwendung von zehn Zyklen mit festen Gaußschen Quantil-Zeitschritten. Für jede Abfrageenergie wird die kontrollierte Evolution separat neu synthetisiert und kompiliert.
2. Factorized Geometric-Time RODEO (FG–R6): Eine komprimierte Implementierung unter Verwendung von sechs Zyklen mit geometrisch abgetasteten Zeitschritten. Diese Methode verwendet eine exakte Abfrage-Phasen-Trennung, bei der die Hamilton-Evolution auf eine Referenzenergie ($E_{ref}$ ) fixiert und einmalig kompiliert wird. Die abfrageabhängige Phase ($\Delta E = E_q - E_{ref}$) wird nach der Kompilierung an das Ancilla-Qubit gebunden.

Beide Methoden nutzen ein Drei-Qubit-Register (zwei System-Qubits und ein recyceltes Ancilla-Qubit) auf dem IBM Aachen Backend. Das physikalische Modell ist eine schematische s-Wellen-np-Quadrattopf-Wechselwirkung. Klassisch konstruierte vierdimensionale Okubo–Lee–Suzuki (OLS) effektive Hamilton-Operatoren werden verwendet, um das System darzustellen. Das Experiment misst vier Positiv-Energie-Niveaus bei fünf verschiedenen Fallenstärken ($\omega = 3,6$ bis $4,0$ MeV), was zwanzig spektrale Eingaben ergibt. Diese Eingaben werden in eine klassische MERE-Analyse eingespeist, um zu den Freiraumwerten zu extrapolieren und Phasenverschiebungen über ein Energiegitter von 0,1–30,0 MeV zu berechnen.

Wesentliche Beiträge

• Schaltkreis-Kompression durch Faktorisierung: Die Arbeit führt den FG–R6 Schaltkreis ein, der die Ancilla-Wiederverwendung mit einer faktorisierten Zeitentwicklungsstrategie kombiniert. Durch die Trennung der Abfragephase von der Hamilton-Entwicklung reduziert der Schaltkreis die Anzahl der kontrollierten Evolutionszyklen von zehn auf sechs und eliminiert die Notwendigkeit, den Evolutionsoperator für jede Abfrageenergie neu zu kompilieren.
• Hardware-Benchmarking: Die Studie liefert ein vollständiges, Hardware-basiertes Benchmark für die Kernstreuung, das über Einzel-Level-Tests hinausgeht und einen vollständigen Satz von zwanzig spektralen Eingaben für die Freiraumrekonstruktion umfasst.
• Analyse des Ressourcen- vs. Genauigkeits-Trade-offs: Die Arbeit zeigt, dass die Reduzierung der Schaltkreistiefe und der Zwei-Qubit-Gatter-Anzahl nicht zwangsläufig die Genauigkeit der endgültigen Streu-Observable verschlechtert, selbst wenn der RMS-Fehler der intermediären Energieniveaus zunimmt.

Ergebnisse

• Ressourcenreduktion: Im Vergleich zur Dynamic Direct R10 Implementierung reduziert die FG–R6 Implementierung die mediane kompilierte Schaltkreistiefe um 33,1 % (von 459 auf 307) und die Zwei-Qubit-Gatter-Anzahl um 37,2 % (von 113 auf 71).
• Energieunsicherheit: FG–R6 erreicht eine kleinere mediane Finite-Shot-Energieunsicherheit (0,534 keV) im Vergleich zu R10 (0,610 keV), obwohl es eine größere mediane gefittete Amplitude und einen größeren medianen RMS-Zentrumsversatz ($\Delta E_c$) relativ zur Referenzenergie aufweist (1,648 keV gegenüber 0,793 keV).
• Genauigkeit der Streuungsphasenverschiebung: Auf dem vollständigen 0,1–30,0 MeV Gitter sinkt die maximale Abweichung der rekonstruierten Phasenverschiebung von der exakten Energie-MERE-Referenz signifikant von 0,842° (R10) auf 0,163° (FG–R6). Die maximale Abweichung von der analytischen Quadrattopf-Lösung sinkt ebenfalls von 0,981° auf 0,156°.
• Energieabhängige Performance: Die Verbesserung wird primär durch den Niedrigenergie-Extrapolationsbereich (unter 10 MeV) getrieben. In der gemeinsamen Interpolationsmenge von 10,0–30,0 MeV liefert die Direct R10 Implementierung tatsächlich etwas kleinere Residuen als FG–R6.
• Fit-Form-Sensitivität: Die Studie stellt fest, dass die Fit-Form-Sensitivität spürbar bleibt. Die Streuung der Ergebnisse über verschiedene Polynomgrade für die Energie- und Fallen-Fits ist erheblich (z. B. eine Halbrange von 20,12° auf dem vollen Gitter für R10 gegenüber 6,75° für FG–R6), was darauf hindeutet, dass die Wahl des Extrapolationsmodells das Endergebnis signifikant beeinflusst.

Bedeutung und Ansprüche
Das Papier behauptet, dass dieses reduzierte Benchmark erfolgreich die Verbindung zwischen der RODEO-Schaltkreiskompression und nuklearen Kontinuumsobservablen herstellt. Die primäre Bedeutung liegt in der Demonstration, dass die Leistung eines Schaltkreises durch die Rekonstruktion der Streu-Observable und deren spezifischen Energiebereich bewertet werden muss, statt allein durch spektrale RMS-Fehler. Die FG–R6 Implementierung wird als validierter Genauigkeits–Ressourcen–Kompromiss für das getestete Spektrum und das Fenster präsentiert, nicht als universelle Minimalzyklus-Lösung.

Die Autoren stellen explizit klar, dass diese Ergebnisse keinen Quantenvorteil oder Vielteilchen-Skalierbarkeit demonstrieren. Die Hardware-Daten bestimmen nicht separat die Effekte der Phasenkompilierung gegenüber der Zeitabtastung, noch etablieren sie die Reproduzierbarkeit unabhängiger Durchläufe. Die Arbeit dient als Ausgangspunkt für die Erweiterung dieser Methoden auf größere Modellräume, realistische Wechselwirkungen und Vielteilchen-Reaktionsdynamiken, wie etwa resonante Streuung höherer Partialwellen (z. B. Neutron–Alpha-Streuung) und Systeme mit langreichweitigen Coulomb-Wechselwirkungen.