| 查看: 1814 | 回复: 4 | |||||
| 【悬赏金币】回答本帖问题，作者xiaomao_16将赠送您 10 个金币 | |||||
xiaomao_16银虫 (小有名气)
|
[求助]
CASTEP计算带电空位缺陷形成能的问题已有1人参与
|
||||
|
如题，从文献中看到有计算公式如下：E=E(defect)-E(perfect)+μ(O)+qEf+ΔFvib(T)；其中μ(O)=E(O)-0.5E(bind_O2)+Δμ(O)；有几个问题请教一下大家： 1. Ef费米能级一般不是取为VBM=零点了吗，那这一项是要考虑为零吗，最后一项说是缺陷形成过程中的振动能，也不知道怎么算，文章里下文完全没提。。可以忽略吗？ 2. E(O2)文章取的实验值没啥说的，E(O)说是电子结构计算得来的，在CASTEP里就是构建一个大立方，放一个氧原子算能量吗？ 3. Δμ(O)怎么确定，文章说Δμ(O)=0定义富氧极限，但氧空位产生不该是缺氧情况吗？ 4. 考虑到带电情况还有修正公式：Ecorr=E(isolated)-E(periodic)+qΔV，最后一项是对齐用的。有个问题是等式右侧不就是我们要求的数吗，那这个修正的意义何在？在CASTEP里计算荷电情况的时候是直接在计算的时候Charge一栏填一个电荷数就ok吗？ |
回复此楼» 收录本帖的淘贴专辑推荐
 缺陷计算 | 计算模拟 | 带电体系形成能 |
» 猜你喜欢
看到大家一直在说代表作NC，不明觉厉
已经有13人回复
-
弱关系
已经有10人回复
-
最浪漫的事
已经有6人回复
-
郑雅婷这种将来估计也会举报自己工作单位领导
已经有17人回复
-
郑雅婷工作后被辞退无业中,于是报复硕导,证据不足失败后可能报复父母
已经有19人回复
-
前些天开会有个人见到人就搞关系，一查此人全是MDPI/Hindawi论文，鄙视!
已经有48人回复
-
蹲一个6月7号江苏省青基公示
已经有21人回复
-
函评意见：既然你有中科院****人才项目了，就不要申请面上了。
已经有20人回复
-
求生物学的博导收留
已经有7人回复
-
等待会评 等待公示
已经有7人回复
【答案】应助回帖
感谢参与，应助指数 +1
|
你可以多读一些paper和tutorial多理解一些吧，这些问题太细太多了！人家回答打字也挺麻烦的。 1. Ef费米能级一般不是取为VBM=零点了吗，那这一项是要考虑为零吗，最后一项说是缺陷形成过程中的振动能，也不知道怎么算，文章里下文完全没提。。可以忽略吗？ Ef 取值一般从VBM到CBM，Ef不同，formation energy不同。 2. E(O2)文章取的实验值没啥说的，E(O)说是电子结构计算得来的，在CASTEP里就是构建一个大立方，放一个氧原子算能量吗？ 是的。 3. Δμ(O)怎么确定，文章说Δμ(O)=0定义富氧极限，但氧空位产生不该是缺氧情况吗？ Δμ(O)很难确定，只能取两个极限值。Δμ(O)=0为富氧环境（O2 单质情况下Δμ(O)=0）。还有一种是O2 deficiency情况，Δμ(O)很小（为一较负值）。氧空位在富氧情况下确实更难形成但这里不是你理解的那样。 4. 考虑到带电情况还有修正公式：Ecorr=E(isolated)-E(periodic)+qΔV，最后一项是对齐用的。有个问题是等式右侧不就是我们要求的数吗，那这个修正的意义何在？在CASTEP里计算荷电情况的时候是直接在计算的时候Charge一栏填一个电荷数就ok吗？ qΔV 因为缺陷带电荷引起的静电势能变化。E(isolated)-E(periodic) 所有DFT计算都采用PBC即周期性边界条件，你计算的是一个单胞，但实际上CASTEP考虑的是无限大周期重复的整个晶体。所以一个缺陷和他旁边相邻的单胞内的缺陷会有相互作用，而这种相互作用是artifacts不是你想要的，所以要减掉。 |
2楼2020-03-19 04:19:05
xiaomao_16
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 373.6
- 散金: 2
- 帖子: 127
- 在线: 23.7小时
- 虫号: 21418430
- 注册: 2020-03-10
- 专业: 凝聚态物性 II ：电子结构
送红花一朵|
感谢回复，主要我不是搞这个的，有些东西一时半会儿搞不明白。我想问的主要是最后一个，E(isolated)我问题里写了，文章说是“用静电理论很容易解出来”，没用DFT 发自小木虫IOS客户端 |
3楼2020-03-19 07:20:09
4楼2020-03-21 01:42:00
5楼2020-03-21 01:57:57
30