计算化学公社

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

[综合交流] 从文献中总结的三种PMF计算方法

[复制链接 Copy URL]

516

帖子

1

威望

4703

eV
积分
5239

Level 6 (一方通行)

本帖最后由 tjuptz 于 2019-7-19 14:22 编辑

我转载的《杂谈自由能计算,PMF,伞形抽样,WHAM》帖子中对MD计算自由能从原理和具体实现上进行了讨论,但当时我没有看具体实例的文献,读完帖子也是模棱两可。最近读了两篇文献,其中涉及了三种方法来计算PMF进而获得分子非键聚合的ΔG。现整理如下:
一、通过径向分布函数RDF计算

二、限制力方法



三、伞状取样技术


随着模拟技术的发展,前两种方法可能因为采样不充分、准确度等问题现在已经少有人用了吧,但方法较简单从这些描述中即可弄懂操作,伞状取样技术比较先进但具体操作就比较复杂了。





评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
15810593730 + 4

查看全部评分 View all ratings

465

帖子

1

威望

2318

eV
积分
2803

Level 5 (御坂)

2#
发表于 Post on 2019-7-19 11:27:47 | 只看该作者 Only view this author
我见过有人用CP2K算化学反应时用第二种方法做thermodynamic integration的。主要是因为这样做不需要重新编译CP2K,因而很方便。如果在CP2K中做umbrella sampling的话,需要plumed。

1149

帖子

6

威望

6635

eV
积分
7904

Level 6 (一方通行)

3#
发表于 Post on 2019-7-19 12:53:05 | 只看该作者 Only view this author
不是很全面
PMF计算有几大类方法:
US类:US,REX-US,自适应US等
MtD类:MtD,WT-MtD,BE-MtD等,
ABF类:ABF,eABF,meta-eABF等,
SMD类:SMD+Jarzynski equality
TAMD类:TAMD+eABF estimator
constrain MD类:constrain MD

还有用广义系综方法(REMD,AMD,ITS等+Reweighting)也可以得到PMF

比较推荐用MtD或者ABF类下面的方法,如果不好选择反应坐标的话就用广义系综类的方法

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
tjuptz + 3

查看全部评分 View all ratings

553

帖子

0

威望

3324

eV
积分
3877

Level 5 (御坂)

4#
发表于 Post on 2019-7-19 13:44:48 | 只看该作者 Only view this author
补充一句,用g(r)计算其实就是在brute-force md下面采样,拿到概率分布,然后直接加log转化成自由能,非常直观,不算近似。

90

帖子

0

威望

330

eV
积分
420

Level 3 能力者

5#
发表于 Post on 2019-8-19 16:49:41 | 只看该作者 Only view this author
楼主,可以请问下关于三种计算方法的文献名字吗?感谢!
生活需要勇气,一定要坚持下去!

516

帖子

1

威望

4703

eV
积分
5239

Level 6 (一方通行)

6#
 楼主 Author| 发表于 Post on 2019-8-20 09:14:42 | 只看该作者 Only view this author
青青青 发表于 2019-8-19 16:49
楼主,可以请问下关于三种计算方法的文献名字吗?感谢!

Energy & Fuels 2009, 23, 1220–1229
J. Phys. Chem. B 2013, 117, 5765−5776

90

帖子

0

威望

330

eV
积分
420

Level 3 能力者

7#
发表于 Post on 2019-8-20 10:03:02 | 只看该作者 Only view this author
tjuptz 发表于 2019-8-20 09:14
Energy & Fuels 2009, 23, 1220–1229
J. Phys. Chem. B 2013, 117, 5765−5776

谢谢!
生活需要勇气,一定要坚持下去!

100

帖子

0

威望

1153

eV
积分
1253

Level 4 (黑子)

8#
发表于 Post on 2025-3-4 16:52:38 | 只看该作者 Only view this author
Daniel_Arndt 发表于 2019-7-19 11:27
我见过有人用CP2K算化学反应时用第二种方法做thermodynamic integration的。主要是因为这样做不需要重新编 ...

怎么用cp2k来实现伞形采样。

8

帖子

0

威望

1275

eV
积分
1283

Level 4 (黑子)

9#
发表于 Post on 2025-3-4 16:58:01 | 只看该作者 Only view this author
fhh2626 发表于 2019-7-19 12:53
不是很全面
PMF计算有几大类方法:
US类:US,REX-US,自适应US等

能否推荐几本书籍和综述,可以直观的补充一下您提到的这些内容的基础。谢谢。

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 08:22 , Processed in 0.183155 second(s), 24 queries , Gzip On.

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