计算化学公社

标题: 乙醇分子H迁移过渡态计算问题（CCSD/def2TZVP） [打印本页]

作者
Author: Machaophy 时间: 2020-1-24 08:59
标题: 乙醇分子H迁移过渡态计算问题（CCSD/def2TZVP）
本帖最后由 Machaophy 于 2020-1-24 09:04 编辑

新手计算乙醇分子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.

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

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

作者
Author: Machaophy 时间: 2020-1-25 00:06

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

而且你写modRedundant试图 ...

我知道错了...

谢谢Sob老师

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