《Advanced Science》：Field-Driven Proton Transport in Water Revealed by Machine Learning Molecular Dynamics
编辑推荐：
过量质子在水中的传输源于质子转移（PT）、氢键重排与溶剂极化之间的耦合。尽管平衡态下的格罗特斯（Grotthuss）扩散已被广泛研究，但强电场下其响应仍难以解析，因为模拟须同时捕捉化学反应性与诱导极化，并在长时间尺度上展开。研究人员为此发展了用于水体系水合过量
过量质子在水中的传输源于质子转移（PT）、氢键重排与溶剂极化之间的耦合。尽管平衡态下的格罗特斯（Grotthuss）扩散已被广泛研究，但强电场下其响应仍难以解析，因为模拟须同时捕捉化学反应性与诱导极化，并在长时间尺度上展开。研究人员为此发展了用于水体系水合过量质子的反应性与可极化机器学习原子间势（MLIP）。经微调的PiNet2模型在零场下为结构与扩散提供高精度参照，而基于AD(OS)的有限场耦合实现了外加电场下的非平衡模拟。在低于0.04 V/?的线性响应区，由漂移速度得到的质子迁移率与平衡扩散系数高度一致；质子漂移速度在高达0.15 V/?场内近似保持线性，而该场强下受约束水合氢离子的离子迁移率已出现下降。模拟显示质子化水表现出独特的场响应：氢键网络基本维持结构涨落与迁移能力，未出现可比拟的介电饱和特征。自由能分析进一步表明，强场可稳定具有四重氢键配位的水合氢离子构型。
研究背景与问题提出
水溶液中过量质子的传输是酸碱化学、电化学与能源转换中的核心过程。与常规离子依靠质心扩散不同，质子通过格罗特斯机制在波动的氢键网络中发生结构扩散，其电荷迁移由共价键重排、氢键拓扑变化、溶剂极化与核量子涨落共同控制。传统第一性原理分子动力学（AIMD）虽能刻画反应性与极化，但计算代价高昂，难以在长时间与大尺度上研究外电场驱动的非平衡质子动力学；经验价键或反应力场又常缺乏普适极化描述。电化学界面、双电层与限域环境中普遍存在强局域电场，而质子是否遵循与普通离子相同的介电饱和规律，此前并不清楚。因此，亟需一种兼具第一性原理精度与高效率、并能一致描述化学反应性与电场耦合的模拟框架。
研究人员开展的 работа 与总体结论
研究人员在《Advanced Science》发表工作，提出反应性与可极化机器学习分子动力学（MLMD）框架：以PiNNAcLe-FT流程将基础模型蒸馏为轻量PiNet2架构，再用收敛的DFT标签微调，得到零场反应性MLIP；随后引入原子偶极加氧化态AD(OS)模型，一致描述总偶极矩与玻恩有效电荷（Born effective charge），实现有限场非平衡模拟。研究表明：零场下模型准确复现水与质子结构及扩散；弱场中由漂移速度得到的质子迁移率与平衡扩散系数相互印证；质子漂移在较宽场强内保持近线性，区别于受约束水合氢离子的介电饱和；强场下纯水出现介电常数下降与结构有序化，而质子化水的氢键网络响应温和，且自由能面被重塑，四重配位水合氢离子构型被稳定。该工作为纳秒尺度电场驱动质子传输提供了实用计算路线。
主要关键技术方法
研究采用三个关键技术：其一，PiNNAcLe-FT工作流，将MACE-MATPES-r2SCAN-0基础模型蒸馏至PiNet2，并基于CP2K生成的收敛DFT标签微调，消除基础模型带来的过结构化偏差；其二，AD(OS)可极化耦合模型，从有限场DFT学习极化与玻恩有效电荷，满足现代极化理论，实现外电场下总偶极演化与受力计算；其三，非平衡MLMD漂移分析，通过监测总偶极位移提取质子漂移速度，并借爱因斯坦关系（Einstein relation）由漂移迁移率反推自扩散系数，同时结合质子转移自由能面（FES）分析氢键配位数与垂直能隙。
研究结果
2.1 质子化水体系的反应性MLIP
研究人员通过PiNNAcLe-FT构建并微调PiNet2势函数。零场径向分布函数显示，先验模型继承基础模型的过结构化液体网络，而微调后模型与匹配级别的AIMD参照一致，能量与受力验证误差分别低于1 meV/原子与31 meV/?，为后续场驱动模拟建立可信零场基准。
2.2 平衡态质子动力学
基于生产轨迹的平均平方位移分析表明，先验模型低估水自扩散，微调模型在330 K下氧自扩散与质子扩散分别接近实验值，质子与水扩散系数之比与热浴环聚合物模拟结果相当。升温近似核量子效应的策略可有效恢复传输动力学。
2.3 有限电场下的质子迁移
在外加电场中监测总偶极随时间增长得到漂移速度。低于0.04 V/?线性响应区内，由漂移速度经爱因斯坦关系得到的质子扩散系数与平衡MD结果吻合；漂移速度对场强近线性可持续至0.15 V/?。对照实验中固定所有O–H键、禁止质子转移，约束水合氢离子漂移速度在高场出现饱和，呈现单价离子典型的介电饱和行为，说明过量质子的结构扩散使其绕开了普通离子迁移率下降机制。
2.4 纯水与质子化水在场中的结构响应
纯水在场强增大时静态介电常数迅速下降，四面体序与O–O结构有序度上升，水自扩散降低，符合常压下电场诱导类冻结趋势。质子化水则不同：O–O径向分布与四面体序仅轻微改变，水自扩散基本维持甚至沿场方向略增，氢键数略减。高场下系统存在长寿命亚稳态导致结构异质性，但MLIP误差评估未见退化。自由能面分析显示，低场与零场相似；高场下沿场方向的质子转移态被优先占据，场耦合项稳定质子位于右侧氧的状态。强场还稳定了零场亚稳的四重氢键配位水合氢构型，水合氢离子与供体水分子偶极沿场对齐，使四配位环境在动力学上可被观测。
讨论与结论翻译
讨论指出，高场结构异质性源于亚稳态，尚需后续研究量化；质子漂移速度误差随场强放大，与系统尺寸和等效质子浓度有关。结论如下：研究人员发展了面向平衡与有限场条件的反应性与可极化MLMD框架，通过微调PiNet2与AD(OS)极化及玻恩有效电荷模型，在同一动力学描述中捕捉质子转移、氢键重排与外场耦合。平衡下微调模型校正过结构化并复现扩散基准；弱场下漂移迁移率与平衡扩散系数定量一致。质子漂移在0.15 V/?内近线性，而约束水合氢离子因介电饱和迁移率下降。质子化水结构响应弱于纯水，强场重塑质子转移自由能面并稳定四重配位构型。该框架对电解液、水电离、电极/电解质双电层、水-气界面及限域水中的场驱动离子传输具有普遍意义。
论文解读已完成，接下来可以选择一个方向继续深化：
