GaAs/AlGaAs量子阱：基于八带闪锌矿结构k·p模型的教学示例
1. 引言
在本教程中，我们将借助OghmaNano所提供的八带闪锌矿结构k·p模型来计算GaAs/AlGaAs量子阱的电子与光学特性。我们首先使用一个现成的示例，查看其层结构，随后运行模拟并分析由此得到的导带与价带子带分布。最后，我们将研究基于这些能态计算出的光学光谱。
该抛物线型量子阱教程采用相对简单的有效质量带模型来处理电子与空穴。在此模型中，哈密顿量包含一个导带基组以及重空穴HH、轻空穴LH和自旋轨道分裂空穴SO这三个价带基组，每个基组又包含两种自旋分量，总计八个基态。这些基态间的耦合作用使得模型能够描述单独使用抛物线近似模型无法呈现的物理现象，尤其是价带混合现象以及平面波矢量变化时能带曲率的变化情况。
2. 打开八带模型示例
在OghmaNano中，点击“文件”功能区上的“新建模拟”按钮。双击“GaAs演示”文件夹（??）。随后从量子阱示例列表中选择“GaAs/AlGaAs量子阱八带k·p模型”（??）并双击。此时系统会自动打开一个预配置好的模拟环境，用户无需手动输入8×8哈密顿量。
3. 量子阱结构与材料参数
主窗口中显示了一个由三层构成的 AlGaAs/GaAs/AlGaAs 结构（??）。点击左侧按钮可打开层编辑器。所提供的示例结构包含两个厚度均为40纳米的 AlGaAs 势垒层，中间夹着一层厚度为10纳米的 GaAs 量子阱。
| 层 | 厚度 | 层类型 |
|---|---|---|
| 铝镓砷 | 40纳米 | 活性势垒 |
| 砷化镓（量子阱） | 10纳米 | 量子阱 |
| 铝镓砷 | 40纳米 | 有源势垒 |
GaAs的能带边缘位于周围AlGaAs的带隙范围之内，从而实现了沿层生长方向的量子限制效应。在层编辑器中，位于中间的GaAs层被明确标记为量子阱 (??)。本示例中，光学材料参数及层类型均已完成配置。
该八带求解器所需的材料参数比标量抛物线带模型更为复杂：除了能带边缘与有效质量外，描述带间耦合的参数、自旋轨道分裂能以及应变信息均需纳入哈密顿量的构建之中。您可通过“电气参数”模块及材料的代码编辑器(??)来查看与砷化镓相关的各项函数。图中所示为价带的形变势参数b与d；在本次初步运算中，您无需对这些参数做任何修改。
八带模型计算过程中有哪些变化？
抛物线模型为各能带分别指定曲率值。而k·p方法则是构建一个矩阵，该矩阵的对角线以外的元素负责将各基础态相互耦合。从示意图上看，受限态问题的处理过程如下
\[ H_{8\times8}(k_x,k_y,-i\partial_z)\,\Psi_n(z) =E_n(k_x,k_y)\,\Psi_n(z). \]
此处\(z\)代表限制方向；\(k_x\)与\(k_y\)则分别表示平行于各层的运动方向。该包络波函数\(\Psi_n\)包含八个分量，分别体现了不同基态对波函数的贡献。在本教程中，我们将重点讨论计算所得的能级值，而非完整的矩阵运算过程。关于该模型的详细阐述可参见八带闪锌矿结构理论章节。
4. 运行仿真
返回“文件”功能区并点击蓝色的“运行模拟”按钮。此时程序将构建八带哈密顿量，并在一系列面内波矢处求解该哈密顿量。您可在“终端”标签页中实时查看计算进程（??）。
示例终端界面显示了稀疏矩阵的构建过程，随后输出各本征态的能量值以及诸如 block converged after 12 iterations 之类的提示信息。对于一系列不同的面内波矢，导带态与价带态均被分别求解。请务必等待模拟运行完毕且按钮重新变为“运行模拟”状态后，方可查看计算结果。
5. 理解计算所得的能带结构
该求解器会在每个面内波矢处计算出一组离散的局域化本征态 \(k_{\parallel}\)。因此，图中呈现的是子带色散关系，而非各层间的空间带边分布。请打开“输出”选项卡并查看 bands 文件夹，或者打开 ?? 中所列出的能带汇总文件。
5.1 导带子带
导带能谱图（??该图显示了两条向上弯曲的曲线。在 \(k_{\parallel}=0\)垂直方向上的能级间隔源于不同的量子限制态所致。随着平面波矢量的增大，能量也相应升高。 \(k_{\parallel}\) 该项反映了载流子在量子阱平面内不受限制的运动情况。
在靠近导带子带极小值处，其色散关系可近似表示为
\[ E_n(k_{\parallel})\simeq E_n(0) +\frac{\hbar^2 k_{\parallel}^{2}}{2m_{n,\parallel}^{*}}. \]
在多带模型中，这仅是一种局部近似：在整个所绘制的波矢范围内，能带曲率乃至表观有效质量未必保持恒定。导带与价带间的耦合作用亦会影响这一特性。
5.2 价电子子带：混合现象与曲率特性
价带图更为耐人寻味（??）。图中所示的两个能带分支起始于不同的能量值，随后向下弯曲，似乎在某一特定的\(k_{\parallel}\)值处相交。在闪锌矿结构的价带中，重空穴态与轻空穴态并非表现为两个完全独立的抛物线形能带：在一定的波矢值下，哈密顿量中的相关项会导致这两种态的特性发生混合。
一种实用的表述方式便是考察两个相邻的能态，它们通过某个矩阵元发生耦合 \(V\)。
\[ H_{\mathrm{illustrative}}= \begin{pmatrix}E_{\mathrm{HH}}(k)&V(k)\\ V^{*}(k)&E_{\mathrm{LH}}(k)\end{pmatrix}. \]
若这些子带之间允许发生耦合，则在相互靠近时它们的能量值会发生分离，且重空穴带与轻空穴带的特性也会在不同分支间发生交换。在因对称性导致相关耦合作用消失的情况下，也可能出现真正的能量交叉现象。仅凭当前截图尚无法区分上述各种情形：为确定各子带的真实特性，需仔细分析对应的本征矢量；必要时还需在看似发生交叉的位置附近计算更多数据点。上文所示的二阶矩阵仅作示例之用，并非完整的八带哈密顿量。
5.3 同时查看导带与价带状态
该综合示意图（??）将电子态与价带电子态置于同一能量标尺上加以呈现。在此示例中，价带态的能量值接近参考零点，而受量子限域效应影响的导带态则位于1.5电子伏特及以上的能量区间。任意两个能级间的能量差值即对应着可能发生的垂直带间跃迁所需的光子能量，该值还取决于相应的光学矩阵元及载流子占据情况。
由于光子动量相对较小，一个实用的初步近似便是认为光学跃迁过程中平面波矢保持不变：
\[ \hbar\omega=E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}), \qquad \lambda=\frac{hc}{\hbar\omega}. \]
不同的电子与空穴态组合便会在不同波长处产生贡献。多带波函数还能决定各带间跃迁的强度与偏振特性，这些信息仅凭能量值无法获知。
6. 量子阱的光学特性
该八带模型计算可得出用于光学模拟的电子能量值及混合态波函数。以本例而言，计算结果分别给出了受激辐射、吸收、净增益以及自发辐射的谱线。各图表中所呈现的多条曲线分别对应模拟过程中所设定的不同载流子状态序列；这些曲线并非代表八个独立的哈密顿量能带。
6.1 受激辐射与吸收
受激辐射现象发生于入射光子促使电子从被占据的导带态跃迁至空置的价带态的过程。吸收则是相反的过程：价带中的电子吸收光子后跃迁至可使用的导带态。这两种过程的电子占据情况可示意如下
\[ \begin{aligned} \text{stimulated emission:}&\quad f_c(1-f_v),\\ \text{absorption:}&\quad f_v(1-f_c). \end{aligned} \]
受激辐射曲线（??）与吸收曲线（??）所涉及的跃迁能量值相同，但由于载流子占据率的不同，两者的强度亦存在差异。在这些计算结果中，大部分光谱特征均出现在790至850纳米范围内；具体形态则取决于所设定的子带结构、矩阵元素、载流子分布以及谱线展宽情况。
6.2 净增益：当受激辐射占主导地位时
净增益曲线（??）综合反映了受激辐射与吸收这两种相互竞争的物理过程。通过计算两种过程的占据概率之差，可得到一个极为简洁的结果：
\[ f_c(1-f_v)-f_v(1-f_c)=f_c-f_v. \]
因此，当参与跃迁的导带态电子占据率高于相应价带态时，该跃迁便会产生净增益。净增益值为负时表明吸收过程占主导地位；正值则意味着受激辐射强度超过了吸收过程。此处所示的净增益系数仅为量子阱层面的光学计算结果，不可直接视作完整激光谐振腔的模态增益，后者还受光场限制效应及其他光学损耗因素的影响。
6.3 自发辐射
自发辐射无需外部入射光子即可触发电子态间的跃迁过程。当某个被占据的电子态通过辐射方式重新结合到空置的价带态时，自发辐射便会发生。此时对应的占据因子仍为\(f_c(1-f_v)\)，不过自发辐射与受激辐射在光谱前因子及物理含义上存在差异。??中所展示的光谱采用了对数形式的垂直坐标轴；因此，垂直方向上的等间距并不意味着发射光子率的绝对值呈等幅增长。
如需详细了解跃迁强度、载流子占据率以及谱线展宽的相关内容，请参阅光学跃迁与增益理论页面。
7. 查找输出文件
bands与gain文件夹内则包含其他计算结果。模拟完成后，在主窗口中选择输出选项（??）。此时可见的结果包括用于绘制能带图的 Ec.csv、Ev.csv 和 Ec_Ev.csv，以及 gain.csv、absorption.csv、net_gain.csv 和 spontaneous.csv。此外，bands 与 gain 文件夹中还存有更多求解器输出数据。双击任意 CSV 文件即可在 OghmaNano 中打开对应图表，也可导出其中的数值数据以供其他程序分析使用。
8. 自行尝试：将本计算结果与抛物线近似模型的结果进行对比。
建议的实验方案：运行抛物线近似下的GaAs/AlGaAs示例程序，同时执行此八带模型计算。需确保GaAs量子阱宽度保持为10纳米，且所用材料参数及其他设置均保持一致。随后对比两者所得的能带分布图与光谱数据。
- 在\(k_{\parallel}=0\)区域附近，最低的导带子带是否近似呈现抛物线形态？
- 价带色散特性与基于抛物线近似的独立带模型相比有何差异？
- 计算所得的光学阈值及相对光谱特征会发生怎样的变化？
作为第二项实验，将GaAs量子阱的厚度从10纳米调整为8纳米，同时保持40纳米厚的AlGaAs势垒层保持不变。随后重新运行模型，并对比各子带的本征能量以及光跃迁所对应的波长范围。
✅ What to look for
在能带极值点附近，基于八带模型的计算结果可近似呈现为有效质量抛物线形式。远离该极值点后，能带混合与非抛物线特性便愈发明显，尤其在价带中表现显著。减小量子阱厚度通常会提升束缚态能量，促使低能量带间跃迁所对应的光子能量升高（即波长变短）；具体变化程度则取决于材料参数及耦合后的本征态特性。