计算化学公社

标题: 求助CCSD（T）模拟氟气键能 [打印本页]

作者
Author: ZZU_SCU 时间: 2021-4-8 22:31
标题: 求助CCSD（T）模拟氟气键能
要讲读书报告了，想给大家讲一些基础知识，想模拟氟气的性质（基态键长和解离能）和实验值对应起来。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

请大神不吝赐教

作者
Author: RAL 时间: 2021-4-8 22:35
都已经自由基了，就不要写r了，直接ccsd(t)就好

作者
Author: 喵星大佬 时间: 2021-4-8 23:05
F2的解离能算得和实验值会差很多

作者
Author: ZZU_SCU 时间: 2021-4-9 12:18

RAL 发表于 2021-4-8 22:35
都已经自由基了，就不要写r了，直接ccsd(t)就好

谢谢您！问题已经解决。

作者
Author: sobereva 时间: 2021-4-10 05:56
顺带一提，写nosymm在此处毫无意义

作者
Author: SuXXXX 时间: 2021-5-29 17:25
想请问一下楼主，CCSD(T)要怎么用,我有点懵，看了帖子，但是可能看不懂。。。。它和平常高斯计算结构优化和单点所用的方法是一个意思吗。。

作者
Author: 北大-陶豫 时间: 2021-5-29 17:31

SuXXXX 发表于 2021-5-29 17:25
想请问一下楼主，CCSD(T)要怎么用,我有点懵，看了帖子，但是可能看不懂。。。。它和平常高斯计算结构优化和 ...

就是一种计算方法的名字。初学的话，你就把它当做一个可以替换你输入文件里的B3LYP的东西就行。注意这东西计算量极大，别妄图用它做优化（除非体系真的极小）。算单点能也算不了太大体系，估计十几个原子基本顶天了吧。

作者
Author: SuXXXX 时间: 2021-5-29 17:40

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

好的好的，谢谢老师！！

作者
Author: SuXXXX 时间: 2021-5-29 17:42

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

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

作者
Author: hdhxx123 时间: 2021-5-29 17:47
本帖最后由 hdhxx123 于 2021-5-29 17:49 编辑

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）

作者
Author: SuXXXX 时间: 2021-5-29 18:28

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

谢谢您的分享！！我看看

作者
Author: hdhxx123 时间: 2021-5-29 22:59

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

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

作者
Author: SuXXXX 时间: 2021-5-30 19:01

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

好哒！！谢谢您！！

作者
Author: fffff 时间: 2021-5-30 23:08

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

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

作者
Author: wzkchem5 时间: 2021-5-31 00:19

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

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

作者
Author: snljty 时间: 2021-5-31 09:28

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

作者
Author: fffff 时间: 2021-5-31 10:02

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

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

作者
Author: 喵星大佬 时间: 2021-5-31 10:12

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

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

作者
Author: snljty 时间: 2021-5-31 10:22
本帖最后由 snljty 于 2021-5-31 10:25 编辑

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

他说的苯的单点

不过我也脑抽，看着苯，数原子个数当成两个苯了...

作者
Author: SuXXXX 时间: 2021-5-31 10:35

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

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

作者
Author: EmperorXiang 时间: 2025-3-9 17:45

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

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

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