光-物质相互作用狄克模型中的复杂性转变

arXiv: 2607.21583v1

论文信息

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

作者: Yicheng Zhang, Erhai Zhao

发布日期: 2026-07-23

arXiv ID: 2607.21583v1

PDF 链接: 下载 PDF

论文背景与研究动机

量子多体系统在远离平衡态的动力学行为一直是现代物理学的前沿热点。传统上，研究者通过 “动力学相变” 这一概念来划分量子淬火后系统的不同演化机制，通常选取某个观测量的长时间平均值作为序参量。然而，这种方法高度依赖于哈密顿量和初始态的具体形式，尤其在系统接近或进入混沌区域时，时间平均观测量往往充满噪声，难以给出清晰的判断。

近年来，一种名为Krylov 复杂度的新工具逐渐进入量子动力学研究的视野。它通过 Lanczos 递归算法构造一组正交基矢（Krylov 基），将复杂的多体时间演化约化为单粒子在 Krylov 晶格上的跳跃问题。初态在 Krylov 空间中的扩散程度衡量了时间演化的复杂度——波包传播得越远、分布得越广，系统动力学就越复杂。

尽管此前已有不少工作尝试将 Krylov 复杂度作为区分可积与混沌动力学的探针，但一个核心问题始终悬而未决：在没有随机无序的物理模型中，改变耦合参数能否导致 Krylov 复杂度或熵出现非解析的突变？ 这正是张一程和赵洱海此篇研究的核心动机。

Dicke 模型与 Krylov 复杂度基础

Dicke 模型的搭建

论文选取的研究对象是描述大量二能级原子与单模腔场耦合的Dicke 模型，其哈密顿量为：

H^=−2gN(a^+a^†)S^x+δa^†a^+ΩS^z\hat{H} = -\frac{2g}{\sqrt{N}}(\hat{a}+\hat{a}^{\dagger})\hat{S}_x + \delta\hat{a}^{\dagger}\hat{a} + \Omega\hat{S}_z

其中 a^†\hat{a}^{\dagger}（a^\hat{a}）是光子的产生（湮灭）算符，S^α\hat{S}_\alpha 为集体自旋算符，gg 为光-物质耦合强度，无量纲化后记为 g~=2g/Ωδ\tilde{g} = 2g/\sqrt{\Omega\delta}。

这个模型之所以值得特别关注，原因有三重：

  1. 实验可及性：冷原子腔 QED 和囚禁离子阵列已经实现了 Dicke 模型的量子模拟，近期甚至观测到了其动力学相变。
  2. 丰富的相图结构：在小耦合极限下系统约化为可积的 Tavis-Cummings 模型，而在大耦合下则呈现出量子混沌特征。
  3. 无限维希尔伯特空间：与有限自旋链不同，Dicke 模型中玻色子模式的存在使得 Krylov 复杂度分析更具普遍意义。

论文聚焦于一个具有明确物理意义的初始态——所有自旋朝上（θ0=0\theta_0=0）且无光子（α0=0\alpha_0=0）的乘积态。在经典极限下，该点对应系统的鞍点，因而是探索动力学不稳定性的理想起点。

Krylov 复杂度的构造流程

从初始态 ∣ψ0⟩|\psi_0\rangle 出发，Lanczos 算法通过逐次施加哈密顿量并 Gram-Schmidt 正交化，生成一组 Krylov 基 {∣k⟩,k=0,1,2,… }\{|k\rangle, k=0,1,2,\dots\}。在这个新基底下，哈密顿量呈现简洁的三对角形式：

H^∣k⟩=ak∣k⟩+bk∣k−1⟩+bk+1∣k+1⟩\hat{H}|k\rangle = a_k|k\rangle + b_k|k-1\rangle + b_{k+1}|k+1\rangle

其中 aka_k 为在位势，bkb_k 为跳跃振幅。原先复杂的多体演化被映射为一维晶格上的单粒子波包动力学，初始波包局域在 k=0k=0 位置。

Krylov 复杂度定义为波包质心的期望位置：

CK(t)=∑kk⋅pk(t)C_K(t) = \sum_k k \cdot p_k(t)

其中 pk(t)=∣⟨k∣ψ(t)⟩∣2p_k(t) = |\langle k|\psi(t)\rangle|^2。相应地，Krylov 熵采用 Shannon 熵形式：

SK(t)=−∑kpk(t)ln⁡pk(t)S_K(t) = -\sum_k p_k(t)\ln p_k(t)

核心发现：复杂性相变的揭示

非单调且陡峭的转变特征

图 1(a)中展示的结果构成了论文最令人震撼的发现：长时间平均 Krylov 复杂度 C‾K\overline{C}_K 随 g~\tilde{g} 并非单调递增，而是先平滑下降，在 g~c∼0.5\tilde{g}_c \sim 0.5 附近突然转折，随后急剧攀升并伴随剧烈波动。

这一趋势在 N=20,30,40N=20,30,40 三种系统尺寸下稳定重现，更关键的是，转变点随 NN 增大而更趋尖锐——这正是相变行为的典型特征。图 1(c)进一步揭示，在 g~<g~c\tilde{g}<\tilde{g}_c 一侧靠近转变点时，C‾K\overline{C}_K 呈现出独特的 N3/2N^{3/2} 标度律，而远在 g~>g~c\tilde{g}>\tilde{g}_c 区域则趋向 N5/2N^{5/2} 标度。标度指数的突变进一步佐证了两个动力学机制之间存在质的差异。

Krylov 熵 S‾K\overline{S}_K（图 1d）给出更富物理意蕴的图景：低 g~\tilde{g} 侧，熵值较低，波包在 Krylov 空间相对局域；越过 g~c\tilde{g}_c 后，熵迅速跳升，特征类似于一级相变中的潜热释放。这不仅是一个定量信号，更直观地揭示了波包构型从 “束缚” 到 “延展” 的突变。

全局相图与两个动力学机制的识别

论文进一步绘出了完整的 C‾K\overline{C}_K 相图（图 1b），覆盖 g~\tilde{g} 和 Ω/δ\Omega/\delta 两个控制参数。蓝色区域对应复杂度极低的状态（自旋锁定于初值、无玻色激发），而红色区域对应混沌的高复杂度区域。Ω=δ\Omega=\delta 共振线的水平切面正对应图 1(a)的曲线。

通过综合传统探针——时间平均磁化强度、平均光子数、自旋-玻色纠缠熵、局域能级间距比和半经典 Lyapunov 指数——论文将两个动力学区域分别定性为 “规则” 和 “混沌” 区。值得注意的是，传统探针在识别精确的转变边界时均存在一定模糊性：能级统计依赖于窗口选择，Lyapunov 指数对截断阈值敏感。相比之下，Krylov 复杂度和熵直接刻画 ∣ψ(t)⟩|\psi(t)\rangle 本身，给出了更为锐利的信号。

物理机制的深度解析

Rindler-Krylov 模型与引力类比

论文最具洞察力的贡献在于通过 Lanczos 系数 (ak,bk)(a_k,b_k) 的渐近行为揭示了转变背后的物理机制。数值结果显示，aka_k 和 bkb_k 在大 kk 下整体呈线性增长（图 2c,d）：

ak∼k,bk∼ηka_k \sim k, \quad b_k \sim \eta k

这种结构具有深刻的时空几何内涵。线性在位势 ak=ka_k = k 代表线性囚禁势，而 bk=ηkb_k = \eta k 的位置依赖跳跃则模拟了Rindler 时空度规——一个匀加速观察者所经历的时空几何。在引力的等效原理视角下，加速观察者感受到指向原点 k=0k=0 的引力，造成向外传播的波包因连续红移而被持续拉伸。

囚禁与解囚禁的竞争：意面化效应

参数 η=bk/ak\eta = b_k/a_k 控制着囚禁与解囚禁的竞争：

  • 当 η<1/2\eta < 1/2 时，囚禁势占据主导，波包维持 “呼吸” 振荡模式，即 CK(t)C_K(t) 呈现周期性脉动。
  • 在临界点 η=1/2\eta = 1/2 处发生解局域化转变，本征态不再局域在原点，而展现幂律衰减，波包持续延展、无法再聚焦。这种极端拉伸现象被形象地类比为黑洞附近的 “意面化” 效应——物体在强引力梯度下被拉成细长条。

Dicke 模型的 η\eta 值随 g~\tilde{g} 单调递增并趋近 0.5，这解释了为何在 g~c\tilde{g}_c 附近系统从振荡模式骤然转入混沌扩散模式。

无序的囚禁作用

仅凭光滑的 Rindler-Krylov 模型尚不足以解释 C‾K\overline{C}_K 在规则区的非单调下降。论文通过引入 “平滑化模型”（以滑动平均消除 ak,bka_k,b_k 的快变涨落）发现，平滑化模型完全丧失了复杂性转变特征——C‾K\overline{C}_K 变得单调递增。这意味着Lanczos 系数的涨落（文中称为 “无序”δak,δbk\delta a_k,\delta b_k）在规则区起着至关重要的额外囚禁作用。

在规则区，δak\delta a_k 和 δbk\delta b_k 呈现单个高大的尖峰（图 2e），该尖峰作为势垒将波包限制在约 k∼450k \sim 450 范围内。随 g~\tilde{g} 增加，尖峰位置向原点移动，囚禁效应增强，C‾K\overline{C}_K 遂平滑下降。达到 g~c\tilde{g}_c 时，尖峰与波包发射区重叠，囚禁达到极限。进入混沌区后，无序结构演变为一系列平缓的宽峰（图 2f），仅对波包产生微弱散射，无法阻止其意面化解体。

数值与解析的交叉验证

论文在 End Matter 中对 Rindler-Krylov 模型的解析推导尤为精妙。在临界点 η=1/2\eta=1/2，本征波函数由合流超几何函数给出：

ϕE(k)=(−1)k1F1(−k;1;2E+1)\phi_E(k) = (-1)^k {}_1F_1(-k;1;2E+1)

其大 kk 渐近行为呈 k−1/4k^{-1/4} 幂律衰减叠加余弦振荡，证实了解局域化的物理本质。

对量子技术与研究的实践启示

尽管本文属于基础理论研究，其成果对量子技术领域具有多重启示意义：

  1. 量子模拟平台的新探针：囚禁离子和超导量子比特实验可直接通过断层扫描技术重构 Krylov 空间中的波包扩散，将 C‾K\overline{C}_K 和 S‾K\overline{S}_K 作为区分动力学机制的 “无模型” 序参量。

  2. 量子控制策略的指引：规则区中波包的周期性呼吸行为意味着系统在特定时间可回到低纠缠态，这对于设计量子存储和相干操控协议具有参考价值。反之，混沌区的意面化效应提示我们，增强耦合并非总导致更高的量子信息扩散，存在复杂度极小点。

  3. 算法复杂度的物理对应：Krylov 复杂度天然与经典计算机模拟量子系统所需的计算资源相关联。C‾K\overline{C}_K 在 g~c\tilde{g}_c 附近达到极小，意味着在中等耦合强度下存在一个 “模拟友好” 窗口，这对经典-量子混合算法的参数选点具有启发意义。

总结与展望

张一程和赵洱海的工作通过绘出 Dicke 模型的 Krylov 复杂度相图，首次在严格意义上回答了 “能否用 Krylov 复杂度捕获量子动力学的非解析转变” 这一开放问题。其关键贡献可凝练为三点：

  • 发现的层面：揭示了复杂度的非单调、陡峭转变和独特的有限尺寸标度律，证明 Krylov 方法能给出比传统探针更锐利的动力学边界。
  • 机制的层面：通过 Rindler-Krylov 模型和意面化引力类比，将两类动力学机制的微观差异归结为波包在合成引力场中的囚禁-解囚禁转变。
  • 方法的层面：通过平滑化模型与原始模型的对比，厘清了 Lanczos 系数涨落在波包囚禁中的决定性作用，为分析其他量子多体系统提供了可迁移的框架。

展望未来，这一框架可推广至更多实验可及的量子多体模型（如 Bose-Hubbard 模型、SYQ 模型等），也可反过来指导量子模拟器的参数设计。尤其值得期待的是，能否在囚禁离子二维晶格中直接观测到 Krylov 空间的呼吸振荡与意面化转变，从而在实验室中实现 “合成引力” 下波包动力学的量子模拟。