计算化学公社

标题: Gaussian计算含过渡金属的结构优化始终报错自洽场收敛失败 [打印本页]

作者
Author: lomon 时间: 2022-10-7 15:47
标题: Gaussian计算含过渡金属的结构优化始终报错自洽场收敛失败
如题，我使用g09对体系进行结构优化，始终不能收敛，目前已经尝试过的方法有：能级提升尝试过500，800，1800；更改泛函m06l为B3LYP；opt=loose配合scf=conver=6 IOp=(7/127=-99)；opt=gdiis等方法，我按照sob老师发的scf自洽场步收敛的解决方法的帖子对命令进行修改仍然不能收敛，包括scf=maxcycle=100也试过，仍然不能收敛。随帖附上最近一次的输入文件和log文件，向各位老师请教问题出现的原因，以及解决办法

作者
Author: wzkchem5 时间: 2022-10-7 15:50
你要算的就是团簇吗，还是本来要算合金，用团簇模型来算的合金

作者
Author: lomon 时间: 2022-10-7 15:51

wzkchem5 发表于 2022-10-7 15:50
你要算的就是团簇吗，还是本来要算合金，用团簇模型来算的合金

算团簇，分析水分子在团簇分子上的吸附和解离行为

作者
Author: chands 时间: 2022-10-7 16:26
团簇在优化过程中变形很大，考虑建模是否合理。
还有你的H2O分成H原子和OH自由基？

作者
Author: lomon 时间: 2022-10-7 17:23

chands 发表于 2022-10-7 16:26
团簇在优化过程中变形很大，考虑建模是否合理。
还有你的H2O分成H原子和OH自由基？

1.之前用m06l的时候就已经发现优化后的结构变化非常大，所以想着添加赝势并更换泛函为m06l，但是更换后目前一直不收敛；这个输入文件用的团簇结构是之前单独优化完成的团簇结构，所以没考虑结构合理性的问题
2.主要目的是分析吸附和解离的问题，所以想着先用这个结构做优化，如果能走的通再推广到其他结构上去

作者
Author: sobereva 时间: 2022-10-8 19:11
犯了低级错误
要用genecp而非gen，仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

另外，opt freq目的用SDD足矣，没必要def2-TZVP那么贵的

本来过渡金属团簇类体系就SCF难收敛，干嘛把SCF迭代次数上限从默认的128减小到100，完全莫名其妙

此类体系用PBE0-D3(BJ)更好，参看
使用量子化学程序基于簇模型计算金属表面吸附问题
http://sobereva.com/540（http://bbs.keinsci.com/thread-16532-1-1.html）

作者
Author: lomon 时间: 2022-10-9 13:37

sobereva 发表于 2022-10-8 19:11
犯了低级错误
要用genecp而非gen，仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入

谢谢Sob老师的回复，我再重新试一下

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