论可调破缺可积性下量子多体混沌的涌现

论文信息

标题: On the emergence of quantum many-body chaos for tunably-broken integrability

作者: Sounak Biswas, Sthitadhi Roy, Roderich Moessner

发布日期: 2026-07-02

arXiv ID: 2607.02506v1

PDF 链接: 下载 PDF

3 分钟速览

• 研究问题：本文研究当可积性通过一个可调参数 λ 被连续打破时，量子多体混沌如何从早期可积行为逐步涌现，并希望定量描述这一跨接过程的时空标度。
• 核心方法：构造一维 qubit 链上的砖墙型量子线路，门集为可积的匹配门与概率 λ 掺杂的非可积 SWAP 门；计算 out‑of‑time‑ordered correlators (OTOC) 并将系综平均后的动力学精确映射到经典马尔可夫过程，进而在连续极限下导出 noisy Fisher‑KPP 方程。
• 关键结果：跨接由特征时间 $t^*(\lambda)\sim\lambda^{-1}$ 和长度 $r^*(\lambda)\sim\lambda^{-1/2}$ 控制；SWAP 门充当局域 “热点”，反复放大 OTOC，最终在时空上累积出弹道传播的混沌前沿，蝴蝶速度 $v_B\sim\sqrt{\lambda}$（带对数修正）。
• 主要局限：模型依赖门系综平均，丢失了单次实验中的量子相干性与环绕效应；未讨论有限温度、能量输运以及真实实验平台中的退相干和门噪声。
• 适合读者：从事量子多体混沌、可积系统破缺、量子线路与 OTOC 理论研究的凝聚态物理和量子信息学者，尤其关注从可积到混沌跨接机制的理论工作者。

论文背景和研究动机

新一代量子模拟平台带来了长程量子相干时间和前所未有的动力学观测窗口，但也对理论模型提出了更好要求：模型应具有有限局域希尔伯特空间（如 qubit）、由局域幺正门构成通用门集，并且可以可调地接近可积性。然而，现有的精确可解模型往往需要精细调节——随机线路或空间随机哈密顿量要求大量局域 Hilbert 空间，双幺正线路的解则对扰动极脆弱，且这些模型通常从初始起就接近最大混沌。

因此，一个 “既可控又通用” 的理论设置十分需要，以便研究早期可积行为到后期充分发展混沌的跨接物理。该跨接不仅涉及基础问题（什么机制支撑跨接？对应的特征时间和长度尺度如何依赖破缺参数？），而且在诸多实验平台上随处可见。本文正是为此提出了一种掺杂 SWAP 门的匹配门线路模型，用 OTOC 作为混沌的探针，定量建立了跨接理论。

核心方法和技术细节

模型与 OTOC

考虑一维 $L$ 个 qubit（或等效的 Majorana 费米子）上的砖墙型线路，时间步 $t\to t+1$ 由先作用在奇数键、后作用在偶数键上的两层门构成。每步中，作用在相邻位点 $x,x+1$ 上的门以概率 $1-\lambda$ 为随机匹配门 $M_{x,x+1}$，以概率 $\lambda$ 为 SWAP 门夹在两个匹配门之间的形式 $S_{x,x+1}=M'\,\text{SWAP}\,\tilde M$。匹配门只包含费米子的二次型，保证 $\lambda=0$ 时系统自由费米子可积；SWAP 门则引入四次费米子相互作用，破坏可积性。

研究的 OTOC 定义为

$$C_{ZZ}(r,t)=\frac12\mathrm{tr}\big([Z_0(t),Z_r][Z_0(t),Z_r]^\dagger\big),$$

也可写为 Majorana 形式。通过将时间演化算符展开为 Majorana 字符串，OTOC 的系综平均值可以表示为这些字符串权重的概率分布。

经典马尔可夫过程与 noisy F‑KPP 方程

对匹配门系综求平均会使字符串权重均匀分配到所有可能的结果上。由此，整个平均 OTOC 的动力学被精确映射成一个经典马尔可夫过程：把 Majorana 字符串看作 1D 链上的硬核粒子占据态，匹配门驱动粒子随机游走，SWAP 门则充当分支（单粒子 → 三粒子，局域放大 OTOC）和湮灭（三粒子 → 单粒子）的 “热点”。

在连续时空极限下，粗粒化 OTOC $C(r,t)$ 满足 noisy Fisher‑KPP 方程

$$\partial_t C = D\,\partial_r^2 C + \lambda H[C] + \sqrt{\lambda K[C]}\,\eta(r,t),$$

其中

$$H[C]=2C(1-C)(1-2C),\qquad K[C]=C(1-C)(1-2C+2C^2),$$

$\eta$ 为白噪声。第一项来自匹配门的扩散（扩散系数 $D=2$），第二项是非线性增益（SWAP 门导致局域 OTOC 放大，$H[C]\sim 2(p_1-p_3)$，$p_1,p_3$ 分别是一个和三个空位的概率），第三项则是由离散性带来的人口统计噪声，导致混沌前沿的扩散展宽。

创新点和贡献

1. 通用可调破缺的可积模型：首次在 qubit 系统中用掺杂 SWAP 门连续调节可积性，模型既保持有限 Hilbert 空间，又形成通用门集，可在整个 λ 范围内得到统计精确的结果。
2. 经典马尔可夫过程的精确映射：将 OTOC 的平均动力学等价于硬核粒子在 SWAP 门处的分支/湮灭过程，使得数值模拟可轻易扩展到任意大系统，并严格推出了非线性微分方程。
3. 跨接尺度的解析识别：明确给出了从可积扩散传播到混沌弹道传播的两个特征标度 $t^*(\lambda)\sim\lambda^{-1}$ 和 $r^*(\lambda)\sim\lambda^{-1/2}$。在 $r\ll r^*$ 区域，OTOC 会先经历可积极限的局部时间极大值，然后才被混沌前沿追上；在 $r\gg r^*$ 区域，到达的 OTOC 前沿已经是混沌的。
4. 微观机制与 F‑KPP 方程的物理联系：把 SWAP 门解释为 “热点”，其累积效果导致全局混沌，而 F‑KPP 方程将 OTOC 的前沿传播与著名的拉涨波等同，自然地给出了蝴蝶速度 $v_B\sim\sqrt{\lambda}$（带对数修正）以及前沿展宽扩散常数几乎不依赖 λ 等性质。

实验结果分析

论文通过模拟经典马尔可夫过程得到以下关键数值结果。

• 可积与混沌的对比（图 1）：λ=0 时，OTOC 呈高斯形扩散，空间剖面 $\propto e^{-r^2/2Dt}$，峰值按 $t^{-1/2}$ 衰减；λ=0.01 时，OTOC 形成弹道光锥，前端以恒定速度 $v_B$ 向外传播，内部饱和。这一鲜明差异直接展示了混沌的涌现。
• 混沌区域的标度形式（图 2）：对于 λ>0，$(r-v_B t)/\sqrt{t}$ 作为标度变量可将不同时间的 OTOC 数据坍塌到同一条曲线上，且同一 $t$ 下不同 λ 的数据也坍塌，表明前沿展宽的扩散常数与 λ 基本无关。提取的 $v_B(\lambda)$ 符合 $\sqrt{\lambda}$ 带对数修正的预测（图 2c）。
• 跨接尺度（图 3、图 4）：以缩放变量 $\tau=\lambda t$ 和 $\rho=r\sqrt{\lambda}$ 表示，OTOC 满足 $C_{ZZ}(r,t)=t^{-1/2}\exp[\tau\,g(\rho/\tau)]$。$r=0$ 处的 OTOC 在 $\tau\approx 1$ 附近明显偏离可积的 $t^{-1/2}$ 衰减；$r=100$ 的数据同样在 $\rho$ 和 $\tau$ 支配下坍塌，证实了 $t^*\sim\lambda^{-1}$ 和 $r^*\sim\lambda^{-1/2}$ 的标度。

noisy F‑KPP 方程的线性化近似（$C\ll 1$ 时 $H[C]\approx 2C$）给出解析解 $C(r,t)\approx e^{2\lambda t(1-r^2/8\lambda t^2)}/\sqrt{4\pi t}$，由此可以推导出 $t^*\sim (1/4\lambda)\ln(1/4\lambda)$，并解释跨接长度 $r^*$ 的出现（式 14–15）。混沌前沿的蝴蝶速度和扩散常数均可从非线性方程的标准理论（Brunet‑Derrida 截断效应）得到解释，与数值非常吻合。

局限与待解决问题

尽管该工作提供了从可积到混沌跨接的完整定量图像，但仍有若干固有局限和待拓展的问题。

• 系综平均与单次实现：OTOC 的解析处理建立在门的系综平均上，而实际实验往往直接测量单次随机电路或特定门序列。单个实现中的量子相干效应和 OTOC 的涨落无法由经典马尔可夫过程描述，也未在文中讨论。
• 忽略退相干与噪声：模型假设完美的幺正门，而现实平台中退相干、门误差等不可避免，这些效应如何修正特征标度目前未知。
• 有限温度与哈密顿体系：模型本质上是一种动力学电路，无限温度假设内建。在具有哈密顿量的可积破缺系统中，能量输运与 OTOC 的关系可能带来新的跨接行为，本文尚未涉及。
• 高阶 OTOC 与设计：论文已展望更高阶 k‑OTOC 可能揭示自由独立性和量子设计的涌现，但未给出具体结果。这类推广是否会引入新的跨接阶次，以及它们如何依赖破缺参数，值得进一步研究。
• 对数修正的可观测性：$v_B$ 的对数修正和前沿扩散常数极弱的对数依赖，在数值上很难与纯幂律明显区分，论文也未给出定量甄别这些修正的实验或数值方案。