| 查看: 2428 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
限制开壳层计算单电能出错
|
||
| 我用WB97XD方法计算H2和一个含Fe的金属有机多孔聚合物的吸附能。由于体系具有磁性，自旋多重度是3，在计算中发现在H2距离聚合物较近的时候，自旋污染比较严重（这个我也不知道为什么，但是从结果看，距离近的时候自旋污染很大，距离远到3A以外，自选污染一下子就减小到可以忽略的程度了），因此我只好换ROWB97XD算，计算都是正常结束的，但是发现算的吸附能明显不对，由近及远，吸附能时大时小，只有两个点看起来是合理的，剩下的六个点从趋势上看明显是错的，而这两个点的输入文件和那六个点的输入文件，除氢分子与聚合物之间的距离以外完全一样。请问版上的各位高人帮忙分析一下问题出在什么地方呢？用限制开壳计算需要注意哪些问题呢？谢谢！~ |
回复此楼» 猜你喜欢
北京师范大学胡德平课题组以及袁萌伟课题组考研招生
已经有0人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有142人回复
-
祈福论文小修
已经有21人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有26人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有40人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
cygwin系统之下安装abinite不知什么原因出错了？
已经有6人回复- vasp跨节点运行出错，mpiexec_node-1 (handle_stdin_input 1089)
已经有5人回复
- td做激发态优化出错
已经有19人回复
- MS discover 计算出错，请高手指点一下
已经有4人回复
- 精修出错，请各位大虫帮忙，谢谢！！！
已经有4人回复
- vasp计算单个氧原子和分子的能量
已经有20人回复
- 为什么算能带时老是出错！
已经有4人回复
- WORD出错了该怎么修复啊？重装也没有用的情况！
已经有4人回复
- 人体尿液能够产生成本低廉电能
已经有15人回复
- 高斯计算单点能出错
已经有7人回复
- MS 能量最小化出错，急求帮助
已经有6人回复
- 请各位大侠看下高斯出错
已经有3人回复
- 大家用castep计算单原子O或N有没有出错？
已经有11人回复
- 关于限制开壳层计算的适用范围
已经有8人回复
- 【求助】【紧急求助】关于开壳层的轨道互换
已经有9人回复
9楼2011-05-16 21:16:02
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