[论文解读] Coarse-graining molecular dynamics models using an extended Galerkin projection

Xiantao Li|arXiv (Cornell University)|Oct 16, 2012

Advanced Mathematical Modeling in Engineering参考文献 44被引用 4

一句话总结

本文提出了一种扩展的伽辽金投影方法，用于晶体固体的分子动力学粗粒化，通过在缺陷附近增强粗粒度子空间，最小化原子/连续体界面处的声子反射。与传统伽辽金方法不同，该方法在不预先计算广义朗之万方程中记忆函数的情况下，保持了变分结构和机械平衡稳定性。

ABSTRACT

We present a new framework for coarse-graining molecular dynamics models for crystalline solids. The reduction method is based on a Galerkin projection to a subspace, whose dimension is much smaller than that of the full atomistic model. The subspace is expanded by adding more coarse-grain variables near the interface between lattice defects and the surrounding regions. This effectively minimizes reflection of phonons at the interface. In this approach, there is no need to pre-compute the memory function in the generalized Langevin equations, a typical model of interface conditions. Moreover, the variational formulation preserves the stability of mechanical equilibria.

研究动机与目标

• 解决粗粒化分子动力学中原子/连续体界面处的人工声子反射问题。
• 克服广义朗之万方程中预先计算记忆核的计算负担。
• 在粗粒化模型中保持原始哈密顿动力学的变分结构和机械平衡稳定性。
• 实现对晶体固体中非平衡动力学（如位错运动和冲击波）的精确模拟。
• 通过引入噪声，提供一种可扩展至有限温度系统的框架。

提出的方法

• 对由基函数张成的维数远小于完整原子系统规模的缩减子空间应用伽辽金投影。
• 通过在晶格缺陷和界面附近增加额外自由度来扩展粗粒度子空间，以抑制声子反射。
• 采用变分公式以保持机械平衡的稳定性，并确保与底层哈密顿动力学的一致性。
• 使用中点求积和有限元离散化近似粗粒度方程，得到具有柯西-伯恩本构关系的弹性动力学模型。
• 在界面处引入重复节点以增强子空间，有效建模过渡区域，而无需依赖广义朗之万方程。
• 将标准伽辽金投影替换为包含缺陷附近局部加密的扩展版本，以改善界面处的波传输性能。

实验结果

研究问题

• RQ1 如何在粗粒化分子动力学模型中最小化原子/连续体界面处的声子反射？
• RQ2 是否可以在保持精度的同时消除广义朗之万方程中预先计算记忆函数的需求？
• RQ3 如何在粗粒化公式中保持原始分子动力学模型的变分结构？
• RQ4 在缺陷附近增强子空间对粗粒化模拟的精度和稳定性有何影响？
• RQ5 该方法能否扩展至有限温度系统，同时满足涨落-耗散定理？

主要发现

• 与传统伽辽金方法相比，扩展伽辽金方法在位错偶极子模拟中显著减少了界面处的声子反射。
• 传统伽辽金方法由于界面反射导致原子区域出现较大数值误差。
• 在Krylov子空间中L=8时，扩展伽辽金方法的误差显著低于L=4，表明子空间增强可提升精度。
• 该方法避免了预先计算广义朗之万方程中的记忆核，简化了实现并降低了计算成本。
• 通过扩展伽辽金投影导出的粗粒化模型结果与全原子分子动力学一致，尤其在捕捉位错核心动力学和波传播方面表现优异。
• 变分公式确保了机械平衡的稳定性，而这一点在许多现有耦合方法中无法保证。