计算化学公社

 找回密码 Forget password
 注册 Register
Views: 320|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 求助Parrinello-Rahman压浴和deform剪切模拟问题

[复制链接 Copy URL]

3

帖子

0

威望

39

eV
积分
42

Level 2 能力者

跑了一个蛋白质配体复合物体系,使用Parrinello-Rahman压浴和deform模拟剪切流,

在跑一段时间后会出现Pressure scaling more than 1%. This may mean your system is not yet equilibrated. Use of Parrinello-Rahman pressure coupling during equilibration can lead to simulation instability, and is discouraged.的提示,没有管,继续跑。
然后提示Step correcting invalid box:
   old box (3×3)
   old box[   0]={2.47868e+01, 0.00000e+00, 0.000000e+00}
   old box[   0]={2.47868e+01, 0.00000e+00, 0.000000e+00}
   old box[   0]={2.47868e+01, 0.00000e+00, 0.000000e+00}
最后显示one or more water can not be settled,最终崩溃
尝试改变了mdp文件中Parrinello-Rahman的tau_p,从1.0,2.0,4.0都改过,还是出现这个情况,请问有什么解决方法吗
PS:正在将步长缩短至0.0005,tau_p改为5.0进行新的尝试。附mpd文件文件!







mdp文件.png (190.75 KB, 下载次数 Times of downloads: 0)

mdp文件.png

shearMD.mdp

3 KB, 下载次数 Times of downloads: 1

849

帖子

4

威望

1920

eV
积分
2849

Level 5 (御坂)

A Student

2#
发表于 Post on 2025-9-2 21:36:16 | 只看该作者 Only view this author
不稳定通常先换回berendsen压浴试试。
敬仰一针见血的指责,厌倦别有用心的赞美。

3

帖子

0

威望

39

eV
积分
42

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-9-2 21:48:44 | 只看该作者 Only view this author
student0618 发表于 2025-9-2 21:36
不稳定通常先换回berendsen压浴试试。

谢谢回答,之前已经在相同体系用berendsen跑过不加deform的,是可以稳定跑完的,在RMSD在50ns可以稳定,我理解的是是要跑一下施加deform但是在berendsen压浴下的吗

849

帖子

4

威望

1920

eV
积分
2849

Level 5 (御坂)

A Student

4#
发表于 Post on 2025-9-2 22:13:42 | 只看该作者 Only view this author
7243007036 发表于 2025-9-2 21:48
谢谢回答,之前已经在相同体系用berendsen跑过不加deform的,是可以稳定跑完的,在RMSD在50ns可以稳定, ...

是。
敬仰一针见血的指责,厌倦别有用心的赞美。

9

帖子

0

威望

633

eV
积分
642

Level 4 (黑子)

5#
发表于 Post on 2025-9-21 11:11:23 | 只看该作者 Only view this author
请问mdp怎么设置的,能发一个文件参考一下吗,我设置的z轴deform得到的应力应变曲线一直有问题

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-25 07:52 , Processed in 0.250503 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list