计算化学公社

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

[GROMACS] 求助:THF/H2O混合溶剂的模拟

[复制链接 Copy URL]

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

本帖最后由 lr8616 于 2026-1-9 00:32 编辑

刚开始学习GROMACS,今天试着模拟THF/H2O的混合溶剂,平衡相之后发现分子在盒子中的分布不均匀了,还想请高手指点一下。

1. 利用packmol生成 10 10 10的盒子,add_box_sides 1.5,填充650个THF和3000个H2O。
2. 利用Sobtop生成THF的itp文件,原子类型用GAFF2,根据老师的帖子将其中电荷替换为RESP2(0.5)H2O用的是opc3.itp。
3. 能量最小化:

4. 平衡相:

然后不知道为啥溶剂分子的分布就变奇怪了呢?

848

帖子

4

威望

1912

eV
积分
2840

Level 5 (御坂)

A Student

2#
发表于 Post on 2026-1-9 00:23:57 | 只看该作者 Only view this author
NPT跑长些。
敬仰一针见血的指责,厌倦别有用心的赞美。

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2026-1-9 00:39:04 | 只看该作者 Only view this author

谢谢回复,我刚才延长平衡相时间到1ns试了,最后得到的溶剂分子不知道为啥分层了呢?


6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

4#
发表于 Post on 2026-1-9 01:20:27 | 只看该作者 Only view this author
lr8616 发表于 2026-1-9 00:39
谢谢回复,我刚才延长平衡相时间到1ns试了,最后得到的溶剂分子不知道为啥分层了呢?

如果两种分子互溶性较差,必然分层
如果已知可以充分互溶,要么模拟条件设定不当(如温度),要么力场/原子电荷不够理想,需要进一步调节,或者看文献里是否有对当前分子优化的专一性的参数并尝试使用
北京科音自然科学研究中心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

94

帖子

0

威望

1005

eV
积分
1099

Level 4 (黑子)

5#
发表于 Post on 2026-1-9 01:22:36 | 只看该作者 Only view this author
试一试OPLS-AA力场,中性溶剂分子用1.2*CM5电荷,离子采用RESP电荷。

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2026-1-9 01:29:38 | 只看该作者 Only view this author
sobereva 发表于 2026-1-9 01:20
如果两种分子互溶性较差,必然分层
如果已知可以充分互溶,要么模拟条件设定不当(如温度),要么力场/ ...

四氢呋喃和水是可以互溶的,我再查查文献,谢谢卢老师

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2026-1-9 01:30:18 | 只看该作者 Only view this author
Dreamer1 发表于 2026-1-9 01:22
试一试OPLS-AA力场,中性溶剂分子用1.2*CM5电荷,离子采用RESP电荷。

我来试试看,感谢!

93

帖子

0

威望

453

eV
积分
546

Level 4 (黑子)

开心小鲨鱼 ORCA的饲养员/食物

8#
发表于 Post on 2026-1-10 00:34:19 | 只看该作者 Only view this author
我查了,THF 和水是部分互溶。从结构看,环烷烃部分就是疏水,虽然有氧原子亲水。如果你想模拟出互溶,可以考虑溶剂分子个数是不是需要调整,参考此贴3楼。有文献的话可以参考文献的数目。我个人感觉,多放水,少点THF,可以弄成互溶效果
如何实现两相分离的模拟结果? - 分子模拟 (Molecular Modeling) - 计算化学公社  http://bbs.keinsci.com/forum.php?mod=viewthread&tid=57520

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2026-1-10 09:05:00 | 只看该作者 Only view this author
SharkYYX2025 发表于 2026-1-10 00:34
我查了,THF 和水是部分互溶。从结构看,环烷烃部分就是疏水,虽然有氧原子亲水。如果你想模拟出互溶,可以 ...

谢谢你的回复!四氢呋喃和水是完全互溶的体系,可以任意比例混溶,主要原因是可以与水形成较强的氢键。因为我想模拟THF和水不同体积比混合的情况,所以先试了1:1。

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2026-1-10 13:22:39 | 只看该作者 Only view this author
sobereva 发表于 2026-1-9 01:20
如果两种分子互溶性较差,必然分层
如果已知可以充分互溶,要么模拟条件设定不当(如温度),要么力场/ ...

卢老师好!有两个问题还想跟您请教一下。

1. 我的初衷是模拟合成的有机分子在四氢呋喃/水 THF/H2O混合溶剂中的动力学行为,有机分子的结构很类似您674博文中的OPP分子结构(环状,有少部分刚性连接和大部分柔性结构,原子数110个),所以本来想根据您的帖子用GAFF力场来模拟,但是失败在了第一步,混合溶剂相平衡时发生了分层(图见1楼和3楼),我猜想是不是THF的电荷小了,不能很好的描述和水之间的氢键作用?如果想用GAFF力场的话不知道还有没有什么尝试的办法?


2. 我后来查了类似方向的文献,发现文献中THF分子用CGenFF力场,水分子用TIP3P模型来模拟两者不同比例混溶。想尝试用这个方法,但是我合成的有机分子在CGenFF网上生成的拓扑文件中芳环连接处的二面角罚分很高,在文献中没搜到报道的参数,这个二面角的位置类似OPP分子标红的部分(很抱歉合成的结构不方便放出来,所以用您的图来类比) 。想请教一下柔性扫描拟合二面角参数的话左边的碳纳米环部分该怎么简化呢?课题组太穷计算资源很有限,不知道这种大分子柔性扫描能算的动吗?



6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

11#
发表于 Post on 2026-1-10 22:37:49 | 只看该作者 Only view this author
lr8616 发表于 2026-1-10 13:22
卢老师好!有两个问题还想跟您请教一下。

1. 我的初衷是模拟合成的有机分子在四氢呋喃/水 THF/H2O混合 ...

1 你之前只说了用RESP2(0.5),但没说在计算时隐式溶剂模型描述的什么溶剂,因此当前电荷是否合理没法判断。最好用THF和水按体积比自定义eps的溶剂环境计算RESP2(0.5)中涉及的溶剂下的RESP电荷。从对范德华相互作用描述来说,我不认为有比GAFF或GAFF2能更好的选择。

2 根本没必要尝试CGenFF。GAFF/GAFF2描述不了的用这个也白搭,这个能描述的用GAFF/GAFF2一样能描述
北京科音自然科学研究中心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

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2026-1-11 19:31:11 | 只看该作者 Only view this author
sobereva 发表于 2026-1-10 22:37
1 你之前只说了用RESP2(0.5),但没说在计算时隐式溶剂模型描述的什么溶剂,因此当前电荷是否合理没法判断 ...

卢老师,最开始计算THF的RESP2(0.5)是用水溶剂,现在按照您的建议把介电常数改成了混合溶剂的再次计算、替换电荷,用GAFF2力场能量最小化、平衡2ns后还是分层了…

我又试了用CGenFF力场,同样平衡后THF和水1:1的体系就混溶了( •̥́ ˍ •̀ू )

根据您的教程判断我要模拟的分子GAFF力场最合适最方便,就是不知道混合溶剂模拟为啥总是不对…

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

13#
发表于 Post on 2026-1-12 04:23:33 | 只看该作者 Only view this author
lr8616 发表于 2026-1-11 19:31
卢老师,最开始计算THF的RESP2(0.5)是用水溶剂,现在按照您的建议把介电常数改成了混合溶剂的再次计算、 ...

你对CGenFF用的什么原子电荷?
北京科音自然科学研究中心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

76

帖子

0

威望

839

eV
积分
915

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2026-1-12 08:10:50 | 只看该作者 Only view this author
sobereva 发表于 2026-1-12 04:23
你对CGenFF用的什么原子电荷?

用的是上传到CGenFF官网上生成top文件里自带的电荷,我上传到附件里了。同时还上传了用Sobtop生成的THF拓扑文件,里边的电荷用自定义eps计算的RESP2(0.5)替换了。

THF_gmx_CGenFF.top

4.6 KB, 下载次数 Times of downloads: 1

THF_RESP2(0.5).itp

10.2 KB, 下载次数 Times of downloads: 1

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

15#
发表于 Post on 2026-1-13 01:11:25 | 只看该作者 Only view this author
lr8616 发表于 2026-1-12 08:10
用的是上传到CGenFF官网上生成top文件里自带的电荷,我上传到附件里了。同时还上传了用Sobtop生成的THF拓 ...

可以尝试用结合混合溶剂算的ADCH电荷,对静电势有较好重现能力,原理也理想。
也可以尝试MK电荷(拟合静电势电荷之一)。THF是比较小的刚性体系,RESP比MK没有额外好处,而由于加了额外限制,可能反倒恶化结果。
ADCH和MK在Multiwfn里都有。
北京科音自然科学研究中心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, 2026-1-23 13:10 , Processed in 0.291028 second(s), 23 queries , Gzip On.

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