24小时热门版块排行榜    

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

gtolv8688

铜虫 (小有名气)

[求助] 限制开壳层计算单电能出错

我用WB97XD方法计算H2和一个含Fe的金属有机多孔聚合物的吸附能。由于体系具有磁性,自旋多重度是3,在计算中发现在H2距离聚合物较近的时候,自旋污染比较严重(这个我也不知道为什么,但是从结果看,距离近的时候自旋污染很大,距离远到3A以外,自选污染一下子就减小到可以忽略的程度了),因此我只好换ROWB97XD算,计算都是正常结束的,但是发现算的吸附能明显不对,由近及远,吸附能时大时小,只有两个点看起来是合理的,剩下的六个点从趋势上看明显是错的,而这两个点的输入文件和那六个点的输入文件,除氢分子与聚合物之间的距离以外完全一样。请问版上的各位高人帮忙分析一下问题出在什么地方呢?用限制开壳计算需要注意哪些问题呢?谢谢!~
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

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

zhou2009

版主 (著名写手)

【答案】应助回帖

gtolv8688(金币+10): 利用浅显易懂的语言解释了复杂的问题,非常感谢! 2011-05-16 11:26:40
2楼专家说,轨道排布很可能有误,我也有同感。

因为使用限制开壳RO计算,是规定安排了所有双电子轨道之后,将单电子轨道放在最高层。这在含有孤对电子的体系时,并不总是对的。
在含有孤对电子的体系,成键的电子对(MO)自然会低于单电子轨道,但是孤对电子轨道(MO)却因为是非键电子,能级可能会高于单电子轨道。
故含有单电子、孤对电子的体系,不能简单地用限制开壳RO计算。

我在算HCl的电子密度差时,Cl由于有一个单电子,计算时可用ROHF或者HF方法。
用ROHF方法算,单电子能级(-0.19956)在两个孤对电子简并能级之上,比H的轨道能级(-0.49982)还高出很多。体现不了Cl的电负性比H大。由此得到的密度差也不合理。
用HF方法算,总能量比ROHF低,单电子能级(-0.57891)在两个孤对电子简并能级之下,也在H的轨道能级之下,这样才合理。

你算题“自旋污染比较严重”,不能走限制开壳RO的路,要另找方法。

[ Last edited by zhou2009 on 2011-5-16 at 09:32 ]
4楼2011-05-16 09:25:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gtolv8688 的主题更新
信息提示
请填处理意见