做美食网站视频,九江网站建设哪家好,最近中国新闻事件,高端网站建站公司大家好#xff0c;我是小马老师。
本文分享一个reaxff反应势的案例。
该案例主要模拟Fe(OH)3在高温下的反应过程#xff0c;主要代码来自lammps自带的案例。
lammps自带案例没有产物输出#xff0c;故在此基础上稍加修改#xff0c;增加了产物输出命令。
反应过程如下图…大家好我是小马老师。
本文分享一个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替换。 如果需要系统学习反应势的模拟可考虑报名培训课程。
