24小时热门版块排行榜    

查看: 3391  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +3 曼111 2026-03-24 4/200 2026-03-24 14:45 by JourneyLucky
[考研] 307求调剂 +4 超级伊昂大王 2026-03-24 4/200 2026-03-24 14:31 by JourneyLucky
[考研] 293求调剂 +6 加一一九 2026-03-24 6/300 2026-03-24 14:29 by JourneyLucky
[考研] 一志愿南航材料专317分求调剂 +4 炸呀炸呀炸薯条 2026-03-23 4/200 2026-03-24 07:32 by wangy0907
[考研] 298-一志愿中国农业大学-求调剂 +11 手机用户 2026-03-17 12/600 2026-03-23 23:51 by 热情沙漠
[考研] 327求调剂 +5 prayer13 2026-03-23 5/250 2026-03-23 22:11 by 星空星月
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +5 niko- 2026-03-22 5/250 2026-03-23 22:01 by fuyu_
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
信息提示
请填处理意见