大家好,我是小马老师。
本文继续介绍lammps常见的错误及解决办法:Non-numeric box dimensions - simulation unstable
最近几天,好几个培训班学员遇到这个问题,在交流群中也有朋友问同样的问题,如下图所示:
这个问题具有一定的代表性,也是初学者经常遇到的一个问题,lammps官网对这个问题给出了清晰的解释:
The box size has apparently blown up.
模拟盒子已经爆炸。
当模拟提示这个问题的时候,可以打开dump输出的轨迹文件,观察一下原子结构是否正常。
多数情况下,看到的轨迹应该如下图所示:
下面分析一下问题产生的原因及解决办法。
1.原因分析
既然盒子已经爆炸,说明原子获得了较高速度,短时间内移动到较远的位置。
根本原因是原子受力过大,已经远远超出正常的受力值。
原子受力过大原因有很多,比如初始模型有原子重叠、力场参数不准确、系综设置不合理等。
2.解决办法
出现这个错误后,首先应该判断模型是否合理,把模型导入到ovito中,查看原子是否重叠。
特别是聚合物或者界面体系,原子容易产生重叠,可使用minimize命令进行初步能量最小化。
如果结构没问题,需要验证力场参数是否准确,最简单的方法是多换几个不同力场进行测试,看能否解决问题。
如果不存在以上问题,需要检查系综设置是否合理,比如被固定的原子是否具有初速度、系综设置是否重复、刚体设置是否合理。
这些是比较常见的原因,但不是全部,当出现错误时,可先按以上的顺序进行排查,以上方法解决不了时,需要综合分析整个模拟过程,能否解决可能就要看个人经验了,毕竟lammps模拟有点复杂。
公众号:lammps加油站