计算化学公社

标题: 随便选的一批分子的T1诊断值 [打印本页]

作者
Author:
sobereva    时间: 2017-3-31 08:48
标题: 随便选的一批分子的T1诊断值
T1诊断对于判断静态相关强弱很有用。随便选了一批五花八门的分子,用g09做了下T1诊断。绝大部分都是B3LYP/def2-TZVP优化的结构,然后在def2-TZVP基组下做的T1诊断。对于T1诊断值大的体系很多还考虑了冻核和不冻核的影响。表格里面带U的是对称破缺计算的结果,R则是限制性闭壳层的结果。若未注明,算的都是体系的基态。对很多体系也测试了在def2-SVP下做诊断,数值和def2-TZVP差别很小,仅对含过渡金属差别略大但也定性一致。

(, 下载次数 Times of downloads: 550)

可见对于含过渡金属的体系冻不冻核往往影响挺大。如预期的,大部分过渡金属确实体现出了比普通有机体系更大的诊断值,但也并非配合物的诊断值一定大,如[Re2Cl8]2-、顺铂。也对比了个别有机反应的反应物、产物、TS的诊断值(HCN异构化、丁二烯闭环),有点TS下的数值明显更大,如前者,但并非都是如此,如后者。对于双自由基体系,如C4H8,可见UHF参考态明显比RHF更好,诊断值小得多。总之,结果基本都和一般理论直觉一致,不过CN自由基数值那么大有点意外。

作者
Author:
我本是个娃娃    时间: 2017-3-31 09:29
恭喜大博士加入测试大军,话说T1诊断是用来定量计算静态相关强弱的吧?如果用Grimme的NFOD值该怎么操作呢?
作者
Author:
sobereva    时间: 2017-3-31 09:54
我本是个娃娃 发表于 2017-3-31 09:29
恭喜大博士加入测试大军,话说T1诊断是用来定量计算静态相关强弱的吧?如果用Grimme的NFOD值该怎么操作呢?


RTFM
作者
Author:
我本是个娃娃    时间: 2017-3-31 09:58
sobereva 发表于 2017-3-31 09:54

RTFM


作者
Author:
moulihui    时间: 2019-8-8 10:24
请问不冻核的关键词是什么?
作者
Author:
sobereva    时间: 2019-8-9 03:59
moulihui 发表于 2019-8-8 10:24
请问不冻核的关键词是什么?

Gaussian里是Full,比如CCSD(T,full),代表做不冻核的CCSD(T)
作者
Author:
beefly    时间: 2019-8-9 05:59
sobereva 发表于 2019-8-9 03:59
Gaussian里是Full,比如CCSD(T,full),代表做不冻核的CCSD(T)

full一定要写在括号里边,有学生曾经用mp2算一个不大不小的体系,误把mp2(full)写成mp2 full,结果撞到一个gaussian的彩蛋,程序把bd(t)及其各种向下兼容的post-HF都算了一遍。原指望两三天出结果,结果半个月都没算完。
作者
Author:
moulihui    时间: 2019-8-9 11:01
sobereva 发表于 2019-8-9 03:59
Gaussian里是Full,比如CCSD(T,full),代表做不冻核的CCSD(T)

谢谢sob老师!
作者
Author:
moulihui    时间: 2019-8-9 11:04
beefly 发表于 2019-8-9 05:59
full一定要写在括号里边,有学生曾经用mp2算一个不大不小的体系,误把mp2(full)写成mp2 full,结果撞到一 ...

多谢提醒!
作者
Author:
Somnus荣小荣    时间: 2020-4-10 10:53
老师,请问一下计算T1诊断,最标准的计算是在什么水平下进行的呢,若使用不同的方法一般的误差能有多大(体系仅包含8个原子,碳氢氧硫)
作者
Author:
sobereva    时间: 2020-4-11 23:14
Somnus荣小荣 发表于 2020-4-10 10:53
老师,请问一下计算T1诊断,最标准的计算是在什么水平下进行的呢,若使用不同的方法一般的误差能有多大(体 ...

不需要大基组,CCSD/cc-pVDZ都足以得到定性正确的结果
作者
Author:
嘤嘤嘤    时间: 2021-7-14 20:48
请问T1诊断值是不是需要在原来的方法基组上进行计算?计算小白表示T1诊断值是什么物理意义都不知道。看了Sob老师说的还是完全不懂,感觉比结构化学还难懂QAQ
作者
Author:
sobereva    时间: 2021-7-14 22:49
嘤嘤嘤 发表于 2021-7-14 20:48
请问T1诊断值是不是需要在原来的方法基组上进行计算?计算小白表示T1诊断值是什么物理意义都不知道。看了So ...

结构化学跟计算量子化学根本没有可比性

我不知道你说的“原来”指的是什么。反正对基组没有特定的要求

参考北京科音中级量子化学培训班(http://www.keinsci.com/KBQC)的幻灯片:

(, 下载次数 Times of downloads: 48)

(, 下载次数 Times of downloads: 56)

没有点基本的量子化学、电子结构方面的知识不可能充分理解。


作者
Author:
嘤嘤嘤    时间: 2021-7-15 10:36
sobereva 发表于 2021-7-14 22:49
结构化学跟计算量子化学根本没有可比性

我不知道你说的“原来”指的是什么。反正对基组没有特定的要求 ...

原来指的是我用DFT方法的结构优化和CCSDT的单电能计算。我这两个都需要T1诊断值吗?还是说只有CCSDT需要?在T1计算的时候,方法是CCSD加关键词T1diage是吗?基组随便?
作者
Author:
biogon    时间: 2021-7-15 10:44
嘤嘤嘤 发表于 2021-7-15 10:36
原来指的是我用DFT方法的结构优化和CCSDT的单电能计算。我这两个都需要T1诊断值吗?还是说只有CCSDT需要 ...

你先搞清楚到底是CCSDT还是CCSD(T),怎么做T1诊断自个看用的程序的手册
作者
Author:
嘤嘤嘤    时间: 2021-7-15 16:01
biogon 发表于 2021-7-15 10:44
你先搞清楚到底是CCSDT还是CCSD(T),怎么做T1诊断自个看用的程序的手册

CCSD(T),好的吧我看了Gaussian的手册搜索T1出现了两个1 Computes the Q1 diagnostic of T. J. Lee and coworkers [195, 654]. Note that Q1 is analogous to the T1
diagnostic for CCSD when it is computed using QCISD instead of the Coupled Cluster method.
2 Computes the T1 diagnostic of T. J. Lee and coworkers [195] (CCSD only).
所以我看一下标注的参考文献就能懂了嘛?
作者
Author:
wzkchem5    时间: 2021-7-15 17:15
嘤嘤嘤 发表于 2021-7-15 09:01
CCSD(T),好的吧我看了Gaussian的手册搜索T1出现了两个1 Computes the Q1 diagnostic of T. J. Lee and ...

对,看参考文献[195],如果有不懂的,找一些量化教材,对着教材的coupled cluster的章节看
作者
Author:
嘤嘤嘤    时间: 2021-7-15 20:39
wzkchem5 发表于 2021-7-15 17:15
对,看参考文献[195],如果有不懂的,找一些量化教材,对着教材的coupled cluster的章节看

好的谢谢老师
作者
Author:
SS123    时间: 2022-1-10 17:52
sobereva 发表于 2019-8-9 03:59
Gaussian里是Full,比如CCSD(T,full),代表做不冻核的CCSD(T)

请问sob老师,对过渡金属团簇而言,在高斯里做冻核和不冻核的CCSD(T),两者算出的能量是不是也差别很大?
作者
Author:
sobereva    时间: 2022-1-11 03:44
SS123 发表于 2022-1-10 17:52
请问sob老师,对过渡金属团簇而言,在高斯里做冻核和不冻核的CCSD(T),两者算出的能量是不是也差别很大?


作者
Author:
糖糖糖豆9988    时间: 2022-11-7 09:14
本帖最后由 糖糖糖豆9988 于 2022-11-7 09:37 编辑

Sob老师您好,请问对于双自由基计算T1诊断值的时候,直接写上关键词,再输入 “读取人为构造对称破缺初猜波函数前提下,优化得到的分子坐标” 就可以吗?计算T1诊断值不需要在重新构造对称破缺态了吧?
作者
Author:
sobereva    时间: 2022-11-7 10:14
糖糖糖豆9988 发表于 2022-11-7 09:14
Sob老师您好,请问对于双自由基计算T1诊断值的时候,直接写上关键词,再输入 “读取人为构造对称破缺初猜波 ...






欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3