| 查看: 5159 | 回复: 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 希望大家指出我的问题所在，感谢大家不吝赐教。大恩不言谢。  |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | 科研啊科研 |
» 猜你喜欢
070305高分子化学与物理 304分求调剂
已经有16人回复
-
esi高被引论文是不是能对中标有所加分和帮助呢
已经有5人回复
-
调剂求助
已经有8人回复
-
315分 085602 求调剂
已经有10人回复
-
求调剂，一志愿郑州大学材料与化工专硕，英二数二342分，求老师收留
已经有7人回复
-
302求调剂
已经有6人回复
-
一志愿北京科技大学材料工程085601，求调剂
已经有6人回复
-
材料与化工306分找调剂
已经有15人回复
-
070300化学求调剂
已经有9人回复
-
一志愿深大085601材料工程专业（专硕）300分可以调剂去哪
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
3楼2015-09-18 04:23:35
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
4楼2015-09-18 22:39:13
hy1987007
木虫 (著名写手)
- 应助: 12 (小学生)
- 金币: 20203.7
- 散金: 2
- 红花: 2
- 帖子: 2555
- 在线: 459小时
- 虫号: 1560286
- 注册: 2012-01-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
niuli5楼2015-09-19 14:38:41
