lammps原子组速度的计算方法

大家好，我是小马老师。
本文介绍lammps模拟过程中原子组速度的计算方法。

在lammps摩擦、冲击、压痕等模拟过程中，涉及到某个原子组的整体移动，如果要输出这个原子组的整体移动速度，可以使用vcm()函数。

vcm()是lammps自带的一个函数，返回原子组质心在设定方向的速度。

对应的命令格式为
vcm(ID,dim)
其中，ID为group-ID，dim为方向，可设置x、y或者z方向。

下面是一个具体的使用案例，金刚石球在盒子内移动，输出移动过程中金刚石球在x方向的质心速度。

对应的代码为：

#初始模拟系统设置
units         metal
boundary      p p p
atom_style    atomic
neighbor	  2 bin
timestep      0.001
#构建模型
lattice       diamond