计算化学公社

标题: npt后改用P-R压浴进行md模拟，盒子剧烈振动，如何解决 [打印本页]

作者
Author: Jus 时间: 2023-11-10 16:41
标题: npt后改用P-R压浴进行md模拟，盒子剧烈振动，如何解决
本帖最后由 Jus 于 2023-11-11 22:04 编辑

使用gromacs进行充分的npt预平衡后，改用 parrinello-rahman压浴进行md模拟，体系盒子不停地剧烈缩放振动，不知如何解决，请求各位大佬的指点(已解决)
Tcoupl = V-rescale
tau_t = 2.0
tc_grps = system
ref_t = 298.15
;
Pcoupl = parrinello-rahman
pcoupltype = semiisotropic
tau_p = 4.0
ref_p = 0 1.01325
compressibility = 0 4.5e-5

作者
Author: wanzongliang 时间: 2023-11-10 17:18
有可能是之前的equilibrium没有真正平衡，或者你可以改变压浴试试

作者
Author: Jus 时间: 2023-11-10 19:26
本帖最后由 Jus 于 2023-11-11 22:04 编辑

wanzongliang 发表于 2023-11-10 17:18
有可能是之前的equilibrium没有真正平衡，或者你可以改变压浴试试

感谢你的回答。
已经跑了10ns，文本文件里说没平衡，后面续跑到13ns了，还是这样，所以我想可能不是没跑平衡的事

作者
Author: sobereva 时间: 2023-11-13 21:43
用c-rescale压浴

作者
Author: science_6 时间: 2025-3-25 08:58

Jus 发表于 2023-11-10 19:26
感谢你的回答。
已经跑了10ns，文本文件里说没平衡，后面续跑到13ns了，还是这样，所以我想可能不是没跑 ...

你好，请问一下怎么看预平衡(也就是eq)怎么查看体系有没有达到平衡？

作者
Author: student0618 时间: 2025-3-25 09:05

science_6 发表于 2025-3-25 08:58
你好，请问一下怎么看预平衡(也就是eq)怎么查看体系有没有达到平衡？

http://bbs.keinsci.com/thread-27122-1-1.html

作者
Author: science_6 时间: 2025-3-25 14:01

student0618 发表于 2025-3-25 09:05
http://bbs.keinsci.com/thread-27122-1-1.html

谢谢你，就是我想知道平衡相怎么样才能达到平衡

作者
Author: Lance先生 时间: 2025-3-25 14:43

science_6 发表于 2025-3-25 14:01
谢谢你，就是我想知道平衡相怎么样才能达到平衡

体系的总能，或者势能提取出来，画模拟时长与总能或者势能变化的图，能量波动上下平稳的话就是平衡了，社长专门写过一个帖子，你可以搜一下。

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