| 查看: 2437 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
限制开壳层计算单电能出错
|
||
| 我用WB97XD方法计算H2和一个含Fe的金属有机多孔聚合物的吸附能。由于体系具有磁性，自旋多重度是3，在计算中发现在H2距离聚合物较近的时候，自旋污染比较严重（这个我也不知道为什么，但是从结果看，距离近的时候自旋污染很大，距离远到3A以外，自选污染一下子就减小到可以忽略的程度了），因此我只好换ROWB97XD算，计算都是正常结束的，但是发现算的吸附能明显不对，由近及远，吸附能时大时小，只有两个点看起来是合理的，剩下的六个点从趋势上看明显是错的，而这两个点的输入文件和那六个点的输入文件，除氢分子与聚合物之间的距离以外完全一样。请问版上的各位高人帮忙分析一下问题出在什么地方呢？用限制开壳计算需要注意哪些问题呢？谢谢！~ |
回复此楼» 猜你喜欢
招贤纳博（已结束）
已经有28人回复
-
表征测试机构后台私聊我
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有274人回复
-
Coordination Chemistry Reviews
已经有1人回复
-
请问现在还有电池材料方向的博导招人吗
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
11楼2011-05-16 21:29:41
mchen10
金虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 1508.2
- 散金: 41
- 红花: 6
- 帖子: 425
- 在线: 159.5小时
- 虫号: 695973
- 注册: 2009-02-04
- 性别: GG
- 专业: 理论和计算化学
2楼2011-05-14 23:12:15
3楼2011-05-15 00:43:19
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42652.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
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