24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4172  |  回复: 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
回复此楼

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

jttjtl

银虫 (小有名气)

引用回帖:
Originally posted by vaspnew at 2010-10-11 12:05:10:
计算氧分子的时候应该让原子驰豫吧

计算总能应该不用驰豫吧?
5楼2010-10-11 15:17:02
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上5B能上会吗? +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 286求调剂 +5 丢掉懒惰 2026-03-27 8/400 2026-03-31 11:27 by Delta2012
[考研] 313求调剂 +6 卖个关子吧 2026-03-31 6/300 2026-03-31 10:58 by Jaylen.
[考研] 求调剂 +4 图鉴212 2026-03-30 4/200 2026-03-31 10:20 by cal0306
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +10 dongfang59 2026-03-27 10/500 2026-03-30 23:42 by 果果妈咪
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 085602化工求调剂(331分) +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 285求调剂 +6 AZMK 2026-03-29 9/450 2026-03-30 21:02 by dophin1985
[考研] 286分调剂 +10 Faune 2026-03-30 11/550 2026-03-30 20:35 by 啊李999
[考研] 材料专硕 085600求调剂 +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
[考研] 0703 化学 求调剂,一志愿山东大学 342 分 +7 Shern—- 2026-03-28 7/350 2026-03-30 16:31 by nothing投稿中
[硕博家园] 求调剂 有机化学考研356分 +10 Nadiums 2026-03-25 11/550 2026-03-30 11:56 by yjolah
[考研] 求调剂 +4 QiMing7 2026-03-25 5/250 2026-03-29 21:10 by 唐沐儿
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 321求调剂 +7 璞玉~~ 2026-03-25 8/400 2026-03-29 06:41 by 544594351
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-03-26 6/300 2026-03-29 01:11 by hanserlol
[考研] 0703化学求调剂,各位老师看看我!!! +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 281求调剂 +6 Koxui 2026-03-24 7/350 2026-03-26 15:37 by 无际的草原
信息提示
请填处理意见