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

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

Basic Information

  • Paper ID: 2510.01983
  • Title: Digital quantum simulation of many-body localization crossover in a disordered kicked Ising model
  • Authors: Tomoya Hayata, Kazuhiro Seki, Seiji Yunoki
  • Classification: quant-ph cond-mat.stat-mech hep-lat hep-th
  • Publication Date: October 15, 2025
  • Paper Link: https://arxiv.org/abs/2510.01983

Abstract

This paper proposes a method for simulating the many-body localization (MBL) crossover phenomenon on near-term quantum devices. The researchers utilized 60 qubits on the IBM Heron r2 superconducting quantum processor ibm_fez to simulate the MBL crossover in a disordered kicked Ising model on a heavy hexagonal lattice. By computing out-of-time-ordered correlators (OTOCs) as indicators of the MBL crossover, quantum chaotic and many-body localized regions are identified from their late-time behavior. The results are verified through two independent error mitigation methods: operator renormalization and zero-noise extrapolation.

Research Background and Motivation

Core Problems

  1. Nonequilibrium dynamics simulation in quantum many-body systems: This represents one of the most promising applications of quantum computers, yet faithful simulation of Hamiltonian evolution on current noisy quantum devices remains highly challenging.
  2. Many-body localization phenomenon: MBL is a mechanism that prevents thermalization in disordered quantum many-body systems, and is crucial for understanding the breaking of ergodicity in quantum systems.
  3. Floquet systems: Periodically driven quantum systems can exhibit exotic phases inaccessible to equilibrium systems, such as time crystals.

Research Motivation

  • Technical motivation: Noise limitations in existing quantum devices restrict accurate Hamiltonian evolution simulation, while Floquet evolution via Trotter decomposition is better suited to near-term quantum devices
  • Scientific motivation: MBL transitions involve nonequilibrium dynamics and complete eigenspectra, which are difficult for classical computers to handle for large systems, making quantum simulation advantageous
  • Practical motivation: Exploring new directions for near-term applications before achieving fault-tolerant quantum computing

Core Contributions

  1. First implementation of digital quantum simulation of MBL crossover on a 60-qubit system, exceeding the capability of classical exact diagonalization methods
  2. Proposed an OTOC-based method for identifying MBL crossover, locating quantum chaotic and many-body localized regions through late-time behavior
  3. Developed a disordered kicked Ising model on heavy hexagonal lattice, fully exploiting the native connectivity and gate operations of IBM quantum devices
  4. Established two independent error mitigation methods: operator renormalization and zero-noise extrapolation, providing mutual verification of result reliability
  5. Demonstrated the feasibility of simulating nonequilibrium dynamics of Floquet many-body systems on near-term quantum devices

Detailed Methodology

Task Definition

Study the many-body localization crossover in a disordered kicked Ising model on a heavy hexagonal lattice. By tuning the disorder strength W, identify the boundary between the quantum chaotic region (W < W_c) and the many-body localized region (W > W_c).

Model Architecture

Disordered Kicked Ising Model

The time-dependent Hamiltonian is defined as:

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}$$ Where: - $X_i$, $Z_i$ are Pauli operators at the i-th lattice site - $B_X(i)$ is a random field uniformly distributed in $[B_{X0}-W, B_{X0}+W]$ - $J$, $B_Z$, $B_{X0}$ are fixed parameters - $W$ is the disorder strength parameter #### Floquet Evolution Operator The evolution operator for a single driving period: $$U_F = e^{-iH_ZT/2}e^{-iH_XT/2}$$ Parameter settings: $JT = π/2$, $B_ZT = 1.3$, $B_{X0}T = π/2$, such that the system is in maximum chaos state when $W=0$. #### Out-of-Time-Ordered Correlators (OTOCs) Defined as: $$\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$$ Normalized form: $$\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}$$ ### Technical Innovations 1. **Circuit optimization**: Manual removal of redundant gates outside the causal cone, significantly improving error mitigation quality 2. **Effective quantum volume**: Utilizing denominator terms to estimate causal cone growth as a complementary MBL indicator 3. **Dual error mitigation**: Combining operator renormalization and ZNE methods for mutual verification 4. **Distance-dependent analysis**: Identifying wavefront propagation characteristics through OTOC behavior at different distances x ## Experimental Setup ### Quantum Device - **Device**: IBM Heron r2 superconducting quantum processor ibm_fez (156 qubits) - **Usage scale**: 60 qubits forming a 2×3 heavy hexagonal lattice - **Gate error rates**: CZ gate ~3.7×10⁻³, SX gate ~2.93×10⁻⁴, readout error ~1.76×10⁻² ### Experimental Parameters - **Number of disorder realizations**: 25 - **Trotter steps**: 10 - **Total circuits**: 250 (25×10) - **Measurements per circuit**: 16,000 - **Pauli string circuits**: 32 ### Evaluation Metrics 1. **Normalized OTOC**: Primary metric reflecting information perturbation propagation 2. **Effective quantum volume**: $V_{eff} = \log F / \log(1-p)$, characterizing causal cone size 3. **Crossover point localization**: $W_c = \arg\max_W \partial\langle\langle\text{OTOC}\rangle\rangle/\partial W$ ### Error Mitigation Methods 1. **Operator renormalization**: Using denominator terms to correct noise effects 2. **Zero-noise extrapolation (ZNE)**: Extrapolation through noise amplification factors f=1.0, 1.5 ## Experimental Results ### Main Results #### OTOC Behavior Characteristics - **Weak disorder region (W≤0.18)**: OTOC drops sharply at distance x=n, exhibits strong late-time decay, showing quantum chaotic behavior - **Strong disorder region (W≥0.18)**: OTOC remains close to 1 or decays slowly, displaying many-body localization characteristics - **Crossover point**: $W_c \sim 0.18$, consistent between two error mitigation methods #### Wavefront Propagation Analysis - **W=0.05**: Clear wavefront propagating near x=n - **W=0.5**: Wavefront disappears, information propagation is suppressed #### Effective Quantum Volume - **Weak disorder**: Superlinear growth at late times, reflecting rapid causal cone expansion - **Strong disorder**: Linear growth, indicating causal cone no longer expands ### Error Mitigation Verification Both error mitigation methods show consistent results in chaotic and MBL regions, with careful analysis needed in the crossover region (0.14≤W≤0.2). The ZNE method makes the phase transition behavior clearer but with relatively larger error bars. ### Experimental Findings 1. **Distance dependence**: OTOC at x=n/2=5 is most suitable as a crossover indicator 2. **System size effects**: The 60-qubit system demonstrates scales inaccessible to classical computation 3. **Noise robustness**: Normalized OTOC shows certain resistance to noise ## Related Work ### MBL Theoretical Research - Theoretical framework based on breakdown of eigenstate thermalization hypothesis - Logarithmic light-cone behavior in one-dimensional systems - Controversies regarding MBL stability in higher dimensions ### Quantum Simulation Experiments - Time crystal implementation on superconducting qubits - Identification of local integral-of-motion operators - Application of excited-state variational quantum eigensolver ### Floquet Systems - Exotic phases in periodically driven systems - Equivalence between Trotter dynamics and Floquet systems - Quantum simulation of prethermalization phenomena ## Conclusions and Discussion ### Main Conclusions 1. **Successfully implemented 60-qubit MBL crossover simulation**, locating the crossover point $W_c \sim 0.18$ 2. **Verified OTOC as an effective MBL indicator**, distinguishing chaotic and localized regions 3. **Demonstrated practical capability of near-term quantum devices**, providing new directions for pre-fault-tolerant applications ### Limitations 1. **Trotter step constraints**: Noise limitations prevent observation of long-time asymptotic behavior 2. **Single system size**: Lack of finite-size scaling analysis prevents determination of true phase transitions 3. **Two-dimensional MBL stability**: Theoretical controversy regarding MBL phase stability in higher dimensions 4. **Causal cone limitations**: Circuit extension becomes impossible when causal cone covers all qubits ### Future Directions 1. **Lower-noise devices**: Exploring longer-time logarithmic light-cone behavior 2. **Finite-size scaling**: Multi-system-size studies to determine phase transition properties 3. **Other nonequilibrium phases**: Extension to exotic phases such as time quasicrystals 4. **Classically intractable problems**: Seeking more quantum advantage applications ## In-Depth Evaluation ### Strengths 1. **Technical innovation**: First MBL simulation on large-scale quantum devices with universally applicable circuit optimization and error mitigation methods 2. **Experimental rigor**: 25 disorder realizations, dual error mitigation, multi-distance analysis ensure result reliability 3. **Physical significance**: Successful identification of MBL crossover provides new experimental platform for quantum many-body physics 4. **Methodological contribution**: Establishes paradigm for studying nonequilibrium dynamics on noisy quantum devices ### Weaknesses 1. **Theoretical depth**: Lacks theoretical predictions and comparisons for crossover point values 2. **Scaling analysis**: Single system size limits understanding of phase transition properties 3. **Noise model**: Depolarizing channel assumption may be oversimplified 4. **Generalizability**: Applicability of methods to other quantum devices and models remains to be verified ### Impact 1. **Academic value**: Provides important experimental benchmark for quantum many-body physics and quantum computing intersection 2. **Technical advancement**: Error mitigation and circuit optimization methods generalizable to other quantum simulation problems 3. **Application prospects**: Demonstrates scientific computing potential of near-term quantum devices 4. **Reproducibility**: Detailed experimental parameters and open data support result reproduction ### Applicable Scenarios 1. **Quantum many-body physics research**: Nonequilibrium dynamics, phase transitions, critical phenomena 2. **Quantum algorithm development**: Error mitigation and circuit optimization techniques 3. **Quantum device evaluation**: Benchmarking and performance characterization 4. **Education and training**: Teaching cases for quantum simulation and many-body physics ## References This paper cites 48 important references covering: - MBL theoretical foundations [10-17, 25, 29-33] - Quantum simulation experiments [19, 26-28, 35] - Floquet systems [1-6, 15-17] - Error mitigation methods [22, 34, 35, 44] - Quantum chaos theory [36-39] Key references include Nandkishore & Huse's MBL review, IBM team's quantum advantage demonstration, and OTOC theoretical foundations. --- **Overall Assessment**: This is an important experimental paper at the intersection of quantum many-body physics and quantum computing. The authors successfully implemented digital simulation of MBL crossover on a large-scale quantum device with innovative technical methods and reliable results. Despite certain limitations in theoretical depth and system size, it opens new directions for scientific applications of near-term quantum devices, possessing significant academic value and technical advancement potential.