计算化学公社

标题: 请教老师，高斯优化分子报错l502 [打印本页]

作者
Author: 不明白 时间: 2021-4-23 15:20
标题: 请教老师，高斯优化分子报错l502
请问是我的结构建的不合适吗，我用小基组能正常结束，换6-31g*就会不收敛，各位老师有好的建议吗Cycle 127  Pass 1  IDiag  1:
RMSU=  2.21D-07 CP:  1.00D+00 -1.02D-02  2.00D-01  1.69D+00  5.51D-03
                  CP: -3.97D-02
E= -2479.61107975939    Delta-E=       0.000164204030 Rises=F Damp=F
DIIS: error= 1.11D-04 at cycle 127 NSaved=  20.
NSaved=20 IEnMin= 4 EnMin= -2479.61159006735    IErMin= 6 ErrMin= 2.91D-05
ErrMax= 1.11D-04  0.00D+00 EMaxC= 1.00D+00 BMatC= 1.74D-04 BMatP= 6.16D-06
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.213D-01-0.234D+00-0.472D+00 0.734D+00-0.316D+00-0.280D+00
Coeff-Com:  0.365D+00-0.640D+00 0.653D+00 0.303D+00 0.494D+00-0.337D+01
Coeff-Com:  0.251D+01 0.855D+00 0.161D+01-0.162D+01 0.438D+00 0.223D+00
Coeff-Com:  0.457D+00-0.695D+00
Coeff:    -0.213D-01-0.234D+00-0.472D+00 0.734D+00-0.316D+00-0.280D+00
Coeff:    0.365D+00-0.640D+00 0.653D+00 0.303D+00 0.494D+00-0.337D+01
Coeff:    0.251D+01 0.855D+00 0.161D+01-0.162D+01 0.438D+00 0.223D+00
Coeff:    0.457D+00-0.695D+00
Gap=    0.026 Goal= None Shift= 0.000
Gap=    0.025 Goal= None Shift= 0.000
RMSDP=9.10D-03 MaxDP=6.54D-01 DE= 1.64D-04 OVMax= 8.99D-02

Cycle 128  Pass 1  IDiag  1:
RMSU=  2.32D-07 CP:  9.99D-01  1.02D+00 -9.52D-02 -2.68D+00  1.44D+00
                  CP: -3.41D-01 -2.44D-01
E= -2479.61123391736    Delta-E=    -0.000154157973 Rises=F Damp=F
DIIS: error= 8.73D-05 at cycle 128 NSaved=  20.
NSaved=20 IEnMin= 4 EnMin= -2479.61159006735    IErMin= 6 ErrMin= 2.91D-05
ErrMax= 8.73D-05  0.00D+00 EMaxC= 1.00D+00 BMatC= 1.03D-04 BMatP= 6.16D-06
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.171D-01-0.204D+00-0.394D+00 0.617D+00-0.273D+00-0.225D+00
Coeff-Com:  0.283D+00-0.563D+00 0.599D+00 0.264D+00 0.487D+00-0.330D+01
Coeff-Com:  0.250D+01 0.748D+00 0.151D+01-0.152D+01 0.415D+00 0.289D+00
Coeff-Com: -0.756D+00 0.534D+00
Coeff:    -0.171D-01-0.204D+00-0.394D+00 0.617D+00-0.273D+00-0.225D+00
Coeff:    0.283D+00-0.563D+00 0.599D+00 0.264D+00 0.487D+00-0.330D+01
Coeff:    0.250D+01 0.748D+00 0.151D+01-0.152D+01 0.415D+00 0.289D+00
Coeff:    -0.756D+00 0.534D+00
Gap=    0.027 Goal= None Shift= 0.000
Gap=    0.026 Goal= None Shift= 0.000
RMSDP=8.31D-04 MaxDP=6.87D-02 DE=-1.54D-04 OVMax= 7.99D-03

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -2479.61123392    A.U. after  129 cycles
         NFock=128  Conv=0.83D-03    -V/T= 2.4408
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.0000 <S**2>= 0.0007 S= 0.0007
<L.S>= 0.000000000000E+00
KE= 1.720943451703D+03 PE=-1.875225745941D+04 EE= 7.692643299948D+03
Annihilation of the first spin contaminant:
S**2 before annihilation    0.0007, after    0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /home-ys/Soft/Gaussian.D01.Linda/g09/l502.exe at Thu Apr 22 03:35:17 2021.
Job cpu time:    3 days  8 hours 32 minutes  0.6 seconds.
File lengths (MBytes):  RWF= 354 Int=    0 D2E=    0 Chk=    30 Scr=    1

作者
Author: hdhxx123 时间: 2021-4-23 15:41
本帖最后由 hdhxx123 于 2021-4-23 17:24 编辑

不知道你比6-31g*的小基组是什么基组，建议你先看这篇博文
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）
收敛问题看这篇博文
解决 SCF 不收敛问题的方法
http://sobereva.com/61
不会搭建初猜用gfn2xtb先优化一下
看博文将 Gaussian 与 Grimme 的 xtb 程序联用搜索过渡态、产生 IRC、做振动分析
http://sobereva.com/421（http://bbs.keinsci.com/thread-10106-1-1.html)基组的写法看这篇博文：

详解 Gaussian 中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60。最后，gjf输入文件最后多打些空行总没错

作者
Author: snljty 时间: 2021-4-23 15:42
本帖最后由 snljty 于 2021-4-23 15:45 编辑

你赝势都没读进去。先把gen换成genecp，不然计算没有意义。先别做优化，先把一个单点收敛了再说。不要说L502，直接说自洽场不收敛就行。6-31G*已经是正常情况能接受的最小基组了，不能再减了，再减几乎完全没有参考价值。金属的LanL2DZ也太小，建议换成SDD或者LanL2TZ（后者要自己拷定义）。B3LYP也不很适合这种体系，建议换成PBE0-D3(BJ)（写作PBE1PBE Em=GD3BJ）

作者
Author: 不明白 时间: 2021-4-23 17:36

hdhxx123 发表于 2021-4-23 15:41
不知道你比6-31g*的小基组是什么基组，建议你先看这篇博文
谈谈量子化学中基组的选择
http://sobereva.co ...

好的，谢谢

作者
Author: 不明白 时间: 2021-4-23 17:37

snljty 发表于 2021-4-23 15:42
你赝势都没读进去。先把gen换成genecp，不然计算没有意义。先别做优化，先把一个单点收敛了再说。不要说L50 ...

好，我记住了，我试一下

作者
Author: 不明白 时间: 2021-4-29 09:31
本帖最后由不明白于 2021-4-29 09:39 编辑

snljty 发表于 2021-4-23 15:42
你赝势都没读进去。先把gen换成genecp，不然计算没有意义。先别做优化，先把一个单点收敛了再说。不要说L50 ...

不收敛的解决方法我试了，3-21g可以收敛，6-31g*不能，这是哪里出现问题了，您知道怎么修改吗

作者
Author: sobereva 时间: 2021-5-3 02:53

不明白发表于 2021-4-29 09:31
不收敛的解决方法我试了，3-21g可以收敛，6-31g*不能，这是哪里出现问题了，您知道怎么修改吗

仔细看
解决SCF不收敛问题的方法
http://sobereva.com/61
没有可额外补充的

小基组绝大多数情况都比大基组更容易收敛
3-21G这么垃圾的基组即便收敛了也没什么实际用处

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