24小时热门版块排行榜    

查看: 5537  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cherrily

木虫 (正式写手)

[交流] 【求助】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文件有很多不明白,问题问的会比较低一点,请大家不吝赐教,谢谢,已经发愁好几天了。
回复此楼

» 收录本帖的淘帖专辑推荐

化学电源之锂离子电池

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

cherrily

木虫 (正式写手)

引用回帖:
Originally posted by ym23 at 2010-05-22 11:45:59:
楼主问了很好的问题,我只能回答第3个问题,Binding E的数值所有原子的结合能。这个值=体系的总能减去所有单原子能量之和,再除以总原子数。

请问是不是整个体系的结合能?
包括所有原子和离子的。
那么为什么是定义成除以原子数?
谢谢
7楼2010-05-24 11:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

ym23

铁杆木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuli8(金币+1): 2010-05-22 13:08:30
楼主问了很好的问题,我只能回答第3个问题,Binding E的数值所有原子的结合能。这个值=体系的总能减去所有单原子能量之和,再除以总原子数。
博观约取、厚积薄发---壁立千仞、无欲则刚
2楼2010-05-22 11:45:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengfengyu

新虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
楼主的问题,我也不懂,同学习中,希望高手们能多多指点一下
3楼2010-05-22 19:14:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水月慧心

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
嗯 随时关注这个帖子 希望有高手可以解答 帮顶
4楼2010-05-23 10:16:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见