乙醇分子H迁移过渡态计算问题（CCSD/def2TZVP）

Machaophy

新手计算乙醇分子H迁移过渡态问题，想利用CCSD/def2TZVP精确计算，但是没有正常结束，我看了输出文件好像也没有明显的报错。
末尾的这一句Form Post-SCF non-separable gradient here不太明白，
不知是因为 自洽场不收敛么？


还请大佬不吝赐教....
输入输出文件部分文本已贴出，输入输出文件也已上传。


输入文件：

%mem=1024MW
%nprocshared=16
#p opt=(ts,noeigen,modRedundant,maxcycle=300,maxstep=5) CCSD/def2TZVP freq pop=npa

1 2
C              
H                  1            B1
H                  1            B2    2            A1
H                  1            B3    2            A2    3            D1    0
C                  1            B4    2            A3    4            D2    0
H                  5            B5    1            A4    2            D3    0
H                  5            B6    1            A5    2            D4    0
O                  5            B7    1            A6    2            D5    0
H                  8            B8    5            A7    1            D6    0

   B1             1.09063538
   B2             1.09194610
   B3             1.09194752
   B4             1.49494918
   B5             1.12566918
   B6             1.14527502
   B7             1.37683408
   B8             0.98577332
   A1           108.86770223
   A2           108.86758789
   A3           108.84948443
   A4           115.92446219
   A5           134.29919408
   A6           113.65098054
   A7           112.36724449
   D1           118.72678175
   D2           120.63635751
   D3            56.53064154
   D4          -115.30664847
   D5           179.99825911
   D6           179.99785947



输出文件：

**********************
DD1Dir will call FoFMem   1 times, MxPair=       262
NAB=    90 NAA=    45 NBB=    36.
Norm of the A-vectors is  9.7405283D-07 conv= 1.00D-06.
RLE energy=       -0.5678213651
DE(Z)  -0.56782136D+00          E(Z)=    -0.15431151308D+03
NORM(A)=   0.10900458D+01
CI/CC converged in   20 iterations to DelEn=-2.80D-09 Conv= 1.00D-08 ErrA1= 9.74D-07 Conv= 1.00D-06

S**2, projected HF & approx projected MPn energies after annihilation of
unwanted spin states (see manual for definitions):

spins       (S**2,0) (S**2,1)    PUHF        PMP2        PMP3        PMP4
annihilated
s+1         0.74987  0.75006 -153.749970 -154.295878 -154.324361
s+1,s+2     0.75000  0.75000 -153.749938 -154.295850 -154.324337
s+1 to s+3  0.75000  0.75000 -153.749938 -154.295850 -154.324337
s+1 to s+4  0.75000  0.75000 -153.749938 -154.295850 -154.324337
s+1 to s+5  0.75000  0.75000 -153.749938
s+1 to s+6  0.75000  0.75000 -153.749938

Discarding MO integrals.
Leave Link  913 at Thu Jan 23 22:58:41 2020, MaxMem=  1073741824 cpu:      5169.5
(Enter /exports/share/G09E01/g09/l1111.exe)
Selected functions:
Compute CCSD 2PDM.
Form Post-SCF non-separable gradient here.
IO1PDM=    0 IOW=    0 IO2PDM=    0.
MDV=1073741824.
Frozen-core window:  NFC=   3 NFV=   0.

sobereva

当前写pop=npa根本毫无意义，你没写density，这样得到的只是HF波函数下的结果

而且你写modRedundant试图避免提供初始Hessian，然而此时收敛到过渡态的难度剧增。这么用CCSD属于瞎折腾。
这么简单的反应没有任何理由用CCSD找过渡态，完全是浪费时间，撑死了用双杂化也足够了。

Machaophy

sobereva 发表于 2020-1-24 17:39
当前写pop=npa根本毫无意义，你没写density，这样得到的只是HF波函数下的结果

而且你写modRedundant试图 ...

我知道错了...

谢谢Sob老师