光与物质相互作用迪克模型中的复杂度转变

论文信息

标题: Complexity transition in the Dicke model of light-matter interaction

作者: Yicheng Zhang, Erhai Zhao

发布日期: 2026-07-23

arXiv ID: 2607.21583v1

PDF 链接: 下载 PDF

3 分钟速览

研究问题：这篇论文试图回答：在光与物质相互作用的 Dicke 模型中，调节耦合强度 $g$ 是否会在 Krylov 复杂度上引发非解析的 “突变”，从而识别出截然不同的量子动力学机制？
核心方法：作者通过对 Dicke 模型进行量子淬火，数值计算态的 Krylov 复杂度 $C_K$ 和 Krylov 熵 $S_K$ ，并将多体动力学简化为 Krylov 空间中的单粒子波包在 “合成引力” 下的跳问题。
关键结果：在 Dicke 模型中首次发现了 Krylov 复杂度随耦合强度的非单调、急剧变化，并识别出一个明确的 “复杂度相变”。在相变点附近， $\overline{C}_K$ 展现出独特的 $N^{3/2}$ 标度律（图 1c），且 Krylov 熵出现跃升。
主要局限：该研究的数值计算受限于有限尺寸系统（最大 $N=40$ 个自旋）；其对 “混沌” 机制的物理解释依赖于一个线性势的简化玩具模型。
适合读者：对量子多体动力学、量子混沌、算子增长、冷原子/离子阱实验以及 Krylov 方法感兴趣的物理学研究者。

论文背景和研究动机

理解孤立量子多体系统如何从非平衡态演化并最终达到热化，是当前物理学的前沿课题。淬火动力学是其主要研究手段，通过突然改变系统哈密顿量的某个参数，观察系统随后的演化。传统上，研究者运用动力学相变来区分不同的动力学机制，但这通常依赖于选择一个恰当的局域观测量作为 “序参量”（Marino et al., 2022; Heyl, 2018）。此方法的成功极大程度上取决于对系统和初态细节的把握，尤其在混沌体系中，长时间平均的观测量会因剧烈涨落而失效。

近年来，Krylov 复杂度（Krylov complexity）作为一种探究算子增长和量子混沌的新视角脱颖而出。它度量了一个初态或算子在一个由 Lanczos 算法构造的特殊正交基（Krylov 基）下的扩散程度。简单的动力学对应着波包在 Krylov 空间中的局域化（复杂度低），而复杂的、混沌的动力学则对应于波包的延展（复杂度高）。此前有研究提出，Krylov 复杂度的长时间平均值 $\overline{C}_K$ 可以作为指示可积-混沌转变的 “序参量”。然而，一个悬而未决的猜想是：在超越可积极限的真实物理系统中，Krylov 复杂度能否真正有效地识别非解析的动力学相变？

本研究的动机正是为了验证这一猜想。作者选择了 Dicke 模型（Dicke model）作为理想的测试平台，该模型描述了 $N$ 个二能级原子（自旋）与一个光腔模式的耦合，是腔量子电动力学和电路量子电动力学的基础模型。该模型本身不含随机无序，但通过调节光与物质耦合强度 $g$ ，可以平滑地从可积的 Tavis-Cummings 极限过渡到量子混沌区域，且拥有无限大的希尔伯特空间，使其成为探索复杂度相变的绝佳候选。

核心方法和技术细节

论文的核心方法是通过量子淬火动力学结合 Krylov 构造，将复杂的多体问题转化为一个直观的单粒子跳问题。

模型与淬火协议：研究对象是 Dicke 模型哈密顿量（公式 1），其特征是自旋-玻色子耦合强度 $g$ 。论文引入无量纲耦合强度 $\tilde{g} = 2g / \sqrt{\Omega\delta}$ ，其中 $\Omega$ 是自旋能级分裂， $\delta$ 是光子能量。作者模拟了一次量子淬火：从初始的纯直积态——所有自旋向上（ $\theta_0=0$ ）且光腔为真空（ $\alpha_0=0$ ）的经典鞍点——出发，在固定哈密顿量下进行厄米时间演化 $|\psi(t)\rangle = e^{-i\hat{H}t}|\psi_0\rangle$ 。
Krylov 空间的构建：为了量化演化过程的复杂度，作者对初始态及其在哈密顿量反复作用下的集合 $\{\hat{H}^k|\psi_0\rangle\}$ 应用递归 Lanczos 算法进行 Gram-Schmidt 正交化，构建出一组 Krylov 基 $\{|k\rangle\}$ 。在此基下，复杂的多体哈密顿量奇妙地化为一个简单的紧束缚模型（公式 2）： $\hat{H}|k\rangle = a_k|k\rangle + b_k|k-1\rangle + b_{k+1}|k+1\rangle$ 其中， $a_k$ 可视为 $k$ 格点上的在位势能， $b_k$ 则是跃迁振幅。初始态 $|\psi_0\rangle$ 位于格点 $k=0$ 上，整个系统的动力学便被编码为 Krylov 晶格上单粒子波包 $p_k(t) = |\langle k|\psi(t)\rangle|^2$ 的扩散过程。
复杂度与熵的定义：Krylov 复杂度 $C_K(t)$ 被定义为波包的 “质心” 位置： $C_K(t) = \sum_k k p_k(t)$ 。Krylov 熵 $S_K(t)$ 则是概率分布 $p_k$ 的香农熵，度量波包在实空间（Krylov 空间）中的弥散程度。波包扩散得越远、越广， $C_K$ 和 $S_K$ 就越大，意味着动力学越复杂。论文关注的指标是它们的长时间平均值 $\overline{C}_K$ 和 $\overline{S}_K$ ，扫描范围是 $gt \in [50, 200]$ 。

创新点和贡献

本研究最主要的创新在于首次在真实物理模型中，展示了 Krylov 复杂度作为动力学相变探针的有效性与独特优势，并提出了极具洞察力的物理解释。

揭示非解析的复杂度相变：这是论文最核心的贡献。作者绘制了 Dicke 模型详尽的 “复杂度相图”（图 1b），并发现 $\overline{C}_K$ 随着耦合强度 $g$ 的增强呈现非单调变化：先平滑下降，在临界点 $\tilde{g}_c \sim 0.5$ 处发生急剧突变，随后迅速上升并伴有剧烈涨落（图 1a）。Krylov 熵 $\overline{S}_K$ 在相变点也表现出从低值到高值的陡峭跃升，类似一级相变中的潜热（图 1d）。这表明 Krylov 复杂度不仅能区分动力学行为，其本身的突变就构成了一个 “相变”。
提出 “振荡” 对 “面条化” 的物理图像：为了解构这一现象的本质，作者将目光投向 Lanczos 系数 $a_k, b_k$ 。他们发现， $a_k$ 和 $b_k$ 整体上随 $k$ 线性增长，类似于一个线性势阱（ $a_k \sim k$ ）与里德勒时空中的 “合成引力” 跳（ $b_k \sim \eta k$ ）的竞争。常规机制（ $g < g_c$ ）对应较小的 $\eta$ 和强无序散射，波包在 Krylov 空间被有效束缚，其质心 $C_K(t)$ 表现出周期性呼吸振荡。混沌机制（ $g > g_c$ ）中，无序景观变得平缓，线性跳导致的连续红移效应使得波包被不可逆地拉伸并撕碎，无法重聚，犹如天文物体的 “面条化”（spaghettification），导致复杂度飙升。
阐明 Lanczos 系数中 “无序” 的关键作用：论文创造性地将 Lanczos 系数相对于平滑平均的涨落（ $\delta a_k, \delta b_k$ ）定义为 “无序”，并证明其在相变中起决定性作用。在常规区，无序在远处形成一个势垒（图 2e），将波包限制在有限区域内，且随 $g$ 增强，势垒向原点移动，导致 $\overline{C}_K$ 下降。在混沌区，此势垒瓦解，波包得以延伸。一个关键的对比是，采用平滑 Lanczos 系数的模型（去除了无序）完全错过了 $\overline{C}_K$ 的非单调行为和相变（图 5c-d），这无可辩驳地证明了 “无序” 是产生复杂度突变的核心要素。

局限与待解决问题

尽管这项研究为理解量子动力学开辟了新路径，但它也留下了一些值得深入探讨的问题。

有限尺寸效应的限制：所有数值模拟均在有限的自旋数 $N \leq 40$ 下进行，同时玻色子激发数也被截断。论文观察到相变在 $N$ 增大时变得更尖锐（图 1a），但无法提供严格的热力学极限（ $N \to \infty$ ）下的证据。其发现的 $C_K$ 在不同区间的标度律 $N^{3/2}$ 和 $N^{5/2}$ 系数的普适性有待在更大尺度的模拟或实验中检验。
Lanczos 无序的来源与表征：论文将 $a_k$ 和 $b_k$ 的局域涨落称为 “无序”，并直观地展示了其作用，但并未对这种由哈密顿量确定性产生的 “无序” 进行系统性表征。例如，它与通常意义上的随机无序有何区别？其统计性质（如关联函数）如何随 $g$ 变化并决定波包动力学？论文未对此给出严格的定量分析，更多是现象学描述。
与其他混沌诊断指标的对应关系：论文指出复杂度相变 “与混沌的出现大致吻合”（broadly aligned with the onset of chaos），并通过能谱统计和李雅普诺夫指数佐证。然而，精确对比揭示（图 3），这些传统指标并不能给出一个清晰、一致的相变边界，突显了 Krylov 复杂度的优越性。但这也引出一个更深层的问题：Krylov 复杂度所定义的这个无需事先选择 “序参量” 的相变，其物理本质究竟是什么？它到底是对热化、混沌，还是对某种尚未被完全认识的动力学行为的一种全新分类？这为未来的理论工作留下了广阔空间。