对称性、BIC 与降阶理论
耦合模理论、微扰观点与模式分裂
读完本章，你应能#
区分时间耦合模理论（Temporal Coupled-Mode Theory, TCMT）、空间耦合模理论与 PCSEL 专用 CWT 的角色差异。
掌握微扰公式如何把“小几何变化”“小折射率漂移”“小损耗改动”转换为频率漂移和损耗变化。
理解模式分裂是偏振选择、模竞争和远场改变的直接来源，而非数学细节。
读前准备#
建议已掌握 第 10 章 、第 11 章 、第 8 章 中关于 CWT、开放系统和对称性破缺的基础内容。
阅读主线#
本章所处模型层级：L1→L2
第 10 章 、第 11 章 给出 PCSEL 专用的耦合波语言，适合描述晶格散射网络。本章引入更通用的耦合模语言，先说明 TCMT 与 CWT 的关系，再用微扰公式描述结构小改动对频率和损耗的影响，最后将这些公式用于模式分裂与对称性破缺。
为什么在 CWT 之外还要学耦合模理论#
很多读者看到这里会自然提问：既然前两章已经有了 CWT，为什么还要再学一个耦合模理论？答案在于两者关注的对象不同。
CWT 的自然基底是由倒格矢连接的内部平面波分量，用于解释光子晶体如何组织反馈和辐射。TCMT 则把结构压缩为一个或几个与端口耦合的共振模；其基底是已形成的共振超模及其输入输出通道。
因此，当问题关注孔形如何改变 x x x y y y
时间耦合模理论的基本语言#
先考虑最简单的单模情形。若一个共振模的复振幅为 a ( t ) a(t) a ( t )
a ˙ = ( i ω 0 − γ i − γ e ) a + k T s + , s − = C s + + d a . \begin{aligned}\dot a &= \left(\ii\omega_0-\gamma_\mathrm{i}-\gamma_\mathrm{e}\right)a + \bm{k}^{\mathsf T}\bm{s}_+,
\\
\bm{s}_- &= \mathbf{C}\bm{s}_+ + \bm{d}\,a.
\end{aligned} a ˙ s − = ( i ω 0 − γ i − γ e ) a + k T s + , = C s + + d a . (12.1)、(12.2) 这里，ω 0 \omega_0 ω 0 γ i \gamma_\mathrm{i} γ i γ e \gamma_\mathrm{e} γ e s + \bm{s}_+ s + s − \bm{s}_- s − C \mathbf{C} C k \bm{k} k d \bm{d} d
为了让这组式子真正能和 RCWA、FDTD/FEM 的输出对接，还必须先说明归一化。教材里最稳妥的约定是把 a ( t ) a(t) a ( t ) 能量振幅 ，于是
∣ a ∣ 2 ∝ U c a v , ∣ s ± ∣ 2 ∝ P p o r t , |a|^2 \propto U_{\mathrm{cav}},
\qquad
|s_\pm|^2 \propto P_{\mathrm{port}}, ∣ a ∣ 2 ∝ U cav , ∣ s ± ∣ 2 ∝ P port , (12.3) 其中 U c a v U_{\mathrm{cav}} U cav P p o r t P_{\mathrm{port}} P port a a a J \sqrt{\mathrm{J}} J s ± s_\pm s ± W \sqrt{\mathrm{W}} W k \bm{k} k d \bm{d} d s − 1 / 2 \mathrm{s}^{-1/2} s − 1/2 ∣ a ∣ 2 |a|^2 ∣ a ∣ 2
实际从数值结果提取 TCMT 参数时，最常用的接口是被动反射/透射谱的共振峰拟合。以单端口近似为例，共振附近的反射系数可写成
r ( ω ) ≈ C 11 + i d 1 k 1 ω − ω 0 + i ( γ i + γ e ) , r(\omega)\approx C_{11}+\frac{\mathrm{i}d_1k_1}{\omega-\omega_0+\mathrm{i}(\gamma_i+\gamma_e)}, r ( ω ) ≈ C 11 + ω − ω 0 + i ( γ i + γ e ) i d 1 k 1 , (12.4) 其中极点线宽由 γ i + γ e \gamma_i+\gamma_e γ i + γ e 2 ( γ i + γ e ) 2(\gamma_i+\gamma_e) 2 ( γ i + γ e ) C 11 C_{11} C 11 γ e \gamma_e γ e γ i \gamma_i γ i Q Q Q Q = ω 0 / 2 ( γ i + γ e ) Q=\omega_0/2(\gamma_i+\gamma_e) Q = ω 0 /2 ( γ i + γ e ) γ e \gamma_e γ e
这组式子的物理意义如下：第一式表示模幅按自身复频率演化并受外部输入驱动；第二式表示输出同时包含非共振直接散射背景和经由共振模再辐射的部分。
TCMT 的最小“端口--腔模--端口”图像。PCSEL 语境下，端口可以是上方自由空间、下方衬底辐射通道，或等效的外部耦合通道；TCMT 正是把内部共振模与这些外部端口之间的能量交换压缩成有限维输入输出模型。图 12.1 · 原始 LaTeX/TikZ 图源独立渲染 对 PCSEL 来说，若我们把“上方自由空间出射”“下方衬底辐射”甚至“侧向端口”都视为外部通道，那么 TCMT 很适合描述谱线、共振峰、临界耦合和端口干涉现象。它尤其适用于分析被动散射谱和端口响应，而不要求我们显式追踪内部每一个倒格波分量。
从单模到多模：模式竞争和分裂如何写进 TCMT#
PCSEL 更常见的是多模或近简并模情形。此时把单个振幅 a a a a \bm{a} a
a ˙ = ( i Ω − Γ ) a + K s + , s − = C s + + D a . \begin{aligned}\dot{\bm{a}}
&=
\left(\ii\bm{\Omega}-\bm{\Gamma}\right)\bm{a}
+
\mathbf{K}\bm{s}_+,
\\
\bm{s}_-
&=
\mathbf{C}\bm{s}_+ + \mathbf{D}\bm{a}.
\end{aligned} a ˙ s − = ( i Ω − Γ ) a + K s + , = C s + + D a . (12.5)、(12.6) 其中 Ω \bm{\Omega} Ω Γ \bm{\Gamma} Γ
因此，TCMT 与 CWT 可视为同一有限维投影思想在不同基底上的表述。实际工作中，可先用 CWT 识别重要内部耦合通道，再把候选超模参数压缩为 TCMT 参数，用于解释散射谱、反射峰、Fano 线形和端口干涉。
微扰观点为什么如此重要#
在研究和设计阶段，常见问题通常是小改动：孔半径轻微变化、刻蚀深度偏差、温度升高导致折射率上升、重掺杂层吸收增加、金属位置偏移等。若每次都依赖全波重算，代价较高且不利于分析机制。微扰理论的作用，是把这些小改动映射为本征频率和损耗的一阶变化。
对于封闭或近似厄米系统，最熟悉的一阶频移公式是
δ ω ω ≈ − 1 2 ∫ V δ ε ( r ) ∣ E ( r ) ∣ 2 d V ∫ V ε ( r ) ∣ E ( r ) ∣ 2 d V . \frac{\delta\omega}{\omega}
\approx
-
\frac{1}{2}
\frac{\int_V \delta\varepsilon(\rvec)\,|\E(\rvec)|^2\,\dd V}
{\int_V \varepsilon(\rvec)\,|\E(\rvec)|^2\,\dd V}. ω δ ω ≈ − 2 1 ∫ V ε ( r ) ∣ E ( r ) ∣ 2 d V ∫ V δ ε ( r ) ∣ E ( r ) ∣ 2 d V . (12.7) 这个公式的意义并不难：如果某个区域的介电常数增加，而该模在这个区域正好场强很强，那么模式更“愿意”把能量放进去，频率通常会向低频方向漂移。反之，若扰动发生在场很弱的区域，频移就会小得多。
式 (12.7) 的价值在于把设计动作改写为场重叠问题。这与前文 CWT 中耦合系数的结构一致：扰动是否有效，取决于它是否作用在目标模式具有显著场强或双正交重叠的区域。
式 (12.7) 严格依赖厄米归一化和小扰动条件。对强开放的 PCSEL 辐射模，真正严格的处理需要使用准正规模（quasinormal mode）或等价的开放系统归一化框架。工程上常把 式 (12.7) 当作第一近似使用，但必须知道它在强辐射、强非厄米情况下只是一种近似，不是普适精确公式。
色散材料的微扰修正
公式式 (12.7) 假设非色散材料 (∂ ε / ∂ ω ≈ 0 \partial\varepsilon/\partial\omega \approx 0 ∂ ε / ∂ ω ≈ 0 [Sauvan et al., 2013] [Alpeggiani et al., 2017] :
δ ω ω ≈ − 1 2 ∫ V [ δ ε + ω ∂ ( δ ε ) ∂ ω ] ∣ E ∣ 2 d V ∫ V [ ε + ω ∂ ε ∂ ω ] ∣ E ∣ 2 d V . \frac{\delta\omega}{\omega} \approx -\frac{1}{2} \frac{\int_V \left[\delta\varepsilon + \omega\frac{\partial(\delta\varepsilon)}{\partial\omega}\right]|\E|^2\dd V}{\int_V \left[\varepsilon + \omega\frac{\partial\varepsilon}{\partial\omega}\right]|\E|^2\dd V}. ω δ ω ≈ − 2 1 ∫ V [ ε + ω ∂ ω ∂ ε ] ∣ E ∣ 2 d V ∫ V [ δ ε + ω ∂ ω ∂ ( δ ε ) ] ∣ E ∣ 2 d V . (12.8) 在 PCSEL 常见的工作波长附近 (近红外),GaAs 和 InP 的材料色散 ∂ n / ∂ ω \partial n/\partial\omega ∂ n / ∂ ω > 10 n m >10\,\mathrm{nm} > 10 nm
先用一个 2 × 2 2\times2 2 × 2 #
在直接写开放 Maxwell 的一般公式之前，先看一个最小非厄米矩阵
M = [ a b c d ] , b ≠ c 或 a , d , b , c ∈ C . \mathbf{M}=
\begin{bmatrix}
a & b\\
c & d
\end{bmatrix},
\qquad b\neq c\ \text{或}\ a,d,b,c\in\mathbb{C}. M = [ a c b d ] , b = c 或 a , d , b , c ∈ C . (12.9) 设其某个本征值为 λ \lambda λ
M v R = λ v R , \mathbf M \bm v^{\mathrm R}=\lambda \bm v^{\mathrm R}, M v R = λ v R , 可取
v R = [ b λ − a ] . \bm v^{\mathrm R}=
\begin{bmatrix}
b\\
\lambda-a
\end{bmatrix}. v R = [ b λ − a ] . (12.10) 左本征矢则满足
( v L ) T M = λ ( v L ) T , (\bm v^{\mathrm L})^{\mathsf T}\mathbf M=\lambda (\bm v^{\mathrm L})^{\mathsf T}, ( v L ) T M = λ ( v L ) T , 等价于 M T v L = λ v L \mathbf M^{\mathsf T}\bm v^{\mathrm L}=\lambda \bm v^{\mathrm L} M T v L = λ v L
v L = [ c λ − a ] . \bm v^{\mathrm L}=
\begin{bmatrix}
c\\
\lambda-a
\end{bmatrix}. v L = [ c λ − a ] . (12.11) 对比 式 (12.10) 、式 (12.11) 可见：除非问题退回厄米或特定 complex-symmetric 情形，左矢量并不等于右矢量的复共轭。在非厄米问题中，关闭投影依赖左/右本征矢的双线性配对，而不是右本征矢自身的厄米内积。
这就是开放 Maxwell 微扰中引入伴随模的最小代数原因。左模可先理解为投影时的测试函数；在开放系统中，E L \E^{\mathrm L} E L E R \E^{\mathrm R} E R
数值演示：
2 × 2 2\times2 2 × 2 非厄米矩阵的左右模
取a = 1 , d = 2 , b = 0.1 , c = 0.3 a=1, d=2, b=0.1, c=0.3 a = 1 , d = 2 , b = 0.1 , c = 0.3 M = 1 0.1 0.3 2 \mathbf{M}=\begin{smallmatrix}1&0.1\\0.3&2\end{smallmatrix} M = 1 0.3 0.1 2
λ ± = 3 ± 1 + 0.12 2 ≈ 3 ± 1.058 2 . \lambda_\pm = \frac{3\pm\sqrt{1+0.12}}{2} \approx \frac{3\pm1.058}{2}. λ ± = 2 3 ± 1 + 0.12 ≈ 2 3 ± 1.058 . 对λ + ≈ 2.029 \lambda_+ \approx 2.029 λ + ≈ 2.029 v + R = [ 0.1 , λ + − 1 ] T ≈ [ 0.1 , 1.029 ] T \bm{v}_+^R=[0.1,\lambda_+-1]^T \approx [0.1,1.029]^T v + R = [ 0.1 , λ + − 1 ] T ≈ [ 0.1 , 1.029 ] T v + L T M = λ + v + L T \bm{v}_+^L{}^T\mathbf{M}=\lambda_+\bm{v}_+^L{}^T v + L T M = λ + v + L T v + L = [ 0.3 , λ + − 1 ] T ≈ [ 0.3 , 1.029 ] T \bm{v}_+^L=[0.3,\lambda_+-1]^T \approx [0.3,1.029]^T v + L = [ 0.3 , λ + − 1 ] T ≈ [ 0.3 , 1.029 ] T v − L ⋅ v + R = 0.3 × 0.1 + ( λ − − 1 ) ( λ + − 1 ) = 0 \bm{v}_-^L{}\cdot\bm{v}_+^R = 0.3\times0.1 + (\lambda_--1)(\lambda_+-1)=0 v − L ⋅ v + R = 0.3 × 0.1 + ( λ − − 1 ) ( λ + − 1 ) = 0 2 × 2 2\times2 2 × 2
用非厄米微扰公式估计孔形误差的影响
设设计好的圆孔 PCSEL 因工艺误差变成椭圆，介电常数扰动为δ ε ( r ) \delta\varepsilon(\bm{r}) δ ε ( r ) E n R ( r ) \bm{E}_n^R(\bm{r}) E n R ( r ) E n L ( r ) \bm{E}_n^L(\bm{r}) E n L ( r )
δ ω n ≈ − ω n 2 ∫ δ ε E n L ⋅ E n R d V ∫ ε E n L ⋅ E n R d V . \delta\omega_n \approx -\frac{\omega_n}{2}\frac{\int\delta\varepsilon\,\bm{E}_n^L\cdot\bm{E}_n^R\dd V}
{\int\varepsilon\,\bm{E}_n^L\cdot\bm{E}_n^R\dd V}. δ ω n ≈ − 2 ω n ∫ ε E n L ⋅ E n R d V ∫ δ ε E n L ⋅ E n R d V . 注意分子中出现的是左模与右模的内积，而不是通常的∣ E ∣ 2 |\bm{E}|^2 ∣ E ∣ 2
为什么非厄米一阶扰动不能直接照搬厄米公式#
上一节只给出了代数例子。开放系统微扰的关键点是：非厄米问题中，右模本身不足以关闭微扰公式。 原因不只是频率变为复数，而是开放 Maxwell 算子不再自伴，右模之间也不再按厄米内积正交。于是，厄米情形默认成立的三件事会失效：
归一化不再能简单用 ∫ ε ∣ E ∣ 2 d V = 1 \int \varepsilon |\E|^2 \dd V=1 ∫ ε ∣ E ∣ 2 d V = 1
投影时不能默认用右模自己的复共轭来消去一阶模形修正；
频移与损耗修正必须被统一看成复频率修正，而不是先算实部、再临时补虚部。
最稳妥的写法，是先把开放系统本征问题看成一般的非线性非厄米本征问题：
L ^ ( ω ~ , p ) E R = 0 , \hat{\mathcal L}(\tilde{\omega},p)\E^{\mathrm{R}} = 0, L ^ ( ω ~ , p ) E R = 0 , (12.12) 其中 p p p 伴随算子 与左模 ：
L ^ A ( ω ~ , p ) E L = 0 , \hat{\mathcal L}^{\mathrm{A}}(\tilde{\omega},p)\E^{\mathrm{L}} = 0, L ^ A ( ω ~ , p ) E L = 0 , (12.13) 这里 L ^ A \hat{\mathcal L}^{\mathrm{A}} L ^ A
⟨ u , L ^ v ⟩ = ⟨ L ^ A u , v ⟩ + 边界项。 \langle u,\hat{\mathcal L}v\rangle = \langle \hat{\mathcal L}^{\mathrm{A}}u, v\rangle + \text{边界项}。 ⟨ u , L ^ v ⟩ = ⟨ L ^ A u , v ⟩ + 边界项 。 (12.14) 直观上，把算子从右边搬到左边，在你选定的内积下整理分部积分项，就得到伴随算子 。在封闭厄米问题里边界项为零且 L ^ A = L ^ \hat{\mathcal L}^{\mathrm{A}}=\hat{\mathcal L} L ^ A = L ^
对一般非厄米问题，真正起投影作用的是左模与右模形成的双正交结构，而不是右模与自身的厄米内积。可把左模先理解为“做投影时的测试函数”，这会比“寻找另一支物理场”更直观。 现在对 式 (12.12) 做一阶微扰。若参数改变 p → p + δ p p\to p+\delta p p → p + δ p
[ ∂ L ^ ∂ ω ∣ n δ ω ~ + δ L ^ ] E n R + L ^ n δ E n R = 0. \left[
\left.\frac{\partial \hat{\mathcal L}}{\partial \omega}\right|_n \delta\tilde{\omega}
+ \delta \hat{\mathcal L}
\right]\E_n^{\mathrm{R}}
+
\hat{\mathcal L}_n\,\delta \E_n^{\mathrm{R}}
=
0. [ ∂ ω ∂ L ^ n δ ω ~ + δ L ^ ] E n R + L ^ n δ E n R = 0. (12.15) 再用左模对它做投影，由于 E n L \E_n^{\mathrm{L}} E n L
δ ω ~ n = − ⟨ E n L , δ L ^ E n R ⟩ ⟨ E n L , ∂ L ^ ∂ ω ∣ n E n R ⟩ . \delta\tilde{\omega}_n
=
-
\frac{\left\langle \E_n^{\mathrm{L}},\,\delta \hat{\mathcal L}\,\E_n^{\mathrm{R}}\right\rangle}
{\left\langle \E_n^{\mathrm{L}},\,\left.\dfrac{\partial \hat{\mathcal L}}{\partial \omega}\right|_n \E_n^{\mathrm{R}}\right\rangle}. δ ω ~ n = − ⟨ E n L , ∂ ω ∂ L ^ n E n R ⟩ ⟨ E n L , δ L ^ E n R ⟩ . (12.16) 这是厄米频移公式在开放系统中的对应形式。分子仍表示扰动与模式的重叠，分母则不再是普通储能积分，而是由伴随投影和频率导数共同定义的广义归一化量 。
若材料存在频率色散，这个分母最好至少写出最小形式。对
L ^ ( ω ) = ∇ × μ r − 1 ∇ × − ω 2 c 2 ε r ( ω ) \hat{\mathcal L}(\omega)=
\nabla\times \mu_r^{-1}\nabla\times
-\frac{\omega^2}{c^2}\varepsilon_r(\omega) L ^ ( ω ) = ∇ × μ r − 1 ∇ × − c 2 ω 2 ε r ( ω ) 这类 Maxwell 算子，有
⟨ E n L , ∂ L ^ ∂ ω E n R ⟩ ≈ − 1 c 2 ∫ ∂ ( ω 2 ε r ( ω ) ) ∂ ω E n L ⋅ E n R d V . \left\langle \E_n^{\mathrm{L}},\,\frac{\partial \hat{\mathcal L}}{\partial \omega}\E_n^{\mathrm{R}}\right\rangle
\approx
-\frac{1}{c^2}
\int
\frac{\partial\!\big(\omega^2\varepsilon_r(\omega)\big)}{\partial \omega}
\,
\E_n^{\mathrm{L}}\cdot \E_n^{\mathrm{R}}
\,
\dd V. ⟨ E n L , ∂ ω ∂ L ^ E n R ⟩ ≈ − c 2 1 ∫ ∂ ω ∂ ( ω 2 ε r ( ω ) ) E n L ⋅ E n R d V . (12.17) 若再退回无色散近似，∂ ( ω 2 ε r ) / ∂ ω → 2 ω ε r \partial(\omega^2\varepsilon_r)/\partial\omega \to 2\omega\varepsilon_r ∂ ( ω 2 ε r ) / ∂ ω → 2 ω ε r
式 (12.16) 比 式 (12.7) 更一般，因为它保留了开放系统中不能省略的两点：左/右模区分，以及算子对频率的显式依赖。只要材料色散、辐射边界或 PML 已进入模型，这两点通常都不能省略。
双正交、伴随模与 QNM 归一化在这里各自扮演什么角色#
到这里很容易产生另一个疑问：左模到底是不是右模的复共轭？答案是，一般并非。对真正的厄米问题，左模与右模在适当内积下确实可以重合；但对开放系统，它们只在更特殊的条件下才会简化。对互易、无磁光非互易的很多 Maxwell 问题，算子常常在某个无复共轭的双线性配对 下呈现 complex-symmetric 结构，此时左模更接近“转置伙伴”，而不是“共轭伙伴”。这也是为什么 QNM 归一化中常出现 E ⋅ E \E\cdot\E E ⋅ E ∣ E ∣ 2 |\E|^2 ∣ E ∣ 2
在实际数值中，常见的一种 QNM 广义归一化可写成
⟨ ⟨ Ψ m , Ψ n ⟩ ⟩ = ∫ Ω + P M L [ E m ⋅ ∂ ( ω ε ) ∂ ω E n − H m ⋅ ∂ ( ω μ ) ∂ ω H n ] d V , \left\langle\!\left\langle \Psi_m,\Psi_n\right\rangle\!\right\rangle
=
\int_{\Omega+\mathrm{PML}}
\left[
\E_m\cdot \frac{\partial(\omega\varepsilon)}{\partial\omega}\E_n
-
\Hf_m\cdot \frac{\partial(\omega\mu)}{\partial\omega}\Hf_n
\right]\dd V, ⟨ ⟨ Ψ m , Ψ n ⟩ ⟩ = ∫ Ω + PML [ E m ⋅ ∂ ω ∂ ( ω ε ) E n − H m ⋅ ∂ ω ∂ ( ω μ ) H n ] d V , (12.18) 其中 Ψ n = ( E n , H n ) \Psi_n=(\E_n,\Hf_n) Ψ n = ( E n , H n )
开放系统归一化通常不再使用普通能量范数；
材料色散会直接进入归一化分母，而不是只在“材料参数表”里出现；
一旦归一化换了，微扰公式的分母也必须随之更换，不能继续沿用封闭腔形式。
更系统的函数空间直觉、QNM 定义和归一化讨论见 第 3 章 、第 28 章 。本章此处的目标，是让读者在使用微扰语言时，不会再把开放系统误当成“只是多了一个复数频率的封闭腔”。
开放系统中微扰不只改频率，还会改损耗#
对 PCSEL 而言，只讨论频移是不够的，因为很多设计动作更直接改写的是损耗，而非共振频率本身。举例来说，轻微打破对称性可能几乎不改变两支模式的频率，却显著改变它们与辐射通道的耦合强度，于是损耗分裂比频率分裂更明显。
在开放系统里，更一般地可以把复频率写成
ω ~ = ω − i γ , \tilde\omega=\omega-\ii\gamma, ω ~ = ω − i γ , 因此，微扰对应复频率变化 δ ω ~ \delta\tilde\omega δ ω ~ 式 (12.16) 在开放系统中输出的是统一的频移与线宽变化。模式工程不仅改变频率，也改变模式与开放通道的耦合；这解释了 第 8 章 中偏振选择常主要体现为损耗分裂而非频率分裂。
从设计角度看，这一结论意味着：只检查能带或共振频率分裂是不够的，还必须检查辐射损耗是否重排；后者常直接决定阈值和单模性。
哪些工程公式可以继续用，哪些地方必须升级到开放系统严格形式#
到这里，读者最需要的是一条清楚的使用边界，而不是更多公式。对 PCSEL 设计，实践上可分成三层。
第一层是可以把厄米公式当作趋势代理 的情形：模式泄露较弱、候选模远离强辐射通道重组区、扰动主要发生在腔内高场区域而不是出光边界附近、并且目标模与竞争模已经明显分离。在这种情况下，式 (12.7) 往往能较可靠地判断“往哪个方向漂移”。
第二层是必须至少用左/右模公式检查一次 的情形：模式本身就是强开放模、设计变量直接作用于出光孔形或上下辐射不对称、材料色散不可忽略、或你已经开始关心虚部变化、偏振损耗分裂和 quasi-BIC 附近的辐射调控。此时，若仍只用 式 (12.7) ，最容易把“频移方向判断正确”误当成“损耗排序也正确”。
第三层是微扰理论本身不宜作为主模型 的情形：近简并模强混合、模式身份交换、扰动显著改变阵列边界或电流孔径、热反馈导致空间分布整体重构，或系统靠近异常点（exceptional point）等强非厄米重组区域。此时应回到完整全波或多物理场本征问题重算，而不是继续扩展一阶公式。
两模耦合为什么足以抓住“分裂”的核心#
尽管真实 PCSEL 可能涉及多模网络，但理解模式分裂时，一个两模模型通常已经足够说明问题。设两个近简并模式的复频率分别为 ω ~ 1 \tilde\omega_1 ω ~ 1 ω ~ 2 \tilde\omega_2 ω ~ 2 κ 12 \kappa_{12} κ 12 κ 21 \kappa_{21} κ 21
M = [ ω ~ 1 κ 12 κ 21 ω ~ 2 ] . \mathbf{M}
=
\begin{bmatrix}
\tilde\omega_1 & \kappa_{12}\\
\kappa_{21} & \tilde\omega_2
\end{bmatrix}. M = [ ω ~ 1 κ 21 κ 12 ω ~ 2 ] . (12.19) 其本征值为
ω ~ ± = ω ~ 1 + ω ~ 2 2 ± κ 12 κ 21 + ( ω ~ 1 − ω ~ 2 2 ) 2 . \tilde\omega_{\pm}
=
\frac{\tilde\omega_1+\tilde\omega_2}{2}
\pm
\sqrt{\kappa_{12}\kappa_{21}+\left(\frac{\tilde\omega_1-\tilde\omega_2}{2}\right)^2}. ω ~ ± = 2 ω ~ 1 + ω ~ 2 ± κ 12 κ 21 + ( 2 ω ~ 1 − ω ~ 2 ) 2 . (12.20) ReciprocityConditionBox[互易性定理的适用条件与非互易机制] 经典 Lorentz 互易定理表述为：对各向同性、无磁光效应的线性介质，若 ( E 1 , H 1 ) (\E_1,\Hf_1) ( E 1 , H 1 ) ( E 2 , H 2 ) (\E_2,\Hf_2) ( E 2 , H 2 )
∫ V ( E 1 ⋅ J 2 − E 2 ⋅ J 1 ) d V = ∮ S ( E 1 × H 2 − E 2 × H 1 ) ⋅ d S . \int_V (\E_1\cdot\Jf_2 - \E_2\cdot\Jf_1)\,\dd V = \oint_S (\E_1\times\Hf_2 - \E_2\times\Hf_1)\cdot\dd\bm{S}. ∫ V ( E 1 ⋅ J 2 − E 2 ⋅ J 1 ) d V = ∮ S ( E 1 × H 2 − E 2 × H 1 ) ⋅ d S . 在无耗散的封闭系统中，这进一步使耦合矩阵满足 κ 12 = κ 21 ∗ \kappa_{12} = \kappa_{21}^\ast κ 12 = κ 21 ∗
互易性成立的条件：
介质张量对称：ε ^ = ε ^ T \hat{\varepsilon} = \hat{\varepsilon}^{\mathsf T} ε ^ = ε ^ T μ ^ = μ ^ T \hat{\mu} = \hat{\mu}^{\mathsf T} μ ^ = μ ^ T
线性响应：无频率转换或非线性混合
封闭系统或适当的辐射边界条件
开放系统中的有效非厄米耦合： 即使介质本身满足互易性，开放系统的辐射通道也会诱导有效耗散耦合，使系统矩阵呈现非厄米结构。这是因为：
模式通过共同辐射通道发生间接耦合，耦合强度依赖于远场干涉条件
辐射损耗矩阵 Γ r a d \bm{\Gamma}_{\mathrm{rad}} Γ rad
在 biorthogonal 基底下，左模与右模不构成简单的复共轭关系
这并不等同于破坏 Lorentz 互易性：在互易条件下，散射矩阵仍满足 S = S T \mathbf S=\mathbf S^{\mathsf T} S = S T κ 12 = κ 21 \kappa_{12}=\kappa_{21} κ 12 = κ 21
式 (12.20) 概括了模式分裂的基本机制：当两模原本简并时，非零耦合会解除简并；若 ω ~ 1 \tilde\omega_1 ω ~ 1 ω ~ 2 \tilde\omega_2 ω ~ 2
对称性破缺怎样导致模式分裂#
现在把上述两模语言落回 第 8 章 的物理图像。设某个理想对称器件在 Gamma 点处支持一对简并或近简并模式，它们对应两种正交偏振或两种不同对称性表示。在完全对称时，系统对这两种模式没有偏好；一旦加入椭圆孔、双晶格扰动、上下一侧更强的辐射出口，或者热、电注入造成横向不均匀，对称性平衡就会被打破，于是出现两类可能的后果。
耦合谐振腔系统中微扰诱导的模式分裂现象。(a) 未微扰情况：两个完全相同的环形谐振腔通过倏逝波耦合，由于系统具有镜像对称性，两个局域模式形成简并对。(b) 施加微扰后：一侧谐振腔发生轻微形变或折射率变化，对称性破缺导致简并解除，能级分裂为ω ± = ω 0 ± Δ \omega_\pm=\omega_0\pm\Delta ω ± = ω 0 ± Δ 图 12.2 · 原始 LaTeX/TikZ 图源独立渲染 第一类后果是频率分裂：某一支模式的有效折射率或光程改变得更多，所以其共振频率偏移更大。
第二类后果是损耗分裂：某一支模式与辐射通道或损耗层的重叠增加得更多，所以其衰减更快或更慢。
在实际 PCSEL 设计中，这两类分裂常同时出现，只是主导程度不同。对于单模和偏振工程，关键在于它们如何改写净增益排序；分裂的重要性不在于谱线分开本身，而在于它决定哪一支模在阈值附近占优。
设方格晶格中一对原本简并的正交偏振 Gamma 点模，在引入轻微椭圆孔后出现分裂。若椭圆长轴沿 x x x x x x
例子
[本书算例] 把椭圆度灵敏度翻译成可测分裂。 在 λ 0 = 940 n m \lambda_0=\SI{940}{nm} λ 0 = 940 nm f 0 = c / λ 0 ≈ 319 T H z f_0=c/\lambda_0\approx\SI{319}{THz} f 0 = c / λ 0 ≈ 319 THz Δ f / f 0 = 5 × 10 − 5 \Delta f/f_0=5\times10^{-5} Δ f / f 0 = 5 × 1 0 − 5
Δ f ≈ 16 G H z , Δ λ ≈ λ 0 2 c Δ f ≈ 0.047 n m . \Delta f\approx\SI{16}{GHz},\qquad
\Delta\lambda\approx\frac{\lambda_0^2}{c}\Delta f\approx\SI{0.047}{nm}. Δ f ≈ 16 GHz , Δ λ ≈ c λ 0 2 Δ f ≈ 0.047 nm . 这里的 5 × 10 − 5 5\times10^{-5} 5 × 1 0 − 5
这个例子说明，微扰分析和两模分裂模型可直接连接对称性破缺与器件表现。
TCMT、CWT 和微扰三者怎样配合使用#
现在可以把这三套语言的分工说得更清楚。
图 12.3 给出了CWT与TCMT两种语言之间的映射关系：CWT聚焦于晶格内部的倒格矢散射网络，适合回答"哪种散射通道主导模式形成"；TCMT 聚焦于共振模与外部端口的能量交换，适合回答"模式如何与自由空间、衬底或波导耦合"。
CWT与TCMT的语言映射。(a)CWT 视角：关注晶格内部由倒格矢驱动的散射过程，基底是行波或驻波分量，耦合系数κ \kappa κ 图 12.3 · 原始 LaTeX/TikZ 图源独立渲染 图 12.3 的用途是区分建模对象：解释孔形如何影响偏振时，用 CWT 追踪散射通道；解释金属厚度如何影响阈值时，用 TCMT 追踪损耗通道；同时涉及两者时，再通过映射关系连接。
如果你的问题是"哪一种倒格矢散射通道在控制Gamma 点模式的重组"，优先使用 CWT。
如果你的问题是"某个共振模式和上方自由空间、下方衬底或旁路端口怎样干涉并形成谱线"，优先使用TCMT。
如果你的问题是"结构改动很小，我想先判断频率和损耗会朝哪个方向变化"，优先使用微扰公式。
三者之间通常是按问题层级配合使用，而非并列选择。常见路径是：先用 CWT 找到主导模和主导耦合；再在该模附近用微扰公式评估小改动；最后将候选参数区间送入 TCMT 或全波散射分析，检查端口响应和输出行为。这种串联方式比无约束参数扫描更便于解释机制。
把 tutorial 中的 Hermitian/非 Hermitian 控制翻译成三种语言
Noda 等人的高功率 PCSEL 路线中经常出现两个参数族：一个是由面内 180 ∘ 180^\circ 18 0 ∘ 90 ∘ 90^\circ 9 0 ∘ [Noda et al., 2023] [Yoshida et al., 2023] 。这套说法可以直接映射到本章三种语言。
在 CWT 语言中，Hermitian 耦合主要来自介电常数 Fourier 分量把四个基本 Bloch 波彼此连接；它决定 Γ \Gamma Γ
在 TCMT 语言中，前者更像共振模之间的内部耦合矩阵 Ω \bm{\Omega} Ω Γ r a d \bm{\Gamma}_{\mathrm{rad}} Γ rad
因此，本书采用的读法是：Hermitian/非 Hermitian 控制是 CWT、TCMT 与开放系统微扰在高功率 PCSEL 上的一次交汇，而非另一套独立理论。它的任务是把目标模与竞争模的复频率排序 做成可设计量，而不是只把某个被动 Q Q Q
一个常见错误是把TCMT 参数、CWT参数和全波仿真输出混在一起直接比较大小，却不检查它们的归一化约定是否一致。对开放系统而言，耦合系数、模幅和端口振幅的定义如果不统一，数值上的"差一个常数"会迅速演变成物理解释上的误判。
为什么本章内容会反复回到后面的器件章节#
从表面上看，微扰与模式分裂像是“纯光学细节”；但在真实 PCSEL 器件中，它们会持续被电注入和温升重新激活。温度上升引起的折射率漂移、本征增益峰漂移、自由载流子吸收上升，都会通过本章的语言转化为频移、损耗变化和模式分裂再排序。也正因为如此，后面的热-电-光耦合章节并非另起炉灶，而是把本章的线性小扰动语言变成空间分布和工作点依赖的器件语言。
因此，本章给出的是一种后续章节会反复使用的分析流程：对小改动，先判断其主要改变频率、损耗还是重叠，而非只适用于冷腔的数学技巧；对模式切换，先判断哪个简并被打破、哪个通道被重新加权。
回看本章#
TCMT、CWT 与微扰理论是互补的三种语言。CWT 用于追踪光子晶体内部耦合网络，TCMT 用于描述共振模与外部端口的交换，微扰理论用于把小结构或材料改动映射为频率和损耗的一阶变化。对开放 PCSEL 问题，严格的一阶扰动应写成左/右模投影形式，而不能简单照搬厄米频移公式；厄米公式在弱开放条件下仍可作为趋势代理，但在强辐射、quasi-BIC、近简并重组或显著色散区，需要使用开放系统归一化甚至完整重算。模式分裂的重要性在于它会改写偏振、远场和阈值排序。有效设计分析应根据问题在这三种语言之间切换。
练习题#
enumerate [基础] 用自己的话解释 TCMT 与 CWT 的区别，并说明为什么它们在 PCSEL 中都不可缺。
[基础] 从 式 (12.16) 的推导思路出发，解释为什么开放系统里必须引入左模或伴随模，而不能只用右模自己关闭投影。
[进阶] 结合 式 (12.7) 与 式 (12.16) ，说明为什么“把材料改在场强最大的地方”通常最有效，但并不一定总是最优设计策略。
[进阶] 从 式 (12.20) 出发，讨论在什么情况下损耗分裂会比频率分裂更直接地影响模式选择。
[综合] 设某 PCSEL 在低温下稳定单偏振工作、升温后出现偏振切换。请用“微扰 + 模式分裂”的语言给出一条可能的解释链。 enumerate
延伸阅读#
itemize 下一章开始进入数值方法部分。PWEM 负责给出候选模目录，本章的分裂与微扰观点则解释这些候选模在后续全波分析中的频率和损耗差异。
一般光学腔与激光模理论可参考 [Siegman, 1986] [Coldren et al., 2012] 中关于耦合模、模竞争和频率漂移的相关章节。
TCMT 与光子晶体 slab 共振的经典文献可参考 [Fan & Joannopoulos, 2002] [Fan et al., 2003] [Suh et al., 2004] 。
开放系统归一化与非厄米微扰可继续阅读 [Leung et al., 1994] [Muljarov et al., 2010] [Sauvan et al., 2013] [Kristensen et al., 2015] ，并与 第 28 章 对照。 itemize
