本文介绍lammps统计一个原子周围不同类型原子数量的方法。
在之前的专栏中,曾介绍过动态统计某一个固定区域内原子数量的方法,也介绍过动态统计某一个原子周围原子数量的方法:
下面介绍第三种类型:动态统计某一个原子周围不同类型原子数量的方法。
以小球的随机碰撞为例,原子类型为1的红色小球(id也为1)周围分布原子类型为2、3的蓝色和黄色小球。
分别统计红色小球周围蓝色原子和红色原子的数量,需要用到count(type,region)命令,改命令可统计区域region内type类型原子的数量。
对应的代码为:
fix 1 all nve
group type2 type 2
group type3 type