大家好,我是小马老师。
在绝大多数的lammps模拟过程中,两个原子之间只要产生了键连接(bond),这个bond默认“永久”存在。
在模拟过程中,这个键是不允许断裂的,当然,也不会有新的键生成。
当键连接的两个原子距离过大,超出势函数允许的范围,就会产生“bond missing“错误,此时,模拟就会中止。
既然键不能断裂也不能重新生成,也就不能模拟体系的化学反应过程。
但反应势是一个例外,这个势允许键的断裂和生成,比如可以模拟聚合物的交联或者裂解。
本文介绍反应势reax/c的设置方法。
- 设置原子类型为charge
atom_style charge
- 设置pair_style
pair_style reax/c lmp_control
其中,reax/c为反应势类型说明,lmp_control为输出控制文件,如果不需要对输出进行控制,可设置为NULL,如: