2025-11-18T17:19:13.403952

Digital quantum simulation of many-body localization crossover in a disordered kicked Ising model

Hayata, Seki, Yunoki
Simulating nonequilibrium dynamics of quantum many-body systems is one of the most promising applications of quantum computers. However, a faithful digital quantum simulation of the Hamiltonian evolution is very challenging in the present noisy quantum devices. Instead, nonequilibrium dynamics under the Floquet evolution realized by the Trotter decomposition of the Hamiltonian evolution with a large Trotter step size is considered to be a suitable problem for simulating in the present or near-term quantum devices. In this work, we propose simulating the many-body localization crossover as such a nonequilibrium problem in the disordered Floquet many-body systems. As a demonstration, we simulate the many-body localization crossover in a disordered kicked Ising model on a heavy-hex lattice using $60$ qubits from $156$ qubits available in the IBM Heron r2 superconducting qubit device named ibm\_fez. We compute out-of-time-ordered correlators as an indicator of the many-body localization crossover. From the late-time behavior of out-of-time-ordered correlators, we locate the quantum chaotic and many-body localized regimes as a function of the disorder strength. The validity of the results is confirmed by comparing two independent error mitigation methods, that is, the operator renormalization method and zero-noise extrapolation.
academic

Simulazione quantistica digitale del crossover di localizzazione a molti corpi in un modello di Ising calciato disordinato

Informazioni Fondamentali

  • ID Articolo: 2510.01983
  • Titolo: Digital quantum simulation of many-body localization crossover in a disordered kicked Ising model
  • Autori: Tomoya Hayata, Kazuhiro Seki, Seiji Yunoki
  • Classificazione: quant-ph cond-mat.stat-mech hep-lat hep-th
  • Data di Pubblicazione: 15 ottobre 2025
  • Link Articolo: https://arxiv.org/abs/2510.01983

Riassunto

Questo articolo propone un metodo per simulare il fenomeno di crossover della localizzazione a molti corpi (Many-body localization, MBL) su dispositivi quantistici di prossima generazione. I ricercatori hanno utilizzato 60 qubit del processore quantistico superconduttore IBM Heron r2 ibm_fez per simulare il crossover MBL nel modello di Ising calciato disordinato su un reticolo esagonale pesato. Calcolando le funzioni di correlazione fuori ordine temporale (Out-of-time-ordered correlators, OTOCs) come indicatore del crossover MBL, hanno identificato le regioni di caos quantistico e localizzazione a molti corpi dal comportamento tardivo. I risultati sono stati verificati mediante due metodi indipendenti di mitigazione degli errori: la rinormalizzazione degli operatori e l'estrapolazione a rumore zero.

Contesto di Ricerca e Motivazione

Problemi Fondamentali

  1. Simulazione della dinamica non-equilibrio di sistemi quantistici a molti corpi: Questa rappresenta una delle applicazioni più promettenti dei computer quantistici, ma la simulazione fedele dell'evoluzione hamiltoniana su dispositivi quantistici rumorosi attuali è estremamente impegnativa.
  2. Fenomeno di localizzazione a molti corpi: MBL è un meccanismo che impedisce la termalizzazione di sistemi quantistici a molti corpi disordinati, ed è fondamentale per comprendere la rottura dell'ergodicità nei sistemi quantistici.
  3. Sistemi Floquet: I sistemi quantistici periodicamente guidati possono esibire fasi singolari non realizzabili in sistemi all'equilibrio, come i cristalli temporali.

Motivazione della Ricerca

  • Motivazione Tecnica: Il rumore nei dispositivi quantistici attuali limita la simulazione precisa dell'evoluzione hamiltoniana, mentre l'evoluzione Floquet realizzata mediante decomposizione di Trotter è più adatta ai dispositivi quantistici di prossima generazione
  • Motivazione Scientifica: La transizione MBL coinvolge dinamiche non-equilibrio e lo spettro completo degli autostati, difficili da trattare per i computer classici su sistemi grandi, rendendo la simulazione quantistica naturalmente vantaggiosa
  • Motivazione Pratica: Esplorare nuove direzioni per applicazioni di prossima generazione prima della realizzazione del calcolo quantistico tollerante ai guasti

Contributi Principali

  1. Prima realizzazione della simulazione quantistica digitale del crossover MBL su un sistema a 60 qubit, con una scala di sistema che supera le capacità dei metodi classici di diagonalizzazione esatta
  2. Proposta di un metodo di identificazione del crossover MBL basato su OTOCs, localizzando le regioni di caos quantistico e localizzazione a molti corpi attraverso il comportamento tardivo
  3. Sviluppo del modello di Ising calciato disordinato su reticolo esagonale pesato, sfruttando pienamente la connettività nativa e le operazioni di gate dei dispositivi quantistici IBM
  4. Istituzione di due metodi indipendenti di mitigazione degli errori: rinormalizzazione degli operatori ed estrapolazione a rumore zero, con verifica incrociata dell'affidabilità dei risultati
  5. Dimostrazione della fattibilità della simulazione della dinamica non-equilibrio di sistemi Floquet a molti corpi su dispositivi quantistici di prossima generazione

Spiegazione Dettagliata dei Metodi

Definizione del Compito

Studio del fenomeno di crossover della localizzazione a molti corpi nel modello di Ising calciato disordinato su reticolo esagonale pesato, identificando il confine tra la regione di caos quantistico (W < Wc) e la regione di localizzazione a molti corpi (W > Wc) mediante variazione dell'intensità del disordine W.

Architettura del Modello

Modello di Ising Calciato Disordinato

L'hamiltoniana dipendente dal tempo è definita come:

H(t)={HX=i=1NBX(i)Xi,t[0,T/2)HZ=J(i,j)ZiZj+BZi=1NZi,t[T/2,T)H(t) = \begin{cases} H_X = \sum_{i=1}^N B_X(i)X_i, & t \in [0,T/2) \\ H_Z = J\sum_{(i,j)} Z_iZ_j + B_Z\sum_{i=1}^N Z_i, & t \in [T/2,T) \end{cases}

dove:

  • XiX_i, ZiZ_i sono gli operatori di Pauli nel sito i-esimo
  • BX(i)B_X(i) è un campo casuale uniformemente distribuito in [BX0W,BX0+W][B_{X0}-W, B_{X0}+W]
  • JJ, BZB_Z, BX0B_{X0} sono parametri fissi
  • WW è il parametro di intensità del disordine

Operatore di Evoluzione Floquet

L'operatore di evoluzione per un singolo ciclo di guida: UF=eiHZT/2eiHXT/2U_F = e^{-iH_ZT/2}e^{-iH_XT/2}

Impostazione dei parametri: JT=π/2JT = π/2, BZT=1.3B_ZT = 1.3, BX0T=π/2B_{X0}T = π/2, assicurando che il sistema sia nello stato di massimo caos quando W=0W=0.

Funzioni di Correlazione Fuori Ordine Temporale (OTOCs)

Definite come: OTOC=0N(UF)nX1(UF)nZm(UF)nX1(UF)n0N\text{OTOC} = \langle 0^N|(U_F)^{\dagger n}X_1(U_F)^nZ_m(U_F)^{\dagger n}X_1(U_F)^n|0^N\rangle

Forma normalizzata: OTOC=0N(UF)nX1(UF)nZm(UF)nX1(UF)n0N0N(UF)nI1(UF)nZm(UF)nI1(UF)n0N\text{OTOC} = \frac{\langle 0^N|(U_F)^{\dagger n}X_1(U_F)^nZ_m(U_F)^{\dagger n}X_1(U_F)^n|0^N\rangle}{\langle 0^N|(U_F)^{\dagger n}I_1(U_F)^nZ_m(U_F)^{\dagger n}I_1(U_F)^n|0^N\rangle}

Punti di Innovazione Tecnica

  1. Ottimizzazione del Circuito: Rimozione manuale di gate ridondanti al di fuori del cono causale, migliorando significativamente la qualità della mitigazione degli errori
  2. Volume Quantistico Effettivo: Utilizzo del termine denominatore per stimare la crescita del cono causale, come indicatore complementare di MBL
  3. Doppia Mitigazione degli Errori: Combinazione dei metodi di rinormalizzazione degli operatori e ZNE per verifica incrociata
  4. Analisi Dipendente dalla Distanza: Identificazione delle caratteristiche di propagazione del fronte d'onda attraverso il comportamento dell'OTOC a diverse distanze x

Configurazione Sperimentale

Dispositivo Quantistico

  • Dispositivo: Processore quantistico superconduttore IBM Heron r2 ibm_fez (156 qubit)
  • Scala Utilizzata: 60 qubit che costituiscono un reticolo esagonale pesato 2×3
  • Tassi di Errore di Gate: Gate CZ3.7×10⁻³, gate SX2.93×10⁻⁴, errore di lettura~1.76×10⁻²

Parametri Sperimentali

  • Numero di Realizzazioni di Disordine: 25
  • Passi di Trotter: 10
  • Numero Totale di Circuiti: 250 (25×10)
  • Misurazioni per Circuito: 16.000 volte
  • Circuiti di Avvolgimento di Pauli: 32

Metriche di Valutazione

  1. OTOC Normalizzato: Indicatore principale, riflette la propagazione della perturbazione di informazione
  2. Volume Quantistico Effettivo: Veff=logF/log(1p)V_{eff} = \log F / \log(1-p), caratterizza la dimensione del cono causale
  3. Localizzazione del Punto di Crossover: Wc=argmaxWOTOC/WW_c = \arg\max_W \partial\langle\langle\text{OTOC}\rangle\rangle/\partial W

Metodi di Mitigazione degli Errori

  1. Rinormalizzazione degli Operatori: Utilizzo del termine denominatore per correggere gli effetti del rumore
  2. Estrapolazione a Rumore Zero (ZNE): Estrapolazione attraverso fattori di amplificazione del rumore f=1.0, 1.5

Risultati Sperimentali

Risultati Principali

Caratteristiche del Comportamento dell'OTOC

  • Regione di Debole Disordine (W≤0.18): L'OTOC diminuisce drasticamente alla distanza x=n, con forte decadimento tardivo, esibendo comportamento di caos quantistico
  • Regione di Forte Disordine (W≥0.18): L'OTOC rimane vicino a 1 o decade lentamente, mostrando caratteristiche di localizzazione a molti corpi
  • Punto di Crossover: Wc0.18W_c \sim 0.18, con risultati coerenti tra i due metodi di mitigazione degli errori

Analisi della Propagazione del Fronte d'Onda

  • W=0.05: Fronte d'onda chiaro che si propaga vicino a x=n
  • W=0.5: Il fronte d'onda scompare, la propagazione dell'informazione è limitata

Volume Quantistico Effettivo

  • Disordine Debole: Crescita superlineare tardiva, riflettendo l'espansione rapida del cono causale
  • Disordine Forte: Crescita lineare, indicando che il cono causale non si espande più

Verifica della Mitigazione degli Errori

I due metodi di mitigazione degli errori producono risultati coerenti nelle regioni di caos e MBL, mentre la regione di transizione (0.14≤W≤0.2) richiede un'analisi più cauta. Il metodo ZNE rende il comportamento della transizione di fase più evidente, ma con barre di errore relativamente più grandi.

Scoperte Sperimentali

  1. Dipendenza dalla Distanza: L'OTOC a x=n/2=5 è il più adatto come indicatore di crossover
  2. Effetti di Dimensione del Sistema: Il sistema a 60 qubit esibisce una scala irraggiungibile dal calcolo classico
  3. Robustezza al Rumore: L'OTOC normalizzato mostra una certa resistenza al rumore

Lavori Correlati

Ricerca Teorica su MBL

  • Quadro teorico basato sulla rottura dell'ipotesi di termalizzazione degli autostati
  • Comportamento del cono di luce logaritmico nei sistemi unidimensionali
  • Controversie sulla stabilità di MBL in sistemi ad alta dimensionalità

Esperimenti di Simulazione Quantistica

  • Realizzazione di cristalli temporali su qubit superconduttori
  • Identificazione di operatori di movimento integrale locale
  • Applicazione di risolutori variazionali di autostati dello stato eccitato

Sistemi Floquet

  • Fasi singolari di sistemi periodicamente guidati
  • Equivalenza tra dinamica di Trotter e sistemi Floquet
  • Simulazione quantistica di fenomeni di preriscaldamento

Conclusioni e Discussione

Conclusioni Principali

  1. Simulazione riuscita del crossover MBL a 60 qubit, localizzando il punto di crossover Wc0.18W_c \sim 0.18
  2. Verifica dell'efficacia dell'OTOC come indicatore di MBL, distinguendo le regioni di caos e localizzazione
  3. Dimostrazione della capacità pratica dei dispositivi quantistici di prossima generazione, fornendo nuove direzioni per applicazioni pre-tolleranza ai guasti

Limitazioni

  1. Limitazione dei Passi di Trotter: I vincoli del rumore impediscono l'osservazione del comportamento asintotico a lungo termine
  2. Singola Dimensione di Sistema: Assenza di analisi di scaling a dimensione finita, impossibilità di determinare l'esistenza di una vera transizione di fase
  3. Stabilità di MBL Bidimensionale: Controversie teoriche sulla stabilità della fase MBL in sistemi ad alta dimensionalità
  4. Limitazione del Cono Causale: Quando il cono causale copre tutti i qubit, non è possibile estendere ulteriormente il circuito

Direzioni Future

  1. Dispositivi a Rumore Inferiore: Esplorazione del comportamento del cono di luce logaritmico a tempi più lunghi
  2. Scaling a Dimensione Finita: Ricerca su più dimensioni di sistema per determinare la natura della transizione di fase
  3. Altre Fasi Non-Equilibrio: Estensione a quasi-cristalli temporali e altre fasi singolari
  4. Problemi Classicamente Difficili: Ricerca di ulteriori applicazioni di vantaggio quantistico

Valutazione Approfondita

Punti di Forza

  1. Innovazione Tecnica: Prima realizzazione della simulazione MBL su dispositivi quantistici su larga scala, con metodi di ottimizzazione del circuito e mitigazione degli errori di applicabilità generale
  2. Completezza Sperimentale: 25 realizzazioni di disordine, doppia mitigazione degli errori, analisi multi-distanza assicurano l'affidabilità dei risultati
  3. Significato Fisico: Identificazione riuscita del crossover MBL, fornendo una nuova piattaforma sperimentale per la fisica quantistica a molti corpi
  4. Contributo Metodologico: Istituzione di un paradigma per lo studio della dinamica non-equilibrio su dispositivi quantistici rumorosi

Insufficienze

  1. Profondità Teorica: Mancanza di previsioni teoriche e confronti per i valori del punto di crossover
  2. Analisi di Scaling: La singola dimensione di sistema limita la comprensione della natura della transizione di fase
  3. Modello di Rumore: L'ipotesi del canale di depolarizzazione potrebbe essere eccessivamente semplificata
  4. Universalità: L'applicabilità del metodo ad altri dispositivi quantistici e modelli rimane da verificare

Impatto

  1. Valore Accademico: Fornisce un importante benchmark sperimentale per il campo dell'intersezione tra fisica quantistica a molti corpi e calcolo quantistico
  2. Progresso Tecnologico: I metodi di mitigazione degli errori e ottimizzazione del circuito possono essere generalizzati ad altri problemi di simulazione quantistica
  3. Prospettive di Applicazione: Dimostra il potenziale del calcolo scientifico dei dispositivi quantistici di prossima generazione
  4. Riproducibilità: I parametri sperimentali dettagliati e i dati aperti supportano la riproduzione dei risultati

Scenari Applicabili

  1. Ricerca sulla Fisica Quantistica a Molti Corpi: Dinamica non-equilibrio, transizioni di fase, fenomeni critici
  2. Sviluppo di Algoritmi Quantistici: Tecniche di mitigazione degli errori e ottimizzazione del circuito
  3. Valutazione di Dispositivi Quantistici: Test di benchmark e caratterizzazione delle prestazioni
  4. Educazione e Formazione: Casi di studio per l'insegnamento della simulazione quantistica e della fisica a molti corpi

Bibliografia

Questo articolo cita 48 importanti riferimenti, coprendo:

  • Fondamenti teorici di MBL 10-17, 25, 29-33
  • Esperimenti di simulazione quantistica 19, 26-28, 35
  • Sistemi Floquet 1-6, 15-17
  • Metodi di mitigazione degli errori 22, 34, 35, 44
  • Teoria del caos quantistico 36-39

I riferimenti chiave includono la rassegna di MBL di Nandkishore & Huse, la prova del vantaggio quantistico del team IBM, e i fondamenti teorici dell'OTOC.


Valutazione Complessiva: Questo è un articolo sperimentale di importanza significativa nel campo dell'intersezione tra fisica quantistica a molti corpi e calcolo quantistico. Gli autori hanno realizzato con successo la simulazione digitale del crossover MBL su un dispositivo quantistico su larga scala, con metodi tecnici innovativi e risultati affidabili. Nonostante alcune limitazioni nella profondità teorica e nella dimensione del sistema, apre nuove direzioni per le applicazioni scientifiche dei dispositivi quantistici di prossima generazione, con significativo valore accademico e impatto tecnologico.