| 查看: 5152 | 回复: 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 希望大家指出我的问题所在，感谢大家不吝赐教。大恩不言谢。  |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | 科研啊科研 |
» 猜你喜欢
275学硕081000服从调剂到其他专业，保不住本专业了
已经有8人回复
-
085600，320分求调剂
已经有4人回复
-
071000生物学调剂
已经有4人回复
-
材料与化工306分找调剂
已经有3人回复
-
一志愿北京科技材料科学与工程288分，求调剂
已经有10人回复
-
一志愿同济大学323分（080500）求调剂
已经有6人回复
-
一志愿085600中科院宁波所276分求调剂
已经有26人回复
-
材料专硕322分
已经有7人回复
-
材料调剂
已经有4人回复
-
294分080500材料科学与工程求调剂
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
magicmonk
至尊木虫 (著名写手)
- 应助: 172 (高中生)
- 金币: 11451.3
- 红花: 24
- 帖子: 1124
- 在线: 3108.8小时
- 虫号: 1191315
- 注册: 2011-01-17
- 专业: 理论和计算化学
7楼2015-09-20 13:33:44
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
4楼2015-09-18 22:39:13
