[论文解读] Implementing Reactivity in Molecular Dynamics Simulations with Harmonic Force Fields

Jordan J. Winetrout, Krishan Kanhaiya+8 more|arXiv (Cornell University)|Jul 30, 2021

Polymer crystallization and properties被引用 7

一句话总结

本文提出了一种反应性 INTERFACE 力场（IFF-R），通过用能量守恒的Morse势能替代经典谐波键势能，实现反应性分子动力学模拟。该方法为每种键类型定义三个可解释的Morse参数，并在解离时设为零能量，从而在保持精度的同时，实现对聚合物、蛋白质和金属等多种材料的键断裂与形成模拟，相较于以往的反应性方法，性能提升约30倍。

ABSTRACT

The simulation of chemical reactions and mechanical properties including failure from atoms to the micrometer scale remains a longstanding challenge in chemistry and materials science. Bottlenecks include computational feasibility, reliability, and cost. We introduce a method for reactive molecular dynamics simulations using a clean replacement of non-reactive classical harmonic bond potentials with reactive, energy-conserving Morse potentials, called the Reactive INTERFACE Force Field (IFF-R). IFF-R is compatible with force fields for organic and inorganic compounds such as IFF, CHARMM, PCFF, OPLS-AA, and AMBER. Bond dissociation is enabled by three interpretable Morse parameters per bond type and zero energy upon disconnect. Use cases for bond breaking in molecules, failure of polymers, carbon nanostructures, proteins, composite materials, and metals are shown. The simulation of bond forming reactions was included via template-based methods. IFF-R maintains the accuracy of the corresponding non-reactive force fields and is about 30 times faster than prior reactive simulation methods.

研究动机与目标

• 解决在原子至微米尺度上模拟化学反应与机械失效的长期挑战。
• 克服反应性分子动力学中的计算瓶颈，特别是高昂的成本与有限的可靠性。
• 开发一种方法，实现键的离解与重新形成，同时保持现有非反应性力场的精度。
• 确保与广泛使用的力场（如 IFF、CHARMM、OPLS-AA、AMBER 和 PCFF）兼容。
• 在不牺牲反应性模拟物理保真度的前提下，实现显著的性能提升。

提出的方法

• 用能量守恒的Morse势能替代经典谐波键势能，以实现键的断裂与形成。
• 为每种键类型定义三个可解释的Morse参数，分别控制键能、平衡键长和解离能。
• 在无限分离时设势能为零（断开时能量为零），以确保能量守恒。
• 与现有力场（如 IFF、CHARMM、OPLS-AA、AMBER、PCFF）集成，实现对有机与无机材料的无缝应用。
• 采用基于模板的方法模拟成键反应，实现动态反应路径的模拟。
• 使用标准分子动力学算法，结合修改后的势能函数，以保持计算效率。

实验结果

研究问题

• RQ1 能否构建一种反应性力场，在保持非反应性谐波力场精度的同时，实现键的断裂？
• RQ2 如何利用Morse势能以物理解释一致且计算高效的方式建模键离解？
• RQ3 该方法在多大程度上能够模拟聚合物、碳纳米结构和复合材料中的复杂失效机制？
• RQ4 该方法能否在包括蛋白质和金属在内的多种化学体系中实现泛化，且仅需极少的重新参数化？
• RQ5 与现有反应性分子动力学方法相比，该方法能实现多大的性能提升？

主要发现

• IFF-R通过用能量守恒的Morse势能干净地替代谐波势能，实现了在分子动力学模拟中键的断裂与形成。
• 该方法相比以往的反应性模拟方法，性能提升约30倍，显著提高了计算可行性。
• 在无限分离时势能为零，准确建模了键离解，确保了能量守恒。
• 该力场在多种体系（包括聚合物、蛋白质和金属）中保持了原始非反应性力场的精度。
• 应用案例成功模拟了聚合物和碳纳米结构中的机械失效，以及通过基于模板的方法实现的动态成键过程。
• IFF-R与主流力场（如 IFF、CHARMM、OPLS-AA、AMBER 和 PCFF）兼容，具备广泛的应用潜力。