计算化学公社

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

[GROMACS] 模拟体系经过nvt过程之后物质糅成一团

[复制链接 Copy URL]

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

老师们模拟的li离子在固态聚合物电解质中的传输,能量最小化后进行了nvt过程,这个过程持续了500ps,然后我检查生成的gro文件,发现如图所示,成为了一团团的,不知道我问题出现在哪?希望老师们帮忙看看,以下是我模拟文件图一是min之后的,图二是nvt之后的

文件.zip

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

25

帖子

0

威望

254

eV
积分
279

Level 3 能力者

2#
发表于 Post on 2025-7-30 20:06:59 | 只看该作者 Only view this author
我看不出你的分子是什么,但看上去是有机分子,那你NVT的情况下盒子体积不变,里面的物质肯定会倾向于团聚在一起啊,说明你一开始分子填充的很稀疏

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-7-30 20:14:19 | 只看该作者 Only view this author
KazusaT 发表于 2025-7-30 20:06
我看不出你的分子是什么,但看上去是有机分子,那你NVT的情况下盒子体积不变,里面的物质肯定会倾向于团聚 ...

是的,一开始建立太小的话,我投放的物质有的放不进去,所以就把盒子尺寸调大了

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2025-7-30 21:55:10 | 只看该作者 Only view this author
KazusaT 发表于 2025-7-30 20:06
我看不出你的分子是什么,但看上去是有机分子,那你NVT的情况下盒子体积不变,里面的物质肯定会倾向于团聚 ...

那老师我接下来是不是可以进行npt平衡,让他密一些呢

25

帖子

0

威望

254

eV
积分
279

Level 3 能力者

5#
发表于 Post on 2025-7-31 04:52:02 | 只看该作者 Only view this author
可以直接进行NPT平衡,让体系逐步弛豫到一个合适的密度

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-7-31 10:52:44 | 只看该作者 Only view this author
本帖最后由 huangzi 于 2025-7-31 10:54 编辑
KazusaT 发表于 2025-7-31 04:52
可以直接进行NPT平衡,让体系逐步弛豫到一个合适的密度

老师我昨天npt跑了10ns,然后得出来的结构还是有很多空隙,因为我想填满我整个盒子,出现这种情况是我npt
时间不够吗还是哪的设置有问题啊,而且图中圈的这一部分物质也没有融合到团簇里一直是一个独立的

微信图片_2025-07-31_105139_446.jpg (63.12 KB, 下载次数 Times of downloads: 0)

微信图片_2025-07-31_105139_446.jpg

npt.mdp

2.39 KB, 下载次数 Times of downloads: 0

86

帖子

0

威望

2954

eV
积分
3040

Level 5 (御坂)

7#
发表于 Post on 2025-7-31 11:02:43 | 只看该作者 Only view this author
huangzi 发表于 2025-7-31 10:52
老师我昨天npt跑了10ns,然后得出来的结构还是有很多空隙,因为我想填满我整个盒子,出现这种情况是我npt ...

如果你直接做npt,还是会团在一起的。如果只是想把空隙消除,把ref_p提高到100,做高温高压的退火处理。

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2025-7-31 11:35:20 | 只看该作者 Only view this author
liuxdhs 发表于 2025-7-31 11:02
如果你直接做npt,还是会团在一起的。如果只是想把空隙消除,把ref_p提高到100,做高温高压的退火处理。

好的老师,体系形成团簇是由于没进行退火处理吗?目前我i的想法是先消除空隙,让体系看起来均匀一些,因为按照我之前模拟的10ns的npt体系是没有平衡的,我猜测可能是由于形成的团簇引起的

6万

帖子

99

威望

5万

eV
积分
120064

管理员

公社社长

9#
发表于 Post on 2025-7-31 12:50:12 | 只看该作者 Only view this author
体系形成团簇是由于没进行退火处理吗?
根本毫无联系
分子间存在弱相互作用,凭什么分子间不会通过吸引作用聚在一起?凭什么会一直保持分散的这种能量较高的状态?
北京科音自然科学研究中心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

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2025-7-31 14:56:06 | 只看该作者 Only view this author
sobereva 发表于 2025-7-31 12:50
根本毫无联系
分子间存在弱相互作用,凭什么分子间不会通过吸引作用聚在一起?凭什么会一直保持分散的这种 ...

那老师我这种算是正常的状态吗,接下来直接进行退火处理将体系消除空隙,然后平衡吗?我在想是不是因为我采用insert-molecules命令插分子的时候,建立的盒子大了,接下来我准备用packmol建立盒子,适当缩小一下,空隙也会小点

6万

帖子

99

威望

5万

eV
积分
120064

管理员

公社社长

11#
发表于 Post on 2025-7-31 14:59:49 | 只看该作者 Only view this author
huangzi 发表于 2025-7-31 14:56
那老师我这种算是正常的状态吗,接下来直接进行退火处理将体系消除空隙,然后平衡吗?我在想是不是因为我 ...

消除空隙靠的是控压,控压把密度弄合理了,自然就没有空隙了
北京科音自然科学研究中心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

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2025-7-31 15:05:41 | 只看该作者 Only view this author
sobereva 发表于 2025-7-31 14:59
消除空隙靠的是控压,控压把密度弄合理了,自然就没有空隙了

那老师接下来还是跑npt控压吧,可能是我初始盒子建立的太大了,npt10ns之后还是分散的簇,接下来用packmol建盒子是不是会好点啊

25

帖子

0

威望

254

eV
积分
279

Level 3 能力者

13#
发表于 Post on 2025-7-31 17:26:07 | 只看该作者 Only view this author
huangzi 发表于 2025-7-31 15:05
那老师接下来还是跑npt控压吧,可能是我初始盒子建立的太大了,npt10ns之后还是分散的簇,接下来用packmo ...

如果你的聚电解质是直链的,你可以用packmol不断增加分子数目,直到难以产生盒子或者堆积成为一个比较合理的状态,如果你知道聚电解质的密度也可以先稍微算一下多少体积大概对应多少分子
直接npt的话肯定会分散开的,分子会由于相互作用吸引团聚在一起,这是降低能量的过程。

39

帖子

0

威望

149

eV
积分
188

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2025-7-31 20:39:50 | 只看该作者 Only view this author
KazusaT 发表于 2025-7-31 17:26
如果你的聚电解质是直链的,你可以用packmol不断增加分子数目,直到难以产生盒子或者堆积成为一个比较合 ...

好的老师,我用packmol重新建盒子了

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 23:04 , Processed in 0.361930 second(s), 23 queries , Gzip On.

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