2022-10-04 Lammps 关于水模型设置的求助

问题描述如下：

模拟选用SPC/E模型进行处理，但是加上用SPC/E对应的bond coeff 设置为0 angle coeff 设置为 0 以后 添加fix shake 命令 之后run 的时候发现运行会被卡住，并不会进行计算，这种现象很怪异

TIP3P长程力作用和SPC/E bond coeff 设置为0以后计算不出结果

同样是使用fix shake 命令的 TIP3P模型短程带能量的结构就可以运行出结果，不知道为什么会出现这种问题，希望有懂得大神帮忙解答一下。

TIP3P短程力模型下可以计算出结果

---

问了问大佬，说是设置fix shake 命令以后就不用管K的值到底是多少了，但是要保证下一秒计算迭代收敛，还是要分配一下K的能量的，不然运算会被一直卡住，之前fix shake 命令被卡住的原因在于计算不收敛，这种不收敛是由多种方面引起的，可能是fix shake 命令不收敛，也有可能是自己的bond atom missing 导致，但是施加fix shake 命令以后会被看不出来，只会显示如上的fix shake 运行界面，建议先去掉fix shake 命令以后把会被盖住的所有命令解决掉，之后再fix shake 。