lammps动态统计中心原子周围原子数量

本文介绍lammps统计一个原子周围一定范围内其它原子数量的方法。

在之前的专栏中，曾介绍过动态统计某一个固定区域内原子数量的方法。

这种方法只能统计一个静态的区域内原子的数量变化，如果要统计某一个运动中的原子周围原子数量的变化，就需要根据该原子的运动状态，实时追踪原子的坐标并更新region区域的范围，然后再统计周围原子的数量。

以小球的随机碰撞为例，原子类型为1的红色小球（id也为1）周围分布原子类型为2、3的蓝色和黄色小球。

动态统计以红色小球为圆心，半径为15的球范围内的数量代码为：

#气体粒子数量
variable npart  equal 200
#模拟基本参数设置
units    lj
dimension 2
atom_styl