直径截断算子演化

论文信息

标题: Diameter truncated operator evolution

作者: Tom Holden-Dye, Max Marvell, Joel Mills, et al.

发布日期: 2026-06-26

PDF 链接: 下载 PDF

研究背景与动机

量子多体系统的非平衡动力学是当代物理学最具挑战性的问题之一。当系统从初始状态开始演化时，局部可观测量通常会弛豫到热平衡，这一过程可以通过粗粒化的流体动力学理论来精确描述。然而，研究系统如何从初始状态逐渐趋近热平衡——即短时间和中等时间尺度的行为——却异常困难。

这种困难根源于量子复杂度的爆炸性增长。随着时间推移，非平衡量子态的纠缠熵快速增加，初始局域的算子会迅速变得高度非局域，“加扰” 到整个系统中。传统的张量网络方法无法有效捕捉完整的多体波函数动力学演化，因为所需的计算资源会随系统规模呈指数增长。

近年来，研究者们发现了一个重要的突破口：虽然算子整体变得极其复杂，但大部分算子的 “重量”（operator weight）对实际关心的局域物理量贡献甚微。特别是，在计算两点相关函数时，只需要保留少数关键的算子演化轨迹即可获得很好的近似。基于这一观察，以 dissipation-assisted operator evolution (DAOE) 和 Pauli Propagation 为代表的一类新方法应运而生，它们通过截断高权重的泡利串项来控制计算复杂度。

本文提出的直径截断算子演化（Diameter Truncated Operator Evolution, DTOE）方法在这一思想基础上另辟蹊径：不再依据算子的泡利重量进行截断，而是根据算子在实空间晶格上的直径——即算子非平凡支持域的空间跨度——来限制模拟的希尔伯特空间。这一选择不仅具有明确的物理动机，而且在计算效率上具有潜在优势。

核心方法：折叠表示与直径截断

DTOE 方法建立在 1+1 维砖墙型量子线路（brickwork circuit）框架之上，处理 $2L$ 个 $q$ 维量子比特组成的一维晶格 $\Lambda = \mathbb{Z}_{2L}$ 。系统的幺正演化算符为：

U(t) = \prod_{t'=1}^t U_{t'} \in U(q^{2L})

其中每一层 $U_{t'}$ 由相邻两比特门 $u \in U(q^2)$ 交错排列构成。

为了高效计算无穷温度下的两点相关函数 $C_{AB}(x,t) = \text{tr}(A^\dagger_x B_0(t) \rho_\infty)$ ，DTOE 使用了巧妙的 “折叠” 算子到态映射技术。具体而言，通过映射 $|a\rangle\langle b| \mapsto |a\rangle \otimes |b\rangle^*$ 将作用于原始希尔伯特空间的算子转化为双倍希尔伯特空间中的向量。折叠门 $v = u \otimes u^*$ 自然地给出了伴随作用 $u A u^\dagger$ 的向量表示。

这个方法的核心创新在于限制演化空间：只保留直径不超过 $d$ 的算子。形式上，定义受限算子空间：

\mathcal{V}_d = \bigoplus_{x \in \Lambda} [\mathcal{P}_d]_x

其中 $\mathcal{P}_d$ 表示以某个固定位置为左端点、长度不超过 $d$ 、且最左端泡利算子非平凡的泡利串集合。关键的是，固定直径 $d$ 的算子数量仅为 $O(L)$ ，而固定重量 $d$ 的算子数量则为 $O(L^d)$ ，二者的比例 $N_{\text{diameter-}d}/N_{\text{weight-}d} \to 0$ （在大 $L$ 极限下，对所有 $d \ge 2$ ）。这意味着直径截断原则上比基于重量的方法更经济。

DTOE 的动力学通过定义转移矩阵 $\mathbf{M}_d$ 来实现，该矩阵由四种子矩阵 $M_d^{0,0}$ 、 $M_d^{0,1}$ 、 $M_d^{1,0}$ 、 $M_d^{1,1}$ 组成，分别对应不同的边界条件组合。在每个时间步，初始向量按以下规则更新：

|B_0(t)\rangle_d = \prod_{t'=1}^t \bigoplus_{x: x\oplus t'=1} [\mathbf{M}_d]_{x,x+1} |B_0\rangle

相关函数则简单地由内积给出： $C_{AB}(x,t)_d = (A_x | B_0(t))_d$ 。

物理洞察：容易区与困难区

DTOE 的有效性建立在深层的物理洞察之上。参考文献[29]对混沌量子系统中两点函数的主导算子轨迹进行了分析，发现在时空晶格中存在一个 “容易” 相和一个 “困难” 相。在容易相中（通常是超出一某个临界速度 $v_c$ 的区域），对两点函数贡献最大的轨迹全程停留在 $O(1)$ 直径的算子空间内。这意味着使用常数值的截断直径 $d$ 即可精确近似这些两点函数，DTOE 的计算复杂度仅为 $O(t^2)$ （多项式级）。

对 1+1 维 Haar 随机线路，临界速度等于蝶形速度 $v_B = (q^2-1)/(q^2+1)$ （对 $q=2$ 的量子比特， $v_B = 3/5$ ）。但在远离 Haar 随机的系统中，容易区的范围可以大大扩展。受扰动对偶幺正线路（perturbed dual-unitary circuits）就是一个典型例子，本文在数值实验中充分验证了 DTOE 在这类系统中的有效性。

守恒荷的存在也会改变容易-困难的划分。对于与守恒密度相关的两点函数，主导轨迹通常局限在与守恒量关联的 “薄” 算子（ $O(1)$ 支持域）中。本文在海森堡 XXZ 模型中验证了这一点——该模型具有 $U(1)$ 总自旋 $S^z$ 守恒——即使在各项同性点（ $\Delta=1$ ），DTOE 依然成功提取出了超扩散输运标度指数。

数值实验与结果分析

受扰动对偶幺正线路：踢动伊辛模型

踢动伊辛模型（kicked Ising model）由以下 Floquet 算符描述：

U_{\text{KI}}[J,B,h] = e^{-iJ Z\otimes Z - i\frac{h}{2}(Z\otimes \mathbf{1} + \mathbf{1}\otimes Z)} \cdot e^{-iB(X\otimes \mathbf{1} + \mathbf{1}\otimes X)} \cdot e^{-iJ Z\otimes Z - i\frac{h}{2}(Z\otimes \mathbf{1} + \mathbf{1}\otimes Z)}

当 $J=B=\pi/4$ 时该模型具有对偶幺正性，两点函数严格限制在光锥边缘。引入扰动 $\eta$ （ $J=B=\pi/4+\eta$ ）后，DTOE 在 $\eta=0.24588$ （远离微扰区）的测试中表现出色。对于 $x = t-21$ 的相关函数 $|C_{XX}|$ ，低直径（ $d=1,2$ ）数据明显偏离真实值，而 $d \ge 8$ 的结果在 $t$ 达到 100 层的长时间窗口内保持 $10^{-2}$ 至 $10^{-3}$ 量级的相对误差。对不同 $\eta \in [0,0.5]$ 的系统性扫描表明，DTOE 的中位相对误差普遍低于 $10^{-2}$ ，只在 $\eta \approx 0.1$ 附近有一个小的性能下降窗口。

海森堡 XXZ 模型：输运标度指数提取

XXZ 线路的门结构为 $U(J,J') = e^{-i[J(X\otimes X + Y\otimes Y) + J'(Z\otimes Z)]}$ 。根据各向异性参数 $\Delta = J'/J$ 的值，系统的自旋输运呈现三种截然不同的标度行为：

$|\Delta| < 1$ ：弹道输运， $\alpha = 1$
$|\Delta| = 1$ （XXX）：超扩散， $\alpha = 2/3$
$|\Delta| > 1$ ：扩散， $\alpha = 1/2$

DTOE 通过计算自旋关联函数的展宽 $\sigma(t) \sim t^\alpha$ 来提取 $\alpha$ 。在弹道区（ $J=0.6501,J'=0.4032,\Delta\approx0.62$ ）， $d=10$ 的数据给出 $\alpha = 1.00 \pm 0.01$ ，与理论预测完美吻合。有趣的是，在极长时间后数据会偏离并趋向 $\alpha=1/2$ ，这反映出有限截断直径不足以捕捉全部准局域守恒荷的长程效应——正是这些长串的贡献保证了真正的弹道输运。DTOE 能保持对正确指数的时标 $t^*$ 与直径呈指数关系（如图 4 所示，通过偶奇直径数据分别拟合得出 $t^* = O(\exp(d))$ ）。

在超扩散区（ $J=J'=0.1163,\Delta=1$ ），不同拟合窗口（ $\delta \in [10,50]$ ）从 $d=10$ 数据提取的 $\alpha$ 落在 $[0.6663, 0.6684]$ 区间，与理论值 $2/3$ 高度吻合。值得强调的是， $\Delta \to 1$ 时所有尺度的磁振子束缚态（准粒子串）都参与输运并对超扩散指数至关重要，这使得该点对 DTOE 尤为困难，但小耦合强度下的结果依然可靠。

在扩散区（ $J=0.4169,J'=0.7281,\Delta\approx1.75$ ），即使是 $d=1$ 的最低截断也能准确复制 $\alpha=1/2$ 。这是因为扩散输运由短关联长度的准粒子主导，直径截断自然保有全部相关知识。

应用前景与未来方向

DTOE 的高效性——在常数 $d$ 下时间复杂度约 $O(t^2)$ ——使其成为在普通计算资源上模拟量子多体动力学的有力工具。本研究所有结果均未使用高性能计算集群，却能达到当前量子计算机难以企及的长时间尺度（100 层以上，而当前量子实验通常只能实现约 25 个 Floquet 周期）。

该方法最具前景的推广是向更高维度的拓展。参考文献[29]论证了在 2D 情形下，Haar 随机线路的整个时空晶格都可能落入容易相，两点函数由穿越 $O(1)$ 直径算子空间的轨迹主导。此时固定直径的算子数量按 $O(L^D)$ （ $D$ 为空间维度）标度，在低维情况下依然保持高效，为研究二维系统的输运性质提供了有希望的途径。

此外，DAOE 框架中使用耗散超算符实现软截断（施加权重衰减而非硬截断）的策略对 DTOE 也有借鉴意义。本文附录 A 给出了构造基于直径的耗散超算符的精确 MPO 表示，其键维数为 $d+1$ （与 DAOE 的 $l^*+1$ 完全对应）。实现这种软截断版本将允许对耗散强度 $\gamma$ 进行外推，从而进一步提升精度。这部分留作未来工作。

总结

直径截断算子演化（DTOE）方法为模拟非平衡量子多体系统的算子动力学提供了一条简洁而有效的新路径。通过基于实空间直径（而非泡利重量）裁剪希尔伯特空间，该方法在保持核心物理内容的同时大幅降低了计算复杂度。对受扰动对偶幺正线路和 XXZ 模型的广泛数值实验证实，DTOE 能在从弹道到扩散的多种输运区精确提取两点相关函数和流体动力学标度指数，时间尺度远超当前量子仿真能力。其向高维系统的潜在推广以及与软耗散外推技术的结合，预示着该方法将成为连接经典模拟与量子物理前沿的重要桥梁。