| 查看: 4938 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】新手求助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朵）
枫魔馨» 猜你喜欢
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
11楼2015-04-03 15:49:17
2楼2010-10-11 12:05:10
lono75
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6963.1
- 散金: 100
- 红花: 2
- 帖子: 2433
- 在线: 536.4小时
- 虫号: 932458
- 注册: 2009-12-23
- 专业: 凝聚态物性 II ：电子结构
3楼2010-10-11 12:15:42
|
哦，发现了，的确写错了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
5