来源论文: https://arxiv.org/abs/2605.10758v1 生成时间: May 17, 2026 00:33

1D监测非相互作用费米子系统中测量诱导相变（MIPT）的缺位：基于超大规模GPU模拟与非线性Sigma模型的深度解析

0. 执行摘要

在监测多体量子系统的动力学演化中，纠缠熵（Entanglement Entropy, EE）随系统尺寸 $L$ 的标度行为是区分不同动力学相的核心指标。近年来，关于一维（1D）受监测非相互作用费米子在存在无序或准周期势时是否会发生从对数法则（Logarithmic-law）到面积法则（Area-law）的测量诱导相变（MIPT）存在显著争议。此前有研究声称在有限的监测强度下观察到了此类相变。

本论文（Yin et al., 2026）通过结合高性能GPU加速的数值模拟和解析的非线性Sigma模型（NLSM）分析，从根本上否定了这一观点。作者指出，先前观察到的所谓“相变”实际上是严重的有限尺寸效应。通过将系统尺寸提升至 $L=18000$（比以往研究大一个数量级），研究表明，只要存在有限强度的监测，系统总是处于面积法则相，其临界监测强度 $\gamma_c$ 严格为零。此外，研究还揭示了无序通过改变系统对称性（从 BDI 类变为 AIII 类）反而会增加相关长度的有趣现象。这一发现对于理解开放量子系统的纠缠动力学以及设计抗噪声量子算法具有重要的理论指导意义。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：MIPT 的真实性与普适性

在传统的随机量子电路中，测量可以通过减少系统算符的空间扩展来诱导纠缠从体积法则（Volume-law）到面积法则（Area-law）的转变。然而，对于非相互作用费米子（即自由费米子），情况更为微妙。在纯净（无背景势）的1D Dirac费米子系统中，已有理论（基于 NLSM）预测系统在任何非零监测强度下都应处于面积法则相。但在引入无序（Disorder）或准周期势（Quasi-periodic potential）后，由于安德森定域化（Anderson Localization）的存在，纠缠标度本应进一步受限。之前的研究（Refs [59, 60]）却得出了相反的结论，认为存在一个有限的监测强度阈值 $\gamma_c > 0$ 才能触发相变。本文的核心使命就是厘清：这种“相变”是真实存在的物理现象，还是由于计算资源限制导致的数值假象？

1.2 理论基础：非线性Sigma模型（NLSM）与场论映射

作者利用 Keldysh 路径积分表述将监测动力学映射到非线性 Sigma 模型上。这是量子相干输运理论中的标准工具，但在监测动力学中涉及分母中的迹（即量子态的归一化），这导致了非幺正性。通过引入 $R$ 份副本（Replicas）并取 $R \to 1$ 的极限，可以将动力学描述为有效玻色子场论。

  • 对称性类（Symmetry Classes）：这是理解本文结论的关键。在没有无序的纯净限下，系统具有手征对称性（CS）和时间反演对称性（TRS），属于 BDI 类。引入对角无序后，手征对称性被打破，系统变为 AIII 类。这种对称性的改变直接导致了重正化群（RG）流方程中系数的变化，从而改变了相关长度 $l_{cor}$ 的标度函数。
  • RG 流方程：一圈（One-loop）重正化群分析表明，质量项 $\mu$（代表监测强度和无序的联合效应）总是会驱动扩散系数 $D$ 减小，使系统流向强耦合区域。只要 $\mu > 0$，RG 流最终都会指向面积法则相。

1.3 技术难点：超大尺度下的纠缠计算

对于自由费米子，纠缠熵可以通过 $L \times L$ 的单粒子相关矩阵 $D_{ij} = \langle \psi | c_i^\dagger c_j | \psi \rangle$ 获得。虽然计算量从指数级降到了多项式级（$O(L^3)$），但要排除相变，需要极大的 $L$（如本文采用的 $L=18000$）。在每一时间步，都需要进行大量的矩阵乘法和正交化（QR分解），这对传统 CPU 构成了巨大的计算挑战。此外，由于相关长度 $l_{cor}$ 在弱监测强度下呈指数级增长，如果系统尺寸 $L$ 不远大于 $l_{cor}$，数值上会表现出虚假的对数增长，误导研究者认为存在临界点。

1.4 方法细节：监测协议与标度分析

  • QSD 协议（量子态扩散）：对应于弱、连续的测量。通过随机薛定谔方程演化态矢量。本文主要用于无序势研究。
  • PM 协议（投影测量）：对应于随机发生的离散测量。本文主要用于准周期势研究。
  • 纠缠熵二阶累积量展开：利用 $S_A \approx \frac{\pi^2}{3} \int C(x,y) dx dy$ 的近似（其中 $C$ 是密度-密度相关函数），将复杂的纠缠熵计算转化为对相关函数衰减特性的分析。通过拟合 $C(r) \propto \exp(-r/l_{cor})$ 提取相关长度。

2. 关键 Benchmark 体系、计算数据与性能表现

2.1 Benchmark 体系：Anderson 定域化与 Aubry-André 模型

  1. 无序势（Disordered Potential）：采用方盒分布 $[-W, W]$。在 $1D$ 限下，无论 $W$ 多小，所有本征态在非监测限下都是定域的。
  2. 准周期势（Quasiperiodic Potential）：采用 $V_i = V \cos(2\pi i / \tau + \theta)$，其中 $\tau$ 为黄金比例。这是著名的 Aubry-André 模型，在 $V>2J$ 时发生金属-绝缘体相变。

2.2 核心数据分析

  • 相关长度 $l_{cor}$ 的演化（图1 & 图2）：数据清晰显示，随着监测强度 $\gamma$ 降低，$l_{cor}$ 呈指数增长。对于不同的无序强度 $W$，拟合得到的临界监测强度 $\gamma_c$ 均在误差范围内等于零。例如，对于 $W=0.75$，$\gamma_c = -0.05 \pm 0.07$。
  • 相关函数 $C(r)$ 的衰减（图1插图）：在小尺度（$r < 500$）时，曲线表现出类似幂律的伪特征（$\alpha \simeq -2.2$），这正是早期研究误判的原因。但在大尺度（$r > 2500$）下，曲线最终回归到指数衰减，证实了面积法则。
  • 互信息（Mutual Information, $\mathcal{I}_2$）的跨越消失（图S1）：在 $L=144$ 时，$\mathcal{I}_2$ 在 $\gamma \sim 0.2$ 处存在伪交叉点。但随着 $L$ 增加到 $2584$，该交叉点彻底消失，曲线单调下移。这进一步排除了热力学极限下的相变。

2.3 性能数据与计算资源

  • GPU 加速效率：计算任务在 NVIDIA A100 GPU 上运行。对于 $L=8192$ 的单个量子轨迹演化到稳态，大约需要 11 小时。如果使用传统单核 CPU，该任务可能需要数周时间。作者通过并行化矩阵运算（CUDA）和利用 MATLAB GPUArray 框架，实现了相对于 CPU 至少 10-50 倍的加速，这使得 $L=18000$ 的采样成为可能。

3. 代码实现细节、复现指南与开源链接

3.1 核心算法：相关矩阵演化法

由于系统保持自由费米子特性（二次哈密顿量），态矢量的全波函数信息编码在相关矩阵 $D$ 中。复现的核心步骤如下：

  1. 初始化：构造 $L/2$ 粒子数的 Neel 态或随机填充态对应的 $D$ 矩阵。
  2. 酉演化：利用步长 $dt$ 更新 $D(t+dt) = U D(t) U^\dagger$，其中 $U = \exp(-i H dt)$。由于 $H$ 是三对角或带状矩阵，可利用稀疏矩阵算法或 FFT 加速。
  3. 监测更新：
    • QSD 协议：引入高斯噪声向量 $dW$，通过非线性的随机项更新矩阵 $D$。关键在于更新后必须通过 QR 分解保持矩阵的秩和正交性，防止数值精度崩溃。
    • PM 协议：根据每个格点的密度 $D_{ii}$ 生成随机数，决定是否坍缩。更新公式涉及矩阵的分块操作（参见论文公式 S1-S2）。
  4. 纠缠熵提取：计算 $D$ 的特征值 $\lambda_i$，带入 $S_A = -\sum [\lambda_i \ln \lambda_i + (1-\lambda_i) \ln(1-\lambda_i)]$。

3.2 软件栈与复现指南

  • 语言：C++ (CUDA) 用于大规模运算；MATLAB (GPUArray) 用于快速原型和准周期势模拟。
  • 核心库：cuBLAS, cuSolver (用于快速特征值分解和矩阵乘法)。
  • 数值稳定性控制：在超大规模模拟中，由于测量过程的非幺正性，相关矩阵的特征值可能越界（超出 $[0,1]$）。必须定期通过重新正交化（Orthonormalization）来重置。建议每 10-20 个步长进行一次 QR 分解。

3.3 开源资源链接建议

虽然论文未直接给出 GitHub 链接，但此类问题的标准实现通常参考以下开源项目：

  • QuSpin (Python): 虽然主要用于小规模精确对角化，但其基础框架可借鉴。
  • TensorNetwork (Google): 可用于实现基于张量网络思想的二次费米子演化。
  • 建议检索关键词: “Monitored free fermions CUDA”, “QSD correlation matrix simulation”。

4. 关键引用文献与局限性评论

4.1 关键引用文献

  1. Refs [59, 60]: 先前声称观察到 MIPT 的研究。本文的主要批判对象。通过对比证明了其结论受限于 $L \le 500$。
  2. Ref [41] (Poboiko et al., 2023): 建立了监测费米子与 NLSM 映射的基础理论框架。
  3. Ref [53] (Fan et al., 2026): 作者的前期工作，展示了 GPU 在此类问题中的巨大威力。
  4. Ref [65] (Liao et al., 2026): 一个几乎并行的独立预印本，得出了相似结论，增强了本文的可信度。

4.2 局限性评论

  • 非相互作用限制：本文的所有结论严格限于二次费米子系统。一旦引入电子-电子相互作用，NLSM 将包含相互作用项，系统可能进入不同的对称性类，MIPT 有可能在相互作用系统中以不同的形式重新出现。
  • 相关长度提取的精度：在极弱监测限下，相关长度 $l_{cor}$ 增长极快，即使是 $L=18000$ 也可能只是“刚刚够用”。RG 流分析虽然提供了定性趋势，但在极强无序限下的高阶修正仍待深入探讨。
  • 实验可观测性：纠缠熵在实验中极难直接测量，尤其是需要大量的“后选择”（Post-selection）量子轨迹。虽然本文证明了其数学上的行为，但在近期量子硬件上复现这一面积法则规律仍面临巨大挑战。

5. 补充内容：从对称性角度看无序的“反直觉”作用

5.1 BDI 与 AIII 对称性类的深层含义

本文最有趣的发现之一是：增加无序强度 $W$ 竟然会增加相关长度 $l_{cor}$（在弱无序限下）。这与通常认为的“无序增加定域化，减小纠缠”背道而驰。

其物理机制在于：

  • 在纯净限下（BDI类），系统的 RG 流受到强约束，相关长度的指数增长受到一个因子 $1/(2\gamma)$ 的限制。
  • 在无序限下（AIII类），对角无序打破了格点间的手征对称性（Chiral Symmetry），使得原本受到保护的某些软模（Soft modes）变得更加活跃。RG 流方程的解析解显示，指数项变为 $1/\gamma$。这个两倍的差异意味着，在相同的 $\gamma$ 下，无序系统的相关长度远大于纯净系统。这也就解释了为什么在以往的小规模模拟中，无序系统看起来更像是一个具有对数标度的“临界”相。

5.2 监测作为有效的“温度”

从统计物理角度看，监测可以被视为一种有效温度，它不断地向系统中注入随机性并破坏量子相干性。安德森定域化本质上是基于波函数相干叠加产生的干涉效应，而监测会诱导波函数坍缩，本质上是在“测量”粒子的位置。这种行为在某种程度上与安德森定域化的机制是竞争的。本文的结果表明，这种竞争在纠缠熵的标度层面上是以“监测”的绝对胜利告终——即只要有测量，相干性所维持的长程对数纠缠就必然坍缩为局域的面积法则纠缠。

5.3 未来展望：迈向高维与相互作用系统

虽然 1D 非相互作用系统没有 MIPT，但作者在引言中提到，在更高维度（2D, 3D）或具有长程跳跃的系统中，MIPT 可能依然存在。本项研究所建立的大规模 GPU 模拟范式，为探索这些更复杂体系的纠缠动力学铺平了道路。量子计算研究者应关注如何利用监测来控制系统纠缠，作为一种除了哈密顿量工程之外的新型操控手段。