计算化学公社

标题: 寻找过渡态时结构不收敛，观察总能量时优化步数到最高点就停止了。 [打印本页]

作者
Author: YXY262626 时间: 2022-9-10 16:04
标题: 寻找过渡态时结构不收敛，观察总能量时优化步数到最高点就停止了。
请问各位老师，我在做过渡态结构的时候，通常都会遇到不收敛（l502错误）的情况，我用Gaussview看输出文件的能量到最高点就结束了，请问这种情况应该怎么解决？
下面是我的输入文件和输出文件。
%nprocshared=28
%chk=1bipycys.chk
%mem=100GB
# opt=(calcfc,ts,noeigen) freq lsda/sdd geom=connectivity

optbipycys

0 1
C                1.76530400 1.51688500 0.80084800
C                0.24324300 0.94782300 2.52715000
C                0.62271500 2.06409400 3.26760000
C                1.61603300 2.91662200 2.75788700
C                2.19511800 2.63487600 1.52274900
C                2.33658700 1.06998600 -0.46381700
C                2.28266600 -0.65309600 -2.08953800
C                3.33833800 -0.06682600 -2.78317500
C                3.92388400 1.10436000 -2.27653900
C                3.42469500 1.66942000 -1.10553500
N                0.81213200 0.69557400 1.33063100
N                1.81528900 -0.09141300 -0.95714800
Pt                0.35339800 -0.88695500 0.17509800
H                2.95679600 3.29175000 1.09324900
H                1.93224300 3.79823300 3.32498200
H                3.84877000 2.59201100 -0.69920900
H                4.76250800 1.57486200 -2.80009300
H                3.69955800 -0.53817600 -3.70130500
H                0.15092000 2.24863900 4.23671500
H                -0.49406100 0.20649400 2.86662700
H                1.78938500 -1.58549800 -2.39567900
Cl             -1.27174800 -1.76732700 1.63878900
Cl                0.29129900 2.79815000 -1.30999500
Cl                2.57229458 -2.46544445 1.63387339
Cl             -0.09601400 -2.65746400 -1.30599500
N                -1.09464400 0.18538800 -0.84623300
H                -1.39620700 -0.46456800 -1.61663100
C                -2.30667200 0.59610500 -0.10216100
H                -2.07745200 0.66183600 0.97595400
C                -2.83964600 1.95766100 -0.60637200
H                -2.28107600 2.27720100 -1.50494200
H                -2.66451100 2.73627200 0.15091500
S                -4.62810100 1.99802100 -1.12053400
H                -5.15850700 1.72731200 0.13368600
O                -3.27322300 -1.32084600 -1.22594700
H                -0.58119400 1.07116200 -1.19985200
C                -3.30930900 -0.52044700 -0.28279100
O                -4.25239600 -0.52532300 0.69937400
H                -4.89352400 -1.27491400 0.56593500

1 5 1.5 6 1.0 11 1.5
2 3 1.5 11 1.5 20 1.0
3 4 1.5 19 1.0
4 5 1.5 15 1.0
5 14 1.0
6 10 1.5 12 1.5
7 8 1.5 12 1.5 21 1.0
8 9 1.5 18 1.0
9 10 1.5 17 1.0
10 16 1.0
11 13 1.0
12 13 1.0
13 22 1.0 24 1.0 25 1.0 26 1.0
14
15
16
17
18
19
20
21
22
23
24
25
26 27 1.0 28 1.0
27
28 29 1.0 30 1.0 37 1.0
29
30 31 1.0 32 1.0 33 1.0
31
32
33 34 1.0
34
35 37 2.0
36
37 38 1.0
38 39 1.0
39

输出文件最后几行和4个NO
SCF Done:  E(RS-VWN) =  -3166.86515351    A.U. after  129 cycles
         NFock=128  Conv=0.34D-03    -V/T= 2.0316
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/g16/l502.exe at Sat Sep 10 14:54:03 2022.
Job cpu time:    0 days  3 hours 20 minutes 14.5 seconds.
Elapsed time:    0 days  0 hours  7 minutes 14.2 seconds.
File lengths (MBytes):  RWF= 663 Int=    0 D2E=    0 Chk=    11 Scr=    1

      Item             Value    Threshold  Converged?
Maximum Force          0.014983    0.000450    NO
RMS    Force          0.001035    0.000300    NO
Maximum Displacement    0.341593    0.001800    NO
RMS    Displacement    0.058256    0.001200    NO
Predicted change in Energy= 2.217827D-03

作者
Author: wzkchem5 时间: 2022-9-10 16:29
在SCF Done后面几行显示Convergence failure -- run terminated，表示是SCF不收敛。
检查结构、电荷、自旋多重度是否合理，如果不合理，改正后重算；如果合理，按照http://sobereva.com/61的方法解决

作者
Author: YXY262626 时间: 2022-9-10 16:33

wzkchem5 发表于 2022-9-10 16:29
在SCF Done后面几行显示Convergence failure -- run terminated，表示是SCF不收敛。
检查结构、电荷、自旋 ...

收到，感谢老师，我这就去改正。

作者
Author: ionexchangeC 时间: 2022-9-10 19:07
不要用SVWN

作者
Author: sobereva 时间: 2022-9-11 00:02
你的体系里除了Pt以外都绝对不要用SDD

把这些都仔细看了再做计算，否则纯粹是白浪费计算资源和时间
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
谈谈量子化学中基组的选择
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）
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

作者
Author: YXY262626 时间: 2022-9-14 09:31

sobereva 发表于 2022-9-11 00:02
你的体系里除了Pt以外都绝对不要用SDD

把这些都仔细看了再做计算，否则纯粹是白浪费计算资源和时间

收到，谢谢社长！

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