Exact chiral symmetry with quantum signal processing

Technische Samenvatting: Exacte Chirale Symmetrie met Quantum Signal Processing

Probleemstelling
De kwantumsimulatie van niet-perturbatieve, niet-evenwichtige observabelen in rooster-QCD (Lattice QCD) staat voor aanzienlijke uitdagingen met betrekking tot fermion-verdubbeling (fermion doubling) en het behoud van chirale symmetrie in Hamiltoniaanse formuleringen. Hoewel de Ginsparg-Wilson (GW) relatie en de overlap-operator elegante oplossingen bieden in Euclidische ruimtetijd, heeft een canonieke Hamiltoniaanse analogon voor overlap-fermionen een gebrek aan een directe vorm gehad. Bestaande benaderingen, zoals domain-wall fermionen, mappen goed naar Hamiltoniaanse settings maar vereisen een expliciete extra dimensie, wat de qubit-kosten verhoogt. Daarentegen bieden overlap-fermionen een formulering in de fysieke ruimtelijke dimensies, maar impliceren zij zeer niet-lokale all-to-all interacties, specifiek het tekenfunctie van de Wilson-Dirac Hamiltonian, $\varepsilon(h_W)$, wat moeilijk efficiënt te implementeren is op kwantumhardware. Het centrale probleem dat wordt geadresseerd is hoe men efficiënte kwantumalgoritmen kan construeren voor chirale-symmetrische rooster-fermionen die de balans tussen geheugen (qubit) kosten en poortcomplexiteit handhaven terwijl ze exacte (of gecontroleerd gebroken) chirale symmetrie behouden.

Methodologie
De auteurs stellen een Quantum Signal Processing (QSP) algoritme voor om de overlap-fermion Hamiltonian te simuleren. De methodologie verloopt via de volgende sleutelstappen:

1. Hamiltoniaanse Formulering: Het werk maakt gebruik van de overlap-Hamiltonian $H_{ov} = \psi^\dagger (\Gamma_0 + \varepsilon(h_W)) \psi$, waarbij $h_W$ de single-particle Wilson-Dirac Hamiltonian is. De tekenfunctie $\varepsilon(X) = X/\sqrt{X^\dagger X}$ wordt benaderd door een polynoom $p_M(x)$ van graad $M$.
2. Block Encoding: De auteurs construeren een block encoding van de single-particle Wilson-Hamiltonian $\hat{h}_W$ (inclusief gaugevelden) met behulp van een "prepare" operator $P$ en een "select" operator $S$. Dit encodeert de Hamiltonian in een uitgebreide Hilbertruimte met ancilla-qubits via een unitaire operator $U_{\hat{h}_W}$.
3. Quantum Signal Processing (QSP): Om de benadering van de tekenfunctie te implementeren, gebruiken de auteurs QSP (specifiek Quantum Singular Value Transformation). Dit maakt het mogelijk om een graad-$M$ polynoom $p_M(U_{\hat{h}_W})$ toe te passen op de block-encoded operator. Het gekozen polynoom benadert de tekenfunctie binnen een spectrale kloof $[\kappa, 1]$ met een fout $\epsilon_e$.
4. Tijdsevolutie: Zodra de overlap-Hamiltonian is block-encoded, wordt de tijdsevolutie-operator $e^{-iHt}$ benaderd met behulp van QSP, waarbij een aantal poorten vereist is dat schaalt met de block-encoding kosten en de evolutietijd.

Belangrijkste Bijdragen en Resultaten

• Algoritmische Constructie: Het artikel levert een concreet QSP-gebaseerd algoritme voor de overlap-Hamiltonian dat de GW-relatie bewaart tot een controleerbare fout $\epsilon_e$. De gemodificeerde chirale operator $\hat{\gamma}_5$ wordt getoond de commutator te hebben met de benaderde Hamiltonian tot een fout van orde $\epsilon_e$.
• Complexiteitsanalyse:
  • Poortcomplexiteit: De kosten om de overlap-Hamiltonian te block-encoderen schalen als $O(Q \kappa^{-1} \log(1/\epsilon_e))$, waarbij $Q$ het aantal roosterplaatsen (vermenigvuldigd met interne vrijheidsgraden) is en $\kappa$ gerelateerd is aan de spectrale kloof. De totale poortcomplexiteit voor tijdsevolutie schaalt als $O(Q \kappa^{-1} \log(1/\epsilon_e) [Qt + \log(1/\epsilon_t)])$.
  • Qubit Kosten: Het algoritme vereist $O(Q + \log Q)$ qubits. Dit is een significante reductie vergeleken met domain-wall fermionen, die $O(Q L_5)$ qubits vereisen, waarbij $L_5$ de omvang van de extra dimensie is.
  • Vergelijking: Hoewel domain-wall fermionen profiteren van geometrische lokaliteit (waardoor een bijna lineaire kosten $O(Q t \text{polylog})$ via Suzuki-Trotter of soortgelijke methoden mogelijk is), brengt de overlap-formulering een hogere poortdiepte met zich mee door de niet-lokale aard van de tekenfunctie-benadering. Echter, de overlap-benadering biedt een superieure asymptotische schaling in geheugen.
• Fysieke Interpretatie van QSP: De auteurs demonstreren dat de vereiste polynoomgraad $M$ om een fout $\epsilon_e$ te bereiken schaalt als $M = O(\kappa^{-1} \log(1/\epsilon_e))$. Zij identificeren deze schaling met de omvang van de extra dimensie $L_5$ in domain-wall fermionen, waar de fout schaalt als $e^{-c L_5}$. De QSP-implementatie "construeert" dus effectief een extra dimensie door de diepte van het circuit, wat de fysieke correspondentie tussen de overlap-operator en de boundary-theorie van domain-wall fermionen weerspiegelt.

Betekenis en Claims
Het artikel claimt dat kwantumsimulaties van Dirac-fermionen met exacte chirale symmetrie "bijna gratis" zijn in de zin dat het toepassen van de overlap-Hamiltonian slechts een logaritmische factor meer kost (in termen van fouttolerantie) dan de Wilson-Dirac Hamiltonian.

Het centrale inzicht is dat QSP een kwantum-algoritmische realisatie biedt van de bekende correspondentie tussen de overlap-operator en de extra dimensie van domain-wall fermionen. De trade-off is expliciet:

• Domain-wall fermionen: Hogere qubit-kosten ($O(Q L_5)$) maar lagere poortdiepte door lokaliteit.
• Overlap-fermionen (via QSP): Lagere qubit-kosten ($O(Q)$) maar hogere poortdiepte ($O(Q^2)$ worst-case schaling voor tijdsevolutie door niet-lokaliteit).

De auteurs concluderen dat de keuze tussen formuleringen afhangt van hardwarebeperkingen (aantal qubits versus poortdiepte). Zij merken op dat hoewel klassieke rooster-QCD vaak rationale benaderingen (zoals Zolotarev) gebruikt voor de tekenfunctie, welke efficiënter zijn, deze geen directe analogen hebben in standaard QSP, wat wijst op een richting voor toekomstig werk. Het artikel verduidelijkt ook dat de schalingresultaten robuust zijn tegen de specifieke encoding van gaugevelden, hoewel constante factoren kunnen variëren.