24小时热门版块排行榜    

查看: 3842  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖
普通回帖

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的回帖

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的回帖

jttjtl

银虫 (小有名气)

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

计算总能应该不用驰豫吧?
5楼2010-10-11 15:17:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)


zxzj05(金币+1):奖励 2010-10-11 20:46:13
要弛豫,不然得不到能量的基态。计算出的键长与实验值总是有偏差
心之所向无惧无悔
6楼2010-10-11 15:42:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


zxzj05(金币+1):奖励 2010-10-11 20:46:39
jttjtl(金币+3): 2010-10-20 10:14:37
必须要弛豫。
7楼2010-10-11 16:41:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vaspstudy

银虫 (小有名气)


jttjtl(金币+5): 2010-10-11 19:15:24
zxzj05(金币+1):奖励 2010-10-11 20:46:26
jttjtl(金币+5): 2010-10-12 08:37:23
引用回帖:
Originally posted by jttjtl at 2010-10-11 09:38:50:
本人刚学vasp,最近用vasp计算了氧原子和氧分子的总能,然后再计算了一下结合能,大小为6.48eV,比实验上的5.2eV大不少,想请大家帮忙看看我的计算有没有问题~
以下是我计算总能时的INCAR文件:
ENCUT = 400
...

下面是个人一些观点,权作交流:
(1)计算单原子能量时不需要relax,计算分子能量时最好relax;
(2)据个人试验,设置相同时,不同的赝势得到的结果可能出入比较大,如果与实验值相差太大,建议换种赝势试一下
8楼2010-10-11 17:53:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jttjtl

银虫 (小有名气)

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

cuicui2397

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
8楼: Originally posted by vaspstudy at 2010-10-11 17:53:41
下面是个人一些观点,权作交流:
(1)计算单原子能量时不需要relax,计算分子能量时最好relax;
(2)据个人试验,设置相同时,不同的赝势得到的结果可能出入比较大,如果与实验值相差太大,建议换种赝势试一下...

想请教一下:
(1)不进行relax,如何计算能量?难道是只进行一步静态自洽计算的得到的能量,就是单原子的能量?
(2)计算体系的能量是看free energy TOTEN的值,还是看energy without entropy的值?
(3)所谓体系能量是relax后的能量,还是先relax再静态自洽计算后的能量?
有点不确定以上问题,请不吝赐教。
10楼2014-04-10 10:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jttjtl 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见