24小时热门版块排行榜    

查看: 19437  |  回复: 66
【奖励】 本帖被评价33次,作者valenhou001增加金币 27.15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 【原创】如何计算中性杂质的形成能 已有3人参与

该贴是在回复“硅晶体内参杂一个铁原子那种情况下最稳定”的问题(http://muchong.com/bbs/viewthread.php?tid=1914787&fpage=0&view=&highlight=&page=1)基础上综合而成的。也就是以硅晶体中掺杂一个Fe原子为例,来介绍如何计算杂质形成能。


将晶胞扩大可能至少需要到64个原子的超原胞,然后根据Si里面的对称性、以及可能的空隙位考虑计算如下的几种情况:
1)将一个Si替换为Fe;
2)将一个Fe放到如下图的三种可能的间隙位置上:


计算步骤:
1)计算完整晶格的晶格常数;
2)构造超原胞,超原胞的大小要大到尽量使得杂质之间的相互作用足够小能足以说明物理问题,以及所能承受的计算资源。
3)计算完整晶格的超原胞的总能;计算参数特别是k点的选取,对在计算完整晶格的晶格常数时所取的k点网格,按超原胞同初基原胞的比例来缩小。
4)优化各种掺杂的原子结构,静态计算出总能。
5)按如下形成能的公式来计算杂质形成能。注意形成能的值,与会掺杂元素的化学势的选取有关的。

Si中Fe杂质的形成能:
delta E(formation) = E_tot(Si带有Fe杂质的超原胞) - E_tot(Si未带有Fe杂质的超原胞,即完整晶格的超原胞) -  n(Fe) * \mu(Fe) +   n(Si) * \mu(Si)

E_tot(...)就是超原胞的总能,有无杂质时,超原胞的大小要一致。
\mu(Fe)是Fe的化学势, \mu(Si)是Si的化学势。n(Fe)是掺杂入超原胞中Fe原子的个数,n(Si)是从超原胞中拿出的Si原子个数。

如果是替位掺杂:n(Fe) =1, n(Si) =1,表示是从超原胞拿出一个Si,然后并放入一个Fe。
如果是间隙掺杂: n(Fe) = 1, n(Si)=0,则表示只有将一个Fe原子放入超原胞,没有从超原胞拿出Si原子。

\mu(Fe)可以取为体心立方(bcc)Fe的总能(单位:能量单位+每原子);
\mu(Si)可以取为金刚石结构Si的总能(单位:能量单位+每原子。)。

[ Last edited by valenhou001 on 2010-5-22 at 07:54 ]
回复此楼

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

仿真建模与计算 第一性原理计算资源总汇 研究生资料集 VASP and MS
第一性原理 VASP 第一原理资料汇编 第一性原理相关文档
Materials-studio vasp计算资料 第一性原理计算经验 VASP使用技巧
MS第一性原理 精华帖-好资源 模拟知识汇总 first principle
MS资源 培训、总结 @个人收集@ 第一性原理杂 Photochemistry
VASP计算资料收集byZF 纳米及第一性原理方向 量子理论见解 杂七杂八
掺杂计算 vasp计算总结 VASP相关主题 MS软件

» 猜你喜欢

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

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

awmc2008

至尊木虫 (文坛精英)


★★★★★ 五星级,优秀推荐

好资料,学习学习,很需要。谢谢楼主。
27楼2010-06-24 23:18:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 67 个回答

freer007

木虫 (正式写手)



小木虫(金币+0.5):给个红包,谢谢回帖交流
多谢楼主分享。
3楼2010-05-22 10:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)



小木虫(金币+0.5):给个红包,谢谢回帖交流
谢 谢, 设为评价贴

[ Last edited by lono75 on 2010-5-22 at 10:21 ]
4楼2010-05-22 10:09:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amynliou

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

写得很好啊  尤其喜欢那三个对称位置  收藏了
7楼2010-05-24 16:52:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 03:26 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-30 01:58 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:49 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:03 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见