高功率掺镱双包层光纤激光器中横向模式不稳定性的研究：基于受激热瑞利散射的三层光学-热学分析
作者、机构与发表信息
本研究由巴西圣保罗大学（USP）微波与光子学功率实验室（LaMP）的Elbis Santos Cardoso和Cláudio Costa Motta，以及核能与能源研究所（IPEN-CNEN）激光与应用中心（CLA）的Ricardo Elgul Samad合作完成。该研究已于2025年发表在学术期刊 Micromachines 上。文章题目为《Transverse Mode Instability in High-Power Yb-Doped Double-Clad Fiber Lasers: A Three-Layer Optical–Thermal Analysis Based on Stimulated Thermal Rayleigh Scattering》，遵循知识共享署名许可（CC BY）协议，可供开放获取。
研究背景与目标
研究领域聚焦于高功率光纤激光器物理，具体关注限制其功率提升的核心瓶颈——横向模式不稳定性（Transverse Mode Instability， TMI）。掺镱双包层光纤激光器（Yb-Doped Double-Clad Fiber Lasers， YDCFLs）因其高效率和功率可扩展性，已成为千瓦级以上连续波高亮度激光系统的首选方案。然而，当功率提升到一定阈值时，会发生TMI现象，表现为激光能量从稳定的基模（LP01模）突然向高阶模（Higher-Order Modes， HOMs）转移，并伴有千赫兹量级的模态含量振荡，导致光束质量急剧恶化。
学术界普遍认为，TMI的物理根源是非线性效应——受激热瑞利散射（Stimulated Thermal Rayleigh Scattering， STRS）。其机制是：传输模式之间的干涉产生纵向调制的光强分布；通过量子缺陷加热（quantum-defect heating），在掺杂纤芯内诱导动态的折射率扰动；径向热扩散在光强调制和折射率响应之间引入了相位延迟，形成一个移动的热光栅，从而实现了模式间的相干耦合，并表现出典型的阈值行为。
尽管STRS理论已被广泛接受，但此前的大量研究在建模双包层光纤（Double-Clad Fiber， DCF）时，常采用等效的两区域近似，即将其简化为仅包含纤芯和单一包层的单包层光纤（Single-Clad Fiber， SCF）。这种简化模型虽然有助于阐明STRS和TMI的许多方面，但无法完全表征内包层的作用。内包层会改变模式的横向本征值结构，并在高功率下允许高阶模的场穿透到纤芯之外，这些效应会改变光学重叠积分，从而影响STRS驱动的耦合效率。
因此，本研究的目标是建立一个定量的、物理一致的分析模型，用于描述三层结构（纤芯、内包层、外包层）双包层光纤中，STRS对TMI的影响。该模型旨在克服单包层近似的局限，通过结合纤芯区域的弱导近似（Weakly Guiding Approximation， WGA）和包层界面处的半弱导近似（Semi-WGA）边界条件，精确计算高阶模在内包层的穿透以及关键的横向本征值（u01和umn）。在此框架内，评估包含增益饱和效应和多层几何结构热本征模的非线性STRS耦合系数，从而定量分析三层几何结构如何通过改变光场和热场的模态结构，最终影响预测的TMI阈值。这项研究旨在为分析掺镱双包层光纤系统中的热光动力学提供一个可靠的工具，并支持设计具有更高模态稳定性的下一代高功率光纤激光器。
研究流程与方法详述
本研究是一个理论建模与数值仿真相结合的工作，不涉及传统意义上的实验样本和湿实验流程，其“流程”核心是理论模型的构建、数值求解与结果分析。具体可分为以下几个主要步骤：
第一步：三层光纤中的光学模态理论构建 研究对象是一个特定几何参数的高功率掺镱双包层光纤激光振荡器。光纤参数为：纤芯半径 a = 12.5 µm， 内包层半径 b = 125 µm， 外包层半径 c = 200 µm；折射率 n1 > n2 > n3。泵浦为双向泵浦（λp = 976 nm），激光发射波长为 λs = 1064 nm， 光纤最佳长度优化为12米。 模型的核心是精确描述三层光纤中的导模。研究采用了Tsao等人的矢量公式，通过德拜势描述电磁传播，并在界面 r=a 和 r=b 处严格强制电场和磁场切向分量的边界条件。即使层间折射率对比度较高（双包层光纤的典型情况），该方法依然有效。 从麦克斯韦方程组和柱坐标变量分离出发，得到归一化参数：u1² = a²(k₀²n1² - β²)， u2² = a²(k₀²n1² - β²)， w3² = b²(β² - k₀²n3²)， 其中k₀为真空波数，β为传播常数。通过数值求解由电磁边界条件导出的特征（本征值）方程，得到β。由于纤芯与内包层间折射率接近（n1 ≈ n2），在r=a界面采用弱导近似（WGA）是有效的；而在r=b界面，n2与n3的对比度显著，WGA不再适用，因此采用Tan等人描述的半弱导近似（Semi-WGA）来连接各区域的场表达式（纤芯用贝塞尔函数Jm， 内包层用Jm和诺依曼函数Nm组合，外包层用修正贝塞尔函数Km）。通过强制界面连续性条件，确定场系数，并最终得到归一化的径向场分布fmn®，用于后续的光学和热学积分计算。
第二步：STRS耦合机制与热学模型建立 此步骤将光学模态与热动力学耦合起来。 1. 模态干涉与移动光栅：总光场表示为所有模式的叠加。LP01模与LPmn模的干涉产生了包含项 cos[(βmn - β01)z - (ωmn - ω01)t] 的移动辐照度（光强）图案，其空间周期由拍长 λ01,mn = 2π/|βmn - β01| 决定。这个移动的光强图案是STRS的根本驱动力。 2. 量子缺陷加热与谐波热源：模式干涉调制了掺杂纤芯内的局部热沉积。加热率由量子缺陷 QD = (λs/λp - 1) 决定。干涉产生的谐波热源项被形式化为 qt(r, φ， z， t) = q̃(r， z) e^(imφ) e^(i(qz-ωt))， 其中径向振幅q̃(r， z) 正比于量子缺陷、饱和有效增益 g̃(r, z) 以及两干涉模场f01®和fmn®的乘积。 3. 谐波热方程与热本征模：由热源qt产生的温度扰动Δt服从热方程。由于热源具有 e^(imφ)e^(i(qz-ωt)) 的依赖关系，求解形式也设定为此形式与径向热本征模 tml® 的线性组合。tml®是三层结构的径向热模，由静态热方程结合 r=b 处的连续性条件和 r=c 处的外边界条件求得。对应的热本征值 ηml 包含了内包层对径向热扩散的影响。通过将热方程投影到tml®上并利用径向正交性，得到温度扰动幅度al(ω, z)的表达式，其中引入了热阻尼参数 γml = (2κ0/ρc)(q² + η²ml)， 它明确了频率依赖性。 4. 热光扰动与非线性耦合系数：温度变化通过热光系数 kt = ∂n/∂T 引起折射率扰动 Δn = kt Δt。该扰动形成一个移动的热光栅，介导STRS。由此导出LP01与LPmn模之间的复非线性耦合系数 χmn(z, ω)。这个系数是模型的核心输出，它综合了量子缺陷、饱和增益、光学场重叠积分、热本征模重叠积分以及包含相位延迟信息的复频率响应函数 fml(ω) = 2(2ω/γml - i) / [1 + (2ω/γml)²]。其实部峰值出现在 ω ≈ γml/2， 虚部则捕捉了光栅与热响应之间的扩散诱导相位延迟，正是这种相位不对称性导致了STRS驱动的非互易耦合。
第三步：功率演化与阈值分析的数值模拟 研究使用MATLAB (R2016a)软件进行数值模拟，所有物理、几何和激光参数总结于论文的Table 1中。 1. 耦合功率方程：在准静态区域，模态功率的轴向演化由耦合方程描述：dP01/dz = g01(z)P01 - g01(z)Re{χmn}P01Pmn； dPmn/dz = [gmn(z) - αmn]Pmn + g01(z)Re{χmn}P01Pmn。其中g01(z)和gmn(z)是饱和模态增益，αmn是高阶模损耗。当Re{χmn}为正且足够大时，高阶模在光纤的大部分长度上获得净增益，导致其功率指数增长，即TMI的发生。 2. STRS微观阈值与TMI宏观阈值：STRI提供了触发模式间相干功率转移的微观机制。在高增益区域，导出一个简单的解析阈值公式：P01_th = (1/Re{χmn}) ln(x P01(0)/Pmn(0))， 其中x是定义阈值的高阶模输出功率分数（通常为1%-5%）。这个阈值定义了高阶模开始因热耦合而指数增长的微观尺度。TMI则是这种增长在放大器输出端表现出的宏观现象。 3. TMI阈值的解析估计：研究还采用了Zervas提出的渐近表达式来估计宏观TMI阈值：P_tmi^thr = [κ0 u²mn (u²mn - u²01)] / [4π neff (dn/dT) α‘s (λs/d0)²]。该公式将阈值与光纤的光学、热学和几何参数联系起来，其中α‘s是包含本底损耗和量子缺陷耗散项的有效损耗。该表达式虽然对复杂几何结构有定量限制，但可作为直接理解各参数如何影响不稳定性的设计参考。 4. 对比分析：为了凸显三层结构的影响，研究对具有相同纤芯半径和数值孔径的单包层光纤（SCF）进行了并行分析。SCF的光学模态采用经典的Gloge标量阶跃折射率色散关系计算，其热学模型也相应简化为两层。
主要研究结果
结果一：光学模态分布的差异 模拟结果显示，在SCF中，LP01和LP11模都紧密局限于纤芯内。而在DCF中，由于存在内包层，导模可以穿透到较低折射率的区域，将光功率重新分布到更大的横向区域，减少了局限于纤芯的功率比例。这种模态局域化的改变影响了有效折射率差和参与TMI阈值模型的物理参数。
结果二：模态本征值与TMI阈值预测 计算并比较了不同模型下的横向本征值及相应的阈值功率（见原文Table 2）。关键发现包括： * 使用Tsao等人的严格本征值时，DCF中LP01-LP11耦合通道的阈值最低（约0.66 kW），而LP01-LP13通道的阈值显著更高。 * 使用Tan的半弱导模型时，出现了u11 < u01的情况，这使得LP01-LP11通道在阈值表达式中失去物理意义，只有LP01-LP13通道可给出有意义的阈值（约0.99 kW）。 * 尽管SCF具有更强的光学限制，但其预测的阈值（约1.37 kW）反而高于DCF。在所使用的渐近框架下，这主要归因于两层和三层结构之间横向光学本征值u01和umn的显著差异。热本征值ηml的贡献是次要但一致的。
结果三：温度分布与STRS诱导的模态耦合演化 * 温度场：模拟显示了沿光纤轴向不对称的温度分布，在输入端（z=0）轴心温度最高（496.6°C），梯度最陡；在光纤中点附近（z≈6 m）轴心温度最低（129.9°C）。这种不对称性反映了泵浦吸收分布。 * 模态功率转移：通过积分耦合微分方程，得到了LP01和LP11模功率沿光纤的演化（原文图6，7）。在无模态耦合时，LP01单调增长。当引入STRS热光耦合后，约在z≈5.0 m处（对应P01≈420 W），LP01功率开始急剧下降，同时LP11功率快速增长并最终成为主导模，这标志着TMI的发生。该阈值与使用典型参数（Re{χ11} ≈ 1.3e-1 W⁻¹， 初始噪声比~10⁻²⁵）通过解析公式估算的值（约430 W）吻合良好，与Dong报道的数值阈值一致。 * 动态行为：超过阈值后，STRS动力学呈现非线性，能量可能间歇性地回流到基模，产生沿光纤的振荡行为，这与Dong描述的“再播种”机制一致。
结果四：横向模式不稳定性阈值的参数化分析 研究进一步将阈值表达为量子缺陷（通过改变λs）和纤芯直径d0的函数进行参数化分析。 * 对于固定的d0=25 µm， TMI阈值功率随量子缺陷增大（即λs增大）而单调降低。在λs=1064 nm， λp=976 nm条件下，计算得到DCF和SCF的阈值分别为0.66 kW和1.37 kW（原文图12）。 * 对于固定的工作波长（λs=1064 nm， λp=976 nm）， TMI阈值功率随纤芯直径d0的增大而升高。在d0=25 µm处得到的结果与基于量子缺陷参数化得到的结果一致（原文图13）。 * 核心结论：DCF表现出比SCF更低的TMI阈值。这一结果并非偶然，它源于一个统一的光-热机制：量子缺陷控制纵向热负载，而横向本征值决定该热分布如何通过STRS耦合LP01模与高阶模。DCF的三层几何结构改变了与模态光强图案相互作用的温度分布，这一效应被编码在横向本征值中，并导致模态因子u²mn(u²mn - u²01)更小，从而增强了负责STRS驱动TMI的热光相位调制效率。
研究结论与价值
本研究的结论是，所观察到的TMI是受激热瑞利散射（STRS）的一种表现：量子缺陷加热产生了一个延迟的折射率调制，形成一个负责相干模态耦合的移动热光栅。这一机制将Dong以及Smith & Smith提出的解释统一在一个一致的模态-热学框架下。
本研究构建并应用了一个耦合光-热公式来分析掺镱双包层光纤中的TMI，明确考虑了多层径向几何结构的影响。研究结果表明，用三层DCF替代两层SCF会改变横向模态本征值和相关的光-热重叠积分。对于所研究的几何结构，内包层施加的本征值结构导致了更强的STRS诱导耦合，因此尽管两种光纤的硅热参数基本相同，DCF的TMI阈值反而更低。本征值分析进一步证实，主导的不稳定性通道敏感地依赖于模态结构。
研究的亮点
其他有价值的发现
研究还通过二维图像直观展示了SCF和DCF中LP01与LP11模的横向光强分布，以及STRS驱动下两模相干叠加产生的动态干涉图样（原文图8-11）。这些可视化结果生动地揭示了TMI发生时光束质量退化的物理图景，即光功率在纤芯横截面内的周期性再分配，为理解不稳定性机制提供了直观辅助。