24小时热门版块排行榜    

查看: 2016  |  回复: 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,以及巨大的单三态能量差。利用电子相关性?

回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

aioroslove

金虫 (小有名气)


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

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的回帖

423444

铜虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
您好,我是做有机的,合了几个双自由基,能给解释一下Open-Shell Broken-Symmetry Singlet 是什么意思吗?具体用高斯怎么计算的吗?我直接用DFT 算出来的闭壳层和单线态开壳层能量结构一直是一样的。真的很着急,实验做完了不知道怎么分析,方便的话可以邮件或者电话联系吗?我的邮箱:dengyj13@lzu.edu.cn
努力就好
6楼2016-01-15 20:44:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chordsword

木虫 (小有名气)

引用回帖:
6楼: Originally posted by 423444 at 2016-01-15 20:44:12
您好,我是做有机的,合了几个双自由基,能给解释一下Open-Shell Broken-Symmetry Singlet 是什么意思吗?具体用高斯怎么计算的吗?我直接用DFT 算出来的闭壳层和单线态开壳层能量结构一直是一样的。真的很着急,实 ...

最近放寒假没上小木冲,你怎么计算的单态开壳层?
7楼2016-02-04 19:35:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chordsword

木虫 (小有名气)

引用回帖:
6楼: Originally posted by 423444 at 2016-01-15 20:44:12
您好,我是做有机的,合了几个双自由基,能给解释一下Open-Shell Broken-Symmetry Singlet 是什么意思吗?具体用高斯怎么计算的吗?我直接用DFT 算出来的闭壳层和单线态开壳层能量结构一直是一样的。真的很着急,实 ...

open-shell broken symmetry singlet 简单的说就是开壳层的单态,阿尔法和贝塔轨道的单电子轨道是不一致的,但是由于自选耦合导致了单态,实际上是两个或多个单电子。  我后来采用的方法是分片段计算,把两个电子设置成相反的自旋态,然后做了NBO和MO分析,还算了自旋密度分布以说明是单态双自由基。

» 本帖已获得的红花(最新10朵)

8楼2016-02-04 19:40:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一盒火柴

金虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by chordsword at 2016-02-04 19:35:42
最近放寒假没上小木冲,你怎么计算的单态开壳层?...

实在抱歉,没想到您给我回了,我是的有机小分子,有一半结构本身就是不对称的,我用opt freq ub3lyp/6-31g(d) nosymm guess=mix 直接计算的,用opt freq ub3lyp/6-31g(d)计算三线态的,用opt freq rb3lyp/6-31g(d)计算闭壳层,能量结果基本是合适的,但是这是我自己捣鼓出来的,不知道对不对。我把log文件发给您,您帮忙看一下可以吗?

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : wuqudaiben-bssos-31-d-1.log
  • 2016-03-15 22:43:59, 1.1 M
9楼2016-03-15 22:44:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一盒火柴

金虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
引用回帖:
8楼: Originally posted by chordsword at 2016-02-04 19:40:33
open-shell broken symmetry singlet 简单的说就是开壳层的单态,阿尔法和贝塔轨道的单电子轨道是不一致的,但是由于自选耦合导致了单态,实际上是两个或多个单电子。  我后来采用的方法是分片段计算,把两个电子设 ...

不好意思,我没看懂,我们实验室是做有机,只是做一点计算解释一下实验,不知道您方便不方便电话或者邮件联系一下?非常感谢!
10楼2016-03-15 22:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chordsword 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 318求调剂 +3 plum李子 2026-03-21 6/300 2026-03-22 11:35 by plum李子
[考研] 求调剂院校信息 +5 CX 330 2026-03-21 5/250 2026-03-22 11:32 by 杨杨杨紫
[考研] 285求调剂 +5 ytter 2026-03-22 5/250 2026-03-22 11:13 by 杨杨杨紫
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
信息提示
请填处理意见