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 。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • 转载自https://zhuanlan.zhihu.com/p/257811414 LAMMPS之约束方法 在分子...
    wuhui1993阅读 8,556评论 0 1
  • 我的图书馆 留言交流 2000个最常用的英语单词 2016-05-15酴羰骀璺 2000个最常用的英语单词 (英语...
    Mr_Wang92阅读 1,074评论 0 0
  • 1 Lammps输入脚本结构 1.1 初始化命令 设置在创建原子或从文件中读取原子之前需要定义的参数。相关命令:u...
    Yuanshuo阅读 27,146评论 1 19
  • Chapter 1 :分类(Classificatio) 1.1 高斯混合模型(GMM:Gaussian-Mixt...
    jdzhangxin阅读 9,266评论 0 5
  • abandon, desert, forsake, leave, give up abandon :强调永远或完全...
    sunxiaohang阅读 3,919评论 0 3

友情链接更多精彩内容