对称性、BIC 与降阶理论
阈值、增益、损耗与约束因子
读完本章,你应能#
严格区分阈值增益、阈值电流、被动 Q Q Q 和净模增益这几个常被混淆的量;
从冷腔损耗和有源区重叠出发推导 PCSEL 的光学阈值条件;
理解为什么高 Q Q Q 不自动等于最低阈值,也不自动等于最低阈值电流;
明确本章结论如何接入后续的 BIC、器件电学和动态学模块。
读前准备#
第 6 章 、第 8 章 中关于面内耦合、辐射通道、偏振与远场的基本概念。
阅读主线#
本章所处模型层级:L2
本章先问“阈值”这个词在 PCSEL 里到底指什么,再把材料增益、模增益、损耗和约束因子串成一条严格的符号链。随后我们说明被动 Q Q Q 与阈值之间的关系、为什么本章只能给出阈值增益而不能直接给出阈值电流,以及为什么阈值之上还必须继续进入动态学、注入和热反馈。
与第 6 章的连接#
第 6 章分析了 PCSEL 的面内反馈机制和垂直辐射通道。现在我们有了完整的损耗图像:α t o t a l = α r a d + α l o s s \alpha_{\rm total} = \alpha_{\rm rad} + \alpha_{\rm loss} α total = α rad + α loss 。本章的自然任务就是回答:需要多大的材料增益才能补偿这些损耗并开始激射?这就引出了阈值条件的严格推导。
阈值首先是净衰减被补偿的条件#
若只看某一候选模式的线性小信号演化,它的复振幅可抽象写成
d a d t = ( i ω 0 − γ r a d − γ i n t + γ o p t ) a , \frac{\dd a}{\dd t}=\left(\ii\omega_0-\gamma_{\mathrm{rad}}-\gamma_{\mathrm{int}}+\gamma_{\mathrm{opt}}\right)a, d t d a = ( i ω 0 − γ rad − γ int + γ opt ) a ,
其中 γ r a d \gamma_{\mathrm{rad}} γ rad 表示辐射损耗,γ i n t \gamma_{\mathrm{int}} γ int 表示内部损耗,γ o p t \gamma_{\mathrm{opt}} γ opt 表示有源区提供的光学增益率。于是模式的净衰减率可写为
γ n e t = γ r a d + γ i n t − γ o p t . \gamma_{\mathrm{net}} = \gamma_{\mathrm{rad}} + \gamma_{\mathrm{int}} - \gamma_{\mathrm{opt}}. γ net = γ rad + γ int − γ opt . (7.1)
阈值的最本质定义就是:净衰减率恰好为零。
γ o p t , t h = γ r a d + γ i n t . \gamma_{\mathrm{opt,th}} = \gamma_{\mathrm{rad}} + \gamma_{\mathrm{int}}. γ opt , th = γ rad + γ int . (7.2)
这是一个平衡条件,而非一个定义式。若把复振幅解成
a ( t ) = a ( 0 ) exp [ ( i ω 0 − γ n e t ) t ] , a(t)=a(0)\exp\!\left[\left(\ii\omega_0-\gamma_{\mathrm{net}}\right)t\right], a ( t ) = a ( 0 ) exp [ ( i ω 0 − γ net ) t ] ,
则模强度或光子数满足
S ( t ) = ∣ a ( t ) ∣ 2 = S ( 0 ) exp [ − 2 γ n e t t ] . S(t)=|a(t)|^2
=
S(0)\exp\!\left[-2\gamma_{\mathrm{net}}t\right]. S ( t ) = ∣ a ( t ) ∣ 2 = S ( 0 ) exp [ − 2 γ net t ] . (7.3)
于是只有当 γ n e t > 0 \gamma_{\mathrm{net}}>0 γ net > 0 时,模式才会衰减;当 γ n e t < 0 \gamma_{\mathrm{net}}<0 γ net < 0 时,线性模型会预言指数增长;而阈值点恰好是这两类行为的分界面。换句话说,式 (7.2) 的地位是“把增长和衰减的边界精确写出来”,而非“把一个符号定义成零”。
从材料增益到模增益:约束因子为什么必须出现#
材料增益并非直接作用在模式上的。模式只在有源区的一部分体积里有效“看到”增益,因此必须通过光学约束因子(optical confinement factor)把材料增益投影成模增益。最常用的近似写法是
g m o d = Γ o p t g m a t , g_{\mathrm{mod}} = \Gamma_{\mathrm{opt}}\, g_{\mathrm{mat}}, g mod = Γ opt g mat , (7.4)
其中 g m a t g_{\mathrm{mat}} g mat 是材料增益,Γ o p t \Gamma_{\mathrm{opt}} Γ opt 是模式与有源区的重叠因子,g m o d g_{\mathrm{mod}} g mod 是模增益。在能量归一化约定下,Γ o p t \Gamma_{\mathrm{opt}} Γ opt 的标准积分定义为
Γ o p t = ∫ V a c t ε r ( r ) ∣ E ( r ) ∣ 2 d V ∫ V t o t a l ε r ( r ) ∣ E ( r ) ∣ 2 d V , \Gamma_{\mathrm{opt}}=
\frac{\int_{V_{\mathrm{act}}}\varepsilon_r(\rvec)|\E(\rvec)|^2\,\mathrm{d}V}
{\int_{V_{\mathrm{total}}}\varepsilon_r(\rvec)|\E(\rvec)|^2\,\mathrm{d}V}, Γ opt = ∫ V total ε r ( r ) ∣ E ( r ) ∣ 2 d V ∫ V act ε r ( r ) ∣ E ( r ) ∣ 2 d V , (7.5)
其中 V a c t V_{\mathrm{act}} V act 为有源区体积,V t o t a l V_{\mathrm{total}} V total 为模式能量主要分布的全结构体积。对单量子阱(厚度约 8 nm)嵌入数百纳米波导层的典型 PCSEL 外延结构,Γ o p t ∼ 1 % \Gamma_{\mathrm{opt}}\sim1\% Γ opt ∼ 1% ;三至五量子阱时约为 3--5
这里的写法刻意保持在“最小接口层”。但它不是无条件成立的定义,而依赖几条近似:材料增益在有源区内可先看作局域标量;目标模式剖面在小信号增益下不发生显著重构;极化选择和增益色散已被压缩成该模式实际看到的一条有效增益曲线。若要把 g m a t g_{\mathrm{mat}} g mat 写成透明载流子密度、微分增益、量子阱数目与温度失谐的函数,要等到 第 18 章 才系统展开;本章先把它当作一个尚未完全展开、但已经明确属于有源区模型的输入量。
若进一步用能量速度 v e n v_{\mathrm{en}} v en 把空间增益换成时间增益率,则
γ o p t = v e n g m o d = v e n Γ o p t g m a t . \gamma_{\mathrm{opt}} = v_{\mathrm{en}}\, g_{\mathrm{mod}} = v_{\mathrm{en}}\Gamma_{\mathrm{opt}}\, g_{\mathrm{mat}}. γ opt = v en g mod = v en Γ opt g mat . (7.6)
同理,若把辐射损耗和内部损耗写成沿传播方向的等效损耗系数,则
γ r a d = v e n α r a d , γ i n t = v e n α i n t . \gamma_{\mathrm{rad}}=v_{\mathrm{en}}\,\alpha_{\mathrm{rad}},
\qquad
\gamma_{\mathrm{int}}=v_{\mathrm{en}}\,\alpha_{\mathrm{int}}. γ rad = v en α rad , γ int = v en α int . (7.7)
把 式 (7.7) 和 式 (7.6) 代回 式 (7.2) ,才得到下面常用的阈值增益公式:
g t h = α r a d + α i n t Γ o p t . g_{\mathrm{th}} = \frac{\alpha_{\mathrm{rad}}+\alpha_{\mathrm{int}}}{\Gamma_{\mathrm{opt}}}. g th = Γ opt α rad + α int . (7.8)
DimensionCheckNote[阈值增益公式的量纲自洽性] Γ o p t \Gamma_{\mathrm{opt}} Γ opt : 无量纲重叠因子; α r a d , α i n t \alpha_{\mathrm{rad}},\alpha_{\mathrm{int}} α rad , α int : 空间损耗系数 [m − 1 \mathrm{m}^{-1} m − 1 ]; g t h g_{\mathrm{th}} g th : 阈值材料增益系数 [m − 1 \mathrm{m}^{-1} m − 1 ]. 右边:[ α ] / [ Γ ] = m − 1 / 1 = m − 1 [\alpha]/[\Gamma] = \mathrm{m}^{-1}/1 = \mathrm{m}^{-1} [ α ] / [ Γ ] = m − 1 /1 = m − 1 ,与左边[ g t h ] = m − 1 [g_{\mathrm{th}}]=\mathrm{m}^{-1} [ g th ] = m − 1 一致. ✓ DimensionCheckNote
严格表述
本书在阈值链中统一使用 v e n v_{\mathrm{en}} v en (能量速度)表示将空间增益/损耗映射为时间增长率的参数。 对弱导模 slab,v e n v_{\mathrm{en}} v en 常可用 c / n e f f c/n_{\mathrm{eff}} c / n eff 作数量级近似; 对 Γ \Gamma Γ 点 standing-wave 模式,带图 Bloch 群速度为零,但 v e n v_{\mathrm{en}} v en 仍然存在且非零, 因为两者描述的是不同层级的物理:前者描述 Bloch 包络的传播速度,后者描述局域光场与增益介质的相互作用速率。 实际建模时,对 GaAs 基 PCSEL 可取 v e n ≈ 0.9 × 10 10 c m / s v_{\mathrm{en}}\approx0.9\times10^{10}\,\mathrm{cm/s} v en ≈ 0.9 × 1 0 10 cm/s 作为参考值, 但最终应与速率方程拟合或传递矩阵法对齐。
更严格的写法可从能量速度而不是带图群速度出发。设模式总储能为 U U U ,有源区提供的净增益功率为 P g a i n P_{\mathrm{gain}} P gain ,则定义时间增长率的直接方式是
γ o p t = P g a i n U . \gamma_{\mathrm{opt}}=\frac{P_{\mathrm{gain}}}{U}. γ opt = U P gain . (7.9)
若再把同一个模式沿传播方向的空间增益写成 g m o d g_{\mathrm{mod}} g mod ,并把能流定义为 P f l o w P_{\mathrm{flow}} P flow ,则
P g a i n = g m o d P f l o w , v e n ≡ P f l o w U , P_{\mathrm{gain}} = g_{\mathrm{mod}} P_{\mathrm{flow}},
\qquad
v_{\mathrm{en}}\equiv \frac{P_{\mathrm{flow}}}{U}, P gain = g mod P flow , v en ≡ U P flow , (7.10)
从而得到
γ o p t = v e n g m o d . \gamma_{\mathrm{opt}} = v_{\mathrm{en}}\, g_{\mathrm{mod}}. γ opt = v en g mod . (7.11)
这里出现的其实是能量速度 v e n v_{\mathrm{en}} v en 。对弱调制 slab 波导,它常可用 c / n e f f c/n_{\mathrm{eff}} c / n eff 作数量级近似;因此本章写作里的 v g v_g v g ,更准确地说应理解为 “把空间增益/损耗换算成时间增长/衰减率的等效能量速度参数”,而不是 exact Γ \Gamma Γ 点 band diagram 的严格 Bloch 群速度。
把 式 (7.2) 、式 (7.6) 联立,就得到常见的阈值增益形式
g t h = α r a d + α i n t Γ o p t , g_{\mathrm{th}} = \frac{\alpha_{\mathrm{rad}}+\alpha_{\mathrm{int}}}{\Gamma_{\mathrm{opt}}}, g th = Γ opt α rad + α int , (7.12)
其中 α r a d \alpha_{\mathrm{rad}} α rad 和 α i n t \alpha_{\mathrm{int}} α int 分别表示与辐射损耗和内部损耗对应的等效损耗系数。
这条式子需要按层级理解:它给出的首先是光学线性阈值 ,不是器件级阈值电流。
例子
[本书算例] 从损耗算到材料阈值增益。 若 α r a d = 5 / c m \alpha_{\mathrm{rad}}=\SI{5}{\per\centi\meter} α rad = 5 /cm 、α i n t = 10 / c m \alpha_{\mathrm{int}}=\SI{10}{\per\centi\meter} α int = 10 /cm ,量子阱光学重叠 Γ o p t = 0.02 \Gamma_{\mathrm{opt}}=0.02 Γ opt = 0.02 ,则阈值处所需模增益是 15 / c m \SI{15}{\per\centi\meter} 15 /cm ,对应材料增益
g t h = 5 + 10 0.02 = 750 / c m . g_{\mathrm{th}}=\frac{5+10}{0.02}=\SI{750}{\per\centi\meter}. g th = 0.02 5 + 10 = 750 /cm .
同一组损耗在较小重叠下会被放大成很高的材料增益要求;这就是为什么只报告冷腔 Q Q Q 仍不足以判断器件阈值。
在本书的语境里,g t h g_{\mathrm{th}} g th 表示“要让某模式净衰减归零所需的模增益阈值”。它依赖该模式的辐射损耗、内部损耗和与有源区的重叠。若把这个量直接改写成 I t h I_{\mathrm{th}} I th ,等于默认注入均匀、材料增益模型已知且热漂移可忽略,而这些默认对真实 PCSEL 往往都不成立。
内部损耗的组成与标定#
很多初学者把内部损耗写成一个符号后就不再追问。但对 PCSEL 而言,这往往正是器件结论出错的地方。更稳妥的写法是把内部损耗继续拆开:
α i n t = α F C A + α a b s + α m e t a l + α s c a t t . \alpha_{\mathrm{int}} = \alpha_{\mathrm{FCA}} + \alpha_{\mathrm{abs}} + \alpha_{\mathrm{metal}} + \alpha_{\mathrm{scatt}}. α int = α FCA + α abs + α metal + α scatt . (7.13)
这里 α F C A \alpha_{\mathrm{FCA}} α FCA 表示自由载流子吸收,α a b s \alpha_{\mathrm{abs}} α abs 表示材料吸收或掺杂层吸收,α m e t a l \alpha_{\mathrm{metal}} α metal 表示金属邻近引入的损耗,α s c a t t \alpha_{\mathrm{scatt}} α scatt 表示与粗糙度、侧壁和工艺缺陷相关的散射损耗。图 7.1 以饼图形式直观展示了典型 PCSEL 器件中各项损耗的占比情况。为避免把经验值误读为“标准常数”,本书统一采用“范围 + 条件 + 来源”写法: [文献值] 对近红外 GaAs 基 PCSEL/半导体激光器的代表性外延与掺杂条件,常见有 α F C A ∼ 1 \alpha_{\mathrm{FCA}}\sim1 α FCA ∼ 1 --10 c m − 1 10\ \mathrm{cm}^{-1} 10 cm − 1 (随掺杂与载流子密度显著变化)、 α a b s ∼ 1 \alpha_{\mathrm{abs}}\sim1 α abs ∼ 1 --5 c m − 1 5\ \mathrm{cm}^{-1} 5 cm − 1 、 α m e t a l ∼ 10 \alpha_{\mathrm{metal}}\sim10 α metal ∼ 10 --100 c m − 1 100\ \mathrm{cm}^{-1} 100 cm − 1 (当模式尾场靠近金属层时)、 α s c a t t ≲ 5 c m − 1 \alpha_{\mathrm{scatt}}\lesssim5\ \mathrm{cm}^{-1} α scatt ≲ 5 cm − 1 [Coldren et al., 2012] [Zhou & Pan, 2023] ; [本书算例] 辐射损耗 α r a d \alpha_{\mathrm{rad}} α rad 在高-Q Q Q 候选模可低至 10 − 1 c m − 1 10^{-1}\ \mathrm{cm}^{-1} 1 0 − 1 cm − 1 量级,在亮模可达 10 c m − 1 10\ \mathrm{cm}^{-1} 10 cm − 1 以上,具体数值必须与结构、边界与提取方法一并报告。
PCSEL 中的损耗机制分解示意图。辐射损耗通常占主导地位(尤其是对于表面发射器件),其次是自由载流子吸收和界面散射损耗。降低总损耗是实现低阈值激射的关键途径之一。注意电流泄漏虽然不直接贡献光学损耗,但会降低注入效率η i n j \eta_{\mathrm{inj}} η inj ,从而间接提高阈值电流密度。图 7.1 · 原始 LaTeX/TikZ 图源独立渲染
PCSEL 之所以不能只看冷腔模场分布,原因就在这里:只要模式在纵向上被拉向重掺杂层、金属邻近区或粗糙界面,原本“在冷腔模型中表现良好”的冷腔候选模,阈值排序就可能被内部损耗完全改写。
为什么高 Q Q Q 有利,但不等于最低阈值#
从被动角度看,品质因子越高,说明能量衰减越慢。若只考虑辐射损耗,这当然通常有利于降低阈值。但对 PCSEL 来说,这个判断绝不能单独使用。高 Q Q Q 模并不必然具备以下优点:
与量子阱有最好的重叠;
对掺杂层、金属层和吸收层有最小重叠;
在电流注入和温升回写之后仍保持最优排序;
在高功率工作时仍保持合适的提取效率和远场形状。
因此,正确的说法是:高 Q Q Q 往往是低阈值的有利因素,但阈值由净模增益条件决定,而不是由 Q Q Q 单独决定。
图 7.2 把这个误区拆成一条更稳妥的判断链。高 Q Q Q 首先只影响辐射损耗入口;它要经过内部损耗、约束因子、注入效率和热回写之后,才可能变成低阈值电流和稳定工作窗口。
从被动高 Q Q Q 到阈值和工作窗口的层级链。高 Q Q Q 有助于降低辐射损耗,但总阈值还要同时检查内部损耗、模场与有源区重叠、注入效率和热回写后的模式排序。图 7.2 · 原始 LaTeX/TikZ 图源独立渲染
从速率方程再看一遍阈值#
若采用单模近似,光子数 S S S 的最简速率方程可写为
d S d t = v g [ Γ o p t g m a t ( N , T ) − α r a d − α i n t ] S + β s p R s p V m o d e . \frac{\dd S}{\dd t} = v_g\bigl[\Gamma_{\mathrm{opt}} g_{\mathrm{mat}}(N,T)-\alpha_{\mathrm{rad}}-\alpha_{\mathrm{int}}\bigr]S + \beta_{\mathrm{sp}}R_{\mathrm{sp}} V_{\mathrm{mode}}. d t d S = v g [ Γ opt g mat ( N , T ) − α rad − α int ] S + β sp R sp V mode . (7.14)
其中 N N N 为等效载流子浓度,T T T 为温度,β s p R s p V m o d e \beta_{\mathrm{sp}}R_{\mathrm{sp}} V_{\mathrm{mode}} β sp R sp V mode 表示耦入该模式的自发辐射项,V m o d e V_{\mathrm{mode}} V mode 为模式体积。
notebox 速率方程的量纲一致性。 式式 (7.14) 中各项必须具有相同量纲:
S S S 是光子数(无量纲计数),故 d S / d t \dd S/\dd t d S / d t 的单位是 s − 1 \mathrm{s}^{-1} s − 1
v g [ ⋯ ] S v_g[\cdots]S v g [ ⋯ ] S 项:v g v_g v g (m/s) × \times × α \alpha α (m− 1 ^{-1} − 1 ) × \times × S S S → \rightarrow → s − 1 \mathrm{s}^{-1} s − 1 ✓ \checkmark ✓
β s p R s p V m o d e \beta_{\mathrm{sp}}R_{\mathrm{sp}} V_{\mathrm{mode}} β sp R sp V mode 项:β s p \beta_{\mathrm{sp}} β sp (无量纲) × \times × R s p R_{\mathrm{sp}} R sp (m− 3 ^{-3} − 3 s− 1 ^{-1} − 1 ) × \times × V m o d e V_{\mathrm{mode}} V mode (m3 ^3 3 ) → \rightarrow → s − 1 \mathrm{s}^{-1} s − 1 ✓ \checkmark ✓
若省略 V m o d e V_{\mathrm{mode}} V mode ,则自发辐射项的量纲将是 m− 3 ^{-3} − 3 s− 1 ^{-1} − 1 ,与光子数变化率的量纲不匹配。 notebox
在单模、小信号、给定模式剖面的近似下,阈值条件就是光子数不再指数衰减:
Γ o p t g m a t ( N t h , T t h ) = α r a d + α i n t . \Gamma_{\mathrm{opt}} g_{\mathrm{mat}}(N_{\mathrm{th}},T_{\mathrm{th}})=\alpha_{\mathrm{rad}}+\alpha_{\mathrm{int}}. Γ opt g mat ( N th , T th ) = α rad + α int . (7.15)
于是 式 (7.2) 与 式 (7.15) 只是同一件事的两种写法:前者强调时间衰减率,后者强调模增益与损耗的平衡。
若希望再从开放模语言看一遍,可把候选模的复频率写成
ω ~ = ω r − i ( γ r a d + γ i n t − γ o p t ) . \tilde{\omega}
=
\omega_r-\ii\bigl(\gamma_{\mathrm{rad}}+\gamma_{\mathrm{int}}-\gamma_{\mathrm{opt}}\bigr). ω ~ = ω r − i ( γ rad + γ int − γ opt ) . (7.16)
在本书采用的时间因子约定下,ℑ ω ~ < 0 \Im\,\tilde{\omega}<0 ℑ ω ~ < 0 对应振幅随时间衰减,ℑ ω ~ > 0 \Im\,\tilde{\omega}>0 ℑ ω ~ > 0 对应线性增长,而阈值则对应
ℑ ω ~ t h = 0. \Im\,\tilde{\omega}_{\mathrm{th}}=0. ℑ ω ~ th = 0. (7.17)
因此,从复频率角度看,阈值条件并非另一条新物理,而只是把“净衰减归零”重写成“复频率虚部穿过零”的开放系统表达。它和 式 (7.2) 、式 (7.15) 必须给出同一个阈值。
这也正是后文 CWT 阈值条件的代数来源。若第 j j j 个 CWT 超模在固定模式归一化下可写成
λ j = v e n , j [ Γ o p t , j g m a t − α r a d , j − α i n t , j ] , \lambda_j
=
v_{\mathrm{en},j}
\bigl[
\Gamma_{\mathrm{opt},j} g_{\mathrm{mat}}
-\alpha_{\mathrm{rad},j}
-\alpha_{\mathrm{int},j}
\bigr], λ j = v en , j [ Γ opt , j g mat − α rad , j − α int , j ] , (7.18)
那么“哪一支超模先达到阈值”就等价于寻找最先满足 Re λ j = 0 \operatorname{Re}\lambda_j=0 Re λ j = 0 的那一支。也就是说, 式 (7.12) 、式 (7.15) 与后文 CWT 的 `max Re( λ j ) = 0 (\lambda_j)=0 ( λ j ) = 0 ` 并非两套 相互竞争的判据,而是单模阈值公式与多模耦合阈值公式的同一条桥。
式 (7.2) 、式 (7.15) 、式 (7.17) 给出的都只是线性光学阈值 :它们回答“哪一个候选模最先补偿净损耗”。它们不是阈上稳态激光解,因为尚未引入饱和和有源自洽场;也不是阈值电流,因为尚未求解注入分布、复合损耗与热回写。
非均匀增益为什么会改写模式排序#
真实器件里,增益往往不是均匀分布的。电流扩展、电流拥挤、热斑和载流子泄漏都会让局域增益依赖空间位置。此时更合适的模增益表达式是
这里分母中的 V a l l V_\mathrm{all} V all 应按“模式实际占据的全部空间”理解,而不是只积到波导芯层或有源层附近一个小盒子为止。对开放 slab 共振模,它至少要覆盖主要场能量所在区域,包括向上下包层延伸的倏逝尾场;否则约束因子和模增益都会被系统性高估。实际截断时,建议与能量归一化所用积分域保持一致:取一个使分母积分值与无穷域积分值之差小于 1 % 1\% 1% 的有限截断面,并在收敛检查中验证该截断不影响最终约束因子结果。
式 (7.19) 的意义非常深刻。它说明:即便两个模式拥有相同的冷腔损耗,若它们在横向上的场分布不同,而注入分布又不均匀,那么它们得到的有效模增益也会不同。这也是后面电注入和热效应能够改写模式排序的根本原因。
为什么本章只能给出阈值增益,不能给出阈值电流#
到这里,读者可能会问:既然阈值条件都写得这么明确了,为什么还不能直接算 I t h I_{\mathrm{th}} I th ?
原因并不复杂。要从 g t h g_{\mathrm{th}} g th 走到 I t h I_{\mathrm{th}} I th ,至少还需要解决三件事:
材料增益 g m a t g_{\mathrm{mat}} g mat 如何依赖载流子浓度和温度;
外加电流如何通过漂移--扩散和复合过程产生该载流子分布;
温升如何反过来改写折射率、损耗和增益谱。
这三步分别对应后面的有源区、电学和热耦合章节。换句话说,本章完成的是从冷腔模式到光学阈值的桥接,而不是从模式直接跳到器件电流。
为什么本章故意停在静态阈值,而不在这里直接讨论动态学#
从教学顺序上看,这里有必要再把边界说清楚。阈值条件回答的是“哪一个模式最先有资格起振”,而不是“起振以后它会怎样演化”。只要器件进入阈值以上区域,载流子抽空、增益压缩、模式竞争和热回写就会一起出现;这已经超出了本章的静态线性阈值范围。
因此,本章故意把问题停在 g t h g_{\mathrm{th}} g th 与阈值电流来源的分界线上。这样做是为了确保读者先把静态层级分清,再进入后面的 第 23 章 去讨论弛豫振荡、调制响应和线宽来源,而非回避动态学。若在这里就把阈值、稳定性和线宽混写在一起,最容易造成的错误恰恰是重新把“能起振”误当成“能稳定单模工作”。
从冷腔候选模到主动激光模:本书采用的理论层级#
到这里还应再补一块总说明。否则读者虽然已经知道“被动候选模、线性阈值模、阈上工作模不是一回事”,却仍可能不知道它们在更高层理论里究竟怎样排位。对本书而言,最有用的分层是:
冷腔候选模 :这是前面能带、BIC、CWT、RCWA、FDTD/FEM 反复追踪的对象。它们是在无饱和、无反转动态、或只含给定损耗的开放 Maxwell 问题里得到的候选共振模或候选模族。
小信号阈值模 :这是在既定冷腔候选模基础上,引入未饱和增益之后得到的线性起振判据。此时我们问的是“哪一支模最先把净损耗补偿到零”,而不是“阈上稳态场到底长什么样”。
主动激光模 :这是更高层的有源自洽对象。它不再把光场视为在固定复介电常数中的被动或线性响应解,而是要求 Maxwell 方程、介质极化、反转分布和饱和非线性共同满足自洽条件。若进一步讨论阈上稳态或随机噪声,还会进入 Maxwell--Bloch 或 steady-state ab initio laser theory(SALT)一类形式化[Türeci et al., 2007] [Türeci et al., 2009] 。
这三层之间的关系,是“前两层是被设计流程反复证明有用的工作近似,第三层则给出更完整的顶层坐标系”,而非“前两层错误、第三层才正确”。本书把冷腔模与小信号阈值模作为主线,并非否认主动激光模理论的重要性,而是基于三点考虑。
第一,对 PCSEL 入门与器件设计最先起决定作用的,通常是模式候选、损耗排序、注入分布和热反馈的主导关系,而非完整阈上自洽解。 这些问题若没先理清,即使直接写出更高层形式化,也很难得到可操作的设计判断。
第二,主动激光模理论要求更完整的材料极化与反转模型。 这会把教学重点迅速推向 Maxwell--Bloch 或更复杂的有源场自洽求解,而偏离本书“从 Maxwell/Bloch/CWT/数值方法走到器件建模链”的主线。
第三,真实 PCSEL 的器件不确定性常常先于顶层主动模理论成为瓶颈。 外延层、注入非均匀、热回写和有限尺寸边界若尚未被清楚建模,那么把全部精力投入阈上自洽场理论,未必比先把冷腔候选模和小信号阈值链条做扎实更有效。
因此,本书采用的立场是:冷腔模理论给出候选模目录,小信号阈值理论给出起振排序,而主动激光模理论提供更高层的理论封顶。前两者构成本书的工作主线,第三者则作为读者应当知道其存在、知道其位置、但本书不把它展开成主教学骨架的上层框架。只要把这三层关系说清楚,就不会再把“冷腔共振频率”“线性阈值模”与“真实阈上激光模”混成同一个概念。
本章应输出的结果#
如果把本章用在研究工作流中,最有价值的输出不应只是一个单独的 g t h g_{\mathrm{th}} g th 数字,而应至少包括:
这些结果比"某模式阈值增益为某数值"更能指导后续设计。图 7.3 以柱状图形式直观展示了不同结构参数下阈值增益各成分的贡献对比,帮助设计者快速识别限制性能的关键因素。
PCSEL 阈值增益分解示意图,展示辐射、吸收、散射与泄漏各项损耗贡献及增益裕度。图 7.3 · 原始 LaTeX/TikZ 图源独立渲染
从注入电流到阈值条件:为什么需要一张 Sankey 图#
公式链条
I i n j → N ( r ) → g m a t → g m o d → g t h I_{\mathrm{inj}}\rightarrow N(\rvec)\rightarrow g_{\mathrm{mat}}\rightarrow g_{\mathrm{mod}}\rightarrow g_{\mathrm{th}} I inj → N ( r ) → g mat → g mod → g th
在文字里看起来很线性,但实际每一跳都会“损耗/折算”一部分可用增益。为把这种逐级折算直观化,图 7.4 给出一张教学型 Sankey 示意图。它不追求某个器件的精确百分比,而是强调:阈值电流上升往往是多个环节的小损失叠加而成,而非由某一项单独造成。
阈值电流密度链路的 Sankey 风格示意。从注入电流到阈值条件并非一步映射,而是经历注入效率、辐射复合效率、增益谱匹配和模场重叠的逐级折算。图中灰色支路表示在各环节被“损耗/失配”消耗的份额。图 7.4 · 原始 LaTeX/TikZ 图源独立渲染
这张图和本章公式可一一对上:注入侧折算对应 式 (7.19) 的空间非均匀项;复合侧折算对应“材料增益不是全都转化为受激增益”;模场侧折算对应 式 (7.5) 中的约束因子;最终只剩能对冲 式 (7.8) 右端总损耗的那一部分。对工程设计而言,这比单独报告一个 g t h g_{\mathrm{th}} g th 数字更有诊断价值,因为它直接告诉你下一步该优先优化“注入”“重叠”还是“损耗”。
回看本章#
本章把阈值问题严格放回了它应在的层级:阈值首先是净模增益补偿总损耗的光学条件。材料增益必须通过约束因子与模式重叠转化为模增益;总损耗则必须拆成辐射与内部各项来源。只有把这些关系讲清楚,后面讨论阈值电流、电注入和热效应才不会偷换概念。与此同时,辐射损耗这一项并非孤立参数:在 第 9 章 中,读者会看到同一类 Γ \Gamma Γ 点辐射通道如何被 BIC / quasi-BIC 语言重新组织;而在 第 23 章 中,读者会看到“谁先达到阈值”还必须继续延伸为“阈上谁更稳定、谁的调制更快、谁的线宽更窄”。这也是 PCSEL 中静态阈值、辐射工程与动态学不可拆开的原因。
练习题#
enumerate [基础] 从复频率定义出发,推导被动 Q Q Q 与等效损耗系数之间的关系。
[基础] 解释为什么同样的材料增益谱,对两个场分布不同的模式会给出不同的模增益。
[进阶] 给出一个 Γ o p t \Gamma_{\mathrm{opt}} Γ opt 增大但阈值电流不降反升的物理场景。
[综合] 将 式 (7.19) 与电注入非均匀联系起来,说明它为什么是后续器件建模的关键接口式。 enumerate
延伸阅读#
itemize [Coldren et al., 2012] [Siegman, 1986] [Henry, 1982] 中关于阈值条件、速率方程与线宽来源的经典章节和论文。
读完 第 18 章 、第 19 章 后回看本章,会更清楚“阈值增益”与“阈值电流”之间相隔哪些物理环节。
若你想继续追问“辐射损耗本身是怎样被结构对称性重排的”,可接着读 第 9 章 ;若你想追问“达到阈值以后会怎样演化”,则应直接进入 第 23 章 。若你还想知道这条链在更高层有源自洽理论中的位置,可把 [Türeci et al., 2007] [Türeci et al., 2009] 作为背景参考,而不必把它当成本书主线的前置条件。 itemize