应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » MP2下波函数的RHF to UHF instability，用UDFT计算却不能自旋破缺
111/2返回列表12下一页
查看: 1976  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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朵）

423444 lzhsh1983

» 猜你喜欢
1楼 2014-08-10 11:17:15
已阅   回复此楼   关注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 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 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
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭