计算化学公社

标题: Gaussian用eomccsd做t2态过渡态搜索求助 [打印本页]

作者
Author: 小西瓜123 时间: 2022-9-17 11:58
标题: Gaussian用eomccsd做t2态过渡态搜索求助
本帖最后由小西瓜123 于 2022-9-17 14:53 编辑

各位老师，我用eomccsd找一个小体系的t2态过渡态，将t1的结构扭转了一下二面角，但是还是出现报错，请问这种应该怎么处理呢？输入文件
%chk=r7-s1-03-ts%nprocshared=4
%mem=8gb
#p EOMCCSD(NStates=5,root=2,triplets,maxcycle=500)/6-31+g* opt=(ts,gediis,noeigentest,maxcycle=200) nosymm

s1-02-ts

0 1
O                1.98022700 -0.35878500 0.29802400
H                -0.00692600 0.10711500 0.18117700
O                -0.85064488 0.76361831 -0.34177421
O                -1.98477518 0.05407884 -0.18894093
H                2.92625400 -0.45382700 0.16535900

--link1--
%chk=r7-s1-03-ts
%nprocshared=4
%mem=8gb
#p EOMCCSD(NStates=5,root=2,triplets,maxcycle=500)/6-311++g** opt=(ts,gediis,noeigentest,maxcycle=200) scf=xqc freq guess=read geom=allcheck nosymm

s1-02-ts

0 1
输出文件
Left  Eigenvector
Alpha Singles Amplitudes
   I SymI A SymA Value
11    1 14    1 0.580394
Beta  Singles Amplitudes
   I SymI A SymA Value
11    1 14    1 -0.580394
Alpha-Beta  Doubles Amplitudes
   I SymI J SymJ A SymA B SymB Value
10    1 13    1 14    1 14    1 0.189920
11    1 13    1 14    1 14    1 0.263576
13    1 10    1 14    1 14    1 -0.189920
13    1 11    1 14    1 14    1 -0.263576
Ground to excited state transition densities written to RWF  633
SavETr:  write IOETrn= 770 NScale= 10 NData=  16 NLR=2 NState= 5 LETran=    180.

EOM-CCSD Z-AMPLITUDE ITERATIONS

  No EOM-CCSD gradients or exc-exc 1PDM for CS triplets
Error termination via Lnk1e in /software//g16/l913.exe at Sat Sep 17 03:27:42 2022.

作者
Author: sobereva 时间: 2022-9-18 04:16
程序功能限制
如果不肯用TDDFT，对这么点体系还不如考虑用CASSCF

作者
Author: 小西瓜123 时间: 2022-9-18 14:00

sobereva 发表于 2022-9-18 04:16
程序功能限制
如果不肯用TDDFT，对这么点体系还不如考虑用CASSCF

谢谢sob老师，用molpro做casscf总是优化不出来，给的活化空间是occ,16,closed，3，用的基组也很小6-31+g*,老师，Gaussian可以用casscf 做激发态优化吗？

作者
Author: sobereva 时间: 2022-9-19 00:28

小西瓜123 发表于 2022-9-18 14:00
谢谢sob老师，用molpro做casscf总是优化不出来，给的活化空间是occ,16,closed，3，用的基组也很小6-31+g* ...

可。没必要带弥散。普通中性体系用2-zeta带弥散还不如用3-zeta不带弥散的

作者
Author: 小西瓜123 时间: 2022-9-19 09:33

sobereva 发表于 2022-9-19 00:28
可。没必要带弥散。普通中性体系用2-zeta带弥散还不如用3-zeta不带弥散的

明白了，谢谢sob老师

作者
Author: zjxitcc 时间: 2022-9-19 10:50
本帖最后由 zjxitcc 于 2022-9-19 10:52 编辑

小西瓜123 发表于 2022-9-18 14:00
谢谢sob老师，用molpro做casscf总是优化不出来，给的活化空间是occ,16,closed，3，用的基组也很小6-31+g* ...

这个不仅跟复杂的计算流程、软件有关，还跟你要找的这个过渡态是否真的存在有关，我没看出来T2下保证能找到这个过渡态。。。顺带说一下这体系有明显至少双自由基特征，用闭壳层EOM-CCSD做T2计算应该是定性错误的，至于是否报错那是次要的

作者
Author: 小西瓜123 时间: 2022-9-19 11:27

zjxitcc 发表于 2022-9-19 10:50
这个不仅跟复杂的计算流程、软件有关，还跟你要找的这个过渡态是否真的存在有关，我没看出来T2下保证能找 ...

老师，怎么判断这个过渡态是否存在呢？我也试过ueomccsd guess=(mix,always)关键词，也算不出来

作者
Author: zjxitcc 时间: 2022-9-19 12:39

小西瓜123 发表于 2022-9-19 11:27
老师，怎么判断这个过渡态是否存在呢？我也试过ueomccsd guess=(mix,always)关键词，也算不出来

没法严格判断是否真的存在，如果你用CASSCF经过合理的计算步骤也找不到，那大概率就是没有你要的过渡态

作者
Author: 小西瓜123 时间: 2022-9-19 15:07

zjxitcc 发表于 2022-9-19 12:39
没法严格判断是否真的存在，如果你用CASSCF经过合理的计算步骤也找不到，那大概率就是没有你要的过渡态

老师，我用molpro找这个反应的t2态的输入文件如下：***,o2

gprint,orbitals=30
memory,6000,m
angstrom
gempty=xyz
geometry={
O                1.98022700 -0.35878500 0.29802400
H                -0.00692600 0.10711500 0.18117700
O                -0.85064488 0.76361831 -0.34177421
O                -1.98477518 0.05407884 -0.18894093
H                2.92625400 -0.45382700 0.16535900
}
basis=6-31+g*

hf

nosym
{multi;maxiter,100;
occ,16
closed,3
wf,26,1,2;state,2
CPMCSCF,GRAD,2.1,spin=1.0}

optg,root=2

{frequencies
print}
您看这对吗？需要把wf,26,1,2;改为wf,26,1,0吗？

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