应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求教关于氧原子化学势和用vasp计算氧分子的化学势的问题
131/2返回列表12下一页
查看: 5776  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

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

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

第一性原理相关文档 科研啊科研

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2015-09-17 18:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wangyan6930

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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吧，也可以我的理解有问题，仅做参考
赞一下(1人) 回复此楼
2楼2015-09-17 22:06:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微6凉

铁虫 (小有名气)
引用回帖:
3楼: Originally posted by stractor at 2015-09-18 04:23:35
根据化学势的定义，O2的化学式就是单个氧分子的能量。除了DFT能量外（与xc的选取相关），还与与温度，压强相关。

您好，我想问一下气体化学势具体的公式是什么呢？和温度压强有关的。我现在算一个东西需要这个公式，若是能给出参考文献就更好了。谢谢
赞一下(1人) 回复此楼
12楼2016-05-06 10:18:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

stractor

金虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-09-18 07:55:25
根据化学势的定义，O2的化学式就是单个氧分子的能量。除了DFT能量外（与xc的选取相关），还与与温度，压强相关。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
3楼2015-09-18 04:23:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aisareru

新虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
氧分子的化学势是氧分子的吉布斯自由能G=U-TS+PV。DFT基态计算得到的只是U中的静态能量部分，U中还含有振动能量。另外，不同的泛函对氧分子化学势的计算影响是不小的。
赞一下 回复此楼
4楼2015-09-18 22:39:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hy1987007

木虫 (著名写手)
晶格参数再取大点，两个原子之间的距离取实验值试试
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

niuli
5楼2015-09-19 14:38:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

【答案】应助回帖

感谢参与，应助指数 +1
注意氧分子的基态是三重态。你虽然加了自旋，但是收敛到的还是单重态的，你可以自己看占据数。
赞一下 回复此楼
6楼2015-09-20 13:19:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)
引用回帖:
6楼: Originally posted by 锐利的碎片 at 2015-09-20 13:19:47
注意氧分子的基态是三重态。你虽然加了自旋，但是收敛到的还是单重态的，你可以自己看占据数。

不是可以通过NUPDOWN设置三重态吗
赞一下 回复此楼
我们最深的恐惧不是我们能力不够，我们最深的恐惧是能力超越了极限。
7楼2015-09-20 13:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher
引用回帖:
7楼: Originally posted by magicmonk at 2015-09-20 13:33:44
不是可以通过NUPDOWN设置三重态吗...

是啊，但是他没设。
回复此楼
8楼2015-09-20 14:32:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuli

金虫 (小有名气)
送红花一朵
引用回帖:
5楼: Originally posted by hy1987007 at 2015-09-19 14:38:41
晶格参数再取大点，两个原子之间的距离取实验值试试

好的好的，我再算算
赞一下 回复此楼
9楼2015-09-21 16:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuli

金虫 (小有名气)
引用回帖:
6楼: Originally posted by 锐利的碎片 at 2015-09-20 13:19:47
注意氧分子的基态是三重态。你虽然加了自旋，但是收敛到的还是单重态的，你可以自己看占据数。

怎样设置三重态呢
赞一下 回复此楼
10楼2015-09-21 16:12:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 niuli 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 关于代码变化问题，想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭