编辑推荐：
本研究推导并实现了量子电动力学哈特里-福克（QED-HF）方法的分子梯度，并应用于水、水二聚体和Corannulene，揭示了不同外部电磁场（包括量子电动力学场）对分子几何结构和反应能垒的影响。通过解析梯度表达式和数值优化，发现静态电场和磁场分别导致分子键长、键角及构型取向的变化，而QED场通过调节偶极矩涨落显著改变分子平衡构型，例如水分子键角减小0.1度，Corannulene的翻转变能垒降低约30%。同时，双场耦合（如正交电场与QED场）可协同调控分子构型，例如水二聚体在QED场中形成更稳定的平面构型。该方法为实验中利用电磁场调控分子反应路径提供了理论支持，适用于光谱学分析和绿色燃料合成。
外部电磁场在化学研究中展现出巨大的潜力，它们可以改变分子结构和反应性，从而影响化学反应的路径和产物。近年来，理论和实验研究不断揭示了电磁场对分子系统的影响，特别是在分子几何结构、反应势垒以及分子取向等方面。这些研究不仅为预测分子性质的变化提供了理论基础，还为开发新的化学反应控制方法提供了实验支持。本文通过引入量子电动力学（QED）和静态电场，探讨了如何在哈特里-福克（Hartree-Fock）框架下推导并实现分子梯度，以研究这些外部场对分子平衡几何结构的影响。为了验证这一方法的可行性，我们选择了具有代表性的分子——水、水二聚体和冠烷（corannulene）——作为研究对象，并观察了它们在不同外部场下的几何结构变化。
在静态电场作用下，分子会表现出显著的取向依赖性。这种现象在许多研究中已经被证实，例如Zhu等人指出，强电场（通常在1.4 V/nm左右）可以显著提高卤素键合复合物的反应速率和选择性。而理论计算则支持这一发现，并解释了分子取向对异裂反应势垒的降低作用。此外，Cassone等人通过计算强静态电场对甲烷和甲醛的生成，提出了电场可能在绿色燃烧燃料的合成中发挥重要作用。而最近的实验则进一步展示了电场在常温下对水的相变，即所谓的“电冷冻”（electrofreezing）。
在磁性场的影响下，分子的电子结构也会发生变化。Wang等人研究了在磁性环境（1 T）中Cu₂O纳米立方体对CO₂还原形成甲酸的催化作用，发现磁性场可以显著提高催化活性。而极端磁性场（如在中子星周围出现的1000 T）则被用于研究分子化学变化，并对观测到的分子光谱进行解释。同时，腔量子电动力学（QED）因其能够产生极化子（polaritons）而受到关注，尤其是在分子与量子场之间有强耦合的情况下。实验研究已经表明，通过在腔镜上施加电压，可以实现对反应选择性和速率的调控。
为了更准确地研究外部场对分子系统的影响，理论计算是不可或缺的。其中，确定分子的平衡几何结构是任何理论研究的基础。外部场能够显著改变分子的平衡几何结构，从而影响反应焓和势垒。例如，Pavošević等人通过计算研究了外部场对分子几何结构的影响，并强调了几何弛豫（relaxation）在获得外部场下的平衡结构中的重要性。为了计算平衡几何结构，需要开发可靠的评估分子梯度的方法。在静态电场下，分子梯度可以使用哈特里-福克理论进行评估，但磁性场的情况更为复杂，因为要获得原点不变的能和梯度，需要使用伦敦原子轨道（LAOs）。
量子电动力学哈特里-福克（QED-HF）方法已经发展出多种理论模型，包括腔Born-Oppenheimer近似哈特里-福克（CBOA-HF）、量子电动力学密度泛函理论（QEDFT）和量子电动力学耦合簇理论（QED-CCSD）。本文展示了QED-HF分子梯度的推导和实现。QED-HF方法等价于CBOA-HF方法，但其中最小化了位移场的影响。为了适应外部场的计算，坐标系统需要进行调整。通常采用笛卡尔坐标进行计算，但这种方法可能会引入数值不稳定性。因此，采用内部坐标系统可以减少几何优化的循环次数并提高计算的稳定性。然而，为了考虑外部场的取向依赖性，必须包含旋转自由度。Wang和Song开发的翻译-旋转-内部坐标（TRIC）系统有效地解决了这一问题，提供了更稳定和高效的优化框架，并被用于上述研究中。
本文的结构如下：第二和第三部分分别展示了QED-HF分子梯度的详细推导和实现，第四部分则扩展了理论以包括外部场的影响，第五部分展示了水、冠烷和水二聚体的结果，第六部分提供了结论。QED-HF理论基于Pauli-Fierz哈密顿量，该哈密顿量在相干态基础上表达为电子哈密顿量与场的相互作用。在QED-HF方法中，分子梯度的推导与实现需要考虑量子场对分子几何结构的修改。通过引入新的积分形式，可以将量子场的影响纳入计算框架，从而实现对分子几何结构的优化。
在静态电场下，分子的几何结构会受到显著影响。例如，电场会使水的平衡几何结构发生偏移，而磁性场则会影响冠烷的翻转势垒。此外，水二聚体中的水分子会呈现正交于量子场方向的取向。这些现象表明，外部场不仅影响分子的几何结构，还影响其反应性。通过使用TRIC系统，可以在计算中同时考虑分子的取向和内部自由度，从而更准确地模拟外部场对分子系统的影响。
在研究水分子的几何结构变化时，我们观察到在QED场和静态电场作用下，水分子的键长和键角都会发生变化。而磁性场对水分子的影响则表现为键角的增加和键长的减少。这些变化在QED场和静态电场的相互作用下更加显著，因为QED场通过减少偶极子的波动来影响分子的几何结构。而静态电场则通过改变分子的取向来影响其几何结构。在QED场和磁性场的共同作用下，水分子的几何结构会进一步改变，表现出不同的取向和键长变化。
在冠烷的研究中，我们观察到在QED场和磁性场的作用下，其翻转势垒都会降低。这表明外部场能够有效改变分子的内部自由度，从而影响其反应路径。而在水二聚体的研究中，我们发现，在低耦合强度下，两个水分子之间的二面角变化较小。然而，在足够高的耦合强度下，全局最小值会发生变化，使得每个水分子在QED场下局部取向趋于平面化。这些结果表明，外部场，特别是与QED场的组合，能够以复杂和系统依赖的方式影响分子的几何结构。
总体来看，这些研究结果表明，外部场对分子系统的影响具有复杂性和多样性。在QED场和静态电场的共同作用下，可以实现对分子几何结构的精细调控，从而影响其反应性和选择性。而在磁性场的作用下，虽然也能改变分子的几何结构，但所需的场强度仍然超出当前实验技术的范围。通过这些研究，我们能够更好地理解外部场对分子系统的影响，并为未来的研究提供理论支持。这些发现不仅揭示了外部场对分子结构的调控作用，还为开发新的化学反应控制方法提供了基础。
