计算化学公社

标题: 过渡金属氧化物团簇SCF不收敛求助 [打印本页]

作者
Author: muyueyan 时间: 2023-10-23 10:24
标题: 过渡金属氧化物团簇SCF不收敛求助
各位老师好，本人在利用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原子成功收敛后才使用的，但是结构非常奇怪。请问是否是我的初始结构有问题呢？请各位老师闲暇之余指点一下，谢谢各位老师！

以下是我的输入文件：

(, 下载次数 Times of downloads: 8)

(, 下载次数 Times of downloads: 3)
(, 下载次数 Times of downloads: 2)

作者
Author: muyueyan 时间: 2023-10-23 10:29
(, 下载次数 Times of downloads: 4)
这是我8个CuO原子优化后的结构，可能是两个氧之间的距离太近了

作者
Author: zjxitcc 时间: 2023-10-23 13:40
不用找报错、不用解决报错、不要试图解决SCF不收敛、不要再用孤立体系量化程序做这个体系计算，应立即改用周期性计算程序（例如QE, VASP, CP2K等）。

作者
Author: muyueyan 时间: 2023-10-23 15:18
好的，因为我们只有Gaussian的版权，所以才没考虑其他，那我用CP2K试试，非常感谢老师的帮助

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

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

作者
Author: muyueyan 时间: 2023-10-24 08:15
我明白了，谢谢sob老师和zjxitcc老师百忙之中的指点！！！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)