来源论文: https://arxiv.org/abs/2606.09106v1 生成时间: Jun 09, 2026 11:33

深度解析 1D Hubbard 模型中的 Family–Vicsek 普适性与动态标度律：无限高温下的量子输运新视角

0. 执行摘要

在凝聚态物理与量子多体动力学领域，理解非平衡态下的输运普适性是一项核心挑战。传统的宏观流体力学描述往往忽略了量子涨落的微观细节，而 Family–Vicsek (FV) 标度律——最初源自经典界面生长理论——为连接微观量子动力学与宏观流体力学提供了一个强有力的标度框架。本文解析的工作针对一维（1D）费米-哈伯德（Hubbard）模型，在无限高温（Infinite Temperature）这一极端且具有代表性的动力学背景下，通过量子生成函数（QGF）方法，系统地研究了电荷、自旋及能量密度的亚系统波动（Roughness）。

核心发现包括：

  1. 普适标度验证：在量子多体系统中，守恒量的转移涨落完全遵循 FV 标度律 $W(\ell, t) \sim \ell^\alpha f(t/\ell^z)$。
  2. 可积性的决定性作用：系统的动力学指数 $z$ 由其可积性决定。在自由极限下，所有扇区呈现弹道输运（Ballistic, $z=1$）；在可积相互作用极限下，电荷与自旋扇区表现为 KPZ（Kardar–Parisi–Zhang）超扩散（$z=3/2$），而能量扇区保持弹道性；一旦引入次近邻相互作用破坏可积性，所有扇区均退化为常规扩散（Diffusive, $z=2$）。
  3. 微观与宏观的交叠：识别出了一个普适的短时间微观机制（$t^2$ 增长），它独立于相互作用强度和可积性，标志着流体力学演化的前奏。

该研究不仅加深了对 Hubbard 模型高能激发态的理解，也为利用冷原子量子模拟器观测非平凡标度行为提供了理论路线图。

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

1.1 核心科学问题

一维量子系统的输运性质往往表现出与高维系统截然不同的奇异性，其根本原因在于一维空间对准粒子散射的严格约束。Hubbard 模型作为关联电子系统的范式模型，在 1D 下具有 Bethe Ansatz 可积性。然而，在无限高温极限下（即所有本征态具有相同的权重），这些守恒律如何从微观层面涌现为宏观的标度行为？具体而言，电荷、自旋和能量这三种不同的守恒量是否遵循相同的普适类？

1.2 理论基础：Family–Vicsek 标度律

FV 标度律描述了长度为 $\ell$ 的亚系统在时间 $t$ 时的粗糙度（Roughness）$W(\ell, t)$：

$$ W(\ell, t) \sim \ell^\alpha f\left(\frac{t}{\ell^z}\right) $$

其中：

  • $\alpha$ (粗糙度指数)：描述饱和时波动随尺寸的增长。在 1D 量子系统中，通常 $\alpha = 1/2$。
  • $z$ (动力学指数)：定义了时间与空间的标度比例，$t \sim \ell^z$。
  • $\beta$ (增长指数)：满足关系 $z = \alpha/\beta$。对于小时间极限，$W(\ell, t) \sim t^\beta$。

本工作的创新之处在于将“界面高度”替换为亚系统中守恒量（如电荷算符 $N_\ell$）的二阶累积量（Cumulant）。

1.3 技术难点：量子生成函数 (QGF) 的构建

在无限高温下，直接计算长时间的相关函数极具挑战，尤其是对于具有非阿贝尔对称性（如 $SU(2)$ 自旋与 $SU_c(2)$ 赝自旋）的系统。传统的 MPS 方法在模拟长时间演化时会面临算符纠缠（Operator Entanglement）快速增长的问题。

本研究采用的 QGF 方案： 定义扭转算符（Twist Operator）：

$$ R_\ell^{(Q)}(\lambda) = e^{i \lambda Q_\ell} $$

通过计算幺正量子生成函数：

$$ G_{\ell, Q}^{(H)}(\lambda, t) = \text{Tr}\left[ R_\ell^{(Q)}(\lambda) U(t) R_\ell^{(Q)\dagger}(\lambda) \rho_\infty U^\dagger(t) \right] $$

其中 $\rho_\infty = 4^{-L} \mathbb{1}$。通过对 $\lambda$ 进行级数展开，可以直接提取出二阶累积量（即粗糙度的平方），而无需重构整个概率分布。为了减少数值误差，作者采用了虚数相位和实数相位的组合提取技术（Eq. 16）。

1.4 方法细节：能量扇区的特殊处理

电荷与自旋的扭转算符在 MPS 形式下是简单的局域算符乘积。但能量密度算符 $h_j$ 涉及相邻格点的相互作用：

$$ h_j = -J \sum_\sigma (c_{j\sigma}^\dagger c_{j+1,\sigma} + h.c.) + \frac{U}{2}(n_{j\uparrow}-1/2)(n_{j\downarrow}-1/2) + \dots $$

由于 $[h_j, h_{j+1}] \neq 0$，能量扭转算符 $R_\ell^{(E)}(\lambda)$ 无法直接分解。作者引入了 二阶 Suzuki-Trotter 分解，将能量算符分为奇偶键贡献（$E_{odd}$ 和 $E_{even}$），从而在 MPO（Matrix Product Operator）框架下高效实现能量涨落的计算。这是处理能量输运涨落的一个重要技术改进。

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

2.1 自由费米子极限 ($U=0, V=0$)

作为数值方法的基准，作者首先解析求解了自由费米子情况。在该极限下：

  • 解析结果：通过 Wick 定理，二阶矩 $\kappa_2^{(N)}(t) \sim \sum |r| J_r^2(2Jt)$。在长时极限下，其表现为线性的 $t$ 增长，意味着 $\beta = 1/2, z = 1$。
  • 数据表现：图 3 展示了完美的标度坍缩（Scaling Collapse）。当使用 $\ell^{1/2}$ 缩放纵轴，使用 $Jt/\ell$ 缩放横轴时，不同长度 $\ell$ 的曲线完美重合。饱和值 WN,sat 正比于 $\ell^{1/2}$，完全符合预测。

2.2 相互作用可积 Hubbard 模型 ($U=2, V=0$)

这是本工作最引人注目的部分。在该体系中，电荷与自旋呈现出 KPZ 普适类：

  • 动力学指数 $z = 3/2$：这标志着超扩散（Superdiffusion）行为。
  • 增长指数 $\beta = 1/3$：计算得到的二阶累积量遵循 $\kappa_2 \sim t^{2/3}$，因此 $W \sim t^{1/3}$。
  • 能量扇区的例外：如图 6(c) 所示，即便在 $U \neq 0$ 时，能量波动依然遵循 $z=1$ 的弹道标度。这是因为能量电流不直接耦合到导致 KPZ 异常的非阿贝尔荷上。

2.3 可积性破坏体系 ($U=2, V=2$)

引入次近邻密度相互作用 $V$ 后，系统不再可积。此时：

  • 全扇区扩散：电荷、自旋和能量均坍缩至 $z=2$ 的曲线。
  • 增长指数 $\beta = 1/4$：二阶累积量 $\kappa_2 \sim t^{1/2}$，对应常规扩散过程。
  • 数据对比：图 5 和 图 6(b, d) 清晰地展示了从 $z=1$ 或 $z=1.5$ 向 $z=2$ 的转变。这证明了 FV 标度律可以作为探测多体系统可积性破坏的灵敏“探针”。

2.4 关键数据汇总表 (Table I 核心内容)

扇区制度 (Regime)参数条件$\alpha$$\beta$$z$
电荷/自旋弹道 (Ballistic)$U=0, V=0$1/21/21
电荷/自旋KPZ$U \neq 0, V=0$1/21/33/2
电荷/自旋扩散 (Diffusive)$V \neq 0$1/21/42
能量弹道 (Ballistic)$V=0, \forall U$1/21/21
能量扩散 (Diffusive)$V \neq 0, \forall U$1/21/42

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

3.1 核心算法：Liouville 空间中的 TEBD

由于研究对象是无限高温下的密度矩阵 $\rho_\infty$，计算并非在 Hilbert 空间进行，而是在 Liouville 算符空间（Operator Space）进行。扭转算符作为初态，在时间演化算符 $L = [H, \cdot]$ 的作用下演化。

复现步骤建议：

  1. 算符表示：将费米子算符通过 Jordan-Wigner 变换（或直接使用费米子张量网络）映射到自旋链形式。
  2. MPO 初始化：根据 Eq. 57 构建作用于指定段 $\ell$ 的扭转算符 $R_\ell^{(Q)}(\lambda)$。对于能量算符，需实现 Suzuki-Trotter 层的 MPO 乘法。
  3. TEBD 演化：
    • 使用二阶或四阶 Trotter 步进。
    • 关键点：保持算符的非阿贝尔对称性。推荐使用支持 $SU(2)$ 对称性保护的库（如 ITensor 或自定义块角矩阵库）。
    • 截断误差控制：由于是无限高温，算符纠缠增长较快，建议 Bond Dimension $\chi$ 至少达到 500-1000 以上。
  4. 累积量提取：执行两次演化（分别对应 $\lambda = r$ 和 $\lambda = ir$），利用 Eq. 16 进行差分计算以抵消低阶截断误差。

3.2 开源工具推荐

虽然论文作者使用的是其开发的专用代码，但以下开源项目具备复现该工作的能力：

  • ITensor (C++/Julia): 强大的张量网络库，非常适合构建自定义 MPO 和处理非阿贝尔对称性。
  • TeNPy (Python): 提供了成熟的 TEBD 和有限尺寸算法接口。
  • QuSpin (Python): 用于小尺寸精确对角化（ED）验证短时间行为。

GitHub 资源链接参考：

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

4.1 关键引用文献

  1. Family & Vicsek (1985): 标度律的奠基之作，定义了非平衡生长的通用框架。[Ref 2]
  2. Lieb & Wu (1968): Hubbard 模型可积性的解析解基础。[Ref 20]
  3. Prosen & Žnidarič (2012): 首次提出 1D 系统中高温扩散输运的数值证据。[Ref 29]
  4. Bulchandani et al. (2021): 综述了可积链中的超扩散与广义流体力学（GHD）。[Ref 33]
  5. Moca et al. (2025/2026): 本团队前期关于自旋链 FV 标度的研究。[Ref 16, 17]

4.2 局限性评论

  1. 无限高温假设：虽然无限高温简化了统计权重，但在真实的冷原子实验中，系统总是处于有限温度。有限温度下对称性破缺对 KPZ 标度的影响尚未在本项目中完全探讨。
  2. 一维局限性：FV 标度的普适性在 1D 表现完美，但在 2D 或准 1D 系统中，由于 $z$ 和 $\alpha$ 的变化，量子涨落的拓扑约束会减弱，结论可能不再成立。
  3. 纠缠限制：TEBD 演化时间受算符纠缠熵（Operator Entanglement Entropy）线性增长的影响，导致无法触及极长的时间窗口。虽然本工作通过标度坍缩部分弥补了这一点，但对于极小频率下的输运系数提取仍存在挑战。

5. 补充：量子涨落与经典界面的深层联系

5.1 为什么量子波动看起来像“界面”？

这是一个物理直觉上的精彩映射。在经典物理中，界面的粗糙度来源于随机沉积导致的局部高度波动。在量子系统中，守恒量（如电荷）在亚系统边界的“流入”与“流出”类似于粒子的沉积与蒸发。由于电荷是守恒的，其波动的累积实际上记录了通过边界的电流的时空积分。这种时间累积效应正是导致 $W \sim t^\beta$ 的动力学根源。

5.2 对称性与赝自旋的物理意义

Hubbard 模型在半填充（Half-filling）处拥有一种特殊的 $SU_c(2)$ 赝自旋对称性（又称 $\eta$-pairing 对称性）。本工作指出，正是这种非阿贝尔对称性的存在，导致了电荷扇区在可积相互作用情形下出现了与自旋扇区相同的 KPZ 行为。这揭示了输运普适类并不取决于算符的具体物理意义（是自旋还是电荷），而取决于其背后的对称性群结构。

5.3 实验观测的可能性

利用现有的量子气体显微镜（Quantum Gas Microscopy）技术，实验学家已经可以实时观测 1D 光晶格中原子的密度涨落。本研究建议，通过测量不同长度采样窗口内的粒子数标准差随时间的演化，即可直接在实验上绘制出 FV 标度曲线。这比测量复杂的高阶关联函数更具操作性，是连接理论与实验的理想桥梁。

5.4 总结与展望

这项工作证明了 Family–Vicsek 标度律是理解关联量子物质动力学的一把钥匙。通过将复杂的量子多体演化映射为简单的标度函数，研究者能够清晰地分辨出可积性、相互作用与输运性质之间的复杂关联。未来的研究方向可能包括：多组分系统的多标度行为（Multiscaling）、无序导致的定位（MBL）对 FV 标度的破坏，以及开放量子系统中的耗散驱动标度行为。