Reshaping quantum annealing landscapes with diagonal catalysts

Technische Samenvatting: Het hervormen van Quantum Annealing Landschappen met Diagonale Katalysatoren

Probleemstelling
Quantum annealing (QA) is een optimalisatiealgoritme dat ontworpen is om de grondtoestand van een kosten-Hamiltoniaan ($H_P$) te vinden die een Quadratic Unconstrained Binary Optimization (QUBO) probleem representeert. De prestaties van QA worden echter vaak beperkt doordat de populatie vast komt te zitten in lokale minima die ver verwijderd zijn van de globale oplossing in termen van Hamming-afstand. Dit probleem wordt verergerd op moeilijke instanties waarbij de energiekloof tussen de grondtoestand en dynamisch toegankelijke geëxciteerde toestanden exponentieel klein wordt. Bij gevolg is de populatie tijdens eindige-tijd anneals vaak diabatisch, wat leidt tot populatieherverdeling over vermeden kruisingen (avoided crossings) tussen toestanden die gelokaliseerd zijn in verschillende bassins van het energie-landschap. Een kritieke uitdaging is dat configuraties met lage energie niet noodzakelijkerwijs overeenkomen met configuraties die dicht bij de oplossing liggen in Hamming-afstand; verre configuraties kunnen nabije configuraties "ondergraven" door lagere energieën te bieden als kandidaten voor concurrerende bassins.

Methodologie en Theoretisch Kader
De auteurs introduceren een wiskundig kader dat energie en Hamming-afstand koppelt om twee-body diagonale katalysatoren te ontwerpen. Deze katalysatoren worden uitsluitend geconstrueerd op basis van de koppelingsstructuur van het probleem (tekens, grootte en graaf-topologie), zonder voorafgaande kennis van de oplossing te vereisen.

1. Hamming-shell Analyse: De auteurs definiëren "Hamming-shells" ($S_{d, z^\star}$) als verzamelingen van configuraties op een specifieke Hamming-afstand $d$ van een referentie-grondtoestand $z^\star$. Zij leiden exacte formules af voor het gemiddelde en de variantie van de energie over deze shells voor elke veldvrije twee-body Ising-operator.

   • Stelling 1 stelt vast dat het gemiddelde energie van een shell een parabolische curve volgt die bepaald wordt door combinatorische factoren, terwijl de variantie (shell-breedte) afhankelijk is van de specifieke koppelingsmatrix.
   • In het kale probleem zijn de gemiddelde energieën van de shells correct geordend, maar de brede variantie veroorzaakt overlap tussen de shells, waardoor verre configuraties lagere energieën kunnen hebben dan nabije configuraties.

2. Katalysatorontwerp ($C$): Het doel is om een diagonale katalysator $C$ te construeren die de scheiding tussen de gemiddelde energie van de shells vergroot ten opzichte van hun fluctuaties binnen de shells, wat effectief een "trechter" (funnel) in het landschap creëert.

   • De ideale katalysator zou proportioneel zijn aan $-z^\star_i z^\star_j$, maar dit vereist kennis van de onbekende grondtoestand.
   • Open-Pad Constructie: De auteurs stellen voor om het grondtoestandpatroon te schatten met behulp van zelf-vermijdende paden (self-avoiding paths) op de interactie-graaf. Door spin-patronen langs deze paden te propageren op basis van koppelings-tekens ($g(p)_{v_{a+1}} = -\text{sgn}(J_{v_a v_{a+1}}) g(p)_{v_a}$), genereren zij lokaal bevredigende patronen.
   • De katalysator $C^{(m)}$ wordt gevormd door de bijdragen van alle $m$-rand zelf-vermijdende paden op te tellen. Elk pad draagt een gewogen rank-one schatting bij van de ideale koppelingsstructuur. De uiteindelijke operator wordt loodrecht op de oorspronkelijke koppelingen geprojecteerd en genormaliseerd.
   • Deze constructie creëert een dichtere interactie-graaf dan het oorspronkelijke probleem, waardoor configuraties die sterke bindingen frustreren worden afgestraft en configuraties die met de gepropageerde patronen in lijn zijn, worden versterkt.

Belangrijkste Resultaten
De auteurs hebben hun methode getest op 200 willekeurige 3-reguliere instanties met $n=20$ spins en koppelingen getrokken uit $[-1, 1]$. Simulaties werden uitgevoerd bij sweep-tijden $T = 2, 6, 10$.

• Hervorming van het Energielandschap: De katalysator slaagt erin het statische landschap te hervormen, waardoor de overlap tussen de Hamming-shells wordt verminderd en de energie-trechter naar de grondtoestand wordt verdiept.
• Waarschijnlijkheidsherverdeling:
  • Energieruimte: De katalysator verschuift de waarschijnlijkheidsmassa van tussenliggende energiebanden naar de grondtoestand-band ($E/E_{GS} \ge 0.95$). Bij $T=6$ nam de mediaan van de massa in de bovenste band toe van 0.067 (niet-gekatalyseerd) naar 0.324, wat een mediane gepaard winst vertegenwoordigt van +351%.
  • Hamming-ruimte: De distributie verschuift naar de grondtoestand-manifold. Bij $T=6$ nam de waarschijnlijkheidsmassa bij Hamming-afstand $\delta \le 2$ toe van 0.10 naar 0.35 (+239%).
• Robuustheid en Vergelijking:
  • De winsten blijven bestaan over verschillende graaf-connectiviteiten (3-, 4-, en 5-reguliere grafen) en op het volledig verbonden Sherrington–Kirkpatrick model, hoewel de mate van verbetering afneemt naarmate de connectiviteit toeneemt.
  • De open-pad katalysator presteert aanzienlijk beter dan een "versterkt-probleem" baseline ($C = H_P$), wat aantoont dat het voordeel voortkomt uit de specifieke structurele herverdeling van koppelingen en niet enkel uit het vergroten van de algehele energieschaal.
  • Er is een afruil tussen de orde van het pad ($m$) en de graaf-dichtheid; matige ordes (bijv. $m=4$ voor 3-reguliere grafen) bieden de beste balans tussen het oplossen van de lokale structuur en het behouden van de getrouwheid aan de ware grondtoestand-volgorde.

Betekenis en Claims
Het artikel beweert een praktische methode te bieden voor het verbeteren van quantum annealing op nabije hardware door gebruik te maken van natieve longitudinale controles (diagonale operatoren) die twee-lokale operatoren zijn en geen complexe gadgets vereisen voor niet-stoquastische katalysatoren.

• Oplossings-agnostisch Ontwerp: De katalysatoren worden geconstrueerd zonder kennis van de oplossing, waarbij enkel vertrouwd wordt op de koppelingsmatrix van het probleem.
• Statische Criteria voor Ontwerp: Het afgeleide shell-moment thema biedt statische criteria (trechterdiepte en shell-overlap) om kandidaat-katalysatoren te screenen voordat er dynamische annealing-simulaties worden uitgevoerd.
• Hardware Compatibiliteit: Omdat ze diagonaal en twee-body zijn, behouden de katalysatoren de transversale-veld structuur en zijn ze compatibel met de connectiviteits- en koppelingsbereik-beperkingen van huidige annealers.
• Toekomstige Richtingen: De auteurs suggereren dat het kader kan worden uitgebreid naar digitale instellingen (bijv. als een fase-separator in QAOA), gecombineerd kan worden met geoptimaliseerde diabatische protocollen, of toegepast kan worden op reverse annealing. Zij merken ook op dat het karakteriseren van de schaling met systeemgrootte en het generaliseren naar veld-inclusieve kostenfuncties een openstaande taak blijft.