大家好,我是小马老师。
本文介绍lammps计算原子组之间相互作用的一个命令:compute group/group。
在lammps模拟中,经常需要计算原子组之间的相互作用,如摩擦球与基体之间、流体与壁面之间等等。
lammps提供compute group/group命令可完成此功能。
compute group/group命令既可以计算原子组之间的作用能,也可以计算两个组之间的作用力。
如:
compute 1 lower group/group upper
lower和upper为两个原子组的id(group-id)。
以上命令共有4个计算结果,其中标量c_1为量原子组之间的作用能,矢量c_1[1]、c_1[2]、c_1[3]分别为lower和upper组在x、y、z方向的作用力。
下面以Cu弛豫为例,简单演示compute group/group使用方法及结果输出。
模型共分上下两部分,计算结果由thermo命令输出。