| 查看: 5823 | 回复: 12 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
niuli金虫 (小有名气)
|
[求助]
求教关于氧原子化学势和用vasp计算氧分子的化学势的问题 已有4人参与
|
|
|
Sample Text 我是一名初学者，在计算缺陷形成能时，有几个有关氧分子化学势的问题想请教大家。  1.文献里有讲 μo=1/2μo2 那么氧分子的化学势就是单个氧分子的能量吗？还是有其他的公式？ 2.如果氧分子的化学势就是单个氧分子的能量，我用vasp对氧分子结构优化和态密度计算后得到的氧分子能量作为氧分子化学势与在小木虫上查到的相差有点 多，我算的结果是-10.5ev。以下是我的输入文件，麻烦大家帮我检查检查。 赝势选择的是LDA 优化结构： INCAR SYSTEM = O2 ISTART = 0 ICHARG = 2 VOSKOWN = 1 PREC = Normal LREAL = .F. IBRION = 2 ISIF = 2 NSW = 100 POTIM = 0.5 EDIFFG = -0.05 ENCUT = 450.00 ISPIN = 2 NELM = 60 EDIFF = 1E-06 LCHARGE = .F. LWAVE = .F. ISMEAR = 1 SIGMA = 0.2 ALGO = Fast KPOINTS Auto generation 0 Gamma 1 1 1 0.0 0.0 0.0 POSCAR O2 1.00000000000000 10.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 10.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 10.0000000000000000 O 2 Direct 0.5039303366104366 0.4904031981679433 0.5158095157155261 T T T 0.4031696553895624 0.4791368018320588 0.5860504712844713 T T T 态密度 INCAR SYSTEM = O2 ISTART = 0 ICHARG = 2 VOSKOWN = 1 PREC = Normal LREAL = .F. IBRION = -1 ISIF = 2 NSW = 0 #POTIM = 0.5 EDIFFG = -0.05 ENCUT = 450.00 ISPIN = 2 NELM = 100 EDIFF = 1E-06 LCHARGE = .T. LWAVE = .T. ISMEAR = 1 SIGMA = 0.2 ALGO = Fast 希望大家指出我的问题所在，感谢大家不吝赐教。大恩不言谢。  |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
-
是这周出结果还是下周出结果？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
12楼2016-05-06 10:18:40
wangyan6930
木虫 (著名写手)
- 应助: 29 (小学生)
- 金币: 1648.4
- 散金: 3258
- 红花: 7
- 帖子: 1382
- 在线: 1040.7小时
- 虫号: 1428980
- 注册: 2011-10-06
- 性别: GG
- 专业: 基础物理学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-09-18 07:55:18
niuli: 金币+20, ★★★很有帮助, 我也用过PBE算，氧原子的化学势结果是4.9。你这样说，我就安心很多呢。 2015-09-21 16:08:23
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-09-18 07:55:18
niuli: 金币+20, ★★★很有帮助, 我也用过PBE算，氧原子的化学势结果是4.9。你这样说，我就安心很多呢。 2015-09-21 16:08:23
| 好像差距也不算太大吧，还有一个问题是这样的，因为我们计算化学势的时候是考虑能量计算的，但是系统的能量值用不同的赝势得到应该是不致的，原因是能量0点的选择不一样，所以用不同软件和计算方法计算出来应该没有可比性。当然VASP GGA-PBE计算出来O原子的大约是4.6-4.8吧，也可以我的理解有问题，仅做参考 |
2楼2015-09-17 22:06:36
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
3楼2015-09-18 04:23:35
4楼2015-09-18 22:39:13
