计算化学公社

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

[GROMACS] 设置了碳纳米管的周期性,为什么聚合物还是会跑到碳管内部?

[复制链接 Copy URL]

199

帖子

0

威望

1165

eV
积分
1364

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 NGC626 于 2025-8-14 22:07 编辑

各位老师好,我在模拟碳纳米管和聚合物吸附的动力学过程时,对碳管设置了周期性,但是在动力学过程后,聚合物仍然会跑到碳管内部,请问是mdp文件中的设置没有生效吗?如果不是的话,请问可能是什么原因导致的呢?体系中有两根碳管,一根弛豫,另一根是完全放开的,分别用AA和BB表示。

CNT.7z

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

社会浮躁

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

2#
发表于 Post on 2025-8-14 22:11:46 | 只看该作者 Only view this author
我没时间具体看这些文件。跑到碳管内部,要么是盒子尺寸不对,导致纳米管实际上不是周期性的,聚合物从两端进入了。VMD里把盒子显示出来看轨迹便知
要么是初始结构本来聚合物就在盒子里
要么是LJ参数不对,导致纳米管没表现出对聚合物位阻作用
没什么其它可能了
北京科音自然科学研究中心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

199

帖子

0

威望

1165

eV
积分
1364

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-8-14 22:20:43 | 只看该作者 Only view this author
本帖最后由 NGC626 于 2025-8-14 22:29 编辑
sobereva 发表于 2025-8-14 22:11
我没时间具体看这些文件。跑到碳管内部,要么是盒子尺寸不对,导致纳米管实际上不是周期性的,聚合物从两端 ...

谢谢sob老师回复,在初始建模的时候,聚合物没有摆放在在碳管中,各分子的拓扑文件也是利用sobtop生成,力场参数应该也没有问题,动力学轨迹里显示,在0.5ns的时候,聚合物的支链就进入了碳管内,体系中的碳管长度大约是2nm,之前看到您的这条回复(盒子尺寸),所以我把盒子的尺寸设得比较大,是7*7*8,请问需要把盒子的尺寸改小一点吗?
社会浮躁

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

4#
发表于 Post on 2025-8-14 22:43:03 | 只看该作者 Only view this author
NGC626 发表于 2025-8-14 22:20
谢谢sob老师回复,在初始建模的时候,聚合物没有摆放在在碳管中,各分子的拓扑文件也是利用sobtop生成, ...

显然这么大盒子会导致这是有限长度的纳米管片段,而不是无限长度的周期性纳米管
顺着纳米管方向的盒子必须和纳米管尺寸一致
北京科音自然科学研究中心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

199

帖子

0

威望

1165

eV
积分
1364

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-8-14 22:59:49 | 只看该作者 Only view this author
sobereva 发表于 2025-8-14 22:43
显然这么大盒子会导致这是有限长度的纳米管片段,而不是无限长度的周期性纳米管
顺着纳米管方向的盒子必 ...

好的,谢谢sob老师
社会浮躁

199

帖子

0

威望

1165

eV
积分
1364

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 7 day ago | 只看该作者 Only view this author
本帖最后由 NGC626 于 2025-8-15 21:29 编辑
sobereva 发表于 2025-8-14 22:43
显然这么大盒子会导致这是有限长度的纳米管片段,而不是无限长度的周期性纳米管
顺着纳米管方向的盒子必 ...

sob老师好,我对碳纳米管和聚合物吸附的体系做动力学模拟时发现,在对碳管设置了周期性后,如果冻结碳管的话,聚合物不会跑到碳管的内部;如果放开碳管的话,聚合物则会跑到碳管的内部,请问这样的现象正常吗?请问是否是因为放开的时候,周期性的设置不会生效呢?
社会浮躁

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

7#
发表于 Post on 6 day ago | 只看该作者 Only view this author
NGC626 发表于 2025-8-15 21:25
sob老师好,我对碳纳米管和聚合物吸附的体系做动力学模拟时发现,在对碳管设置了周期性后,如果冻结碳管 ...

VMD里显示出盒子,从头到尾观看轨迹了解分子是怎么运动的,自然就知道怎么回事了
北京科音自然科学研究中心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

199

帖子

0

威望

1165

eV
积分
1364

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 6 day ago | 只看该作者 Only view this author
本帖最后由 NGC626 于 2025-8-16 14:26 编辑
sobereva 发表于 2025-8-16 02:35
VMD里显示出盒子,从头到尾观看轨迹了解分子是怎么运动的,自然就知道怎么回事了

谢谢sob老师,我仔细查看了各分子的运动轨迹,发现碳管在运动的时候,当它开始在盒子里倾斜,也就是它的长度远小于某一条边的时候,就会有小分子进入内部,而如果碳管固定不动,也就是其长度方向一直和盒子最短的那条边平行的时候,小分子就没法进去,请问这是因为运动导致碳管在盒子内的某个方向上变成有限长度了吗?
社会浮躁

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

9#
发表于 Post on 5 day ago | 只看该作者 Only view this author
NGC626 发表于 2025-8-16 14:14
谢谢sob老师,我仔细查看了各分子的运动轨迹,发现碳管在运动的时候,当它开始在盒子里倾斜,也就是它的 ...

根据你的描述,貌似是盒子尺寸不对,或者没用periodic-molecule=yes,或者bond项描述的拓扑关系本身就不对,导致gromacs没考虑跨盒子的C-C键,因此纳米管才会像孤立分子一样自发旋转

PS: 如果参加过北京科音分子动力学与GROMACS培训班(http://www.keinsci.com/KGMX)的话,仔细看里面“纳米球、管、板的模拟” 部分的讲解和例子,了解怎么正确模拟跨盒子的无限延展的体系
北京科音自然科学研究中心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, 2025-8-22 21:08 , Processed in 0.304556 second(s), 23 queries , Gzip On.

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