24小时热门版块排行榜    

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

独孤笑愚

铁虫 (初入文坛)

[求助] 计算氢原子能量 已有1人参与

VASP计算单个氢原子吸附在金属表面的吸附能,计算的氢原子能量为-0.007eV,请教各位大神指点
INCAR如下:
PREC=Accurate
LWAVE=.FALSE.
LCHARG=.FALSE.
ISMEAR=0
SIGMA=0.01
NSW=0
EDIFF=1E-5
ENCUT=400
LUSE_VDW=.TRUE.
AGGAC=0.0000
IBRION=-1
谢谢!
回复此楼

» 猜你喜欢

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

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

独孤笑愚

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by hakuna at 2014-12-02 18:55:37
1.自旋要考虑
2.盒子要足够大,且去掉对称性
3.gama点

加上自旋得到的氢原子能量为-1.07eV
如果计算氢原子时加自旋,计算氢吸附在金属表面时要加自旋吗
4楼2014-12-03 09:07:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

yahxi66

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流! 2014-12-03 21:15:25
独孤笑愚: 金币+1 2015-01-08 10:41:07
The result from Turbomole at TZVPP basis set level /PBE functional. When the H has single positive charge.

Total energy =     -0.4996173348 H =       -13.5952871 eV

HOMO-LUMO Separation
   HOMO:     1. a   1 a1g   -0.27759341 H =     -7.55371 eV
   LUMO:     2. a   2 a1g   +0.21685392 H =     +5.90090 eV
   Gap :                    +0.49444733 H =    +13.45460 eV

Number of MOs=     14, Electrons=      1.00, Symmetry: oh   

   Nr.   Orbital    Occupation       Energy
   14. b   3 a1g                 +763.877043 H =    +20786.164 eV
   13. b   2 t1u                 +714.921510 H =    +19454.015 eV
   12. b   1 t2g                 +552.702424 H =    +15039.807 eV
   11. b   1 eg                  +552.702424 H =    +15039.807 eV
   10. b   2 a1g                 +364.996525 H =     +9932.066 eV
    9. b   1 t1u                 +222.262431 H =     +6048.072 eV
    8. b   1 a1g                  +20.678526 H =      +562.692 eV
    7. a   2 t1u                   +3.424311 H =       +93.180 eV
    6. a   1 eg                    +3.063591 H =       +83.365 eV
    5. a   1 t2g                   +3.063591 H =       +83.365 eV
    4. a   3 a1g                   +2.084403 H =       +56.720 eV
    3. a   1 t1u                   +0.531815 H =       +14.471 eV
    2. a   2 a1g                   +0.216854 H =        +5.901 eV
    1. a   1 a1g      1.000        -0.277593 H =        -7.554 eV
2楼2014-12-02 17:49:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

★ ★ ★
ljw4010: 金币+1, 谢谢交流! 2014-12-03 21:15:46
独孤笑愚(franch代发): 金币+2, 谢谢回帖交流。。。 2015-01-08 11:02:50
1.自旋要考虑
2.盒子要足够大,且去掉对称性
3.gama点
3楼2014-12-02 18:55:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

引用回帖:
4楼: Originally posted by 独孤笑愚 at 2014-12-03 09:07:40
加上自旋得到的氢原子能量为-1.07eV
如果计算氢原子时加自旋,计算氢吸附在金属表面时要加自旋吗...

H吸附计算加自旋终归没错,只是考虑自旋计算时间会长些,对结果有多大影响得看计算结果,一般会去测试一下。多数情况因为表面很大,会把吸附物自旋的影响给“淹灭”了,为节省时间,可以考虑不计自旋影响,慎重一点的话,建议先测试。
5楼2015-01-09 11:30:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见