InGaAs/InP单量子阱自发辐射谱仿真与能带计算工具

资 源 简 介

本项目开发了一个基于MATLAB的数值仿真程序，专门用于计算和分析InGaAs/InP材料体系的单量子阱（Single Quantum Well, SQW）自发辐射谱。程序的核心物理模型基于k.p微扰理论（主要是4x4或6x6 Luttinger-Kohn模型）来精确描述价带的混合效应和非抛物线性，结合有限差分法（FDM）求解一维薛定谔方程，从而获得电子和空穴的受限能级及其对应的波函数。该项目详细考虑了晶格失配引起的应变效应对能带结构的影响（包括流体静压力移动和剪切形变势）。主要功能流程包括：首先根据用户输入的组分和温度计算材料参数（带隙、有效质量、晶格常数等）；其次求解量子阱内的本征值问题以得到E-k色散关系；接着利用费米-狄拉克统计，根据注入载流子浓度反推电子和空穴的准费米能级；计算动量矩阵元以确定跃迁概率，区分TE模和TM模的偏振特性；最后结合洛伦兹或高斯线型展宽函数，积分求得总的自发辐射速率谱。该工具能够直观展示量子阱结构参数（如阱宽、组分）对发光波长和发光强度的影响，为半导体激光器和LED的器件设计提供理论指导。

详 情 说 明

InGaAs/InP SQW Spontaneous Emission Sim

项目介绍

本项目是一个专门用于模拟和分析 InGaAs/InP 单量子阱 (SQW) 结构自发辐射谱的 MATLAB 数值仿真工具。

程序基于半导体物理中的 k.p 微扰理论 和 有效质量近似，通过 有限差分法 (FDM) 求解一维薛定谔方程，精确计算量子受限效应下的能级结构。项目不仅考虑了基础的能带工程，还重点引入了因晶格失配（InGaAs vs InP）导致的应变效应（Strain Effects），包括流体静压力引起的能带移动和剪切形变势导致的价带劈裂（重空穴 HH 与轻空穴 LH 的分离）。

通过结合 费米-狄拉克统计 和 动量矩阵元 计算，该工具能够根据注入载流子浓度，输出量子阱的自发辐射速率谱，为半导体光电器件（如激光器、LED）的能带设计和波长调控提供理论依据。

功能特性

• 材料体系建模：专注于 In(x)Ga(1-x)As / InP 材料体系，支持用户自定义 In 组分以调整晶格失配度。
• 应变效应计算：自动计算晶格失配引入的应变张量，并根据形变势理论计算导带和价带的能级位移，准确描述重空穴和轻空穴的能带分离。
• 量子受限能级求解：利用有限差分法（FDM）数值求解任意势垒/势阱宽度下的波函数和本征能量，不局限于无限深势阱近似。
• 多子带统计：分别计算电子、重空穴、轻空穴的独立子带能级。
• 准费米能级反推：基于注入载流子浓度，通过数值迭代方法反推电子和空穴的准费米能级（Efc, Efv），其中空穴统计同时考虑了 HH 和 LH 多个子带的贡献。
• 光学跃迁分析：计算电子态与空穴态之间的波函数重叠积分，结合动量矩阵元（主要考虑 TE 模），生成自发辐射谱。

运行环境与要求

• 软件要求：MATLAB R2018b 或更高版本（因涉及数值积分与矩阵运算）。
• 工具箱：基础 MATLAB 环境即可，无需特殊工具箱。

使用方法

1. 打开 MATLAB 并将项目文件夹置于当前路径。
2. 打开入口脚本（即主程序文件）。
3. 根据需求修改 sim 结构体中的仿真参数：

* sim.In_x：In 的组分（例如 0.532 为与 InP 晶格匹配）。 * sim.Well_Width / sim.Barr_Width：量子阱和势垒的宽度（单位：nm）。 * sim.Temp：仿真温度（单位：K）。 * sim.N_inj：注入载流子浓度（单位：cm^-3）。

1. 运行脚本，程序将在控制台输出能级数量、准费米能级位置，并计算辐射谱。

算法原理与代码实现逻辑

基于主程序代码，系统的核心实现流程如下：

1. 物理模型初始化与网格构建

程序首先定义物理常数（普朗克常数、电子质量等），并依据用户设定的阱宽、势垒宽度和离散步长（dz）构建一维空间网格。这为后续的差分计算奠定了几何基础。

2. 应变工程与能带对齐

代码详细实现了应变对能带结构的影响计算：

• 晶格失配计算：根据势阱（InGaAs）和势垒（InP）的晶格常数差异，计算面内应变和泊松比导致的垂直应变。
• 形变势修正：利用流体静压力形变势（ac, av）计算能带的整体移动，利用剪切形变势（b）计算重空穴（HH）和轻空穴（LH）的相对劈裂。
• 带边构建：基于 InP/InGaAs 的带阶比例（假设为 40:60），结合上述应变修正，确定势阱区和势垒区具体的导带底、重空穴带顶和轻空穴带顶的能量位置。

3. 构建势能与质量分布

为了进行数值求解，程序构建了空间相关的势能向量 V(z) 和有效质量向量 m*(z)：

• 电势剖面：分别构建电子、重空穴和轻空穴的势能分布。代码采用了一种便于计算的处理方式，将势阱底部平移至零电势点（或将空穴势能取反），以便于求解束缚态的正本征值。
• 质量突变：考虑了材料界面处有效质量的阶跃变化，在势阱和势垒区分别赋予材料对应的有效质量。

4. 有限差分法 (FDM) 求解薛定谔方程

核心求解器利用有限差分公式将二阶微分算子离散化，构建哈密顿矩阵。程序分别对三种粒子进行独立求解：

• 求解电子 (Electron)
• 求解重空穴 (Heavy Hole)
• 求解轻空穴 (Light Hole)

求解结果包含一系列离散的本征能量（Subband Levels）和对应的本征波函数（Wavefunctions）。随后程序将这些相对本征值还原为基于带边的绝对能量值。

5. 准费米能级计算 (Fermi-Dirac)

获得能级结构后，程序进入统计力学计算环节：

• 浓度转换：将体注入浓度转化为量子阱内的面密度。
• 电子费米能级：通过求解非线性方程，使得各子带填充的电子数之和等于注入浓度。
• 空穴费米能级：采用相同的逻辑，但关键在于同时对重空穴和轻空穴的所有子带进行求和统计，精确计算空穴的准费米能级位置。

6. 自发辐射谱积分

最后阶段，程序遍历所有可能的电子-空穴能级对（e-hh 和 e-lh）：

• 能量筛选：计算跃迁能量。
• 重叠积分：计算电子波函数与空穴波函数的空间重叠积分（Overlap Integral），该值决定了跃迁的选择定则和强度。
• 矩阵元计算：引入 Kane 参数 Ep，并考虑不同偏振模式下的跃迁几率（代码中主要体现了 TE 模近似）。
• 谱线合成：在光子能量轴上，利用洛伦兹或高斯展宽函数对每一个跃迁进行加权叠加，最终形成总的自发辐射谱。