计算化学公社

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

[GROMACS] 求助:用alchemical_analysis计算联苯的溶剂化自由能得到诡异的重叠分布图该如何改进

[复制链接 Copy URL]

40

帖子

0

威望

261

eV
积分
301

Level 3 能力者

如题,我的自由能计算参数设置如下:
free-energy              = yes
couple-moltype           = BIP
couple-lambda0           = none
couple-lambda1           = vdw-q
couple-intramol          = no
separate-dhdl-file       = yes
sc-alpha                 = 0.5
sc-power                 = 1
sc-sigma                 = 0.3
init-lambda-state        = $LAMBDA$
coul-lambdas             = 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00
vdw-lambdas              = 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
nstdhdl                  = 100
calc-lambda-neighbors    = -1

用alchemical_analysis分析得到采样窗口的重叠分布图和相应的重叠矩阵如下,这样的采样窗口该如何改进



cfm.png (249.52 KB, 下载次数 Times of downloads: 46)

重叠分布图

重叠分布图

O_MBAR.png (253.8 KB, 下载次数 Times of downloads: 29)

重叠矩阵

重叠矩阵

12

帖子

0

威望

97

eV
积分
109

Level 2 能力者

2#
发表于 Post on 2022-2-28 18:32:29 | 只看该作者 Only view this author
我计算的时候也遇到过这种情况,是因为所取的窗口太密集了,减少几个窗口应该可以解决这个问题

40

帖子

0

威望

261

eV
积分
301

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-3-14 16:04:24 | 只看该作者 Only view this author
王家小姐 发表于 2022-2-28 18:32
我计算的时候也遇到过这种情况,是因为所取的窗口太密集了,减少几个窗口应该可以解决这个问题

请问,如何判断哪些个窗口应该减少呢?

12

帖子

0

威望

97

eV
积分
109

Level 2 能力者

4#
发表于 Post on 2022-5-16 16:59:09 | 只看该作者 Only view this author
活性炭 发表于 2022-3-14 16:04
请问,如何判断哪些个窗口应该减少呢?

您好,我判断的标准是对角线面积很大,可以减少些窗口,我是隔一个减少一个

40

帖子

0

威望

261

eV
积分
301

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-6-8 15:13:56 | 只看该作者 Only view this author
王家小姐 发表于 2022-5-16 16:59
您好,我判断的标准是对角线面积很大,可以减少些窗口,我是隔一个减少一个

好的,谢谢回复!

20

帖子

0

威望

441

eV
积分
461

Level 3 能力者

6#
发表于 Post on 2023-7-28 16:09:08 | 只看该作者 Only view this author
您好,想问下有解决吗?我看到有人说增加窗口,有人说减少窗口,我不知道该怎么搞了

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 05:36 , Processed in 0.236713 second(s), 23 queries , Gzip On.

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