求助CCSD（T）模拟氟气键能

ZZU_SCU

要讲读书报告了，想给大家讲一些基础知识，想模拟氟气的性质（基态键长和解离能）和实验值对应起来。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

请大神不吝赐教

EmperorXiang

fffff 发表于 2021-5-30 23:08
老师，请问一下算出的解离能和实验差很多可以解释吗？是理论方法的缺陷还是说CCSD(T)方法没问题但算出来F ...

Full-CCSD(T)/cc-pVQZ//CCSD/def2-TZVP算出来的键焓是152kJ·mol-1 和实验值相差不多。此外实验本身也有误差。

SuXXXX

snljty 发表于 2021-5-31 09:28
你没做Freq原理上肯定不可能有热力学校正量。看卢老师博文《使用Shermo结合量子化学程序方便地计算分子的 ...

好的好的，谢谢老师！！我再去学习一下这些基础知识

snljty

喵星大佬 发表于 2021-5-31 10:12
对于F2，有没有DLPNO计算量应该是没啥区别

他说的苯的单点 不过我也脑抽，看着苯，数原子个数当成两个苯了...

喵星大佬

snljty 发表于 2021-5-31 09:28
你没做Freq原理上肯定不可能有热力学校正量。看卢老师博文《使用Shermo结合量子化学程序方便地计算分子的 ...

对于F2，有没有DLPNO计算量应该是没啥区别

fffff

wzkchem5 发表于 2021-5-31 00:19
CCSD(T)应该不至于差很多。问题在于F2有一定的多参考态性质，所以HF以及很多DFT方法算出来的键能都特别不 ...

好的，谢谢老师。我去看看单参考态和多参考态的知识

snljty

SuXXXX 发表于 2021-5-31 09:24
各位老师，我还想请教一下，上周六我问完我就用苯试着算了一下气相单点：# CCSD(T)/6-31G* geom=connectivi ...

你没做Freq原理上肯定不可能有热力学校正量。看卢老师博文《使用Shermo结合量子化学程序方便地计算分子的各种热力学数据》 http://bbs.keinsci.com/thread-17494-1-1.html 相关内容。另外用6-31G*搭配CCSD(T)是完全荒谬的，至少要用到cc-pVTZ。如果算不动，就考虑DLPNO-CCSD(T)之类的。看卢老师博文 《谈谈量子化学中基组的选择》 http://bbs.keinsci.com/thread-3545-1-1.html

wzkchem5

fffff 发表于 2021-5-30 16:08
老师，请问一下算出的解离能和实验差很多可以解释吗？是理论方法的缺陷还是说CCSD(T)方法没问题但算出来F ...

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

fffff

喵星大佬 发表于 2021-4-8 23:05
F2的解离能算得和实验值会差很多

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

SuXXXX

hdhxx123 发表于 2021-5-29 22:59
对了，记得这篇也一起看了。泛函：http://sobereva.com/272

好哒！！谢谢您！！

hdhxx123

SuXXXX 发表于 2021-5-29 18:28
谢谢您的分享！！我看看

对了，记得这篇也一起看了。泛函：http://sobereva.com/272

SuXXXX

hdhxx123 发表于 2021-5-29 17:47
cc-pVTZ起步。基组的选择：谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com ...

谢谢您的分享！！我看看

hdhxx123

SuXXXX 发表于 2021-5-29 17:42
老师，我还想问一下，如果选用这个计算方法的话，基组应该如何选择，可以和B3LYP这些方法一样用def2tzvp ...

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

SuXXXX

北大-陶豫 发表于 2021-5-29 17:31
就是一种计算方法的名字。初学的话，你就把它当做一个可以替换你输入文件里的B3LYP的东西就行。注意这东 ...

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

SuXXXX

北大-陶豫 发表于 2021-5-29 17:31
就是一种计算方法的名字。初学的话，你就把它当做一个可以替换你输入文件里的B3LYP的东西就行。注意这东 ...

好的好的，谢谢老师！！