2025-11-12T03:25:30.825484

Dissecting Superradiant Phase Transition in the Quantum Rabi Model

Yang, Luo
The phase transition is both thermodynamically and quantum-mechanically ubiquitous in nature or laboratory and its understanding is still one of most active issues in modern physics and related disciplines. The Landau's theory provides a general framework to describe \textit{phenomenologically} the phase transition by the introduction of order parameters and the associated symmetry breakings; and is also taken as starting point to explore the critical phenomena in connection with phase transitions in renormalization group, which provides a complete theoretical description of the behavior close to the critical points. In this sense the microscopic mechanism of the phase transition remains still to be uncovered. Here we make a first attempt to explore the microscopic mechanism of the superradiant phase transition in the quantum Rabi model (QRM). We firstly perform a diagonalization in an operator space to obtain three fundamental patterns involved in the QRM and then analyze explicitly their energy evolutions with increasing coupling strengths. The characteristic behaviors found uncover the microscipic mechanism of the superradiant phase transition: one is active to drive the happening of phase transition, the second responses rapidly to the change of the active pattern and wakes up the third pattern to stablize the new phase. This kind of dissecting mechanism explains for the first time why and how happens the superradiant phase transition in the QRM and paves a way to explore the microscopic mechanism of the phase transitions happening popularly in nature.
academic

Dissecting Superradiant Phase Transition in the Quantum Rabi Model

基本信息

  • 论文ID: 2212.05186
  • 标题: Dissecting Superradiant Phase Transition in the Quantum Rabi Model
  • 作者: Yun-Tong Yang, Hong-Gang Luo (兰州大学物理科学与技术学院，兰州理论物理中心)
  • 分类: quant-ph (量子物理)
  • 提交时间: 2022年12月10日
  • 论文链接: https://arxiv.org/abs/2212.05186

摘要

相变现象在自然界和实验室中普遍存在，其理解是现代物理学最活跃的研究课题之一。朗道理论通过引入序参量和对称性破缺提供了描述相变的唯象框架，重整化群理论则完整描述了临界点附近的行为。然而，相变的微观机制仍未被充分揭示。本文首次尝试探索量子Rabi模型(QRM)中超辐射相变的微观机制。通过在算符空间进行对角化，获得QRM涉及的三个基本模式(patterns)，并明确分析它们随耦合强度增加的能量演化。研究发现：第一个模式主动驱动相变发生，第二个模式快速响应并唤醒第三个模式来稳定新相。这种解剖机制首次解释了QRM中超辐射相变为何及如何发生，为探索自然界普遍存在的相变微观机制开辟了道路。

研究背景与动机

要解决的核心问题

本文旨在回答量子相变领域的根本性问题：相变为何发生(why)以及如何发生(how)。尽管朗道理论和重整化群方法在唯象和临界行为描述上取得巨大成功，但相变的微观动力学机制仍是未解之谜。

问题的重要性

  1. 理论意义：相变理解是现代物理学核心问题之一，从1822年de la Tour观察超临界液体开始，历经200年发展，但微观机制仍不清楚
  2. 普适性：量子Rabi模型是描述光与物质相互作用的基本模型，理解其相变机制对广泛物理系统具有指导意义
  3. 实验进展：超辐射相变已在单个囚禁离子38,39和核磁共振平台40中实现，急需理论解释

现有方法的局限性

  1. 朗道理论：仅提供唯象描述，基于平均场处理，无法揭示微观机制
  2. 重整化群：完整描述临界点附近行为，但本质上是"描述性"而非"解释性"理论
  3. 传统理解：将相变归因于热力学或量子涨落及不同序之间的竞争，但未回答"为何"和"如何"的问题

研究动机

作者在先前工作37中指出，超辐射相变源于耦合强度增加导致的量子双阱势，本文进一步深入探索相变发生的动力学过程和微观机制。

核心贡献

  1. 方法论创新：提出算符空间对角化方法，将量子Rabi模型哈密顿量精确分解为三个基本模式(patterns)，每个模式由特征能量和波函数分量标记
  2. 微观机制揭示：首次明确了超辐射相变的三模式动力学机制：
    • 模式λ₁(主动驱动者)：主动降低系统能量，启动相变
    • 模式λ₂(敏感探测者)：快速响应λ₁的变化，携带超辐射光子，唤醒λ₃
    • 模式λ₃(被动稳定者)：被唤醒后平衡λ₁，稳定新相
  3. 定量分析：通过能量及其一阶、二阶导数的详细分析，定量描述了三个模式在相变过程中的角色转换和相互作用
  4. 理论验证：模式分解方法与数值精确对角化结果完全一致，证明方法的严格性和有效性
  5. 普适性框架：为探索其他复杂系统(如Dicke模型、自旋-玻色子模型)的相变微观机制提供了新范式

方法详解

任务定义

研究目标：解析量子Rabi模型中超辐射量子相变的微观动力学机制

输入：标准量子Rabi模型哈密顿量 H^=ℏω(a^†a^+Δ2σ^x+g(a^+a^†)σ^z)\hat{H} = \hbar\omega\left(\hat{a}^\dagger\hat{a} + \frac{\Delta}{2}\hat{\sigma}_x + g(\hat{a}+\hat{a}^\dagger)\hat{\sigma}_z\right)

其中：

  • a^†(a^)\hat{a}^\dagger(\hat{a})：单光子模式的产生(湮灭)算符
  • σ^x,σ^z\hat{\sigma}_x, \hat{\sigma}_z：泡利矩阵，描述二能级原子
  • Δ\Delta：能级间隔(无量纲)
  • gg：耦合强度(无量纲)
  • 临界点：gc=1+1+Δ2/16g_c = \sqrt{1+\sqrt{1+\Delta^2/16}}

输出：三个基本模式的特征及其在相变过程中的动力学演化

模型架构：算符空间对角化

1. 哈密顿量重构

利用关系σ^yσ^z=iσ^x\hat{\sigma}_y\hat{\sigma}_z = i\hat{\sigma}_x，将哈密顿量改写为二次型：

H^=(−iσ^yσ^za^†)(0Δ40Δ40g0g1)(iσ^yσ^za^)\hat{H} = \begin{pmatrix} -i\hat{\sigma}_y & \hat{\sigma}_z & \hat{a}^\dagger \end{pmatrix} \begin{pmatrix} 0 & \frac{\Delta}{4} & 0 \\ \frac{\Delta}{4} & 0 & g \\ 0 & g & 1 \end{pmatrix} \begin{pmatrix} i\hat{\sigma}_y \\ \hat{\sigma}_z \\ \hat{a} \end{pmatrix}

2. 模式分解

对3×3矩阵进行本征分解，得到三个本征值{λn}\{\lambda_n\}和本征矢{un}\{u_n\} (n=1,2,3)，哈密顿量可表示为：

H^=∑n=13λnA^n†A^n\hat{H} = \sum_{n=1}^3 \lambda_n \hat{A}_n^\dagger \hat{A}_n

其中模式算符： A^n=un,1(iσ^y)+un,2σ^z+un,3a^\hat{A}_n = u_{n,1}(i\hat{\sigma}_y) + u_{n,2}\hat{\sigma}_z + u_{n,3}\hat{a}

3. 完备基展开

在完备基∣σz,m⟩|\sigma_z, m\rangle中求解(其中σ^z∣σz⟩=±∣σz⟩\hat{\sigma}_z|\sigma_z\rangle = \pm|\sigma_z\rangle，a^†a^∣m⟩=m∣m⟩\hat{a}^\dagger\hat{a}|m\rangle = m|m\rangle，m=0,1,⋯ ,Nm=0,1,\cdots,N为截断Fock基)：

⟨σz,m∣H^∣m′,σz′⟩=∑n=13λn∑σz′′,q[A^n†]m,σz;q,σz′′[A^n]q,σz′′;m′,σz′\langle\sigma_z,m|\hat{H}|m',\sigma_z'\rangle = \sum_{n=1}^3 \lambda_n \sum_{\sigma_z'',q} [\hat{A}_n^\dagger]_{m,\sigma_z;q,\sigma_z''}[\hat{A}_n]_{q,\sigma_z'';m',\sigma_z'}

采用N=200N=200的截断足以描述所关注的物理。

技术创新点

1. 算符空间对角化 vs 传统希尔伯特空间对角化

  • 传统方法：直接在Fock基上数值对角化哈密顿量矩阵，得到能级和本征态，但无法揭示内部结构
  • 本文方法：在算符空间先行对角化，将哈密顿量分解为三个"基本激发模式"，每个模式保持特定的算符组合关系

2. 模式的物理意义

三个模式具有明确的物理特征(图1)：

  • 模式λ₁：负本征能量，iσ^yi\hat{\sigma}_y和a^\hat{a}同相，与σ^z\hat{\sigma}_z反相
  • 模式λ₂：接近零的本征能量，主要由光子模式a^\hat{a}主导
  • 模式λ₃：正本征能量，iσ^yi\hat{\sigma}_y和a^\hat{a}反相，与σ^z\hat{\sigma}_z同相

3. 相位关系的不变性

关键发现：各模式内部的相对相位在整个耦合强度范围保持不变，这是标记模式的充要条件，反映了模式的本征特性

4. 精确性保证

该分解在给定Fock基下是精确的——哈密顿量矩阵可完全重构，与数值精确对角化结果完全一致(图1c-e)

实验设置

参数选择

  • 能级间隔：Δ=50\Delta = 50 (固定，对应大失谐情况)
  • 耦合强度范围：g/gc∈[0,2]g/g_c \in [0, 2]，覆盖正常相(g<gcg<g_c)、临界点(g∼gcg\sim g_c)和超辐射相(g>gcg>g_c)
  • 临界耦合强度：gc=1+1+Δ2/16≈2.55g_c = \sqrt{1+\sqrt{1+\Delta^2/16}} \approx 2.55 (对于Δ=50\Delta=50)
  • Fock基截断：N=200N=200

分析方法

  1. 能量分析：计算基态和第一激发态的总能量及三个模式的能量贡献
  2. 导数分析：计算能量关于耦合强度的一阶和二阶导数，揭示动力学过程
  3. 波函数分析：在Fock基中展示波函数的模式分解
  4. 物理量计算：光子数⟨a^†a^⟩\langle\hat{a}^\dagger\hat{a}\rangle和自旋翻转⟨σ^x⟩\langle\hat{\sigma}_x\rangle

验证策略

通过以下三个方面验证模式分解的有效性：

  1. 能谱对比：模式能量求和与精确对角化结果完全一致(图1c)
  2. 光子数对比：各模式光子数求和与精确结果一致(图1d)
  3. 自旋翻转对比：各模式⟨σ^x⟩\langle\hat{\sigma}_x\rangle求和与精确结果一致(图1e)

实验结果

主要结果：三模式动力学机制

1. 基态能量演化(图2a)

正常相 (g/gc<1g/g_c < 1)：

  • 三个模式能量几乎不变
  • 总能量保持稳定

临界区 (g/gc∼1g/g_c \sim 1)：

  • λ₁首先下降：主动降低能量，启动相变
  • λ₂快速补偿：以相反方式响应，暂时稳定总能量
  • λ₃被唤醒：开始响应但较缓慢

超辐射相 (g/gc>1g/g_c > 1)：

  • λ₃权重超越λ₂：开始主导平衡λ₁
  • λ₁和λ₃达到新平衡：两者大致抵消
  • 总能量平滑下降：进入稳定的超辐射相

2. 二阶导数分析(图2a2)

揭示了更清晰的动力学过程：

  • λ₁: 在g/gc∼0.8g/g_c \sim 0.8处出现显著负峰，表明快速降能
  • λ₂: 立即出现正峰响应，随后在λ₃介入后回落
  • λ₃: 在g/gc∼1.0g/g_c \sim 1.0处被激活，随后持续增长
  • 总能量: λ₂的二阶导数与总能量高度一致，证实其"探测者"角色

3. 第一激发态(图2b)

展现完全类似但稍微平滑的过程，证实机制的普适性

模式特性分析(图3)

本征能量演化(图3 a1-a3)

  • λ₁: 从-0.25平滑下降至-0.30
  • λ₂: 始终接近0，轻微下降
  • λ₃: 从+0.25上升至+0.30

波函数分量(图3 b1-b3)

模式λ₁和λ₃(二能级主导)：

  • un,2(σ^z)u_{n,2}(\hat{\sigma}_z)：保持主导地位且几乎不变(~0.9)
  • un,1(iσ^y)u_{n,1}(i\hat{\sigma}_y)：随耦合强度缓慢减小
  • un,3(a^)u_{n,3}(\hat{a})：光子分量持续增加

模式λ₂(光子主导)：

  • u2,3(a^)u_{2,3}(\hat{a})：始终主导(~1.0)
  • u2,1(iσ^y)u_{2,1}(i\hat{\sigma}_y)：随耦合强度增加
  • u2,2(σ^z)u_{2,2}(\hat{\sigma}_z)：始终可忽略

这解释了为何λ₂携带大部分超辐射光子(图5a)，而对自旋翻转贡献极小(图5b)

波函数空间分布(图4)

正常相(g/gc=0.5g/g_c=0.5)

  • λ₂主导总波函数，几乎完全重合
  • λ₁和λ₃被轻微激发但反相，相互抵消

临界点(g/gc=1.0g/g_c=1.0)

  • 波函数开始变得不对称
  • 三个模式权重发生显著变化

超辐射相(g/gc=1.5g/g_c=1.5)

  • 总波函数和模式分量高度不对称
  • λ₂仍然主导，但λ₁和λ₃权重大幅增加
  • 关键特征：λ₁与λ₃始终反相，λ₃与λ₂同相

物理量的模式分解(图5)

光子数⟨a^†a^⟩\langle\hat{a}^\dagger\hat{a}\rangle(图5a)

  • λ₂贡献：主导超辐射光子，从0增至~45
  • λ₁和λ₃贡献：极小且相互抵消
  • 与λ₁和λ₃主要由二能级构成一致

自旋翻转⟨σ^x⟩\langle\hat{\sigma}_x\rangle(图5b)

  • λ₁和λ₃贡献：主导自旋翻转，反相关系
  • λ₂贡献：几乎为零
  • 与λ₂主要由光子模式构成一致

定量发现总结

  1. 相变启动：λ₁在g/gc≈0.8g/g_c \approx 0.8开始主动降能
  2. 快速响应：λ₂在g/gc≈0.9g/g_c \approx 0.9立即响应
  3. 模式唤醒：λ₃在g/gc≈1.0g/g_c \approx 1.0被激活
  4. 角色转换：λ₃权重在g/gc≈1.1g/g_c \approx 1.1超越λ₂
  5. 新相稳定：g/gc>1.2g/g_c > 1.2后系统进入稳定超辐射相

相关工作

量子Rabi模型研究历史

  1. 基础理论 16,17：I.I. Rabi于1936-1937年提出，描述二能级系统与单模光场相互作用
  2. 早期分析 18,19：发现强耦合、超强耦合、深强耦合区域的压缩效应，是超辐射的前兆
  3. 相变理论 26-37：
    • De Liberato & Ciuti 26：多模玻色场与平带电子耦合的量子相
    • Hwang等27：量子相变和普适动力学
    • Ying等28：基态相图，提出临界点公式gc=1+1+Δ2/16g_c=\sqrt{1+\sqrt{1+\Delta^2/16}}
    • Liu等29：普适标度和临界指数

实验进展

  1. 超强耦合实现 20-25：
    • 2004年Wallraff等在超导量子比特中实现强耦合
    • 2010年Niemczyk等实现超强耦合
    • 2017-2018年Yoshihara等实现深强耦合和能级反转
    • 2020年Mueller等在等离激元纳米晶体中实现深强光-物质耦合
  2. 相变观测 38-40：
    • 2021年Cai等在单个囚禁离子中观测到量子相变
    • 2021年Chen等通过反压缩在NMR平台模拟超辐射相变

本文与前人工作的关系

  • 作者先前工作37：指出相变源于量子双阱势的形成
  • 本文突破：首次揭示相变的微观动力学机制——三模式协同演化过程
  • 方法论创新：算符空间对角化提供了全新分析视角，超越传统唯象理论

相关模型

论文指出该方法可推广至：

  • Dicke模型41：多原子集体与光场耦合
  • 自旋-玻色子模型42：研究退相干和耗散机制

结论与讨论

主要结论

  1. 微观机制揭示：首次从微观动力学角度回答了超辐射相变"为何"和"如何"发生：
    • 为何：λ₁模式主动降能突破能量势垒
    • 如何：λ₂快速响应→唤醒λ₃→λ₃平衡λ₁→新相稳定
  2. 三模式角色定位：
    • λ₁(initiator)：相变主动驱动者
    • λ₂(inspector)：敏感探测者和超辐射光子载体
    • λ₃(stabilizer)：被动稳定者
  3. 方法论贡献：算符空间对角化为研究量子相变提供了新范式，超越朗道理论和重整化群的唯象/描述性框架
  4. 普适性：该解剖机制可能适用于自然界广泛存在的相变现象

局限性

  1. 模型简化：
    • 仅考虑单模光场和单个二能级原子
    • 忽略了环境耗散和退相干效应
    • 零温近似(基态分析)
  2. 参数限制：
    • 主要分析了大失谐情况(Δ=50\Delta=50)
    • 未系统研究不同Δ\Delta值的影响
    • Fock基截断(N=200N=200)可能在极强耦合区域不足
  3. 机制验证：
    • 缺乏直接实验验证三模式动力学的方案
    • 未提供测量各模式贡献的实验可观测量
  4. 理论深度：
    • 未建立与重整化群理论的联系
    • 缺少对临界指数的解释
    • 未讨论有限温度效应

未来方向

论文明确指出：

  1. 模型推广：扩展至Dicke模型和自旋-玻色子模型
  2. 复杂物理：探索退相干和耗散机制的微观起源
  3. 实验设计：发展测量模式分解的实验方案
  4. 普适机制：寻找其他物理系统相变的类似模式结构

深度评价

优点

1. 方法论创新性(★★★★★)

  • 原创性：算符空间对角化是研究量子相变的全新视角，首次将哈密顿量分解为具有明确物理意义的基本模式
  • 严格性：分解在给定Fock基下精确，与数值精确对角化完全一致
  • 普适性：方法可推广至其他复杂量子系统

2. 物理洞察深刻(★★★★★)

  • 微观机制：首次回答相变"为何"和"如何"发生，超越传统唯象理论
  • 动力学图像：清晰描绘了三模式的时序演化和相互作用
  • 角色定位：initiator-inspector-stabilizer的三角色机制具有很强的物理直觉

3. 分析充分性(★★★★☆)

  • 多层次分析：从能量、导数、波函数、物理量多角度验证
  • 定量描述：提供了详细的数值结果和演化曲线
  • 自洽性：所有分析相互支持，形成完整图景

4. 表达清晰度(★★★★☆)

  • 逻辑严密：从方法介绍→验证→应用→结论，层次分明
  • 图表丰富：5个主图包含多个子图，信息密集但清晰
  • 物理图像：用"主动/被动"、"驱动/稳定"等直观语言描述抽象概念

不足

1. 理论完备性

  • 缺少解析理论：主要依赖数值计算，未建立解析理论框架
  • 临界行为：未讨论临界指数、标度律等重整化群关心的问题
  • 普适类：未明确该相变属于哪个普适类

2. 实验可验证性

  • 间接验证：三模式分解是理论构造，缺乏直接实验测量方案
  • 可观测量：未提出能区分三模式贡献的实验可观测量
  • 实验参数：Δ=50\Delta=50对应的实验条件未讨论

3. 适用范围

  • 单一参数：仅分析了Δ=50\Delta=50的情况
  • 零温限制：未考虑有限温度效应
  • 理想模型：忽略耗散、退相干等实际因素

4. 技术细节

  • Fock基收敛性：未系统研究NN的收敛性
  • 数值精度：未报告数值误差和精度分析
  • 计算效率：未讨论与传统方法的计算成本对比

影响力评估

对领域的贡献(★★★★☆)

  1. 范式转变：从"描述相变"到"解剖相变"，开辟新研究方向
  2. 方法工具：算符空间对角化可成为量子多体理论的新工具
  3. 理论深化：为量子Rabi模型研究提供了更深层次的理解

实用价值(★★★☆☆)

  1. 理论指导：为设计和控制量子相变实验提供微观图像
  2. 模型推广：方法可应用于量子光学、腔QED等领域
  3. 局限性：需要进一步发展才能直接指导实验

可复现性(★★★★☆)

  1. 方法清晰：算符空间对角化步骤明确
  2. 参数明确：所有计算参数均已给出
  3. 代码缺失：未提供开源代码，需自行实现

适用场景

适合应用的领域

  1. 量子光学：腔QED系统的相变研究
  2. 超导量子比特：电路QED中的超强耦合现象
  3. 冷原子系统：光晶格中的量子相变
  4. 量子化学：分子系统的激发态相变

适合研究的问题

  1. 相变机制：其他量子模型的微观动力学
  2. 临界现象：模式演化与临界指数的关系
  3. 非平衡动力学：淬火动力学中的模式演化
  4. 拓扑相变：拓扑不变量的模式分解

不适合的情况

  1. 强耗散系统：开放量子系统需要扩展方法
  2. 多体纠缠：大量粒子系统计算成本过高
  3. 非线性系统：超出二次型哈密顿量的系统

参考文献(精选)

相变理论基础

  • 8 Landau, L.D. (1937). On the theory of phase transitions. Phys. Z. Sowjet. 11.
  • 12 Wilson, K.G. & Kogut, J. (1974). The renormalization group and the ε expansion. Phys. Rep. 12, 75.

量子Rabi模型关键文献

  • 16-17 Rabi, I.I. (1936-1937). 量子Rabi模型原始论文
  • 27 Hwang, M.-J. et al. (2015). Quantum phase transition and universal dynamics in the Rabi model. Phys. Rev. Lett. 115, 180404.
  • 28 Ying, Z.-J. et al. (2015). Ground-state phase diagram of the quantum Rabi model. Phys. Rev. A 92, 053823.

实验观测

  • 38 Cai, M.-L. et al. (2021). Observation of a quantum phase transition in the quantum Rabi model with a single trapped ion. Nature Commun. 12, 1126.
  • 40 Chen, X. et al. (2021). Experimental quantum simulation of superradiant phase transition beyond no-go theorem. Nature Commun. 12, 6281.

总体评价：这是一篇在量子相变研究领域具有重要方法论创新的工作。通过算符空间对角化这一新颖方法，作者首次揭示了超辐射相变的微观动力学机制，回答了长期悬而未决的"为何"和"如何"问题。尽管存在理论完备性和实验可验证性方面的不足，但其开创性的思想和清晰的物理图像为量子相变研究开辟了新方向，具有较高的学术价值和潜在影响力。该工作适合对量子光学、量子相变和理论方法发展感兴趣的研究者深入研读。