Non-adiabatic Ehrenfest dynamics with norm-conserving and ultra-soft pseudo-potentials with nuclear velocity corrections on the atomic orbitals within the Projector Augmented Wave Method framework

本文在投影增强波（PAW）方法框架内，推导出了一个具有伽利略不变性的第一性原理埃伦费斯特分子动力学框架，该框架通过在原子轨道上引入与核速度相关的相位，消除了对于规范守恒和超软伪势产生的虚假非绝热耦合。

要点

想象一下，你正试图预测一群人（电子）如何在围绕着一群舞者（原子核）的房间里移动。在量子物理的世界里，我们使用复杂的数学来模拟这场舞蹈。通常情况下，科学家们假设舞者是静止不动的，而人群在他们周围移动。但实际上，舞者一直在移动、旋转并改变位置。

这篇论文探讨了一个特定的问题，即当我们尝试模拟那些舞者开始移动时会发生的情况。

问题所在：“幽灵”移动

当科学家模拟原子移动时，他们经常使用一种被称为“伪势”（pseudo-potential）的捷径。你可以把它想象成使用一张简化的地图，而不是详细的卫星照片。这可以节省大量的计算能力。

然而，使用这些地图的旧方法存在一个缺陷。当“舞者”（原子核）以恒定速度移动时，模拟有时会错误地显示“人群”（电子）突然跳到新的能量状态或改变其行为。

论文称这种现象为违反了伽利略不变性（Galilean invariance）。用日常生活的语言来说，这就像是如果你在一列匀速行驶的火车上，你杯中的咖啡相对于你应该是静止的。但旧的模拟却说，仅仅因为火车在移动，咖啡就会突然溅出来。这在现实世界中是不合理的，但由于数学逻辑出现了偏差，从而产生了本不该存在的“幽灵”运动。

解决方案：“移动步道”

作者通过改变描述电子的方式解决了这个问题。

在旧的方法中，他们将电子视为被“粘”在舞者的位置上。如果舞者移动，电子的“家”只是僵硬地平移到了新位置。

在这一新方法中，作者为电子增加了一个特殊的“速度因子”。想象一下，电子不仅仅是坐在舞者身上，它们还乘坐着一条与舞者速度完全一致的移动步道。

• 相位偏移（The Phase Shift）： 他们增加了一种取决于原子核移动速度的数学“相位”（一种时间调整）。
• 结果： 现在，当原子核移动时，电子会完美地随之移动，就像移动步道上的乘客一样。这消除了“幽灵”运动。现在的模拟遵循了这样一个规则：恒定的运动不应该引起系统突然且无法解释的变化。

两种类型的地图

论文研究了制作这些简化地图（伪势）的两种不同方式：

1. 保范数型（Norm-Conserving，标准地图）： 这是较简单的版本。作者发现，通过添加“移动步道”的速度因子，完全解决了这个问题。数学变得简洁清晰，且“幽灵”作用力消失了。
2. 超软型（Ultra-Soft，灵活地图）： 这是用于重原子的更复杂、更灵活的版本。在这里，修复过程更为棘手。作者发现，他们不仅需要考虑原子核的速度，还需要考虑加速度（即原子核加速或减速的快慢）。
   • 他们发现，如果一个原子核在加速，它会对电子产生一个微小的“推力”（就像汽车加速时，你会感觉到被推向椅背的感觉）。
   • 旧的数学忽略了这个推力。新的数学包含了它，确保即使在原子加速或减速时，模拟依然保持准确。

根据论文，为什么这很重要

作者不仅仅是修复了一个漏洞；他们恢复了模拟中基本的物理定律。

• 不再有悖论： 他们证明了，如果你移动整个系统，电子不应该突然跳转到新的状态。他们的新方法确保了这种情况不会发生。
• 更高的准确度： 通过包含这些速度和加速度的调整，这个“简化的地图”（伪势）现在的表现与“详细的卫星照片”（全电子计算）完全一致，且不需要消耗那么多计算资源。

总结

这篇论文为模拟移动原子提供了一套新的规则。这就像是升级了电子游戏的软件，使得当角色奔跑时，物理引擎不会出现故障。通过为电子添加“速度调整”并考虑“加速度推力”，作者确保了关于原子和电子如何相互作用的模拟在物理上是正确的——无论原子是在匀速巡航，还是在加速或减速。

技术摘要

技术摘要：PAW 框架下包含核速度修正的非绝热 Ehrenfest 动力学

问题陈述
标准的基于第一性原理的非绝热分子动力学，特别是在 Ehrenfest 框架内，通常存在违反伽利略不变性（Galilean invariance）的问题。这是因为有效电子哈密顿量通常使用随原子刚性平移的类原子轨道构建，忽略了与核速度相关的相位（即电子-平移因子，ETFs）。当原子核运动时，这种忽略会导致伪非绝热耦合，其表现形式是悖论性的：由原子的恒定速度运动诱发了电子跃迁。此外，这种方法破坏了涉及 Born 有效电荷、原子间力常数与频率相关电磁极化率之间的求和规则。虽然这些问题在全电子计算以及特定的上下文（如振动圆二色性 VCD）中已得到解决，但目前使用伪势实现的基于第一性原理的非绝学分子动力学通常忽略了这些与速度相关的相位。

方法论
作者在投影缀加平面波（PAW）方法框架内，推导了一种第一性原理 Ehrenfest 分子动力学形式体系，其中明确将核速度相关的相位纳入了原子轨道基组。该方法论流程如下：

1. 包含速度的原子轨道： 作者为运动的原子核建立了定域原子轨道基组。通过求解运动原子核的孤立原子薛定谔方程，他们推导出原子轨道必须包含一个相位因子 $e^{i\alpha_s(\hat{r})}$，其中 $\alpha_s(\hat{r}) = \frac{m}{\hbar}\dot{R}_s(t) \cdot (\hat{r} - R_s(t))$。该相位解释了瞬时核速度 $\dot{R}_s(t)$。
2. PAW 变换推广： 标准 PAW 变换算符 $\hat{T}$（用于将伪波函数映射为全电子波函数）被推广为 $\hat{T}_{\dot{R}}$。这一新算符通过使用上述包含速度的原子轨道替换了静态原子轨道。
3. 拉格朗日形式体系： 作者采用了半经典 Ehrenfest 拉格朗日形式体系。至关重要的是，由于拉格朗日量依赖于核加速度（源于速度相关相位的时导数），其运动方程是使用适用于高阶导数的广义 Euler-Lagrange 方程推导出来的。
4. 有效哈密顿量的推导： 通过对全电子哈密顿量应用广义 PAW 变换来构建有效哈密顿量。作者仔细计算了变换算符 $\hat{T}_{\dot{R}}$ 的时间导数，这引入了依赖于核速度和加速度的附加项。

主要贡献与结果
本文为范数守恒（norm-conserving）和超软（ultra-soft）伪势提供了明确的推导：

• 范数守恒情况： 推导出的有效哈密顿量在非定域部分具有类似于 Peierls 相位的项，该项依赖于核速度（$e^{i\frac{m}{\hbar}\dot{R}_s \cdot \hat{r}} v^{nl}_s e^{-i\frac{m}{\hbar}\dot{R}_s \cdot \hat{r}}$）。在这种情况下，核加速度项消失或可以忽略不计，且结果与之前的范数守恒势研究结论一致（文献 [14]）。原本存在于刚性平移轨道方案中的伪耦合被消除，从而恢复了伽利略不变性。
• 超软情况： 作者识别并推导了超软伪势特有的额外项。这些项同时依赖于核速度和核加速度。它们源于 PAW 变换及增广电荷（$Q_{ij}$）的时间导数。具体而言，哈密顿量包含了涉及位置算符与非定域势对易项以及与加速度相关力（$m\ddot{R}_s \cdot \hat{r}$）相关的项。这些贡献在以往的文献中被忽视了。
• 运动方程： 本文给出了电子动力学（薛定谔方程）和核动力学（力）的完整方程集。核力包括标准的 Hellmann-Feynman 项、速度相关修正项以及新的加速度相关项。
• 守恒能量： 利用适用于依赖于加速度的拉格朗日量的 Ostrogradsky 构造，推导出了一个守恒能量表达式。该表达式包含了与哈密顿量及重叠矩阵对核速度和加速度的导数相关的修正项。

意义与主张
作者声称，其工作实现了使用伪势描述伽利略不变的非绝热 Ehrenfest 分子动力学。通过在原子轨道基组中引入核速度相关相位，该方法：

1. 消除了伪耦合： 它消除了当整个系统进行恒定速度全局平移时发生的虚假电子跃迁，而这种缺陷是标准刚性平移方法的固有缺陷。
2. 恢复了物理一致性： 它确保了伪势哈密顿量能够正确重现运动原子核下的全电子哈密顿量性质，从而恢复了频率相关的求和规则和正确的振动响应。
3. 推广至超软势： 该工作将对速度相关效应的理解从范数守恒势扩展到了超软势，强调了在超软情况下，核加速度项是不可忽略的，并且必须包含这些项才能获得准确的动力学振动响应（例如，力常数矩阵和 Born 有效电荷）。

文章结论指出，虽然哈密顿量的定域部分不受影响，但在 PAW 方法中，将这些相位纳入非定域部分和增广区域对于准确描述凝聚态系统中的非绝热过程至关重要。