计算化学公社

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

[GROMACS] 求共享固液界面或三相体系的模型

[复制链接 Copy URL]

6

帖子

0

威望

39

eV
积分
45

Level 2 能力者

本帖最后由 panyongcai 于 2022-4-19 17:36 编辑

初入gromacs,哪位有固液界面或三相体系的模型输入文件参考 一下,谢谢。


想模拟如下体系:盒子里下为固体表面,上为液体,液体中的气体分子自动形成气泡的过程。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2022-4-19 17:19:48 | 只看该作者 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

6

帖子

0

威望

39

eV
积分
45

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-4-19 17:35:41 | 只看该作者 Only view this author
本帖最后由 panyongcai 于 2022-4-19 17:37 编辑
sobereva 发表于 2022-4-19 17:19
说具体输入文件指什么

sob老师好,就是固体和液体,气体的pdb和对应生成的gro和top文件。类似图片中的体系,其中固体是固定的,液体填充满,一定个数的气体分子随机分布。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2022-4-19 18:15:06 | 只看该作者 Only view this author
构建初始结构,packmol就可以
分子动力学初始结构构建程序Packmol的使用
http://sobereva.com/473http://bbs.keinsci.com/thread-12549-1-1.html

但packmol填充液体区域不可能填充得非常致密。如果填充水,可以用gmx solvate基于自带的比如spc216.gro来填充,会填得非常致密。如果填充别的液体也想填充得特别紧密的话,需要先对相应的液体做NPT得到纯溶剂盒子的gro,再让gmx solvate依据它来填充别的体系的特定的空白区域。

top文件可以用sobtop
http://sobereva.com/soft/Sobtop
对石墨烯、每类分子都得到itp文件,恰当地include到主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

6

帖子

0

威望

39

eV
积分
45

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2022-4-19 18:51:04 | 只看该作者 Only view this author
sobereva 发表于 2022-4-19 18:15
构建初始结构,packmol就可以
分子动力学初始结构构建程序Packmol的使用
http://sobereva.com/473(http: ...

谢谢Sob老师详细的指导,我尽快学习尝试,非常感谢。

55

帖子

0

威望

536

eV
积分
591

Level 4 (黑子)

6#
发表于 Post on 2023-8-28 19:40:33 来自手机 | 只看该作者 Only view this author
sobereva 发表于 2022-4-19 18:15
构建初始结构,packmol就可以
分子动力学初始结构构建程序Packmol的使用
http://sobereva.com/473(http: ...

sob老师,我构建了一个命名为box.pdb的3*3*10的盒子,(0 0 0)->(3 3 3)是固体二氧化硅的界面模型。现在想在这个盒子的 (0 0 3)->(3 3 6) 填充水从而建立固液界面模型。sob老师,如何用
packmol进行对已有盒子的填充呢,如果packmol不行,如何建立这样的盒子呢。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

7#
发表于 Post on 2023-8-28 19:53:58 | 只看该作者 Only view this author
关于时间DE 发表于 2023-8-28 19:40
sob老师,我构建了一个命名为box.pdb的3*3*10的盒子,(0 0 0)->(3 3 3)是固体二氧化硅的界面模型。 ...

如果盒子当前(0 0 3)->(3 3 6)是真空区,直接用gmx solvate填充就完了,比用packmol填充得更致密
北京科音自然科学研究中心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

55

帖子

0

威望

536

eV
积分
591

Level 4 (黑子)

8#
发表于 Post on 2023-8-28 23:37:09 来自手机 | 只看该作者 Only view this author
sobereva 发表于 2023-8-28 19:53
如果盒子当前(0 0 3)->(3 3 6)是真空区,直接用gmx solvate填充就完了,比用packmol填充得更致密

收到,谢谢sob老师,但用solvate有一个问题,有个别水分子会插入到二氧化硅层里。
然后我发现建立的二氧化硅不是周期性结构,在盒子边缘会有断键,这个怎么解决呢。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

9#
发表于 Post on 2023-8-29 01:35:42 | 只看该作者 Only view this author
关于时间DE 发表于 2023-8-28 23:37
收到,谢谢sob老师,但用solvate有一个问题,有个别水分子会插入到二氧化硅层里。
然后我发现建立的二氧 ...

把scale设大一些,使得水插不进去孔洞缝隙





弄清楚是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

55

帖子

0

威望

536

eV
积分
591

Level 4 (黑子)

10#
发表于 Post on 2023-8-29 10:00:20 | 只看该作者 Only view this author
sobereva 发表于 2023-8-29 01:35
把scale设大一些,使得水插不进去孔洞缝隙

应该是盒子问题,跑em的时候会报错

106

帖子

0

威望

719

eV
积分
825

Level 4 (黑子)

11#
发表于 Post on 2024-9-1 13:21:07 | 只看该作者 Only view this author
关于时间DE 发表于 2023-8-28 19:40
sob老师,我构建了一个命名为box.pdb的3*3*10的盒子,(0 0 0)->(3 3 3)是固体二氧化硅的界面模型。 ...

请问你是如何使用packmol构建固体二氧化硅的界面模型的?我根据体积和密度计算二氧化硅分子数构建出来的好像是液相而非固体二氧化硅了

55

帖子

0

威望

536

eV
积分
591

Level 4 (黑子)

12#
发表于 Post on 2024-9-2 21:15:44 | 只看该作者 Only view this author
zhangfaxue 发表于 2024-9-1 13:21
请问你是如何使用packmol构建固体二氧化硅的界面模型的?我根据体积和密度计算二氧化硅分子数构建出来的 ...

二氧化硅晶体去构建

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 01:10 , Processed in 0.207964 second(s), 23 queries , Gzip On.

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