高斯做Bi掺杂二氧化硅团簇的结构优化收敛问题求助

尛瞐尛 · 发表于 Post on 2022-11-4 16:43:00

老师们，我用高斯16正在做Bi掺杂二氧化硅团簇的结构优化，遇到了不收敛问题，看了sob老师的教程后，加了SCF(novaracc,noincfock) int=acc2e=12 SCF=maxcyc=500等关键词后仍然不收敛，我猜测是结构问题，但是不知道问题在哪，能请老师帮我看看输入文件吗？感激不尽！

ABetaCarw · 发表于 Post on 2022-11-4 16:59:25

弱弱问一句这个模型怎么来的？Bi为啥孤零零在那里挂着，你想优化一个什么结构才这么建模的？

wzkchem5 · 发表于 Post on 2022-11-4 17:02:28

实验上已知Bi是掺杂在那个地方吗？实验上已知Bi的价态吗？

尛瞐尛 · 发表于 Post on 2022-11-4 17:04:54

ABetaCarw 发表于 2022-11-4 16:59
弱弱问一句这个模型怎么来的？Bi为啥孤零零在那里挂着，你想优化一个什么结构才这么建模的？

我是跟着这篇论文做的，他做的是铋的间隙掺杂

尛瞐尛 · 发表于 Post on 2022-11-4 17:07:09

wzkchem5 发表于 2022-11-4 17:02
实验上已知Bi是掺杂在那个地方吗？实验上已知Bi的价态吗？

实验上铋的价态有三价，二价，零价都有

尛瞐尛 · 发表于 Post on 2022-11-4 17:07:48

尛瞐尛发表于 2022-11-4 17:07
实验上铋的价态有三价，二价，零价都有

至于位置我是根据计算形成能确定

chands · 发表于 Post on 2022-11-4 17:17:14

你把Bi去掉后剩下的结构收敛吗？

wzkchem5 · 发表于 Post on 2022-11-4 17:24:10

又仔细看了一下结构，41号氢和55号氢太近了，此外Bi和Si有可能也有点太近了。
其实建议用GFN2-xTB先做一个构象搜索，或者简单地退一下火也可以，可以让所有羟基的取向以及Bi的位置都是合理的。另外这个团簇有点太小了，没有把Bi完全包起来，可能不能很好地反映Bi掺杂在体相里的性质，如果大团簇DFT算不动，最起码用ONIOM(DFT:GFN2-xTB)。
用赝势粗优化以后，建议用全电子相对论方法cross-check一下，赝势对于这么重的元素未必靠谱。单点能建议用全电子相对论方法，并且加SOC。
还有一个问题（虽然和你的SCF收敛问题无关）是必须加色散校正。

尛瞐尛 · 发表于 Post on 2022-11-4 18:36:12

chands 发表于 2022-11-4 17:17
你把Bi去掉后剩下的结构收敛吗？

好的，我试试，感谢老师！

尛瞐尛 · 发表于 Post on 2022-11-4 18:36:55

wzkchem5 发表于 2022-11-4 17:24
又仔细看了一下结构，41号氢和55号氢太近了，此外Bi和Si有可能也有点太近了。
其实建议用GFN2-xTB先做一个 ...

好的，谢谢老师！

zhu359 · 发表于 Post on 2023-8-7 11:03:15

问一下这个结构为什么有氢存在啊

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 高斯做Bi掺杂二氧化硅团簇的结构优化收敛问题求助