lammps运行原子散开解决方法

——作者水平有限，仅供学习交流使用，如有问题还请指出

Lammps在运行某个in文件时发现某部分原子在初始建模是正常的，但在之后的时间步就会散开（还有一种在初始建模时原子就已经散开），本人通过平时与ERROR斗智斗勇过程并加以思考总结出该问题出现的可能原因，照例，我们根据一个实例去说这个问题，in文件如下；

每隔1000步输出在dump文件中，总共运行5000步也就是说共保存了5次。我们可以看到在初始（0）时，模型是我们所需要的，但在此之后金刚石球部分的原子散开。

造成这个现象的原因是金刚石球有了初始速度，可以回看in文件区域划分以及温度（速度）初始化的group。

注意看10-16行命令，可能会觉得我们对mobile层温度初始化，怎么金刚石（ball）就有速度了，仔细看12行命令（mobile=all-lower）是不是把金刚石（ball）也化为mobile区了，是不是对金刚石也温度初始化并赋予其速度了？

好，既然找到问题所在，就很好去改正，命令修改为

修改之后运行尝试一下

运行正常，修改完成。