| 查看: 5874 | 回复: 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 希望大家指出我的问题所在，感谢大家不吝赐教。大恩不言谢。  |
回复此楼» 猜你喜欢
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
13楼2016-09-12 23:55:54
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
10