计算化学公社

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

[GROMACS] 求助MD模拟自组装5ns后分子间相互作用能为何大于0

[复制链接 Copy URL]

82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

我只选取了4个有机分子进行自组装模拟,将有机分子吸附质划分为MOL组、石墨基底分为GRA组,分子也是经过优化后的,分子力场是从http://zarbi.chem.yale.edu/ligpargen/这个网站获取的oplsaa力场,进行了能量最小化以及位置限制的平衡相模拟后做了无限制的产生相模拟,各个过程中自组装的构象都保持的很好,但是对于吸附质分子间的相互作用能 LJ-SR(MOL-MOL)一直大于零。之前做过类似分子的模拟,两者间构象的区别是产生相跑完后此分子间O—H氢键的键长较小,平均大概1.8埃。请问是什么原因导致分子间相互作用能大于0的呢?感谢!

力场文件.zip

137.05 KB, 下载次数 Times of downloads: 1

system.gro

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

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2022-4-17 03:37:02 | 只看该作者 Only view this author
MOL-MOL不仅体现分子间,也同时体现分子内作用
你可以先模拟单个被吸附的分子看是什么情况(或者把单个分子定义为能量组然后rerun,看组内的各项能量)
为正没有必然的异常,只要模拟结果合理完全不用管它
北京科音自然科学研究中心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

82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-4-17 21:21:14 | 只看该作者 Only view this author
sobereva 发表于 2022-4-17 03:37
MOL-MOL不仅体现分子间,也同时体现分子内作用
你可以先模拟单个被吸附的分子看是什么情况(或者把单个分 ...

谢谢sob老师,我用rerun计算了单个分子的能量(ONE-ONE),与之前4个分子的MOL-MOL进行了比较,请问想得到4个分子间的范德华相互作用大小是取LJ-SR(MOL-MOL)与4倍LJ-SR(ONE-ONE)的差值吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2022-4-17 22:31:23 | 只看该作者 Only view this author
Rosefinch 发表于 2022-4-17 21:21
谢谢sob老师,我用rerun计算了单个分子的能量(ONE-ONE),与之前4个分子的MOL-MOL进行了比较,请问想得 ...

当前模拟中,四个单元的非键相互作用能不可能恰好精确相同。应当把每个单元都定义成能量组,用整体的减去这四个的
北京科音自然科学研究中心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

82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-4-18 00:58:13 | 只看该作者 Only view this author
sobereva 发表于 2022-4-17 22:31
当前模拟中,四个单元的非键相互作用能不可能恰好精确相同。应当把每个单元都定义成能量组,用整体的减去 ...

请问老师,我又计算了另外三个分子的LJ-SR,但是发现它们四个的LJ-SR加起来才69.8714,比整体的LJ-SR(MOL-MOL)=87.1813还要小,也就是减出来值依旧是正的,请问这是什么原因呀?难道是范德华作用对于这四个分子的结合起副作用吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2022-4-18 02:03:08 | 只看该作者 Only view this author
Rosefinch 发表于 2022-4-18 00:58
请问老师,我又计算了另外三个分子的LJ-SR,但是发现它们四个的LJ-SR加起来才69.8714,比整体的LJ-SR(MOL ...

有可能
因为形成氢键会导致相应的原子间距离较近,交换互斥作用明显

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
Rosefinch + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 06:16 , Processed in 0.184860 second(s), 26 queries , Gzip On.

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