[资源] 【原创】如何计算中性杂质的形成能 已有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软件

» 猜你喜欢

• 求助:我三月中下旬出站，青基依托单位怎么办？ 已经有9人回复
• Cas 72-43-5需要30g，定制合成，能接单的留言 已经有8人回复
• 北京211副教授，35岁，想重新出发，去国外做博后，怎么样？ 已经有8人回复
• 磺酰氟产物，毕不了业了！ 已经有5人回复
• 论文终于录用啦！满足毕业条件了 已经有25人回复
• 2026年机械制造与材料应用国际会议 （ICMMMA 2026） 已经有3人回复
• 自荐读博 已经有3人回复
• 不自信的我 已经有5人回复
• 投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球 已经有8人回复

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

• 掺杂化学势及形成能计算 已经有6人回复
• g_energy里面的一些能量项 已经有8人回复
• 如何根据LD50值确定杂质限度？ 已经有3人回复
• 掺杂后形成能的计算 已经有35人回复
• 求教一个小白问题 已经有10人回复
• 如何求元素的化学势 已经有6人回复
• 怎样设置有缺陷体系电荷问题 已经有18人回复
• 掺杂的超晶胞体浓度怎么算？？ 已经有16人回复
• 材料计算中的不等价掺杂问题 已经有10人回复
• vasp计算掺杂体系的形成能 已经有7人回复
• 【求助】化学势的计算相关 已经有13人回复
• 【求助】杂质的百分含量怎么计算啊 详细点，最好举例子 已经有5人回复
• 【其他】缺陷形成能与电中性要求 已经有12人回复
• 【求助】体心立方结构间隙问题 已经有6人回复
• 【求助】重金求计算ZnO或GaN等材料方面形成能，结合能，激活能的计算过程 已经有29人回复
• 【分享】石蜡切片固定，大家认为中性甲醛和10%福尔马林哪个更好？ 已经有6人回复