计算化学公社

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

[GROMACS] 用MS构建高分子聚合物周期性体系如何用于Groamcs模拟

[复制链接 Copy URL]

43

帖子

0

威望

175

eV
积分
218

Level 3 能力者

各位老师前辈,我最近需要做一个PET(聚对苯二甲酸乙二醇酯)晶面与蛋白质界面的分子动力学模拟。成功按照网上教程在MS中构建出来了PET的晶面结构,目前遇到的问题是导出为PDB格式,整体的键连关系显得就会很混乱,在帖子中看有人用在VMD中用Dynamic Bond,发现没有周期性键连关系,没法像Sobtop老师在Sobtop教程中例子那样有两端连接的键。
在网上看到有前辈说导成为cif格式在在Gaussview中转换成mol2文件,用Gaussview打开后氢原子与其他原子如苯环碳没有键连关系,并且苯环的键连关系错误,出现三键和双键。手动更改耗时很长,总改不全。
请问一下我该怎么改才能将这个结构应用到Gromacs运行呢,Sobtop老师开发的sobtop软件的一些例子我也在看,想着先构建出mol2文件,使用PDB文件重新定义各个类型的原子还没太搞明白哪个原子和哪个原子相对应。
这个结构整了好几天还没整明白,请求各位老师前辈帮助,希望各位老师前辈能帮我理清一下基本的应用思路,谢谢大家



202406191548023741..png (971.95 KB, 下载次数 Times of downloads: 15)

用gauss转换的键连关系整体错误

用gauss转换的键连关系整体错误

202406191547022117..png (252.71 KB, 下载次数 Times of downloads: 12)

202406191547022117..png

Autosave of PET (1 0 0) (2).pdb

230.73 KB, 下载次数 Times of downloads: 3

Autosave of PET (1 0 0) (2).cif

213.78 KB, 下载次数 Times of downloads: 3

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

2#
发表于 Post on 2024-6-19 22:19:05 | 只看该作者 Only view this author
sobtop通常是对单个分子建立拓扑文件,而不是对整个盒子里所有分子一起创建拓扑文件。先把单个完整的链弄出来在sobtop里创建对应的[moleculetype]和相关信息,模拟整个体系时直接在[molecules]里指定聚合物链的分子数就完了(注意这要求体系里所有聚合物分子里原子顺序都是相同的,而且不同分子里原子的出现不能互相交错)
北京科音自然科学研究中心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

43

帖子

0

威望

175

eV
积分
218

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-6-20 09:05:47 | 只看该作者 Only view this author
sobereva 发表于 2024-6-19 22:19
sobtop通常是对单个分子建立拓扑文件,而不是对整个盒子里所有分子一起创建拓扑文件。先把单个完整的链弄出 ...

老师,一个链的构建是对自己构建的晶面体系转换成PDB格式然后只取一条链并在PDB文件中加入晶胞信息,在像Sobtop教程里PDB文件生成相应的gro和top文件那样做吗,因为我是根据晶胞做的晶面,不同层间链的原子顺序还是有差别的,如果只是按照取一条链来计算,能反应实际的结构吗。
还有一个困扰我的问题是我MS生成的PDB格式键连关系有些乱,其实我有些搞不清在VMD或者Gauss打开后这种混乱的键连关系对我后续做模拟有影响吗,因为我发现不同软件打开后,键连关系都不太一样,可能是判断标准不一样,有人说可以调整判断的标准,但是在这个软件中调完,在其他软件中会适用吗,还是根据调整后的键连关系直接用于Gromacs模拟。这些混乱连键的关系如果有影响我是需要一个个在Gauss中作校正吗,然后在生成单链的top、itp等文件吗。
后续生成聚合物体系是需要用Packmol构建吗,对链加以限制,使其能形成片层结构,我到时候在研究一下怎么加限制才能实现这个建模。
老师我还有一个问题是我自己画的PET三聚体结构,做了能量优化后计算RESP电荷,我发现我自己画的这个结构的原子顺序和用MS构建的晶面结构各个单元的原子顺序完全不一样,我如果只做单链的话可能会好替换一些直接进行手动替换。之前对片层结构的替换说实话很难对应上各个原子电荷。
老师,我问的问题可能相当无知,之前做蛋白配体复合物模拟也是跟着教程做的我对很多基础性的东西确实不太了解,真的麻烦老师您了

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

4#
发表于 Post on 2024-6-21 06:03:24 | 只看该作者 Only view this author
雨中宇 发表于 2024-6-20 09:05
老师,一个链的构建是对自己构建的晶面体系转换成PDB格式然后只取一条链并在PDB文件中加入晶胞信息,在像 ...

一条链就是一个分子,没晶胞信息的事

每个分子里原子顺序不同就必须想办法调整成相同

没有叫Gauss的程序,Gaussian还是GaussView分清楚

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

2

帖子

0

威望

67

eV
积分
69

Level 2 能力者

5#
发表于 Post on 2025-7-25 10:49:58 | 只看该作者 Only view this author
请问您后面怎么解决的啊?

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

6#
发表于 Post on 2025-7-25 11:49:42 | 只看该作者 Only view this author
本帖最后由 牧生 于 2025-7-25 11:53 编辑

最简单,对电荷要求不高的方法:

MS画出的聚合物结构式,加氢饱和后,先简单优化一下,保存。用gv打开,分子链的首尾端各去掉一个氢原子,然后把首尾段那两个碳原子用形式化学键连接起来,存为mol2格式。

用sobtop得到拓扑文件,电荷就直接调用open babel来得到MMFF94电荷。



加上periodic-molecules = yes
又菜又爱玩

2

帖子

0

威望

67

eV
积分
69

Level 2 能力者

7#
发表于 Post on 2025-7-25 12:56:53 | 只看该作者 Only view this author
牧生 发表于 2025-7-25 11:49
最简单,对电荷要求不高的方法:

MS画出的聚合物结构式,加氢饱和后,先简单优化一下,保存。用gv打开, ...

请问可以用PET晶体的cif文件在MS构建PET晶区的结构,然后用gromacs做分子动力学模拟吗?还有别的构建PET表面的方法吗?

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 03:52 , Processed in 0.162428 second(s), 23 queries , Gzip On.

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