计算化学公社

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

[GROMACS] 求助:gmx跑em时emtol无法达到100KJ/(mol*nm)

[复制链接 Copy URL]

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

本帖最后由 confidence 于 2022-4-19 16:58 编辑

为模拟水包油乳液在氨基化碳纳米管和亲水型二氧化硅中的分离情况,构建了如图所示模型。


其中1是树脂,2是表面接-Si(OCH3)3的二氧化硅,3是水包油乳液,乳液所用油为异辛烷,表面活性剂为十二烷基硫酸钠,4是表面接-NH3的碳纳米管。

异辛烷、十二烷基硫酸钠、树脂的itp文件从LigParGen网站获取,二氧化硅、碳纳米管的top文件使用x2top获取,采用力场皆为OPLS-AA。

跑能量最小化采用的mdp文件如下图所示:


跑em时不论固定二氧化硅+碳纳米管与否,都无法达到收敛标准,看过sob老师发的帖子,感觉应该是模型不合理或者其他因素,但本人尝试其他模型搭建也未能成功。

希望各位老师能给我提点意见,看看究竟是哪里出了问题,谢谢各位老师
拓扑,mdp文件.rar (1.06 MB, 下载次数 Times of downloads: 12) 这是模型用到的所有itp,mdp,top及结构文件。

在跑em时用的是冻结,后跑npt时改为位置限制。



6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

2#
发表于 Post on 2022-3-23 14:44:07 | 只看该作者 Only view this author
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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-3-23 16:13:08 | 只看该作者 Only view this author
sobereva 发表于 2022-3-23 14:44
MD能跑就行,没必要强求达到当前的最大受力收敛标准

感谢sob老师的回答,目前我这建模emol连5000都收敛不了,跑npt时会出现很多的LINCS提示,然后就停止跑了。看过您给的原因估计是建模不合理,但本人看不出究竟是哪里出的问题,能麻烦sob老师看一下可能是哪里有问题么

6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

4#
发表于 Post on 2022-4-12 01:29:31 | 只看该作者 Only view this author
confidence 发表于 2022-3-23 16:13
感谢sob老师的回答,目前我这建模emol连5000都收敛不了,跑npt时会出现很多的LINCS提示,然后就停止跑了 ...

http://sobereva.com/soft/Sobtop#FAQ8仔细看了
当前这模型很怪,水都朝着一个方向
北京科音自然科学研究中心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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-4-15 16:22:49 | 只看该作者 Only view this author
sobereva 发表于 2022-4-12 01:29
把http://sobereva.com/soft/Sobtop#FAQ8仔细看了
当前这模型很怪,水都朝着一个方向

谢谢sob老师的回答,这两天我把模型各部分拆开跑发现。乳液+水能跑npt,二氧化硅在真空条件下也能跑npt。但一旦二氧化硅+水/二氧化硅+树脂就跑不了npt了,请问老师这有可能是什么原因呢?

6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

6#
发表于 Post on 2022-4-15 19:33:56 | 只看该作者 Only view this author
confidence 发表于 2022-4-15 16:22
谢谢sob老师的回答,这两天我把模型各部分拆开跑发现。乳液+水能跑npt,二氧化硅在真空条件下也能跑npt。 ...

检查拓扑文件和结构文件对应关系,检查top文件字段顺序,检查初始结构合理性等等
北京科音自然科学研究中心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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-4-18 18:51:04 | 只看该作者 Only view this author
sobereva 发表于 2022-4-15 19:33
检查拓扑文件和结构文件对应关系,检查top文件字段顺序,检查初始结构合理性等等

嗯嗯谢谢sob老师,字段顺序和对应关系应该是没有问题的,有仔细检查过。但是这个初始结构合理性我不知道怎么确认,老师您有什么标准么?

6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

8#
发表于 Post on 2022-4-19 00:56:02 | 只看该作者 Only view this author
confidence 发表于 2022-4-18 18:51
嗯嗯谢谢sob老师,字段顺序和对应关系应该是没有问题的,有仔细检查过。但是这个初始结构合理性我不知道 ...

http://sobereva.com/soft/Sobtop#FAQ8里已经说了,主要检查是否有不合理接触等
北京科音自然科学研究中心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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2022-4-19 15:23:40 | 只看该作者 Only view this author
sobereva 发表于 2022-4-19 00:56
http://sobereva.com/soft/Sobtop#FAQ8里已经说了,主要检查是否有不合理接触等

嗯嗯,谢谢sob老师,我对着您的链接中所提出的方法自己检查对照了,但是还是未能找出究竟问题出在哪里。我把所用到的所有文件都上传了,能否麻烦老师您帮忙看看是哪里有问题呢?

14

帖子

0

威望

223

eV
积分
237

Level 3 能力者

10#
发表于 Post on 2024-3-29 19:44:22 | 只看该作者 Only view this author
您好楼主,想请教一个问题,我看您水包油乳液这部分,是用油球做为油滴放在水相中间了,但是一般水包油乳液中的油滴直径在微米级别,而这里的油球直径在纳米级别,请问这部分尺度换算应该怎么考虑呢?可以用这个模型来代表水包油乳液吗?

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2024-3-30 14:50:18 | 只看该作者 Only view this author
doh88330 发表于 2024-3-29 19:44
您好楼主,想请教一个问题,我看您水包油乳液这部分,是用油球做为油滴放在水相中间了,但是一般水包油乳液 ...

我这里油球的直径是几纳米,主要是由我整个盒子的大小决定,而且在这个尺度下就已经能模拟出油滴的分离情况了。

14

帖子

0

威望

223

eV
积分
237

Level 3 能力者

12#
发表于 Post on 2024-3-31 10:41:28 | 只看该作者 Only view this author
confidence 发表于 2024-3-30 14:50
我这里油球的直径是几纳米,主要是由我整个盒子的大小决定,而且在这个尺度下就已经能模拟出油滴的分离情 ...

但和实际的体系中油滴的大小尺寸数量级相差较大呀,这个尺度上的换算需要考虑吗?我现在模拟也遇到这个问题

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 01:24 , Processed in 0.171306 second(s), 23 queries , Gzip On.

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