计算化学公社

标题: 求助gauss计算单点能出现l502报错 [打印本页]

作者
Author:
vvkki    时间: 2024-5-20 12:05
标题: 求助gauss计算单点能出现l502报错
各位老师好,我在跑完IRC对两边的产物进行频率优化和单点能计算,频率优化使用b3lyp/6-31g* em=gd3bj,单点能使用m062x/def2tzvp em=gd3。单点能计算出现l502报错,原因是SCF不收敛,学习sob老师的帖子解决SCF不收敛问题的方法 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com),先后尝试添加关键词int=superfine、scf(nodiis,fermi)、scf(novaracc,noincfock)、scf(novaracc,noincfock)、scf(novaracc,conver=6,nodiis,vshift=400)、scf(novaracc,conver=6) iop=(7/127=-99)、逐步提高优化精度b3lyp/def2tzvp优化后m062x/def2tzvp em=gd3计算单点能,均显示l502报错,附上输入文件和输出文件 (, 下载次数 Times of downloads: 10) (, 下载次数 Times of downloads: 1) ,请各位老师帮忙看一下,非常感谢!

作者
Author:
wzkchem5    时间: 2024-5-20 17:05
尝试先用def2svp收敛以后,得到的波函数作为实际的def2tzvp计算的初猜
作者
Author:
wcy鱼饵    时间: 2024-5-20 18:01
如果用DEF2TZVP收敛较慢,可以用SVP优化完再用DEF2TZVP优化,亲测有效
作者
Author:
vvkki    时间: 2024-5-20 20:09
wzkchem5 发表于 2024-5-20 17:05
尝试先用def2svp收敛以后,得到的波函数作为实际的def2tzvp计算的初猜

老师您好,我尝试了使用b3lyp/def3svp进行振动分析和频率优化,后使用m062x/deftzvp计算单点能,结果还是显示未收敛,输出文件显示
SCF Done:  E(UM062X) =  -431.313137799     A.U. after  129 cycles
            NFock=128  Conv=0.87D-02     -V/T= 2.0062
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7520 S= 0.5010
<L.S>= 0.000000000000E+00
KE= 4.286607791839D+02 PE=-2.371639671103D+03 EE= 8.204693990312D+02
Annihilation of the first spin contaminant:
S**2 before annihilation     0.7520,   after     0.7500
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/software/g16/l502.exe at Mon May 20 19:46:09 2024.
除此之外还尝试了使用m062x/def2svp计算单点能,显示不收敛;m062x/def2svp进行频率和振动分析,也是显示不收敛
作者
Author:
vvkki    时间: 2024-5-20 20:11
wcy鱼饵 发表于 2024-5-20 18:01
如果用DEF2TZVP收敛较慢,可以用SVP优化完再用DEF2TZVP优化,亲测有效

您好,我现在def2tzvp计算单点能不能收敛,且svp优化后用def2tzvp优化也不可以,您可以说一下具体的操作步骤吗,可能是我没有领悟清楚
作者
Author:
wzkchem5    时间: 2024-5-20 21:38
vvkki 发表于 2024-5-20 13:09
老师您好,我尝试了使用b3lyp/def3svp进行振动分析和频率优化,后使用m062x/deftzvp计算单点能,结果还是 ...

检查你现在m062x/def2tzvp迭代后的轨道是否符合化学常识
作者
Author:
sobereva    时间: 2024-5-21 02:49
def2-SVP下优化如果能收敛,最后def2-TZVP算单点的时候读取优化任务的chk里的波函数当初猜。如果还失败,上传此时的输入输出文件
作者
Author:
vvkki    时间: 2024-5-23 09:42
老师您好,我读取了chk的波函数作为初猜,依旧显示l502报错,这是输入和输出文件,请您帮忙看看 (, 下载次数 Times of downloads: 8) (, 下载次数 Times of downloads: 0)




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3