Worldline-Susceptibility Scheduling for Quantum Annealing Beyond Local-Adiabatic Evolution

Technische Samenvatting: Worldline-Susceptibiliteit Planning voor Quantum Annealing

Probleemstelling
De prestaties van quantum annealing (QA) zijn kritiek afhankelijk van het annealing-schema, dat bepaalt hoe de totale looptijd wordt verdeeld over het evolutiepad. Hoewel het lokale adiabatische Roland–Cerf-schema theoretisch optimaal is—waarbij de tijd omgekeerd evenredig wordt toegewezen aan het kwadraat van de instantane spectrale kloof ($\Delta(s)^{-2}$)—is dit in de praktijk onhaalbaar voor grote optimalisatieproblemen. Het construeren van een dergelijk schema vereist volledige kennis van het instantane laag-energetische spectrum, wat herhaalde diagonalisatie van een exponentieel grote Hamiltoniaan vereist. Bijgevolg is er behoefte aan een computationeel goedkope surrogaat die de kritieke regio's van de anneal kan identificeren zonder expliciete spectrale informatie.

Methodologie
De auteurs stellen een surrogaat-planningframework voor op basis van de worldline magnetisatie-susceptibiliteit ($\chi_m$) gemeten tijdens Simulated Quantum Annealing (SQA). De methodologie verloopt als volgt:

1. SQA-simulatie: Het quantum annealing-proces wordt gemapt naar een equivalent klassiek statistisch-mechanisch systeem met behulp van de Suzuki–Trotter-decompositie. Dit transformeert het transversale veld Ising-model naar een $(d+1)$-dimensionaal klassiek systeem met $M$ gekoppelde replica's (Trotter-slices).
2. Observabele Meting: Evenwichtsconfiguraties worden gesampled met behulp van klassieke Monte Carlo-technieken. Hieruit wordt de worldline magnetisatie $m$ berekend, en de susceptibiliteit wordt afgeleid via de fluctuatie–dissipatie-stelling: $\chi_m = NM(\langle m^2 \rangle - \langle |m| \rangle^2)$.
3. Surrogaat Constructie: De gemeten susceptibiliteit $\chi_m(s)$ dient als een surrogaat voor de inverse-kwadraat spectrale kloof. Een tijdtoewijzingsgewicht $w(s) = \chi_m(s) + \chi_0$ (waarbij $\chi_0$ een regularisatieconstante is) wordt gedefinieerd. De cumulatieve tijdtoewijzingsfunctie $\tau(s)$ wordt geconstrueerd door deze gewichten te integreren, en het uiteindelijke schema $s(t)$ wordt verkregen door de inversie van $\tau(s)$.
4. Validatie: Het framework, geïmplementeerd in de open-source software Qanneal, valideert de surrogaat-schema's tegen exact diagonalisatie (ED) resultaten voor Sherrington–Kirkappert spin-glas instanties. Prestaties worden gekwantificeerd door de uiteindelijke grondtoestands-waarschijnlijkheid ($P_{GS}$) verkregen door het oplossen van de tijd-afhankelijke Schrödinger-vergelijking voor lineaire, exacte Roland–Cerf, en de voorgestelde surrogaat-schema's.

Belangrijkste Resultaten
De studie toont aan dat het worldline-susceptibiliteits-schema consequent beter presteert dan zowel conventionele lineaire annealing als, contra-intuïtief, het exacte Roland–Cerf-schema voor een aanzienlijk deel van de SK-instanties.

• Identificatie van Kritieke Regio's: De piek in de worldline-susceptibiliteit identificeert betrouwbaar de brede regio waar de spectrale kloof minimaal is, ook al reproduceert het niet de gedetailleerde gap-profiel kwantitatief.
• Superieure Eindtijd-Prestaties: Voor representatieve instanties ($n=10$ en $n=12$) levert het surrogaat-schema hogere grondtoestands-waarschijnlijkheden dan het exacte Roland–Cerf-schema over verschillende totale annealing-tijden.
• Ontdekking van Foutmodi: De auteurs identificeren twee onderscheidende eindtijd-foutmodi van het exacte lokale-adiabatische schema die verklaren waarom het onderpresteert:
  1. Boundary-Gap Trap: Wanneer de minimale kloof aan het einde van de anneal optreedt ($s \approx 1$), wijst het Roland–Cerf-schema de meerderheid van de looptijd toe aan een regio waar het transversale veld-driver effectief is verdwenen. Dit maakt de extra tijd dynamisch ineffectief, aangezien quantumfluctuaties daar onderdrukt zijn. Het susceptibiliteit-gebaseerde schema vermijdt dit door de tijd te concentreren in de bredere crossover-regio waar quantumfluctuaties nog actief zijn.
  2. Oscillerende Instabiliteit: Wanneer de minimale kloof intern ligt, concentreert het exacte schema de looptijd zo scherp rond de kloof dat het coherente multilevel-interferentie induceert. Dit resulteert in een niet-monotone, oscillerende afhankelijkheid van de succeswaarschijnlijkheid ten opzichte van de totale annealing-tijd, wat het vlakkere susceptibiliteit-gebaseerde schema voorkomt.
• Statistische Robuustheid: Grootschalige numerieke experimenten over verschillende wanorde-realisaties en systeemgroottes ($n=10$ tot $n=20$) bevestigen dat de superioriteit van het surrogaat-schema een generiek kenmerk is, en geen artefact van specifieke instanties. Naarmate de systeemgrootte toeneemt, neemt de prevalentie van boundary-gap instanties toe, wat de prestaties van het exacte Roland–Cerf-schema verder verslechtert ten opzichte van het surrogaat.

Betekenis en Claims
Het artikel betoogt dat onder realistische eindtijd-condities, robuustheid in planning belangrijker is dan de precieze reproductie van het instantane spectrale kloofprofiel. De auteurs claimen dat:

• Goedkope evenwichtsobservabelen (specifiek worldline-susceptibiliteit uit SQA) effectieve surrogaten kunnen dienen voor spectrale-gap-gebaseerde strategieën.
• De exacte lokale-adiabatische conditie, hoewel theoretisch optimaal in de asymptotische limiet, contraproductief kan zijn in eindtijd-regimes vanwege de "boundary-gap trap" en "oscillerende instabiliteit".
• De voorgestelde methodologie een schaalbaar pad biedt voor de optimalisatie van schema's die de computationele onhandelbaarheid van exacte spectrale berekeningen omzeilt, waardoor het toepasbaar is op probleemgroottes die buiten het bereik van exact diagonalisatie liggen.

Het werk is geïmplementeerd in het open-source Qanneal framework, wat een verenigd platform biedt voor SQA, susceptibiliteitsmeting en schema-constructie, waardoor reproduceerbaar onderzoek naar de dynamica van quantum annealing wordt gefaciliteerd.