非绝热动力学混合半经典模拟的最佳实践：可观测物理量决定量子化自由度和方法选择

《Physical Chemistry Chemical Physics》：Best practices for hybrid semiclassical simulations of nonadiabatic dynamicsOpen Access

　　

混合半经典（SC）方法通过选择性量子化关键自由度，为非绝热动力学模拟提供高效途径，并在降低计算成本的同时保留量子极限SC方法的大部分精度。本工作建立基于系统特征和关注的可观测物理量选择合适混合SC方法的实用指南。研究人员比较若干用于非绝热动力学的混合SC方法，包括混合量子-经典初始值表示（MQC-IVR及其变体）、部分线性化密度矩阵（PLDM）方法以及近期提出的非绝热混合Wigner动力学（NHWD）。使用一维（1D）散射模型和跨越弱与强系统-浴耦合区的自旋–玻色模型，研究人员以数值精确量子动力学和经典极限线性化半经典（LSC）近似为基准对这些SC方法进行基准测试。研究人员表明，混合方法的最佳选择及需量子化的自由度取决于关注的可观测物理量：电子可观测物理量，如态布居，仅量子化电子态变量即可准确捕获，而核可观测物理量则需要对相关核模式进行显式量子化。研究人员进一步发现，NHWD对强电子-核耦合系统特别有效，此情形需量子化电子和核自由度。最后，研究人员证明使用Meyer–Miller–Stock–Thoss（MMST）映射哈密顿量的对称化形式可一致地提高混合SC模拟的精度。

非绝热过程的高效理论模拟对理解光化学反应、金属表面通过电子–空穴对激发的振动能量耗散、生物体系电子转移以及锥形交叉附近动力学十分重要。数值精确方法如准绝热路径积分（QUAPI）、层级运动方程（HEOM）和多组态含时Hartree（MCTDH）仍受限于系统–浴模型哈密顿量和较少的系统自由度；近似虚时路径积分和表面跳跃方法虽有用，但并非由实时相关函数的精确量子表达式导出，不能通过受控纳入量子效应系统提高精度。基于半经典（SC）理论的方法利用扩展相空间中经典轨迹系综计算实时量子相关函数，可捕获非绝热跃迁、深隧穿和核相干。然而，量子极限SC方法如双Herman–Kluk初始值表示（DHK-IVR）计算昂贵，仅适用于低维系统；经典极限的LSC和对称准经典（SQC）方法虽广泛应用，但不能保持零点能（ZPE）并捕获干涉效应。混合方法通过对子集自由度采用量子极限处理，在保持计算效率的同时系统改进经典极限模拟，但尚缺乏针对给定可观测物理量选择合适混合SC方法及量子化自由度的系统框架。为此，研究人员比较量子极限、经典极限和混合方法，使用一维散射模型和凝聚相自旋–玻色模型，所有比较方法均采用MMST映射哈密顿量描述电子态动力学。

关键技术方法概述：研究人员采用一维散射模型（电子态与单一核自由度强耦合）和自旋–玻色模型（不同温度、对称/非对称电子态、弱/强系统–浴耦合；5组参数，4或50个浴模），比较MQC-IVR、AMQC-IVR、可分离前因子近似AMQC SP、PLDM、NHWD(E)与NHWD(V_n)，并以DHK-IVR、Husimi-IVR、LSC及精确量子结果为基准；MMST映射哈密顿量采用常规与对称化形式；精确基准采用离散变量表示（DVR）及小矩阵路径积分（SMatPI）/QUAPI。

4.1 Choosing the quantum subsystem in MQC（MQC中量子子系统的选择）：通过高入射能一维散射模型比较最终核动量分布与电子态布居，发现核动量分布需量子化核自由度，而电子态相关函数仅量子化电子变量即可准确获得；NHWD(E)能准确计算态布居，但因仅量子化电子变量而不能捕获核动量分布。由此得出，混合量子–经典模拟中选择量子化子系统由可观测物理量决定，而非系统变量的固有属性。

4.2 Choosing the form of the electronic state operators（电子态算符形式的选择）：通过比较σ_z算符的单激发振子（SEO）形式与连续映射变量（MV）形式在MQC-IVR中的收敛性，发现二者收敛至相同数值结果，但SEO形式在相同轨迹数下最大标准差更小、所需轨迹显著更少；研究人员将其归因于SEO形式更好地将映射变量约束在时间零点和后续时间的正确量子相空间。

4.3 Choosing a mixed quantum-classical SC method（混合量子–经典SC方法的选择）：通过比较高/低入射能下σ_z(t)相关函数，发现MQC-IVR（电子c_el=0.01、核c_N=5）和AMQC-IVR与精确结果吻合较好，而AMQC SP在强电子–核耦合模型中显著不准确；LSC精度可比，短时甚至更好，说明强耦合模式宜在同一理论水平处理。低入射能态分辨透射/反射布居中，MQC-IVR适度核量子化和AMQC-IVR与精确结果最接近；仅量子化电子变量的PLDM和NHWD(E)因近似电子–核相互作用而不准确，LSC显著较差，符合可观测物理量性质决定量子化程度。

4.4 Choosing the form of the MMST Hamiltonian（MMST哈密顿量形式的选择）：通过改变MMST哈密顿量形式比较PLDM和NHWD(E)的σ_z相关函数，发现二者结果显著不同，对称化形式一致产生更准确结果；LSC使用对称化形式略有改善，DHK-IVR结果不变。

4.5 Choosing a method for high-dimensional systems（高维系统方法的选择）：在自旋–玻色模型(a)–(d)中比较PLDM、LSC、NHWD(E)和AMQC SP，发现NHWD(E)与PLDM最准确，LSC一般较差，AMQC SP最不准确；弱系统–浴耦合且可观测为电子算符时，NHWD(E)或PLDM（结果相同）为最佳选择。强系统–浴耦合模型(e)中，准确电子态相关函数需采用NHWD(V)同时量子化强耦合浴模与电子变量。投影位置相关函数中，NHWD(V₁)选择性量子化被观测核模式较NHWD(E)和LSC-IVR显著改善；观察强耦合模式1时，NHWD(V₂)量子化模式1（Rabi频率模）和模式3（半Rabi频率模）最佳；观察弱耦合模式4时，量子化模式4与模式1仅边际改善，完全捕获其结构需量子化所有更强耦合模式及被观测模式。50浴模时，NHWD(V₁)和NHWD(V₂)捕获振荡行为，LSC快速阻尼，NHWD(E)较长时振荡频率不准确；10 a.u.、10⁵轨迹下NHWD(V₂)运行约7.38 min、标准差0.05，DHK约3139 min、标准差2.42且不能收敛50浴模，表明NHWD(V₂)兼具DHK精度和低计算成本，使高维系统模拟可行。

结论部分表明，研究人员建立非绝热混合量子–经典SC模拟的最佳实践：量子化模式的选择强烈依赖可观测物理量的性质；电子可观测物理量通过量子化电子态变量即可准确描述，而核可观测物理量或涉及电子与核算符的可观测物理量需显式量子化相关核自由度。NHWD(E)和PLDM是弱电子–核耦合且可观测为电子算符系统的最有效方法；NHWD(E)一般更优，因其系统推导且可按需量子化额外模式改进。NHWD(V_n)对显著振动耦合系统最有效，并证明对高维系统准确高效。在MMST框架内，对称化MMST哈密顿量提高精度且无额外计算成本；电子态算符采用离散SEO形式较MV形式改善数值收敛。未来工作将聚焦消除NHWD模拟量子化子系统高维时的剩余瓶颈，包括SC前因子的额外近似、探索电子自由度的不同映射变量方案以及纳入MQC-IVR的相位过滤思想。