| 查看: 1869 | 回复: 10 | ||
[求助]
限制开壳层计算单电能出错
|
| 我用WB97XD方法计算H2和一个含Fe的金属有机多孔聚合物的吸附能。由于体系具有磁性，自旋多重度是3，在计算中发现在H2距离聚合物较近的时候，自旋污染比较严重（这个我也不知道为什么，但是从结果看，距离近的时候自旋污染很大，距离远到3A以外，自选污染一下子就减小到可以忽略的程度了），因此我只好换ROWB97XD算，计算都是正常结束的，但是发现算的吸附能明显不对，由近及远，吸附能时大时小，只有两个点看起来是合理的，剩下的六个点从趋势上看明显是错的，而这两个点的输入文件和那六个点的输入文件，除氢分子与聚合物之间的距离以外完全一样。请问版上的各位高人帮忙分析一下问题出在什么地方呢？用限制开壳计算需要注意哪些问题呢？谢谢！~ |
回复此楼» 猜你喜欢
清华大学化工系张吉松老师和杭州高等研究陈兴坤老师招聘联合培养博士生
已经有6人回复
-
法国 波尔多大学 CSC
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有179人回复
-
法国 波尔多大学 （化学） CSC
已经有0人回复
-
工业萃取剂有机物含量怎么测
已经有5人回复
-
江湖救急，关于镍钼催化剂
已经有3人回复
-
光催化反应器冷点的问题
已经有0人回复
-
草酸盐沉淀-热分解法制备金属氧化物
已经有0人回复
-
酶联生物|ELISA试剂盒实验小贴士
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
- 【求助】MS写perl脚本计算 the energy of the nonbond 时出错
已经有4人回复
- 【求助】【紧急求助】关于开壳层的轨道互换
已经有9人回复
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
- 金币: 42635.7
- 红花: 153
- 帖子: 1728
- 在线: 1274.1小时
- 虫号: 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
5楼2011-05-16 11:32:20
|
一般情况下ROHF包括多个Slater行列式，它的轨道能并不是唯一的，可以有也可以没有Koopmann定理。不同的程序包可以给出不同的轨道能和相同的总能量，参见JCP 125, 204110 (2006) 比如Cl的基态，ROHF/6-31G* 5d, Gaussian 09A01给出 Alpha occ. eigenvalues -- -104.86600 -10.60162 -8.06641 -8.06641 -8.06384 Alpha occ. eigenvalues -- -1.05128 -0.48445 -0.48445 -0.19785 SCF Done: E(ROHF) = -459.442780764 A.U. after 7 cycles Convg = 0.1424D-08 -V/T = 1.9999 molpro 2010. 1 给出 CONVERGENCE THRESHOLDS: 1.00E-05 (Density) 1.00E-07 (Energy) ... !RHF STATE 1.3 Energy -459.44278076 Orbital energies: 1.1 2.1 3.1 4.1 5.1 -104.867755 -10.603340 -1.070927 0.500364 1.047729 1.2 2.2 3.2 -8.067214 -0.489407 0.608471 1.3 2.3 3.3 -8.071657 -0.568175 0.569678 1.4 1.064504 1.5 2.5 3.5 -8.067214 -0.489407 0.608471 1.6 1.115252 1.7 1.064504 欢迎重复计算 |
6楼2011-05-16 14:44:56
7楼2011-05-16 20:42:12
8楼2011-05-16 20:51:25
9楼2011-05-16 21:16:02
10楼2011-05-16 21:25:14
20