计算化学公社

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

[量化理论] 为什么GAUSSIAN中的restricted和unrestricted结果一样,而其他人却不是?

[复制链接 Copy URL]

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

我算了一个小体系,算甲烷在一个C-H键变化过程中的能量变化(the potential energy curve was obtained by altering the length of a single C–H bond while keeping the remaining three C–H bonds at their equilibrium bond length ~1.086 Å and the HCH angles at the tetrahedral value,Full configuration interaction potential energy curves for breaking bonds to hydrogen:
An assessment of single-reference correlation methods,
Antara Dutta and C. David Sherrill
),但是发现了一个奇怪的现象。
如下图所示,无论是B3LYP/MP2/CCSD/CCSD(T),我用gaussian算的restricted和unrestricted结果都一模一样,并且和文献中的restricted是一样的,但是他的unrestricted却有更好的结果,不知道是gaussian自己的处理有问题,还是我哪里设置得不对?贴上了我的一个log文件供检查。

另外,好像QCISD(T)的计算结果比CCSD(T)还好?这下为难了。。。






CH3-H.png (143.38 KB, 下载次数 Times of downloads: 140)

CH3-H.png

CH4_6_UCCSDT.log

1.2 MB, 下载次数 Times of downloads: 8

6万

帖子

99

威望

5万

eV
积分
120180

管理员

公社社长

34#
发表于 Post on 2015-2-5 19:25:43 | 只看该作者 Only view this author
那你看看RHF/6-31G*和对称破缺的UHF/6-31G*谁低。理应肯定是后者低,这样的话就等于RCCSD(T)的相关能比UCCSD(T)相关能更大才导致了这种情况。
北京科音自然科学研究中心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

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

33#
 楼主 Author| 发表于 Post on 2015-2-5 18:39:03 | 只看该作者 Only view this author
sobereva 发表于 2015-2-5 16:53
我不清楚你说的是什么偏离更远。

CBS-QB3里CCSD(T)/6-31+G(d')不是直接作为能量的,而只是用来获得高 ...

这样啊,好的,谢谢sob!

PS: 我指的偏离更远是因为之前在CCSD(T)/6-31G(d)下RCCSD(T)比UCCSD(T)的能量偏低,而按照“一般能量更低的作为基态或者稳定态”的说法,事实上应该偏高才对,但如果使用stable=opt之后RCCSD(T)的结果更低,就会与偏高的事实相去更远。

6万

帖子

99

威望

5万

eV
积分
120180

管理员

公社社长

32#
发表于 Post on 2015-2-5 16:53:25 | 只看该作者 Only view this author
诸葛壹次心 发表于 2015-2-5 16:42
这。。。更低的话岂不是偏得更远了?看来ccsd(t)/6-31g(d)这样的搭配果然很有问题呀。但我想起来好像CBS- ...

我不清楚你说的是什么偏离更远。

CBS-QB3里CCSD(T)/6-31+G(d')不是直接作为能量的,而只是用来获得高阶的相关能的校正。这和经常用的一种近似有一定类似:
CCSD(T)/大基组 ≈ CCSD(T)/小基组 - MP2/小基组 + MP2/大基组
北京科音自然科学研究中心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

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

31#
 楼主 Author| 发表于 Post on 2015-2-5 16:42:10 | 只看该作者 Only view this author
sobereva 发表于 2015-2-5 13:46
不稳定的话,用stable=opt,通常会产生出真正基态波函数,给出的能量应当比先前的更低。

这。。。更低的话岂不是偏得更远了?看来ccsd(t)/6-31g(d)这样的搭配果然很有问题呀。但我想起来好像CBS-QB3的方法中有一步就是CCSD(T)/6-31+G(d'),貌似两个基组也差不了多少,是不是说明键长较长时使用CBS-QB3的方法也要比较小心啊?

6万

帖子

99

威望

5万

eV
积分
120180

管理员

公社社长

30#
发表于 Post on 2015-2-5 13:46:39 | 只看该作者 Only view this author
诸葛壹次心 发表于 2015-2-5 13:29
这样子啊,我之前算的那个RCCSD(T)/6-31g(d)的case检验了stability,发现是不稳定的,我想这时虽然能量更 ...

不稳定的话,用stable=opt,通常会产生出真正基态波函数,给出的能量应当比先前的更低。
北京科音自然科学研究中心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

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

29#
 楼主 Author| 发表于 Post on 2015-2-5 13:29:56 | 只看该作者 Only view this author
helpme 发表于 2015-2-5 09:30
能量是计算化学里面的“黄金规则”,绝大多数情况下,以能量低的为准,能量低的为基态。

这样子啊,我之前算的那个RCCSD(T)/6-31g(d)的case检验了stability,发现是不稳定的,我想这时虽然能量更低应该也不能用作基态是吗?

296

帖子

1

威望

2588

eV
积分
2904

科音成员

28#
发表于 Post on 2015-2-5 09:30:16 | 只看该作者 Only view this author
诸葛壹次心 发表于 2015-2-4 14:22
为什么想要这个呢?能详细说明一下原因么?谢谢~

能量是计算化学里面的“黄金规则”,绝大多数情况下,以能量低的为准,能量低的为基态。
华北电力大学数理学院,理论与计算化学,团簇、表面的结构与反应机理。(招第一性原理计算,量子化学计算方向的教师、硕士/博士研究生)

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

27#
 楼主 Author| 发表于 Post on 2015-2-5 01:49:54 | 只看该作者 Only view this author
sobereva 发表于 2015-2-4 19:13
之所以有著名的kasha rule,正因为激发态能量相近,交叉频繁,所以...太普遍了

PS:最近恰有一篇用 ...

好的,谢谢sob~

6万

帖子

99

威望

5万

eV
积分
120180

管理员

公社社长

26#
发表于 Post on 2015-2-4 19:13:56 | 只看该作者 Only view this author
诸葛壹次心 发表于 2015-2-4 18:43
好的,谢谢sob~我赶紧去看看相关的例子~

PS: 看到“很普遍”这几个字顿时有种欲哭无泪的感觉,这简直 ...


之所以有著名的kasha rule,正因为激发态能量相近,交叉频繁,所以...太普遍了

PS:最近恰有一篇用了Multiwfn的研究含铂化合物的单-三重态交叉的文章有兴趣可以看看:http://www.sciencedirect.com/science/article/pii/S1566119915000312
北京科音自然科学研究中心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

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

25#
 楼主 Author| 发表于 Post on 2015-2-4 18:43:16 | 只看该作者 Only view this author
sobereva 发表于 2015-2-4 17:47
如果是自旋多重度不同的态的交叉,DFT实际上就可以。
对于自旋多重度相同的态的交叉(一般是避免交叉 ...

好的,谢谢sob~我赶紧去看看相关的例子~

PS: 看到“很普遍”这几个字顿时有种欲哭无泪的感觉,这简直就是一个无底洞啊。。。我最开始只是想研究下barrierless反应的。。。

6万

帖子

99

威望

5万

eV
积分
120180

管理员

公社社长

24#
发表于 Post on 2015-2-4 17:47:34 | 只看该作者 Only view this author
诸葛壹次心 发表于 2015-2-4 17:13
经典啊!原来是这样!谢谢sob~
那么这篇文中提到的“势能面交叉”如果存在的话,是否就需要用态平均(st ...


如果是自旋多重度不同的态的交叉,DFT实际上就可以。
对于自旋多重度相同的态的交叉(一般是避免交叉),离交叉点比较远的地方,用TDDFT之类普通研究激发态的方法都可以研究,但离交叉点近的区域,态平均CASSCF是标准方法。如果在CASSCF基础上进一步用二阶微扰,如CASPT2,就应该用其多态版本,单态版本只适合研究某个态。

势能面交叉是很普遍的现象,特别是对有机分子来说。随便找找光化学的资料都会有很多例子。
北京科音自然科学研究中心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

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

23#
 楼主 Author| 发表于 Post on 2015-2-4 17:13:19 | 只看该作者 Only view this author
sobereva 发表于 2015-2-4 15:03
RHF不能正确表现解离曲线,拉远后能量会比两个片段之和要高,见此贴讨论
小谈几种量子化学方法对氢分子 ...

经典啊!原来是这样!谢谢sob~
那么这篇文中提到的“势能面交叉”如果存在的话,是否就需要用态平均(state average)和多态(multi state)的方法来解决呢?(我前几天刚看到这个,但是还没搞太明白)另外一般什么体系下可能会出现势能面交叉呢?

6万

帖子

99

威望

5万

eV
积分
120180

管理员

公社社长

22#
发表于 Post on 2015-2-4 15:03:27 | 只看该作者 Only view this author
诸葛壹次心 发表于 2015-2-4 14:22
为什么想要这个呢?能详细说明一下原因么?谢谢~

RHF不能正确表现解离曲线,拉远后能量会比两个片段之和要高,见此贴讨论
小谈几种量子化学方法对氢分子解离过程的描述
http://hi.baidu.com/sobereva/item/fa2772ad3197ac17a9cfb7d8
相应地,由于RHF波函数作为参考态不正确,RCCSD(T)也没法正确描述解离曲线。
北京科音自然科学研究中心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

55

帖子

1

威望

208

eV
积分
283

Level 3 能力者

21#
 楼主 Author| 发表于 Post on 2015-2-4 14:22:30 | 只看该作者 Only view this author
helpme 发表于 2015-2-3 19:34
“现在是偏高”-------我想要的就是这句话。

为什么想要这个呢?能详细说明一下原因么?谢谢~

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

GMT+8, 2025-8-17 08:07 , Processed in 0.197733 second(s), 25 queries , Gzip On.

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