太赫兹量子级联激光器中有源区上激发态电子向高能级泄漏的研究*

李金锋 万婷 王腾飞 周文辉 莘杰 陈长水2)†

1) (华南师范大学信息光电子科技学院, 广东省微纳光子功能材料与器件重点实验室,广州市特种光纤光子器件与应用重点实验室, 广州 510006)

2) (江门珠西激光智能科技有限公司, 江门 529000)

(2018 年10 月22日收到; 2018 年11 月21日收到修改稿)

利用热力学统计理论和激光器输出特性理论, 建立了太赫兹量子级联激光器(THz QCL)有源区中上激发态电子往更高能级电子态泄漏的计算模型, 以输出功率度量电子泄漏程度研究分析了晶格温度和量子阱势垒高度对电子泄漏的影响. 数值仿真结果表明, 晶格温度上升会加剧电子泄漏, 并且电子从上激发态泄漏到束缚态的数量大于泄漏到阱外连续态, 同时温度的上升也会降低激光输出功率. 增加量子阱势垒高度能抑制电子泄漏, 并且有源区量子阱结构中存在一个最优量子阱势垒高度. THz QCL经过最优量子阱势垒高度优化后, 工作温度得到提升, 其输出功率相比于以往的结果也有所提高. 研究结果对优化THz QCL有源区结构、抑制电子泄漏和改善激光器输出特性有指导作用.

1 引 言

太赫兹(terahertz, THz)波是指频率为0.1—10 THz的电磁波, 具有安全性、宽带性、指纹谱特性和穿透性等特点, 被广泛应用于医学、通信和雷达等领域[1-6]. 太赫兹量子级联激光器(terahertz quantum cascade laser, THz QCL)是一种有效获得THz波的半导体激光器[7-10], 其利用半导体异质结在外加电压的情况下, 具有量子态级联的效果, 进而声子辅助电子共振使得电子从上一级到达下一级, 外部注入的单个电子可以发射出多个光子. 然而, 存在一部分电子在运输过程中会偏离输运路径, 直接降低激光器系统的粒子数反转, 使得激光器输出功率受限, 这部分电子被称为泄漏电子. 为了更好地提高激光器输出效率, 有必要对电子泄漏进行研究.

电子泄漏主要分为三个方式: 一是电子从上激发态通过长纵声子散射到低能态; 二是电子从下激发态散射到高能束缚态和连续态; 三是电子从上激发态散射到高能束缚态和连续态. 在THz QCL中, 处于上激发态的电子向下激发态跃迁发光占据主导, 因此电子从上激发态通过长纵声子散射到低能态相对较小; 电子下激发态与高能束缚态和连续态的耦合强度低, 因此电子从下激发态散射到高能束缚态和连续态跃迁概率低; 而第三种泄漏方式中上激发态的电子是由上一周期的基态注入, 该状态下电子具有较高的不稳定性, 当温度升高时, 电子向更高能级泄漏概率大, 因此在高温的环境下, 电子从上激发态散射到高能束缚态和连续态是电子泄漏的主要途径. 目前, 大多数研究者都是通过间接测量来解释电子泄漏[11-16], 通过建立理论模型来分析的研究较少.

本文利用热力学理论结合激光器输出特性, 建立电子从上激发态散射到高能束缚态和连续态的电子泄漏计算模型, 重点研究该模型中晶格温度和量子阱势垒高度对电子泄漏的影响, 进而找到抑制电子泄漏的有效措施来提高激光器输出功率.

2 理论模型

本节利用如图1所示的THz QCL能级结构[17]研究电子从上激发态泄漏到高能束缚态和阱外连续态的泄漏方式. 图中1, 2, 3, 4分别表示基态、3态的孪生态、下激发态、上激发态, 5和6表示近邻束缚态和次近邻束缚态, 7表示阱外连续态,1′表示上一周期的基态, 4′′表示下一周期的上激发态. 该能级结构由材料系统建造的量子阱形成, 发射的激光频率3.9 THz,h为普朗克常量,Eph=15.6 meV为光子能量; 有源区材料的声子能量ELO=36.8 meV;能级参数39 meV分别表示7态与4态、6态与4态、5态与4态之间的能量差;J45,J46,J47分别表示电子从4能态泄漏到5, 6, 7能态的泄漏电流密度. 该激光器系统是通过优化振子强度, 并采用斜跃迁形式发光, 能形成更大的粒子数反转, 具有一定的典型性. 在THz波的产生过程中, 电子从上一周期的1′态弛豫到达4态, 然后从4态跃迁到3态发射出THz光子, 接着3态和2态上的电子通过声子辅助共振到达1态, 随后进入下一级的 4′′态. 在此过程中, 4态上的电子泄漏到5态、6态和7态造成电子泄漏, 导致电子利用率降低, 不利于激光器的粒子数反转, 进而影响激光器输出功率. 我们用激光器的输出功率值衡量电子泄漏程度, 半导体激光器输出功率Pout可以描述为[18]

图1 THz QCL级联的能级结构Fig.1. Energy level structure of a single THz QCL cascade.

式中hυ为发射光子的能量;q为电子电荷;αm和αw分别为镜面损失和波导损失;A为接触面积;η为内量子效应;J和Jth分别为注入电流密度和阈值电流密度.

THz QCL发射的光子能量通常在2—40 meV之间[19], 这远比中红外QCL要小, 这种结构上的差异使得THz QCL中4态波函数与5态、6态和7态的波函数耦合更强, 4态上的电子不稳定性更大, 造成4态上的电子向更高能级泄漏成为主要的泄漏部分. 在激光的阈值条件下, 我们主要分析电子从4态泄漏到5态、6态和7态的泄漏模型,(1)式改为

式中Jesc为总泄漏电流密度, 由4态电子泄漏到5态、6态和7态的泄漏电流密度J45 ,J46 ,J47 相加得到, 即在没有电子泄漏且电子势能被完全转化为光能的情况下, 激光器的最大输出功率为为了更加清晰地看出电子泄漏对激光器输出功率的影响, 将(2)式进行归一化:

由(3)式可知, 归一化的输出功率与泄漏电流密度J45,J46,J47有关, 接下来将对这三个参数进行分析. 先对4态上的电子泄漏到5态和6态对应的泄漏电流密度J45,J46进行分析. 设j为高能束缚态(j=4,5 ), 当温度升高时, 4态上的电子能吸收足够的能量跃迁到子带j, 在有限的能量空间中, 这种散射机制形成的电流密度可以表示为

式中n4为4态上的电子浓度;分别为4态上电子泄漏到5态和6态的时间;kB为玻尔兹曼常数.

接下来对4态上的电子泄漏到7态对应的泄漏通电流密度J47进行分析. 由电子在能态上占据的时间与泄漏到该能态的电流密度之间的相互关系得

式中τ47为电子从4态散射到7态的时间. 电子的电荷q乘以4态上的电子浓度n4等于4态上的电荷浓度, 1 /τ47表示电子从4态跃迁到7态的概率,因此(7)式的物理意义为4态上的电子泄漏到7态对应的泄漏电流密度. 它是一个与温度有关的函数, 该函数关系可由热激发模型得到[22]

将(5)式, (6)式和(9)式代入(3)式可以得到电子泄漏模型下的归一化输出功率表达式

由(10)式可知, 激光器的输出功率与晶格温度和量子阱势垒高度密切相关, 在接下来的工作中将详细讨论这两个参数对激光器输出功率的影响.在激光系统中, 7态位于量子阱口, 为激光器系统最高的势能态, 我们用7态与4态之间的能量差E47作为量子阱势垒参数, 并在第3节中对参数E47最优值进行讨论与分析.

3 结果与讨论

本节将依据第2节建立的理论模型研究晶格温度和量子阱势垒高度对THz QCL电子泄漏的影响. 选取的仿真分析参数如下: 激光器发射频率为3.9 THz, 最高运行温度为186 K, 限制因子Tc=0.86[23], 电子-晶格耦合常数散射时间注入电流密度J=1000 A/cm2[17,24], 电子浓度等于掺杂浓度

经过第2节的研究可发现, 晶格温度是影响电子泄漏的一个重要参数, 为解释电子泄漏的内在机制, 首先分析泄漏电流密度和晶格温度TL的关系,其关系曲线如图2所示. 由图2可知, 在186 K的晶格温度下, 电子从4态泄漏到5态、6态和7态的总泄漏电流密度为810 A/cm2, 模拟数值与文献报道的实验值吻合[17,24]. 从图2还可以看出, 当晶格温度从80 K增加到200 K时, 电子泄漏随着温度的增加而不断增大, 且在晶格温度增加的过程中电子主要是从4态泄漏到5态, 还可以看出电子从4态泄漏到5态和7态对温度的变化比较敏感,而泄漏到6态的数目相对较少. 出现该现象的原因是4态和5态之间的能级差等于39 meV, 该能级差与THz QCL有源区选用的GaAs材料的辅助共振声子能量接近, 在这种环境下, 4态上的电子容易吸收一个声子的能量而跃迁到5态[25]. 由能级结构可知, 6态处于较高能级, 在量子阱中能级越高, 量子阱对电子所施加的势能就越低, 因此电子被该能态束缚的概率就越低, 当电子向量子阱上部泄漏时, 泄漏到量子阱较高束缚态的概率相对较低, 出现了电子泄漏到6态的数量少的现象.

图2 泄漏电流密度与晶格温度的关系Fig.2. Relationship between leakage current density and temperature.

然后研究在不同的泄漏通道下, 晶格温度TL的变化对输出功率Pout/Pmax的影响, 结果如图3所示, 红色、绿色和蓝色曲线分别表示4态电子泄漏到5, 6态(束缚态), 7态(连续态)以及5, 6, 7(束缚态以及连续态)态对输出功率的影响. 在考虑不同的泄漏情况下, 只把J45,J46表达式代入(3)式,不考虑J47值影响, 便可得到4态电子泄漏到5,6态对功率的影响; 同理, 只将J47表达式代入(3)式, 不考虑J45与J46值的影响, 可得到 4 态电子泄漏到7态对功率的影响; 将J45 ,J46 与J47 同时代入(3)式, 可得到4态电子泄漏到5, 6, 7态对功率的影响. 从图3可以看出, 随着晶格温度的升高, 归一化输出功率将不断降低, 而且电子从4态泄漏到5, 6态对输出功率的影响要比泄漏到7态显著, 由此也可反映出晶格温度的升高将加剧电子泄漏. 另外, 处于5, 6态和7态的电子不是稳定状态的电子, 随着温晶格度的增加, 电子热运动加剧,5, 6, 7态的电子有一部分继续向更高能级跃迁或是逃离有源区, 在这些过程中, 电子因运动与晶格产生碰撞, 致使晶格温度升高, 进一步恶化系统的温度特性和输出功率特性.

图3 归一化输出功率与晶格温度的关系Fig.3. Relationship between normalized output power and lattice temperature.

最后研究量子阱势垒高度参数E47对电子泄漏的影响. 图4模拟了不同晶格温度条件下输出功率随势垒参数E47的变化关系. 由仿真结果可知,在186 K温度附近,E47最优值为130 meV, 该值与报道的原始结构值吻合[17,24], 进一步证明了本文计算模型的正确性. 从图4还得出, 在势垒高度较小的情况下, 归一化输出功率也相对较小; 当势垒在某一区间内增加时, 归一化输出功率呈现线性递增; 当势垒高度增加到某一临界值时, 归一化输出功率趋于稳定. 这是因为势垒高度过低时, 电子受到的束缚力低, 电子泄漏的数目多, 造成输出功率相对较低. 当势垒在某一区间内增加, 归一化输出功率呈现线性递增, 这说明电子泄漏得到了抑制,电子的利用率得到提高, 因此输出功率上升. 通过解量子阱有源区中的薛定谔方程, 可以得到能级的解, 薛定谔方程为

式中E为能态的解;ψ为对应能态的波函数;m∗(z)为有效质量,z为量子阱生长方向;V(z) 为量子阱势能参数;F(z) 为施加的电压参数. 当施加电压不变时, 量子阱势垒高度增加后, 通过数值求解薛定谔方程发现量子阱中的束缚态数目增多, 而从图3讨论的结果可知, 电子从4态泄漏到5, 6态(阱内束缚态)对激光器造成的负面影响比电子从4态泄漏到7态要大. 因此, 需要计算出一个最优的量子阱势垒高度, 该值应满足量子阱中的束缚态数目少、此时的势垒高度能有效地对电子泄漏进行抑制. 由图4可知, 当势垒高度增加到一定值时,归一化输出功率趋于稳定值, 定义当有激光输出时, 归一化输出功率首次到达这个稳定值时对应的势垒参数即为所求. 对文献[17]报道的THz QCL进行势垒高度的优化, 如图4所示, 将晶格温度设置在210 K下, 模拟得出了对应的最优阱高参数E47= 185 meV, 同时用(2)式模拟得出该条件下的输出功率值为8 mW, 优化后取得的效果与最高186 K操作温度下实验报道的输出功率5 mW相比有所改善.

图4 归一化输出功率与势垒高度参数的关系Fig.4. Relationship between normalized output power and barrier height parameters.

4 结 论

利用热力学统计理论和激光输出特性理论, 建立了THz QCL有源区中上激发态电子往更高能级电子态泄漏的计算模型, 以输出功率来衡量电子泄漏程度, 研究了晶格温度和量子阱势垒高度对电子泄漏的影响. 研究发现, 晶格温度上升, 电子从上激发态泄漏到近邻束缚态和连续态的数量增加相对较多, 泄漏到次近邻束缚态的数量相对较小.同时, 电子泄漏导致电子向外部散发热量, 进一步加剧温度升高, 形成温度-电子泄漏的恶性循环, 相应地, THz QCL的输出功率会因温度的升高而降低, 影响THz QCL的正常运行. 适当增加量子阱势垒高度能够有效地抑制电子泄漏, 进而改善输出功率. 本文对THz QCL系统进行有源区势垒优化, 当量子阱势垒高度提高到185 meV时, 升高温度至210 K能得到THz QCL激光输出功率为8 mW, 改善了以往在最高186 K操作温度下得到的输出功率. 这些研究结果对研究THz QCL的电子泄漏温度特性和THz QCL有源区结构优化设计提供理论依据.