计算化学公社

标题: ccsd做irc求助 [打印本页]

作者
Author: 小西瓜123 时间: 2022-1-22 19:01
标题: ccsd做irc求助
本帖最后由小西瓜123 于 2022-1-23 19:52 编辑

各位老师，我有一个很小的体系（O2H2）用ccsd/avtz找到了过渡态，用CCSD/avtz IRC(gradientonly,euler)做IRC的时候只算了5个点就Normal结束了，于是我用freq，先做frequency,存chk文件，然后再做irc，输入文件如下
%chk=r3-ts-irc-test
%nprocshared=8
%mem=16gb
# ccsd(maxcycle=500)/aug-cc-pvtz freq

freq

0 1
O                0.05120100 1.24334900 0.00000000
H                0.06339300 0.07509000 0.00000000
O                0.05120100 -1.08448300 0.00000000
H                -0.88260300 -1.34602500 0.00000000

--Link1--
%chk=r3-ts-irc-test
# CCSD/aug-cc-pVtZ Geom=Check Guess=Read IRC=(Maxpoints=40,readcartesianFC,lqa)

irc

0 1

可是这个时候就会报错Analytic second derivatives necessary but not available.需要二阶解析导数，我的二阶解析导数已经存到了chk里，为啥就读不了呢？请问老师这种情况应该怎么办呢？

作者
Author: LittlePupil 时间: 2022-1-22 23:50
"freq=numer"是做数值二阶梯度。

作者
Author: 小西瓜123 时间: 2022-1-23 08:13

LittlePupil 发表于 2022-1-22 23:50
"freq=numer"是做数值二阶梯度。

嗯嗯对，已修正，我直接用的freq做的，但还是报错

作者
Author: LittlePupil 时间: 2022-1-23 16:26

小西瓜123 发表于 2022-1-23 08:13
嗯嗯对，已修正，我直接用的freq做的，但还是报错

因为G16里没有CCSD的二阶解析导数。

Analytic frequencies are available for the AM1, PM3, PM3MM, PM6, PDDG, DFTB, DFTBA, HF, DFT, MP2, CIS, TD and CASSCF methods.

Numerical frequencies are available for MP3, MP4(SDQ), CID, CISD, CCD, CCSD, EOM-CCSD and QCISD.

作者
Author: LittlePupil 时间: 2022-1-23 19:28
另外，你的描述“O1H3”和你贴的输入文件里的不一样，而且我按你的输入文件内容做波函数稳定性测试，结果显示不稳定。

作者
Author: 小西瓜123 时间: 2022-1-23 19:53
本帖最后由小西瓜123 于 2022-1-23 19:58 编辑

LittlePupil 发表于 2022-1-23 16:26
因为G16里没有CCSD的二阶解析导数。

原来如此，我用g09试了也不行，那哪个版本有ccsd的二阶解析导数呢

作者
Author: 小西瓜123 时间: 2022-1-23 19:54
本帖最后由小西瓜123 于 2022-2-17 10:10 编辑

LittlePupil 发表于 2022-1-23 19:28
另外，你的描述“O1H3”和你贴的输入文件里的不一样，而且我按你的输入文件内容做波函数稳定性测试，结果显 ...

同时在试几个体系，已修正，那我还需要先stable=OPT，稳定波函数吗

作者
Author: sobereva 时间: 2022-1-23 21:54

小西瓜123 发表于 2022-1-23 19:53
原来如此，我用g09试了也不行，那哪个版本有ccsd的二阶解析导数呢

哪个版本都没有
CFOUR才支持CCSD的解析Hessian

作者
Author: LittlePupil 时间: 2022-1-23 23:52

小西瓜123 发表于 2022-1-23 19:54
同时在试几个体系，已修正，那我还需要先opt=stable，稳定波函数吗

是“stable=opt”.
不过我不清楚这样是否有意义，如果你确定体系是自旋单重态，或许应该用非限制方法结合“guess=(mix,always)”跑IRC.

作者
Author: 小西瓜123 时间: 2022-1-24 10:33

LittlePupil 发表于 2022-1-23 23:52
是“stable=opt”.
不过我不清楚这样是否有意义，如果你确定体系是自旋单重态，或许应该用非限制方法结 ...

好，谢谢老师

作者
Author: 小西瓜123 时间: 2022-1-24 10:33

sobereva 发表于 2022-1-23 21:54
哪个版本都没有
CFOUR才支持CCSD的解析Hessian

谢谢老师

作者
Author: Hamilton1758 时间: 2024-6-14 11:02
我目前也做了一些单重态计算可以交流吗

作者
Author: mizu-bai 时间: 2024-6-14 12:55
第一步 freq=num，第二步里 irc=gradientonly

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