目录

随机合金 InGaAs 的有效能带结构

版本:2015.1

随机合金正在成为学术上非常重要的材料。例如,SiGe 是最著名的热电材料之一,InGaAs 则是未来 CMOS(互补金属氧化物半导体)器件中替代或补充硅的主要候选材料。热电元件和晶体管器件的模拟通常都涉及到能带结构参数,如带隙和有效质量。

在本教程中,您将了解如何设置一个随机合金的结构并计算所谓的有效能带结构。本教程将聚焦于 In0.53Ga0.47As,并采用 MGGA 交换关联函数进行 DFT 计算以获得精确的带隙。

注意

在正常的能带结构中,每条能带在 E 和 k 方向上都是具有零宽度/加宽的线:应用 Bloch’s 定理,波数 k 为好量子数。在无序合金中,Bloch’s 定理不是很严格地适用,人们也不能预知 k 是否仍然是一个好量子数或者整个能带结构的概念是否定义准确。

在 QuantumATK 2015 中,可以按照参考文献 [PZ10] [PZ12] [HES11] 中描述的方法计算随机合金的“有效能带结构”。目的为研究随机合金的能带结构是否被很好地定义,以及提取能带参数是否合理。

方法

在计算有效能带结构时,人们会考虑用一个包含许多原子的较大超胞。本例中,我们采用 3×3×3 的原始 FCC 晶胞,包含 64 个原子。我们以 InAs 为 In0.53Ga0.47As 的初始模型,将其中一些铟原子由镓原子随机取代。在所有的计算中,In/Ga 的比率固定保持在 0.53/0.47。

包含超过 500 条电子能带的超胞的能带结构中,我们可以“展开”能带结构使其仅包含那些对应于初始晶胞的能带。在展开过程中,可以计算光谱权重 $\mid \langle e^{ i \bf{k} \cdot \bf{r}} \mid \psi_{j,K} \rangle \mid$,其中 $\mid \psi_{j,K} \rangle$ 为超胞在 k 点 $K$ 处的本征态。如果超胞只是单纯性单胞的简单复制(如 InAs 的超胞),光谱权重为 0 或 1,能带展开可以明确完成。然而,在无序合金的情况下,如 In0.53Ga0.47As,能带权重可以是 0 和 1 之间的任意实数。

提示

本教程使用特定版本的QuantumATK创建,因此涉及的截图和脚本参数可能与您实际使用的版本略有区别,请在学习时务必注意。

  • 不同版本的QuantumATK的py脚本可能不兼容;
  • 较新的版本输出的数据文件默认为hdf5;
  • 老版本的数据文件为nc文件,可以被新版本读取。

InAs 的能带结构

您将采用初始(1x1x1)晶胞首先计算 InAs 的正常能带结构。然后创建一个 3x3x3 的 InAs 构型,您将展示该重复超胞的有效带结构复制初始晶胞的有效带结构。

1×1×1 InAs

exchange_correlation = MGGA.TB09LDA

并更改为

exchange_correlation = MGGA.TB09LDA(c=0.91)

注意

你可以从教程 Accurate bandstructure calculations for 2D confined InAs 中了解更过关于 MGGA 交换关联势和 MGGA c参数的信息。

3×3×3 InAs

结果

1×1×1(左图)和 3×3×3(右图)InAs 构型的能带结构如下所示。3x3x3 的能带结构在区域边界处折叠,使得难以定位与 1x1x1 构型的 X 和 L 波谷最小值处相对应的波谷。

3×3×3 InAs 的有效能带结构

打开 Scripter 窗口,按照以下步骤:

脚本至此已完成,保存并运行计算。在笔记本上大约需要 3 分钟完成。

注意

在计算运行期间,以下警告会出现在 log 文件中:

Warning: The calculation did not converge to the requested tolerance!

这是完全没问题的。执行单个 SCF 迭代以获得用于有效带结构计算的初始晶胞的基函数。

分析工具显示了 k,E(k) 平面的带密度图。每个能带都将在带密度上增加 0 到 1 之间的权重。在 $\Gamma$ 点有三个简并价带,带密度变为 3(上图中为亮黄色)。导带的权重为 1。

为了查看所有权重接近 0 的能带,您可以在对数尺度下绘制能带密度。在图上单击鼠标右键,然后在菜单底部选择 Log scale,制得的图如下左图所示。很多能带都是红色的。通过检查颜色条,您会注意到这些条带的权重大约为10-9,这基本上只是数字噪声。可以尝试调整“Data range”(右上角)并将最小值设置为 -6。“噪声能带”现在已经消失了,如下右图所示。

重要

到目前为止,您已经看到有效能带结构起到了作用,并给出我们期望的结果:3x3x3 InAs 的有效能带结构与简单的 1x1x1 InAs的正常能带结构相同。

在下一节中,您将学习如何设置随机 InGaAs 合金 - 在该系统中,无法计算简单的 1x1x1 构型的正常能带结构,而且能带结构的整个概念是否被定义明确也并不明显。

随机合金 In0.53Ga0.47As

您现在将设置一个随机合金 In0.53Ga0.47As 的结构并计算它的有效能带结构。我们还应该考虑样本平均有效能带结构,分析有限能带增宽的结果。

随机合金

打开 Builder ,创建合金构型:

在新构型中,大约有 47% 的铟原子被镓原子随机取代。将此 In0.53Ga0.47As 随机合金构型发送到 Script Generator

有效能带结构计算

结果

使用 Effective Bandstructure Analyzer 研究计算的有效带结构。它应该类似于左下图。能带的加权将不再是 0 或 1,而是介于两者之间。特别是在 L- 点和 X- 点导带最小值附近,与可见的 InAs 的简单能带结构存在明显的偏差。尝试调整数据范围并将最大值设置为0.5。然后如下右图所示,一些能带变得更加明显。在 $\Gamma$- 点附近,价带和导带都是明确定义的,随机无序的影响似乎相当有限。

注意

由于随机合金的产生方式,得到的有效能带结构不一定与上图完全相同。

有效能带结构的样本平均

通过多次重复上述过程产生不同的随机构型应该可以获得数量上更可靠的有效能带结构。一些 python 脚本可以优雅高效地帮助我们完成任务。

下载 Python 脚本 ↓ingaas_ebs.py,它将创建许多随机的 In0.53Ga0.47As 合金,并执行每一个的有效能带结构计算。然后,您就可以对结果进行平均。生成合金的默认数量(number_of_runs)为 10。通过增加该值,您将获得更好的结果。

运行脚本,观察出现在 LabFloor 上文件 InGaAs_3x3x3.nc 中的 EffectiveBandstructure 数据块。标记所有生成的数据对象,然后单击 Effective Bandstructure Analyzer。分析器将显示所有有效带结构密度的平均值。可以尝试通过更改数据范围和最大值来调整颜色:

追踪具有最高权重的最重要的能带是可能实现的:如上图所示,用鼠标画一个方框。窗口底部将出现一个进度条,表示正在分析用于能带追踪的有效能带结构数据。片刻之后会出现一个新窗口,如下图所示。该图显示了用有限值绘制的追踪能带。线的中心被识别为一条能带。右键单击图片并选择 Export data 可以将数据导出为 text 文件。

提示

按以下方法计算能带:

  1. 对于每个 k,首先通过进行 Lorentzian 函数的卷积(能量的展宽为 0.2 eV)使能带密度平滑。
  2. 找到平滑后能带密度的峰值。
  3. 通过选择最接近前一个 k 点(或者,对于第一个 k 点,最接近您所绘制方框的中心)处找到的峰值能量形成能带。
  4. 继续下一个 k 点,重复该过程。

在峰值能量附近的 -0.15-0.15 eV 区间里计算宽度作为原始数据(即未平滑的)的二阶矩。最大宽度为 0.3 eV。请注意,如果两条或以上能带彼此非常靠近,由于其他能带的密度将会影响到宽度,计算的宽度将会有略微的误差。尝试在最高价带周围画个方框。能带的追踪图将会很快更新出现,如下图所示。靠近 $\Gamma$- 点,第二高的价带会错误地对第一个的宽度做出贡献。

有效增宽

与正常的能带结构相反,有效能带将有一个有限的增宽。如果增宽很小,则能带结构是明确的;如果增宽很大,则能带的概念将无效。如上图所示,随机合金 In0.53Ga0.47As 的平均有效能带结构在接近于 $\Gamma$- 点处定义明确。然而,导带在 L- 点和 X- 点的波谷最小值周围没有明确定义,该处能带有一个 $w \approx$ 0.06 eV 的有限增宽。

在 $\Gamma$- 点周围定义明确的有效能带结构可以根据此处的长波长理解。局部随机变化不会被长波长电子“看见”。相反,在 X 和 L 点,波函数在原子尺度上的变化受到了无序的影响。

有限能带增宽的分析和结果

人们可以就一个周期 $\tau = \hbar / w =$ 11 fs jieshi 解释有效能带增宽。与都从虚拟晶体近似(VCA)计算得到的近似纵向有效质量 $m_l$ = 1.5 $m_e$ 和横向有效质量 $m_t$ = 0.2 $m_e$ 一起,我们这样可以估计有效平均自由程(MFP):

$$\mathrm{MFP} = v_k \tau ,$$

在这里,能带速度为$v_k = 1/\hbar \cdot \partial E(k)/\partial k$。假设一个在波谷极小值处的抛物线能带 $E(k) = \hbar^2 k^2/2m$,相对于 X- 或 L- 波谷极小值,关于能量函数的平均自由程为:

$$\mathrm{MFP} = \frac{\hbar}{\omega}\sqrt{\frac{2E}{m}} .$$

在波谷极小值点,平均自由程消失,在 X- 或 L- 波谷最小值之上 0.01 eV 的恒定能量表面处,我们分别在横向和纵向估算得到 1.5 nm 和 0.5 nm 的平均自由程。明显地,随机无序是在这些波谷间输运的限制因素。平均自由程 vs 能量图可以借助 python 脚本 ↓ mean_free_path.py 绘制。

最终评论

对于科学用途,我们建议使用更大的超胞,如 4x4x4,甚至更大,这将更像真正的随机合金。如若选择本教程中使用的相对较小的 3x3x3 的晶胞将使计算相对较快。小的晶胞尺寸可能会以定量的方式影响结果,但我们确实希望主要结论保持不变。

由于接近于 $\Gamma$- 点的有效能带结构具有较窄的增宽,可以合理地假设布洛赫定理适用于此且随机无序对输运性质的影响很小。在本例中,有关采用虚拟晶体近似进行更有效地随机合金能带结构研究的详解可参见 Virtual Crystal Approximation for InGaAs random alloy simulations

参考