CaO-Al2O3-SiO2熔体中的自扩散系数：一种高通量分子动力学方法

《Materials & Design》：Self-diffusion coefficients in CaO-Al 2O 3-SiO 2 melts: a high-throughput molecular dynamics approach

　　

氧化物熔体中准确的扩散数据对于模拟冶金和地球科学中的高温过程至关重要，然而现有的测量仍然稀少且不一致。本研究提出了对CaO–Al₂O₃–SiO₂熔体中自扩散的高通量分子动力学（MD）研究。研究人员在液相区1400°C至1600°C之间，对40种成分进行了总计119次模拟，得到了Ca、Al、Si和O的扩散系数。研究人员分析了这些趋势与熔体结构的关系，并使用Thibodeau–Jung模型进行拟合，得到两种参数化：纯MD拟合和结合MD与选定实验数据的混合拟合。纯MD拟合表现出极好的内部一致性，并再现了物理上有意义的扩散系数排序。混合拟合改善了对实验数据的吻合度，特别是对于Al、Si和O。两种拟合均与Kwon的完全实验参数化进行了比较，并通过扩散偶实验和Al₂O₃及SiO₂颗粒溶解的相场模拟进行了验证。与Kwon的参数化不同，基于MD的参数化捕捉了结构趋势，并再现了实验观察到的上坡扩散。这项工作为氧化物熔体中的自扩散建立了一个经过验证的、可迁移的框架，并展示了MD-informed策略在数据稀缺的多组分氧化物系统中进行预测性输运建模的潜力。

**论文解读：CaO-Al₂O₃-SiO₂熔体自扩散系数的高通量分子动力学研究**

**研究背景与问题**
在冶金和地球科学领域，高温过程（如钢水精炼中的耐火氧化物溶解）的精确建模依赖于氧化物熔体中可靠的扩散数据。然而，CaO–Al₂O₃–SiO₂熔体的自扩散系数实验测量极为稀缺且不一致，现有模型（如Kwon等基于实验数据的参数化）在预测扩散系数排序时出现非物理结果（如D_Al > D_O > D_Ca > D_Si），与实验和分子动力学（MD）模拟结果矛盾。为填补这一数据空白并建立物理上一致的扩散描述，本研究采用高通量MD方法，系统获取温度-成分依赖的自扩散系数，并通过Thibodeau–Jung模型进行参数化，最终通过扩散偶实验和相场模拟进行验证。论文发表在《Materials》。

**主要技术方法**
研究人员采用以下关键技术：
1. **分子动力学模拟**：使用Bouhadja等开发的经验力场（Born–Mayer–Huggins势），通过LAMMPS软件（版本27 June 2024）在NPT系综下进行模拟，系统包含约5000个原子，生产运行时长50 ns，涵盖1400°C–1600°C（50°C间隔）和40种成分（液相区，共119个输入集）。
2. **自扩散系数计算**：基于Einstein关系，从均方位移（MSD）的线性扩散区提取Ca、Al、Si、O离子的自扩散系数，采用滑动窗口分析确定统计不确定性。
3. **Thibodeau–Jung模型拟合**：将扩散系数表示为温度与网络形成阳离子对分数（由Modified Quasi-Chemical Model（MQM）计算）的函数，通过非线性最小二乘拟合获得MD-only和混合（MD+实验）两种参数化。
4. **相场模拟**：采用三维相场模型模拟Al₂O₃和SiO₂颗粒在熔体中的溶解，验证扩散参数化的动力学预测能力（无样本队列来源，所有模拟数据基于MD和实验文献）。

**研究结果**
**3.1 MD-derived self-diffusion coefficients**
通过MD模拟获得的自扩散系数在所有成分和温度下遵循D_Ca > D_Al > D_O > D_Si的排序，与Ca作为网络修饰剂、Al和Si作为网络形成剂以及O与网络耦合的结构角色一致。氧扩散系数可通过Si和Al扩散的组成加权平均准确描述，而算术平均会高估D_O。扩散系数随CaO含量增加而增大（因桥氧转化为非桥氧），随温度降低而减小，且活化能随网络形成度（NBO/T）增加而升高。

**3.2 Fitting of the self-diffusion coefficients**
采用Thibodeau–Jung模型对MD数据进行拟合，得到MD-only和混合（MD+实验，权重因子w=1）参数化。MD-only拟合在内部一致性上表现优异，对所有元素的平均绝对百分比误差（MAPE）低于13.7%。混合拟合显著改善了Al、Si、O与实验数据的吻合度（MD-only对Al、Si、O的扩散系数低估，混合拟合校正了这一偏差），同时保留了正确的扩散系数排序。与Kwon的实验参数化相比，MD-based参数化避免了非物理的排序问题。

**3.3 Validation of diffusion profiles in CaO-Al₂O₃-SiO₂ melts**
通过扩散偶实验（Sugawara和Oishi的数据）验证，MD-based参数化（特别是混合拟合）成功再现了实验中的上坡扩散现象，而Kwon参数化因人为降低Ca迁移率以抑制假性上坡扩散，导致其他区域上坡扩散被削弱。混合拟合的均方误差（MSE）为0.394，优于MD-only的0.539，但略高于Kwon的0.269（后者针对这些实验调参）。

**3.4 Validation of Al₂O₃ and SiO₂ dissolution kinetics**
通过相场模拟对比Al₂O₃（Liu等实验）和SiO₂（Feichtinger等实验）颗粒溶解实验，所有参数化均预测了正确的溶解时间尺度。MD-based参数化（特别是混合拟合）在考虑对流缺失和有限域饱和效应后，给出的溶解时间合理；Kwon参数化在1500°C时预测溶解路径进入两相区（slag–anorthite），但实验未观察到二次相形成。MD-only拟合因低估Al、Si扩散而预测更长的溶解时间，混合拟合在实验数据丰富区域接近Kwon参数化，在数据稀疏区域更接近MD-only。

**讨论与结论**
讨论部分指出，Bouhadja力场对Al、Si、O扩散系数的系统性低估可通过混合拟合部分校正，同时揭示了未来力场改进（如机器学习势）的必要性。在实验数据稀缺的多组分氧化物体系中，MD-based策略提供了物理上合理的可扩展路径。研究结论翻译如下：
“本研究证明了分子动力学（MD）作为一种强大且物理上合理的方法，可用于确定CaO–Al₂O₃–SiO₂熔体中温度-成分依赖的扩散系数。使用Bouhadja经验力场，在液相区对40种成分和5个温度（1400°C–1600°C）进行了119次MD模拟，获得了Ca、Al、Si和O自扩散系数的密集且一致的数据集。预测的扩散趋势反映了各元素在熔体网络中的已知结构角色，并在所有条件下一致再现了正确的扩散系数排序D_Ca > D_Al > D_O > D_Si，与先前的实验和计算研究一致。为将MD结果集成到介观模拟中，采用Thibodeau–Jung模型（该模型考虑了熔体中网络形成的影响）对自扩散系数进行了拟合，开发了MD-only和混合（等权重结合MD与实验数据）两种参数化。MD-only模型以高内部精度忠实捕捉了成分和温度趋势，而混合模型改善了对现有实验扩散系数的拟合，特别是对于Al、Si和O。通过实验基准的验证强化了MD-based方法的稳健性：1773 K下的扩散偶模拟再现了上坡扩散等关键特征，而Al₂O₃和SiO₂颗粒溶解的相场模拟预测了正确的溶解时间尺度。这些结果巩固了MD-based方法在模拟氧化物熔体中扩散控制转变的可靠性和可迁移性。本研究建立了一个物理一致、高保真的计算框架，用于模拟CaO–Al₂O₃–SiO₂熔体的动力学。由于CaO–Al₂O₃–SiO₂体系在实验上相对表征充分，得以与实验进行详细比较并开发混合参数化方法。所得结果桥接了原子尺度模拟与介观输运模型，提供了一个经过验证、可迁移的框架，用于模拟氧化物熔体中的扩散。同时，Bouhadja力场对熔体动力学的轻微低估凸显了未来使用机器学习势进行改进的机会。然而，在实验数据稀缺或不可用的更复杂多组分氧化物体系中，本研究提出的MD-based策略为冶金和地球科学应用中的预测性扩散建模提供了一条物理上合理且可扩展的道路。”