A near-term quantum simulation of the transverse field Ising model hints at Glassy Dynamics

Sintesi Tecnica: Una Simulazione Quantistica a Breve Termine del Modello di Ising a Campo Trasverso Suggerisce una Dinamica Vetrosa

Definizione del Problema
Il documento affronta la sfida computazionale di simulare sistemi di spin glass su larga scala, che sono critici per comprendere le fasi esotiche della materia e le transizioni vetrose nella fisica molti-corpo quantistica. Gli approcci computazionali classici, come la dinamica molecolare e le simulazioni Monte Carlo, faticano a causa della crescita esponenziale delle funzioni d'onda elettroniche nei sistemi multi-elettronici e dell'elevata complessità computazionale richiesta per modellare sistemi disordinati su tempi di rilassamento estesi. Nello specifico, i sistemi vetrosi esibiscono paesaggi energetici accidentati con numerosi minimi locali ed tempi di rilassamento estremamente lunghi, rendendoli difficili da simulare accuratamente utilizzando l'informatica classica ad alte prestazioni (HPC). Gli autori mirano a esplorare se le tecniche di simulazione quantistica digitale (DQS), specificamente gli algoritmi Variational Quantum Eigensolver (VQE), possano indagare efficacemente queste configurazioni di spin disordinate sull'attuale hardware quantistico a breve termine.

Metodologia
Gli autori hanno impiegato un approccio ibrido quantistico-classico utilizzando l'algoritmo Variational Quantum Eigensolver (VQE) per simulare il Modello di Ising a Campo Trasverso (TFIM) con l'aggiunta di un campo longitudinale. Lo studio si è concentrato su due geometrie specifiche:

1. Catena di Ising 1D: Una catena di lunghezza $n=25$ (25 qubit).
2. Reticolo di Ising 2D: Un reticolo quadrato $6 \times 6$ (36 qubit).

Componenti Tecniche Chiave:

• Costruzione dell'Hamiltoniana: L'Hamiltoniana del sistema include interazioni tra vicini più prossimi lungo l'asse z ($J \sum \sigma_i^z \sigma_{i+1}^z$), un campo trasverso nella direzione x ($-B_x \sum \sigma_i^x$) e un campo longitudinale nella direzione z ($-B_z \sum \sigma_i^z$). Le interazioni diagonali sono state escluse.
• Selezione dell'Ansatz: Gli autori hanno utilizzato un ansatz "TwoLocal" composto da blocchi di rotazione $R_x$ e $R_y$. Fondamentalmente, questo ansatz è stato progettato senza entanglement tra i qubit. Questa scelta è stata fatta per ridurre la profondità del circuito e mitigare la suscettibilità al rumore sui dispositivi a breve termine, pur catturando le fluttuazioni quantistiche necessarie per la struttura del problema.
• Ottimizzazione: L'ottimizzatore classico COBYLA è stato utilizzato per minimizzare il valore di aspettazione dell'Hamiltoniana, con un massimo di 25.000 iterazioni (sebbene la maggior parte delle esecuzioni sia confluita prima).
• Implementazione Software: Gli autori hanno sviluppato un pacchetto Python personalizzato, Qiskit Glassy Dynamics (successivamente riferito come QuGlassyIsing), per costruire modularmente le Hamiltoniane per entrambi i reticoli 1D e 2D. Questo pacchetto gestisce l' generazione dei termini di interazione ($\sigma^z\sigma^z$) e dei termini di campo ($\sigma^x, \sigma^z$) in modo efficiente.
• Ambiente di Simulazione: Le simulazioni sono state condotte principalmente utilizzando il backend Qiskit Aer GPU. Una dimostrazione di prova di concetto è stata inoltre eseguita sul dispositivo reale a 7 qubit IBM Q Oslo utilizzando un modello ridotto a 7 spin.

Risultati Chiave

• Paesaggi Energetici: Le simulazioni hanno rivelato che l'aumento del campo longitudinale ($B_z$) abbassa l'energia del sistema allineando gli spin, mentre l'aumento del campo trasverso ($B_x$) induce fluttuazioni quantistiche e inversioni di spin senza abbassare l'energia in modo altrettanto significativo.
• Fasi Disordinate: Lo studio ha identificato un regime specifico in cui emergono configurazioni di spin disordinate. Questa "fase disordinata" è caratterizzata da una miscela di domini ferromagnetici e antiferromagnetici.
  • In assenza di un campo longitudinale ($B_z=0$), il sistema ha esibito un perfetto ordinamento antiferromagnetico.
  • Uno stato disordinato è stato osservato solo per combinazioni specifiche di $B_x$ e $B_z$ (ad esempio, $B_z=1, B_x=2$), dove la competizione tra i campi interrompe l'ordinamento uniforme, creando momenti locali di ferromagnetismo all'interno di uno sfondo antiferromagnetico.
  • Gli autori notano che le tre configurazioni di spin più probabili nelle simulazioni 2D erano quasi identiche, suggerendo un ordine coerente e strutturato anche all'interno della fase disordinata.
• Dimostrazione Hardware: L'esecuzione sul dispositivo IBM Q Oslo (7 qubit) ha mostrato un profilo di convergenza dell'energia con fluttuazioni significative e un valore di energia finale ($-4.0139$) superiore al minimo teorico. Questo risultato è attribuito agli errori di gate, agli errori di lettura e alla decoerenza dei qubit inerenti ai dispositivi NISQ (Noisy Intermediate-Scale Quantum).

Significatività e Rivendicazioni
Il documento sostiene di fornire una dimostrazione di "prova di principio" che i simulatori quantistici a breve termine possono rilevare fasi disordinate nei modelli di Ising che sono difficili da accedere tramite metodi classici. Gli autori sostengono che le loro scoperte:

• Validano l'Utilità del VQE: Dimostrano che il VQE, anche con ansatz non entangled, può modellare funzioni d'onda complesse in tempo polinomiale, riducendo i costi computazionali rispetto ai metodi classici.
• Elucidano la Dinamica Vetrosa: I risultati suggeriscono che l'interazione tra i campi trasverso e longitudinale crea un paesaggio energetico "accidentato" con fasi meta-stabili, riflettendo il comportamento degli spin glass e dei processi catalitici dove il disordine porta a meccanismi di reazione inefficienti.
• Sviluppo di Strumenti: L'introduzione del pacchetto QuGlassyIsing offre un'interfaccia modulare e user-friendly per la costruzione di Hamiltoniane per la dinamica vetrosa, facilitando la ricerca futura in sistemi quantistici complessi.
• Prospettive Future: Gli autori dichiarano con modestia che, sebbene i limiti attuali del NISQ (tempo di coerenza, fedeltà dei gate) impediscano la scalabilità verso sistemi più grandi e complessi, il loro lavoro pone le basi per future applicazioni di calcolo quantistico tollerante ai guasti (FTQC). Propongono che, man mano che la correzione degli errori quantistici migliorerà, questi metodi potrebbero sbloccare le simulazioni di comportamento vetroso in sistemi 3D e geometrie più complesse, potenzialmente informando la progettazione di nuovi materiali di memoria magnetica e sistemi catalitici.

Il documento conclude che, sebbene l'implementazione attuale sia limitata dai vincoli hardware, essa evidenzia con successo il potenziale degli strumenti di elaborazione dell'informazione quantistica per portare a una comprensione sfumata della dinamica e della struttura dei sistemi vetrosi.