2023.01.15 rgd32 会议文章 accepted

今天rgd32的会议文章 accepted with minor revisions 啦！

ReaxFF自带的输出命令

根据 pair_style reaxff command 的描述，ReaxFF自带14个能量有关的输出命令，在in文件中插入：

compute reax all pair reaxff
variable eb      equal c_reax[1]
variable ea      equal c_reax[2]
variable elp      equal c_reax[3]
variable ew    equal c_reax[11]

variable emol      equal c_reax[4]
variable ev      equal c_reax[5]
variable epen      equal c_reax[6]
variable ecoa      equal c_reax[7]
variable ehb      equal c_reax[8]
variable et      equal c_reax[9]
variable eco      equal c_reax[10]
variable ep      equal c_reax[12]
variable efi      equal c_reax[13]
variable eqeq      equal c_reax[14]

thermo_style custom step temp epair v_eb v_ea v_elp v_ew v_emol v_ev v_epen v_ecoa v_ehb v_et v_eco v_ep v_efi v_eqeq

其中，epair是lammps自己输出的，后面是reax输出的。

对于1个O2，结果上显示为 epair = v_eb+v_ea+v_elp+v_ew，其他量是0。

对于1个O，结果为 epair = v_ea+v_elp，其他为0。

对于整个C(0001)，结果为 epair = v_eb+v_ea+v_elp+v_ew+v_ev+v_epen+v_et+v_eco+v_ep+v_eqeq，其他为0。

根据从输出得到的规律，用ReaxFF计算时，lammps自己认为的epair=ep=reaxff全部输出之和。epair的大小正好是每个原子compute pe/atom 命令输出的pe之和。

另外，按手册的说法，epair = evdwl + ecoul。