计算化学公社

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

[GROMACS] 磷酸氢根离子采取RESP2电荷、用sobtop生成拓扑文件后,能量最小化出错

[复制链接 Copy URL]

6

帖子

0

威望

105

eV
积分
111

Level 2 能力者

本帖最后由 照夜清 于 2023-7-4 21:10 编辑

我的需求是实现HPO4 2-离子在水溶液环境下模拟


前期结构拓扑文件准备:
1.首先通过Gaussview画出HPO4 2-的结构,得到HPO42-.mol2文件。
2.然后运行 ./RESP2.sh HPO42-.mol2 -2 1,得到HPO42-.chg。
3.在Sobtop中依次载入HPO42-的mol2、chg文件,使用GAFF原子类型,分别产生HPO42-的gro、itp、top文件。

gromacs模拟,尝试了如下方案:
1.先是用gmx insert-molecules -box 12 12 12 -ci HPO42-.gro -o  test.gro -nmol 54 产生盒子

然后直接用em.mdp开始能量最小化(见文件test.gro HPO42-.itp test.top em.mdp)
程序很快进行不下去

2.后来又尝试在盒子中添加54个CA离子平衡,同样失败;

3.在盒子中加满水,同样失败

发现对少量(20个左右)HPO42-离子能量最小化可以成功,感觉是拓扑文件的问题,但不知道问题具体在哪,请各位大佬、老师指教!


ps:也尝试将优化后的chg中bonds键长改写进itp文件中,同样出错

202307021852241544..png (129.46 KB, 下载次数 Times of downloads: 15)

54个纯HPO42-离子能量最小化

54个纯HPO42-离子能量最小化

HPO42-.mol2

443 Bytes, 下载次数 Times of downloads: 8

Gaussview画出

HPO42-.chg

354 Bytes, 下载次数 Times of downloads: 1

RESP2.sh脚本产生

RESP2.sh

3.2 KB, 下载次数 Times of downloads: 4

HPO42-.itp

3.14 KB, 下载次数 Times of downloads: 10

HPO42-.gro

384 Bytes, 下载次数 Times of downloads: 5

HPO42-.top

342 Bytes, 下载次数 Times of downloads: 4

em.mdp

390 Bytes, 下载次数 Times of downloads: 4

test.top

719 Bytes, 下载次数 Times of downloads: 5

test.gro

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

1376

帖子

0

威望

3984

eV
积分
5360

Level 6 (一方通行)

2#
发表于 Post on 2023-7-2 21:25:09 | 只看该作者 Only view this author
本帖最后由 牧生 于 2023-7-3 07:41 编辑

我使用你的结构式,同样的提示

Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 500 (which may not be possible for your system).
It stopped because the algorithm tried to make a new step whose size was too
small, or there was no change in the energy since last step. Either way, we
regard the minimization as converged to within the available machine
precision, given your starting configuration and EM parameters.


虽然我查到的磷酸氢根的结构式也是这样的

   但是试了好多次,都无法正常跑MD


但是,如果我把氢原子画在和磷原子相连,那么则一切都非常正常

又菜又爱玩

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

3#
发表于 Post on 2023-7-3 01:02:03 | 只看该作者 Only view this author
照夜清 发表于 2023-7-2 19:01
实在是尝试了多种方法,三天就耗在这一个小离子上,求助!!!

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

6

帖子

0

威望

105

eV
积分
111

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2023-7-4 21:13:14 | 只看该作者 Only view this author
牧生 发表于 2023-7-2 21:25
我使用你的结构式,同样的提示

Energy minimization has stopped, but the forces have not converged t ...

感觉可能是GAFF力场参数的问题,不能直接通过指认原子类型来生成拓扑文件。

6

帖子

0

威望

105

eV
积分
111

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2023-7-4 21:14:12 | 只看该作者 Only view this author
sobereva 发表于 2023-7-3 01:02
有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员 ...

好的,学生看的不够认真,已删除并重新编辑!

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

6#
发表于 Post on 2023-7-5 02:37:01 | 只看该作者 Only view this author
照夜清 发表于 2023-7-4 21:13
感觉可能是GAFF力场参数的问题,不能直接通过指认原子类型来生成拓扑文件。

不是。关键是你给sobtop载入的mol2文件里连接关系必须正确
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

6

帖子

0

威望

105

eV
积分
111

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2023-7-8 16:22:53 | 只看该作者 Only view this author
sobereva 发表于 2023-7-5 02:37
不是。关键是你给sobtop载入的mol2文件里连接关系必须正确

我再次检查了sobtop载入磷酸氢根的mol2文件中的连接关系,应该是正确的

L[N5Y2QVCS~4G}W)D)Z]Q3V.png (15.9 KB, 下载次数 Times of downloads: 12)

L[N5Y2QVCS~4G}W)D)Z]Q3V.png

265

帖子

0

威望

1209

eV
积分
1474

Level 4 (黑子)

8#
发表于 Post on 2023-7-9 19:03:35 | 只看该作者 Only view this author
照夜清 发表于 2023-7-8 16:22
我再次检查了sobtop载入磷酸氢根的mol2文件中的连接关系,应该是正确的

你把1-4项pairs删除了试试。

63

帖子

0

威望

1087

eV
积分
1150

Level 4 (黑子)

9#
发表于 Post on 2023-7-31 23:52:40 | 只看该作者 Only view this author
slxc920113 发表于 2023-7-9 19:03
你把1-4项pairs删除了试试。

大佬你好,我也有遇到类似的情况。我把1-4的pair删去。发现浓度一高,出现了聚集。不知道这个是否合理呢?

4

帖子

0

威望

97

eV
积分
101

Level 2 能力者

10#
发表于 Post on 2024-6-24 15:22:09 | 只看该作者 Only view this author
请问这个问题最后有比较好的解决办法了吗

38

帖子

0

威望

1159

eV
积分
1197

Level 4 (黑子)

11#
发表于 Post on 2024-6-25 00:23:21 | 只看该作者 Only view this author
本帖最后由 yxdd98 于 2024-6-25 00:51 编辑
xxx1025 发表于 2024-6-24 15:22
请问这个问题最后有比较好的解决办法了吗

我最近做了一个硫酸分子的模拟,也是分子数多了跑能量最小化显示力很大,达不到收敛限,分子数少了虽然能量最小化能跑通,但一跑动力学就崩溃。
说一下我的解决方法,在用sobtop时,选择键长和键角是基于Hessian矩阵获得,但二面角等是用力场预定义的 (好像是对应选项7?),之后就能正常跑起来了,而且观察了下轨迹,硫酸整体四面体维持较好,两个OH也能自由转动。
或者也可以让全部力常数都基于Hessian矩阵获得,但这样就刚性过强,那些氢就不能自由转了

1632

帖子

4

威望

4088

eV
积分
5800

Level 6 (一方通行)

喵星人

12#
发表于 Post on 2024-6-25 01:31:09 | 只看该作者 Only view this author
首先这种简单离子的Amber力场参数你搜一下文献就有了,而且专门拟合的不需要自己造轮子

4

帖子

0

威望

97

eV
积分
101

Level 2 能力者

13#
发表于 Post on 2024-6-27 10:34:54 | 只看该作者 Only view this author
yxdd98 发表于 2024-6-25 00:23
我最近做了一个硫酸分子的模拟,也是分子数多了跑能量最小化显示力很大,达不到收敛限,分子数少了虽然能 ...

感谢!我试试~

2

帖子

0

威望

13

eV
积分
15

Level 1 能力者

14#
发表于 Post on 2024-11-11 16:58:13 | 只看该作者 Only view this author
yxdd98 发表于 2024-6-25 00:23
我最近做了一个硫酸分子的模拟,也是分子数多了跑能量最小化显示力很大,达不到收敛限,分子数少了虽然能 ...

你好,可以发我一份硫酸根的拓扑文件吗

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 05:16 , Processed in 0.176516 second(s), 24 queries , Gzip On.

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