lammps统计一个原子周围不同类型原子数量的方法

本文介绍lammps统计一个原子周围不同类型原子数量的方法。
在之前的专栏中，曾介绍过动态统计某一个固定区域内原子数量的方法，也介绍过动态统计某一个原子周围原子数量的方法：
下面介绍第三种类型：动态统计某一个原子周围不同类型原子数量的方法。
以小球的随机碰撞为例，原子类型为1的红色小球（id也为1）周围分布原子类型为2、3的蓝色和黄色小球。

分别统计红色小球周围蓝色原子和红色原子的数量，需要用到count(type,region)命令，改命令可统计区域region内type类型原子的数量。
对应的代码为：

fix             1 all nve
group type2 type 2
group type3 type