超快半导体激光器中的有效拉比频率：自启动谐波频率梳

《Light-Science & Applications》：Effective Rabi frequency in ultrafast semiconductor lasers: self-starting harmonic frequency combs

　　

光学频率梳已成为光学与光子学中的关键研究课题。一种特殊的梳态是谐波频率梳（HFC），其中光学线间隔为腔自由光谱范围（FSR）的整数倍。最近在具有快速增益恢复的半导体激光器中，如量子级联激光器（QCLs），观察到HFC的自发形成，尽管其潜在物理机制尚不清楚。在这项工作中，研究人员基于有效拉比频率（ERF）与激光腔模式之间的共振现象，为QCLs中HFC的形成提供了物理解释。这通过用于描述多模激光器动力学的有效半导体麦克斯韦-布洛赫方程（ESMBEs）的数值积分结果以及连续波（CW）发射的线性稳定性分析（LSA）得到证实。

**研究背景与问题**
光学频率梳（OFC）因其在计量学、光谱学和光通信中的精密测量作用而成为光学与光子学的核心课题。近年来，片上集成梳源的发展催生了两种重要器件：光泵浦克尔微谐振器和快速半导体激光器，如量子级联激光器（QCL）和量子点激光器（QDL）。QCL作为超快半导体激光器，其载流子寿命短（皮秒量级），且具有非零线宽增强因子（LEF，α），已在中红外和太赫兹（THz）区域实现频率梳生成。然而，QCL中一种特殊梳态——自启动谐波频率梳（HFC），其光学线间隔为腔自由光谱范围（FSR）的整数倍——虽在2016年首次被观测到，但物理机制长期未明。早期研究基于法布里-珀罗（FP）激光器的空间烧孔（SHB）效应，但仅适用于近阈值条件，且无法解释高泵浦下的HFC现象。因此，揭示HFC的自发形成机制成为该领域的关键未解问题。

**研究内容与结论**
研究人员采用有效半导体麦克斯韦-布洛赫方程（ESMBEs）模型，对单向环腔QCL的多模动力学进行数值模拟，并结合线性稳定性分析（LSA）和数值积分，提出了HFC形成的物理解释。研究发现，HFC源于一种有效拉比频率（ERF）与激光腔模式之间的共振现象。ERF依赖于场强度的平方根，类似于两能级系统中的拉比频率，但受非零α因子的调制。通过LSA，研究人员指出，当ERF接近某一腔模时，连续波（CW）发射失稳，触发自启动HFC，其阶数等于ERF与失稳CW频率之间的腔模数（或FSR倍数）。该结论在α<1和α>1两种情况下均得到验证，并揭示了α因子对失稳类型（RNGH型或Benjamin-Feir型）的影响。此外，通过改变泵浦参数和腔长，研究人员成功预测并数值模拟了不同阶数的HFC，且与实验观察（如THz环腔QCL中的HFC）一致。

**重要意义**
论文发表在《Light-Science》，其核心贡献在于：首次为QCL中自启动HFC的形成提供了清晰的物理图像，即ERF与腔模的共振机制。这一发现不仅解答了领域内长期悬而未决的问题，还为HFC的生成与控制提供了新思路——通过调节偏置条件或腔设计（如机械或电光调谐）来定制HFC特性，无需依赖外部注入或缺陷。该工作所发展的ERF估计和LSA方法可作为QCL结构设计的工具，加速HFC在片上集成光子学中的应用，包括中红外/THz脉冲生成、微波-太赫兹光子学及下一代无线通信中的载波生成。

**主要关键技术方法**
研究人员采用以下关键方法：1. 有效半导体麦克斯韦-布洛赫方程（ESMBEs）模型，用于描述环腔QCL的多模动力学，该模型包含非零α因子、非对称增益和折射率轮廓等半导体特性；2. 线性稳定性分析（LSA），对CW解进行小扰动，通过求解雅可比矩阵特征值获得参数增益，并识别ERF；3. 时域行波法数值积分，结合有限差分格式在时间和空间域求解ESMBEs。样本队列未涉及，但模拟参数基于典型THz QCL的实验值（如τe=7 ps，τd=0.1 ps）。

**研究结果**
**有效拉比频率（ERF）与阻尼**
研究人员通过分析ESMBEs中材料变量（极化P和载流子密度D）的短时扰动态，推导出ERF，其表达式为复数特征值的虚部。ERF随场强度X=|F|²的变化曲线显示：当α≠0时，低X区呈下降趋势，随后转换为√X增长；当α=0时，完全遵循两能级系统的√X趋势。同时，阻尼系数（实部绝对值）随X增大而显著降低，尤其在α=1.05（典型QCL值）时，高X区阻尼近乎消失，这为HFC的实验观测提供了依据。此外，快载流子动力学（大b值）导致阻尼减小，而量子阱激光器（QWL）因载流子寿命长（b小）阻尼高，无法支持ERF振荡和HFC。

**线性稳定性分析（LSA）与谐波梳形成**
对CW₀解（k=0）进行LSA，计算参数增益（特征值实部最大值）。对于典型THz QCL参数（μ=8.01），参数增益峰位于距CW频率8倍FSR的腔模处，且与ERF（红色标记）高度吻合。数值模拟证实该参数下系统稳定运行在8阶HFC。当泵浦参数μ增至11.03时，参数增益峰移至9倍FSR，ERF同步偏移，数值模拟得到9阶HFC。这表明HFC形成源于ERF与腔模的共振。对于α<1，失稳属RNGH型（参数增益气球下限非零）；对于α>1，失稳属Benjamin-Feir型（下限为零，即相位失稳），且更接近QCL实际（尤其在近阈值时）。

**腔长对谐波梳的影响**
固定α=0.92，改变腔长L：L=1 mm时，参数增益峰和ERF位于距CW频率2倍FSR处，模拟得2阶HFC；L=1.5 mm时，得3阶HFC；L=2 mm时，得4阶HFC。当L=1.75 mm（介于1.5和2 mm之间），参数增益峰值位于两相邻腔模之间，无正增益模，模拟得到稳定CW而非HFC。这一反例进一步证实了共振机制。对于α>1，类似现象被观察到，但参数增益线型不同。

**讨论与结论**
讨论部分指出，LSA在预测HFC形成中具有显著能力，且将自启动HFC解释为ERF与腔模共振的结果。与两能级系统不同，QCL中非零LEF使ERF阻尼更低，从而在超快激光器中触发CW失稳。此外，简化的复金兹堡-朗道方程（CGLE）仅能描述α>1时的长波（BF型）失稳，而ESMBEs模型揭示了更全面的失稳图景。研究结论翻译如下：
总之，研究人员为QCL中自启动HFC的形成提供了物理解释，这被认为是QCL频率梳领域的主要未解问题之一。理解HFC自发形成的物理机制，为生成和控制这些梳态带来了范式转变。以往研究主要依赖外部干预（如光注入、电注入或反射率不连续/缺陷），而本工作提出了替代方案，包括通过偏置条件或腔设计（如机械或电光调谐）来调节谐振器模式间距。此外，所发展的ERF估计和LSA方法可作为设计QCL结构的补充工具，以专门生成HFC态。由于QCL适合片上集成，这些发现有望加速基于HFC的光子学应用，包括中红外和THz脉冲生成、用于射频波形合成的微波和THz光子学，以及为下一代无线网络生成THz载波——这些领域因实验平台复杂性而长期处于起步阶段。