Benchmarking Open-Access Quantum Hardware Through Many-Body Dynamics in the 1D Tilted-Field Ising Model

Riepilogo Tecnico: Benchmarking di Hardware Quantistico Open-Access attraverso la Dinamica Many-Body nel Modello di Ising a Campo Inclinato 1D

Definizione del Problema
La simulazione di Hamiltoniane è una delle principali applicazioni proposte per i computer quantistici a breve termine, tuttavia la sua utilità pratica dipende dalla capacità dell'hardware fisico di preservare la dinamica many-body prevista tra i limiti di campionamento, il rumore dei gate, la decoerenza e i vincoli di profondità del circuito. Mentre il modello di Ising a campo trasversale è un sistema integrabile standard, l'aggiunta di un campo longitudinale crea il modello di Ising a campo inclinato, che rompe l'integrabilità e introduce dinamiche complesse e non banali, inclusi il caos quantistico e la termalizzazione. Il problema centrale affrontato è determinare l'estensione con cui l'attuale hardware quantistico ad accesso aperto possa riprodurre questi osservabili many-body fisicamente significativi rispetto a simulazioni ideali e prive di rumore.

Metodologia
Lo studio utilizza una catena di Ising a campo inclinato unidimensionale a 15 qubit come sistema di benchmark, eseguito sul backend IBM Quantum ibm_marrakesh.

• Inizializzazione del Sistema: Il sistema viene inizializzato con tutti i qubit nello stato spin-up ($|0\rangle$), eccetto per il qubit centrale ($q8$), che viene invertito in spin-down ($|1\rangle$) tramite un gate X per creare un'eccitazione localizzata.
• Hamiltoniana e Parametri: La dinamica è governata dall'Hamiltoniana $H = -J \sum Z_i Z_{i+1} - h_x \sum X_i - h_z \sum Z_i$. Lo studio confronta un caso di riferimento integrabile ($h_z = 0$) contro il caso a campo inclinato non integrabile ($h_z = h_x = 0.35625$) attraverso tre regimi di interazione: $J < h_x$, $J = h_x$, e $J > h_x$.
• Tecnica di Simulazione: L'evoluzione temporale è approssimata utilizzando una decomposizione di Trotter simmetrica del secondo ordine con un passo temporale $\Delta t = 0.1$. Questa scelta bilancia l'accuratezza (l'errore locale scala come $O(\Delta t^3)$) con la profondità del circuito per rimanere fattibile su hardware ad accesso aperto.
• Osservabili: Tre osservabili chiave vengono tracciate per valutare la dinamica:
  1. Magnetizzazione Locale: $\langle Z_i(t) \rangle$ per tracciare l'orientamento dello spin individuale.
  2. Magnetizzazione Totale: $M_z(t) = \sum \langle Z_i(t) \rangle$ per misurare la memoria della polarizzazione globale dello spin.
  3. Correlazioni di Spin Connesse: $C_{ij}^{conn}(t) = \langle Z_i Z_j \rangle - \langle Z_i \rangle \langle Z_j \rangle$ per analizzare la diffusione delle correlazioni e le strutture a cono di luce.
• Esecuzione Sperimentale: I circuiti sono stati eseguiti su ibm_marrakesh con 2048 shot per circuito. Lo studio ha impiegato il livello di resilienza 1 con deaccoppiamento dinamico (sequenza XpXm), ma non ha utilizzato l'estrapolazione del rumore zero implementata dall'utente o la piena correzione degli errori quantistici. I risultati dell'hardware sono stati campionati ogni 0.5 unità temporali fino a $t=35$.

Risultati Chiave
Lo studio confronta le simulazioni classiche Trotterizzate ideali con i risultati dell'esecuzione sull'hardware:

• Magnetizzazione Locale: L'hardware riproduce l'andamento generale del decadimento della magnetizzazione nei primi tempi. Tuttavia, i valori dell'hardware sono "compressi", mostrando una riduzione della separazione tra i qubit rispetto alla simulazione ideale. Nei casi a campo inclinato ($h_z = h_x$), la simulazione ideale mantiene valori di magnetizzazione locale più elevati e una struttura oscillatoria, mentre i dati dell'hardware mostrano un collasso più rapido verso lo zero e una ridotta variazione sito-sito.
• Magnetizzazione Totale: Similmente alla magnetizzazione locale, i dati dell'hardware mostrano un rapido calo iniziale seguito da uno stato stazionario inferiore rispetto alla simulazione ideale. Nel regime a campo inclinato, la simulazione ideale mantiene una significativa magnetizzazione positiva (memoria dello stato iniziale), mentre i dati dell'hardware mostrano una riduzione molto più forte, indicando una perdita della memoria della magnetizzazione globale.
• Correlazioni Connesse: La simulazione ideale esibisce una chiara struttura a "cono di luce" dove le correlazioni si diffondono verso l'esterno dall'eccitazione centrale ($q8$) in modo strutturato. Al contrario, i risultati dell'hardware mostrano un pattern significativamente più debole e rumoroso. Il fronte di diffusione distinto è in gran parte oscurato, sostituito da fluttuazioni sparse a bassa intensità, particolarmente ai tempi lunghi e alle distanze maggiori dal centro.

Contributi Chiave

• Protocollo di Benchmarking: Il documento stabilisce un protocolo di benchmarking riproducibile a livello di osservabile per la fedeltà della diffusione delle correlazioni in un sistema a catena di spin non integrabile su hardware IBM di classe Heron.
• Confronto Quantitativo: Fornisce un confronto diretto e quantitativo tra dinamiche ideali e rumorose per un sistema a 15 qubit, isolando specificamente gli effetti della rottura dell'integrabilità sulle prestazioni dell'hardware.
• Caratterizzazione del Rumore: Lo studio dimostra che, sebbene l'hardware ad accesso aperto possa catturare caratteristiche grossolane della dinamica nei tempi iniziali, la diffusione delle correlazioni a lungo termine e le strutture oscillatorie dettagliate sono sempre più oscurate dal campionamento finito, dagli errori dei gate a due qubit, dagli errori di lettura e dalla decoerenza.

Significato e Rivendicazioni
Il documento afferma modestamente che lo studio serve come un benchmark a breve termine piuttosto che come una simulazione completamente corretta dagli errori. Il suo significato primario risiede nel dimostrare i limiti attuali dell'hardware quantistico ad accesso aperto nel preservare un comportamento many-body fisicamente significativo.

• Preservazione della Dinamica: L'hardware preserva alcune caratteristiche dei tempi iniziali (tendenze generali della magnetizzazione) ma non riesce a preservare in modo affidabile la diffusione delle correlazioni a lungo termine.
• Limitazioni: I risultati evidenziano che, senza mitigazione o correzione degli errori, la struttura dettagliata delle dinamiche non integrabili (come specifici pattern di correlazione o memoria oscillatoria) è difficile da estrarre dai dati rumorosi.
• Direzioni Future: Gli autori suggeriscono che, sebbene la mitigazione degli errori (ad esempio, mitigazione dell'errore di lettura, estrapolazione del rumore zero) potrebbe migliorare l'accuratezza, la piena correzione degli errori quantistici rimane un overhead significativo per i dispositivi a breve termine. Propongono come lavoro futuro confronti tra diversi backend, catene di spin più piccole per isolare il rumore e l'inclusione di osservabili più complessi come gli OTOC (out-of-time-order correlators) e l'entropia di entanglement per sondare ulteriormente il caos quantistico e la diffusione dell'informazione.