7Feb2024 版本的LAMMPS新增一条关于反应力场的命令：compute reaxff/atom command。它可以生成一个per-atom类型的数组和一个local类型的数组。

compute reaxatom all reaxff/atom bonds yes

dump reax all local 23 reax.bond c_reaxatom[*]
dump_modify reax colname 1 id_i colname 2 id_j colname 3 bo

dump data some custom 23 atom.txt mol type id c_reaxatom[1] c_reaxatom[2] c_reaxatom[3] x y z vx vy vz 
dump_modify data colname c_reaxatom[1] abo colname c_reaxatom[2] nlp colname c_reaxatom[3] nb

per-atom类型的数组包含三列，分别是该原子的键序之和、孤对电子数和化学键数目。

local类型的数组也包含三列，分别是形成化学键的原子i、原子j和该化学键的键序。

如上所示的脚本输出为：

  1. local（由 dump local 输出）：


1

  1. per-atom（由 dump custom 输出）：


2

大佬论坛