计算化学公社

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

[量化理论] CIVECTOR里列出的权重比接近1的两个组态DFT能量和CCSD(T)能量差很大

[复制链接 Copy URL]

43

帖子

0

威望

212

eV
积分
255

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
算了一个态平均CASSCF作业,第二个根里的CIVECTOR里列出来的两个组态权重比接近1,可是这个组态在HF, B3LYP,CCSD(T)算出来能差都很大。是否说明这个态和其他态平均到一起有问题导致对第二个根的波函数描述不好。如果是,要想state-specific地算第二根只需要把其它态权重设成0就行了吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2018-6-15 07:41:46 | 只看该作者 Only view this author
某个根里如果俩组态函数贡献都特别大,说明这个根本身不可能很好通过单个组态来描述,此时正体现出CASSCF的价值
“HF, B3LYP,CCSD(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

43

帖子

0

威望

212

eV
积分
255

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2018-6-16 16:57:58 | 只看该作者 Only view this author
sobereva 发表于 2018-6-15 07:41
某个根里如果俩组态函数贡献都特别大,说明这个根本身不可能很好通过单个组态来描述,此时正体现出CASSCF的 ...

有种说法说如果某个组态能量很高,那么波函数里它的贡献会很小,比如H2的sigma*轨道都双占的组态能量就属于高能组态,H2波函数这个组态贡献也很小。是否这种说法就不合理?另外,我这个体系不存在alpha,beta单子分开那种对称破缺情况。从CIVECTOR列出来的组态来看,选取任何一个当参考态,CCSD(T)通过激发也能构造出其它的组态形式,其实也考虑了其它几种组态,为何通过以第二个根中某个组态的HF为基础做的CCSD(T)与基态上相同手续做出来的CCSD(T)能量差与实验值比较大。是否CC的级别不够高?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2018-6-17 13:08:43 | 只看该作者 Only view this author
对于算稳定体系的基态来说(大多是单参考特征),若某组态自身能量显著高于作为参考态的HF波函数,那么它掺入基态波函数的比率应该较小,否则肯定会把这个态的能量拉高。

第二个问题看不明白。
至于为什么第二个根里两个组态比率甚至达到1:1,看看组态的构成以及赝正则轨道,试图分析原因。
北京科音自然科学研究中心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

43

帖子

0

威望

212

eV
积分
255

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2018-6-20 23:24:48 | 只看该作者 Only view this author
sobereva 发表于 2018-6-17 13:08
对于算稳定体系的基态来说(大多是单参考特征),若某组态自身能量显著高于作为参考态的HF波函数,那么它掺 ...

我想了一下,可能正是因为第二根两个混入50%成分HF能量很高的组态才是激发态(第二个根), 第一个根(基态)那个HF能量高的组态组成成分极小以至于没有在CIVECTOR里列出. 我用能量很高的组态做参考态做CC试图得到激发态能量干脆就是错误的,得到是基态波函数很差的基态能量

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 00:31 , Processed in 0.151521 second(s), 20 queries , Gzip On.

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