过渡金属氧化物团簇SCF不收敛求助

muyueyan · 发表于 Post on 2023-10-23 10:24:39

各位老师好，本人在利用Gaussian09进行CuO团簇的限制性优化时，算了四天也没算完。考虑到可能是金属原子太多（共96个原子，冻结），就先优化半层（共16个原子）以检验方法的正确性，但遇到了SCF不收敛（L502出错）的问题，精读了sob老师的解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)这篇博文，尝试了PBE0、TPSSH等泛函，Cu使用SDD，O使用6-31G（d）基组并使用scf=conver=6，vshift=350，int=ultrafine、opt=gdiss等方法后，依旧难以收敛。这些泛函和基组是最开始用8个CuO原子成功收敛后才使用的，但是结构非常奇怪。请问是否是我的初始结构有问题呢？请各位老师闲暇之余指点一下，谢谢各位老师！

以下是我的输入文件：

CuO111.gjf (6.12 KB, 下载次数 Times of downloads: 3)

CuO3-4.gjf (1.23 KB, 下载次数 Times of downloads: 2)

muyueyan · 发表于 Post on 2023-10-24 08:15:28

我明白了，谢谢sob老师和zjxitcc老师百忙之中的指点！！！

sobereva · 发表于 Post on 2023-10-24 04:33:41

如果用当前的簇模型是为了近似表现固体，这么算纯属瞎折腾。SCF总是很难收敛，而且还不能不考虑边界效应（需要加背景电荷、capped ECP，很折腾），而且还没法优化晶胞参数。

CP2K优化这种晶体又快又容易又不花钱

muyueyan · 发表于 Post on 2023-10-23 15:18:00

好的，因为我们只有Gaussian的版权，所以才没考虑其他，那我用CP2K试试，非常感谢老师的帮助

zjxitcc · 发表于 Post on 2023-10-23 13:40:59

不用找报错、不用解决报错、不要试图解决SCF不收敛、不要再用孤立体系量化程序做这个体系计算，应立即改用周期性计算程序（例如QE, VASP, CP2K等）。

muyueyan · 发表于 Post on 2023-10-23 10:29:22

这是我8个CuO原子优化后的结构，可能是两个氧之间的距离太近了

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

[新手求助] 过渡金属氧化物团簇SCF不收敛求助