《Computer Physics Communications》：BORAT-EM: High Frequency Electromagnetic Solver for Atmospheric Entry Blackout Analysis
编辑推荐：
为解决大规模原子系统分子动力学（MD）模拟的计算瓶颈，同时保持关键区域（如接触面、裂纹区）的原子级精度，研究人员提出了一种新型混合多尺度方法。该方法结合了全原子模拟和由优化势参数控制的“准原子”连续介质建模，并利用机器学习框架进行高效参数优化。通过与LAMMPS软件集成，该方法在保持高精度的前提下，显著提升了计算速度，为固体力学研究开辟了新途径。
分子动力学（MD）模拟已成为研究原子尺度材料的不可或缺工具，然而，尽管计算能力呈指数级增长，大规模系统的全原子模拟在计算上仍然难以承受。许多物理过程，如固体中的裂纹扩展、纳米粒子聚集、烧结、自组装、枝晶生长或相分离动力学，都需要同时进行原子和连续介质描述，以捕捉局部化学性质和长程形态。传统的纯连续介质或纯粗粒化方法在此类情况下往往失效。因此，发展一种能够无缝衔接原子与连续尺度、且计算高效的混合方法，是材料模拟领域的一大挑战。
针对这一难题，Skolkovo科学技术研究所的Artem Chuprov、Egor Nuzhin、Alexey Tsukanov和Nikolay Brilliantov团队在《Computer Physics Communications》上发表研究，提出了一种新颖、快速且精确的混合MD–介观尺度方法。该方法的核心是准原子（quasi-atoms） 概念，它将连续介质区域视为由不同尺寸的“准原子”构成的复合材料。研究人员通过一种概念上等同于在线机器学习（Machine Learning， ML）的优化方法，自动校准准原子间的相互作用势，使其构成的复合介质的弹性特性与真实原子介质相匹配。这种方法允许直接应用标准分子动力学软件包（如LAMMPS），并辅以基于机器学习的优化器，从而实现了对大规模系统的高效模拟，同时保持了关键区域的原子精度。
为了开展这项研究，作者主要采用了以下几项关键技术方法：首先，提出了“准原子”作为连续介质区域动态离散节点的概念，其尺寸可以从原子尺度平滑过渡到宏观尺度。其次，构建了连接Python优化器与LAMMPS MD引擎的桥梁框架。第三，采用了一种基于Regis和Shoemaker建立的代理优化（surrogate optimization） 方法（一种高效的“高成本黑箱优化”技术），通过定义并最小化一个损失函数（loss function） 来优化准原子势参数，该函数量化了复合介质与真实介质弹性张量分量之间的相对误差。最后，使用优化后的参数，在LAMMPS中对不同尺寸和组成的粒子碰撞进行了全原子与混合模拟的对比验证，并与经典接触理论（赫兹理论和JKR理论）进行比较。
研究结果主要分为以下几个部分：
1. 有效势参数的优化结果
优化程序针对包含N=2、3和4种原子/准原子类型的多尺度介质实施。结果表明，经过30-40次迭代后，目标参数（弹性常数）的相对误差可收敛至低于0.94%，同时计算时间大幅减少，验证了该方法的数值效率。
2. 粒子碰撞：原子尺度与多尺度建模对比
为验证多尺度方法，研究人员对球形粒子碰撞进行了复合结构与全原子结构的对比模拟。结果显示，无论是纯排斥性碰撞还是具有吸引力的碰撞，两种模拟得到的碰撞时间演化和相互作用力都表现出显著的一致性，证实了该方法在保持碰撞过程基本物理规律的同时，大幅降低了计算复杂度。
3. 与宏观接触理论的比较
进一步，研究模拟了微米级（半径0.1微米，约含2亿原子）大粒子的碰撞，并与经典接触理论——非粘性接触的赫兹（Hertz）理论和粘性接触的约翰逊-肯达尔-罗伯特（JKR）理论——进行了比较。对于纯排斥碰撞，模拟结果与赫兹理论在正重叠区域（即机械接触时）吻合良好。对于具有吸引力的碰撞，在正重叠区域也与JKR理论基本一致，但在负重叠区域（无机械接触时）则观察到显著偏差，这归因于JKR理论未考虑非接触的范德华吸引力。研究还使用多体Tersoff势对硅粒子进行了模拟，得出了类似结论。
4. 多尺度方法的计算效率
所开发的多尺度方法显著减少了大型系统的计算时间。模拟结果表明，对于半径为30埃的粒子，复合粒子模拟比全原子模拟快十倍；对于40埃的粒子，加速效果更为显著。此外，准原子类型越多，计算时间减少越明显。对于半径35埃的粒子，使用四种原子/准原子类型的模拟比使用两种类型的快约四倍。
5. 与现有其他方法的比较
研究将本方法与准连续（QC）方法、原子到连续介质（AtC）方法、粗粒化/自适应分辨率方案（AdResS）、异质多尺度方法（HMM）以及近场动力学等现有多尺度方法进行了对比。分析指出，本方法在保持原子模拟能力、适用于固体、开源实现以及计算速度方面具有综合优势。特别是与LAMMPS中的AtC包相比，在相同的模拟设置下，本方法即使使用更高分辨率（更多准原子），其计算速度仍然更快。
研究结论与讨论部分强调，该工作成功开发了一种新颖的混合多尺度模拟方法，能够精确描述原子尺度相互作用和连续介质动力学。其核心创新在于引入了“准原子”概念，并利用基于机器学习的代理优化技术高效确定其相互作用势参数，确保了复合介质与真实原子介质在弹性特性上的一致性。通过与LAMMPS集成，该方法实现了对标准分子动力学软件的便捷应用。
该方法的意义重大。首先，它有效克服了全原子模拟的计算限制，使得模拟包含数亿原子的微米级系统成为可能，为固体力学中的多尺度问题（如裂纹增长、压痕等）研究开辟了新途径。其次，通过与经典连续介质理论的对比验证，不仅确认了该方法在机械接触区域的准确性，还揭示了在存在吸引力时，连续理论（如JKR）需补充非接触原子力才能完整描述微观相互作用，这对理解微纳尺度接触力学具有指导意义。最后，与现有其他混合方法相比，该方法在计算速度、实现便利性和精度方面展现出显著优势，尤其适用于对计算效率和局部原子精度均有高要求的复杂固体系统模拟。该框架具有通用性，可扩展至其他势函数（如金属系统的Morse势、机器学习驱动的势函数）和更多物理现象（如粘弹性）的建模，展现了广阔的应用前景。
