Feynman's clock and hierarchy-informed sampling for quantum error mitigation

Technische Samenvatting: Feynman's Klok en Hiërarchie-geïnformeerde Sampling voor Quantumfoutmitigatie

Probleemstelling
Nabije quantumalgoritmen worden geconfronteerd met aanzienlijke beperkingen door de ruisgevoelige aard van de huidige NISQ-hardware (Noisy Intermediate-Scale Quantum), die de fysieke uitvoering van quantumcircuits verstoort. Hoewel foutcorrectie vereist dat grootschalige apparaten beschikbaar zijn die momenteel nog niet aanwezig zijn, zijn efficiënte technieken voor quantumfoutmitigatie (QEM) noodzakelijk voor toepassingen in de nabije toekomst. Bestaande methoden zijn fundamenteel beperkt, maar blijven de primaire methode om quantumruis te verminderen. Dit artikel behandelt de behoefte aan een algemene QEM-techniek die in staat is om ruis te mitigeren die voortkomt uit de uitvoering van willekeurige quantumcircuits, waarmee de toepassing wordt uitgebreid voorbij eerdere methoden die beperkt waren tot specifieke spinketen-simulaties.

Methodologie
De voorgestelde techniek, een uitbreiding van het BBGKY-geïnformeerde sampling-mitigatieschema (BBGKY-ISM), werkt via een proces van drie fasen:

1. Circuit-naar-Hamiltoniaan Mapping via Feynman's Klok:
   De auteurs mappen een willekeurig $N_Q$-qubit circuit $U$, bestaande uit $N_G$ poorten, op de tijdsevolutie van een quantummechanisch systeem dat wordt beheerst door een tijdonafhankelijke Hamiltoniaan. Hiervoor wordt gebruikgemaakt van de Feynman quantumcomputer-constructie, waarbij een klokregister $H_C$ met $N_C = \lceil \log_2(N_G + 1) \rceil$ ancilla-qubits wordt geïntroduceerd. Het systeem evolueert onder de Feynman-klok-Hamiltoniaan $H$, waarbij de kloktoestanden de voortgang van de circuitberekening volgen. De volledige dynamica wordt beschreven door $|\psi(t)\rangle = e^{-itH}|\psi(0)\rangle$.

2. Reconstructie via Hadamard-testen:
   In plaats van de tijdsevolutie te simuleren via Trotterisatie (wat de ruis zou versterken), reconstrueert de methode de tijdsafhankelijke verwachtingswaarden van belang via Hadamard-testen op gedeeltelijke uitvoeringen van het oorspronkelijke circuit.

• De Z-metingen op het data-register worden geherinterpreteerd als tijdsafhankelijke verwachtingswaarden van enkelvoudige Pauli-Z operatoren, geprojecteerd op de finale kloktoestand $|N_G\rangle_C$.
• Deze grootheden worden gedecomposeerd in een set "grootheden van belang" die Pauli-strings omvatten.
• De reële en imaginaire delen van de noodzakelijke matrixelementen worden gemeten met behulp van Hadamard-testen, waardoor de diepe circuitdiepte die bij Trotter-segmenten hoort, wordt vermeden.
• De uiteindelijke Z-meetwaarden worden geëxtraheerd door de gereconstrueerde tijdsevolutie-grootheden op te tellen, gewogen door de waarschijnlijkheidsamplitude van de kloktoestand, wat de tijdsafhankelijkheid effectief verwijdert.

3. BBGKY-Geïnformeerde Sampling Mitigatie:
   De ruisgevoelige schattingen verkregen uit de Hadamard-testen worden gemitigeerd met behulp van het BBGKY-ISM-schema.

• Fysica-geïnformeerde Verdeling: De methode bemonster (samples) mitigatiekandidaten uit een waarschijnlijkheidsverdeling die gecodeerd is met de ruisvrije dynamica. Deze verdeling wordt gedefinieerd door een actie $S(\vec{x})$ die wordt geminimaliseerd wanneer de configuratie $\vec{x}$ de ideale, ruisvrije evolutie benadert.
• Hiërarchische Constraints: De constraints worden afgeleid van een BBGKY-achtige hiërarchie van bewegingsvergelijkingen voor Pauli-string correlatoren. De hiërarchie verbindt de grootheden van belang met hogere-orde correlatoren.
• Simulated Annealing: Een Markov Chain Monte Carlo-procedure (Simulated Annealing) wordt gebruikt om configuraties te bemonsteren die de actie minimaliseren. Een parameter $r$ (straal) bepaalt hoeveel lagen van de hiërarchie in de constraints worden opgenomen. Een grotere $r$ codeert een groter deel van de ideale dynamica, wat de constraints op de sampling aanscherpt.

Belangrijkste Bijdragen

• Generalisatie van BBGKY-ISM: Het artikel breidt het BBGKY-ISM-schema uit van specifieke spinketen-simulaties naar willekeurige quantumcircuits door gebruik te maken van de Feynman-klokconstructie.
• Efficiënte Schaling van de Overhead: De auteurs demonstreren dat zowel de klassieke als de quantum overhead van de voorgestelde methode polynomiaal schalen met de omvang van het circuit ($N_G$), het aantal qubits ($N_Q$) en het aantal tijdstappen ($N_T$). Specifiek schaalt het aantal vereiste Hadamard-testen als $N_Q \cdot \text{poly}(N_G)$.
• Post-processing Karakter: De techniek is volledig een post-processing methode, waardoor deze gekoppeld kan worden aan bestaande mitigatieschema's of toegepast kan worden na de data-acquisitie zonder de uitvoering op de quantumhardware te wijzigen.
• Controleerbare Foutreductie: De methode introduceert een instelbare parameter $r$ (de hiërarchie-straal) die gebruikers in staat stelt de intensiteit van de mitigatie te controleren.

Resultaten
De methode werd numeriek getest op instelbare Bell-toestand preparatie-circuits ($N_Q=2, N_G=2$) onder een ruismodel dat de IBM Fez-processor nabootst.

• Systematische Reductie: De resultaten tonen een systematische reductie van de quantumfout ($\Delta' < \Delta$) over alle geteste circuitparameters en hiërarchie-radii.
• Instelbaarheid: Het vermogen tot foutreductie neemt toe naarmate de hiërarchie-straal $r$ groter wordt. Dit bevestigt dat het coderen van een groter deel van de BBGKY-hiërarchie de sampling effectiever dwingt richting de ruisvrije dynamica.
• Zelfconsistentie: Bij een specifieke straal ($r=4$ in de geteste casus) wordt het stelsel van vergelijkingen volledig zelfconsistent ($z=1$). Op dat punt reduceert de methode tot een puur klassieke numerieke berekening die de exacte ruisvrije dynamica herstelt, waarbij de bijdrage van de quantumruis effectief wordt genegeerd.

Betekenis en Claims
Het artikel claimt een nieuwe, algemene quantumfoutmitigatietechniek te bieden die zowel efficiënt (polynomiale schaling) als controleerbaar is. Door circuits te mappen naar Hamiltoniaan-dynamica en gebruik te maken van physics-informed sampling, biedt de methode een systematische manier om ruis te verminderen zonder extra quantumbronnen nodig te hebben buiten de initiële circuituitvoering en standaard post-processing. De auteurs benadrukken dat de effectiviteit van de techniek instelbaar is via de hiërarchie-straal, wat het aanpasbaar maakt aan verschillende ruisniveaus en computationele beperkingen. Het werk suggereert potentiële toekomstige toepassingen voor complexere algoritmen (bijv. Grover's, Shor's) en diepere variatie-ansatze, waarbij wordt opgemerkt dat het compatibel is als een post-processing laag met andere mitigatiestrategieën.