24小时热门版块排行榜    

查看: 1980  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chordsword

木虫 (小有名气)

[交流] MP2下波函数的RHF to UHF instability,用UDFT计算却不能自旋破缺 已有4人参与

RT,最近一个体系中DFT的单三态能量差是10 kcal/mol,MP2下单三太能量差为30 kcal/mol,单态低。 审稿人提出MP2下可能存在RHF instability,测试了一下,MP2下果真有RHF to UHF instability,现在正在等待UMP2+guess=mix的结果,估计要好几天。但是DFT下测试比起哦啊曾是stable的,同时用gaussian09的分片段模块来打破对称性得到的依然是S2=0的闭壳层单态。根据symmetry broken理论,混合singlet和triplet,所以S2=0为闭壳层单态,S2=2是三态,数值越接近1,则单态双自由基的特性就越明显。而S2=0也就是说DFT方法无法打破对称性,这个单态应该是闭壳层无疑。 现在的问题是如何去解释在MP2下产生的RHF to UHF的instability,以及巨大的单三态能量差。利用电子相关性?

回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chordsword

木虫 (小有名气)

引用回帖:
2楼: Originally posted by aioroslove at 2014-08-11 21:33:56
如果用 stable=opt 找到了比基态更低的三重态能级,然后应该怎么办呢?
如果我只看单重态的跃迁,基态应该不能以这个三重态为参考态吧??

结果证实stable=opt得到的结构是不好的,优化之后又回到原结构。用gauss=mix和UMP2可以得到稳定的波函数,但是s2=0,并且能量和结构依然与闭壳层一致,说明这个结构应当的闭壳层。
3楼2014-08-11 23:50:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

aioroslove

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
如果用 stable=opt 找到了比基态更低的三重态能级,然后应该怎么办呢?
如果我只看单重态的跃迁,基态应该不能以这个三重态为参考态吧??
2楼2014-08-11 21:33:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aioroslove

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
似乎用stable=opt 优化后,得到的是一个单重态的开壳层,这样s2也是等于0.但是这样不符合电子排布的原理
4楼2014-08-12 10:47:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chordsword

木虫 (小有名气)

引用回帖:
4楼: Originally posted by aioroslove at 2014-08-12 10:47:47
似乎用stable=opt 优化后,得到的是一个单重态的开壳层,这样s2也是等于0.但是这样不符合电子排布的原理

我这个例子用guess=mix,scf=nosymm和UMP2优化opt后得到了和闭壳层一样几何构型和能量的单态,但wavefunction是稳定的。说明检测波函数的时候HF倾向于打破对称性,而MP2优化后表明这依然是闭壳层,S2=0.又采用了gaussian09的分片段,同样也无法打破对称性,可见是闭壳层无疑,只是在MP2下RHF是不稳定的。
5楼2014-08-12 12:54:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] [复试调剂]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见