Pauli-Sparse regularised Counterdiabatic Shortcuts for Linear-Ramp QAOA

Technische Samenvatting: Pauli-ijle Geregulariseerde Counterdiabatische Afkortingen voor Lineaire-Ramp QAOA

Probleemstelling
Combinatorische optimalisatie via het Quantum Approximate Optimization Algorithm (QAOA) staat voor een kritieke afweging tussen parameterefficiëntie en dynamische nauwkeurigheid. Lineaire-ramp QAOA (LR-QAOA) vermindert de klassieke optimalisatiebelasting door de lagen-hoeken te beperken tot een adiabatisch geïnspireerd schema met weinig parameters. Echter, voor Hamiltonians met kleine of exponentieel kleine spectrale kloven, lijdt de eindige-tijd adiabatische evolutie onder sterke diabatische transities, wat ervoor zorgt dat het systeem uit het gewenste laag-energetische oplossingsmanifold lekt. Hoewel counterdiabatische (CD) sturing een principieel mechanisme biedt om deze transities te onderdrukken door een correctieterm toe te voegen die is afgeleid van de Adiabatic Gauge Potential (AGP), vereisen standaardbenaderingen vaak dichte operatoren of vooraf gespecificeerde ansatzes die onhaalbaar zijn onder eindige gate-budgetten. Bovendien vereist het berekenen van exacte AGP's voor harde optimalisatieproblemen vaak het oplossen van microscopisch kleine, exponentieel kleine energie-splitsingen binnen een laag-energetisch manifold, wat computationeel duur is en potentieel onnodig.

Methodologie
De auteurs stellen een raamwerk voor het construeren van een Pauli-ijle, geregulariseerde counterdiabatische uitbreiding van LR-QAOA. De kernmethodologie bestaat uit het direct oplossen van de geregulariseerde AGP-vergelijking in de Pauli-basis zonder de volledige dichte Liouvilliaanse matrix te vormen.

1. Geregulariseerde AGP-vergelijking: De methode lost het lineaire systeem $(\mathcal{L}_H^2 + \eta I) A_\lambda^{(\eta)} = -i \mathcal{L}_H(G)$ op, waarbij $\mathcal{L}_H(X) = [H, X]$ de commutator-superoperator is, en $G = \partial_\lambda H$, met $\eta > 0$ als een regularisatieparameter. Deze vergelijking is afgeleid van het minimaliseren van een geregulariseerd functionaal dat de annulering van diabatische transities balanceert tegen de norm van de potentiaal.
2. Energie-resolutie Mechanisme: De parameter $\eta$ fungeert als een energie-resolutieschaal. Het onderdrukt transities geassocieerd met energieverschillen $|\omega| \ll \sqrt{\eta}$ terwijl het transities met $|\omega| \gg \sqrt{\eta}$ behoudt. Dit stelt de methode in staat om exponentieel kleine splitsingen binnen een doelmatig laag-energetisch manifold te negeren (het vermijden van de noodzaak voor exacte adiabatische gevolg), terwijl het wel lekken naar hogere-energie excitatietoestanden voorkomt.
3. Inexacte Conjugate Gradient (CG) Solver: In plaats van een dichte AGP te berekenen en te comprimeren, gebruikt de auteur een inexacte conjugate-gradient methode in Pauli-coördinaten.
   • Matrix-vrije Toepassing: De operator $(\mathcal{L}_H^2 + \eta I)$ wordt symbolisch toegepast met behulp van Pauli-commutatoren, waardoor de constructie van de $4^n$-dimensionale matrix wordt vermeden.
   • Truncatie: Tijdens de CG-iteraties wordt de Pauli-expansie afgekapt om alleen de $m$ grootste coëfficiënten te behouden (magnitude-truncatie). Dit transformeert de lineaire solver in een support-discovery procedure, die een gate-budget-bewuste set implementeerbare Pauli-rotaties identificeert.
   • Residu-herberekening: Om foutaccumulatie door truncatie te voorkomen, wordt het residu in elke stap exact herberekend in plaats van recursief bijgewerkt.
4. Galerkin Refit en Certificering: Zodra een ijle support is ontdekt, wordt een Galerkin refit uitgevoerd om de coëfficiënten op deze specifieke support te optimaliseren. De kwaliteit van de benadering wordt a posteriori gecertificeerd met een residu-grens: $\|A - A_\lambda^{(\eta)}\|_{HS} \leq \|R\|_{HS} / \eta$, waarbij $R$ het residu van het lineaire systeem is. Dit biedt een berekenbare foutgrens zonder kennis van de exacte oplossing te vereisen.
5. Circuit Integratie: De resulterende ijle Pauli-generator wordt in de LR-QAOA-circuit ingevoegd als een counterdiabatische laag. De variationele optimalisatie wordt uitgevoerd over een gereduceerde set parameters: de ramp-amplitudes $(\Delta\beta, \Delta\gamma)$ en de counterdiabatische amplitude $\Delta_{CD}$.

Kernbijdragen

• Directe Oplossing in Pauli-coördinaten: In tegenstelling tot eerdere CD-QAOA werken die ansatzes voorschrijven (bijv. lokale termen, geneste commutatoren) of dichte AGP's berekenen, lost dit werk de geregulariseerde AGP-vergelijking direct in Pauli-coördinaten op, waarbij simultaan de potentiaal wordt benaderd en een ijle set implementeerbare rotaties wordt geselecteerd.
• Regularisatie als Fysische Filter: Het artikel herkadert de regularisatieparameter $\eta$ niet enkel als een numerieke stabilisator, maar als een fysische filter die transities onder een specifieke energieschaal onderdrukt. Dit is bijzonder effectief voor optimalisatieproblemen waar het doel is om binnen een laag-energetisch manifold te blijven in plaats van een enkele eigentoestand te volgen door exponentieel kleine kloven.
• Sparse Lineaire Algebra Pipeline: De auteurs ontwikkelen een volledige pipeline die matrix-vrije symbolische commutatoren, getrunceerde CG voor support-discovery, Galerkin refitting voor coëfficiëntoptimalisatie en residu-gebaseerde foutcertificering combineert. Deze aanpak importeert ideeën uit de inexacte Krylov-methoden naar de context van kwantumcontrole.
• Theoretische Garanties: Het artikel biedt strikte grenzen aan de afstand tussen de geregulariseerde en ongeregulariseerde potentialen, de conditiegetal van het lineaire systeem, en de fout van de ijle benadering ten opzichte van de exacte geregulariseerde oplossing.

Resultaten
Numerieke experimenten werden uitgevoerd op Ferromagnetische Keten (FC) instanties en geperturbeerde FC–MaxCut/MarketSplit instanties:

• FC Validatie: Op FC-instanties vertoonde standaard LR-QAOA zwakke prestaties naarmate de systeemgrootte toenam, waarbij het faalde om het geordende manifold te bereiken. De voorgestelde LR-CD-QAOA met exacte AGP bereikte hoge geordende-manifold waarschijnlijkheden ($P_{ord} \approx 0.9$) zelfs bij lage diepte. Cruciaal is dat de methode robuust bleef over een reeks $\eta$-waarden, wat bevestigt dat het niet de exponentieel kleine grondtoestand-splitting hoeft te resolveren om lekken naar geëxciteerde toestanden te onderdrukken.
• Sparse CG Convergentie: Experimenten toonden aan dat een klein aantal CG-iteraties ($K \approx 4$) en een matig Pauli-support-budget ($m=64$) voldoende waren om de meeste voordelen van de dichte AGP te herstellen. De benaderingsratio verbeterde monotoon met zowel het aantal iteraties als de support-grootte.
• Blended Instanties: Op geperturbeerde MaxCut en MarketSplit instanties (waar de Hamiltonian wordt gedomineerd door de FC-component), presteerde LR-CD-QAOA significant beter dan de ongecorrigeerde LR-QAOA. Terwijl LR-QAOA zelfs bij diepte $p=8$ ver van het optimum bleef, bereikte de CD-gecorrigeerde versie benaderingsratio's die dicht bij het exacte optimum liggen bij matige dieptes.

Betekenis
Het artikel claimt dat dit raamwerk de praktische toepasbaarheid van QAOA op QUBO-optimalisatie substantieel vergroot. Door de robuustheid te verbeteren over heterogene probleemlandschappen — specifiek die met bijna-gedegenereerde laag-energetische structuren en kleine spectrale kloven — pakt de methode een primaire bottleneck in nabije-term kwantumoptimalisatie aan. Het biedt een pad voor het implementeren van effectieve counterdiabatische sturing op kwantumhardware met beperkte gate-budgetten, waarmee de beperkingen van parameterefficiënte maar dynamisch inefficiënte lineaire ramps worden overstegen. Het werk vestigt een hybride optimalisatieparadigma waarbij klassieke ijle lineaire algebra-technieken worden gebruikt om kwantumcircuits te construeren die dynamisch zijn aangepast aan de specifieke spectrale eigenschappen van de probleem-Hamiltonian.