大家好,我是小马老师。
本文分享一个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<