来源论文: https://arxiv.org/abs/2603.16709v1 生成时间: Mar 18, 2026 03:20

开放量子拉比模型中的 Kibble-Zurek 机制：非马尔可夫记忆效应下的普适标度深度解析

0. 执行摘要

在非平衡统计物理与量子模拟的前沿领域，理解系统穿越相变点时的缺陷产生规律是核心课题之一。Kibble-Zurek 机制（KZM）为预测这种非平衡动力学中的普适幂律标度提供了基本框架。然而，当系统与环境耦合形成开放量子系统时，传统的 KZM 往往因马尔可夫耗散（无记忆环境）而失效。本文深度解析了一项突破性研究：在开放量子拉比模型（Open Quantum Rabi Model, OQRM）中，引入非马尔可夫欧姆库（Ohmic bath）不仅诱导了Berezinskii-Kosterlitz-Thouless (BKT) 相变，更出人意料地保持了 KZM 的普适标度完整性。

本研究利用先进的矩阵乘积态（MPS）数值技术，结合密度矩阵重整化群（DMRG）与时变变分原理（TDVP），证明了在非马尔可夫环境下，耗散并不与绝热动力学竞争，而是重新定义了相变的普适类。这一发现确立了 KZM 作为开放系统中相变关键性的可靠见证者，为量子态工程与绝热量子计算提供了理论支撑。

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

1.1 核心科学问题：开放系统中的 KZM 悖论

Kibble-Zurek 机制的核心思想是“冻结（Freeze-out）”假设。当一个系统缓慢扫过临界点时，由于临界减速（Critical Slowing Down），其弛豫时间 $\tau$ 会趋于发散。当扫掠速率超过系统响应能力时，系统状态会停止随瞬时基态演化而“冻结”，从而产生激发态（缺陷）。

在孤立系统中，激发能 $E_{exc}$ 与扫掠时间 $t_Q$ 遵循普适幂律：$E_{exc} \propto t_Q^{-\mu}$。但在开放系统中，马尔可夫环境引入的随机噪声通常会破坏这一标度。本研究的问题在于：如果环境具有记忆效应（非马尔可夫性），KZM 的普适性是否能够幸存？

1.2 理论基础：量子拉比模型与欧姆库

量子拉比模型（QRM）描述了一个两能级系统（Qubit）与一个谐振子模式的耦合：

$$H_{Q-O} = -\frac{\Delta}{2}\sigma_x + \omega_0 a^\dagger a + g\sigma_z(a^\dagger + a)$$

其中 $\Delta$ 是隧道能，$g$ 是耦合强度。当该模型与欧姆库耦合时，总哈密顿量变为 $H = H_{Q-O} + H_I + H_B$。库的特征由欧姆谱密度定义：

$$J(\omega) = \alpha\omega\Theta(\omega_c - \omega)$$

通过积分掉库自由度，系统可以映射到一个具有长程相互作用的一维经典自旋链，该链在特定耦合强度下经历 BKT 相变。BKT 相变的特征是弛豫时间具有指数型发散，而非普通的幂律发散：

$$\tau \sim \exp\left(\frac{B}{\sqrt{|g - g_c|}}\right)$$

1.3 技术难点：多体动力学与连续谱模拟

模拟开放量子拉比模型的难点在于：

  1. 波色激发发散：在临界点附近，谐振子模式和环境的波色子激发数量极大，传统的基底截断法失效。
  2. 长程相互作用：非马尔可夫环境在时间维度上引入了复杂的自相关，在 MPS 表示中对应于空间上的长程耦合。
  3. 计算成本：为了捕获 BKT 相变的微妙细节，需要极高的数值精度和足够大的环境自由度截断。

1.4 方法细节：MPS, DMRG 与 TDVP

作者采用了基于张量网络（Tensor Networks）的方案：

  • MPS 几何结构：使用“星型结构（Star Geometry）”离散化欧姆库。这种方法允许将环境模式直接耦合到中心系统，避免了链式结构中能量尺度的过快衰减。
  • DMRG（基态寻优）：用于确定不同耦合强度 $g$ 下的瞬时基态 $|\psi_{GS}\rangle$，特别是为了寻找临界点 $g_c$。
  • TDVP（时间演化）：这是模拟非平衡淬火过程的关键。TDVP 允许在 MPS 流形上投影薛定谔方程，能够高效处理具有大量环境模式的非马尔可夫动力学。相比于传统的 TEBD，TDVP 在处理长程力和非单位演化算符时更具优势。

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

2.1 模拟参数设置

研究选定的参数体系如下：

  • 耗散强度 $\alpha = 0.2$（弱耗散区）。
  • 截止频率 $\omega_c = 10\Delta$。
  • 谐振子频率 $\omega_0 = 0.75\Delta$。
  • 临界耦合强度通过蒙特卡洛与有限尺寸标度分析确定为 $g_c/\Delta \approx 0.9165$。

2.2 弛豫动力学数据 (Benchmark 1)

为了验证 BKT 性质，研究分析了纵向自旋磁化强度的弛豫函数 $\Sigma_z(t) = \langle\sigma_z(t)\rangle/\langle\sigma_z(0)\rangle$。

  • 结果：在 $g < g_c$ 时，系统表现出相干拉比振荡；接近 $g_c$ 时，进入过阻尼状态。
  • 拟合结果：弛豫时间 $\tau$ 完美符合 BKT 标度公式 $\tau \propto \exp(B/\sqrt{g_c - g})$，拟合得到的 $g_c/\Delta \approx 0.918$，与平衡态预测误差仅为 $10^{-3}$。这证明了动态过程与平衡态临界性的一致性。

2.3 淬火动力学与 KZM 标度 (Benchmark 2)

执行线性淬火协议 $g(t) = (g_f/t_Q)t$。

  • 冻结时间 $t_f$：通过超越方程 $Ae^{\sqrt{B/|g-g_c|}} = \frac{t_Q}{g_f}|g-g_c|$ 计算。计算利用了 Lambert W 函数。
  • 激发能 $E_{exc}$：在 $t_f$ 时刻评估系统能量与瞬时基态能量之差。
  • 标度指数：数值模拟显示 $E_{exc} \propto t_f^{-0.992}$。
  • 激发概率 $P_{exc}$：表现出 $P_{exc} \propto t_f^{-1.07}$ 的幂律。

2.4 性能数据说明

MPS 模拟中，环境模式被离散化为多达 100-200 个谐振子，单个演化路径在高性能计算集群上耗时数十至数百 CPU 小时。TDVP 的步长保持在 $0.01/\Delta - 0.05/\Delta$ 以保证收敛性。

3. 代码实现细节，复现指南与开源链接

3.1 软件包推荐：ITensor

论文明确使用了 ITensor 库（Ref [42]），这是一个基于 Julia 和 C++ 的张量网络开发利器，特别适合处理 MPS 与多体哈密顿量演化。

3.2 复现指南

  1. 离散化库：使用欧姆谱密度 $J(\omega)$。推荐使用线性或对数离散化方案。在 ITensor 中，定义一个 SiteType 为 "Oscillator" 的数组来表示库。
  2. 构建 MPO：量子拉比模型的项是局域的，但由于星型结构，你需要构建一个特殊的 MPO (Matrix Product Operator)。对于星型结构，建议直接使用 AutoMPO 功能，手动指定所有 $\sigma_z$ 与环境模式 $x_i$ 之间的相互作用。
  3. 基态寻找：运行 DMRG 算法获取不同 $g$ 下的基态。确保键维（Bond Dimension）足够大（例如 $D \geq 500$），并监控截断误差（Truncation Error）。
  4. 淬火演化：
    • 初始化 MPS 为 $t=0$ 的基态。
    • 使用 tdvp 函数进行演化。注意 tdvp 在 ITensor 中通常作为插件存在（如 ITensorTDVP.jl）。
    • 在每个步长更新 Hamiltonian 中的 $g(t)$。
  5. 后处理：计算 $\langle \psi(t) | H(g(t)) | \psi(t) \rangle - E_{GS}(g(t))$。

3.3 开源实现参考

虽然作者未直接放出本论文的私有仓库，但可以参考以下类似的非马尔可夫演化项目：

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

4.1 关键引用文献

  1. Kibble (1976) & Zurek (1985)：KZM 的奠基性工作（Ref [7, 8]）。
  2. Kosterlitz & Thouless (1973)：BKT 相变的理论源头（Ref [23]）。
  3. Hwang, Puebla, & Plenio (2015)：量子拉比模型中的量子相变研究（Ref [14]），本研究将其扩展到了开放系统。
  4. Polkovnikov (2005)：绝热微扰论视角下的 KZM 标度推导（Ref [16]）。
  5. Damski (2005)： Landau-Zener 映射到 KZM 的经典论文（Ref [32]）。

4.2 局限性评论

尽管本工作极具创新性，但仍存在以下局限：

  • 空间维度的缺失：研究集中在 $d=0$ 的系统。虽然文章论证了零维系统也能通过有效映射展现 KZM，但在真实的高维凝聚态系统中，空间关联长度 $\xi$ 的竞争可能会引入更复杂的效应。
  • 欧姆谱密度的特异性：结论高度依赖于非马尔可夫记忆（欧姆库）。对于亚欧姆（Sub-Ohmic）或超欧姆（Super-Ohmic）库，BKT 相变可能会演变为其他类型的相变，KZM 是否依然鲁棒尚待验证。
  • 计算截断效应：波色子能级的截断始终是 MPS 模拟中的“幽灵”。虽然作者宣称收敛，但在极强耦合区，光子数的非线性增长可能导致细微的系统偏差。

5. 补充说明：为什么这对于量子化学和计算物理很重要？

5.1 从原子尺度到量子模拟器

对于量子化学研究者而言，量子拉比模型不仅仅是一个抽象的物理模型，它是腔量子电动力学 (QED) 和 光物质相互作用 的基石。在极强耦合（Deep Strong Coupling）体制下，分子振动与腔模的强耦合会导致极化激元（Polariton）的形成，从而改变分子的化学反应路径。

5.2 绝热制备策略的优化

本论文给出的 KZM 标度规律直接指导了我们如何通过“时间换精度”。如果你在实验室中试图制备一个高度纠缠的量子基态，非马尔可夫记忆不仅不是敌人，反而是某种程度上的“守护者”——它确保了即使在耗散存在的情况下，如果你遵循正确的 KZM 扫掠速率，你依然能获得普适的、可预测的结果。

5.3 动力学临界性作为诊断工具

传统上，我们通过测量静态关联函数来判定相变。本文展示了一种全新的思路：通过观察系统对外界干扰的弛豫速度（弛豫动力学），我们可以精确锁定相变点。这种动态探测技术在嘈杂的实验环境（如超导量子比特或捕获离子平台）中具有极高的实用价值。

5.4 总结

这项工作弥合了非平衡动力学标度理论与真实开放量子系统之间的鸿沟。它告诉我们，记忆效应是维持物理学对称性和普适性的纽带。在未来的量子模拟实验中，利用非马尔可夫环境工程，我们或许能设计出更具鲁棒性的量子相变路径。

关键词：#量子拉比模型 #Kibble-Zurek #非马尔可夫 #张量网络 #BKT相变 #量子模拟