来源论文: https://arxiv.org/abs/2603.15744v1 生成时间: Mar 18, 2026 12:21

测量诱导相变中的后选择临界性：深度解析 P-RQC 与随机张量网络的普适类共性

0. 执行摘要

测量诱导相变（Measurement-Induced Phase Transitions, MIPT）是近年来量子多体动力学领域的研究热点，它描述了单一量子轨迹在幺正演化与局部测量竞争下，纠缠熵从体积律（Volume-law）向面积律（Area-law）转换的过程。然而，传统的 MIPT 依赖于符合 Born 定律的测量结果采样，这在实验观测中面临指数级的“后选择”难题。本文通过数值模拟和理论映射，深入探讨了强制后选择（Forced Post-selection）下的临界行为。研究发现：

  1. 后选择根本性地改变了 MIPT 的普适类，其关联长度指数 $\nu \approx 2.1$，显著大于标准 MIPT（$\nu \approx 1.3$）。
  2. 后选择 MIPT 与随机张量网络（RTN）的纠缠相变表现出高度一致的临界指数，证明两者属于同一普适类。
  3. 该普适类具有负的有效中心荷 $c_{\text{eff}} \approx -0.4$，这在失序系统中极为罕见。
  4. 在平移对称的准周期电路中，诱导相变需要至少 3 维的局部希尔伯特空间（qutrits）。

本解析将从理论基础、数值基准、实现细节及科学局限性四个维度对该工作进行深度剖析。

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

1.1 科学背景与核心问题

量子动力学的鲁棒性可以通过局部监测来刻画。在标准 MIPT 中，测量结果根据 Born 定律 $p(m) = \langle \psi | P_m | \psi \rangle$ 随机出现。若要研究特定的量子轨迹，实验者必须在指数多次的重复中寻找相同的测量序列，即“后选择难题”。

本文的核心问题是：如果我们手动“强制”所有的测量结果为特定值（例如 $|0\rangle$ 态），这种显式后选择会如何改变相变的性质？更进一步，这种受迫演化是否与随机张量网络中的几何纠缠相变具有相同的物理起源？

1.2 理论基础：统计力学映射与 Replica Trick

理解 MIPT 的强有力工具是将其映射到二维经典统计力学模型。对于标准 MIPT，我们需要处理分母（归一化因子），这导致了 $k \to 1$ 的 replica 极限。而在后选择情况下（P-RQC），测量结果不再受 Born 概率加权，这相当于在 replica 极限中取 $k \to 0$。

  • 标准 MIPT (Born Rule)：对应于 $S_{n} \sim \lim_{k\to 1} \frac{1}{1-n} \log \text{Tr}(\rho^n)$。在统计力学中，这对应于一种特殊的具有相互作用的自旋模型，其对称性通常与排列群相关。
  • 后选择 MIPT (P-RQC)：对应于 $k=0$ 的极限。在这种情况下，配分函数 $Z_X$ 仅由演化的非幺正算符序列决定，而不需要考虑波函数的投影概率。这使得它在数学结构上直接等同于随机张量网络（RTN）的收缩问题。

1.3 技术难点：非幺正动力学下的临界提取

后选择过程引入了非幺正性，波函数的模长会随时间指数衰减。数值上，必须频繁重归一化以防止溢出，但这会掩盖系统的总权重信息。为了提取有效中心荷 $c_{\text{eff}}$，研究者必须精确计算准自由能密度 $f(L)$ 的有限尺寸修正，这要求极高的数值稳定性。

1.4 方法细节：模型构建

  • P-RQC 模型：采用一维砖墙式（brick-layer）几何结构的随机 Haar 幺正门。在每个时空点，以概率 $p$ 进行强力测量（投影到 $|0\rangle$ 态）。
  • RTN 模型：构建一个二维张量网络，其张量元素从高斯分布中独立同分布采样。通过调节边缘键维数（Bond Dimension）或耦合强度来驱动相变。
  • 量热指标：使用三体互信息（Tripartite Mutual Information, $I_3$）作为主要的探测手段。$I_3$ 在临界点具有尺度不变性，是确定 $p_c$ 的利器。

2. 关键基准体系，计算所得数据，性能数据

2.1 临界指数对比 (Table I 分析)

通过对不同系统的有限尺寸放缩（Scaling Collapse），论文给出了下表所示的核心数据（简化版）：

模型$p_c / \chi_c$$\nu$ (关联长度)$\eta$ (反常维度)$z$ (动力学)$c_{\text{eff}}$ (中心荷)
P-RQC0.24(1)2.1(4)0.17(2)0.96(5)-0.40(6)
RTN0.9(9)2.2(5)0.16(2)1.05(5)-0.35(2)
Standard MIPT0.168(5)1.2(2)0.19(1)1.06(4)0.25(3)
Percolation0.51.330.2110.2914

数据解读：

  1. $\nu$ 的漂移：P-RQC 的 $\nu \approx 2.1$ 远大于标准 MIPT 的 $1.2$。根据 Harris 判据，$\nu > 2/d$（此处 $d=1$），这意味着后选择相变对静态失序更加鲁棒。
  2. 负中心荷：这是该工作最令人惊讶的发现。通常 CFT 的中心荷 $c > 0$。负值暗示了非幺正 CFT 的存在，这与随机系统中的 replica 处理有关。
  3. 普适类一致性：P-RQC 与 RTN 在 $\nu, \eta, z, c_{\text{eff}}$ 上几乎完全吻合，严格证明了后选择相变在本质上就是张量网络的几何相变。

2.2 净化动力学（Purification Dynamics）

论文研究了将一个辅助比特（Ancilla）与系统最大纠缠后的纯化过程。在面积律相（Disentangling phase），辅助比特的熵 $S_a$ 随时间指数衰减；而在体积律相，纯化时间随系统尺寸指数增长。数值结果显示，$S_a$ 的放缩曲线在 $p_c$ 处展现出完美的重叠。

2.3 平移对称性的影响

实验探讨了如果门电路不再随机，而是固定的、平移对称的情况。结果发现：

  • 对于 Qubit (d=2) 系统，平移对称的后选择电路不表现出相变，弱测量强度 $p_w$ 趋于零时系统依然处于某种平庸态。
  • 对于 Qutrit (d=3) 系统，即便没有时空随机性，相变依然存在。这表明局部希尔伯特空间的维度是驱动此类非随机相变的“阈值资源”。

3. 代码实现细节，复现指南，所用的软件包

3.1 核心算法：张量收缩与矩阵乘积态 (MPS)

复现该工作的核心在于模拟非幺正演化下的纠缠熵。对于小尺寸系统（$L < 24$），建议使用全状态矢量演化；对于大尺寸，必须使用 MPS 方法。

复现步骤建议：

  1. Haar 随机门生成：
    • 使用 QR 分解：生成一个随机复高斯矩阵 $A$，进行 $A = QR$ 分解，得到的 $Q$ 即为 Haar 随机门。
  2. 演化算符应用：
    • 按照砖墙结构交替应用两比特门。
    • 测量算符应用：$|\psi\rangle \leftarrow (P_0 \otimes I) |\psi\rangle$。注意，每次测量后必须手动记录波函数的范数以计算自由能密度 $f(L)$。
  3. $I_3$ 计算：
    • 将系统分为四个相等的子区 A, B, C, D，分别计算 $S_A, S_{A\cup B}$ 等组合。

3.2 软件包推荐

  • ITensors (Julia)：处理随机张量网络和 MPS 演化的首选。其 random_itensor 和 apply 接口非常适合此类任务。
  • QuTiP (Python)：适合快速实现小尺寸的 Qubit 链模拟，但在处理 $L > 20$ 的非幺正演化时速度受限。
  • JAX (Python)：利用 GPU 加速大规模的随机矩阵运算，尤其是在需要大量样本平均（Disorder average）时效率极高。

3.3 开源资源链接（类似项目参考）

虽然论文未提供直接代码库，但可以参考以下社区资源：

4. 关键引用文献，以及对这项工作局限性的评论

4.1 关键引用文献

  1. Skinner et al. (2019) [1]：MIPT 的开创性论文，首次提出了测量驱动的纠缠转换。
  2. Li et al. (2018) [2]：独立发现了量子芝诺效应相关的纠缠相变。
  3. Vasseur et al. (2019) [15]：建立了 MIPT 与随机张量网络统计力学映射的理论框架。
  4. Gullans & Huse (2020) [5]：定义了通过辅助比特纯化来探测相变的方法。

4.2 工作局限性评论

  • 有限尺寸效应：尽管理论上 $\nu \approx 2.1$，但在 $L \le 20$ 的模拟中，放缩函数（Scaling functions）在 P-RQC 和 RTN 之间仍存在细微差异。这暗示可能存在更高阶的非普适修正，或者当前的系统尺寸尚未完全进入渐近区。
  • 计算复杂度：后选择虽然在理论上简化了统计映射（$k=0$），但在实验复现上仍然是“不可能的任务”。该研究更多是理论上的“存在性证明”，而非直接指导实验观测。
  • 负 $c_{\text{eff}}$ 的物理图景：虽然数值上得到了负中心荷，但对应的共形场论（CFT）的具体拉格朗日量尚不明确。这种非幺正 CFT 是否具有稳定的重整化群流仍有待研究。
  • 维度限制：目前的研究集中在一维系统（1+1D 时空）。在更高维度下，后选择是否会诱导出更复杂的拓扑序或不同的普适类，目前尚无结论。

5. 其他必要补充：各向异性因子与 CFT 修正

5.1 各向异性因子 $v$ 的计算

在 RQC 中，由于门电路在时间轴和空间轴上的不对称性（时间轴是单向演化，空间轴是双向耦合），系统表现出各向异性。为了正确拟合 Casimir 形式的自由能：

$$F(L, t) / A = f(L = \infty) - \frac{\pi c_{\text{eff}}}{6L^2}$$

必须先确定 $v$。论文通过比较空间关联函数与时间关联函数的衰减率（见 Fig. 6），计算出 $v_{\text{RQC}} \approx 0.70$。相比之下，RTN 在构造上是对称的，$v_{\text{RTN}} = 1$。这种细节处理显示了作者在临界数据分析上的严谨性。

5.2 为什么 Qutrit 很重要？

在平移对称电路中，$d=2$（Qubits）提供的状态空间太窄，导致测量操作太容易将系统推向单一的 product state。而 $d=3$（Qutrits）增加了状态的“冗余度”，使得系统能够维持一种竞争平衡。这一发现对于设计鲁棒的量子纠错码具有潜在启发：在没有随机性的情况下，增加局部编码维度是维持长程纠缠的一种策略。

5.3 结语

这项工作打通了“量子线路后选择”与“静态随机张量网络”之间的壁垒，为理解测量诱导的非平衡态相变提供了统一的几何视角。对于量子化学或量子物理研究者而言，这种从动态演化到静态网络映射的思维方式，是处理多体关联问题的有力武器。