来源论文: https://arxiv.org/abs/2604.23192v1 生成时间: May 02, 2026 12:03

0. 执行摘要

在本研究中，科研团队系统性地探讨了全局守恒律（Conservation Laws）如何塑造多体量子系统中的相干性（Coherence）展宽过程。研究涵盖了三类典型动力学模型：U(1)对称随机线路、电荷-偶极矩守恒线路（分形子动力学）以及遍历性哈密顿量动力学（混合场伊辛模型）。

研究通过全局探针（参与熵 $S_d$）和局部探针（相干性相对熵 $C_d$）揭示了守恒律的深远影响。主要结论包括：

  1. 全局演化： 守恒律打破了无约束随机线路中常见的对数级快速饱和，代之以受流体动力学制约的幂律弛豫。在有限尺寸系统中，这一过程表现为从幂律衰减到指数衰减的交叉演化，且弛豫时间随系统尺寸呈代数增长。
  2. 局部演化： 对称约束的随机线路展现出清晰的“上升-峰值-下降”剖面。峰值时间 $\tau_c^m$ 随子系统尺寸 $L_A$ 呈代数缩放。相比之下，遍历性哈密顿量在较大子系统时会将尖锐的峰展宽为长延展的平台，暗示了本质上不同的局部弛豫机制。
  3. 理论解释： 研究结合了复本张量网络（RTN）、矩阵乘积态（MPS）模拟及解析分析（如SSEP有效模型和稀疏区域分析），证明了相干性可以作为对称约束热化过程的敏感探针，直接联系了量子资源动力学与多体输运现象。

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

1.1 核心科学问题

量子动力学的核心问题之一是系统如何失去对初始状态的记忆。在没有守恒律的系统（如Haar随机演化）中，信息展宽极快，相干性迅速耗散。然而，当物理定律要求某些物理量（如电荷、动量或偶极矩）守恒时，这些“慢模”会成为演化的瓶颈。本研究聚焦于：守恒律究竟如何在定量和定性上改变相干性的消失速度？ 这种约束在全局和局部层面分别有何表现？

1.2 理论基础：相干性的资源理论

研究采用了相干性资源理论（Resource Theory of Coherence）的视角。在给定的参考基底（通常是计算基底）下，相干性代表了态处于基矢叠加的能力。

  • 参与熵（Participation Entropy, $S_k$）： 全局测度，定义为 $S_k(|\Psi\rangle) = \frac{1}{1-k} \log_2 \sum_x p_x^k$，其中 $p_x = |\langle x|\Psi\rangle|^2$。本文主要研究 Renyi-2 参与熵 ($k=2$)，它直接关联到逆参与比（IPR）和碰撞概率。$S_d$ 测量波函数在希尔伯特空间中的非集中程度。
  • 相干性相对熵（Relative Entropy of Coherence, $C_d$）： 局部测度，定义为 $C_d(\rho_A) = S_d(\rho_A) - S_R(\rho_A)$。这里 $S_d$ 是对角熵，$S_R$ 是冯·诺依曼（或Renyi-2）纠缠熵。$C_d$ 量化了子系统在扣除掉与外界纠缠贡献后，自身保留的纯相干性。

1.3 技术难点

  1. 希尔伯特空间爆炸： 随着自旋数 $L$ 增加，基矢数量呈指数级 $q^L$ 增长，常规精确对角化（ED）难以处理 $L > 20$ 的系统。
  2. 随机线路平均： 随机线路的结果依赖于具体实现，需要进行大量的样本平均以获取统计特性。
  3. 守恒律导致的希尔伯特空间破碎（Fragmentation）： 尤其在偶极矩守恒系统中，系统会分裂成无数相互独立的Krylov子空间，导致动力学极其复杂，需要针对特定片段（Fragment）进行解析和模拟。

1.4 方法细节

  • 复本张量网络（RTN）： 用于模拟大尺寸系统的 Renyi-2 量。通过将算符映射到双倍希尔伯特空间中的向量，将电路平均转化为局域转移矩阵（Transfer Operator）的计算。对于 $q=2$ 的情况，转移矩阵为 $36 \times 36$；对于 $q=3$，则达到 $225 \times 225$。这允许模拟高达 $L=128$ 的系统。
  • 矩阵乘积态与TDVP： 对于哈密顿量动力学，使用时间依赖变分原理（TDVP）结合MPS。针对 $S_d$ 的计算，研究团队开发了一种创新的算法：通过构造概率分布向量的MPS表示，并将其范数作为对角纯度。这种方法避免了对 $q^{L_A}$ 个项的暴力求和，将计算复杂度从指数级降至 $O(L_A d \chi_p \chi^3)$。
  • 对称简易排斥过程（SSEP）映射： 利用稀疏区域理论，将 $U(1)$ 对称动力学映射为两个独立的扩散模。解析推导出全局相干性偏离值 $\Delta S_d(t) \sim t^{-1}$ 的演化律。

2. 关键 Benchmark 体系与计算数据

2.1 U(1) 对称随机线路 (Spin-1/2 & Spin-1)

  • 初始态： 采用不相干的 Néel 态（如 $|+1, -1, +1, -1, \dots\rangle$）。
  • 全局数据：
    • 偏离值 $\Delta S_d(t) = S_d(\infty) - S_d(t)$ 在中间时间段表现为明显的幂律衰减 $t^{-\beta_p}$，其中 $\beta_p \approx 1.04$。
    • 弛豫时间 $\tau_p$ 随 $L$ 呈代数增长，$\tau_p \sim L^{1.95}$，符合电荷扩散的流体动力学特征。
  • 局部数据：
    • 相干性峰值时间 $\tau_c^m \sim L_A^{0.52}$。这源于 $S_d$ 的代数上升与 $S_R$ 的线性增长（弹道增长）之间的竞争。

2.2 电荷与偶极矩守恒线路 (Fractonic Circuits)

  • 特征： 极其严重的希尔伯特空间破碎。演化被局限在特定的 Krylov 片段中。
  • 数据：
    • 全局 $\Delta S_d(t)$ 同样呈现幂律衰减，但指数 $\beta_p \approx 1.91$，明显大于 $U(1)$ 情况，反映了偶极矩约束下更复杂的输运模式。
    • 局部峰值时间缩放指数 $\alpha_m \approx 0.313$，显示出更快的相干性耗散速度（相对于 $U(1)$）。

2.3 混合场伊辛模型 (MFIM)

  • 哈密顿量参数： $b = (5+\sqrt{5})/8$, $h_i = (\sqrt{5}+1)/4$。这种配置确保了系统的非积性与量子遍历性。
  • 关键发现：
    • 峰值到平台的交叉： 在小型子系统 $L_A=2, 3$ 时，相干性呈现清晰的峰值；但随着 $L_A$ 增大，峰值展宽并形成稳定的平台。这说明哈密顿量动力学中的局部相干性具有极强的鲁棒性，与随机线路的快速衰减形成鲜明对比。
    • 全局 $\Delta S_d(t) \sim t^{-0.83}$，局部 $\Delta S_d(t) \sim t^{-0.45}$，证明即便在遍历哈密顿量中，慢扩散模依然主导后期弛豫。

3. 代码实现细节与复现指南

3.1 环境要求与软件包

复现本工作的数值结果推荐以下技术栈：

  • 核心框架： C++ 开发，利用多线程加速随机样本采样。
  • 张量网络库： ITensor (用于 MPS, TDVP) 或 Quimb (用于高效的张量收缩)。
  • 矩阵运算： Eigen 或 Armadillo。
  • 精确对角化： QuSpin 库可用于快速构建 $U(1)$ 和偶极矩守恒的哈密顿量基底。

3.2 复现关键步骤（以RTN为例）

  1. 构造转移算符： 根据论文附录 B，根据 Weingarten 函数或 Haar 测度推导 $n=2$ 的局部转移矩阵 $T_{i,i+1}$。注意区分不同电荷扇区 $Z_1, Z_2$ 的配对方式。
  2. 构建基底： 对于 $q=2$，定义六个局部单点复本态（见 Eq. B10）。
  3. 电路收缩： 模拟砖墙结构（Brick-wall circuit）。每层演化相当于在 MPS 上作用一层 $T$ 算符。为了保持 bond dimension $\chi$ 可控，必须在每步引入 SVD 截断（通常保留前 800-1024 个奇异值）。
  4. 观测量提取： 参与熵通过边界向量 $\langle\langle \Gamma_A |$ 与演化后的复本态 $\mathcal{T}^t |\rho_0^{\otimes 2}\rangle\rangle$ 的内积获得。

3.3 开源资源链接

  • ITensor Repo: https://github.com/ITensor/ITensor (推荐用于哈密顿量模拟)
  • 作者提供的代码（拟定）： 论文提到 “Code and data will be publicly shared at publication”，通常会托管在 GitHub 组织的个人主页（如 https://github.com/sreemayee-aditya/coherence-dynamics）。

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

4.1 关键引用

  1. [65] Baumgratz et al., Phys. Rev. Lett. 113, 140401 (2014): 定义了量子相干性资源理论的公理化框架。
  2. [16] von Keyserlingk et al., Phys. Rev. X 8, 021013 (2018): 奠定了带守恒律随机算符动力学（算符氢动力学）的基础。
  3. [58] Rakovszky et al., Phys. Rev. Lett. 122, 250602 (2019): 揭示了扩散导致的纠缠熵亚线性增长机制。
  4. [61] Turkeshi et al., Phys. Rev. B 102, 064305 (2020): 开发了复本张量网络模拟量子线路中 Renyi 指数的方法。

4.2 工作局限性评价

  • 维度限制： 本工作完全基于一维系统。在高维系统中，扩散常数和流体动力学长尾会发生显著变化，相干性的幂律指数是否依然稳健尚存疑问。
  • 复本指数 $n$： 研究主要集中在 $n=2$（Renyi-2），虽然这在数值上易于处理且实验可测，但冯·诺依曼极限 ($n \to 1$) 的某些精细行为可能被掩盖。
  • 破碎效应的非普适性： 论文提到分形子线路的结果高度依赖于选定的 Krylov 片段。这意味着在真实的物理系统中，不同的初态可能会导致截然不同的相干性寿命。

5. 补充解析：稀疏区域效应与初态选择

5.1 稀疏区域（Rare Regions）的物理图景

这是理解本文 $S_R$ 与 $S_d$ 竞争的关键。在带电荷守恒的线路中，如果某个区域的所有点电荷都处于极值态（全向上或全向下），该区域会被“冻结”，导致电荷传输无法跨越。只有当外界的波动扩散到该区域并将其“激活”时，纠缠才能继续增长。这就是为什么 $S_R$ 在存在守恒律时会变慢的原因。然而，本工作选取的 Néel 态本身含有极高密度的移动电荷，有效抑制了稀疏区域的形成，使得 $S_R$ 在初期仍保持弹道增长。这解释了为何相干性峰值出现在极短的时间内。

5.2 对角熵 $S_d$ 的扩散模型

附录 C 中详细讨论了 $\Delta S_d(t) \sim t^{-1}$ 的起源。在 $n=2$ 复本形式下，对角纯度的演化可以看作是两个复本中的粒子相互独立扩散的过程。由于每个复本贡献一个 $t^{-1/2}$ 的衰减因子，两者乘积自然导出了 $t^{-1}$。这一简洁的解析关系与数值模拟中 $\beta_p \approx 1.04$ 的结果惊人吻合，展现了统计力学模型在量子多体问题中的强大预测力。

5.3 实验可观测性

论文特别强调了 Renyi-2 观测量在冷原子和超导量子比特平台上的可实现性。利用影子层析成像（Shadow Tomography）或随机测量方案，可以直接测量对角纯度 $P_{A,diag}$，这意味着本研究提出的“相干性峰”可以作为验证量子处理器是否遵循对称约束热化的实验标尺。