自适应QM/MM不再在结构与能量间二选一：神经网络抵消切换力，水合结构无伪峰且能量守恒
核心概要
自适应多尺度模拟长期面临一个两难：溶剂分子在不同分辨率间切换会引入非物理的切换力，要么牺牲水合结构，要么破坏哈密顿量守恒。集成分子模型（IMM）用一个只以多尺度权重函数为输入的多层感知机学习抵消这些切换力，在纯水、K+、Na+、NH4+、Ca2+、Cl–水溶液和H3O+体系中给出无伪峰的水合结构，并在NVE系综下长期稳定守恒能量。修正力每步仅占约1.8%的计算时间，且在K+上训练的模型可直接用于Na+和NH4+而不需重新训练。
深度剖析
自适应QM/MM模拟中，同一溶剂分子在不同分辨率间切换会带来非物理的切换力，使模拟只能在结构失真与能量漂移之间取舍。 此前的势导向方法保留切换力，水合结构出现约6.0 Å处的非物理尖峰和过深的间隙极小值，模拟约20 ps即崩溃；力导向方法人为去掉切换力，结构较合理但需要非保守的记账项，能量持续漂移。 在纯水体系（2138个分子）中，势导向SCMP的径向分布函数出现上述伪峰，力导向SCMP的记账项单调上升，而集成分子模型在NVE系综下长时间运行无能量漂移。
修正势由一个多层感知机给出，其输入只是各分区的多尺度权重函数，而不是原子坐标。 权重函数在纯QM和纯MM体相区梯度为零，因此修正力被严格限制在过渡区，不会泄漏到体相，并自动满足牛顿第三定律与动量守恒。 网络为两层隐藏层（512和256个节点）、SiLU激活，用Huber损失和迭代“训练—模拟”主动学习采样，纯水4轮、K+与H3O+体系5轮后相邻轮次RDF的Jensen-Shannon散度低于1.0×10⁻⁴。
在K+水溶液上训练好的修正势可直接用于其他一价阳离子，无需重新训练。 此前每换一种溶质就要重新训练或重新构造修正函数；这里Na+和NH4+的RDF在零样本条件下消除了边界伪影，并与力导向参考一致。 同一模型用于Ca2+和Cl–时边界伪影重新出现，说明零样本迁移要求目标溶质与训练数据具有相似的宏观溶剂化结构。
消除切换力后，H3O+的扩散系数可以在不依赖恒温器的NVE轨迹上直接测得。 力导向SCMP的扩散系数随Langevin摩擦系数变化，γ0=100 ps⁻¹时为(5.7±1.0)×10⁻¹ Ų/ps，γ0=10 ps⁻¹时为(10.5±1.1)×10⁻¹ Ų/ps。 集成分子模型在纯NVE下给出(6.3±1.3)×10⁻¹ Ų/ps，作者将其略低于实验值9.4±0.1×10⁻¹ Ų/ps归因于原子核的经典处理和DFTB3对氢键的过强结合。
启示与展望
验证集中在纯水、含单个K+、Na+、NH4+、Ca2+、Cl–或H3O+的水溶液，高精度区用DFTB3/3OB，低精度区用SPC-Fw柔性水模型，SCMP取N=80个分区、每分区32个QM溶剂分子。纯水体系为2138个分子，离子体系为2137个水分子、边长40.00 Å的立方盒，H3O+体系为2047个水分子、边长39.48 Å。作者指出N=80是保守选择，此前工作表明N=40或48也可用于水合体系。
网络只以权重函数为输入，是对完整原子坐标的降维投影，瞬时残余力不为零，作者论证其在给定描述符下的条件系综平均为零。K+体系未做纯QM参考，因为单个离子的RDF需要数十皮秒采样，而DFTB3水模型在NVT/NVE下会出现空洞、在NPT下密度偏高。H3O+扩散系数略低于实验值，作者归因于经典核处理和DFTB3的过强氢键。零样本迁移对Ca2+和Cl–不成立。