这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
atk:带电点缺陷的形成能_手动方法 [2019/10/22 22:18] – [GaAs 的带电 As 空位] xie.congwei | atk:带电点缺陷的形成能_手动方法 [2019/10/30 11:25] (当前版本) – [带电点缺陷的形成能(手动方法)] dong.dong | ||
---|---|---|---|
行 24: | 行 24: | ||
<WRAP center important 100%> | <WRAP center important 100%> | ||
=== 注意 === | === 注意 === | ||
- | 在 O-2018.6 版本中,我们引入了一个所谓的 //Study Object// 将本教程中描述的工作流程自动化。现在是 QuantumATK 中研究带电点缺陷的推荐方法,如果您可以访问 | + | 在 O-2018.6 版本中,我们引入了一个所谓的 //Study Object// 将本教程中描述的工作流程自动化。现在是 QuantumATK 中研究带电点缺陷的推荐方法,如果您可以使用 |
</ | </ | ||
行 128: | 行 128: | ||
第一项与我们刚刚对中性点缺陷的计算类似,而第二项和第三项与之前完全相同,不需要重新计算。价带最大值由无缺陷块体的态密度(DOS)计算确定,$\mu_e$ 将被视为自由参数。最后一项 $E_{corr}$ 根据 Freysoldt 等人 <color # | 第一项与我们刚刚对中性点缺陷的计算类似,而第二项和第三项与之前完全相同,不需要重新计算。价带最大值由无缺陷块体的态密度(DOS)计算确定,$\mu_e$ 将被视为自由参数。最后一项 $E_{corr}$ 根据 Freysoldt 等人 <color # | ||
==== 带电缺陷 ==== | ==== 带电缺陷 ==== | ||
+ | |||
+ | |||
+ | As 空位总共可以维持 5 个电荷态,包括上文研究过的中性状态。我们现在将计算这些点缺陷的形成能。对于非常大的系统,使用弛豫过的中性点缺陷作为电荷态的初始猜测可能是值得的,但对于本例使用的系统尺寸来说并不重要。 | ||
+ | |||
+ | 可以按照上述步骤在计算器中添加电荷,也可以直接复制脚本并在计算器模块中添加 '' | ||
+ | |||
+ | <WRAP center important 100%> | ||
+ | === 注意 === | ||
+ | 谨记在脚本 [[https:// | ||
+ | </ | ||
==== 形成能 ==== | ==== 形成能 ==== | ||
+ | |||
+ | |||
+ | 现在我们已经准备好了所有需要的部分,可以通过从保存的数据文件中提取总能量和计算修正项得到完整的形成能。为此,您可以使用脚本 [[https:// | ||
+ | 您还需要如教程:[[https:// | ||
+ | |||
+ | //表 7 含有 As 空位的 2x2x2 GaAs 超胞形成能(包括周期性和能带偏移校正)。// | ||
+ | |||
+ | ^ 电荷态 | ||
+ | | +1 | ||
+ | | 0 | 3.22 eV | 0.00 eV | 0.00 eV | | ||
+ | | -1 | ||
+ | | -2 | ||
+ | | -3 | ||
+ | |||
+ | |||
+ | 上面给出的形成能是根据 VBM 处的电子化学势计算的,即 $\mu_e = 0$ eV。考虑到 $\mu_e$ 作为自由参数能使费米能级移位,我们也可以为每个缺陷画出形成能与电子化学势的函数关系图。这也相当于随材料的掺杂水平而变化。结果如下图所示,在给定 $\mu_e$ 下,最稳定的电荷态即最低的线。由图可知,对于 p 型掺杂的 GaAs,在 $\mu_e = 0$ eV 时,带正电荷的 As 空位是最稳定的。随着电子化学势的逐渐增大,最稳定的电荷态从 $\mu_e > 0.35$ 开始变为中性,然后是负电荷态。尤其是本征(无掺杂) GaAs 在 $\mu_e = E_{gap}/ | ||
+ | |||
+ | |||
+ | {{ : | ||
+ | |||
+ | |||
+ | 图97 GaAs 中 As 空位各种电荷态的形成能随 $\mu_e$ 函数的变化(相对于 VBM 的定义),包含 VBM 为零的情况,计算采用了 2x2x2 正交晶系超胞。 | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
==== 体系大小的影响 ==== | ==== 体系大小的影响 ==== | ||
+ | |||
+ | |||
+ | 如果我们将超胞尺寸增加到 3x3x3,k 点网格更改为 2x2x2,我们会得到下表列出的形成能,以及使用相同大小超胞的的文献值< | ||
+ | |||
+ | //表 8 计算形成能(包括周期性和能带偏移校正)与文献值的比较// | ||
+ | |||
+ | ^ 电荷态 | ||
+ | | +1 | ||
+ | | 0 | 3.22 eV | 3.26 eV | 3.25 eV | | ||
+ | | -1 | ||
+ | | -2 | ||
+ | | -3 | ||
+ | |||
+ | |||
+ | 我们可以看到形成能有怎样的细微变化,中性点缺陷的偏移约为 0.04 eV。带电缺陷的差异略大,这是由于 k 点密度的小更改引起 VBM 的轻微偏移。本例没有与参考文献完全一致的,但应该注意的是,参考文献中没有使用鬼原子,采用了较陈旧的静电校正方法和 LDA 交换关联。 | ||
+ | |||
+ | 这些变化也会影响形成能和 $\mu_e$ 的函数关系图,尽管线的斜率相同,但截距略有不同,导致缺陷稳定性的整体图像发生显着变化。具体而言,中性状态现在仅在非常窄的区域稳定,且从 $\mu_e = 0.4$ 起带负电的状态最稳定。 | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | 图98 GaAs 中 As 空位各种电荷态的形成能随 $\mu_e$ 函数的变化(相对于 VBM 的定义),包含 VBM 为零的情况,计算采用了 3x3x3 正交晶系超胞。 | ||
+ | |||
+ | |||
+ | |||
+ | |||
===== 附录 ===== | ===== 附录 ===== | ||
行 137: | 行 206: | ||
==== 周期校正方法 ==== | ==== 周期校正方法 ==== | ||
+ | 正如上文已经提到的,我们会采用由 Freysoldt 等人 <color # | ||
- | ===== 参考 ===== | + | $$E_{corr} |
+ | $E_{lat}$ 修正了电荷周期图像间的静电相互作用,$q \Delta V_{q/b}$ 则修正了能带偏移。$E_{lat}$ 通过考虑模型电荷分布的静电作用和对比周期晶胞中边界条件为非周期性的势计算得到: | ||
+ | |||
+ | $$E_{lat} = \int_\Omega \left[\frac{1}{2} q (\tilde{V}_q^{lr}-V_q^{lr}) \right] d^3 r .$$ | ||
+ | |||
+ | $V_q^{lr}$ 是模型电荷分布的长程势,$\tilde{V}_q^{lr}$ 周期晶胞的修正量。$\Delta V_{q/b}$ 由引入点缺陷造成的静电势移位计算得来: | ||
+ | |||
+ | $$ \Delta V_{q/b} = \tilde{V}_{q/ | ||
+ | |||
+ | 这里的 $\tilde{V}_{q/ | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | 图 99 能带偏移校正图。绿色是缺陷的模型电位,红色是由于引入缺陷引起的变化,而黑色则是绝对能量零点的偏移。 | ||
+ | |||
+ | 注意,所获得的差异是在两个屏蔽势间定义的,参见文献 <color # | ||
+ | |||
+ | |||
+ | |||
+ | ==== GaAs 的介电常数 ==== | ||
+ | |||
+ | 为了说明半导体材料(围绕带电缺陷)对能带偏移校正的屏蔽影响,应该在脚本 [[https:// | ||
+ | |||
+ | 将优化的 GaAs 块体结构拖放到 **Builder** 并发送到 **Script Generator**。添加一个 //New Calculator// | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | 运行计算,应该不会超过几分钟。您现在应该在 **LabFloor** 上看到 // | ||
+ | |||
+ | |||
+ | {{ : | ||
+ | |||
+ | |||
+ | |||
+ | 本例中,我们只对上边那张图的截距感兴趣,看起来大约为 13。要获得精确计算的值,请打开带有 Text Representation 的 // | ||
+ | |||
+ | |||
+ | {{ : | ||
+ | |||
+ | ===== 参考 ===== | ||
+ | * 英文原文:https:// | ||
+ | * [EMM05] Fedwa El-Mellouhi and Normand Mousseau. Self-vacancies in gallium arsenide: An ab initio calculation. //Phys. Rev. B//, 71(12): | ||
+ | * [FGH+14] Christoph Freysoldt, Blazej Grabowski, Tilmann Hickel, Jörg Neugebauer, Georg Kresse, Anderson Janotti, and Chris G. Van De Walle. First-principles calculations for point defects in solids. //Rev. Mod. Phys//., 86(1): | ||
+ | * [FNVanDWalle09] (1, | ||
+ | * [PRC+08] J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. //Phys. Rev. Lett//., 100(13): | ||