Quantum-inspired dynamical models on quantum and classical annealers

Titel: Quantum-geïnspireerde dynamische modellen op quantum- en klassieke annealers

Auteurs: Philipp Hanussek, Jakub Pawłowski, Zakaria Mzaouali en Bartłomiej Gardas.

---

1. Het Probleem

Het simuleren van de real-time evolutie van veel-deeltjes quantum-systemen is fundamenteel voor de gecondenseerde materie-fysica, quantumchemie en opkomende quantumtechnologieën. Klassieke simulaties worden echter onuitvoerbaar door de exponentiële groei van de Hilbert-ruimte, zelfs bij slechts enkele tientallen qubits.

Hoewel quantumhardware (zoals quantum annealers) beloofd wordt om dit probleem op te lossen, ontbreekt het aan een rigoureuze, eerlijke benchmarking. Bestaande studies kampen met twee hoofdbeperkingen:

1. Er is geen efficiënte, schaalbare methode om de Schrödinger-vergelijking (een reeks unitaire operatoren) om te zetten naar een statische optimalisatieprobleem die op annealers kan worden uitgevoerd.
2. Er ontbreken gestandaardiseerde benchmark suites om echte quantum-versnelling te onderscheiden van artefacten veroorzaakt door specifieke probleeminstellingen of toevallige inbeddingen op de hardware.

2. Methodologie

De auteurs introduceren een nieuw, physics-geïnspireerd benchmarking framework dat real-time quantumdynamica mapt naar een gemeenschappelijk optimalisatieformaat: Quadratic Unconstrained Binary Optimization (QUBO).

• Parallel-in-time encoding: De kern van de methode is het discretiseren van de tijdsinterval $[t_0, t_f]$ in $N$ stappen. De tijdsontwikkeling wordt gecodeerd in een "geschiedenis-toestand" (history state).
• Van Hamiltoniaan naar QUBO: De dynamica van een $n$-qubit systeem (mogelijk met niet-Hermitische Hamiltonianen) wordt omgezet in een kwadratische vorm $f(x) = \frac{1}{2}\langle x|A|x\rangle - \langle x|\phi\rangle$. Door variabelen in vast-punt binaire notatie te schrijven, ontstaat een QUBO-instantie met $O(LNR)$ spins.
• Solver-agnostisch: Deze QUBO-instanties kunnen zonder aanpassing worden uitgevoerd op zowel quantum annealers als klassieke optimalisatieoplossers.

Gebruikte Solvers:

1. Quantum Annealers: Twee generaties D-Wave-systemen:
   • Advantage (Pegasus-topologie).
   • Advantage2 (Zephyr-topologie, met 33% hogere connectiviteit).
   • Er werden zowel standaard voorwaartse annealing als cyclische annealing (een verfijningsstrategie) toegepast.
2. Klassieke Solvers:
   • VeloxQ: Een GPU-versnelde, physics-geïnspireerde heuristiek die QUBO's oplost zonder hardware-embeddingsbeperkingen.
   • Simulated Annealing (SA): Uitgevoerd op zowel CPU als GPU als referentiepunt.

Testgevallen:
Er werden acht representatieve dynamische modellen getest, variërend van enkel-qubit rotaties en multi-qubit verstrengeling (Bell, GHZ, cluster toestanden) tot PT-symmetrische en andere niet-Hermitische generatoren. Sommige problemen zijn "native" voor de hardware-topologie, andere vereisen "minor-embedding".

3. Belangrijkste Bijdragen

• Unificatie van Dynamica en Optimalisatie: Het artikel biedt een praktische, open-source suite om quantumdynamica direct te vertalen naar QUBO-problemen, waardoor een "like-for-like" vergelijking mogelijk is tussen fundamenteel verschillende rekenparadigma's.
• Uitgebreide Benchmark: Het is de eerste studie die systematisch prestaties vergelijkt tussen twee generaties D-Wave-processors en state-of-the-art klassieke heuristieken op een reeks fysiek gemotiveerde dynamische modellen.
• Schaalbaarheidstest: De auteurs testen het framework tot zeer grote instanties ($N \simeq 10^5$), waarmee een klassieke baseline wordt gevestigd voor toekomstige hardware.

4. Resultaten

De resultaten tonen een gemengd maar instructief beeld van de huidige staat van quantum versus klassiek rekenen:

• Hardware Vooruitgang (Advantage2 vs. Advantage):

  • De nieuwste Advantage2-processor presteert consistent beter dan zijn voorganger.
  • Voor native problemen is de succeskans (ground-state success probability) 4 tot 5 keer hoger.
  • Voor niet-native problemen (die embedding vereisen) is de verbetering een orde van grootte (tot 26 keer hoger).
  • Dit wijst op snelle hardwarevoortgang in qubit-coherentie, controleprecisie en connectiviteit.

• Klassieke Dominantie:

  • VeloxQ (de GPU-versnelde klassieke solver) behaalt de kortste absolute runtimes en behoudt de leiding op de momenteel toegankelijke probleemgroottes.
  • VeloxQ toont een schalingsexponent ($\beta$) van ongeveer 100, wat aanzienlijk beter is dan de quantum annealers (die rond $\beta \approx 9-33$ liggen, afhankelijk van de configuratie).
  • Klassieke solvers hebben geen last van de beperkte connectiviteit van de quantumhardware en kunnen problemen oplossen die voor de quantumchip onoplosbaar zijn (bijv. systeem H3 en H8).

• Cyclische Annealing:

  • Cyclische annealing verbetert de zoekopdracht lokaal en helpt om uit lokale minima te ontsnappen, maar levert geen fundamentele verandering in de exponentiële schaling van de workload op. Het is vooral nuttig als een gerichte verfijningsstap.

• Grootschalige Schaling ($N \approx 10^5$):

  • Bij zeer grote problemen (tot 80.000 variabelen) vertonen klassieke solvers een machts-wet schaling (power-law), waarbij VeloxQ Simulated Annealing consistent verslaat.
  • Dit bevestigt dat klassieke heuristieken momenteel nog superieur zijn voor het oplossen van deze specifieke dynamische problemen op grote schaal.

5. Betekenis en Conclusie

Dit werk positioneert het "parallel-in-time QUBO framework" als een veelzijdig en fysiek gemotiveerd testplatform om de voortgang richting quantum-competitieve simulatie van dynamische systemen kwantitatief te volgen.

• Conclusie: Hoewel D-Wave's nieuwste hardware aanzienlijke verbeteringen laat zien ten opzichte van eerdere generaties, overtreffen geoptimaliseerde klassieke algoritmen (vooral op GPU's) momenteel nog steeds de quantumannealers in snelheid en schaalbaarheid voor deze specifieke taken.
• Toekomstperspectief: De studie benadrukt dat een categorisch quantumvoordeel nog niet is bereikt. Echter, met de toenemende dichtheid van qubit-lattices en bredere koppelingstelsels, zouden medium-schaal dynamische simulaties (tientallen tot honderden effectieve vrijheidsgraden) binnenkort haalbaar kunnen worden binnen één annealing-run.
• Bijdrage aan de gemeenschap: De open-source benchmark suite biedt een transparante maatstaf waarmee toekomstige hardware- en algoritmische doorbraken kunnen worden gemeten.