lammps教程：不同原子组势能计算方法

大家好，我是小马老师。
本文如何使用lammps命令计算原子组的势能。
模拟吸附、拉拔、界面等过程时，需要计算不同原子组的势能，根据势能差计算吸附能或者界面能。
关于势能计算，lammps只提供了单原子势能计算命令（compute pe/atom)和一个默认的总势能（pe）。
对于不同原子组的势能，可以通过组内原子势能求和的方式得到。
如铜、水润湿体系，分别根据原子类型设置两个原子组：water、cu：

group water type 1 2
group cu type 3

计算每个原子的势能

compute mype all pe/atom

把两个组内的原子势能相加，得到原子的势能，compute reduce sum命令实现单原子量的求和计算。

compute