24小时热门版块排行榜

返回列表

niuli

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1805.5
散金: 103
红花: 1
帖子: 219
在线: 121.6小时
虫号: 484506
注册: 2007-12-27
性别: MM
专业: 凝聚态物性 II ：电子结构

[求助] 求教关于氧原子化学势和用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

      希望大家指出我的问题所在，感谢大家不吝赐教。大恩不言谢。

回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档

科研啊科研

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
论文投稿，期刊推荐已经有4人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼2015-09-17 18:38:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

wangyan6930

木虫 (著名写手)

应助: 29 (小学生)
金币: 1639.4
散金: 3251
红花: 7
帖子: 1381
在线: 1039.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

好像差距也不算太大吧，还有一个问题是这样的，因为我们计算化学势的时候是考虑能量计算的，但是系统的能量值用不同的赝势得到应该是不致的，原因是能量0点的选择不一样，所以用不同软件和计算方法计算出来应该没有可比性。当然VASP GGA-PBE计算出来O原子的大约是4.6-4.8吧，也可以我的理解有问题，仅做参考