LAMMPS级联碰撞中出现杂乱原子及势函数力值与一阶导不一致？

zhaoyanx4

我计算了10w原子的玻璃（charge类型的非晶结构）模型用于级联碰撞模拟。但计算出的碰撞过程中在模拟盒边缘出现了大量离位原子（如下图一）。推测可能是势函数的问题，并且lammps运行过程中也对势值报了如下warning：WARNING: 1 of 5000 force values in table Si-OO are inconsistent with -dE/dr.
WARNING:  Should only be flagged at inflection points (../pair_table.cpp:465)
势参形式如下图二。我对势函数和势参反复确认、提高精度都没找到产生warning的原因，曾尝试定位报错的势参值，但这个值不在分段点周围并且重新计算也无法解决。想请教各位老师有没有曾遇见过类似情况，指点一下该如何解决级联碰撞或者势函数的问题。