变分极化子理论：强耦合光-物质与电子-声子体系基态的变分处理方法

《The Journal of Physical Chemistry Letters》：Variational Polaron Theory for Ground States of Strongly Coupled Light–Matter and Electron–Phonon Systems

　　

强光-物质耦合与电子-声子耦合会产生由虚玻色子激发所修饰的基态，这使得在超强耦合区间内，裸态截断（bare-state truncation）与微扰处理方法不再可靠。研究人员提出了一种基于态依赖极化子变换（state-dependent polaron transformation）的非微扰变分基态框架，该框架结合了乘积态（product-state）拟设与针对残余物质-玻色子纠缠的二阶微扰修正。研究表明，在无穷耦合极限下，优化后的变换参考系渐近退耦，因为主导的线性耦合项被抵消，同时非对角物质跃迁被位移振子重叠积分（displaced-oscillator overlaps）所抑制。该方法在弱耦合与强耦合极限下均渐近精确，并且在固定极化子变换（fixed polaron transformation）最不可靠的中间区间仍保持较高精度。Dicke模型基准测试复现了基态能量、保真度（fidelity）以及超辐射相变（superradiant transition），二阶能量误差低于0.2%。Holstein模型基准测试的误差低于0.5%，并阐明了平移对称性对波函数质量的影响。该缀饰基（dressed-basis）框架能够对强耦合光-物质与电子-声子体系进行非微扰建模。

研究背景与问题

强耦合是分子科学、凝聚态物理与腔量子电动力学（cavity quantum electrodynamics, cavity QED）中的核心主题。在光学腔中，电子与振动激发可与受限电磁模式杂化形成分子极化子（molecular polaritons），从而改变分子光谱、能量流与化学反应活性。在分子与有机材料中，载流子可通过电子-声子耦合被局域振动模式强烈修饰，形成极化子（polarons）及极化子输运。这些体系的共同物理特征是：相关低能态并非裸物质态与玻色子态的简单乘积，而是包含虚光子或虚振动激发的缀饰态（dressed states）。这种缀饰是许多腔致与声子致分子和材料性质改变的微观起源，但也使基态问题本质上具有非微扰性。

在超强耦合与深强耦合区间，耦合强度与特征玻色子频率或物质激发能相当甚至超过后者，旋转波近似（rotating-wave approximation, RWA）不再可靠，反旋波过程直接贡献于基态，裸态截断可能低效或失控。弱耦合微扰理论因玻色子缀饰不再微小而失效，而传统强耦合极化子图像在中间区间可能定量不准确。因此，一种有用的基态理论应在弱耦合与强耦合之间插值，在任一极限描述都不充分的区间保持精确，并提供物质-玻色子缀饰的紧凑物理表示。

已有研究发展了指数拟设与变分Lang-Firsov变换、分子量子电动力学与腔QED分子轨道理论、变分极化子变换等方法。然而，许多变换参考系方法使用变换后的玻色子真空或相干态作为零阶态，未显式扫描残余玻色子希尔伯特空间，可能遗漏不完整缀饰变换后存留的非真空玻色子激发。此外，已有研究未系统利用变换哈密顿量在无穷耦合极限下的渐近退耦性质作为变分基态理论的组织原则。

研究内容与结论

研究人员引入了一种基于态依赖极化子变换的非微扰变分基态框架。该框架一般性地证明：优化后的变换参考系在无穷耦合极限下渐近退耦——优化位移抵消主导的对角线性物质-玻色子耦合，而非对角物质区跃迁被位移振子重叠积分抑制。渐近退耦不仅作为形式化的强耦合陈述，还作为紧凑缀饰基描述何时可靠的物理基础与定量诊断。有限耦合体系仍保留残余物质-玻色子关联，但强耦合极限表明适当优化的缀饰参考系应使这些关联微扰化。因此，研究人员从物质与玻色子区因子化的零阶态出发，通过二阶微扰修正纳入残余关联。在有限耦合下，位移参数通过最小化基态能量进行变分优化，使缀饰基在裸弱耦合描述与渐近退耦强耦合描述之间插值。

乘积态参考仅在物质与玻色子区之间因子化，玻色子分量不限于变换真空，也不假设在独立模式上因子化，而是在保留的多模玻色子希尔伯特空间中展开，并与物质分量及位移参数共同优化。残余物质-玻色子纠缠通过二阶微扰理论纳入。

研究人员用两个范式模型进行基准测试。Dicke模型作为集体光-物质耦合的最小模型，检验反旋波耦合、集体增强与超辐射量子相变。Holstein模型描述巡游电子自由度与局域振动模式的耦合，直接关联分子与有机材料中的极化子形成。通过与精确对角化（exact diagonalization）对比，研究人员量化了有效退耦区间的实际起点以及中间耦合区间的残余能量与保真度误差。

关键技术方法

研究人员采用的核心方法包括：态依赖变分极化子变换（state-dependent variational polaron transformation），通过位移参数λ_α,μ吸收主导玻色子缀饰；物质-玻色子乘积态拟设（product-state ansatz），其中玻色子分量在保留的多模Fock空间中展开并显式优化非真空分量；变分参数（物质系数、玻色子系数与位移参数）的联合梯度下降优化；以及围绕优化缀饰参考系的二阶微扰修正（second-order perturbative correction），用于纳入残余物质-玻色子纠缠。位移振子重叠积分通过Laguerre多项式公式解析求值。

研究结果

Dicke模型基准（N = 1）：零阶变分计算的最大相对能量误差约为1.2%，二阶修正将最大误差降至约0.2%。保真度从变分水平的最小约0.995提升至微扰修正后的超过0.999。误差在g→0与g→∞极限下均消失，仅在中间区间达到峰值。截断收敛数据显示，N_cut从16增至20时二阶最大误差仅从2.34×10^–3变为2.33×10^–3，表明结果已收敛。

Holstein模型基准：误差低于0.5%，并阐明了平移对称性对波函数质量的影响。非真空玻色子分量优化相比变换真空限制，将两种基准的最大二阶能量误差降低了一个数量级以上。

讨论与结论

该框架结合了渐近退耦缀饰表示、变分优化、残余玻色子希尔伯特空间的显式扫描以及残余纠缠的微扰修正，既避免了裸态截断，也避免了变换真空限制。乘积态结构具有计算优势：物质区原则上可用针对关联物质体系的问题特定方法独立处理，而玻色子缀饰在优化变换参考系中处理。该框架为分子极化子与极化子基态提供了紧凑且定量可靠的起点，在弱耦合与强耦合极限下均渐近正确，并在固定极化子变换最不可靠的中间区间保持精确。