计算化学公社

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

[新手求助] 高斯计算CCSD(T)计算单点能和频率的续算问题

[复制链接 Copy URL]

9

帖子

0

威望

110

eV
积分
119

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
用高斯的ccsd(t)做单点能和频率计算的时候,计算了很多天了,然后断电了,怎么进行续算

516

帖子

1

威望

5054

eV
积分
5590

Level 6 (一方通行)

2#
发表于 Post on 2023-9-21 10:15:55 | 只看该作者 Only view this author
多大的体系?单点建议转投ORCA,用DLPNO-CCSD(T), 频率不知道你做的目的是什么?

9

帖子

0

威望

110

eV
积分
119

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-9-21 13:09:50 | 只看该作者 Only view this author
niobium 发表于 2023-9-21 10:15
多大的体系?单点建议转投ORCA,用DLPNO-CCSD(T), 频率不知道你做的目的是什么?

体系不大,体系一共八九个原子,但是计算资源不是很多。也算了一周多了。频率计算是想要体系的自由能。就想知道怎么续算,学习一下,以后也好用。

3811

帖子

4

威望

8000

eV
积分
11891

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2023-9-21 13:35:08 | 只看该作者 Only view this author
danr 发表于 2023-9-21 13:09
体系不大,体系一共八九个原子,但是计算资源不是很多。也算了一周多了。频率计算是想要体系的自由能。就 ...

问一下,为啥没提到结构优化,你结构优化收敛了,已经进入了freq阶段?那具体是CCSD振幅迭代阶段,还是(T)能量计算阶段,抑或是计算数值导数中某一阶段?
自动做多参考态计算的程序MOKIT

5万

帖子

99

威望

5万

eV
积分
112436

管理员

公社社长

5#
发表于 Post on 2023-9-21 14:12:31 | 只看该作者 Only view this author
opt freq撑死了用CCSD也就够了,用解析梯度都没有的CCSD(T)做等于自虐,而且连红外强度都得不到,而且结果比CCSD也几乎不会有可察觉的进步。尤其是算自由能的目的,热力学校正量带来的误差相对于电子能量的误差来说只是个零头,用CCSD(T)完全是费力不讨好
用CCSD(T)做振动分析一定要真正懂得自己在做什么,不是真懂就绝对不要做
北京科音自然科学研究中心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!

516

帖子

1

威望

5054

eV
积分
5590

Level 6 (一方通行)

6#
发表于 Post on 2023-9-21 20:15:03 | 只看该作者 Only view this author
zjxitcc 发表于 2023-9-21 13:35
问一下,为啥没提到结构优化,你结构优化收敛了,已经进入了freq阶段?那具体是CCSD振幅迭代阶段,还是(T ...

目测他不知道opt和freq要在同一级别。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
zjxitcc + 1 那样的话一周机时就白费了。。

查看全部评分 View all ratings

9

帖子

0

威望

110

eV
积分
119

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2023-9-22 10:31:20 | 只看该作者 Only view this author
sobereva 发表于 2023-9-21 14:12
opt freq撑死了用CCSD也就够了,用解析梯度都没有的CCSD(T)做等于自虐,而且连红外强度都得不到,而且结果 ...

谢谢sob老师

9

帖子

0

威望

110

eV
积分
119

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2023-9-22 10:34:31 | 只看该作者 Only view this author
niobium 发表于 2023-9-21 20:15
目测他不知道opt和freq要在同一级别。

所以就是说freq的话,必须和opt的方法一致是吗。不能就是说我用DFT优化的结构,然后用CC来计算freq得到自由能,这样是不准确的。但是DFT的结构做cc的单点是可以的吧

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
niobium + 1 看来我目测测对了

查看全部评分 View all ratings

9

帖子

0

威望

110

eV
积分
119

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2023-9-22 10:35:35 | 只看该作者 Only view this author
zjxitcc 发表于 2023-9-21 13:35
问一下,为啥没提到结构优化,你结构优化收敛了,已经进入了freq阶段?那具体是CCSD振幅迭代阶段,还是(T ...

额,想着用DFT的优化的结构来着

3811

帖子

4

威望

8000

eV
积分
11891

Level 6 (一方通行)

MOKIT开发者

10#
发表于 Post on 2023-9-22 10:47:40 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2023-9-22 10:49 编辑
danr 发表于 2023-9-22 10:34
所以就是说freq的话,必须和opt的方法一致是吗。不能就是说我用DFT优化的结构,然后用CC来计算freq得到自 ...

单点是可以。直接freq不行,因为DFT优化的几何结构在CCSD(T)下力不为零,未经CCSD(T)的opt、直接算freq得到热力学量无意义。
自动做多参考态计算的程序MOKIT

9

帖子

0

威望

110

eV
积分
119

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2023-11-25 11:22:18 | 只看该作者 Only view this author
zjxitcc 发表于 2023-9-22 10:47
单点是可以。直接freq不行,因为DFT优化的几何结构在CCSD(T)下力不为零,未经CCSD(T)的opt、直接算freq得 ...

我想问一下,我用DFT优化,用CCSD(T)计算单点能,那么CCSD(T)的单点能的零点能修正是加上DFT的零点能进行修正吗

3811

帖子

4

威望

8000

eV
积分
11891

Level 6 (一方通行)

MOKIT开发者

12#
发表于 Post on 2023-11-25 13:24:08 | 只看该作者 Only view this author
danr 发表于 2023-11-25 11:22
我想问一下,我用DFT优化,用CCSD(T)计算单点能,那么CCSD(T)的单点能的零点能修正是加上DFT的零点能 ...

是的
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 23:01 , Processed in 0.856236 second(s), 21 queries , Gzip On.

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