来源论文: https://arxiv.org/abs/2602.15933v1 生成时间: Feb 21, 2026 02:03

邻近集成性的长程相互作用自旋链中的 KPZ 类输运解析

0. 执行摘要

在一维量子多体系统中，输运性质的分类是统计物理的核心课题。通常认为，通用的非集成系统在高温下表现为扩散输运，而具有特殊对称性的集成模型（如 XXX Heisenberg 链）则表现出属于 Kardar-Parisi-Zhang (KPZ) 普适类的超扩散行为（动态临界指数 $z=3/2$）。然而，本研究（Anand et al., 2026）通过最先进的张量网络数值模拟发现，即便是在非集成的长程幂律相互作用模型中，KPZ 类超扩散依然能极其顽强地存在于极长的时间尺度内。该发现挑战了“非集成必扩散”的直觉，并提出了“邻近集成性”（Proximity to Integrability）的物理图景：这些长程模型可以被视为对 Inozemtsev 集成模型族的微扰。这一结论对于利用 Rydberg 原子阵列和极性分子开展量子模拟实验具有重要的指导意义。

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

1.1 核心科学问题：输运的普适性边界

本工作的核心在于探索：长程相互作用是否会立即破坏 KPZ 类的动力学规律？ 在最近邻（NN）Heisenberg 链中， spin 输运已被证实属于 KPZ 普适类。但在 Haldane-Shastry (HS) 模型（$1/r^2$ 相互作用）中，输运是弹性的（ballistic, $z=1$）。这引出了一个关键矛盾：介于最近邻（$\alpha o \infty$）和 HS（$\alpha = 2$）之间的幂律模型（$1/r^\alpha$）是非集成的，理应在宏观时间下表现为扩散（$z=2$），但实验和数值观测往往在可观测时间内看到非扩散行为。本文旨在厘清这些行为背后的微观机制。

1.2 理论基础：Inozemtsev 模型与 KPZ 缩放

  • Inozemtsev 模型：这是一个关键的集成族，其相互作用强度定义为 $f_{Inoz}(r) = \sinh^2(\kappa)/\sinh^2(\kappa r)$。当 $\kappa o 0$ 时，退化为 HS 模型（弹道输运）；当 $\kappa o \infty$ 时，退化为 NN Heisenberg 模型（KPZ 输运）。
  • KPZ 动力学：特征在于电荷转移 $P_Q(t) \sim t^{1/z}$。对于 KPZ 类，$z=3/2$。其双点关联函数应符合特定的缩放函数 $\langle q_r(t) q_0(0) angle = bt^{-1/z} f(br/t^{1/z})$，其中 $f$ 为 KPZ 缩放函数。

1.3 技术难点：长程作用下的混合态演化

  1. 混合态表征：在无限高温下研究输运需要处理密度矩阵 $ ho$。为了利用张量网络，必须对密度矩阵进行“向量化”（Vectorization），将其转变为加倍希尔伯特空间（Doubled Hilbert Space）中的纯态 $| ho angle angle$。
  2. 长程 MPO 构建：幂律相互作用 $1/r^\alpha$ 无法直接用有限键维度的矩阵乘积算符（MPO）精确表示。传统做法是将其拟合为多个指数函数之和 $f(r) \approx \sum_i c_i e^{-\lambda_i r}$。这引入了模型近似误差与计算复杂度的权衡。
  3. 计算复杂度：长程作用意味着 MPO 的键维度显著增加，且演化算符不再是局域的，传统的 TEBD 算法失效，必须使用时间相关变分原理（TDVP）。

1.4 方法细节：加倍空间 TDVP 与弱对称性保护

  • 向量化演化：演化算符变为 $U \otimes U^*$。构建加倍空间的哈密顿量 $H^D = H \otimes I - I \otimes H^*$。本研究通过复杂的代数推导（见补充材料 SM2），给出了长程 XXZ 模型下自旋流和能量流的算符表达式。
  • 弱 $S^z$ 对称性：尽管初态是混合态，但在加倍空间中通过保守 $Q^D = Q \otimes I - I \otimes Q$ 电荷，可以获得至少 5 倍的计算加速。这使得模拟 $L=2000$ 个格点、时间长达 $t=1000/J$ 成为可能。
  • 初态制备：采用弱极化领域墙态 $ ho_\lambda \propto e^{\lambda \hat{R}}$，其中 $\hat{R}$ 是左右电荷差。这种方法避免了传统线性响应理论在大尺寸下的数值噪声。

2. 关键 Benchmark 体系，计算所得数据与性能数据分析

2.1 幂律 Heisenberg 模型的输运指数分析

论文考察了 $\alpha \in [2.2, 6.0]$ 的模型。

  • 发现 1：对于所有 $\alpha > 2$，自旋输运均表现出向 $z_S = 3/2$ 的收敛倾向。即使是强非集成的 $\alpha=3$ 模型，在 $t=400/J$ 范围内也未观察到向 $z=2$（扩散）的交叉。
  • 数据对比：在 $t=400$ 时，$\alpha=4.5$ 的自旋关联函数与 KPZ 缩放函数的拟合误差比 Gaussian 拟合低一个数量级。具体而言，相对误差 $\epsilon$ 在 $t > 200$ 后持续下降，证明了 KPZ 描述的稳健性。

2.2 Inozemtsev 模型作为“锚点”

研究揭示了 Inozemtsev 模型在整个 $\kappa$ 范围内（除 $\kappa=0$ 极端情况外）都呈现 KPZ 输运。

  • 关键数据：对于 $\kappa=0.4$ 的 Inozemtsev 模型，模拟得到的动态指数 $z_S \approx 1.503$，极度接近理想值 $1.5$。
  • 能量输运：与自旋不同，能量输运在所有考察的模型中均为弹道的（$z_E = 1$）。这反映了能量流作为准守恒荷的特殊地位。

2.3 实验相关模型：Rydberg 与极性分子

  • 极性分子模型 ($H_{pm}$)：$1/r^3$ 各向异性。发现随着各向异性参数 $\Delta$ 的增加，系统从弹道输运（$\Delta < 0.5$）转向超慢扩散甚至亚扩散（$\Delta > 2$）。
  • Rydberg 原子模型 ($H_{ryd}$)：$1/r^3$ 交换项与 $1/r^6$ Ising 项混合。数据显示其能量输运呈现显著的超扩散性质，即使自旋输运已表现为类扩散。这为实验观测非平凡流体动力学提供了直接目标。

2.4 数值收敛性 Benchmark (SM1/S4)

  • 键维度（Bond Dimension）：模拟中使用的 $\chi$ 最高达到 384。通过对 $\chi=128, 192, 256, 384$ 的对比，确认了自旋转移 $P_S$ 的收敛精度在 1% 以内。
  • SVD 截断：设置阈值为 $10^{-6}$，在长程演化中有效平衡了精度与计算资源。

3. 代码实现细节，复现指南与开源工具

3.1 核心软件包：TeNPy

该研究主要基于 TeNPy (Tensor Network Python) 库实现。TeNPy 因其高效的对称性处理和矩阵运算优化，已成为量子多体动力学领域的标准工具。

3.2 关键算法复现逻辑

复现本工作的关键步骤如下：

  1. 构建指数求和算符：利用 tenpy.networks.mpo.MPO.from_exp_sums 方法。建议针对 $\alpha=3$ 使用 $K=4$ 个指数项进行拟合，这足以在 $t < 600$ 的时间内保持极低误差（$\epsilon < 10^{-4}$）。
  2. 加倍希尔伯特空间定义：
    • 定义 local site，每个 site 有 4 个能级（对应原来的算符空间）。
    • 手动构建 $H^D$。注意项之间的符号，尤其是 $I \otimes H^*$ 部分。
  3. TDVP 演化配置：
    • 使用 tenpy.algorithms.tdvp.TwoSiteTDVP。
    • 步长策略：初期使用极小步长 $dt=0.04$ 以允许键维度自然增长，随后可增加至 $dt=0.4$ 或 $4.0$（取决于稳定性）。
  4. 算符测量：
    • 自旋密度：$\langle \langle I | (s_i^z \otimes I) | ho angle angle$。
    • 自旋流：需根据论文 SM2 中的三体算符公式手动编写测量函数。

3.3 模拟参数建议

  • 系统尺寸 $L$: 100 到 2000。
  • 初始极化 $\mu$: 0.05（弱领域墙，确保线性响应）。
  • 边界条件: Open Boundary Conditions (OBC)。

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

4.1 关键文献

  1. Kardar, Parisi, Zhang (1986): 提出了经典的 KPZ 方程，本工作的动力学普适类起源。
  2. Inozemtsev (1990): 定义了连接 HS 和 NN 模型的集成链。
  3. Ljubotina et al. (2017/2019): 首次在数值和实验上确立了 NN Heisenberg 链中的 KPZ 输运。
  4. Bulchandani et al. (2024): 为长程系统的准弹道输运提供了广义流体动力学（GHD）框架。

4.2 局限性评论

尽管该工作展示了令人信服的数值证据，但仍存在以下局限：

  1. 渐近扩散收敛问题：从理论上讲，非集成系统在无限长时间后终将回归扩散。本工作观测到的 KPZ 行为可能是一个“长寿命的中间态交叉”。虽然时间尺度对于实验已足够长，但严格的数学渐近极限仍待探讨。
  2. 三体算符测量的开销：在计算电流关联时，SM2 导出的公式涉及大量三体算符，这使得在大规模 MPS 上的测量变得极其耗时（需要数天），限制了采样频率。
  3. 指数拟合的尾部效应：尽管 4 项指数拟合在短程表现良好，但在极远距离下会失效。对于需要极高红外（IR）精度的物理量，这种近似可能引入偏差。

5. 补充内容：从量子化学到量子模拟的跨界启示

5.1 量子模拟器的“结构不稳定性”

对于量子化学和物理研究者来说，一个重要的启示是：完美的集成性是不稳定的，但集成性的特征是极其稳健的。 在实际的实验装置中，我们永远无法实现完美的 $1/r^2$ 或纯粹的最近邻相互作用。本研究证明了，即便存在这些“扰动”，系统依然会沿着邻近的集成模型路径演化。这意味着，我们可以在不追求完美精度的前提下，利用量子模拟器探索基础的普适性原理。

5.2 输运性质与量子糾纏的关系

在加倍空间的 TDVP 演化中，算符空间的算符纠缠（Operator Entanglement）增长速度决定了模拟的难度。KPZ 行为通常伴随着对数增长的算符纠缠，这正是张量网络方法能够成功的根本原因。如果系统表现为纯粹的热化扩散，纠缠增长会更快，模拟难度将呈指数上升。

5.3 结论与展望

该工作不仅完善了一维自旋输运的拼图，还通过 Inozemtsev 模型族提供了一个统一的视角。未来的研究方向可能包括：

  • 具有 $U(1)$ 破坏项（如 $h_x$ 场）的长程链中的输运转变。
  • 二维长程系统中的类 KPZ 行为探索。
  • 利用本文开发的电流测量 MPO 技术，精确计算各向异性材料的 Drude 权重。

本文为技术深度解析，旨在为同行提供研究背景与落地参考。