这是本文档旧的修订版!
版本:2017.2
由电子-声子相互作用产生的非弹性散射效应在确定通过器件的载流子传输中起到了基本作用。在本教程中,你将研究在反向偏压情况下,这种效应对通过硅 $p$-$n$ 结中电流的影响。在这个系统中,载流子的传输主要是由 $p$ 区价带和 $n$ 区导带之间的带间隧穿控制。非弹性散射效应可大大增强该过程。
具体来说,您将使用 ATK 中可应用的 InelasticTransmissionSpectrum 模块来计算在电子-声子相互作用的情况下器件的透射光谱。您还将使用 QuantumATK 中的 Inelastic Transmission Spectrum Analyzer 插件对电流进行详细分析,并确定哪种声子模是造成电流增强的主要原因。
QuantumATK 中采用的方法基于 Lowest Order Expansion (LOE) [FPBJ07] 法和 eXtended LOE (XLOE) [LCF + 14] 法。在本教程中,您将使用 XLOE 法,该方法处理 $\mathbf{k} \to \mathbf{k}\pm\mathbf{q}$ 过渡中初始态和最终态间的有限能量差时会由于电子-声子的耦合而更精确。然而也可以使用更简单的 LOE 法执行本教程,将会获得定性相似的结果。您可以在 ATK 手册 InelasticTransmissionSpectrum 模块的“注释”部分中找到有关这两种方法的更多详细理论信息。
为构建器件,您可以按照 Silicon p-n junction 教程中“构建器件”部分的说明进行操作。但在本教程中,您将创建一个较短的器件以加快计算速度。在 Builder,按照以下变动创建器件:
您可以在此处下载生成的器件构型:↓ Si_pn_junction.py。
在本节,您将计算分析 $p$-$n$ 结在反向偏压 $V_\mathrm{bias} = -0.4\ \mathrm{V}$ 下的电子结构,以及使用 ATK 中 TransmissionSpectrum 模块里常规 Landauer 公式获得的没有电子-声子相互作用的相关传输光谱。
点击 按钮,将结构从 Builder 发送到 Script generator。按照以下添加模块并设置参数:
最后,在 Global IO 选项,更改 Default output file 名称为 transmission_V-0.4.hdf5
。
利用 按钮将脚本发送到 Job manager,保存为 transmission_V-0.4.py
,然后按下 运行计算。使用 4 个 CPU 只需要不到 2 分钟就能完成。您可以在此处下载完整的脚本:↓ transmission_V-0.4.py。
在 LabFloor,从文件 transmission_V-0.4.hdf5
中选择 ProjectedLocalDensityOfStates 数据块,点击 Projected Local Density of States 插件。在插件窗口,从左上角的下拉菜单选择 Spectral Current,设置 Data range 的最大值为 0.1。在屏幕左侧的面板中,使用 Zoom 按钮确保电流的最小值约为 $10^{-24} \mathrm{A/eV}$。低于此值的特性可能与数字噪声相关。
从上图的左侧面板可以看出,谱电流的最大值出现在偏压窗口内($-0.2\ \mathrm{eV} \leq \mathrm{Energy} \leq 0.2\ \mathrm{eV}$)。此外,从右侧面板显示的器件态密度可以看出,在偏压窗口内,器件的中心区域(20\ \mathrm{Å} \leq \mathrm{z} \leq 45\ \mathrm{Å})没有电子态。因此,左右电极之间的隧穿是唯一可能的传输机制。这表明反向偏压下器件的电子输运受隧穿控制。
为分析在 $\mathbf{k}$ 空间隧穿的可能性。选择同一文件中 Transmission Spectrum 数据块,点击 Transmission Analyzer 插件。在插件窗口,设置 Data range 的最大值为 0.1。
从上图可以看出,费米能级处布里渊区的 $\Gamma$ 点是隧穿的隧道概率的最高峰值。
在本节,您将采用 XLOE 近似 [LCF+14] 计算硅 $p$-$n$ 结在反向偏压 $V_\mathrm{bias} = -0.4\ \mathrm{V}$ 和电子-声子相互作用下的透射谱。
为计算 InelasticTransmissionSpectrum,您首先需要计算系统的 DynamicalMatrix 和 HamiltonianDerivatives。
在 LabFloor,将文件 transmission_V-0.4.hdf5
中的 DeviceConfiguration 数据块拖拽到 Script generator。删除 DeviceSemiEmpiricalCalculator 并替换为 ForceFieldCalculator。
按照以下添加模块并设置其参数:
最后,在 Global IO 选项处,将 Default output file 的名称更改为 dynmat.hdf5
。
点击 按钮将脚本发送到 Job manager,保存为 dynmat.py
。按下 按钮运行计算。您也可以在此处下载完整的脚本:↓ dynmat.py。
或者,您也可以采用 special thermal displacement - Landauer method (STD-Landauer)将电声耦合效应包含在传输计算内,计算成本会降低很多。The STD-Landauer case study 中探讨的系统与本教程讨论的相似。
在 LabFloor,将文件 transmission_V-0.4.hdf5
中的 DeviceConfiguration 数据块拖拽到 Script generator,按照以下修改 New Calculator 模块:
添加一个 Study Objects HamiltonianDerivatives 模块,参照以下设置参数:
在 Global IO 的选项里,更改 Default output file 的名称为 dHdR_V-0.4.hdf5
。
点击 按钮将脚本发送到 Job manager,保存为 dHdR_V-0.4.py
,按下 按钮运行计算。本次计算使用 8 个 CPU 将耗时 20 分钟。您可以在此处下载完整的脚本:↓ dHdR_V-0.4.py。
在此,为加快计算速度,对哈密顿量导数做了非自洽计算。
为计算电流的非弹性部分,点击 Script generator,添加以下模块:
打开 Analysis from File 模块,选择 transmission_V-0.4.hdf5
文件,加载文件中包含的 DeviceConfiguration 数据块。
您将会注意到还添加了两个其他的模块:
在我们的示例中,动力学矩阵和哈密顿量导数已经在之前计算出来了,可以重复使用。删除这两个模块,再添加两个 Analysis from File 模块,按照下图排列:
点击从上数第二个 Analysis from File 模块。选择 dynmat.hdf5
文件,加载 DynamicalMatrix 数据块。取消勾选这个文件中包含的 DeviceConfiguration 数据块。该面板应如下图所示:
然后,点击从上数第三个 Analysis from File 模块。选择 dHdR_V-0.4.hdf5
文件,加载 HamiltonianDerivatives 数据块。该面板应如下图所示:
最后,为 Analysis InelasticTransmissionSpectrum 模块设置如下参数:
在 Global IO 选项里,更改 Default output file 的名称为 xloe_V-0.4.hdf5
。
点击 按钮将脚本发送到 Job manager,保存为 xloe_V-0.4.py
,按下 按钮运行计算。本次计算使用 24 个 CPU 将耗时一个半小时。您可以在此处下载完整的脚本:↓ xloe_V-0.4.py。
利用 Inelastic Transmission Spectrum Analyzer 可以详细地分析非弹性透射谱的结果。
在 LabFloor,选择文件 xloe_V-0.4.hdf5
的 InelasticTransmissionSpectrum 数据块,然后点击屏幕右侧插件面板中的 Inelastic Transmission Spectrum Analyzer 插件。
首先,您将分析电流对 $\mathbf{k}$ 和 $\mathbf{q}$ 的依赖性。在分析器的主窗口,设置以下参数:
您将获得的图与下图相似,展示了在取样的 k 点范围内电流对 $\mathbf{k}$ 依赖性。从图中可以明显看出,对电流的主要贡献来源于由 $k_\mathrm{A}$ 和 $k_\mathrm{B}$ 定义的二维布里渊区的 $\Gamma$ 点,也不可忽视有限 $\mathbf{k}$ 矢量的贡献。
接下来,把 Plot type 更改为 “Current vs. q-points”。您将获得与以下相似的图,展示了电流对 $\mathbf{q}$ 的依赖性。从图中可以明显看出,对电流的主要贡献来源于由 $q_\mathrm{A}$ 和 $q_\mathrm{B}$ 定义的二维布里渊区的 $\Gamma$ 点,也不可忽视有限 $\mathbf{q}$ 矢量的贡献。
上面的两个图都表明电流主要与发生在二维布里渊区 $\Gamma$ 点的传输右端,无论是 $\mathbf{k}$ 空间和 $\mathbf{q}$ 空间。为分析得出 $\Gamma$ 点的哪一种声子模贡献最大,按照下面设置分析器的主窗口:
从上图可以看出,对电流贡献最大的声子指数为 66,能量为 $\hbar \omega = 63.01\ \mathrm{meV}$。
选择在 “设置动力学矩阵计算” 部分得到的文件 dynmat.hdf5
中的 VibrationalMode 数据块,使用 Vibration Visualizer 可以将声子模可视化。
在上面显示的图中,电流值是负的,因为模拟是在反向偏压条件下执行的。