Non-numeric box dimensions - simulation unstable 错误原因及解决办法

大家好，我是小马老师。
本文继续介绍lammps常见的错误及解决办法：Non-numeric box dimensions - simulation unstable
最近几天，好几个培训班学员遇到这个问题，在交流群中也有朋友问同样的问题，如下图所示：

这个问题具有一定的代表性，也是初学者经常遇到的一个问题，lammps官网对这个问题给出了清晰的解释：

The box size has apparently blown up.
1

模拟盒子已经爆炸。
当模拟提示这个问题的时候，可以打开dump输出的轨迹文件，观察一下原子结构是否正常。
多数情况下，看到的轨迹应该如下图所示:

下面分析一下问题产生的原因及解决办法。
1.原因分析
既然盒子已经爆炸，说明原子获得了较高速度，短时间内移动到较远的位置。
根本原因是原子受力过大，已经远远超出正常的受力值。
原子受力过大原因有很多，比如初始模型有原子重叠、力场参数不准确、系综设置不合理等。
2.解决办法
出现这个错误后，首先应该判断模型是否合理，把模型导入到ovito中，查看原子是否重叠。
特别是聚合物或者界面体系，原子容易产生重叠，可使用minimize命令进行初步能量最小化。
如果结构没问题，需要验证力场参数是否准确，最简单的方法是多换几个不同力场进行测试，看能否解决问题。
如果不存在以上问题，需要检查系综设置是否合理，比如被固定的原子是否具有初速度、系综设置是否重复、刚体设置是否合理。
这些是比较常见的原因，但不是全部，当出现错误时，可先按以上的顺序进行排查，以上方法解决不了时，需要综合分析整个模拟过程，能否解决可能就要看个人经验了，毕竟lammps模拟有点复杂。
公众号：lammps加油站