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