| 查看: 2937 | 回复: 4 | |||||
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吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 缺陷计算 | 计算模拟 | 带电体系形成能 |
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有11人回复
-
什么时候能放榜呀？
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有8人回复
-
UV压敏胶开发
已经有5人回复
-
面上提前没消息，有中的吗
已经有16人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有25人回复
-
2026年国自然面上资助率
已经有22人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有12人回复
-
微信指数没变化，科研之友没阅读
已经有19人回复
【答案】应助回帖
感谢参与，应助指数 +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不是你想要的，所以要减掉。 |
» 本帖已获得的红花（最新10朵）
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