计算化学公社

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

[Lammps] 伞状抽样如何做出看窗口重叠情况的图?

[复制链接 Copy URL]

33

帖子

0

威望

229

eV
积分
262

Level 3 能力者

各位老师和同学,我在lammps中用fix colvars伞状抽样,然后用WHAM得到PMF,那如何做出窗口overlap的图以判断窗口重叠情况呢?

微信图片_20190227204004.jpg (189.74 KB, 下载次数 Times of downloads: 46)

微信图片_20190227204004.jpg

71

帖子

0

威望

6801

eV
积分
6872

Level 6 (一方通行)

2#
发表于 Post on 2019-2-27 22:08:24 | 只看该作者 Only view this author
wham的hist选项也可以输出下面那个图

561

帖子

0

威望

3410

eV
积分
3971

Level 5 (御坂)

3#
发表于 Post on 2019-2-28 15:00:13 | 只看该作者 Only view this author
Count your CV and draw histogram...

33

帖子

0

威望

229

eV
积分
262

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2019-3-1 14:01:59 | 只看该作者 Only view this author
xiaowu759 发表于 2019-2-27 22:08
wham的hist选项也可以输出下面那个图

可以吗?我仔细看了WHAM手册,也加进去了,但是输出是这样的(如下图),好像一个窗口就得好多点吧?我的为什么只有这些呢?

捕获.JPG (114.48 KB, 下载次数 Times of downloads: 44)

捕获.JPG

71

帖子

0

威望

6801

eV
积分
6872

Level 6 (一方通行)

5#
发表于 Post on 2019-3-1 16:46:11 | 只看该作者 Only view this author
wasngsimin 发表于 2019-3-1 14:01
可以吗?我仔细看了WHAM手册,也加进去了,但是输出是这样的(如下图),好像一个窗口就得好多点吧?我的 ...

看错了,以为你用的是gromacs里的g_wham呢

33

帖子

0

威望

229

eV
积分
262

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2019-3-1 19:07:56 | 只看该作者 Only view this author
xiaowu759 发表于 2019-3-1 16:46
看错了,以为你用的是gromacs里的g_wham呢

嗯嗯 我用的是Lammps

316

帖子

0

威望

1253

eV
积分
1569

Level 5 (御坂)

7#
发表于 Post on 2025-10-6 17:25:38 | 只看该作者 Only view this author
可以写个小脚本,把反应坐标序列输入到一个文件里,然后复制到origin中绘制分布图,再把分布图中的直方图部分隐藏就是这个了
######
顺便请教下各位,这个重叠什么程度比较好呢?如果重叠部分非常大,甚至接近90%的话是不是不太合理?是否要增加力常数

156

帖子

0

威望

1637

eV
积分
1793

Level 5 (御坂)

8#
发表于 Post on 2025-10-6 17:40:17 | 只看该作者 Only view this author
ljh123 发表于 2025-10-6 17:25
可以写个小脚本,把反应坐标序列输入到一个文件里,然后复制到origin中绘制分布图,再把分布图中的直方图部 ...

重叠90%,说明你可以减少一些windows,
当然提高力学常数也可以。。。

316

帖子

0

威望

1253

eV
积分
1569

Level 5 (御坂)

9#
发表于 Post on 2025-10-6 17:51:38 | 只看该作者 Only view this author
beyond 发表于 2025-10-6 17:40
重叠90%,说明你可以减少一些windows,
当然提高力学常数也可以。。。

可以请教下您重叠度过高的影响是什么吗?除了增加了不必要的计算量以外,会影响自由能曲线吗

156

帖子

0

威望

1637

eV
积分
1793

Level 5 (御坂)

10#
发表于 Post on 2025-10-6 18:10:56 | 只看该作者 Only view this author
ljh123 发表于 2025-10-6 17:51
可以请教下您重叠度过高的影响是什么吗?除了增加了不必要的计算量以外,会影响自由能曲线吗

不会
30%是比较合理的

35

帖子

0

威望

360

eV
积分
395

Level 3 能力者

11#
发表于 Post on 2025-11-16 23:54:56 | 只看该作者 Only view this author
伞形采样的窗口重叠图其实就是被“弹簧”牵引的物质的密度分布。所以我直接用lammps计算这个物质的密度分布,然后输出其密度分布和坐标关系的图就好了。

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 05:59 , Processed in 0.211238 second(s), 23 queries , Gzip On.

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