计算化学公社

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

[Gaussian/gview] 求助CCSD(T)模拟氟气键能

[复制链接 Copy URL]

146

帖子

0

威望

3272

eV
积分
3418

Level 5 (御坂)

要讲读书报告了,想给大家讲一些基础知识,想模拟氟气的性质(基态键长和解离能)和实验值对应起来。Gaussian16的CCSD(T)优化了基态F2的构型,键长1.4129与文献报道实验值1.4131一致。通过键能=两个F自由基的能量-基态F2分子的能量  的方法计算键能时,计算单个氟原子总是报错。关键字# rccsd(t) cc-pvqz nosymm
Title Card Required

0 2
F                  0.00000000    0.00000000    0.20000000


报错信息
NBasis=    55 RedAO= T EigKep=  2.97D-02  NBF=    55
NBsUse=    55 1.00D-06 EigRej= -1.00D+00 NBFU=    55
RHF requested, but multiplicity is not 1.
Error termination via Lnk1e in /public3/home/sc50101/cloud/business/gaussian/gaussian-g16/g16/l401.exe at Thu Apr  8 18:42:32 2021.
Job cpu time:       0 days  0 hours  1 minutes 18.7 seconds.
Elapsed time:       0 days  0 hours  0 minutes  2.1 seconds.
File lengths (MBytes):  RWF=      6 Int=      0 D2E=      0 Chk=      1 Scr=      1

请大神不吝赐教

142

帖子

0

威望

3602

eV
积分
3744

Level 5 (御坂)

2#
发表于 Post on 2021-4-8 22:35:56 | 只看该作者 Only view this author
都已经自由基了,就不要写r了,直接ccsd(t)就好

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
ZZU_SCU + 4 谢谢

查看全部评分 View all ratings

1657

帖子

5

威望

4562

eV
积分
6319

Level 6 (一方通行)

喵星人

3#
发表于 Post on 2021-4-8 23:05:43 | 只看该作者 Only view this author
F2的解离能算得和实验值会差很多

146

帖子

0

威望

3272

eV
积分
3418

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2021-4-9 12:18:57 | 只看该作者 Only view this author
RAL 发表于 2021-4-8 22:35
都已经自由基了,就不要写r了,直接ccsd(t)就好

谢谢您!问题已经解决。

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

5#
发表于 Post on 2021-4-10 05:56:18 | 只看该作者 Only view this author
顺带一提,写nosymm在此处毫无意义
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

87

帖子

0

威望

724

eV
积分
811

Level 4 (黑子)

6#
发表于 Post on 2021-5-29 17:25:35 | 只看该作者 Only view this author
想请问一下楼主,CCSD(T)要怎么用,我有点懵,看了帖子,但是可能看不懂。。。。它和平常高斯计算结构优化和单点所用的方法是一个意思吗。。

602

帖子

1

威望

4788

eV
积分
5410

Level 6 (一方通行)

7#
发表于 Post on 2021-5-29 17:31:40 | 只看该作者 Only view this author
SuXXXX 发表于 2021-5-29 17:25
想请问一下楼主,CCSD(T)要怎么用,我有点懵,看了帖子,但是可能看不懂。。。。它和平常高斯计算结构优化和 ...

就是一种计算方法的名字。初学的话,你就把它当做一个可以替换你输入文件里的B3LYP的东西就行。注意这东西计算量极大,别妄图用它做优化(除非体系真的极小)。算单点能也算不了太大体系,估计十几个原子基本顶天了吧。
现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

87

帖子

0

威望

724

eV
积分
811

Level 4 (黑子)

8#
发表于 Post on 2021-5-29 17:40:22 | 只看该作者 Only view this author
北大-陶豫 发表于 2021-5-29 17:31
就是一种计算方法的名字。初学的话,你就把它当做一个可以替换你输入文件里的B3LYP的东西就行。注意这东 ...

好的好的,谢谢老师!!

87

帖子

0

威望

724

eV
积分
811

Level 4 (黑子)

9#
发表于 Post on 2021-5-29 17:42:54 | 只看该作者 Only view this author
北大-陶豫 发表于 2021-5-29 17:31
就是一种计算方法的名字。初学的话,你就把它当做一个可以替换你输入文件里的B3LYP的东西就行。注意这东 ...

老师,我还想问一下,如果选用这个计算方法的话,基组应该如何选择,可以和B3LYP这些方法一样用def2tzvp这样的基组吗

442

帖子

0

威望

2410

eV
积分
2852

Level 5 (御坂)

娃娃儿鱼

10#
发表于 Post on 2021-5-29 17:47:15 | 只看该作者 Only view this author
本帖最后由 hdhxx123 于 2021-5-29 17:49 编辑
SuXXXX 发表于 2021-5-29 17:42
老师,我还想问一下,如果选用这个计算方法的话,基组应该如何选择,可以和B3LYP这些方法一样用def2tzvp ...

cc-pVTZ起步。基组的选择:谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html)。
这里也有提到:详谈 Multiwfn 产生 ORCA 量子化学程序的输入文件的功能
http://sobereva.com/490http://bbs.keinsci.com/thread-13560-1-1.html

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
北大-陶豫 + 5 我很赞同

查看全部评分 View all ratings

真·探

87

帖子

0

威望

724

eV
积分
811

Level 4 (黑子)

11#
发表于 Post on 2021-5-29 18:28:44 | 只看该作者 Only view this author
hdhxx123 发表于 2021-5-29 17:47
cc-pVTZ起步。基组的选择:谈谈量子化学中基组的选择
http://sobereva.com/336(http://bbs.keinsci.com ...

谢谢您的分享!!我看看

442

帖子

0

威望

2410

eV
积分
2852

Level 5 (御坂)

娃娃儿鱼

12#
发表于 Post on 2021-5-29 22:59:14 | 只看该作者 Only view this author
SuXXXX 发表于 2021-5-29 18:28
谢谢您的分享!!我看看

对了,记得这篇也一起看了。泛函:http://sobereva.com/272
真·探

87

帖子

0

威望

724

eV
积分
811

Level 4 (黑子)

13#
发表于 Post on 2021-5-30 19:01:24 | 只看该作者 Only view this author
hdhxx123 发表于 2021-5-29 22:59
对了,记得这篇也一起看了。泛函:http://sobereva.com/272

好哒!!谢谢您!!

95

帖子

0

威望

2445

eV
积分
2540

Level 5 (御坂)

14#
发表于 Post on 2021-5-30 23:08:04 | 只看该作者 Only view this author
喵星大佬 发表于 2021-4-8 23:05
F2的解离能算得和实验值会差很多

老师,请问一下算出的解离能和实验差很多可以解释吗?是理论方法的缺陷还是说CCSD(T)方法没问题但算出来F2解离能就是和实验值不一样。

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

15#
发表于 Post on 2021-5-31 00:19:18 | 只看该作者 Only view this author
fffff 发表于 2021-5-30 16:08
老师,请问一下算出的解离能和实验差很多可以解释吗?是理论方法的缺陷还是说CCSD(T)方法没问题但算出来F ...

CCSD(T)应该不至于差很多。问题在于F2有一定的多参考态性质,所以HF以及很多DFT方法算出来的键能都特别不准。CCSD(T)原则上也是不适用于多参考态体系,但是当多参考态性质不大的时候比HF、DFT还是准很多的,具体到F2我觉得CCSD(T)还是能用的

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
fffff + 5 谢谢

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2025-8-17 01:43 , Processed in 0.167177 second(s), 21 queries , Gzip On.

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