来源论文: https://arxiv.org/abs/2602.22331v1 生成时间: Feb 27, 2026 02:14

非相互作用系统中的算符 Rényi 熵与纠缠增长：基于 Schwinger-Keldysh 场论的深度解析

0. 执行摘要

量子多体系统中的动力学演化一直是凝聚态物理和高能物理研究的核心。近年来，算符增长（Operator Growth）与纠缠熵（Entanglement Entropy）的生成被认为是理解量子信息传播、热化（Thermalization）以及输运性质的关键纽带。本文深度解析了 Priesh Roy 与 Sumilan Banerjee 的最新研究成果。该工作提出了一种衡量算符增长的新度量——子系统算符 Rényi 熵（Subsystem Operator Rényi Entropy, $S_{OA}^{(2)}$），并构建了一个统一的 Schwinger-Keldysh (SK) 场论框架，用以计算非相互作用系统中的算符和状态纠缠增长。

研究表明，该方法能够精确捕捉包括一维/二维 Aubry-André (AA) 模型和二维 Anderson 模型在内的多种体系中的输运特性。关键结论包括：

  1. 算符 Rényi 熵提供了一个状态无关（State-independent）的探测手段，能直接反映守恒荷的输运行为。
  2. 通过 SK 场论导出的关联矩阵方法，可以高效计算大尺寸系统的演化。
  3. 研究揭示了弹道、扩散、亚扩散（Sub-diffusive）以及超弹道（Super-ballistic）增长模式与系统输运指数之间的严格对应关系。
  4. 非相互作用系统中的 Schmidt 值分布呈现出与相互作用系统截然不同的均匀性，反映了其独特的准粒子动力学特征。

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

1.1 核心科学问题：量子信息传播与输运的脱钩与耦合

在传统的量子多体物理中，输运性质通常通过线性响应理论（如 Kubo 公式）或关联函数来定义。然而，在量子信息领域，人们更倾向于使用量子纠缠增长和算符演化（如 Out-of-time-order correlators, OTOC）来描述动力学。长期以来，物理学界试图寻找这两种描述之间的统一桥梁：纠缠的增长速度是否受限于物理量的输运速度？

以往的研究发现：

  • 在扩散性相互作用系统中，$n>1$ 的 Rényi 熵 $S_A^{(n)}(t) \sim \sqrt{t}$，而 von Neumann 熵 $S_A(t) \sim t$。
  • 在可积系统中，所有熵通常呈弹道式线性增长 $S \sim t$。
  • 在多体定位（MBL）系统中，纠缠呈对数增长。

然而，对于具有**次弹道输运（Sub-ballistic transport）**特性的非相互作用系统（如准晶模型、Anderson 定位边缘），算符增长的细节仍不明朗。本文的核心使命在于定义一种新的度量，并证明其在这些复杂输运背景下的鲁棒性。

1.2 理论基础：算符 Rényi 熵的定义

作者定义了子系统算符 Rényi 熵 $S_{OA}^{(2)}(t)$。考虑一个局部算符 $\mathcal{O}$（如某点处的费米子数算符 $\hat{n}_l$），随时间演化为 $\mathcal{O}(t) = e^{iHt}\mathcal{O}e^{-iHt}$。定义一个基于部分迹（Partial Trace）的“等效密度矩阵”：

$$\rho_{OA}(t) = \text{Tr}_B[\mathcal{O}(t)] / Z_{\mathcal{O}}$$

其中 $Z_{\mathcal{O}} = ext{Tr}[\mathcal{O}]$。由此定义的第二 Rényi 熵为：

$$e^{-S_{OA}^{(2)}(t)} = ext{Tr}_A [ ho_{OA}^2(t)]$$

与传统的算符纠缠熵（Operator Entanglement Entropy）不同，此定义不需要构造双倍希尔伯特空间（Doubled Hilbert space），且包含了丰富的空间位置信息，使其能够直接映射到局域物理量的扩散过程。

1.3 技术难点：Schwinger-Keldysh 路径积分的构建

计算 $ ext{Tr}_A[( ext{Tr}_B \mathcal{O})^2]$ 是极其困难的，因为它涉及到部分迹的平方。传统方法（如精确对角化 ED）受限于希尔伯特空间的维度（$2^L$），难以处理超过 20 个格点的系统。技术突破在于：

  1. 引入费米子位移算符（Fermionic Displacement Operators）：利用等式 $ ext{Tr}_A[( ext{Tr}_B O)^2] = \int f_A(\xi_1, \xi_2) \prod_{\alpha} d^2\xi_\alpha ext{Tr}[O D_A(\xi_\alpha)]$，将部分迹的平方转化为费米子 Wigner 特征函数的乘积。
  2. Keldysh 闭时轮廓（Closed-time Contour）：算符 $\mathcal{O}(t)$ 的演化涉及 $e^{-iHt}$ 和 $e^{iHt}$，这在 SK 场论中自然对应于向前（+）和向后（-）分支。作者通过对 Grassmann 场进行路径积分，成功导出了关联矩阵 $\mathbb{C}_O(t)$ 的闭合表达式。

1.4 方法细节：非相互作用系统的简化

对于费米子非相互作用系统，二次型哈密顿量使得路径积分可以解析积出。最终，熵的计算被简化为关联矩阵的行列式：

$$S_{OA}^{(2)}(t) = - ext{tr} \ln [(I - \mathbb{C}_O(t))^2 + \mathbb{C}_O^2(t)]$$

对于初态为纯积态（Product State）的情况，类似的推导给出了状态纠缠熵的公式。这使得计算复杂度从指数级降到了多项式级 $O(L^3)$，从而允许模拟数百个格点的长时间演化。

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

2.1 一维 Aubry-André (AA) 模型

AA 模型具有经典的定位-非定位相变。当准周期势能强度 $V < 1$ 时为金属相， $V > 1$ 时为定位相， $V = 1$ 为临界点。

  • 饱和时间缩放（Saturation Time Scaling）： 作者定义 $t_{op}$ 为熵达到饱和值一半所需的时间。在金属相，$t_{op} \sim L^1$（弹道输运）；在临界点 $V=1$，$t_{op} \sim L^2$（扩散输运）。 数据点：在 $V=1$ 时，拟合得到的指数 $\alpha_{op} \approx 2.0$。这完美契合了该模型在临界点处的亚扩散性质。

  • 增长曲线分析： 算符 Rényi 熵在达到饱和前表现出 $\ln t$ 的对数增长特征。与之对应，状态 von Neumann 熵 $S^{(1)}$ 在临界点表现为 $\sqrt{t}$ 增长（扩散性），而在金属相表现为 $t$ 的线性增长（弹道性）。

2.2 二维 Aubry-André 模型与超弹道现象

二维 AA 模型（无次近邻跳跃 $t'=0$）在 $V < 1$ 时表现出有趣的**异常超弹道（Super-ballistic）**行为。

  • 异常指数：在 $V=0.25$ 时，饱和时间的缩放指数 $\alpha_{vN} \approx 0.6$，这意味着纠缠传播速度超过了线性。作者通过引入次近邻跳跃 $t'=1/3$ 破坏了模型的特殊对称性（van-Hove 奇异性），使得系统恢复了预期的弹道（$\alpha \approx 1$）或亚扩散（$\alpha \approx 2$）行为。

2.3 二维 Anderson 模型

研究了在弱定位 regime（格点尺寸 $L \ll \xi$ 定位长度）下的扩散行为。

  • 性能数据：通过对 100 个随机实现取平均，计算得到 $\alpha_{op} \approx 2.2$。这表明即使在定位效应尚未显现的尺度下，算符 Rényi 熵也能精准探测到由随机势引起的扩散延迟。

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

3.1 算法流程

要复现本文结果，研究者应遵循以下数值计算步骤：

  1. 构造单粒子格林函数： 计算单粒子哈密顿量的本征值 $\epsilon_\alpha$ 和本征态 $\phi_\alpha(i)$。构造 Keldysh 轮廓上的格林函数 $G_{ij}(t, t')$：

    $$iG^<_{ij}(t, t') = \sum_\alpha \phi_\alpha(i)\phi^*_\alpha(j) [-n_F(\epsilon_\alpha) e^{-i\epsilon_\alpha(t-t')}]$$

    对于算符增长，使用无限温分布 $n_F = 1/2$；对于状态增长，使用真空或特定填充分布。

  2. 构建关联矩阵 $\mathbb{C}(t)$： 利用文中公式 (15) 或 (23)。注意算符位移带来的 “kick” 项处理。矩阵维度为子系统大小 $L_A imes L_A$。

  3. 计算熵值： 对 $\mathbb{C}(t)$ 进行奇异值分解（SVD）或直接计算迹对数。对于 von Neumann 熵，使用公式：

    $$S_A(t) = - ext{tr} [(I - \mathbb{C}) \ln(I - \mathbb{C}) + \mathbb{C} \ln \mathbb{C}]$$

3.2 软件包建议

  • 编程语言：推荐使用 Julia (因其高性能线性代数库和对量子算符的良好支持) 或 Python (with NumPy/SciPy)。
  • 并行化：由于需要进行大量的随机势实现（Disorder Realizations）取平均，建议使用 Distributed.jl 或 multiprocessing 在集群上运行。
  • 核心库：LinearAlgebra (用于 eigen 和 logdet)。

3.3 开源资源参考

虽然作者未提供官方 repo，但此类计算可参考以下开源框架：

  • QuantumOptics.jl: 可用于构建基础哈密顿量。
  • ITensors.jl: 虽然本文主要讨论非相互作用系统，但 ITensor 是对比相互作用系统纠缠增长的标准工具。

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

4.1 关键引用文献

  1. Banerjee et al. [38, 40]: 奠定了 SK 场论计算纠缠熵的理论基石。
  2. Maldacena et al. [5]: 建立了算符增长与混沌（Chaos）之间的普适界限。
  3. Rakovszky et al. [10]: 探讨了扩散系统中的 Rényi 熵增长，是本文的重要对比对象。
  4. Aubry & André [34]: AA 模型的原始文献。

4.2 工作局限性评价

尽管该工作在理论和数值上非常出色，但仍存在以下局限：

  • 非相互作用限制：目前的结果主要集中在自由费米子。虽然场论框架具有普适性，但一旦引入相互作用，Grassmann 积分将不再是二次型，必须引入 DMFT 或 RPA 等近似方法。这会导致计算复杂度大幅上升且引入系统误差。
  • 算符依赖性：本文主要研究粒子数算符 $\hat{n}_l$。对于非守恒荷的算符增长，其与输运的耦合关系可能不如文中展现的那么直接。
  • Schmidt 值的均匀性：作者指出非相互作用系统中所有 Schmidt 值表现一致，但这可能是自由粒子体系的一个“人工产物”。在现实的电子关联材料中，Schmidt 谱的“尾巴”通常包含关键的物理信息。

5. 其他补充：物理直觉与深度思考

5.1 准粒子图像（Quasiparticle Picture）的失效与修正

在弹道系统中，纠缠增长可以用“纠缠对”向外传播的图像完美解释。但在亚扩散或定位系统中，准粒子图像需要修正。本文的研究表明，在临界 AA 模型中，由于单粒子波函数的**多重分形（Multifractality）**特性，信息并不是由单一速度的准粒子携带，而是经历了一种能量依赖的扩散过程。这在算符 Rényi 熵的对数增长阶段得到了体现。

5.2 对量子计算的启示

理解算符增长的界限对于设计容错量子计算机至关重要。本文提出的 $S_{OA}^{(2)}$ 可以作为评估量子门电路中噪声传播和错误扩散的灵敏度量。如果我们能控制算符增长处于亚扩散状态，可能有助于延长系统的量子相干时间。

5.3 总结视角

本项研究不仅提供了一个强大的计算框架（SK 场论），更重要的是，它将量子信息理论中抽象的“熵增长”具象化为凝聚态物理中熟悉的“输运过程”。对于量子化学家而言，这种将复杂动力学演化转化为格林函数计算的思路，在处理大分子激子动力学或电子转移过程时具有重要的借鉴价值。