计算化学公社

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

[GROMACS] 求助,做PMF时,直方图做不出来,什么原因

[复制链接 Copy URL]

48

帖子

0

威望

451

eV
积分
499

Level 3 能力者

本帖最后由 wgg1181543722 于 2021-1-29 17:15 编辑

我最近还是在做pmf,可无论怎么做,做出来的直方图都是一个峰,PMF图也是乱糟糟的,实在是没办法,在这里求助一下大佬,麻烦你们抽空帮我看一下问题出在哪里,
注:我是严格按照附件的文献做的,附件还包含了模拟用到的mdp文件

拉伸示意图.png (6.51 KB, 下载次数 Times of downloads: 52)

拉伸示意图

拉伸示意图

拉伸的物质.png (17.82 KB, 下载次数 Times of downloads: 43)

拉伸的物质

拉伸的物质

文献的步骤.png (177.97 KB, 下载次数 Times of downloads: 38)

文献的步骤

文献的步骤

文献的PMF.png (83.19 KB, 下载次数 Times of downloads: 39)

文献的PMF

文献的PMF

文献的直方图.png (56.98 KB, 下载次数 Times of downloads: 37)

文献的直方图.png

我做的-力随时间的变化.png (11.22 KB, 下载次数 Times of downloads: 39)

我做的-力随时间的变化.png

我做的-PMF.png (10.08 KB, 下载次数 Times of downloads: 42)

我做的-PMF.png

我做的-直方图.png (10.05 KB, 下载次数 Times of downloads: 42)

我做的-直方图.png

Molecular insights into multilayer 18-crown-6-.pdf

2.34 MB, 下载次数 Times of downloads: 22

SI.pdf

303.73 KB, 下载次数 Times of downloads: 2

pull.mdp

1.07 KB, 下载次数 Times of downloads: 9

npt_umbrella.mdp

1.07 KB, 下载次数 Times of downloads: 7

md_umbrella.mdp

1.03 KB, 下载次数 Times of downloads: 8

561

帖子

0

威望

3410

eV
积分
3971

Level 5 (御坂)

2#
发表于 Post on 2021-1-29 18:41:35 | 只看该作者 Only view this author
这文献里什么破data,error搞这么大……

513

帖子

0

威望

4388

eV
积分
4901

Level 6 (一方通行)

truffle

3#
发表于 Post on 2021-6-2 12:08:59 | 只看该作者 Only view this author
k64_cc 发表于 2021-1-29 18:41
这文献里什么破data,error搞这么大……

想请教下您,对于PMF,约1 kJ/mol的error bar算大吗
No problem is insoluble in all conceivable circumstances.

561

帖子

0

威望

3410

eV
积分
3971

Level 5 (御坂)

4#
发表于 Post on 2021-6-2 16:08:42 | 只看该作者 Only view this author
naoki 发表于 2021-6-2 12:08
想请教下您,对于PMF,约1 kJ/mol的error bar算大吗

对于这个文献来讲,能垒也就2kJ/mol上下,1kJ的误差就有些不可接受了。如果降低误差有难度还则罢了,但是对于US来说,降低误差就是一个力大砖飞的活,多sample就好了。需要做,可以做,但是没做,就不太行。

513

帖子

0

威望

4388

eV
积分
4901

Level 6 (一方通行)

truffle

5#
发表于 Post on 2021-6-2 16:49:43 | 只看该作者 Only view this author
k64_cc 发表于 2021-6-2 16:08
对于这个文献来讲,能垒也就2kJ/mol上下,1kJ的误差就有些不可接受了。如果降低误差有难度还则罢了,但是 ...

明白了,谢谢~
No problem is insoluble in all conceivable circumstances.

513

帖子

0

威望

4388

eV
积分
4901

Level 6 (一方通行)

truffle

6#
发表于 Post on 2021-6-3 14:17:06 | 只看该作者 Only view this author
k64_cc 发表于 2021-6-2 16:08
对于这个文献来讲,能垒也就2kJ/mol上下,1kJ的误差就有些不可接受了。如果降低误差有难度还则罢了,但是 ...

我还想请教一下您,我跑三个平行的拉伸动力学准备做error bar时发现,三组平行的PMF曲线的横坐标距离不是完全一样的,比如第一组在0.25 nm处有个PMF值,第二组在0.26处有值,这种情况怎么做error bar呢,谢谢!
No problem is insoluble in all conceivable circumstances.

33

帖子

0

威望

1087

eV
积分
1120

Level 4 (黑子)

7#
发表于 Post on 2024-7-24 09:23:21 | 只看该作者 Only view this author
您好,能发一下您的gro,top,itp文件嘛?我想做一下复现

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 17:39 , Processed in 4.525913 second(s), 23 queries , Gzip On.

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