Exact chiral symmetry with quantum signal processing

Technische Zusammenfassung: Exakte chirale Symmetrie mit Quantum Signal Processing

Problemstellung
Die Quantensimulation von nichtperturbativen, Nichtgleichgewichts-Observablen in der Gitter-QCD (Quantenchromodynamik) steht vor erheblichen Herausforderungen hinsichtlich des Fermion-Verdopplungsproblems (Fermion Doubling) und der Bewahrung der chiralen Symmetrie in Hamiltonian-Formulierungen. Während die Ginsparg-Wilson (GW)-Relation und der Overlap-Operator elegante Lösungen im euklidischen Raumzeit-Kontext bieten, fehlte bisher ein kanonischer Hamiltonian-Analog für Overlap-Fermionen. Bestehende Ansätze, wie etwa Domain-Wall-Fermionen, lassen sich gut auf Hamiltonian-Settings abbilden, erfordern jedoch eine explizite zusätzliche Dimension, was die Qubit-Kosten erhöht. Im Gegensatz dazu bieten Overlap-Fermionen eine Formulierung in den physischen Raumdimensionen, beinhalten jedoch hochgradig nichtlokale All-zu-Alle-Wechselwirkungen, spezifisch die Vorzeichenfunktion des Wilson-Dirac-Hamiltonians, $\varepsilon(h_W)$, die auf Quantenhardware schwierig effizient zu implementieren ist. Das zentrale Problem, das hier adressiert wird, ist die Konstruktion effizienter Quantenalgorithmen für chirale symmetrische Gitter-Fermionen, die die Kosten für Speicher (Qubits) gegen die Gatter-Komplexität abwägen und gleichzeitig eine exakte (oder kontrolliert gebrochene) chirale Symmetrie beibehalten.

Methodik
Die Autoren schlagen einen Quantum Signal Processing (QSP)-Algorithmus zur Simulation des Overlap-Fermion-Hamiltonians vor. Die Methodik umfasst folgende wesentliche Schritte:

1. Hamiltonian-Formulierung: Die Arbeit nutzt den Overlap-Hamiltonian $H_{ov} = \psi^\dagger (\Gamma_0 + \varepsilon(h_W)) \psi$, wobei $h_W$ der Single-Particle Wilson-Dirac-Hamiltonian ist. Die Vorzeichenfunktion $\varepsilon(X) = X/\sqrt{X^\dagger X}$ wird durch ein Polynom $p_M(x)$ des Grades $M$ approximiert.
2. Block-Kodierung (Block Encoding): Die Autoren konstruieren eine Block-Kodierung des Single-Particle Wilson-Hamiltonians $\hat{h}_W$ (einschließlich der Eichfelder) unter Verwendung eines „Prepare“-Operators $P$ und eines „Select“-Operators $S$. Dies kodiert den Hamiltonian in einen unitären Operator $U_{\hat{h}_W}$, der auf einem erweiterten Hilbert-Raum mit Ancilla-Qubits wirkt.
3. Quantum Signal Processing (QSP): Zur Implementierung der Vorzeichenfunktion-Approximation verwenden die Autoren QSP (speziell Quantum Singular Value Transformation). Dies ermöglicht die Anwendung eines Polynoms vom Grad $M$, $p_M(U_{\hat{h}_W})$, auf den block-kodierten Operator. Das Polynom wird gewählt, um die Vorzeichenfunktion innerhalb einer Spektrallücke $[\kappa, 1]$ mit einem Fehler $\epsilon_e$ zu approximieren.
4. Zeitevolution: Sobald der Overlap-Hamiltonian block-kodiert ist, wird der Zeitevolutionsoperator $e^{-iHt}$ mittels QSP approximiert, wobei die Anzahl der Gatter mit den Block-Kodierungs-Kosten und der Evolutionszeit skaliert.

Zentrale Beiträge und Ergebnisse

• Algorithmische Konstruktion: Die Arbeit liefert einen konkreten QSP-basierten Algorithmus für den Overlap-Hamiltonian, der die GW-Relation bis zu einem kontrollierbaren Fehler $\epsilon_e$ bewahrt. Es wird gezeigt, dass der modifizierte chirale Operator $\hat{\gamma}_5$ mit dem approximierten Hamiltonian bis auf einen Fehler der Ordnung $\epsilon_e$ kommutiert.
• Komplexitätsanalyse:
  • Gatter-Komplexität: Die Kosten für die Block-Kodierung des Overlap-Hamiltonians skalieren als $O(Q \kappa^{-1} \log(1/\epsilon_e))$, wobei $Q$ die Anzahl der Gitterplätze (multipliziert mit den internen Freiheitsgraden) ist und $\kappa$ mit der Spektrallücke zusammenhängt. Die gesamte Gatter-Komplexität für die Zeitevolution skaliert als $O(Q \kappa^{-1} \log(1/\epsilon_e) [Qt + \log(1/\epsilon_t)])$.
  • Qubit-Kosten: Der Algorithmus benötigt $O(Q + \log Q)$ Qubits. Dies ist eine signifikante Reduktion gegenüber Domain-Wall-Fermionen, die $O(Q L_5)$ Qubits benötigen, wobei $L_5$ die Ausdehnung der zusätzlichen Dimension darstellt.
  • Vergleich: Während Domain-Wall-Fermionen von geometrischer Lokalität profitieren (was eine nahezu lineare Kostenstruktur $O(Q t \text{polylog})$ via Suzuki-Trotter oder ähnlichen Methoden erlaubt), verursacht die Overlap-Formulierung eine höhere Gatter-Tiefe aufgrund der nichtlokalen Natur der Vorzeichenfunktion-Approximation. Dennoch bietet der Overlap-Ansatz eine überlegene asymptotische Skalierung im Speicherbedarf.
• Physikalische Interpretation von QSP: Die Autoren zeigen, dass der benötigte Polynom-Grad $M$, um einen Fehler $\epsilon_e$ zu erreichen, als $M = O(\kappa^{-1} \log(1/\epsilon_e))$ skaliert. Sie identifizieren diese Skalierung mit der Größe der zusätzlichen Dimension $L_5$ bei Domain-Wall-Fermionen, wobei der Fehler als $e^{-c L_5}$ skaliert. Somit „konstruiert“ die QSP-Implementierung effektiv eine zusätzliche Dimension durch die Tiefe des Schaltkreises und spiegelt damit die physikalische Korrespondenz zwischen dem Overlap-Operator und der Randtheorie von Domain-Wall-Fermionen wider.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass Quantensimulationen von Dirac-Fermionen mit exakter chiraler Symmetrie im Sinne von „fast kostenlos“ sind, da das Anwenden des Overlap-Hamiltonians nur einen logarithmischen Faktor mehr (in Bezug auf die Fehlertoleranz) kostet als der Wilson-Dirac-Hamiltonian.

Die zentrale Erkenntnis ist, dass QSP eine quantenalgorithmische Realisierung der bekannten Korrespondenz zwischen dem Overlap-Operator und der zusätzlichen Dimension von Domain-Wall-Fermionen bereitstellt. Es besteht ein Trade-off:

• Domain-Wall-Fermionen: Höhere Qubit-Kosten ($O(Q L_5)$), aber geringere Gatter-Tiefe aufgrund von Lokalität.
• Overlap-Fermionen (via QSP): Niedrigere Qubit-Kosten ($O(Q)$), aber höhere Gatter-Tiefe ($O(Q^2)$ Worst-Case-Skalierung für die Zeitevolution aufgrund von Nichtlokalität).

Die Autoren kommen zu dem Schluss, dass die Wahl zwischen den Formulierungen von den Hardware-Beschränkungen abhängt (Anzahl der Qubits vs. Schaltungstiefe). Sie merken an, dass die klassische Gitter-QCD oft rationale Approximationen (z. B. Zolotarev) für die Vorzeichenfunktion verwendet, die effizienter sind, diese jedoch keine direkten Analoga im Standard-QSP besitzen, was einen Weg für zukünftige Arbeiten aufzeigt. Das Paper stellt zudem klar, dass die Skalierungsergebnisse robust gegenüber der spezifischen Kodierung der Eichfelder sind, wenngleich konstante Faktoren variieren können.