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文件有很多不明白，问题问的会比较低一点，请大家不吝赐教，谢谢，已经发愁好几天了。

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

1楼2010-05-21 15:22:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

381233029

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 274.1
散金: 20
帖子: 17
在线: 33.1小时
虫号: 2412176
注册: 2013-04-11
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

9楼: Originally posted by study163 at 2010-06-10 09:48:01
1) 分子轨道波函数由各个原子的波函数线性组合而成，这里的值是对应的分子轨道波函数的本征值
2)你说的这个值是体系的总能
3）结合能已经有人回答了
4）这个数据里的charge说明了这个原子的电荷得失情况，正的 ...

Charge只能表示带了多少电荷，不能反映电荷得失的能力，这个是我帖子上的一位同仁说的，关于这点我也很迷茫，正在努力寻找答案