计算化学公社

 找回密码 Forget password
 注册 Register
Views: 677|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] Gaussian过渡态优化报错,计算意外停止

[复制链接 Copy URL]

6

帖子

0

威望

43

eV
积分
49

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
在gaussian过渡态优化报错,计算意外停止,不知道是什么原因,还望各位老师帮忙解决,下面是输入文件的开头


#b3lyp/def2TZVP em=gd3bj opt=(ts,calcfc,noeigen,recalc=10) freq

Title Card Required

0 2





以下是log输出文件的最后一部分


Stoichiometry    C12H18O4PS2(2)
Framework group  C1[X(C12H18O4PS2)]
Deg. of freedom   105
Full point group                 C1      NOp   1
Largest Abelian subgroup         C1      NOp   1
Largest concise Abelian subgroup C1      NOp   1
                         Standard orientation:                        
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
      1          6           0       -2.879607    1.088881   -0.888650
      2          6           0       -1.795839   -0.282840    0.791303
      3          6           0       -2.475658    0.405482    1.776337
      4          6           0       -3.360153    1.435166    1.446946
      5          6           0       -3.556412    1.767905    0.104841
      6          6           0       -1.972551    0.040764   -0.579223
      7          6           0       -1.288995   -0.611267   -1.625813
      8          1           0       -3.028185    1.356398   -1.928026
      9          1           0       -1.106357   -1.065144    1.072247
     10          1           0       -2.316292    0.144713    2.815585
     11          1           0       -4.238162    2.567104   -0.159023
     12          1           0       -1.480486   -0.269930   -2.636515
     13          6           0       -0.294243   -1.680583   -1.471906
     14         16           0        1.472130   -1.127557   -1.838144
     15         15           0        2.010412   -0.019210   -0.169454
     16         16           0        3.869057    0.520542   -0.310183
     17          8           0        1.608660   -0.968271    1.037812
     18          8           0        0.974535    1.168537    0.042000
     19          6           0        1.615612   -0.497624    2.420535
     20          6           0        0.882268    2.309439   -0.849061
     21          6           0        0.115754    3.392842   -0.135316
     22          6           0        1.225156   -1.655807    3.299874
     23          1           0        2.617923   -0.135051    2.652405
     24          1           0        0.909185    0.327512    2.502484
     25          1           0        1.890904    2.634615   -1.107477
     26          1           0        0.371228    1.982982   -1.755721
     27          1           0        0.008119    4.254123   -0.798069
     28          1           0       -0.879700    3.046587    0.144594
     29          1           0        0.645217    3.713427    0.763705
     30          1           0        1.227178   -1.330100    4.342253
     31          1           0        0.224044   -2.013883    3.053796
     32          1           0        1.931505   -2.481060    3.195159
     33          1           0       -0.280137   -2.150510   -0.493690
     34          1           0       -0.421614   -2.457266   -2.227765
     35          1           0       -3.887521    1.971481    2.225322
     36          8           0       -2.799286   -1.972145   -1.099478
     37          8           0       -2.920663   -2.524340   -0.041012
---------------------------------------------------------------------
Rotational constants (GHZ):           0.3942314           0.2224682           0.2119664
Standard basis: def2TZVP (5D, 7F)
There are   810 symmetry adapted cartesian basis functions of A   symmetry.
There are   715 symmetry adapted basis functions of A   symmetry.
   715 basis functions,  1167 primitive gaussians,   810 cartesian basis functions
    85 alpha electrons       84 beta electrons
       nuclear repulsion energy      2056.9488520158 Hartrees.
NAtoms=   37 NActive=   37 NUniq=   37 SFac= 1.00D+00 NAtFMM=   60 NAOKFM=F Big=F
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
Nuclear repulsion after empirical dispersion term =     2056.8570815788 Hartrees.
One-electron integrals computed using PRISM.
   1 Symmetry operations used in ECPInt.
ECPInt:  NShTT=   41328 NPrTT=  130309 LenC2=   37459 LenP2D=   82813.
LDataN:  DoStor=T MaxTD1= 6 Len=  172
NBasis=   715 RedAO= T EigKep=  6.40D-06  NBF=   715
NBsUse=   715 1.00D-06 EigRej= -1.00D+00 NBFU=   715
ExpMin= 9.04D-02 ExpMax= 6.07D+04 ExpMxC= 2.07D+03 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
Harris functional with IExCor=  402 and IRadAn=       5 diagonalized for initial guess.
HarFok:  IExCor=  402 AccDes= 0.00D+00 IRadAn=         5 IDoV= 1 UseB2=F ITyADJ=14
ICtDFT=  3500011 ScaDFX=  1.000000  1.000000  1.000000  1.000000
FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=           0
         NFxFlg=           0 DoJE=T BraDBF=F KetDBF=T FulRan=T
         wScrn=  0.000000 ICntrl=       500 IOpCl=  0 I1Cent=   200000004 NGrid=           0
         NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
Petite list used in FoFCou.
Initial guess <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7500 S= 0.5000
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Integral accuracy reduced to 1.0D-05 until final iterations.
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Problem detected with inexpensive integrals.
Switching to full accuracy and repeating last cycle.
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -1907.12243816     A.U. after  129 cycles
            NFock=128  Conv=0.11D-03     -V/T= 2.0029
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.9810 S= 0.6095
<L.S>=  0.00000000000   
Annihilation of the first spin contaminant:
S**2 before annihilation     0.9810,   after     0.7557
Convergence failure -- run terminated.
Error termination via Lnk1e in /mnt/app/g16avx2/g16/l502.exe at Thu Apr 11 11:33:14 2024.
Job cpu time:       0 days  8 hours 24 minutes 49.0 seconds.
Elapsed time:       0 days  1 hours  4 minutes  1.4 seconds.
File lengths (MBytes):  RWF=    488 Int=      0 D2E=      0 Chk=     17 Scr=      1
Error: segmentation violation
   rax 0x0000000000000000, rbx 0x0000000041a45b58, rcx 0x0000148bc7328127
   rdx 0x000000000000000f, rsp 0x00007ffe55ad5ab8, rbp 0x000000000326cd30
   rsi 0x000000000000000b, rdi 0x0000000000004f0d, r8  0x00000000018a23d0
   r9  0x00007ffe55ad5294, r10 0x0000148bc72e4468, r11 0x0000000000000206
   r12 0x00007ffe55ad5b70, r13 0x00007ffe55ad5be8, r14 0x00007ffe55ad5be0
   r15 0x0000000041a45b78
  /lib64/libpthread.so.0(+0x16910) [0x148bc7638910]
  /lib64/libc.so.6(kill+0x9) [0x148bc7328127]


5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

6#
发表于 Post on 2024-4-12 05:40:49 | 只看该作者 Only view this author
-ZQ- 发表于 2024-4-11 17:32
请问#p是这样添加吗

#P b3lyp/def2TZVP em=gd3bj opt=(ts,calcfc,noeigen,recalc=10) freq

显然是,手册里也都明确说了
这种常识性知识都没有的话建议系统学一遍扫扫盲
北京科音初级量子化学培训班(http://www.keinsci.com/workshop/KEQC_content.html

另外,opt freq目的用def2TZVP严重浪费,6-311G*足矣,看
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387http://bbs.keinsci.com/thread-6600-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

6

帖子

0

威望

43

eV
积分
49

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2024-4-11 17:32:18 | 只看该作者 Only view this author
Kamistry 发表于 2024-4-11 17:19
文章为什么不看?

请问#p是这样添加吗

#P b3lyp/def2TZVP em=gd3bj opt=(ts,calcfc,noeigen,recalc=10) freq

35

帖子

3

威望

794

eV
积分
889

Level 4 (黑子)

4#
发表于 Post on 2024-4-11 17:19:15 | 只看该作者 Only view this author
Gaussian中对应SCF不收敛的错误提示是L502出错,并伴随Convergence criterion not met和Convergence failure -- run terminated.的具体提示(但从G16 B.01开始,仅提示Convergence criterion not met而不报错中断)。

文章为什么不看?

6

帖子

0

威望

43

eV
积分
49

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-4-11 17:11:19 | 只看该作者 Only view this author
Kamistry 发表于 2024-4-11 15:58
SCF不收敛,看http://sobereva.com/61。另外建议用#p输出更详细的信息。

请问会有其他的原因吗,好像不是SCF不收敛,只计算了一步就停住了

我把前面一部分输出补上,如下
# opt=(calcfc,ts,noeigen,recalc=10) freq b3lyp def2tzvp em=gd3bj
----------------------------------------------------------------
1/5=1,10=4,11=1,18=20,26=3,38=1,71=10/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=44,7=101,11=2,25=1,30=1,71=2,74=-5,124=41,140=1/1,2,3;
4//1;
5/5=2,38=5/2;
8/6=4,10=90,11=11/1;
11/6=1,8=1,9=11,15=111,16=1/1,2,10;
10/6=1,13=1/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7/10=1,25=1/1,2,3,16;
1/5=1,10=4,11=1,18=20,26=3,71=10/3(3);
2/9=110/2;
7/9=1,25=1,44=-1/16;
99//99;
2/9=110/2;
3/5=44,7=101,11=2,25=1,30=1,71=2,74=-5,124=41,140=1/1,2,3;
4/5=5,16=3,69=1/1;
5/5=2,38=5/2;
8/6=4,10=90,11=11/1;
11/6=1,8=1,9=11,15=111,16=1/1,2,10;
10/6=1,13=1/2;
7/10=1,25=1/1,2,3,16;
1/5=1,11=1,18=20,26=3,71=10/3(-8);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
7/9=1,25=1,44=-1/16;
99/9=1/99;
-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge =  0 Multiplicity = 2
C                     3.20566  -1.10962  -0.31364
C                     1.91582   0.91737  -0.63923
C                     2.79031   1.59379   0.18779
C                     3.87596   0.93542   0.77003
C                     4.07546  -0.42198   0.50903
C                     2.09271  -0.46364  -0.91454
C                     1.22002  -1.2109   -1.7317
H                     3.35919  -2.16475  -0.50617
H                     1.07564   1.44513  -1.06572
H                     2.62791   2.64575   0.38889
H                     4.91305  -0.94076   0.95895
H                     1.43422  -2.26587  -1.85767
C                     0.00885  -0.70426  -2.38976
S                    -1.59187  -1.33131  -1.61223
P                    -1.76321  -0.2378    0.14208
S                    -3.44481  -0.65566   1.01575
O                    -1.53637   1.25682  -0.34194
O                    -0.46492  -0.45057   1.0352
C                    -1.35901   2.35326   0.60667
C                    -0.15048  -1.7048    1.69263
C                     0.88793  -1.42589   2.7484
C                    -1.20998   3.62451  -0.18629
H                    -2.23278   2.38473   1.25888
H                    -0.47327   2.1449    1.20555
H                    -1.06486  -2.11194   2.12623
H                     0.22021  -2.39462   0.93335
H                     1.16266  -2.36189   3.23938
H                     1.78563  -0.99143   2.30731
H                     0.49758  -0.74168   3.50393
H                    -1.08008   4.46232   0.50207
H                    -0.33704   3.57816  -0.83962
H                    -2.09686   3.8119   -0.79397
H                    -0.0656    0.37763  -2.4335
H                    -0.09182  -1.09983  -3.40178
H                     4.55564   1.47239   1.41912
O                     2.38852  -0.13882  -3.10838
O                     2.43008   1.05453  -3.22748

Add virtual bond connecting atoms O36        and C6         Dist= 4.23D+00.
Add virtual bond connecting atoms O36        and C7         Dist= 3.97D+00.

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Initialization pass.
                           ----------------------------
                           !    Initial Parameters    !
                           ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,5)                  1.3806         calculate D2E/DX2 analytically  !
! R2    R(1,6)                  1.4202         calculate D2E/DX2 analytically  !
! R3    R(1,8)                  1.0835         calculate D2E/DX2 analytically  !
! R4    R(2,3)                  1.3807         calculate D2E/DX2 analytically  !
! R5    R(2,6)                  1.4193         calculate D2E/DX2 analytically  !
! R6    R(2,9)                  1.08           calculate D2E/DX2 analytically  !
! R7    R(3,4)                  1.3968         calculate D2E/DX2 analytically  !
! R8    R(3,10)                 1.0833         calculate D2E/DX2 analytically  !
! R9    R(4,5)                  1.3966         calculate D2E/DX2 analytically  !
! R10   R(4,35)                 1.0824         calculate D2E/DX2 analytically  !
! R11   R(5,11)                 1.0831         calculate D2E/DX2 analytically  !
! R12   R(6,7)                  1.4099         calculate D2E/DX2 analytically  !
! R13   R(6,36)                 2.2374         calculate D2E/DX2 analytically  !
! R14   R(7,12)                 1.0838         calculate D2E/DX2 analytically  !
! R15   R(7,13)                 1.4686         calculate D2E/DX2 analytically  !
! R16   R(7,36)                 2.1            calculate D2E/DX2 analytically  !
! R17   R(13,14)                1.8868         calculate D2E/DX2 analytically  !
! R18   R(13,33)                1.0853         calculate D2E/DX2 analytically  !
! R19   R(13,34)                1.0912         calculate D2E/DX2 analytically  !
! R20   R(14,15)                2.0743         calculate D2E/DX2 analytically  !
! R21   R(15,16)                1.9405         calculate D2E/DX2 analytically  !
! R22   R(15,17)                1.5873         calculate D2E/DX2 analytically  !
! R23   R(15,18)                1.5901         calculate D2E/DX2 analytically  !
! R24   R(17,19)                1.4606         calculate D2E/DX2 analytically  !
! R25   R(18,20)                1.4506         calculate D2E/DX2 analytically  !
! R26   R(19,22)                1.5057         calculate D2E/DX2 analytically  !
! R27   R(19,23)                1.0908         calculate D2E/DX2 analytically  !
! R28   R(19,24)                1.0893         calculate D2E/DX2 analytically  !
! R29   R(20,21)                1.5069         calculate D2E/DX2 analytically  !
! R30   R(20,25)                1.0908         calculate D2E/DX2 analytically  !
! R31   R(20,26)                1.0908         calculate D2E/DX2 analytically  !
! R32   R(21,27)                1.0921         calculate D2E/DX2 analytically  !
! R33   R(21,28)                1.0905         calculate D2E/DX2 analytically  !
! R34   R(21,29)                1.0915         calculate D2E/DX2 analytically  !
! R35   R(22,30)                1.0921         calculate D2E/DX2 analytically  !
! R36   R(22,31)                1.0913         calculate D2E/DX2 analytically  !
! R37   R(22,32)                1.0913         calculate D2E/DX2 analytically  !
! R38   R(36,37)                1.2            calculate D2E/DX2 analytically  !
! A1    A(5,1,6)              121.2829         calculate D2E/DX2 analytically  !
! A2    A(5,1,8)              120.1084         calculate D2E/DX2 analytically  !
! A3    A(6,1,8)              118.6067         calculate D2E/DX2 analytically  !
! A4    A(3,2,6)              120.9289         calculate D2E/DX2 analytically  !
! A5    A(3,2,9)              119.3336         calculate D2E/DX2 analytically  !
! A6    A(6,2,9)              119.727          calculate D2E/DX2 analytically  !
! A7    A(2,3,4)              120.7344         calculate D2E/DX2 analytically  !
! A8    A(2,3,10)             119.4731         calculate D2E/DX2 analytically  !
! A9    A(4,3,10)             119.792          calculate D2E/DX2 analytically  !
! A10   A(3,4,5)              119.4216         calculate D2E/DX2 analytically  !
! A11   A(3,4,35)             120.2776         calculate D2E/DX2 analytically  !
! A12   A(5,4,35)             120.3007         calculate D2E/DX2 analytically  !
! A13   A(1,5,4)              120.3611         calculate D2E/DX2 analytically  !
! A14   A(1,5,11)             119.7459         calculate D2E/DX2 analytically  !
! A15   A(4,5,11)             119.8918         calculate D2E/DX2 analytically  !
! A16   A(1,6,2)              117.2698         calculate D2E/DX2 analytically  !
! A17   A(1,6,7)              119.2893         calculate D2E/DX2 analytically  !
! A18   A(1,6,36)             112.166          calculate D2E/DX2 analytically  !
! A19   A(2,6,7)              123.4374         calculate D2E/DX2 analytically  !
! A20   A(2,6,36)              93.7508         calculate D2E/DX2 analytically  !
! A21   A(6,7,12)             117.4471         calculate D2E/DX2 analytically  !
! A22   A(6,7,13)             125.97           calculate D2E/DX2 analytically  !
! A23   A(12,7,13)            116.5337         calculate D2E/DX2 analytically  !
! A24   A(12,7,36)            108.1043         calculate D2E/DX2 analytically  !
! A25   A(13,7,36)             89.3711         calculate D2E/DX2 analytically  !
! A26   A(7,13,14)            113.6013         calculate D2E/DX2 analytically  !
! A27   A(7,13,33)            114.7421         calculate D2E/DX2 analytically  !
! A28   A(7,13,34)            111.506          calculate D2E/DX2 analytically  !
! A29   A(14,13,33)           106.8393         calculate D2E/DX2 analytically  !
! A30   A(14,13,34)           100.57           calculate D2E/DX2 analytically  !
! A31   A(33,13,34)           108.5209         calculate D2E/DX2 analytically  !
! A32   A(13,14,15)           104.0872         calculate D2E/DX2 analytically  !
! A33   A(14,15,16)           109.8055         calculate D2E/DX2 analytically  !
! A34   A(14,15,17)           103.1024         calculate D2E/DX2 analytically  !
! A35   A(14,15,18)           109.6967         calculate D2E/DX2 analytically  !
! A36   A(16,15,17)           117.6376         calculate D2E/DX2 analytically  !
! A37   A(16,15,18)           115.2035         calculate D2E/DX2 analytically  !
! A38   A(17,15,18)           100.4037         calculate D2E/DX2 analytically  !
! A39   A(15,17,19)           121.7445         calculate D2E/DX2 analytically  !
! A40   A(15,18,20)           123.178          calculate D2E/DX2 analytically  !
! A41   A(17,19,22)           107.6842         calculate D2E/DX2 analytically  !
! A42   A(17,19,23)           108.2223         calculate D2E/DX2 analytically  !
! A43   A(17,19,24)           108.1922         calculate D2E/DX2 analytically  !
! A44   A(22,19,23)           111.7045         calculate D2E/DX2 analytically  !
! A45   A(22,19,24)           111.7499         calculate D2E/DX2 analytically  !
! A46   A(23,19,24)           109.155          calculate D2E/DX2 analytically  !
! A47   A(18,20,21)           107.8711         calculate D2E/DX2 analytically  !
! A48   A(18,20,25)           108.734          calculate D2E/DX2 analytically  !
! A49   A(18,20,26)           107.7584         calculate D2E/DX2 analytically  !
! A50   A(21,20,25)           111.6067         calculate D2E/DX2 analytically  !
! A51   A(21,20,26)           111.7505         calculate D2E/DX2 analytically  !
! A52   A(25,20,26)           108.9974         calculate D2E/DX2 analytically  !
! A53   A(20,21,27)           109.251          calculate D2E/DX2 analytically  !
! A54   A(20,21,28)           110.9548         calculate D2E/DX2 analytically  !
! A55   A(20,21,29)           110.766          calculate D2E/DX2 analytically  !
! A56   A(27,21,28)           108.4346         calculate D2E/DX2 analytically  !
! A57   A(27,21,29)           108.423          calculate D2E/DX2 analytically  !
! A58   A(28,21,29)           108.9443         calculate D2E/DX2 analytically  !
! A59   A(19,22,30)           109.1196         calculate D2E/DX2 analytically  !
! A60   A(19,22,31)           111.0132         calculate D2E/DX2 analytically  !
! A61   A(19,22,32)           110.965          calculate D2E/DX2 analytically  !
! A62   A(30,22,31)           108.3474         calculate D2E/DX2 analytically  !
! A63   A(30,22,32)           108.4167         calculate D2E/DX2 analytically  !
! A64   A(31,22,32)           108.9045         calculate D2E/DX2 analytically  !
! A65   A(6,36,37)            104.2566         calculate D2E/DX2 analytically  !
! A66   A(7,36,37)            126.3005         calculate D2E/DX2 analytically  !
! D1    D(6,1,5,4)              0.0424         calculate D2E/DX2 analytically  !
! D2    D(6,1,5,11)           179.6449         calculate D2E/DX2 analytically  !
! D3    D(8,1,5,4)           -179.4457         calculate D2E/DX2 analytically  !
! D4    D(8,1,5,11)             0.1568         calculate D2E/DX2 analytically  !
! D5    D(5,1,6,2)              0.2755         calculate D2E/DX2 analytically  !
! D6    D(5,1,6,7)           -179.0694         calculate D2E/DX2 analytically  !
! D7    D(5,1,6,36)           107.102          calculate D2E/DX2 analytically  !
! D8    D(8,1,6,2)            179.7712         calculate D2E/DX2 analytically  !
! D9    D(8,1,6,7)              0.4262         calculate D2E/DX2 analytically  !
! D10   D(8,1,6,36)           -73.4024         calculate D2E/DX2 analytically  !
! D11   D(6,2,3,4)              0.2879         calculate D2E/DX2 analytically  !
! D12   D(6,2,3,10)          -179.4576         calculate D2E/DX2 analytically  !
! D13   D(9,2,3,4)            179.1006         calculate D2E/DX2 analytically  !
! D14   D(9,2,3,10)            -0.6449         calculate D2E/DX2 analytically  !
! D15   D(3,2,6,1)             -0.4393         calculate D2E/DX2 analytically  !
! D16   D(3,2,6,7)            178.876          calculate D2E/DX2 analytically  !
! D17   D(3,2,6,36)          -117.773          calculate D2E/DX2 analytically  !
! D18   D(9,2,6,1)           -179.2474         calculate D2E/DX2 analytically  !
! D19   D(9,2,6,7)              0.068          calculate D2E/DX2 analytically  !
! D20   D(9,2,6,36)            63.419          calculate D2E/DX2 analytically  !
! D21   D(2,3,4,5)              0.0433         calculate D2E/DX2 analytically  !
! D22   D(2,3,4,35)          -179.8568         calculate D2E/DX2 analytically  !
! D23   D(10,3,4,5)           179.7879         calculate D2E/DX2 analytically  !
! D24   D(10,3,4,35)           -0.1122         calculate D2E/DX2 analytically  !
! D25   D(3,4,5,1)             -0.2078         calculate D2E/DX2 analytically  !
! D26   D(3,4,5,11)          -179.8098         calculate D2E/DX2 analytically  !
! D27   D(35,4,5,1)           179.6922         calculate D2E/DX2 analytically  !
! D28   D(35,4,5,11)            0.0903         calculate D2E/DX2 analytically  !
! D29   D(1,6,7,12)             0.7227         calculate D2E/DX2 analytically  !
! D30   D(1,6,7,13)           178.0718         calculate D2E/DX2 analytically  !
! D31   D(2,6,7,12)          -178.5795         calculate D2E/DX2 analytically  !
! D32   D(2,6,7,13)            -1.2304         calculate D2E/DX2 analytically  !
! D33   D(1,6,36,37)         -114.6432         calculate D2E/DX2 analytically  !
! D34   D(2,6,36,37)            6.858          calculate D2E/DX2 analytically  !
! D35   D(6,7,13,14)         -107.8448         calculate D2E/DX2 analytically  !
! D36   D(6,7,13,33)           15.4833         calculate D2E/DX2 analytically  !
! D37   D(6,7,13,34)          139.3674         calculate D2E/DX2 analytically  !
! D38   D(12,7,13,14)          69.5258         calculate D2E/DX2 analytically  !
! D39   D(12,7,13,33)        -167.1461         calculate D2E/DX2 analytically  !
! D40   D(12,7,13,34)         -43.262          calculate D2E/DX2 analytically  !
! D41   D(36,7,13,14)         179.5167         calculate D2E/DX2 analytically  !
! D42   D(36,7,13,33)         -57.1552         calculate D2E/DX2 analytically  !
! D43   D(36,7,13,34)          66.7289         calculate D2E/DX2 analytically  !
! D44   D(12,7,36,37)        -178.1354         calculate D2E/DX2 analytically  !
! D45   D(13,7,36,37)          64.0617         calculate D2E/DX2 analytically  !
! D46   D(7,13,14,15)          75.2408         calculate D2E/DX2 analytically  !
! D47   D(33,13,14,15)        -52.3083         calculate D2E/DX2 analytically  !
! D48   D(34,13,14,15)       -165.5173         calculate D2E/DX2 analytically  !
! D49   D(13,14,15,16)        175.5421         calculate D2E/DX2 analytically  !
! D50   D(13,14,15,17)         49.3821         calculate D2E/DX2 analytically  !
! D51   D(13,14,15,18)        -56.8933         calculate D2E/DX2 analytically  !
! D52   D(14,15,17,19)       -169.9995         calculate D2E/DX2 analytically  !
! D53   D(16,15,17,19)         69.032          calculate D2E/DX2 analytically  !
! D54   D(18,15,17,19)        -56.7605         calculate D2E/DX2 analytically  !
! D55   D(14,15,18,20)        -69.5858         calculate D2E/DX2 analytically  !
! D56   D(16,15,18,20)         54.9022         calculate D2E/DX2 analytically  !
! D57   D(17,15,18,20)       -177.677          calculate D2E/DX2 analytically  !
! D58   D(15,17,19,22)       -178.3929         calculate D2E/DX2 analytically  !
! D59   D(15,17,19,23)        -57.4875         calculate D2E/DX2 analytically  !
! D60   D(15,17,19,24)         60.6646         calculate D2E/DX2 analytically  !
! D61   D(15,18,20,21)       -163.4154         calculate D2E/DX2 analytically  !
! D62   D(15,18,20,25)        -42.2241         calculate D2E/DX2 analytically  !
! D63   D(15,18,20,26)         75.781          calculate D2E/DX2 analytically  !
! D64   D(17,19,22,30)        179.7509         calculate D2E/DX2 analytically  !
! D65   D(17,19,22,31)        -60.905          calculate D2E/DX2 analytically  !
! D66   D(17,19,22,32)         60.3509         calculate D2E/DX2 analytically  !
! D67   D(23,19,22,30)         61.0546         calculate D2E/DX2 analytically  !
! D68   D(23,19,22,31)       -179.6014         calculate D2E/DX2 analytically  !
! D69   D(23,19,22,32)        -58.3455         calculate D2E/DX2 analytically  !
! D70   D(24,19,22,30)        -61.5631         calculate D2E/DX2 analytically  !
! D71   D(24,19,22,31)         57.781          calculate D2E/DX2 analytically  !
! D72   D(24,19,22,32)        179.0369         calculate D2E/DX2 analytically  !
! D73   D(18,20,21,27)       -178.3995         calculate D2E/DX2 analytically  !
! D74   D(18,20,21,28)        -58.8999         calculate D2E/DX2 analytically  !
! D75   D(18,20,21,29)         62.2317         calculate D2E/DX2 analytically  !
! D76   D(25,20,21,27)         62.2162         calculate D2E/DX2 analytically  !
! D77   D(25,20,21,28)       -178.2843         calculate D2E/DX2 analytically  !
! D78   D(25,20,21,29)        -57.1527         calculate D2E/DX2 analytically  !
! D79   D(26,20,21,27)        -60.1265         calculate D2E/DX2 analytically  !
! D80   D(26,20,21,28)         59.373          calculate D2E/DX2 analytically  !
! D81   D(26,20,21,29)       -179.4954         calculate D2E/DX2 analytically  !
--------------------------------------------------------------------------------
Trust Radius=3.00D-01 FncErr=1.00D-07 GrdErr=1.00D-06 EigMax=2.50D+02 EigMin=1.00D-04
Number of steps in this run=    195 maximum allowed number of steps=    222.
Search for a saddle point of order  1.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

                          Input orientation:                          
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
      1          6           0        3.205661   -1.109621   -0.313638
      2          6           0        1.915820    0.917375   -0.639233
      3          6           0        2.790308    1.593787    0.187791
      4          6           0        3.875955    0.935423    0.770034
      5          6           0        4.075455   -0.421984    0.509026
      6          6           0        2.092705   -0.463641   -0.914541
      7          6           0        1.220015   -1.210903   -1.731704
      8          1           0        3.359190   -2.164752   -0.506175
      9          1           0        1.075640    1.445128   -1.065716
     10          1           0        2.627914    2.645748    0.388890
     11          1           0        4.913046   -0.940762    0.958954
     12          1           0        1.434224   -2.265865   -1.857674
     13          6           0        0.008849   -0.704258   -2.389764
     14         16           0       -1.591867   -1.331308   -1.612229
     15         15           0       -1.763207   -0.237797    0.142076
     16         16           0       -3.444810   -0.655659    1.015747
     17          8           0       -1.536366    1.256822   -0.341939
     18          8           0       -0.464919   -0.450572    1.035201
     19          6           0       -1.359008    2.353259    0.606666
     20          6           0       -0.150482   -1.704796    1.692628
     21          6           0        0.887932   -1.425886    2.748395
     22          6           0       -1.209982    3.624512   -0.186295
     23          1           0       -2.232781    2.384726    1.258879
     24          1           0       -0.473273    2.144902    1.205553
     25          1           0       -1.064862   -2.111938    2.126229
     26          1           0        0.220206   -2.394622    0.933354
     27          1           0        1.162659   -2.361885    3.239376
     28          1           0        1.785633   -0.991431    2.307314
     29          1           0        0.497580   -0.741680    3.503930
     30          1           0       -1.080079    4.462323    0.502065
     31          1           0       -0.337043    3.578158   -0.839616
     32          1           0       -2.096856    3.811899   -0.793975
     33          1           0       -0.065604    0.377629   -2.433498
     34          1           0       -0.091825   -1.099827   -3.401778
     35          1           0        4.555640    1.472395    1.419118
     36          8           0        2.388517   -0.138825   -3.108385
     37          8           0        2.430081    1.054527   -3.227480
---------------------------------------------------------------------
                    Distance matrix (angstroms):
                    1          2          3          4          5
     1  C    0.000000
     2  C    2.424544   0.000000
     3  C    2.780713   1.380663   0.000000
     4  C    2.409532   2.414226   1.396810   0.000000
     5  C    1.380638   2.788621   2.412079   1.396596   0.000000
     6  C    1.420227   1.419256   2.436144   2.824031   2.441224
     7  C    2.442122   2.491426   3.743871   4.233134   3.714405
     8  H    1.083487   3.405955   3.864173   3.392173   2.140306
     9  H    3.410184   1.079960   2.129194   3.386966   3.868438
    10  H    3.863953   2.133398   1.083252   2.151299   3.394229
    11  H    2.136157   3.871723   3.394797   2.152050   1.083109
    12  H    2.618957   3.442318   4.573817   4.807821   3.997157
    13  C    3.833305   3.054602   4.434108   5.256176   5.001994
    14  S    4.975114   4.278687   5.567738   6.380464   6.119243
    15  P    5.065314   3.934477   4.908288   5.794042   5.853080
    16  S    6.797213   5.826641   6.679986   7.495700   7.540939
    17  O    5.299781   3.481551   4.371987   5.534709   5.919043
    18  O    3.965712   3.216042   3.936243   4.564481   4.570850
    19  C    5.802991   3.786625   4.238995   5.426029   6.102860
    20  C    3.955127   4.072218   4.668343   4.902462   4.572206
    21  C    3.853302   4.245403   4.392511   4.291617   4.022797
    22  C    6.475039   4.159854   4.501789   5.832025   6.692791
    23  H    6.652813   4.792373   5.196560   6.297309   6.945052
    24  H    5.141441   3.258500   3.462734   4.535229   5.269249
    25  H    5.019458   5.070403   5.687906   5.961321   5.647490
    26  H    3.481258   4.039488   4.802988   4.947765   4.351356
    27  H    4.285549   5.134628   5.254403   4.932734   4.438739
    28  H    2.983260   3.513208   3.490717   3.231941   2.966714
    29  H    4.694988   4.682912   4.659161   4.658356   4.676841
    30  H    7.076687   4.779605   4.827749   6.088769   7.101833
    31  H    5.899379   3.492179   3.843646   5.227272   6.106563
    32  H    7.250432   4.950130   5.456029   6.811366   7.597434
    33  H    4.172157   2.727041   4.062811   5.109759   5.142587
    34  H    4.517756   3.966283   5.333597   6.106516   5.755005
    35  H    3.389940   3.393156   2.155760   1.082414   2.155810
    36  O    3.069329   2.726850   3.745418   4.290526   4.001450
    37  O    3.711542   2.642404   3.476298   4.252630   4.341523
                    6          7          8          9         10
     6  C    0.000000
     7  C    1.409874   0.000000
     8  H    2.159751   2.643448   0.000000
     9  H    2.168104   2.742059   4.308007   0.000000
    10  H    3.413749   4.620915   4.947405   2.442728   0.000000
    11  H    3.419350   4.577237   2.461549   4.951509   4.290671
    12  H    2.138015   1.083835   2.354202   3.811463   5.531355
    13  C    2.564495   1.468555   4.111658   2.740622   5.079661
    14  S    3.849115   2.816993   5.141104   3.888813   6.134167
    15  P    4.004435   3.654805   5.511110   3.514264   5.259058
    16  S    5.867449   5.442185   7.133584   5.401876   6.940477
    17  O    4.056849   3.952056   5.974999   2.716963   4.450220
    18  O    3.216070   3.327592   4.465206   3.221933   4.423859
    19  C    4.707793   4.982227   6.626630   3.090158   4.003564
    20  C    3.656458   3.721323   4.167024   4.362781   5.324138
    21  C    3.974229   4.497531   4.152742   4.777593   5.017271
    22  C    5.305759   5.628003   7.382089   3.278287   4.002282
    23  H    5.616652   5.813255   7.421813   4.151173   4.944833
    24  H    4.228893   4.770339   6.015888   2.836808   3.245789
    25  H    4.683310   4.573420   5.148260   5.236696   6.268215
    26  H    3.263381   3.082749   3.460969   4.412678   5.612382
    27  H    4.660831   5.102910   4.346581   5.747585   5.945471
    28  H    3.279208   4.084331   3.430522   4.221166   4.197482
    29  H    4.705807   5.306029   5.127854   5.098818   5.071140
    30  H    6.028138   6.516567   8.040011   4.025982   4.130613
    31  H    4.716506   5.114231   6.837714   2.568387   3.342093
    32  H    5.987253   6.091757   8.097631   3.967389   5.008246
    33  H    2.770066   2.160733   4.680544   2.076733   4.512790
    34  H    3.370942   2.126596   4.629048   3.646499   5.982932
    35  H    3.906437   5.315464   4.285684   4.276160   2.480777
    36  O    2.237401   2.100000   3.437744   2.899146   4.476839
    37  O    2.787176   2.972168   4.316539   2.580757   3.955914
                   11         12         13         14         15
    11  H    0.000000
    12  H    4.668136   0.000000
    13  C    5.943147   2.180236   0.000000
    14  S    7.005526   3.176612   1.886808   0.000000
    15  P    6.762677   4.282012   3.125377   2.074297   0.000000
    16  S    8.362910   5.886789   4.850529   3.285751   1.940541
    17  O    6.936614   4.850895   3.229104   2.883597   1.587329
    18  O    5.400797   3.907782   3.466871   3.009085   1.590125
    19  C    7.093190   5.933921   4.494222   4.307407   2.663230
    20  C    5.173137   3.928198   4.206232   3.624798   2.675291
    21  C    4.431590   4.713796   5.262529   5.017309   3.902947
    22  C    7.723004   6.669473   5.007901   5.171003   3.915522
    23  H    7.887437   6.692371   5.280087   4.739508   2.888836
    24  H    6.212452   5.698839   4.612641   4.612503   2.910699
    25  H    6.202384   4.705381   4.850630   3.855280   2.817260
    26  H    4.912953   3.046351   3.734316   3.300642   3.035118
    27  H    4.613601   5.105182   5.980487   5.673409   4.760852
    28  H    3.406076   4.369759   5.030107   5.185156   4.224983
    29  H    5.100282   5.652189   5.914042   5.557745   4.082538
    30  H    8.082056   7.560331   6.020133   6.188564   4.763128
    31  H    7.156740   6.190835   4.567460   5.125853   4.190367
    32  H    8.648657   7.109093   5.232228   5.232316   4.169838
    33  H    6.167158   3.093399   1.085327   2.434015   3.145503
    34  H    6.640032   2.464290   1.091230   2.346528   4.011922
    35  H    2.482502   5.869854   6.318125   7.405515   6.669589
    36  O    4.853823   2.645610   2.549304   4.416328   5.273718
    37  O    5.260469   3.727346   3.107645   4.947455   5.532421
                   16         17         18         19         20
    16  S    0.000000
    17  O    3.023748   0.000000
    18  O    2.987003   2.441251   0.000000
    19  C    3.683952   1.460644   2.973971   0.000000
    20  C    3.522989   3.851145   1.450572   4.371231   0.000000
    21  C    4.729479   4.756502   2.390915   4.890570   1.506897
    22  C    4.975864   2.395142   4.318968   1.505683   5.749297
    23  H    3.282083   2.078408   3.348773   1.090803   4.609588
    24  H    4.087689   2.076917   2.601072   1.089314   3.893791
    25  H    3.003010   4.202703   2.076153   4.725841   1.090809
    26  H    4.057477   4.247936   2.063759   5.014281   1.090764
    27  H    5.393005   5.762429   3.340736   5.960080   2.132729
    28  H    5.398002   4.807168   2.641171   4.895704   2.152968
    29  H    4.662713   4.787657   2.665665   4.628143   2.151370
    30  H    5.661233   3.346009   4.979879   2.129999   6.349405
    31  H    5.570080   2.659823   4.445442   2.153258   5.861456
    32  H    5.005112   2.654601   4.917089   2.152637   6.356531
    33  H    4.938017   2.703838   3.588487   3.849494   4.622620
    34  H    5.563661   4.123482   4.499724   5.440329   5.130536
    35  H    8.288457   6.345103   5.389917   6.034821   5.684799
    36  O    7.162631   5.000576   5.040693   5.835754   5.652301
    37  O    7.446077   4.909173   5.368128   5.544779   6.203277
                   21         22         23         24         25
    21  C    0.000000
    22  C    6.206462   0.000000
    23  H    5.145704   2.161415   0.000000
    24  H    4.121137   2.160839   1.776578   0.000000
    25  H    2.161287   6.186735   4.726131   4.395260   0.000000
    26  H    2.163026   6.287211   5.381940   4.600248   1.776024
    27  H    1.092077   7.293942   6.162929   5.208055   2.502682
    28  H    1.090491   6.041417   5.352131   4.019094   3.068167
    29  H    1.091488   5.966336   4.719061   3.815422   2.493370
    30  H    6.602278   1.092081   2.493571   2.496708   6.771931
    31  H    6.278124   1.091329   3.069492   2.501099   6.457798
    32  H    6.992263   1.091299   2.503899   3.068063   6.684643
    33  H    5.569014   4.111176   4.728503   4.065973   5.290328
    34  H    6.236254   5.823142   6.200607   5.648121   5.703510
    35  H    4.859948   6.360138   6.851327   5.078173   6.703543
    36  O    6.181436   5.970804   6.840884   5.658207   6.574206
    37  O    6.651448   5.394784   6.605992   5.410191   7.134657
                   26         27         28         29         30
    26  H    0.000000
    27  H    2.491391   0.000000
    28  H    2.511428   1.770588   0.000000
    29  H    3.068715   1.771266   1.775766   0.000000
    30  H    6.992456   7.687169   6.419865   6.211429   0.000000
    31  H    6.255242   7.360122   5.940544   6.182555   1.770299
    32  H    6.846411   8.062752   6.911128   6.777770   1.771047
    33  H    4.370670   6.418338   5.270361   6.068202   5.131414
    34  H    4.535110   6.875423   6.010851   6.940063   6.866888
    35  H    5.829732   5.433909   3.812121   5.071136   6.445307
    36  O    5.111345   6.836578   5.515451   6.903752   6.799784
    37  O    5.838897   7.422834   5.935925   7.229994   6.151732
                   31         32         33         34         35
    31  H    0.000000
    32  H    1.775855   0.000000
    33  H    3.585739   4.313726   0.000000
    34  H    5.339321   5.911500   1.766672   0.000000
    35  H    5.785712   7.390992   6.115313   7.173304   0.000000
    36  O    5.137309   6.409624   2.597096   2.676136   5.271692
    37  O    4.441574   5.832513   2.705003   3.321390   5.126743
                   36         37
    36  O    0.000000
    37  O    1.200000   0.000000

35

帖子

3

威望

794

eV
积分
889

Level 4 (黑子)

2#
发表于 Post on 2024-4-11 15:58:35 | 只看该作者 Only view this author
SCF不收敛,看http://sobereva.com/61。另外建议用#p输出更详细的信息。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 04:01 , Processed in 0.175936 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list