计算化学公社

 找回密码 Forget password
 注册 Register

二氧化硅与水、钠离子的弱相互作用优化收敛问题

查看数: 18376 | 评论数: 17 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-2-28 12:29

正文摘要:

各位老师们好,初学gaussian,关于二氧化硅与水、钠离子的弱相互作用优化收敛问题想请假一下。我用了二氧化硅团簇、一个钠离子、一个水分子(MS导入为mol文件),总共55各原子。在gaussian上设置相关参数后,用的40 ...

回复 Reply

北大-陶豫 发表于 Post on 2021-3-2 12:30:20
sobereva 发表于 2021-3-1 16:00
PM6对这个体系不怎么靠谱,而且还是没带色散校正的情况

3-21G就太糟糕了,因为还达不到能定性正确描述 ...

嗷,是我无知了,谢谢指正
zzzco 发表于 Post on 2021-3-1 22:11:24
wzkchem5 发表于 2021-3-1 21:52
这种体系不建议用clean,clean只对比较小的、常见的分子可靠,对较大或者不常见的分子,可能clean后的结 ...

好的 谢谢 我再试试 感谢化学公社每一位热心解答的老师
zzzco 发表于 Post on 2021-3-1 22:09:46
wzkchem5 发表于 2021-3-1 21:52
这种体系不建议用clean,clean只对比较小的、常见的分子可靠,对较大或者不常见的分子,可能clean后的结 ...

好的 谢谢 我再试试 感谢化学公社每一位解答的老师
wzkchem5 发表于 Post on 2021-3-1 21:52:38
zzzco 发表于 2021-3-1 18:03
sob老师,我的课题是二氧化硅与阳离子(盐溶液)之间的相互作用过程,用了分子模拟后打算继续用量化加点 ...

这种体系不建议用clean,clean只对比较小的、常见的分子可靠,对较大或者不常见的分子,可能clean后的结构比clean之前还不合理。
建议在MS里用UFF优化结构,然后导入GaussView,不clean直接算。如果初始团簇足够合理,连UFF优化都可以省。
zzzco 发表于 Post on 2021-3-1 18:03:38
本帖最后由 zzzco 于 2021-3-1 18:26 编辑
sobereva 发表于 2021-3-1 16:01
已经说了,当前模型根本没意义,说耗时多少也没意义

顺带一提,恰当模型下,可以把底层(不与Na和水直 ...

sob老师,我的课题是二氧化硅与阳离子(盐溶液)之间的相互作用过程,用了分子模拟后打算继续用量化加点深度,然后结合宏观实验做一个讨论。量化模型的合理性和可靠性是我一直在努力的,我是用的MS软件中自带的无定形二氧化硅构建的4.5埃的团簇,本来是类似于圆形的,先MS的clean,保存成mol文件,再Gaussian上用clean,就变成了不是圆形状的二氧化硅团簇了。请问老师有什么好的建议吗
sobereva 发表于 Post on 2021-3-1 16:01:55
zzzco 发表于 2021-3-1 14:37
好的,sob老师,昨天的收敛结果我看错了,是电脑出bug了,看了fchk文件后,还没计算成功,而且计算freq本 ...

已经说了,当前模型根本没意义,说耗时多少也没意义

顺带一提,恰当模型下,可以把底层(不与Na和水直接接触的)二氧化硅的原子用相对小一些的6-31G*
sobereva 发表于 Post on 2021-3-1 16:00:10
北大-陶豫 发表于 2021-3-1 14:51
不要一上来就 B3LYP/6-311G** opt freq,你应该先预优化一下,比如先 opt pm6,然后 opt B3LYP/3-21G,逐 ...

PM6对这个体系不怎么靠谱,而且还是没带色散校正的情况

3-21G就太糟糕了,因为还达不到能定性正确描述体系的程度,拿它预优化比起不做预优化,最终总耗时甚至更高
zzzco 发表于 Post on 2021-3-1 15:28:28
北大-陶豫 发表于 2021-3-1 14:51
不要一上来就 B3LYP/6-311G** opt freq,你应该先预优化一下,比如先 opt pm6,然后 opt B3LYP/3-21G,逐 ...

明白了 我现在改 谢谢指导!!
北大-陶豫 发表于 Post on 2021-3-1 14:51:26
zzzco 发表于 2021-3-1 14:37
好的,sob老师,昨天的收敛结果我看错了,是电脑出bug了,看了fchk文件后,还没计算成功,而且计算freq本 ...

不要一上来就 B3LYP/6-311G** opt freq,你应该先预优化一下,比如先 opt pm6,然后 opt B3LYP/3-21G,逐步加大基组,这样可以节约大量时间。
zzzco 发表于 Post on 2021-3-1 14:37:30
sobereva 发表于 2021-2-28 20:11
你都说是钠离子了,显然整体应当带一个正电荷

好的,sob老师,昨天的收敛结果我看错了,是电脑出bug了,看了fchk文件后,还没计算成功,而且计算freq本身也比较耗时,白高兴了。。。另外请教一下sob老师,这个体系计算大概需要多久,假设构型没出错,也按照老师给的关键词写输入文件(如果我能早点知道不收敛就可以早点改,像之前算了两三天等不来)
sobereva 发表于 Post on 2021-2-28 20:11:17
zzzco 发表于 2021-2-28 19:20
您好,是这样的,我弄了两个文件,一个是有电荷,一个是无电荷。我是看了sob老师的文章才知道怎么加电荷

你都说是钠离子了,显然整体应当带一个正电荷
zzzco 发表于 Post on 2021-2-28 19:20:36
北大-陶豫 发表于 2021-2-28 14:50
楼上说得很好了,我再补充一点。你为什么是无电荷 doublet,不应该是 +1 电荷 singlet 吗?doublet 的话就 ...

您好,是这样的,我弄了两个文件,一个是有电荷,一个是无电荷。我是看了sob老师的文章才知道怎么加电荷
zzzco 发表于 Post on 2021-2-28 17:23:29
sobereva 发表于 2021-2-28 14:14
有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员 ...

按照freeman老师和sob老师的修改建议,我已改正,并且收敛了。非常感谢老师们。sob老师让看的文章我之前都有看过,是学习不到位,我会继续努力。临近毕业,我需要用gaussian算很少的一部分数据,很多操作上确实不熟练,见笑了。还有,sob老师,我报了您这届的班,但要写论文,所以去不了,我想单独请教我体系的问题,不知私聊您可好?
北大-陶豫 发表于 Post on 2021-2-28 14:50:29
楼上说得很好了,我再补充一点。你为什么是无电荷 doublet,不应该是 +1 电荷 singlet 吗?doublet 的话就不得不用 ub3lyp,比 rb3lyp 又慢一大截

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

GMT+8, 2026-2-25 21:09 , Processed in 0.184468 second(s), 25 queries , Gzip On.

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