lammps案例：反应力场reax/c设置方法及案例解析

大家好，我是小马老师。
在绝大多数的lammps模拟过程中，两个原子之间只要产生了键连接（bond），这个bond默认“永久”存在。
在模拟过程中，这个键是不允许断裂的，当然，也不会有新的键生成。
当键连接的两个原子距离过大，超出势函数允许的范围，就会产生“bond missing“错误，此时，模拟就会中止。
既然键不能断裂也不能重新生成，也就不能模拟体系的化学反应过程。
但反应势是一个例外，这个势允许键的断裂和生成，比如可以模拟聚合物的交联或者裂解。
本文介绍反应势reax/c的设置方法。

1. 设置原子类型为charge

atom_style	charge

2. 设置pair_style

pair_style	reax/c lmp_control

其中，reax/c为反应势类型说明，lmp_control为输出控制文件，如果不需要对输出进行控制，可设置为NULL，如：