24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4177  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jttjtl

银虫 (小有名气)

[交流] 【求助】新手求助vasp计算氧分子结合能的问题 已有4人参与

本人刚学vasp,最近用vasp计算了氧原子和氧分子的总能,然后再计算了一下结合能,大小为6.48eV,比实验上的5.2eV大不少,想请大家帮忙看看我的计算有没有问题~
以下是我计算总能时的INCAR文件:
ENCUT = 400
GGA = PE
ISMEAR = 0
ISIF = 0
NSW = 0
EDIFF = 1E-05
LREAL = .FALSE.
LDIAG=.TRUE.
PREC=Accurate
NELMDL=5
ISPIN=2
NUPDOWN=2

以下是氧原子的POSCAR:
SYSTEM
1.000000
15.000000000000000       0.0000000000000000       0.0000000000000000
0.0000000000000000       15.000000000000000       0.0000000000000000
0.0000000000000000       0.0000000000000000       15.000000000000000
O
1
Direct
   7.5000000000000000   7.5000000000000000   7.500000000000000

氧分子直接是取了课本中的健长1.208A,然后放到15*15*15的盒子里面的;
氧分子的POSCAR:
1.000000
15.00000000000000      0.000000000000000      0.000000000000000
0.000000000000000      15.00000000000000      0.000000000000000
0.000000000000000      0.000000000000000      15.00000000000000
O
2
Direct
   7.5000000000000000   7.5000000000000000   7.5000000000000000
   8.5805000000000007   7.5000000000000000   7.5000000000000000

KPOINTS均为:
KPOINTS
0
gamma
1  1  1
0. 0. 0.

计算出来的结果是:
氧原子的总能:
  free energy    TOTEN  =        -1.67175867 eV

  energy without entropy =       -1.51749635  energy(sigma->0) =       -1.59462751

氧分子的总能:

  free energy    TOTEN  =        -9.82983769 eV

  energy without entropy =       -9.82983769  energy(sigma->0) =       -9.82983769
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jttjtl

银虫 (小有名气)

谢谢zxzj05的提示,我换了HSE06来计算,跟实验值非常接近了,5.18eV.
9楼2010-10-12 08:38:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

vaspnew

金虫 (小有名气)


jttjtl(金币+2): 2010-10-11 15:17:16
zxzj05(金币+1):奖励 2010-10-11 20:45:52
计算氧分子的时候应该让原子驰豫吧
2楼2010-10-11 12:05:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)


jttjtl(金币+5):谢谢提醒~ 2010-10-11 15:15:53
zxzj05(金币+1):奖励 2010-10-11 20:46:01
坐标错了,
D 为分数坐标,
由天命,成人事!
3楼2010-10-11 12:15:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jttjtl

银虫 (小有名气)

引用回帖:
Originally posted by lono75 at 2010-10-11 12:15:42:
坐标错了,
D 为分数坐标,

哦,发现了,的确写错了sweat:~
不过貌似不影响结果,这是氧原子的OUTCAR中的原子位置
position of ions in fractional coordinates (direct lattice)
   0.50000000  0.50000000  0.50000000
   0.58050000  0.50000000  0.50000000

position of ions in cartesian coordinates  (Angst):
   7.50000000  7.50000000  7.50000000
   8.70750000  7.50000000  7.50000000:
4楼2010-10-11 15:15:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕求调剂 +10 hyl3153942 2026-03-29 10/500 2026-03-31 16:31 by hypershenger
[考研] 343求调剂 +8 爱羁绊 2026-03-28 8/400 2026-03-31 16:12 by 不吃魚的貓
[考研] 085601 材料工程 313分 求调剂 +6 Ong3 2026-03-27 6/300 2026-03-31 15:18 by yanflower7133
[考研] 求收留 +8 1943443204 2026-03-28 8/400 2026-03-31 15:00 by -迷了路啊路
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +13 辰沐5211314 2026-03-26 13/650 2026-03-31 14:37 by 记事本2026
[考研] 085600材料与化工调剂 +16 kikiki7 2026-03-30 16/800 2026-03-31 10:03 by 氯化亚硝酰
[考研] 材料专硕 085600求调剂 +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
[考研] 332求调剂 +6 Lyy930824@ 2026-03-29 6/300 2026-03-30 10:53 by 1172367218
[考研] 一志愿中南大学化学0703总分337求调剂 +6 niko- 2026-03-27 6/300 2026-03-30 10:25 by herarysara
[考研] 求调剂 +10 张zz111 2026-03-27 11/550 2026-03-30 09:17 by 无际的草原
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 一志愿华理,数一英一285求A区调剂 +8 AZMK 2026-03-25 12/600 2026-03-28 18:15 by AZMK
[有机交流] 高温高压反应求助 10+4 chibby 2026-03-25 4/200 2026-03-27 21:08 by BT20230424
[考研] 复试调剂,一志愿南农083200食品科学与工程 +5 XQTJZ 2026-03-26 5/250 2026-03-27 14:49 by 狂炫麦当当
[考研] 085600,材料与化工321分调剂 +4 大馋小子 2026-03-27 6/300 2026-03-27 14:11 by 松花缸1201
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
[考研] 0703化学338求调剂! +6 Zuhui0306 2026-03-26 7/350 2026-03-27 10:35 by shangxh
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 341求调剂 +7 青柠檬1 2026-03-26 7/350 2026-03-27 00:19 by wxiongid
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
信息提示
请填处理意见