小木虫

当前位置：首页 > 分子模拟 >请教势函数的书写哪出错了？

请教势函数的书写哪出错了？

作者 huanghehua: 来源: 小木虫 150 3 举报帖子

+关注

之前有一个写法正确但准确性并不好的势函数，其余不动只替换成下图中的势函数之后报错，应该是势函数的写法上出了问题，希望可以得到老师们的指导

in.png

报错.png

势参数.png

势函数.png 返回小木虫查看更多

今日热帖

精华评论

华强

Lost atom, 原子丢失了。一开始有376个，在第一步的时候，LAMMPS发现只有127个了，然后就报错告诉你原子丢失了。这个错误通常是对原子做某种操作造成的原子丢失，这部分在你设置势函数之后。就你截图的来看，你应该是对原子做了某个热学处理，因为报错是由thermo这个cpp文件反馈回来的。你检查一下是不是对原子的操作那部分出了问题。
huanghehua

引用回帖:
2楼: Originally posted by 华强 at 2019-09-09 20:36:34
Lost atom, 原子丢失了。一开始有376个，在第一步的时候，LAMMPS发现只有127个了，然后就报错告诉你原子丢失了。这个错误通常是对原子做某种操作造成的原子丢失，这部分在你设置势函数之后。就你截图的来看，你应该 ...

感谢，检查了thermo之后，不出现原子丢失，方便再请教一下，在上图中也可以看到，能量最小化过程中的etotal数量级极大且无法收敛，之后在fix 1 all npt 1.0 1.0 0.1 iso 0.0 0.0 1.0 drag 0.2过程中dump原子相对位置没有发生改变，但盒子的体积在不断缩小，导致最终报错近邻表太大，我该往哪个方向修改
，
oxox6085

引用回帖:
3楼: Originally posted by huanghehua at 2019-09-10 07:37:34
感谢，检查了thermo之后，不出现原子丢失，方便再请教一下，在上图中也可以看到，能量最小化过程中的etotal数量级极大且无法收敛，之后在fix 1 all npt 1.0 1.0 0.1 iso 0.0 0.0 1.0 drag 0.2过程中dump原子相对位 ...

多贴点input可能会知道，感觉是没加上电荷