|
|
请问各位老师,我在做过渡态结构的时候,通常都会遇到不收敛(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
|
|