Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

Technische Samenvatting: Het verkennen van de grondtoestand van het antiferromagnetische Heisenberg-model op een Kagome-rooster met behulp van een geometrisch geïnformeerde Variational Quantum Eigensolver

Probleemstelling
Het artikel behandelt de uitdaging van het vinden van de grondtoestand van het antiferromagnetische Heisenberg-model op een Kagome-rooster, een systeem dat bekend staat om geometrische frustratie en potentiële quantum spin-vloeistoffasen. Het bepalen van de grondtoestand van dergelijke $k$-lokale Hamiltonians is een QMA-hard probleem, wat klassieke simulatie voor grote systemen onuitvoerbaar maakt vanwege de exponentiële groei van de Hilbertruimte. Hoewel Variational Quantum Eigensolvers (VQE) een hybride quantum-klassieke benadering bieden die geschikt is voor Noisy Intermediate-Scale Quantum (NISQ) apparaten, worden ze geconfronteerd met aanzienlijke hindernissen: de noodzaak voor diepe circuits, de gevoeligheid voor ruis (coherente en incoherente fouten), barren plateaus waar gradiënten verdwijnen, en de moeilijkheid van het ontwerpen van ansatzen die zowel expressief als hardware-efficiënt zijn. Specifiek beogen de auteurs de eigenschappen van de grondtoestand te karakteriseren, zoals spin-spin correlaties en statische structuurfactoren, op echte quantumhardware zonder te vertrouwen op zware probleemspecifieke reducties of uitgebreide foutmitigatie tijdens de optimalisatielus.

Methodologie
De auteurs stellen een geometrisch geïnformeerde VQE-aanpak voor die specifiek is afgestemd op het Kagome-rooster, waarbij de focus ligt op twee fundamentele geometrieën: een enkele driehoek en een Kagome-ster (zes hoekgedeelde driehoeken).

1. Geometrisch Geïnformeerde Ansatz:
   De kerninnovatie is een aangepaste, ondiepe, hardware-efficiënte ansatz. In tegenstelling tot de Hamiltonian Variational Ansatz (HVA), die te diep kan zijn voor NISQ-apparaten, is deze ansatz geconstrueerd met een patroon van geparametriseerde single-qubit rotaties ($R_y$) en entangling gates die beperkt zijn tot de natuurlijke hardware-connectiviteit (het vermijden van SWAP-gates).

   • Euclidische Parameterruimte: Een cruciaal kenmerk van dit ontwerp is dat de Fubini-Study metriek-tensor analytisch bewezen diagonaal en constant is. Dit zorgt ervoor dat de parameterruimte van nature Euclidisch, singulariteitsvrij en vrij van overparameterisatie is. Bijgevolg is de quantum natural gradient equivalent aan de standaard gradiënt, wat de optimalisatielandschappen vereenvoudigt.
   • Hardware Compatibiliteit: De circuitdiepte wordt geminimaliseerd (één entangling laag) en het aantal parameters wordt verminderd door de plaatsing van $R_y$ gates af te wisselen.

2. Optimalisatiestrategie (AQNGD/AGD):
   De auteurs maken gebruik van een Adaptive Gradient Descent (AGD) optimizer. Omdat de ansatz een diagonale en constante Fubini-Study metriek oplevert, reduceert de Adaptive Quantum Natural Gradient Descent (AQNGD) tot AGD. De leersnelheid wordt dynamisch aangepast met behulp van de Armijo backtracking regel. Dit maakt grote initiële stappen mogelijk om convergentie te versnellen en verkleint automatisch de stapgrootte naarmate de optimizer de minimum nadert, waardoor stagnatie in nauwe geometrieën wordt voorkomen. Dit wordt vergeleken met de Simultaneous Perturbation Stochastic Approximation (SPSA).

3. Foutmitigatie (Post-verwerking):
   Foutmitigatie wordt toegepast na de VQE-optimalisatie om de uiteindelijke energie-inschattingen en structurele eigenschappen te verfijnen, in plaats van tijdens de iteratieve training.

   • Readout Error Mitigation (REM): Een gepartitioneerde responsmatrix-aanpak wordt gebruikt om SPAM-fouten (State Preparation and Measurement) te corrigeren, wat efficiënter schaalt dan volledige matrixinversie voor ijle metingen.
   • Zero-Noise Extrapolation (ZNE): De auteurs passen ZNE toe met kwadratische extrapolatie (en Bayesiaanse polynoomregressie) om resultaten te extrapoleren naar de zero-noise limiet. Ze merken op dat hoewel ZNE de Rayleigh-Ritz variational principe niet strikt naleeft, het hier wordt gebruikt om de nauwkeurigheid te verbeteren.

Belangrijkste Resultaten
De studie werd uitgevoerd op echte quantumhardware (IBM Oslo, IBM Kyoto, IBM Torino, en IBM Algiers) en gevalideerd via state-vector simulaties.

• Convergentie en Optimalisatie:

  • Op een 3-qubit driehoeksysteem convergeerde de AGD-geoptimaliseerde VQE soepel naar de grondtoestandsenergie ($E \approx -3$) met een hoge getrouwheid (99,89% op simulatie) met slechts 19 iteraties. De uiteindelijke parameters kwamen nauw overeen met de theoretische exacte waarden.
  • Op een 12-qubit Kagome-ster demonstreerde AGD een superieure convergentie vergeleken met SPSA, waarbij minder iteraties nodig waren om de grondtoestandsenergie ($E \approx -18$) te bereiken. De uiteindelijke getrouwheid op simulatie was 96,37%.
  • De stapgrootte-adaptatie via de Armijo-regel bleek effectief te zijn in het navigeren door het kostenlandschap, waarbij gestart werd met grote stappen en de verfijning toen de convergentie naderde.

• Ruisresistentie van Structurele Eigenschappen:

  • Een belangrijke bevinding is dat spin-spin correlaties en de statische structuurfactor ($S(q)$) resistent zijn tegen ruis. Zelfs zonder foutmitigatie behielden de experimentele $S(q)$ kaarten de kwalitatieve signatuur van de grondtoestand (bijv. de onderdrukking bij nul momentum die geen ferromagnetische orde aangeeft).
  • De Pearson correlatiecoëfficiënt tussen de ruizige experimentele $S(q)$ en de exacte theoretische $S(q)$ bleef extreem hoog (>99,9%).
  • Simulatie-experimenten onthulden een "stap-achtige" afhankelijkheid van de statische structuurfactor van de ruisniveaus: de correlatiestructuur blijft robuust tot een kritische ruisdrempel, waarna deze snel degradeert, wat lijkt op een faseovergang. Dit contrasteert met energie, die monotoon afneemt met ruis.

• Effectiviteit van Foutmitigatie:

  • Het toepassen van REM en ZNE na de optimalisatie verbeterde de energie-inschattingen aanzienlijk. Voor de Kagome-ster was de ruwe energie $\approx -14,9$, terwijl REM alleen de waarde verbeterde naar $\approx -17,3$, en ZNE (kwadratisch) de waarde dichter bij de exacte waarde van $-18$ bracht.
  • De auteurs observeerden dat het combineren van REM en ZNE soms leidde tot het overschrijden van de grondtoestandsenergie (het schenden van de variational bound), een gedrag dat consistent is met de niet-variationele aard van ZNE. ZNE alleen bood vaak de beste nauwkeurigheid onder de geteste configuraties.

Betekenis en Claims
Het artikel claimt dat een geometrisch geïnformeerde ansatz-ontwerp, dat een natuurlijk Euclidische parameterruimte garandeert, efficiënte training van VQE op huidige NISQ-hardware mogelijk maakt zonder de noodzaak van complexe, probleemspecifieke reducties.

• Robuustheid: Het werk demonstreert dat betekenisvolle fysieke eigenschappen, specifoor spin-spin correlaties en de statische structuurfactor, met voldoende nauwkeurigheid uit ruizige quantumhardware kunnen worden geëxtraheerd om de grondtoestand te karakteriseren, zelfs zonder foutmitigatie tijdens de optimalisatielus.
• Optimalisatie-efficiëntie: Het gebruik van AGD met adaptieve stapgroottes bleek efficiënter dan SPSA voor deze specifieke hardware-efficiënte ansatz, waarbij snellere convergentie in termen van iteraties werd bereikt.
• Schaalbaarheid: Door diepe circuits en zware foutmitigatie tijdens de training te vermijden, biedt de aanpak een schaalbaar pad voor het verkennen van gefrustreerde quantummagneten. De auteurs suggereren dat dit framework kan worden uitgebreid naar grotere Kagome-roosters en andere gefrustreerde quantumsystemen.
• Bescheidenheid: De auteurs erkennen dat hoewel hun methode succesvol kleine fundamentele cellen (driehoek en ster) kan karakteriseren, de algemene grondtoestand van het oneindige Kagome-rooster een open vraag blijft in de literatuur. Ze beweren niet de oplossing voor het oneindige roosterprobleem te hebben gevonden, maar hebben eerder een levensvatbare methodologie gedemonstreerd om deze systemen op huidige apparaten te verkennen. Ze merken ook op dat hoewel ZNE de energie-inschattingen verbetert, het de variational principe niet behoudt, en dat de combinatie van REM en ZNE zorgvuldige interpretatie vereist.