ReaxFF 势，全称 reactive force field，即反应力场，由 Adri van Duin 和 William A. Goddard,

III 于 2001 年提出[22]。ReaxFF 势包含以下 9 项：

Esystem  Ebond Eover Eunder Eval  Epen Etors  Econj EvdWaals  ECoulomb

（2.2）

其中 Ebond 为键能， Eover 为过饱和键能， Eunder 为不饱和键能， Eval 为键角能， Epen 为共价 键修正能， Etors 为二面角能， Econj 为共轭能， EvdWaals 为范德华相互作用力能， ECoulomb 为静电

相互作用能。ReaxFF 是一种依靠键级的力场，它考虑了键距与键级之间、键级与键能之间的 关系以及一些化学价理论，能够很好的描述化学键的断裂与生成。另外，ReaxFF 还包含了范 德华相互作用以及静电相互作用来描述非键相互作用。来!自~751论-文|网www.751com.cn

由于 ReaxFF 势精确考虑了分子的复杂性（涉及到键长、键角、二面角等等），因此一 般情况下使用 ReaxFF 势计算比 AIREBO 势准确。本文第一部分研究高能脉冲下的缺陷转化 时采用了两种势函数分别进行了计算，在第二部分概率统计主要是在 ReaxFF 势下进行，第 三部分加载应力模拟使用 AIREBO 势。文献综述

2.3.3 边界条件

Lammps 中可设置的边界条件有四种：p,f,s,m。

p 即为周期边界条件。除了 p 之外另外三种均为非周期边界调节。其中 f 为 fixed，即固 定边界，如果一个原子跨过了边界，则这个原子视为“丢失”。s 为 shrink-wrapped，即无论 原子运动到哪，边界都要调整到包裹所有原子，就像一层塑料膜包裹一样。m 为 shrink-wrapped with a minimum value，在 s 边界条件的基础上有一个最小值，比如 z 方向的上边界被定为 50 或以上，则即使所有的原子都在 50 以内，z 方向上边界始终大于或等于 50。

本文第一、第二部分工作边界条件设置为 m，第三部分工作设置为周期性边界 p。