乙烯双自由基计算问题

atom · 发表于 Post on 2017-4-27 23:16:55

各位好：最近几天在看自由基的计算问题，发现自己还是很疑惑。Q1：看了sob其中一篇博文（http://sobereva.com/353），我自己做了乙烯双自由基的扭转，最后得出的wfn文件我用高斯是可以加载自旋密度的，但是multiwfn缺看不到？文件见附件。

Q2：双自由基的计算中涉及到自旋是1和3两种，我做了乙烷C-C键拉伸，发现乙烷一直是单态是基态，两个C-C距离拉到10埃，发现3态和自旋极化单重态能量是一致的，这是发生了简并？对于任意的两个自由基，怎么判断单态和三态是基态？
Q3：一个技术问题，在计算双自由基的时候，计算单态是不是直接uhf（或者其他泛函）+guess mix，计算三态是不是直接hf（或者其他泛函）？
由于鄙人对自由基这里很晕，所述可能存在欠缺，希望有人指条明路，兄弟谢过了！

liyuanhe211 · 发表于 Post on 2017-4-28 01:23:23

无限远时两个自由基的状态互不影响，当然1和3就一样了。
没有特别一般的规则能简单的判断“任意自由基”如何，即使简单有机物的情况单重态稳定在稠环芳烃里也有意外，过渡金属就更复杂了。
A3: 后者是。前者有时可以，有时也要用其他初猜才行。

冰释之川 · 发表于 Post on 2017-4-28 08:56:33

本帖最后由冰释之川于 2017-4-28 08:59 编辑

把 isovalue 设小一点，比如0.003，你就能看到等值面了！（一般做电子密度差，包括自旋密度，isovalue基本都在0.001~0.005量级的）

atom · 发表于 Post on 2017-4-28 11:15:38

liyuanhe211 发表于 2017-4-28 01:23
无限远时两个自由基的状态互不影响，当然1和3就一样了。
没有特别一般的规则能简单的判断“任意自由基”如 ...

嗯嗯好的谢谢李老师了

atom · 发表于 Post on 2017-4-28 11:16:15

冰释之川发表于 2017-4-28 08:56
把 isovalue 设小一点，比如0.003，你就能看到等值面了！（一般做电子密度差，包括自旋密度，isovalue基本 ...

的确是isovalue值的问题谢谢你了呢

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 乙烯双自由基计算问题

评分 Rate

评分 Rate

浏览过的版块