Dynamic heterogeneity in sodium silicate melts via machine-learning potential

技术摘要：通过机器学习势函数研究钠硅酸盐熔体中的动态非均匀性

问题陈述
碱金属掺杂二氧化硅玻璃中的离子扩散是一个关键现象，其应用领域涵盖了从化学物理学到全固态二次电池的广泛范围。虽然“优先路径”（preferential pathway）图景表明碱金属离子是通过空间上截然不同的通道进行迁移，而非进行均匀扩散，但从弛豫动力学角度对这些过程进行详细表征的研究仍然有限。以往的研究主要依赖于具有各向同性相互作用的简单填充型玻璃模型，或缺乏捕捉网络玻璃复杂局部环境精度的经验势函数。相反，第一性原理计算（从头算分子动力学，AIMD）虽然具有高精度，但在观察玻璃转变附近的动态非均匀性和结构弛豫时，由于需要较大的系统尺寸和较长的时间尺度（纳秒级），其计算成本随 $O(N^3)$ 比例缩放，因此在计算上是难以实现的。

方法论
为了克服这些计算障碍，作者采用了机器学习原子间势函数（MLIPs）对钠硅酸盐熔体进行分子动力学（MD）模拟。

• 势函数构建： 构建了一个基于 PaiNN 型等变图神经网络的势函数，该模型是在由从头算分子动力学（AIMD）熔体淬火轨迹中提取的 11,184 个结构组成的训练集上训练而成的。训练集覆盖了四种 Na$_2$O 组分（从 Na$_6$Si$_{57}$O$_{117}$ 到 Na$_{48}$Si$_{36}$O$_{96}$）和四种质量密度，确保模型能够捕捉多样化的局部环境。损失函数高度加权了原子受力（$w_F=0.97$），以确保轨迹的高保真度。
• 模拟设置： 生产运行针对包含 3,000 个原子的系统，涉及三种不同的组分：Na$_2$O $\cdot$ 2SiO$_2$ (NS2)、Na$_2$O $\cdot$ 4SiO$_2$ (NS4) 和 Na$_2$O $\cdot$ 6SiO$_2$ (NS6)。
• 热力学方案： 系统在 4000 K 下平衡，并以 $10^{12}$ K s$^{-1}$ 的速率冷却至 300 K。生产运行在 NVT 系综下进行，温度范围从 4000 K 降至 1500 K。至关重要的是，模拟时长被延长到了纳秒量级（在 1500 K 时高达 5.1 ns），这一量级在之前是无法通过从头算方法实现的。
• 分析方法： 研究利用均方位移（MSD）、键断裂相关函数（$C_B(t)$）、范·霍夫分布函数的自部分（$G_s(r,t)$）以及非高斯参数（$\alpha_2(t)$）来表征弛豫动力学和空间非均匀性。

关键结果

1. 弛豫动力学与解耦： 模拟成功捕捉到了网络形成元素 Si 和 O 特有的、具有玻璃转变特征的两步弛豫过程，这由 MSD 和 $C_B(t)$ 在中间时间尺度上的明显平台期所证实。相比之下，Na 离子的弛豫显著更快。即使在 1500 K 时，Si 和 O 的相关性仍能持续约 1 ns，而 Na 的相关性在 10–100 ps 内就已衰减。这证实了一种动力学解耦现象，即碱金属的传输是由离散事件介导的，使得 Na 离子能够独立于迟缓的硅酸盐基质进行运动。
2. 跳跃机制： 钠的范·霍夫函数自部分在长时间尺度下表现出清晰的双峰特性。位于约 3.5 Å 处的第二峰对应于在相邻间隙位点之间的离散跳跃。这为碱金属传输通过跳跃机制而非连续扩散提供了直接证据。
3. 动态非均匀性： 非高斯参数 $\alpha_2(t)$ 表明，尽管所有物种都表现出动态非均匀性，但氧原子的行为最为显著。氧的 $\alpha_2(t)$ 峰值强度始终高于硅和钠。这归因于氧位点的结构多样性（桥联氧与非桥联氧）及其作为 SiO$_4$ 四面体顶点的角色，这些顶点会经历罕见的、随机的结构重排。
4. 组分依赖性： 增加钠含量（从 NS6 到 NS2）系统地加速了所有物种的动力学过程。然而，随着钠含量的降低（NS6），非高斯行为变得更加剧烈，这表明更刚性、更连通的网络增强了弛豫事件的间歇性以及移动粒子与不移动粒子之间的分离。

意义与主张
作者声称，这项工作通过弥合从头算精度与观察动态非均匀性所需的巨大时空尺度之间的差距，为研究现实玻璃系统建立了新的基础。

• 方法论进展： 该研究证明了 MLIPs 可以追踪在纳秒量级下的硅酸盐网络弛豫和扩散离子运动，这一量级此前在具有第一性原理精度的研究中是无法触及的。
• 物理洞察： 结果为“优先路径”图景提供了强有力的动力学证据，表明碱金属离子通过离散跳跃事件从基质中解耦。
• 非均匀性起源： 本文阐明了钠硅酸盐熔体中的动态非均匀性源于两种不同的机制：共价网络的罕见、随机重排（主导氧动力学）以及碱金属改性剂的快速、通道介导的跳跃。
• 未来应用价值： 通过实现接近从头算的精度并在结构弛豫区间内进行研究，该方法为理解局部约束如何决定宏观传输性质提供了严谨的基础，为离子传导材料的计算设计提供了新的标准。