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

Sintesi Tecnica: Sondare lo Stato Fondamentale del Modello di Heisenberg Antiferromagnetico su Reticolo Kagome utilizzando un Variational Quantum Eigensolver Geometricamente Informato

Definizione del Problema
Il documento affronta la sfida di trovare lo stato fondamentale del modello di Heisenberg antiferromagnetico su un reticolo kagome, un sistema noto per la frustrazione geometrica e i potenziali stati di liquido di spin quantistico. Determinare lo stato fondamentale di tali hamiltoniane $k$-locali è un problema QMA-hard, il che rende la simulazione classica intrattabile per sistemi di grandi dimensioni a causa della crescita esponenziale dello spazio di Hilbert. Sebbene i Variational Quantum Eigensolvers (VQE) offrano un approccio ibrido quantistico-classico adatto ai dispositivi Noisy Intermediate-Scale Quantum (NISQ), essi affrontano ostacoli significativi: la necessità di circuiti profondi, la suscettibilità al rumore (errori coerenti e incoerenti) e i plateau sterili (barren plateaus) dove i gradienti svaniscono, oltre alla difficoltà di progettare ansatz che siano sia espressivi che efficienti dal punto di vista dell'hardware. Nello specifico, gli autori mirano a caratterizzare le proprietà dello stato fondamentale, come le correlazioni spin-spin e i fattori di struttura statici, su hardware quantistico reale senza fare affidamento su pesanti riduzioni specifiche del problema o su un'estesa mitigazione dell'errore durante il ciclo di ottimizzazione.

Metodologia
Gli autori propongono un approccio VQE geometricamente informato adattato per il reticolo kagome, concentrandosi su due geometrie fondamentali: un singolo triangolo e una stella di kagome (sei triangoli che condividono gli angoli).

1. Ansatz Geometricamente Informato:
   L'innovazione principale è un ansatz personalizzato, poco profondo e hardware-efficient. A differenza dell'Hamiltonian Variational Ansatz (HVA), che può essere troppo profondo per i dispositivi NISQ, questo ansatz è costruito utilizzando un pattern di rotazioni a singolo qubit parametrizzate ($R_y$) e gate di entanglement limitati alla connettività nativa dell'hardware (evitando i gate SWAP).

   • Spazio dei Parametri Euclideo: Una caratteristica critica di questo design è che il tensore metrico di Fubini-Study è dimostrato analiticamente essere diagonale e costante. Ciò assicura che lo spazio dei parametri sia naturalmente euclideo, privo di singolarità e privo di sovra-parametrizzazione. Di conseguenza, il gradiente naturale quantistico è equivalente al gradiente standard, semplificando il panorama dell'ottimizzazione.
   • Compatibilità con l'Hardware: La profondità del circuito è minimizzata (un solo strato di entanglement) e il numero di parametri è ridotto alternando la disposizione dei gate $R_y$.

2. Strategia di Ottimizzazione (AQNGD/AGD):
   Gli autori impiegano un ottimizzatore Adaptive Gradient Descent (AGD). Poiché l'ansatz produce un metrico di Fubini-Study diagonale e costante, l'Adaptive Quantum Natural Gradient Descent (AQNGD) si riduce ad AGD. Il tasso di apprendimento viene adattato dinamicamente utilizzando la regola di backtracking di Armijo. Ciò consente passi iniziali ampi per accelerare la convergenza, riducendo automaticamente la dimensione del passo man mano che l'ottimizzatore si avvicina al minimo, evitando la stagnazione in geometrie strette. Questo viene confrontato con la Simultaneous Perturbation Stochastic Approximation (SPSA).

3. Mitigazione dell'Errore (Post-Processing):
   La mitigazione dell'errore viene applicata dopo l'ottimizzazione VQE per raffinare le stime finali dell'energia e delle proprietà strutturali, piuttosto che durante l'addestramento iterativo.

   • Mitigazione dell'Errore di Lettura (REM): Viene utilizzato un approccio a matrice di risposta partizionata per invertire gli errori SPAM (State Preparation and Measurement), scalando in modo più efficiente rispetto alla piena inversione della matrice per misurazioni sparse.
   • Zero-Noise Extrapolation (ZNE): Gli autori applicano ZNE utilizzando l'estrapolazione quadratica (e la regressione polinomiale bayesiana) per estrapolare i risultati al limite di rumore zero. Notano che, sebbene la ZNE non rispetti strettamente il principio variazionale di Rayleigh-Ritz, viene utilizzata qui per migliorare l'accuratezza.

Risultati Chiave
Lo studio è stato condotto su hardware quantistico reale (IBM Oslo, IBM Kyoto, IBM Torino e IBM Algiers) ed è stato validato tramite simulazioni di stato vettoriale.

• Convergenza e Ottimizzazione:

  • Su un sistema a triangolo a 3 qubit, il VQE ottimizzato tramite AGD è convergito fluidamente allo stato fondamentale dell'energia ($E \approx -3$) con alta fedeltità (99,89% nella simulazione) utilizzando solo 19 iterazioni. I parametri finali corrispondevano strettamente ai valori teorici esatti.
  • Su una stella di kagome a 12 qubit, l'AGD ha dimostrato una convergenza superiore rispetto a SPSA, richiedendo meno iterazioni per raggiungere l'energia dello stato fondamentale ($E \approx -18$). La fedeltà finale nella simulazione è stata del 96,37%.
  • Si è dimostrato che l'adattamento della dimensione del passo tramite la regola di Armijo permette di navigare efficacemente il panorama del costo, partendo da passi ampi e raffinandoli man mano che la convergenza si avvicina.

• Resilienza al Rumore delle Proprietà Strutturali:

  • Un risultato significativo è che le correlazioni spin-spin e il fattore di struttura statico ($S(q)$) sono resilienti al rumore. Anche senza mitigazione dell'errore, le mappe sperimentali di $S(q)$ hanno mantenuto la firma qualitativa dello stato fondamentale (ad esempio, la soppressione allo zero momentum che indica l'assenza di ordine ferromagnetico).
  • Il coefficiente di correlazione di Pearson tra il $S(q)$ sperimentale rumoroso e il $S(q)$ teorico esatto è rimasto estremamente alto (>99,9%).
  • Gli esperimenti di simulazione hanno rivelato una dipendenza "a gradini" del fattore di struttura statico rispetto ai livelli di rumore: la struttura di correlazione rimane robusta fino a una soglia critica di rumore, dopo la quale degrada rapidamente, somigliando a una transizione di fase. Questo contrasta con l'energia, che decade monotonicamente con il rumore.

• Efficacia della Mitigazione dell'Errore:

  • L'applicazione di REM e ZNE post-ottimizzazione ha migliorato significativamente le stime dell'energia. Per la stella di kagome, l'energia grezza era $\approx -14,9$, mentre la REM da sola l'ha migliorata a $\approx -17,3$, e la ZNE (quadratica) l'ha portata più vicina al valore esatto di $-18$.
  • Gli autori hanno osservato che combinare REM e ZNE a volte portava a superare l'energia dello stato fondamentale (violando il limite variazionale), un comportamento coerente con la natura non variazionale della ZNE. La ZNE da sola ha spesso fornito l'accuratezza migliore tra le configurazioni testate.

Significato e Rivendicazioni
Il documento afferma che un design di ansatz geometricamente informato, che garantisce uno spazio dei parametri naturalmente euclideo, consente un addestramento efficiente del VQE su attuali hardware NISQ senza la necessità di complesse riduzioni specifiche del problema.

• Robustezza: Il lavoro dimostra che proprietà fisiche significative, specificamente le correlazioni spin-spin e il fattore di struttura statico, possono essere estratte da hardware quantistico rumoroso con sufficiente accuratezza per caratterizzare lo stato fondamentale, anche senza mitigazione dell'errore durante il ciclo di ottimizzazione.
• Efficienza di Ottimizzazione: L'uso di AGD con dimensioni del passo adattive si è dimostrato più efficiente rispetto a SPSA per questo specifico ansatz hardware-efficient, ottenendo una convergenza più rapida in termini di iterazioni.
• Scalabilità: Evitando circuiti profondi e una pesante mitigazione dell'errore durante l'addestramento, l'approccio offre una via scalabile per sondare magneti quantistici frustrati. Gli autori suggeriscono che questo framework può essere esteso a reticoli kagome più grandi e altri sistemi quantistici frustrati.
• Modestia: Gli autori riconoscono che, sebbene il loro metodo caratterizzi con successo piccole celle fondamentali (triangolo e stella), lo stato fondamentale generale del reticolo kagome infinito rimane una questione aperta nella letteratura. Non pretendono di aver risolto il problema del reticolo infinito, ma hanno dimostrato una metodologia valida per sondare tali sistemi su dispositivi attuali. Notano inoltre che, sebbene la ZNE migliori le stime dell'energia, non preserva il principio variazionale, e che la combinazione di REM e ZNE richiede un'interpretazione attenta.