24小时热门版块排行榜    

查看: 7401  |  回复: 12

bingmou

金虫 (著名写手)

[交流] 【求助】空位形成能计算的讨论 已有7人参与

到底应该怎么计算空位形成能呢?

stractor在这里(http://muchong.com/bbs/viewthread.php?tid=1728484&fpage=0&view=&highlight=&page=2)曾经说过:

这样计算形成能: 设E1为完整晶格的能量,E2为有一个空位的晶格的能量,E3为挖走这个Al原子的能量,这形成能 E=(E2+E3)-E1. 对E3的处理不同研究者会给出不同的能量,有的是给出FccAl晶格中一个Al原子的能量,有的给出是原子状态的能量,有的还会考虑其它环境因素而给出另外的E3表达式(Phys. Rev. B 80, 184110 (2009))。

而这个E3的确有很多中处理方法,例如第三种其实就是这里(http://www.materialssimulation.com/node/258)所采用的化学势的方法;第二种方法就是侯老师在他的手册中提到的采用单原子的方法,用VASP计算的同学都应该很熟悉。第一种是比较特殊的,采用所谓的最稳定晶体结构的方法。

但是我觉得这后面两个种都有点问题

首先空位形成应该是这样的一个过程:晶体中的空位是在不断形成又同时消失的,原子有可能跳入间隙位而产生了空位,也有可能是位错而产生了空位。所以这个少掉的原子还在这体系中,而不是回到单质中(第一种情况),因为这时候已经不是在生成反应中了;也不是回到真空基态(第二种情况),那是对表面原子才有可能的。因此我觉得我们在计算空位形成能的时候必须采用第一种方法,确定这个离开位置的原子的化学势

但是怎么确定化学势呢?

有人认为应该采用E_F=E_{defect}-\frac{N-1}{N}E_{no defect}的方式,也就是认为\frac{E_{no defect}}{N}是这个原子的化学势。如果对于单个成分的体系来说,这样是可以的;但是如果是对于GaN这种体系还可以吗?在Ga和N这两种成分的化学势相等吗?我觉得应该是可以的,因为一个体系稳定的要求就是两种元素的化学势是相等的,所以这样的公式是比较可靠的。这个公式来源:ab initio calculation of the formation energy of charged vacancies in germanium (感谢enola,http://muchong.com/bbs/viewthread.php?tid=2440380)

不知各位对此有什么看法?是否同意我的看法呢?

这里讨论的是自发缺陷,而非摻杂缺陷。摻杂缺陷我是认同采用单质(第一种情况)的做法的+适当考虑环境的因素的。

[ Last edited by bingmou on 2010-10-20 at 15:55 ]
回复此楼

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

研究生资料集 VASP 第一原理资料汇编 vasp计算资料
带电掺杂体系形成能 模拟知识汇总 杂七杂八 MS学习
我的分享

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bingmou

金虫 (著名写手)

★ ★
aylayl08(金币+2):谢谢讨论 2010-10-22 15:12:49
在黄昆《固体物理学》(1988年10月第一版 2005年11月第21次印刷)中第544页有:
对\omgea(形成能)的含义,容易产生误解,认为它代表吧一个在格点上的原子从晶体中拿走所需要的能量。实际上,根据上面的讨论,n改变时原子书目N没有变,也就是说,并没有吧原子拿出晶体以外,所以应该把\omega理解为将晶格内部一个格点上的原子放到晶体表面所需要的能量,如图12-17所示。

这里说的移到表面的能量应该怎么计算呢?
6楼2010-10-21 10:28:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hgye

金虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢参与热点问题讨论 2010-10-21 10:10:21
采用化学势计算是合理的,用单质态确定E3只是化学势的一个极限取值。在GaN这样的化合物体系中,原子化学势是有一个变化范围的,不能简单认为两个组份的化学势相等。详见J. Appl. Phys. 95,3851
2楼2010-10-20 20:39:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

引用回帖:
Originally posted by hgye at 2010-10-20 20:39:09:
采用化学势计算是合理的,用单质态确定E3只是化学势的一个极限取值。在GaN这样的化合物体系中,原子化学势是有一个变化范围的,不能简单认为两个组份的化学势相等。详见J. Appl. Phys. 95,3851

如果不能默认为相等,那么这时候Ga的形成能计算化学势应该怎么计算呢?这时候是否应该采用单质的情况给出范围,然后画曲线呢?

[ Last edited by bingmou on 2010-10-20 at 23:39 ]
3楼2010-10-20 21:22:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
bingmou(金币+2):感谢详细回复 2010-10-20 23:36:48
zzy870720z(金币+2):谢谢参与热点讨论 2010-10-21 10:10:49
先补充一点吧。
“首先空位形成应该是这样的一个过程:晶体中的空位是在不断形成又同时消失的,原子有可能跳入间隙位而产生了空位,也有可能是位错而产生了空位。所以这个少掉的原子还在这体系中,而不是回到单质中(第一种情况)……”。
这种情况不是孤立的空位了,而是缺陷形成的complex。bingmou你提到的这种空位-间隙对是Frenkel 缺陷。当然从研究角度讲,如果对complex感兴趣,可以研究一下。

考虑化学势就是因为缺陷的产生受环境的影响的,比如,间隙或者空位会与其它元素结合生产其它相,会对化学势有一定的限制。

至于哪些条件来限制化学势,这很复杂,一般要考虑各种可能的相,得出一些化学势的可能值的范围或者定出几个特定的点。
我即将做半导体ZnO的研究,在这方面看了一些文献。推荐几篇,大家参考一下:
1. Choi et al., "First-principles study of native defects and lanthanum impurities in NaTaO3", PHYSICAL REVIEW B 78, 014115 (2008).
2. Chris G. Van de Walle and J. Neugebauer, "First-principles calculations for defects and impurities: Applications to III-nitrides", J. Appl. Phys. 95, 3851 (2004).

[ Last edited by stractor on 2010-10-20 at 23:30 ]
4楼2010-10-20 23:09:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

★ ★
aylayl08(金币+2):谢谢探讨 2010-10-22 15:12:40
哦,你不说我都忘了还有Schottky缺陷就是刚好对应第一和第二种情况的缺陷了,谢谢!对于空位也还仔细考虑阿
谢谢你的文献,hgye推荐的文章就是你的第二篇,非常不错!
所谓特殊的点就是那几个交界点是吗?我以前也计算过缺陷的形成能,今天再计算的时候突然想到这个问题,对空位的形成重新考虑了一下

[ Last edited by bingmou on 2010-10-21 at 09:35 ]
5楼2010-10-20 23:36:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xx4751

铜虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送鲜花一朵
在一个晶体中放入一个原子的形成能是怎么算啊,放入两个原子或两个以上呢。谢谢
7楼2012-06-07 16:55:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

引用回帖:
2400078楼: Originally posted by xx4751 at 2012-06-07 16:55:59
在一个晶体中放入一个原子的形成能是怎么算啊,放入两个原子或两个以上呢。谢谢

一般是在晶体的空隙位放入一个原子,然后计算能量差

» 本帖已获得的红花(最新10朵)

8楼2012-06-08 13:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xx4751

铜虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送鲜花一朵
引用回帖:
2400186楼: Originally posted by bingmou at 2012-06-08 13:04:03
一般是在晶体的空隙位放入一个原子,然后计算能量差...

公式是什么样的,一个和两个(或以上)
9楼2012-06-08 13:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
有空位又有掺杂的话,形成能又该如何结算呢?第二是形成能和结合能又有什么区别呢
希望有一天成为大神
10楼2013-11-18 17:04:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bingmou 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见