计算化学公社

标题: GAMESS-US计算过渡金属能量分解异常且SCF不收敛问题 [打印本页]

作者
Author: genglj 时间: 2020-3-21 13:10
标题: GAMESS-US计算过渡金属能量分解异常且SCF不收敛问题
各位好，请问在使用GAMESS-US计算能量分解时遇到能量分解结果异常的问题怎么处理。
我的体系是Co团簇与NH3的复合物，用高斯（BPW91/def2tzvp for Co,N,H）优化后的稳定结构作为inp文件的结构,
刚开始用的MP2，报错911；改为DFT后（BPW91/def2tzvp for Co,N,H）正常结束后能量分解结果异常并且SCF不收敛（能量震荡），
尝试更换泛函基组为 Bpw91/Lanl2TZ(f) for Co, 6-311+g** for N,H 及 B3lyp/Lanl2TZ(f) for Co, 6-311+g** for N,H , 均不收敛
也尝试了 $SCF DIRSCF=.T. SHIFT=.T. $END 及 $SCF DIIS=.TRUE. ETHRSH=2.0 SWDIIS=0.005 $END，均不收敛
请教各位大神们这种情况该如何修正使其收敛保证能量分解结果正常
谢谢
(, 下载次数 Times of downloads: 27)

作者
Author: sobereva 时间: 2020-3-21 13:21
就用def2-TZVP。本来这种过渡金属成簇的体系SCF就是难收敛，再给NH加弥散函数更难收敛。

确保GAMESS-US算单点能顺利结束之前不要尝试做能量分解。

用Multiwfn载入Gaussian之前同级别算完的fch文件，用主功能100的子功能2导出成GAMESS-US输入文件，把里面记录的初猜当做当前计算的初猜，若单点任务能顺利收敛，再去做能量分解

另外可考虑改用PBE0-D3(BJ)。纯泛函普遍不如杂化泛函容易收敛。另外就算用纯泛函，也没什么理由用BPW91，既不流行也不是什么好选择。

作者
Author: genglj 时间: 2020-3-21 13:42

sobereva 发表于 2020-3-21 13:21
就用def2-TZVP。本来这种过渡金属成簇的体系SCF就是难收敛，再给NH加弥散函数更难收敛。

确保GAMESS-US ...

谢谢sob老师，我去试试
另外，这个体系的高斯计算分析部分都用的BPW91泛函，若能量分解分析部分用不同的泛函能够被认可吗

作者
Author: sobereva 时间: 2020-3-22 15:24

genglj 发表于 2020-3-21 13:42
谢谢sob老师，我去试试
另外，这个体系的高斯计算分析部分都用的BPW91泛函，若能量分解分析部分用不同的 ...

如果BPW91你怎么弄都没法成功使用LMOEDA，那你就有充分理由在LMOEDA这部分换泛函，即便审稿人问也可以回复得有理有据。并且这点应当在文章中注明。

作者
Author: genglj 时间: 2020-3-23 11:15
好的，谢谢您的解答

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