lammps案例：反应势reax模拟氢气燃烧过程

大家好，我是小马老师。

本文根据Lammps官方案例写了一个氢气燃烧的lammps in文件，使用fix reax/c/species对氢气燃烧后的反应物进行分析。

主要模拟过程：

（1）建模

使用第三方软件建立H2和O2混合模型，其中氢气分子共200个，氧气分子100个。

（2）读取data文件，并设置反应势参数

pair_style  reax/c lmp_control
pair_coeff  * * ffield.reax.cho H O
fix         q all qeq/reax 1 0.0 10.0 1e-6 param.qeq

（3）设置反应物输出

fix 2 all reax/c/species 1 100 100 species.out element H O

每100步输出一次反应物的种类和各反应物的数量，数据保存到species.out文件中，可编程统计各反应物随时间的变化情况。

以下为in文件代码

units   real
atom_style  charge
read_data H2_O2.data

pair_style  reax/c lmp_control
pair_coeff  * * ffield.reax.cho H O
fix         q all qeq/reax 1 0.0 10.0 1e-6 param.qeq
neighbor  2 bin
neigh_modify  every 10 delay 0 check no
timestep  0.25

dump        1 all atom 100 mini.xyz
minimize    1e-8 1e-8 10000 10000
undump      1
reset_timestep  0

velocity    all create 1000 9898
fix         2 all reax/c/species 1 100 100 species.out element H O 
dump        1 all atom 1000 nvt.xyz
fix       1 all nvt temp 1000 1000 25
run         10000

燃烧效果如下图所示：

更多lammps教程，请关注微信公众号：lammps加油站