24小时热门版块排行榜    

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

hpfamily1988

木虫 (著名写手)

[求助] 氢气能量计算求助 已有2人参与

我查了一下氢气的能量大概在-6.7eV左右,我结构优化然后静态计算结果为:
free  energy   TOTEN  =         0.837274 eV
结果相差很大,不知道哪儿出了问题,麻烦大家帮我看看,多谢!
下面是我计算的文件:
结构优化:
INCAR:
SYSTEM = H2 dimer in a box
ISMEAR = 0
SIGMA = 0.01
PREC = Accurate
EDIFF  = 1E-5
ISPIN = 2
IBRION = 2
ENCUT = 500
ISTART = 0
ICHARG = 2
ISIF = 2
ISYM = 0
LWAVE = .FALSE.

KPOINTS:
Automatic mesh
0
G
1 1 1
0 0 0

POSCAR:
H2 molecule in a box
  1.0
  10 0.0 0.0
  0.0 11 0.0
  0.0 0.0 12
2
cart
  0 0 0
  0 0 1

静态计算:
INCAR:
SYSTEM = H2 dimer in a box
ISMEAR = 0
SIGMA = 0.01
PREC = Accurate
EDIFF  = 1E-5
ISPIN = 2
IBRION = 2
ENCUT = 500
ISTART = 0
ICHARG = 2
ISYM = 0
LWAVE = .FALSE.

KPOINTS:
Automatic mesh
0
G
1 1 1
0 0 0

POSCAR:
H2 molecule in a box                    
   1.00000000000000     
    10.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000   11.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000   12.0000000000000000
   H
   2
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.0833333333333357

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

hy1987007

木虫 (著名写手)

【答案】应助回帖

我计算的结果为  free  energy   TOTEN  =  -6.77073721 eV
不对的地方请指正!

SYSTEM = test
PREC   = A
ENCUT  = 500
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA  = 0.2
IBRION = 2
ISIF   = 2
NSW    = 500
EDIFFG = -0.01
LCHARG = F
LWAVE  = F
GGA    = PE

Au
0
G
1 1 1
0 0 0


H2
  1.00
15.0     0      0
0        15.0   0
0        0      15.0
2
S
C
0      0     0 T T T
0.75   0     0 T T T
11楼2016-04-08 14:04:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

hpfamily1988

木虫 (著名写手)

有没有人一起交流交流。。。
2楼2016-04-06 23:13:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvchaojia

金虫 (小有名气)

不懂帮顶,我才哈,不同的温压条件和结构下,H2的能力应该是不同的。
再说,这个能量只是个相对量,能直接和文献里的数值对比吗? 个人观点。

» 本帖已获得的红花(最新10朵)

3楼2016-04-07 10:49:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xyj880506

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
hpfamily1988: 金币+20, ★★★很有帮助, 多谢交流 很有帮助 赞! 2016-04-07 22:56:19
hpfamily1988: 金币+40, ★★★★★最佳答案 2016-04-08 00:11:05
氢气优化时必须用ISPIN=1,否则优化完就会自动解离。这是vasp的一个bug
4楼2016-04-07 11:34:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YLHlRHNwYkce 2026-09-02 7/350 2026-09-03 17:02 by T0rGB46095mJ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 10/500 2026-09-03 17:00 by T0rGB46095mJ
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 rFEsUBKXRll0 2026-09-02 6/300 2026-09-03 16:49 by T0rGB46095mJ
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 16:38 by T0rGB46095mJ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 5/250 2026-09-03 16:38 by T0rGB46095mJ
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 14:48 by T0rGB46095mJ
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:03 by T0rGB46095mJ
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +8 gdd2018 2026-08-28 13/650 2026-09-03 06:22 by 羽毛枫f
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 4/200 2026-09-03 03:14 by rM1TE0WVDIIY
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
请填处理意见