求助gauss计算单点能出现l502报错

vvkki · 发表于 Post on 2024-5-20 12:05:22

各位老师好，我在跑完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报错，附上输入文件和输出文件

2-CIS-RE-sp-1.gjf (1.98 KB, 下载次数 Times of downloads: 10)

2-CIS-RE-sp-1.log (176.64 KB, 下载次数 Times of downloads: 1) ，请各位老师帮忙看一下，非常感谢！

wzkchem5 · 发表于 Post on 2024-5-20 17:05:41

尝试先用def2svp收敛以后，得到的波函数作为实际的def2tzvp计算的初猜

wcy鱼饵 · 发表于 Post on 2024-5-20 18:01:27

如果用DEF2TZVP收敛较慢，可以用SVP优化完再用DEF2TZVP优化，亲测有效

vvkki · 发表于 Post on 2024-5-20 20:09:10

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进行频率和振动分析，也是显示不收敛

vvkki · 发表于 Post on 2024-5-20 20:11:42

wcy鱼饵发表于 2024-5-20 18:01
如果用DEF2TZVP收敛较慢，可以用SVP优化完再用DEF2TZVP优化，亲测有效

您好，我现在def2tzvp计算单点能不能收敛，且svp优化后用def2tzvp优化也不可以，您可以说一下具体的操作步骤吗，可能是我没有领悟清楚

wzkchem5 · 发表于 Post on 2024-5-20 21:38:26

vvkki 发表于 2024-5-20 13:09
老师您好，我尝试了使用b3lyp/def3svp进行振动分析和频率优化，后使用m062x/deftzvp计算单点能，结果还是 ...

检查你现在m062x/def2tzvp迭代后的轨道是否符合化学常识

sobereva · 发表于 Post on 2024-5-21 02:49:18

def2-SVP下优化如果能收敛，最后def2-TZVP算单点的时候读取优化任务的chk里的波函数当初猜。如果还失败，上传此时的输入输出文件

vvkki · 发表于 Post on 2024-5-23 09:42:03

老师您好，我读取了chk的波函数作为初猜，依旧显示l502报错，这是输入和输出文件，请您帮忙看看

2-CIS-RE-3-sp-1.gjf (1.92 KB, 下载次数 Times of downloads: 8)

2-CIS-RE-3-sp-1.log (173.19 KB, 下载次数 Times of downloads: 0)

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

[新手求助] 求助gauss计算单点能出现l502报错