求助：lammps中高熵合金在进行短程有序优化过程中特定原子向边界上挤压该如解决？

墨茶茶

各位老师好！我最近在用lammps做5元高熵合金的分子动力学短程有序优化，当我优化到一定步数的时候特定原子发生了向边界挤压的情况，疑似在短程有序优化过程中盒子的边界条件发生了改变？短程有序优化过程中我使用的是NPT系综，三个方向的边界条件都是周期性边界。请问各位老师出现这种情况是什么原因呢？该如何解决呢？以下附上我的模型和短程有序优化代码，恳请各位老师帮助！

原子排布

相分布

短程有序部分脚本

墨茶茶

自我纠错，经检查发现在建模过程中region命令使用的是 units box为单位，导致模型从一开始就不具有周期性(在边界上原子重叠)，改为units lattice之后就没有再发生相同问题

yhz301415

同学你好，我也是在做短程有序，可以跟你求一个in文件吗。我们做的体系不一样，我是难熔中熵合金，我刚接触短程有序，我想请教一个问题，就是我看到文献上都是给MC一个温度进行迭代，从而根据aij认为产生了短程有序，这适用于所有高熵合金吗，我也不知道我前面对于产生短程有序的理解是否正确。