光腔强耦合量子电动力密度泛函理论：分子结构与相互作用调控新范式

《Journal of Chemical Theory and Computation》：Quantum-Electrodynamical Time-Dependent Density Functional Theory Description of Molecules in Optical Cavities

本文系统介绍了量子电动力时间依赖密度泛函理论（QED-TDDFT）的张量积表示方法，为研究光学腔强耦合分子的电子-光子关联效应提供了计算高效的理论工具。研究通过结合实空间电子波函数与截断Fock空间光子态，成功预测了腔场对分子结合能、几何结构和极化子谱的显著调控作用，为光控化学反应和分子组装等前沿领域奠定了理论基础。

　　

研究背景与方法概述

光学腔量子电动力学（QED）是通过将光与物质限制在特定电磁环境中，精确调控两者相互作用的基础物理框架。腔内强耦合光-物质相互作用已发展成为调节物理与化学性质的重要研究方向，并在实验中受到广泛关注。描述这类耦合体系在理论与计算上充满挑战，因为需要同时处理电子的多体问题和光子自由度。为此，多种计算方案应运而生。

其中，基于密度泛函理论的量子电动力学密度泛函理论（QED-DFT）成为连接量子光学和电子结构理论的桥梁，为精确描述光腔中的光-物质强耦合现象提供了可行路径。然而，标准电子交换关联泛函（XC）无法充分描述强耦合下的电子-光子关联，这导致对极化子能级、基态修正以及光子分布统计的预测不准确。因此，开发专门的QED交换关联（QED-XC）势至关重要。

本工作提出的量子电动力学时间依赖密度泛函理论的张量积表示方法，其核心是将耦合的电子-光子波函数构建在实空间网格与Fock态表示的张量积上，形成了所谓的QED-DFT-TP及QED-TDDFT-TP框架。这种方法采用与波函数方法类似的耦合电子-光子波函数，并利用Fock态的正交性保持了实空间DFT哈密顿量的稀疏结构，从而能够采用传统实空间DFT中成熟且计算高效的迭代对角化技术。

形式主义

本研究考虑的所有体系均为非相对论性的，因此光-物质耦合可以通过Pauli-Fierz（PF）非相对论QED哈密顿量一致地描述。在偶极近似下，该哈密顿量由动能、Kohn-Sham势和光子哈密顿量三部分组成。其中，Kohn-Sham势包含了描述光与物质相互作用的势能项，而一个独立的类麦克斯韦方程则用于描述光子。通过引入与耦合强度λ相关的参数，可以清晰地展现顺磁耦合项和抗磁耦合项（又称海鸥项）对体系的影响。

耦合系统的轨道是在实空间和Fock空间的张量积上定义的。在KS框架下，每个分子轨道都与代表量化光子模式的特定Fock基态配对。哈密顿量的光-物质相互作用部分控制着不同光子态之间轨道元素的耦合。基态计算遵循传统的DFT方法，而为了计算吸收光谱，则采用在基态波函数上施加弱瞬时扰动（δ-kick），然后进行实时传播，通过记录偶极矩的演化并进行傅里叶变换来获取。

结果分析：基准测试与应用示例

研究的首要目标是在多个分子和分子二聚体上，将QED-DFT-TP方法与其他已建立的方法进行基准测试。计算确保了所有方法的参数一致性，包括耦合强度λ、腔频率ω和分子几何结构。结果显示，对于LiH分子，QED-DFT-TP与光子数QED-FCI（PN-QED-FCI）的计算结果非常吻合，仅在解离极限处存在微小差异。研究进一步探讨了耦合强度和腔频率对LiH基态能量以及|1?光子Fock态占据数的影响，发现增强耦合会导致基态能量因抗磁项贡献而升高，而增加腔频率则呈现相反的效果。

对于BH₃分子，通过实时传播计算吸收光谱并提取Rabi分裂，与QED-FCI方法进行了比较。虽然QED-TDDFT-TP始终低估了分裂的大小且表现出更强的非线性行为，但两种方法都显示出分裂随λ值增加而增大的趋势，整体一致性令人鼓舞。

针对分子间相互作用，研究了(H₂)₂、Ar₂和(H₂O)₂等二聚体体系。结果表明，腔限制能够以依赖于偏振方向的方式显著改变分子间的结合能和几何结构。例如，对于(H₂)₂，沿分子间轴偏振的腔模式比垂直偏振的模式产生更低的结合能。Ar₂的计算显示，当腔偏振垂直于二聚体轴且耦合较弱时，结合能几乎没有变化；而当偏振平行于轴时，结合能随耦合强度增加而显著降低。对于水二聚体，腔效应显著削弱了氢键，结合能下降幅度远超耦合簇CCSD方法的预测，突显了在强耦合区光对氢键作用的显著调控潜力。

最后，考察了HF二聚体在不同分子排列（平行与反平行）和腔偏振方向下的行为。对于反平行排列的HF-FH构型，光-物质耦合降低了两种偏振方向下的结合强度，并且平行偏振缩短了平衡距离，而垂直偏振则增大了平衡距离。|1?光子态的占据数在反平行构型中显著增强，且其分布展现出与偏振相关的特性。

总结

综上所述，本研究发展的QED-DFT-TP/QED-TDDFT-TP方法为在光学腔强耦合条件下模拟分子系统提供了一个准确且可扩展的计算工具。通过对一系列小分子体系的基准测试和应用展示，该方法在预测基态能量、极化子光谱以及腔修饰的分子间相互作用方面表现出良好的性能。研究工作证实了腔约束能够深刻影响分子的结合性质，其效应强烈依赖于腔模的偏振方向，这为未来利用光学腔主动调控化学反应路径、分子组装和材料特性开辟了新的理论途径。