| 查看: 2449 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
限制开壳层计算单电能出错
|
||
| 我用WB97XD方法计算H2和一个含Fe的金属有机多孔聚合物的吸附能。由于体系具有磁性，自旋多重度是3，在计算中发现在H2距离聚合物较近的时候，自旋污染比较严重（这个我也不知道为什么，但是从结果看，距离近的时候自旋污染很大，距离远到3A以外，自选污染一下子就减小到可以忽略的程度了），因此我只好换ROWB97XD算，计算都是正常结束的，但是发现算的吸附能明显不对，由近及远，吸附能时大时小，只有两个点看起来是合理的，剩下的六个点从趋势上看明显是错的，而这两个点的输入文件和那六个点的输入文件，除氢分子与聚合物之间的距离以外完全一样。请问版上的各位高人帮忙分析一下问题出在什么地方呢？用限制开壳计算需要注意哪些问题呢？谢谢！~ |
回复此楼» 猜你喜欢
多相催化-2026年CSC博士及洪堡合作博士后招生！
已经有0人回复
-
请问现在还有电池材料方向的博导招人吗
已经有5人回复
-
物理化学论文润色/翻译怎么收费?
已经有56人回复
-
大连工业大学招收储能电池方向博士1名
已经有0人回复
-
中国地质大学（武汉）—国家级青年人才杨明教授组-招收博士-新能源材料化学及催化材料
已经有25人回复
-
《把心放慢，世界就会温柔起来》
已经有0人回复
-
求光催化产过氧化氢的PPT ！！！毕业答辩需要
已经有1人回复
-
钠离子硬碳负极扣式半电池组装都没有电流！！！
已经有2人回复
-
Ni元素XPS分析
已经有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人回复
- 【求助】【紧急求助】关于开壳层的轨道互换
已经有9人回复
|
一般情况下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
8楼2011-05-16 20:51:25
10楼2011-05-16 21:25:14