计算化学公社

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

[GROMACS] Gromacs如何计算同组分子内的相互作用能?

[复制链接 Copy URL]

2

帖子

0

威望

327

eV
积分
329

Level 3 能力者

本帖最后由 yuancheng 于 2020-5-17 23:02 编辑

作为一名初学者,目前刚学会了gmx mdrun方法计算两个不同能量组之间的相互作用能,即如图所示*ndx中a与b之间的。

gmx energy命令后选取:21 Coul-SR:a-b、22 LJ-SR:a-b 两相数据作为短程库伦相互作用和短程范德华相互作用进行分析。

然而如果我还想继续分析a组中(310个原子代表5个分子)分子间的相互作用能时就遇到了一定的困难,之前我以为就是 17  Coul-SR:a-a   18  LJ-SR:a-a 两个数据,但是看到算出的结果明显感觉不对(正负号感觉都是错的),对于分子数少的我可以采取这个5个分子分别建一个ndx中的组,然后分别算出每一对的Coul-SR:a-b、LJ-SR:a-b最后进行加和得到总的。

但是如果我的a组中分子数是上百个的话,这种手动计算就显得非常耗时且容易出错,显得非常不现实。我也想到了可能可以用一些计算脚本来算,但是苦于能力有限还编写不了脚本。

我想请教的是:1、Coul-SR:a-a 和LJ-SR:a-a这种选项代表什么意思?2、如何便捷地计算同一个组内分子间的相互作用能?有没有类似的脚本可以参考呢?


1.PNG (26.55 KB, 下载次数 Times of downloads: 26)

ndx

ndx

2.PNG (35.08 KB, 下载次数 Times of downloads: 22)

2.PNG

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

2#
发表于 Post on 2020-5-17 23:59:12 | 只看该作者 Only view this author
Coul-SR:a-a 和LJ-SR:a-a就是a组内的实空间计算的库仑作用和范德华作用,分子内和分子间的非键作用都包括

做MD的不会写脚本寸步难行,不会就赶紧学
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1

帖子

0

威望

13

eV
积分
14

Level 1 能力者

3#
发表于 Post on 2021-8-25 16:19:25 | 只看该作者 Only view this author
请问您的问题解决了吗?我也遇到了相同的问题。如果挨个建索引,工作内容太大了

                                               

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

4#
发表于 Post on 2021-9-8 06:15:07 | 只看该作者 Only view this author
lyj980317 发表于 2021-8-25 16:19
请问您的问题解决了吗?我也遇到了相同的问题。如果挨个建索引,工作内容太大了

看此文了解点脚本编写常识
详谈Multiwfn的命令行方式运行和批量运行的方法
http://sobereva.com/612http://bbs.keinsci.com/thread-24929-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

5

帖子

0

威望

37

eV
积分
42

Level 2 能力者

5#
发表于 Post on 2022-9-26 08:07:29 | 只看该作者 Only view this author
学姐好,我想问下有AB的能量文件和A和B各自的能量文件,怎用使用gmx mdrun方法计算两个不同能量组之间的相互作用能啊,谢谢啦!

109

帖子

0

威望

499

eV
积分
608

Level 4 (黑子)

6#
发表于 Post on 2025-9-13 10:59:00 | 只看该作者 Only view this author
sobereva 发表于 2020-5-17 23:59
Coul-SR:a-a 和LJ-SR:a-a就是a组内的实空间计算的库仑作用和范德华作用,分子内和分子间的非键作用都包括
...

请问sob老师这个实空间是考虑到组分和其周围环境的影响吗,因为我计算的离子液体,计算组分内正离子之间的相互作用出现了负值的情况,不知道这是不是有可能出现的情况?

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

7#
发表于 Post on 2025-9-13 23:42:27 | 只看该作者 Only view this author
wanzongliang 发表于 2025-9-13 10:59
请问sob老师这个实空间是考虑到组分和其周围环境的影响吗,因为我计算的离子液体,计算组分内正离子之间 ...

Ewald及其变体诸如PME,算的静电作用是考虑到无穷远距离的,是实空间部分和倒易空间部分之和,只有实空间部分才能做组的分解。“实空间”和“周围环境”概念上没直接联系
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

109

帖子

0

威望

499

eV
积分
608

Level 4 (黑子)

8#
发表于 Post on 2025-9-15 08:42:15 | 只看该作者 Only view this author
sobereva 发表于 2025-9-13 23:42
Ewald及其变体诸如PME,算的静电作用是考虑到无穷远距离的,是实空间部分和倒易空间部分之和,只有实空间 ...

好的感谢sob老师

1

帖子

0

威望

17

eV
积分
18

Level 1 能力者

9#
发表于 Post on 2025-10-25 16:50:27 | 只看该作者 Only view this author
wanzongliang 发表于 2025-9-13 10:59
请问sob老师这个实空间是考虑到组分和其周围环境的影响吗,因为我计算的离子液体,计算组分内正离子之间 ...

你好,请问这个问题解决了吗?我也是在溶液环境中,整体带正电的离子组分的Coul-SR 、 LJ-SR 都是负值,这是正常的吗?

109

帖子

0

威望

499

eV
积分
608

Level 4 (黑子)

10#
发表于 Post on 2025-11-5 14:45:57 | 只看该作者 Only view this author
小熊的立可白 发表于 2025-10-25 16:50
你好,请问这个问题解决了吗?我也是在溶液环境中,整体带正电的离子组分的Coul-SR 、 LJ-SR 都是负值, ...

我没有找到这个问题的原因,我自己写了分析脚本计算电荷之间的相互作用得到的结果是预想的结果

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 07:31 , Processed in 0.309259 second(s), 29 queries , Gzip On.

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