计算化学公社

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

[Gaussian/gview] ccsd做irc求助

[复制链接 Copy URL]

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
#
本帖最后由 小西瓜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里,为啥就读不了呢?请问老师这种情况应该怎么办呢?

68

帖子

7

威望

1901

eV
积分
2109

Level 5 (御坂)

12#
发表于 Post on 2024-6-14 12:55:12 | 只看该作者 Only view this author
第一步 freq=num,第二步里 irc=gradientonly

23

帖子

0

威望

215

eV
积分
238

Level 3 能力者

11#
发表于 Post on 2024-6-14 11:02:40 | 只看该作者 Only view this author
我目前也做了一些单重态计算 可以交流吗

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2022-1-24 10:33:53 | 只看该作者 Only view this author
sobereva 发表于 2022-1-23 21:54
哪个版本都没有
CFOUR才支持CCSD的解析Hessian

谢谢老师

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2022-1-24 10:33:31 | 只看该作者 Only view this author
LittlePupil 发表于 2022-1-23 23:52
是“stable=opt”.
不过我不清楚这样是否有意义,如果你确定体系是自旋单重态,或许应该用非限制方法结 ...

好,谢谢老师

99

帖子

0

威望

3299

eV
积分
3398

Level 5 (御坂)

8#
发表于 Post on 2022-1-23 23:52:56 | 只看该作者 Only view this author
小西瓜123 发表于 2022-1-23 19:54
同时在试几个体系,已修正,那我还需要先opt=stable,稳定波函数吗

是“stable=opt”.
不过我不清楚这样是否有意义,如果你确定体系是自旋单重态,或许应该用非限制方法结合“guess=(mix,always)”跑IRC.
「狩りってのは先に焦せった方の負けだ、ハンターも獲物もな」

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

7#
发表于 Post on 2022-1-23 21:54:47 | 只看该作者 Only view this author
小西瓜123 发表于 2022-1-23 19:53
原来如此,我用g09试了也不行,那哪个版本有ccsd的二阶解析导数呢

哪个版本都没有
CFOUR才支持CCSD的解析Hessian
北京科音自然科学研究中心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!

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2022-1-23 19:54:41 | 只看该作者 Only view this author
本帖最后由 小西瓜123 于 2022-2-17 10:10 编辑
LittlePupil 发表于 2022-1-23 19:28
另外,你的描述“O1H3”和你贴的输入文件里的不一样,而且我按你的输入文件内容做波函数稳定性测试,结果显 ...

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

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2022-1-23 19:53:18 | 只看该作者 Only view this author
本帖最后由 小西瓜123 于 2022-1-23 19:58 编辑
LittlePupil 发表于 2022-1-23 16:26
因为G16里没有CCSD的二阶解析导数。

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

99

帖子

0

威望

3299

eV
积分
3398

Level 5 (御坂)

4#
发表于 Post on 2022-1-23 19:28:18 | 只看该作者 Only view this author
另外,你的描述“O1H3”和你贴的输入文件里的不一样,而且我按你的输入文件内容做波函数稳定性测试,结果显示不稳定。
「狩りってのは先に焦せった方の負けだ、ハンターも獲物もな」

99

帖子

0

威望

3299

eV
积分
3398

Level 5 (御坂)

3#
发表于 Post on 2022-1-23 16:26:43 | 只看该作者 Only view this author
小西瓜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.
「狩りってのは先に焦せった方の負けだ、ハンターも獲物もな」

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

2#
 楼主 Author| 发表于 Post on 2022-1-23 08:13:24 | 只看该作者 Only view this author
LittlePupil 发表于 2022-1-22 23:50
"freq=numer"是做数值二阶梯度。

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

99

帖子

0

威望

3299

eV
积分
3398

Level 5 (御坂)

楼主
发表于 Post on 2022-1-22 23:50:13 | 只看该作者 Only view this author
"freq=numer"是做数值二阶梯度。
「狩りってのは先に焦せった方の負けだ、ハンターも獲物もな」

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

GMT+8, 2024-11-24 14:25 , Processed in 0.176267 second(s), 22 queries , Gzip On.

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