Long-Range Electrostatics and Born Effective Charges from Invariant Moment Tensor Potentials with Latent Ewald Summation (MTP-LES)

技术摘要：基于具有潜在 Ewald 求和（MTP–LES）的旋转不变矩张量势的远程静电相互作用与 Born 有效电荷研究

问题陈述

机器学习原子间势（MLIPs）已实现了高精度和低计算成本，但其本质上受限于对局部短程描述符（通常为 4–6 Å 切断半径）的依赖。这种局部性假设无法捕捉对于模拟极性液体、铁电氧化物及有限分子系统至关重要的长程静电相互作用和极化响应。虽然最近的框架（如 Latent Ewald Summation, LES）已经证明，仅通过能量和力训练即可推导出每个原子的潜在电荷（per-atom latent charges），从而在没有显式电荷监督的情况下恢复 Born 有效电荷（BECs）和偶极矩等性质，但现有的实现方式仍局限于基于图或等变架构（例如 CACE-LR, MACE-LES, Nequip-LES）。这些架构显式地编码了方向信息，这引发了一个疑问：这类特征是否对于恢复具有物理意义的潜在静电响应是绝对必要的，或者关键在于如何约束和稳定电荷自由度。

方法论

作者引入了 MTP–LES 框架，该框架将 Latent Ewald Summation (LES) 形式论与利用旋转不变描述符的 矩张量势 (MTP) 相结合。与以往依赖于等变消息传递的 LES 实现不同，MTP–LES 采用矩张量的标量收缩。为了在这一不变描述符空间内实现潜在电荷的稳定学习，作者实施了三种特定的架构修改：

1. 独立分支：模型使用用于短程能量 ($E_{SR}$) 和潜在电荷预测 ($q_i$) 的独立神经网络分支。这防止了电荷自由度被“奴役”于局部能量表示中。
2. 物种相关偏置与中性投影：一个可学习的、物种相关的电荷偏置 ($b_{s_i}$) 为电荷通道提供具有化学合理排序的初始化。至关重要的是，模型通过将平均电荷投影为零 ($\sum_i q_i = 0$) 来强制执行每个构型的显式电荷中性。这消除了在无约束模型中出现的非物理残余晶胞电荷。
3. 可选偶极损失：对于有限分子系统（非周期性系统），由于不存在长程 Ewald 贡献，且仅靠能量/力训练可能不足以充分约束绝对电荷尺度，因此添加了一个可选的辅助损失项，用于将总分子偶极矩与参考值（通过 MBIS 电荷计算）进行匹配。

总能量定义为 $E_{total} = E_{SR} + E_{LR}$，其中 $E_{LR}$ 在周期性系统中通过 Ewald 求和计算，在非周期性系统中通过阻尼成对库仑求和计算。力与 BECs 分别通过对总能量和极化强度的自动微分获得。

关键结果

MTP–LES 框架在四种不同的材料系统中进行了评估：

• 铁电体 PbTiO$_3$：在基于 SCAN-DFT 数据的严格基准测试下，MTP–LES 实现了 $R = 0.994$ 的对角 BEC 相关性，与在相同数据上的 CACE-LR 基准线（$R = 0.993$）相匹配。模型成功恢复了异常的 Ti 和 O 电荷响应（例如，$Z^*_{Ti} \approx +6$），且无需任何电荷或 BEC 监督。虽然 BEC 量级的绝对 RMSE 高于 CACE-LR，但保留了张量相关性和异常的物种排序。
• 液态水：模型重现了与实验一致的结构相关性（氧-氧径向分布函数）和红外吸收光谱（librational, bending, and stretching modes）。添加 LES 模块后，与仅使用 MTP 的基准相比，力精度提高了约 21%。恢复的 BECs 处于其他已发表的 LES 增强基准（MACE-LES, Nequip-LES 等）范围内。
• 极性二肽：对于真空中的有限系统，引入偶极匹配损失稳定了潜在电荷尺度。模型在预测分子偶极矩（$R^2 = 0.996$）和 BECs（$R^2 = 0.93$）方面表现出高精度，优于缺乏偶极监督的消融实验。
• 铪氧化物 (HfO$_2$)：应用于一个新构建的单斜 HfO$_2$ 数据集，模型保持了结构保真度（重现了七配位的 Hf–O 结构），并恢复了异常的 Hf Born 电荷（$Z^* \approx +5$），此前并无针对该系统的 LES 基准。
• 计算效率：在液态水系统（高达 15,360 个原子）上的计时基准测试显示，LES 校正相对于仅使用 MTP 的基准仅引入了 3.4% 的墙钟时间开销，同时保持了随系统规模近乎线性的扩展性。

意义与主张

本文声称，准确的潜在电荷恢复并不完全取决于是否使用基于图或等变的架构。相反，关键因素在于电荷自由度是如何被暴露、约束和稳定的。通过解耦电荷与能量分支、强制电荷中性并利用具有化学合理先验的初始化，MTP–LES 证明了旋转不变的、基于描述符的势函数可以有效地模拟长程静电相互作用和极化响应。

这项工作扩展了 LES 框架在计算高效的描述符型架构中的适用性，为在 MLIPs 中引入长程静电相互作用提供了一条途径，且无需等变特征传播的高昂计算成本，也无需在训练期间使用显式的电荷标签。作者将 MTP–LES 定位为模拟极性、离子型和铁电材料的一种可行替代方案，在这些材料中，长程相互作用具有重要的物理意义。