24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jttjtl

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 420.9
散金: 62
红花: 1
帖子: 106
在线: 18.9小时
虫号: 608704
注册: 2008-09-21
专业: 半导体晶体与薄膜材料

[交流] 【求助】新手求助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