| 查看: 5155 | 回复: 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 希望大家指出我的问题所在，感谢大家不吝赐教。大恩不言谢。  |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | 科研啊科研 |
» 猜你喜欢
调剂
已经有6人回复
-
280求调剂
已经有6人回复
-
302求调剂一志愿华中师范大学
已经有5人回复
-
346求调剂
已经有4人回复
-
材料专硕322分
已经有9人回复
-
372分材料与化工（085600）一志愿湖南大学求调剂
已经有5人回复
-
318求调剂，计算材料方向
已经有10人回复
-
310求调剂
已经有14人回复
-
315求调剂
已经有3人回复
-
265求调剂
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
4楼2015-09-18 22:39:13
wangyan6930
木虫 (著名写手)
- 应助: 29 (小学生)
- 金币: 1638.4
- 散金: 3258
- 红花: 7
- 帖子: 1382
- 在线: 1040.3小时
- 虫号: 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
hy1987007
木虫 (著名写手)
- 应助: 12 (小学生)
- 金币: 20203.7
- 散金: 2
- 红花: 2
- 帖子: 2555
- 在线: 459小时
- 虫号: 1560286
- 注册: 2012-01-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
niuli5楼2015-09-19 14:38:41
