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

基本信息

论文ID: 2510.01983
标题: Digital quantum simulation of many-body localization crossover in a disordered kicked Ising model
作者: Tomoya Hayata, Kazuhiro Seki, Seiji Yunoki
分类: quant-ph cond-mat.stat-mech hep-lat hep-th
发表时间: 2025年10月15日
论文链接: https://arxiv.org/abs/2510.01983

摘要

本文提出了一种在近期量子设备上模拟多体定域化(Many-body localization, MBL)交叉现象的方法。研究者利用IBM Heron r2超导量子处理器ibm_fez上的60个量子比特，在重六角格子上模拟了无序踢Ising模型中的多体定域化交叉现象。通过计算时间无序关联函数(Out-of-time-ordered correlators, OTOCs)作为MBL交叉的指标，从其晚期行为中识别出量子混沌和多体定域化区域。研究结果通过算符重整化方法和零噪声外推两种独立的误差缓解方法得到验证。

研究背景与动机

核心问题

量子多体系统非平衡动力学模拟：这是量子计算机最有前景的应用之一，但在当前噪声量子设备上忠实模拟哈密顿量演化极具挑战性。
多体定域化现象：MBL是一种阻止无序量子多体系统热化的机制，对理解量子系统的遍历性破缺具有重要意义。
Floquet系统：周期性驱动的量子系统可以展现平衡系统无法实现的奇异相，如时间晶体。

研究动机

技术动机：现有量子设备噪声限制了精确哈密顿量演化的模拟，而Floquet演化通过Trotter分解实现，更适合近期量子设备
科学动机：MBL转变涉及非平衡动力学和完整本征谱，经典计算机难以处理大系统，量子模拟具有天然优势
实用动机：为容错量子计算实现前的近期应用探索新方向

核心贡献

首次在60量子比特系统上实现MBL交叉的数字量子模拟，系统规模超出经典精确对角化方法的能力范围
提出基于OTOCs的MBL交叉识别方法，通过晚期行为定位量子混沌和多体定域化区域
开发了适用于重六角格子的无序踢Ising模型，充分利用IBM量子设备的原生连接性和门操作
建立了两种独立的误差缓解方法：算符重整化和零噪声外推，相互验证结果可靠性
证明了近期量子设备模拟Floquet多体系统非平衡动力学的可行性

方法详解

任务定义

研究无序踢Ising模型在重六角格子上的多体定域化交叉现象，通过调节无序强度W，识别量子混沌区域(W < Wc)和多体定域化区域(W > Wc)的边界。

模型架构

无序踢Ising模型

时间依赖的哈密顿量定义为：

$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}$

其中：

$X_i$ , $Z_i$ 是第i个格点的Pauli算符
$B_X(i)$ 是随机场，均匀分布在 $[B_{X0}-W, B_{X0}+W]$
$J$ , $B_Z$ , $B_{X0}$ 为固定参数
$W$ 为无序强度参数

Floquet演化算符

单个驱动周期的演化算符： $U_F = e^{-iH_ZT/2}e^{-iH_XT/2}$

参数设置： $JT = π/2$ , $B_ZT = 1.3$ , $B_{X0}T = π/2$ ，使得 $W=0$ 时系统处于最大混沌状态。

时间无序关联函数(OTOCs)

定义为： $\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$

归一化形式： $\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}$

技术创新点

电路优化：手动移除因果锥外的冗余门，显著提高误差缓解质量
有效量子体积：利用分母项估算因果锥增长，作为MBL的补充指标
双重误差缓解：结合算符重整化和ZNE方法，提供相互验证
距离依赖分析：通过不同距离x的OTOC行为识别波前传播特征

实验设置

量子设备

设备：IBM Heron r2超导量子处理器ibm_fez (156量子比特)
使用规模：60量子比特构成2×3重六角格子
门误差率：CZ门~~3.7×10⁻³, SX门~~2.93×10⁻⁴, 读出误差~1.76×10⁻²

实验参数

无序实现数：25个
Trotter步数：10步
总电路数：250个 (25×10)
每电路测量次数：16,000次
Pauli缠绕电路数：32个

评价指标

归一化OTOC：主要指标，反映信息扰动传播
有效量子体积： $V_{eff} = \log F / \log(1-p)$ ，表征因果锥大小
交叉点定位： $W_c = \arg\max_W \partial\langle\langle\text{OTOC}\rangle\rangle/\partial W$