计算化学公社

标题: gromacs钠电池电解液NPT平衡阶段压力波动极大 [打印本页]

作者
Author: Tonue 时间: 2025-11-13 22:22
标题: gromacs钠电池电解液NPT平衡阶段压力波动极大
各位老师好，我在模拟电解液体系时，在NPT阶段压力波动极大（RMSD达到数千）且运行时一直伴随如下提示：Warning: pressure scaling more than 1%。体系总原子数九千左右，分子数在三百左右。已经试过了C-rescale，PR，bredesen三种方式，PR会在运行一段时间后直接崩溃，另外两种可以产生结果。希望各位老师能够指点一下 (, 下载次数 Times of downloads: 3)

作者
Author: student0618 时间: 2025-11-13 22:33
pressure看的是平均，gmx energy出的波动大正常。

要看收敛看Density或volume或box x/y/z

以下帖文已经有个Pressure的图及解说：
谈谈怎么判断分子动力学模拟是否达到了平衡
http://bbs.keinsci.com/forum.php ... 2&fromuid=64740
(出处: 计算化学公社)

作者
Author: Tonue 时间: 2025-11-14 00:03

student0618 发表于 2025-11-13 22:33
pressure看的是平均，gmx energy出的波动大正常。

要看收敛看Density或volume或box x/y/z

好的，感谢您的解答和版主的文章，受益良多

作者
Author: sobereva 时间: 2025-11-15 06:18
我在论坛里说过无数次，体系达到平衡之前不要用PR压浴，这压浴不是什么好东西

那叫Berendsen不叫bredesen

作者
Author: yuzc 时间: 2025-12-24 15:05
PR压浴也是扩展相空间法，对体系的初态要求还挺严格的，应该充分平衡体系。
所以一般就是先用berendsen平衡完再用PR采样。不过现在可以直接从头到尾用C-rescale了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)