应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】新手求助vasp计算氧分子结合能的问题
51/1返回列表
查看: 4888  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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朵）

枫魔馨

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-11 09:38:50
已阅   回复此楼   关注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的回帖
查看全部 13 个回答

vaspnew

金虫 (小有名气)
★
jttjtl(金币+2): 2010-10-11 15:17:16
zxzj05(金币+1):奖励 2010-10-11 20:45:52
计算氧分子的时候应该让原子驰豫吧
赞一下(1人) 回复此楼
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 为分数坐标,
赞一下(2人) 回复此楼
由天命,成人事！
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的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 青基如何破局 +5 lyz123lyz 2026-08-27 6/300 2026-08-27 17:08 by Cathyau
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +7 Kittylucky 2026-08-27 7/350 2026-08-27 16:31 by msuhlx
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
关闭
请填处理意见
关闭