On the electronic path integral normal modes of the Meyer-Miller-Stock-Thoss representation of nonadiabatic dynamics

技术摘要：论 Meyer–Miller–Stock–Thoss 表示法中电子路径积分正规模的性质

问题陈述
准确且高效地模拟非绝热动力学对于理解复杂系统中的电荷和能量转移至关重要。验证此类方法有效性的一个关键准则在于是否能保持量子玻尔兹曼分布（QBD）的守恒。虽然 Matsubara 动力学通过截断虚时路径积分的高阶正规模，成功实现了单表面核动力学的 QBD 守恒，但将其扩展到非绝热系统（NA-Mats）时却面临挑战。现有的 NA-Mats 方法对核变量应用了 Matsubara 截断，但对 Meyer–Miller–Stock–Thoss (MMST) 表示下的电子变量则保留了完整的表示而不进行正规模截断。因此，该方法在处理一般系统时无法实现 QBD 守恒。这提出了一个根本性的问题：在 MMST 表示法中截断高阶路径积分电子正规模，是否会导致 QBD 守恒，还是说这些模式在本质上与其对应的核模式不同？

方法论
作者利用一个仅含电子系统的模型来研究 MMST 框架内电子正规模的性质，以隔离电子自由度的行为。

• 理论框架： 研究采用了 MMST 哈密顿量，该哈密顿量将多能级量子系统映射为耦合谐振子，并引入电子位置 ($q$) 和动量 ($p$) 变量。作者使用与核 Matsubara 正规模相同的变换矩阵 ($T$) 定义了电子正规模 ($\check{q}_k, \check{p}_k$)，从而使自由环聚合物弹簧常数矩阵对角化。
• Liouvillian 分析： 在正规模基组下推导了精确的电子 Liouvillian。对于核与电子自由度解耦的系统，研究表明正规模的演化方程是可分离的，且其形式与 bead 演化方程一致。
• 相关函数： 作者利用 bead 变量和正规模变量构建了广义 Kubo 变换 (KT) 相关函数。他们特别考察了观测量（如状态布居）是否可以表示为有限个低阶正规模的函数。
• 截断协议： 计算测试通过两种方式进行截断：
  1. Bead 截断： 仅传播前 $j$ 个 bead，而将其余 bead 保持在 $t=0$ 时刻。
  2. 正规模截断： 仅传播最低的 $k$ 个正规模，而将其余高阶模保持在 $t=0$ 时刻。
• 系统模型： 使用具有非对称势矩阵 ($V$) 的二能级电子系统，并在约化单位制 ($\beta=1$) 下进行 8-bead 环聚合物模拟。

核心贡献与结果
本研究是对 MMST 表示法中电子正规模的首次分析，得出了与传统核 Matsubara 动力学相比具有反直觉性的发现：

1. 非有限观测量： 不同于核观测量通常依赖于有限数量的低频 Matsubara 模式，MMST 表示中的电子观测量是所有正规模的函数。这意味着在进行截断时，变换矩阵无法被消除，因此观测量依赖于全套模式。
2. 缺乏分布收敛性： 在核 Matsubara 动力学中，高阶正规模的分布会变窄，从而允许对其进行积分处理。而在 MMST 电子表示中，高阶正规模并不受分布约束；它们的分布不会变窄，这导致在不损失精度的情况下无法对高频模式进行积分。
3. 单轨迹 QBD 守恒失效： 代数分析和计算结果表明，在 MMST 正规模中进行动力学截断并不能使单条轨迹实现 QBD 守恒。在多 bead 情况下的 QBD 守恒要求完成整个路径积分回路；截断破坏了这一条件。
4. 系综平均伪影： 虽然通过对轨迹进行系综平均时 QBD 似乎得到了保持（由于某种平均效应），但它在单个轨迹上并不守恒。这使得该方法区别于真正的 QBD 守恒动力学。
5. 相关函数不准确： 在 MMST 正规模中截断动力学会导致相关函数不准确，因为必须包含所有模式才能获得合理的动力学行为。

意义与主张
论文得出结论，MMST 路径积分正规模并非开发准确且满足 QBD 守恒的非绝热动力学方法的理想选择。作者认为，Matsubara 动力学实现 QBD 守恒所需的特定属性——即观测量需为有限数个模式的函数，以及高阶模式分布变窄——在 MMST 电子变量中并不成立。因此，简单地将 Matsubara 截断应用于 NA-Mats 框架中的电子变量，并不能产生一种能使单条轨迹实现 QBD 守恒的方法。这项工作表明，可能需要其他的映射表示或度量标准，才能同时实现 QBD 守恒、拉比振荡重现以及经典标度律的要求。