“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

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

[GROMACS] 含二氧化硅的盐溶液冻结模拟与冰拼接后能量最小化不收敛

[复制链接 Copy URL]

13

帖子

0

威望

149

eV
积分
162

Level 3 能力者

利用其他大神帖子中提到的amber14sb_OL15-Silica.ff中力场参数,并用sobtop生成了二氧化硅的itp文件,现含二氧化硅的盐溶液能量平衡后无问题,但与冰结构拼接后能量最小化不收敛,目出现边界处出现
发散情况,确定界面处无原子重叠情况,恳请各位老师指点

冰与含二氧化硅盐溶液拼接后的结构图.png (43.11 KB, 下载次数 Times of downloads: 36)

冰与含二氧化硅盐溶液拼接后的结构图

冰与含二氧化硅盐溶液拼接后的结构图

能量最小化不收敛截图.png (17.89 KB, 下载次数 Times of downloads: 35)

能量最小化不收敛截图

能量最小化不收敛截图

能量最小化后的结构图.png (101.21 KB, 下载次数 Times of downloads: 32)

能量最小化后的结构图

能量最小化后的结构图

206

帖子

0

威望

1157

eV
积分
1363

Level 4 (黑子)

终身学习

2#
发表于 Post on 2025-5-13 14:57:11 | 只看该作者 Only view this author
看一下atom 1491是什么原子,这么大的受力可能是拓扑文件或者起始结构不合理
另外说句不重要的,盒子周围的真空区有些太大了,能量最小化完冰的结构已经明显看到散了
Open source enables open science.

1580

帖子

0

威望

5119

eV
积分
6699

Level 6 (一方通行)

3#
发表于 Post on 2025-5-13 16:10:00 | 只看该作者 Only view this author
你的二氧化硅只是孤立的一个小块,那么,有使用羟基饱和边缘吗?

又菜又爱玩

13

帖子

0

威望

149

eV
积分
162

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2025-5-13 17:26:17 | 只看该作者 Only view this author
牧生 发表于 2025-5-13 16:10
你的二氧化硅只是孤立的一个小块,那么,有使用羟基饱和边缘吗?

在ms里已经用OH饱和过了,净电荷为0,但是在vmd中显示就是这样的,之前您在另一篇帖子回复我的推荐二氧化硅力场参数已经使用了,但这次就是不收敛,请大神费心再帮忙看看

本版积分规则 Credits rule

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

GMT+8, 2026-4-13 18:31 , Processed in 0.210225 second(s), 23 queries , Gzip On.

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