24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cherrily

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6364.8
散金: 80
帖子: 402
在线: 202.4小时
虫号: 416393
注册: 2007-06-30
性别: MM
专业: 无机纳米化学

[交流] 【求助】Dmol3计算outmol分析电子密度，谢谢！已有8人参与

几何优化模型之后，我计算了轨道，得到结果outmol文件
第一个问题：
里面显示如下：
Fermi Energy 0.15390Ha    4.188eV -kTS_e= -0.014745Ha
  HOMO is orbital number 40
  LUMO is orbital number 41

state                      eigenvalue       occupation
                           (au)       (ev)

1 - 1 a       -18.607907    -506.347    1.000
2 + 1 a       -18.607857    -506.346    1.000
3 - 2 a       -18.607726    -506.342    1.000
4 + 2 a       -18.607658    -506.340    1.000
。
。
。
34 +  17 a       -0.079579       -2.165    1.000
35 -  18 a       -0.079369       -2.160    1.000
36 +  18 a       -0.079285       -2.157    1.000
37 -  19 a       -0.077259       -2.102    1.000
38 +  19 a       -0.077139       -2.099    1.000
39 +  20 a       -0.076044       -2.069    1.000
40 -  20 a       -0.075971       -2.067    1.000
41 -  21 a       0.161169       4.386    0.189
42 +  21 a       0.161405       4.392    0.182
43 -  22 a       0.162940       4.434    0.141
44 +  22 a       0.163105       4.438    0.137
我查看了以前的老贴子，说电子密度就是轨道的平方，在import之后我是可以得到HOMO的图的，这里显示的特征值是不是就是轨道函数呢？

第二个问题：
      Total E (au)    Binding E (au)    Convergence    Time (m) Iter
-----------------------------------------------------------------------------
Ef    -305.5343502    -1.2624257       0.83E-04       0.623    26

binding energy  -1.262426Ha  -34.35236eV -792.200kcal/mol -3314.564kJ/mol
这个Ef的能量是不是就是费米能级的能量？那我理解为HOMO的能量，对不对呢？

第三个问题：
上一个问题的数据中，Binding E的数值是费米能级的什么结合能？

第四个问题：
计算布居结果如下
Mulliken population analysis

Mulliken atomic charges:
      charge spin
  O(  1)  -0.592 0.001
  H(  2) 0.181  -0.004
  H(  3) 0.174  -0.004
  O(  4)  -0.595  -0.001
  H(  5) 0.162 0.005
  H(  6) 0.172 0.005

请问这个spin的数值能说明什么问题？这个例子里面的数据有说明了什么?

第五个问题：
Plotting output:
  property:                         file name:
  DMol3 total electron density                               4H2O-_density.grd
  orbital 40  20 -a    -0.0760  1.00 4H2O-_homo.grd
  orbital 41  21 -a    0.1614  0.18 4H2O-_lumo.grd
这个里面显示的 -0.0760和0.1614 和第一个问题里面的轨道的特征值很相近，是不是就是特征值？如果是，为什么有一点差别？这是同一次收敛计算的结果啊。
如果不是，那是什么呢？

刚刚学习MS，对于outmol文件有很多不明白，问题问的会比较低一点，请大家不吝赐教，谢谢，已经发愁好几天了。

回复此楼