大家好，我是小马老师。

本文分享一个reaxff反应势的案例。

该案例主要模拟Fe(OH)3在高温下的反应过程，主要代码来自lammps自带的案例。

lammps自带案例没有产物输出，故在此基础上稍加修改，增加了产物输出命令。

反应过程如下图所示：

in文件代码已经注释，以供参考。

#模型基本设置
units           real
atom_style      full
read_data       FeOH.data
neighbor        2 bin
neigh_modify    every 10 delay 0 check no
timestep        0.25
thermo 100
#势参数设置
pair_style      reaxff lmp_control
pair_coeff      * * ffield.reax.Fe_O_C_H H O Fe
#系综控温
fix             1 all nve
fix             2 all qeq/reax 1 0.0 10.0 1e-6 param.qeq
fix             3 all temp/berendsen 1000.0 1000.0 100.0
#产物命令
fix     4 all reaxff/species 1 100 100 species.out element  H O Fe
fix     5 all reaxff/bonds 100 bonds.reaxff
#输出轨迹文件
dump 1 all custom 100 dump.xyz id mol type q x y z
run             5000

所有的reaxff模拟，如：氧化、腐蚀、热裂解、燃烧等模拟，均可采用此模拟框架，只需要更换data文件、势文件，代码中的control文件可以不写，电荷平衡文件也可以用reaxff替换。
如果需要系统学习反应势的模拟，可考虑报名培训课程。