附录与索引

数学附录

读完本章，你应能

读前准备

建议已掌握 第 3 章、第 4 章、第 7 章、第 12 章 中关于 Maxwell 方程、Bloch 展开、阈值条件和微扰观点的内容。

阅读主线

本章所处模型层级：跨层数学支撑（L1→L4）

本附录的数学对象关系图：从函数空间与算符出发，连接到谱结构、模态展开与格林张量，最终落到可直接使用的公式。
本附录的数学对象关系图：从函数空间与算符出发，连接到谱结构、模态展开与格林张量，最终落到可直接使用的公式。图 28.1 · 原始 LaTeX/TikZ 图源独立渲染

图 图 28.1 先把本章数学对象之间的关系固定下来，再进入各节公式推导，可明显降低“符号会看、结构不清”的阅读负担。

首次阅读建议路径

若是第一次阅读本附录，建议按以下顺序进行：先读时间因子与傅里叶约定，再读周期/Bloch 展开与卷积结构，再读 QNM 与归一化，最后读互易与微扰公式。其余段落可在需要时回查。

本附录是把正文反复调用、却不宜每次都从头展开的数学骨架集中讲清楚，而非把若干常见恒等式堆在一起的公式仓库。我们先统一时间因子和频域约定，因为一切损耗、增益和复频率的符号都依赖它；随后建立周期介质中的傅里叶与 Bloch 展开，解释卷积结构为什么在 第 5 章 的能带概念、第 13 章 的 PWEM、以及后续 RCWA 和 CWT 中一再出现；再讨论内积、归一化和开放系统中的复频率与 QQ；最后给出互易性与一阶微扰公式，并明确说明它们在哪些条件下可靠、在哪些条件下只能作为工程近似。

为什么数学附录不是装饰，而是主结论的一部分

PCSEL 的建模链横跨电磁场、半导体输运和热传导。只要这些模块之间的时间因子、损耗符号、傅里叶约定或归一化不一致，最终就会出现一种很危险的情况：看起来像是“物理结论不同”，实际却只是数学约定变了。最典型的例子就是时间因子的选择。一旦一个章节使用 e−iωt\ee^{-\ii\omega t}，另一个局部推导却默默改成 e+iωt\ee^{+\ii\omega t}，那么介质虚部、损耗符号、复频率虚部和 QQ 的定义就会全部翻转。若读者没有意识到这种符号翻转，很容易把纯符号问题误读成腔损耗差异或增益符号差异。

因此，本附录的目标是把正文各章默认共享的数学地基显式写出来，而非补充“额外数学”。只要地基统一，后面的开放系统、耦合波理论、RCWA、FDTD/FEM 提取和器件级回写才可能真正接成一条链。

时间因子与频域约定：先把增益和损耗的符号说清楚

全书统一采用时间因子

E(r,t)=ℜ ⁣{E(r)e−iωt},H(r,t)=ℜ ⁣{H(r)e−iωt}.\E(\rvec,t) = \Re\!\left\{\E(\rvec)\ee^{-\ii\omega t}\right\}, \qquad \Hf(\rvec,t) = \Re\!\left\{\Hf(\rvec)\ee^{-\ii\omega t}\right\}.(28.1)

在这一约定下，频域 Maxwell 方程写成

∇×E=+iωB,∇×H=J−iωD.\begin{aligned}\nabla\times \E &= +\ii\omega\B,\\ \nabla\times \Hf &= \J - \ii\omega\D.\end{aligned}(28.2)、(28.3)

若采用线性各向同性本构关系

D=ε0εE,B=μ0μH,\D = \varepsilon_0\varepsilon\E, \qquad \B = \mu_0\mu\Hf,

则在无外加电流时可得到熟悉的波动方程。重要的是它对损耗和增益符号的约束，而非方程形式本身。

在约定 式 (28.1) 下，若介质为被动介质，则通常要求

ℑ ε(ω)>0,\Im\,\varepsilon(\omega) > 0,(28.5)

因为这对应正的平均耗散功率。这个结论可以从时均功率密度看出来。对无外加源的谐波场，时均耗散功率密度满足

Pabs=ωε02ℑ ε ∣E∣2,P_{\mathrm{abs}} = \frac{\omega\varepsilon_0}{2}\Im\,\varepsilon\,|\E|^2,(28.6)

其中假设磁损耗可以忽略或已单独处理。于是，在本书采用的符号体系里，正的虚部表示吸收，若要以等效介电常数虚部的方式表示增益，则会对应负的有效虚部。这条符号约定看似基础，却是所有“阈值增益符号到底该写成正还是负”的根源。

式 (28.5)、式 (28.6) 都依赖于本书统一采用的时间因子 \(^-ω t\)。若改用 \(^+ω t\)，则介质虚部和复频率虚部的物理判据会整体翻转。任何引用外部文献中的损耗或增益公式时，都必须先检查其时间因子约定。

为了避免正文不同章节在符号上“看起来一样、其实约定不同”，这里把本附录与主文完全一致的傅里叶变换约定写清楚。时间傅里叶变换采用

f^(ω)=∫−∞+∞f(t) e+iωt dt,f(t)=12π∫−∞+∞f^(ω) e−iωt dω,\hat f(\omega) = \int_{-\infty}^{+\infty} f(t)\,\ee^{+\ii\omega t}\,\dd t, \qquad f(t) = \frac{1}{2\pi}\int_{-\infty}^{+\infty}\hat f(\omega)\,\ee^{-\ii\omega t}\,\dd\omega,(28.7)

它与时间因子 式 (28.1) 保持同一符号体系。平面上的空间傅里叶变换采用

f~(kx,ky)=∬f(x,y) e−i(kxx+kyy) dx dy,\tilde f(k_x,k_y) = \iint f(x,y)\,\ee^{-\ii(k_x x+k_y y)}\,\dd x\,\dd y,(28.8)
f(x,y)=1(2π)2∬f~(kx,ky) e+i(kxx+kyy) dkx dky.f(x,y) = \frac{1}{(2\pi)^2} \iint \tilde f(k_x,k_y)\,\ee^{+\ii(k_x x+k_y y)}\,\dd k_x\,\dd k_y.(28.9)

把这一步写明的价值在于：后面近远场变换、角谱传播、频率噪声谱以及 RCWA/FDTD/FEM 的谱量定义，都能回到同一组符号，而不会在章节切换时悄悄换掉正负号。

周期介质中的傅里叶展开：为什么卷积会反复出现

正文的 PWEM、RCWA 和 CWT 之所以看起来形式不同，却总会出现相似的矩阵耦合项，根源在于周期介质中的傅里叶卷积。设平面内周期单胞为 AcA_c，倒格矢为 G\gvec。任意满足该周期的标量场 f(r∥)f(\rvec_\parallel) 都可展开为

f(r∥)=∑GfG eiG⋅r∥,fG=1Ac∫Acf(r∥)e−iG⋅r∥ d2r.f(\rvec_\parallel)=\sum_{\gvec} f_{\gvec}\,\ee^{\ii\gvec\cdot\rvec_\parallel}, \qquad f_{\gvec}=\frac{1}{A_c}\int_{A_c} f(\rvec_\parallel)\ee^{-\ii\gvec\cdot\rvec_\parallel}\,\dd^2r.(28.10)

若场满足 Bloch 形式，则还应再乘上相位因子 eik⋅r∥\ee^{\ii\kvec\cdot\rvec_\parallel}。于是任意 Bloch 场都可写为

Ek(r∥,z)=eik⋅r∥∑GEk+G(z)eiG⋅r∥.\E_{\kvec}(\rvec_\parallel,z) = \ee^{\ii\kvec\cdot\rvec_\parallel} \sum_{\gvec} \E_{\kvec+\gvec}(z) \ee^{\ii\gvec\cdot\rvec_\parallel}.(28.11)

关键的是乘积项的展开。若介电函数也按同样方式展开，

ε(r∥)=∑GεGeiG⋅r∥,\varepsilon(\rvec_\parallel) = \sum_{\gvec} \varepsilon_{\gvec} \ee^{\ii\gvec\cdot\rvec_\parallel},

那么 εE\varepsilon\E 的傅里叶系数不再是简单乘法，而是卷积：

(εE)G=∑G′εG−G′EG′.(\varepsilon E)_{\gvec} = \sum_{\gvec'} \varepsilon_{\gvec-\gvec'} E_{\gvec'}.(28.13)

因此周期结构中的本征问题总会变成“一个分量与许多倒格矢分量耦合”的矩阵问题。PWEM 把这个卷积直接写成平面波矩阵；RCWA 把它写进层内本征模矩阵；CWT 则主动只保留少数最重要的耦合项。

若单胞几何只保留最强的四个一级傅里叶分量，那么 CWT 中保留的四波耦合，本质上就是把式 (28.13) 从无限维卷积截断成最小闭合子空间。它的物理解释更强，但严格性当然不如 RCWA 或大基底 PWEM。

开放 Maxwell 问题的函数空间直觉：为什么普通 Hilbert 空间不够

在封闭、无损、厄米问题中，我们习惯把本征模视为某个 Hilbert 空间里的正交基。这个直觉很有用，但它并不能原封不动地搬到 PCSEL 的开放模问题上。根源在于开放系统的共振模满足的是纯出射边界条件，而不是封闭边界条件。对一个复频率模

ω~n=ωn−iγn,γn>0,\tilde{\omega}_n=\omega_n-\ii\gamma_n, \qquad \gamma_n>0,(28.14)

其远场渐近形式可写成

En(r)∼eik~nrrfn(r^),k~n=ω~nc.\E_n(\rvec)\sim \frac{\ee^{\ii \tilde{k}_n r}}{r}\bm{f}_n(\hat{\rvec}), \qquad \tilde{k}_n=\frac{\tilde{\omega}_n}{c}.(28.15)

由于

eik~nr=eiωnr/ce+γnr/c,\ee^{\ii \tilde{k}_n r} = \ee^{\ii \omega_n r/c}\ee^{+\gamma_n r/c},(28.16)

模场若按纯出射条件延拓到无穷远，空间上通常不是衰减而是增长。这并非说模式“物理发散”，而是说它不再属于普通平方可积的束缚态函数空间。开放模更自然的数学身份，是 Green 张量或散射矩阵在复频率平面上的极点态 [Leung et al., 1994] [Alpeggiani et al., 2017]。

这条函数空间直觉非常关键，因为它直接解释了后面三个常见困惑：

封闭腔本征模更像“盒子里的驻波”；开放共振模更像“散射振幅在复频率平面上的极点”。只要把这一步想通，QNM、伴随模和广义归一化就不再是后贴的数学技巧，而是开放问题自然逼出来的语言。

QNM / resonant states 的最小定义

本书中所说的准正规模（quasinormal mode, QNM）或共振态（resonant state），最小可以定义为满足无源 Maxwell 方程且只允许出射波的非平凡解：

L^(ω~n)EnR=0,outgoing boundary condition only.\hat{\mathcal L}(\tilde{\omega}_n)\E_n^{\mathrm{R}}=0, \qquad \text{outgoing boundary condition only}.(28.17)

这里上标 R\mathrm{R} 表示右模（right mode），是为了和稍后出现的左模区分。对开放系统，这个定义与“ω~n\tilde{\omega}_n 是 Green 张量或散射矩阵极点”是等价的 [Sauvan et al., 2013] [Alpeggiani et al., 2017]。两种说法分别对应两种直觉：

这也是为什么 PCSEL 中的开放模分析天然应使用复频率语言。只要目标问题涉及法向辐射、上下出光不对称、quasi-BIC、远场整形或辐射损耗排序，把模式写成 QNM / resonant states 就比把它硬塞回封闭腔本征模更贴近真实物理。

QNM 理论的三个核心数学问题（进阶选读）

准正规模理论在实际应用中面临三个相互关联的数学挑战：(i) 模场的空间增长导致传统L2L^2内积发散；(ii) 非厄米算符的本征模不构成正交系而需采用双正交基；(iii) QNM集合的完备性问题。本节逐一解决这些问题。

1. 发散积分的正则化处理

考虑一维 Fabry-Perot 腔的 QNM解。在腔外区域(∣x∣>L/2|x|>L/2)，模场具有如下渐近形式:

En(x)∼E0eik~n∣x∣,k~n=ω~nc=kn−iκn,E_n(x) \sim E_0 \ee^{\ii \tilde{k}_n |x|}, \qquad \tilde{k}_n = \frac{\tilde{\omega}_n}{c} = k_n - \ii \kappa_n,(28.19)

其中κn>0\kappa_n > 0对应模式的辐射衰减率。将复波矢代入模场表达式可得，当x→∞x\to\infty时，∣En(x)∣2∼e2κnx|E_n(x)|^2 \sim \ee^{2\kappa_n x}指数增长，导致传统的L2L^2范数积分发散:

∫−∞∞∣En(x)∣2dx=∞.\int_{-\infty}^{\infty} |E_n(x)|^2 \dd x = \infty.(28.20)

为解决这一问题，现代 QNM理论采用以下几种等价的正则化方案:

(a) 体积极限正则化 (Volume regularization): 先将系统限制在有限体积VL=[−L,L]V_L=[-L,L]内计算，最后取L→∞L\to\infty极限并提取有限部分(finite part)。具体而言，将发散积分写为:

∫−LL∣En(x)∣2dx=Ae2κnL+B(L),\int_{-L}^{L} |E_n(x)|^2 \dd x = A \ee^{2\kappa_n L} + B(L),

其中AA是发散的系数，B(L)B(L)是随LL缓变的部分。定义正则化后的内积为去掉指数发散项后的剩余部分:

⟨En,Em⟩reg:=lim⁡L→∞[∫−LLEn∗(x)Em(x)dx−δnmAn2κne2κnL].\langle E_n, E_m \rangle_{\rm reg} := \lim_{L\to\infty} \left[ \int_{-L}^{L} E_n^*(x) E_m(x) \dd x - \delta_{nm} \frac{A_n}{2\kappa_n}\ee^{2\kappa_n L} \right].(28.22)

(b) 解析延拓法 (Analytic continuation): 将 QNM视为复波数平面上的亚纯函数，通过将实轴积分路径旋转到复平面的适当位置来实现收敛。设k=k~n+seiθk = \tilde{k}_n + s \ee^{\ii\theta}，选择旋转角θ\theta使得被积函数沿新路径指数衰减。这种方法与留数定理自然衔接。

(c) PML 正则化: 引入完美匹配层坐标变换x→x~(x)=x+if(x)x \to \tilde{x}(x) = x + \ii f(x)，其中f(x)f(x)在计算域外部单调递增。在此变换下，原发散的被积函数变为:

∫−∞∞∣En(x)∣2dx⟶∫−∞∞∣En(x~)∣2dx~dxdx,\int_{-\infty}^{\infty} |E_n(x)|^2 \dd x \quad \longrightarrow \quad \int_{-\infty}^{\infty} |E_n(\tilde{x})|^2 \frac{\dd\tilde{x}}{\dd x} \dd x,

由于ℑ[k~nx~]<0\Im[\tilde{k}_n \tilde{x}] < 0在外区成立，积分自动收敛。这是商业软件中最常用的实现方式。

这三种方法在数学上是等价的，都给出了相同的广义归一化结果 [Kristensen et al., 2015] [Sauvan et al., 2013]。

2. 双正交基的严格构造

对于非厄米 Maxwell 算符Θ^\hat{\Theta}，其右本征模{∣Rn⟩}\{|R_n\rangle\}满足:

Θ^∣Rn⟩=ω~n2∣Rn⟩,\hat{\Theta}|R_n\rangle = \tilde{\omega}_n^2 |R_n\rangle,

但这些右模之间并不正交:⟨Rm∣Rn⟩≠0\langle R_m|R_n\rangle \neq 0 (m≠nm\neq n)。为了构建有效的模态展开，必须引入伴随算符Θ^†\hat{\Theta}^\dagger的左本征模{∣Ln⟩}\{|L_n\rangle\}:

Θ^†∣Ln⟩=(ω~n2)∗∣Ln⟩.\hat{\Theta}^\dagger|L_n\rangle = (\tilde{\omega}_n^2)^* |L_n\rangle.(28.25)

关键结论是：左右模构成双正交系(biorthogonal system):

⟨Lm∣Rn⟩=Nnδmn,\langle L_m|R_n\rangle = N_n \delta_{mn},(28.26)

其中NnN_n是依赖于归一化约定的常数。这里的内积应采用前述正则化方案之一来计算。

双正交系的几何图像

在非厄米系统中，左模和右模不再是彼此的共轭转置，而是两个独立的向量集。右模{∣Rn⟩}\{|R_n\rangle\}张成状态空间，而左模{∣Ln⟩}\{|L_n\rangle\}作为其对偶基底(dual basis)。任意物理态∣Ψ⟩|\Psi\rangle按右模展开时，展开系数由对应的左模投影给出:cn=⟨Ln∣Ψ⟩/Nnc_n=\langle L_n|\Psi\rangle/N_n。这种关系类似于协变矢量和逆变矢量——通过配对缩并得到标量。

3. 完备性关系与模态展开

若 QNM集合{∣Rn⟩}\{|R_n\rangle\}是完备的，则任意物理态∣Ψ(t)⟩|\Psi(t)\rangle可以展开为:

∣Ψ(t)⟩=∑ncn(t)∣Rn⟩,cn(t)=⟨Ln∣Ψ(t)⟩⟨Ln∣Rn⟩.|\Psi(t)\rangle = \sum_n c_n(t) |R_n\rangle, \qquad c_n(t) = \frac{\langle L_n|\Psi(t)\rangle}{\langle L_n|R_n\rangle}.(28.27)

完备性的数学表述是闭包关系(closure relation):

∑n∣Rn⟩⟨Ln∣⟨Ln∣Rn⟩=I^,\sum_n \frac{|R_n\rangle\langle L_n|}{\langle L_n|R_n\rangle} = \hat{I},(28.28)

其中I^\hat{I}是单位算符。

然而，对于开放系统，QNM的完备性并非总是成立。Leung等人 [Leung et al., 1994]证明了以下重要结论:

从 Green 函数的角度看，完备性等价于说 Green 算符可以用 QNM展开:

G↔(r,r′;ω)=∑nEn(r)⊗EnL(r′)2ω~n(ω−ω~n),\overleftrightarrow{G}(\rvec,\rvec';\omega) = \sum_n \frac{\E_n(\rvec)\otimes\E_n^{\rm L}(\rvec')}{2\tilde{\omega}_n(\omega-\tilde{\omega}_n)},(28.29)

其中EnL\E_n^{\rm L}是左模的电场分量。这个展开式在远离分支切割(branch cut)的区域一致收敛，而在接近连续谱阈值时需要额外小心处理 [Sauvan et al., 2013] [Kristensen et al., 2015]。

实际应用建议

对于 PCSEL仿真，推荐采用以下工作流程:

  1. 使用 COMSOL/Lumerical 的 QNM求解器获取前NN个主导模式的复频率ω~n\tilde{\omega}_n和场分布En(r)\E_n(\rvec)

  2. 在同一平台上计算对应的左模EnL\E_n^{\rm L}(或通过互易性关系间接获得)

  3. 验证模态截断误差：先计算双正交投影系数

    cn=⟨Ln∣Ψtarget⟩⟨Ln∣Rn⟩,c_n=\frac{\langle L_n|\Psi_{\rm target}\rangle}{\langle L_n|R_n\rangle},

    再以重构残差定义误差

    η=∥Ψtarget−∑n=1NcnRn∥W∥Ψtarget∥W\eta= \frac{\left\|\Psi_{\rm target}-\sum_{n=1}^{N} c_n R_n\right\|_{W}} {\|\Psi_{\rm target}\|_{W}}

  4. 若η<1%\eta < 1\%,则认为当前 QNM基组足够用于后续的微扰分析和耦合计算

光子学中的 Green 函数方法简介

前面已经从本征模的角度讨论了开放系统的 QNM 框架。但在现代纳米光子学中，还有一个更加统一的语言——Green 函数方法。它不仅包含了本征模的所有信息，还能自然地处理光源激发、散射问题和自发辐射耦合等物理过程。

从标量波动方程的 Green 函数出发

考虑最简单的标量 Helmholtz 方程：

[∇2+k02ε(r)]ψ(r)=−s(r),\left[\nabla^2 + k_0^2\varepsilon(\rvec)\right]\psi(\rvec) = -s(\rvec),(28.30)

其中s(r)s(\rvec)是源项。我们希望找到一个"单位点源的响应"G(r,r′)G(\rvec,\rvec')，使得对任意源分布都能通过叠加原理写出解：

ψ(r)=∫G(r,r′)s(r′) d3r′.\psi(\rvec) = \int G(\rvec,\rvec') s(\rvec')\,\dd^3 r'.(28.31)

这就要求GG满足：

[∇2+k02ε(r)]G(r,r′)=−δ(r−r′).\left[\nabla^2 + k_0^2\varepsilon(\rvec)\right]G(\rvec,\rvec') = -\delta(\rvec-\rvec').(28.32)

这就是 Green 函数的定义：它是算符L^=∇2+k02ε(r)\hat{L}=\nabla^2+k_0^2\varepsilon(\rvec)的逆算符的核（kernel）。

因果性与推迟 Green 函数

对于时谐问题，还需要指定边界条件来确定 Green 函数的唯一性。物理上合理的要求是因果性：响应不能超前于激励。在频域中，这等价于要求当r→∞r\to\infty时只有出射波：

G(r,r′)∼eik0nr4πr,r→∞.G(\rvec,\rvec') \sim \frac{\ee^{\ii k_0 n r}}{4\pi r}, \qquad r\to\infty.(28.33)

这个条件被称为 Sommerfeld 辐射条件。满足它的 Green 函数称为推迟 Green 函数（retarded Green's function）。

Dyson 方程及其物理含义

设我们把介电常数写成背景加扰动：ε(r)=εbg(r)+Δε(r)\varepsilon(\rvec)=\varepsilon_{\mathrm bg}(\rvec)+\Delta\varepsilon(\rvec)。相应地，Green 函数也可以分解为背景部分G0G_0和散射修正ΔG\Delta G。将它们代入定义式并整理，可得到著名的 Dyson 方程：

G(r,r′)=G0(r,r′)+∫G0(r,r′′) k02Δε(r′′) G(r′′,r′) d3r′′.G(\rvec,\rvec') = G_0(\rvec,\rvec') + \int G_0(\rvec,\rvec'')\,k_0^2\Delta\varepsilon(\rvec'')\,G(\rvec'',\rvec')\,\dd^3 r''.(28.35)

这个方程的物理图像极其清晰：总响应等于直接响应加上所有可能的单次/多次散射路径的总和。如果把这个方程迭代展开，就得到 Born 级数：

G=G0+G0(k02Δε)G0+G0(k02Δε)G0(k02Δε)G0+⋯=G0+G0TG0,\begin{aligned}G &= G_0 + G_0(k_0^2\Delta\varepsilon)G_0 + G_0(k_0^2\Delta\varepsilon)G_0(k_0^2\Delta\varepsilon)G_0 + \cdots \\ &= G_0 + G_0 T G_0,\end{aligned}

其中TT算符总结了所有阶次的散射过程。

Dyson 方程的核心思想是"递归"：要知道总的散射响应，可以先问"第一次散射发生在哪里"，然后再问"散射之后的传播又如何"。这种自相似的结构正是 Green 函数方法强大之处的源泉。

矢量 Green 张量与 Maxwell 方程

对真实的电磁场问题，我们需要处理的是矢量波动方程。此时 Green 函数推广为一个二阶张量G‾(r,r′)\overline{\bm{G}}(\rvec,\rvec')（也称为 dyadic Green's function），满足：

∇×∇×G‾(r,r′)−k02ε(r)G‾(r,r′)=I‾ δ(r−r′),\nabla\times\nabla\times\overline{\bm{G}}(\rvec,\rvec') - k_0^2\varepsilon(\rvec)\overline{\bm{G}}(\rvec,\rvec') = \overline{\bm{I}}\,\delta(\rvec-\rvec'),(28.38)

其中I‾\overline{\bm{I}}是单位并矢（unit dyad）。电场与电流源的关系为：

E(r)=iωμ0∫G‾(r,r′)⋅J(r′) d3r′.\E(\rvec) = \ii\omega\mu_0 \int \overline{\bm{G}}(\rvec,\rvec')\cdot\J(\rvec')\,\dd^3 r'.(28.39)

注意这里的因果关系：偶极子源J\J产生场E\E，中间的传递函数就是 Green 张量。

周期结构中的 dyadic Green 函数与 Bloch-Floquet 展开

对 PCSEL 这样的周期结构，Green 张量也必须满足 Bloch 定理。具体来说，可以将其展开为：

G‾(r,r′;k)=∑nEnR(r;k)⊗EnL(r′;k)2ω~n(k)[ω~−ω~n(k)],\overline{\bm{G}}(\rvec,\rvec';\kvec) = \sum_n \frac{\E_n^{\mathrm R}(\rvec;\kvec)\otimes\E_n^{\mathrm L}(\rvec';\kvec)}{2\tilde{\omega}_n(\kvec)[\tilde{\omega}-\tilde{\omega}_n(\kvec)]},(28.40)

其中EnR/L\E_n^{\mathrm{R/L}}分别是右模和左模，ω~n\tilde{\omega}_n是复本征频率。这个表达式有几个关键特征：

  1. 分母中的极点结构清楚地表明：Green 函数的奇点对应系统的本征模；

  2. 分子中的双线性形式反映了非厄米系统的双正交特性；

  3. 求和覆盖了所有离散模和连续谱贡献（后者在某些情况下需要写成积分）。

式 (28.40) 的严格推导需要小心处理收敛性和完备性问题。对开放系统，由于 QNM 的空间增长特性，这个展开在某些区域可能条件收敛甚至发散。实用的解决方案包括使用 PML 正则化或采用适当的解析延拓技术 [Sauvan et al., 2013] [Kristensen et al., 2015]。

Green 函数极点与准正则模的对应关系

现在我们回到之前留下的问题：为什么说 QNM 是 Green 函数的极点？从 式 (28.40)可以看出，当探测频率ω~\tilde{\omega}趋近某个本征值ω~n\tilde{\omega}_n时，对应的项会发散。在复频率平面上，这些发散的点就是极点（pole）。

更严格地说，可以在极点附近做 Laurent 展开：

G‾(ω~)≈Resnω~−ω~n+regular part,\overline{\bm{G}}(\tilde{\omega}) \approx \frac{\text{Res}_n}{\tilde{\omega}-\tilde{\omega}_n} + \text{regular part},(28.41)

其中留数Resn\text{Res}_n正好正比于EnR⊗EnL\E_n^{\mathrm R}\otimes\E_n^{\mathrm L}。这就建立了两种语言之间的精确对应：

tabularll

QNM 语言 & Green 函数语言

本征频率 ω~n\tilde{\omega}_n & 极点位置

模场 EnR/L\E_n^{\mathrm{R/L}} & 留数的因子分解

归一化条件 & 留数的归一化

tabular

Purcell 因子的 Green 函数表达

Green 函数方法的一个重要应用是计算自发辐射增强因子（Purcell factor）。考虑一个位于r0\rvec_0的电偶极子p\bm{p}，其辐射功率可通过 Green 张量表示为：

P=ω2ℑ[p∗⋅E(r0)]=ω2μ02∣p∣2 p^⋅ℑ[G‾(r0,r0)]⋅p^,P = \frac{\omega}{2}\Im\left[\bm{p}^*\cdot\E(\rvec_0)\right] = \frac{\omega^2\mu_0}{2}|\bm{p}|^2\,\hat{\bm{p}}\cdot\Im\left[\overline{\bm{G}}(\rvec_0,\rvec_0)\right]\cdot\hat{\bm{p}},(28.42)

其中p^\hat{\bm{p}}是偶极矩方向的单位矢量。相对于真空中的辐射功率P0P_0，Purcell 因子定义为：

FP(r0,p^,ω)=PP0=6πcω p^⋅ℑ[G‾(r0,r0)]⋅p^.F_P(\rvec_0,\hat{\bm{p}},\omega) = \frac{P}{P_0} = \frac{6\pi c}{\omega}\,\hat{\bm{p}}\cdot\Im\left[\overline{\bm{G}}(\rvec_0,\rvec_0)\right]\cdot\hat{\bm{p}}.(28.43)

这个公式的强大之处在于：它不需要显式地知道任何本征模，只需要计算源点处的 Green 张量即可。当然，如果将 式 (28.40)代入，也能得到用 QNM 展开的形式：

FP(r0,p^,ω)≈∑nℜ ⁣[Anω−ω~n],F_P(\rvec_0,\hat{\bm p},\omega) \approx \sum_n \Re\!\left[\frac{A_n}{\omega-\tilde{\omega}_n}\right],(28.44)

其中

An=6πcω (p^ ⁣⋅ ⁣EnR(r0))(p^ ⁣⋅ ⁣EnL(r0))2ω~n⟨EnL,EnR⟩A_n= \frac{6\pi c}{\omega}\, \frac{\bigl(\hat{\bm p}\!\cdot\!\E_n^{\mathrm R}(\rvec_0)\bigr)\bigl(\hat{\bm p}\!\cdot\!\E_n^{\mathrm L}(\rvec_0)\bigr)} {2\tilde{\omega}_n\langle \E_n^{\mathrm L},\E_n^{\mathrm R}\rangle}

是由左右模重叠与归一化共同决定的复残数。仅在单模、弱重叠、高 QQ 近似下，式 (28.44) 才可进一步退化为 Lorentzian 峰形；一般开放系统可出现非 Lorentz 线形和模间干涉项。

伴随模、双正交与广义归一化

一旦算子不再自伴，右模本身通常不足以完成投影与微扰分析。对第一次接触这套语言的读者，最稳妥的直觉是：先把左模看成做投影时的测试函数，而不要先把它想成“另一支也要单独观察到的物理场”。

伴随算子是什么？

选定一个内积 ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle 后，伴随算子由“分部积分后的移项”定义：

⟨u,L^v⟩=⟨L^Au,v⟩+边界项。\langle u,\hat{\mathcal L}v\rangle = \langle \hat{\mathcal L}^{\mathrm{A}}u, v\rangle + \text{边界项}。(28.45)

若边界项为零，则 L^A\hat{\mathcal L}^{\mathrm{A}} 就是把 L^\hat{\mathcal L} “搬到左边”得到的算子；封闭厄米问题中更进一步有 L^A=L^\hat{\mathcal L}^{\mathrm{A}}=\hat{\mathcal L}。

为什么开放系统必须区分左右模？

开放边界使边界项不再自动消失，因此左模与右模必须分开定义。更一般的做法是同时引入右模与伴随左模：

L^(ω~n)EnR=0,L^A(ω~n)EnL=0.\hat{\mathcal L}(\tilde{\omega}_n)\E_n^{\mathrm{R}}=0, \qquad \hat{\mathcal L}^{\mathrm{A}}(\tilde{\omega}_n)\E_n^{\mathrm{L}}=0.(28.46)

其中 L^A\hat{\mathcal L}^{\mathrm{A}} 是伴随算子，EnL\E_n^{\mathrm{L}} 是左模。对一般非厄米问题，真正起作用的是左模和右模形成的双正交结构，而不是右模与自身的厄米内积。也就是说，开放系统中的“正交”更接近

⟨EmL,EnR⟩∝δmn,\left\langle \E_m^{\mathrm{L}},\E_n^{\mathrm{R}}\right\rangle \propto \delta_{mn},(28.47)

而不是封闭腔里熟悉的 ⟨Em,En⟩ε=δmn\langle \E_m,\E_n\rangle_\varepsilon=\delta_{mn}。

对 Maxwell 开放问题，进一步的细节取决于介质是否互易、是否色散，以及数值上是否通过 PML 处理开放边界。对许多互易系统，一个常用的广义 QNM 归一化形式是

⟨ ⁣⟨Ψm,Ψn⟩ ⁣⟩=∫Ω+PML[Em⋅∂(ωε)∂ωEn−Hm⋅∂(ωμ)∂ωHn]dV,\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,(28.48)

其中 Ψn=(En,Hn)\Psi_n=(\E_n,\Hf_n)，且积分区域包含物理区与 PML [Sauvan et al., 2013] [Kristensen et al., 2015]。这个表达式的意义不在于要求读者把它背下来，而在于明确三件事： enumerate 广义归一化里经常不出现普通的复共轭； 材料色散会直接进入归一化，而不是只出现在“材料模型”章节； PML 并不只是吸收边界的数值技巧，在 QNM 归一化里它还承担把开放问题转成可积配对的数学角色。 enumerate

在封闭、无损、无色散极限下，式 (28.48) 会退回更熟悉的能量型归一化；但对真实 PCSEL 的开放模，若继续坚持只用

⟨Em,En⟩ε=∫Vε(r)Em∗(r)⋅En(r) dV=δmn.\langle \E_m,\E_n\rangle_{\varepsilon} = \int_V \varepsilon(\rvec)\E_m^*(\rvec)\cdot\E_n(\rvec)\,\dd V = \delta_{mn}.(28.49)

往往就会在归一化、模耦合和微扰分母上埋下系统错误。也正因此，本书正文只有在明确处于封闭或弱开放近似下，才会继续使用式 (28.49) 这类厄米归一化。

这里有两条最容易被忽略、却又最关键的原因。第一，对纯出射 QNM 取复共轭后，会把出射边界条件翻成入射边界条件，因此 E∗\E^\ast 本身通常不是同一个开放本征问题的合法测试函数。第二，开放模在无穷远处的空间延拓并不平方可积，若仍执意把它塞进普通厄米内积，分子、分母和归一化常数都可能一起失真。正因为这两点，开放系统一阶微扰不能机械照搬封闭腔公式。

这意味着一个非常重要的实务结论：若两个不同求解器采用了不同归一化方式，不能直接比较它们输出的“场强大小”，必须比较与归一化无关的不变量，或者先把它们转换到同一内积和归一化框架下。正文中大量出现的重叠因子之所以用分式定义，也是为了让归一化常数在分子分母中相消。

例如，量子阱重叠因子常写成

ΓQW=∫VQWε∣E∣2 dV∫Vε∣E∣2 dV,\Gamma_{\mathrm{QW}} = \frac{\int_{V_{\mathrm{QW}}} \varepsilon |\E|^2\,\dd V} {\int_V \varepsilon |\E|^2\,\dd V},(28.50)

它之所以稳健，就在于不依赖模整体振幅的选取。

复频率、衰减率与 QQ：不要混淆振幅寿命和能量寿命

开放腔模常写成复频率形式

ω~=ωr−iγ,γ>0.\tilde{\omega}=\omega_r-\ii\gamma, \qquad \gamma>0.(28.51)

代回时间因子后，场的时间依赖为

e−iω~t=e−iωrte−γt.\ee^{-\ii\tilde{\omega} t} = \ee^{-\ii\omega_r t}\ee^{-\gamma t}.

因此场振幅的衰减时间常数是 1/γ1/\gamma。但能量与场振幅平方成正比，所以腔内能量按

U(t)∝e−2γtU(t) \propto \ee^{-2\gamma t}

衰减，对应的能量寿命为

τp=12γ.\tau_p=\frac{1}{2\gamma}.(28.54)

这时品质因子定义为每个周期储能与损耗的比值，可写成

Q=ωrτp=ωr2γ.Q = \omega_r\tau_p = \frac{\omega_r}{2\gamma}.(28.55)

这个推导之所以需要在附录里明确写出，是因为很多实际文稿会把 1/γ1/\gamma 直接称为“光子寿命”，或者把振幅衰减率和能量衰减率混用。对被动光学问题，这种混淆可能只差一个因子 2；但一旦进入阈值条件、线宽估计和速率方程接口，这个因子 2 就会在不同模块之间制造系统偏差。

互易性：为什么许多耦合公式其实来自同一个定理

正文中的耦合模理论、微扰公式和灵敏度分析，背后都有一个共同工具，即 Lorentz 互易性。设两组场 (E1,H1)(\E_1,\Hf_1) 和 (E2,H2)(\E_2,\Hf_2) 在同一线性、时不变、局域介质中分别满足

ε=εT,μ=μT,\varepsilon=\varepsilon^{\mathsf T}, \qquad \mu=\mu^{\mathsf T},(28.56)

所对应的 Maxwell 方程。式 (28.56) 是这一节最关键的适用范围前提：它排除了磁光、旋光、显式时变调制和其他非互易介质。只有在这一前提下，下面的 Lorentz 互易性推导才按最熟悉的对称形式成立。

本节默认介质张量满足 式 (28.56)。若器件包含磁光偏置、时空调制或其他非互易机制，则互易关系必须改写，不能直接把本节公式平移过去。

∇×E1=+iωμH1,∇×H1=J1−iωεE1,∇×E2=+iωμH2,∇×H2=J2−iωεE2.\begin{aligned}\nabla\times\E_1 &= +\ii\omega\mu\Hf_1, & \nabla\times\Hf_1 &= \J_1-\ii\omega\varepsilon\E_1, \\ \nabla\times\E_2 &= +\ii\omega\mu\Hf_2, & \nabla\times\Hf_2 &= \J_2-\ii\omega\varepsilon\E_2.\end{aligned}

对第一式与第二组场做混合乘积，再用向量恒等式

∇⋅(E1×H2−E2×H1)=H2⋅(∇×E1)−E1⋅(∇×H2)−H1⋅(∇×E2)+E2⋅(∇×H1),\begin{aligned} \nabla\cdot(\E_1\times\Hf_2-\E_2\times\Hf_1) &= \Hf_2\cdot(\nabla\times\E_1) - \E_1\cdot(\nabla\times\Hf_2) \\ &\quad - \Hf_1\cdot(\nabla\times\E_2) + \E_2\cdot(\nabla\times\Hf_1), \end{aligned}(28.59)

可得到

∇⋅(E1×H2−E2×H1)=E2⋅J1−E1⋅J2.\nabla\cdot(\E_1\times\Hf_2-\E_2\times\Hf_1) = \E_2\cdot\J_1-\E_1\cdot\J_2.

对体积积分并使用散度定理，就得到互易性形式

∫∂V(E1×H2−E2×H1)⋅n^ dS=∫V(E2⋅J1−E1⋅J2) dV.\int_{\partial V}(\E_1\times\Hf_2-\E_2\times\Hf_1)\cdot\hat{\bm n}\,\dd S = \int_V (\E_2\cdot\J_1-\E_1\cdot\J_2)\,\dd V.(28.61)

式 (28.61) 的价值在于，它为“某个局域扰动为什么会以某种重叠积分方式耦合到目标模”提供了统一解释。很多看起来不同的耦合系数公式，本质上都能追溯到互易性与边界条件。

互易性并非无条件成立。若介质存在显著非互易效应、时变调制、磁光非对称张量或非局域响应，则式 (28.61) 需要修改，甚至可能失效。将其直接套用到所有有源介质问题上是不严谨的。

一阶微扰公式：什么时候它是定理，什么时候它只是工程近似

现在可以回到正文中最关键的开放系统数学缺口：非厄米一阶扰动到底应怎样写。若先从最一般的角度出发，把本征问题写成

L^(ω~,p)ER=0,\hat{\mathcal L}(\tilde{\omega},p)\E^{\mathrm{R}}=0,(28.62)

其中 pp 是几何、材料或工作点参数，那么对小扰动 p→p+δpp\to p+\delta p 做一阶展开，可得

[∂L^∂ω∣nδω~+δL^]EnR+L^n δEnR=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.(28.63)

再引入左模投影，就得到开放系统的一阶扰动骨架

δω~n=−⟨EnL, δL^ EnR⟩⟨EnL, ∂L^∂ω∣nEnR⟩.\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 }.(28.64)

这就是开放系统里真正普适的一阶公式。它比封闭腔公式更“硬”，原因不在于形式更复杂，而在于它保留了两件开放系统中不能被省略的信息：左/右模的双正交结构，以及算子对频率的显式依赖。

若系统退回封闭、无损、弱开放极限，且左模与右模在厄米内积下重合，同时忽略显式色散，那么式 (28.64) 就会退化为熟悉的频移公式

δωω≈−12∫Vδε(r)∣E(r)∣2 dV∫Vε(r)∣E(r)∣2 dV.\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}.(28.65)

它的物理意义仍然非常直观：若正的电场能量主要落在折射率增大的区域，则共振频率通常下降。正文中关于热漂移、材料修正和几何微扰的许多定性判断，都可用这一公式快速给出方向感。

但这里必须准确区分三种层级。

第一层，严格定理层级：式 (28.64)。只要左模、右模、伴随算子和归一化框架是自洽的，它就是开放系统一阶扰动的正确起点。

第二层，弱开放工程层级：式 (28.65)。当泄露较弱、模式局域较强、色散不主导，且扰动主要发生在腔内高场区域时，它通常是非常有用的趋势公式。

第三层，超出微扰适用层级：近简并模强混合、模式身份交换、扰动直接重构辐射边界、材料色散和吸收显著、或电热回写已经把问题从“局部小改动”推成“空间分布重构”。此时无论厄米还是非厄米一阶公式都不应继续扛主问题，必须回到完整本征问题或全耦合重算。

若扰动是附加小吸收 δε′′>0\delta\varepsilon''>0，则同样的骨架会给出复频率虚部修正，也就是额外损耗，而非实部折射率。于是，从数学上看，“热漂移改频率”和“掺杂吸收改线宽”其实是同一套开放系统微扰公式的不同投影：前者更看重 ℜ δω~\Re\,\delta\tilde{\omega}，后者更看重 ℑ δω~\Im\,\delta\tilde{\omega}。

若某个 PCSEL 候选模已经处于 quasi-BIC 附近，那么轻微孔形扰动往往首先重排的是辐射通道相消条件，而不是腔内平均折射率。此时若只用式 (28.65) 看实部频移，可能会得到“改动很小”的错觉；而真正决定阈值排序的，恰恰是虚部损耗的放大变化。这里就必须至少回到式 (28.64) 的层级，甚至直接做全波重算。

这些数学工具在 PCSEL 仿真中怎么实际使用

附录的价值，最终还是要回到建模实践。对本书关心的 PCSEL 问题，这些工具至少有四类直接用途。

第一，用时间因子和复频率关系检查损耗与增益符号是否在不同模块间保持一致。若 RCWA、FDTD 和器件增益模型对介质虚部的正负含义不一致，必须立即停下来修正。

第二，用 Fourier-Bloch 展开理解数值方法之间的联系。这样在比较 PWEM、RCWA 和 CWT 时，读者就知道它们差别主要在截断和投影，而不是在“互相矛盾的物理定律”。

第三，用归一化无关的重叠积分比较不同求解器结果，避免把场强数值本身的差异误判为物理差异；若必须比较开放模的绝对模幅或灵敏度，则应先确认所用的是同一种左/右模与广义归一化框架。

第四，用复频率与一阶微扰公式快速做 sanity check。例如，温度升高若使波导折射率整体增加，那么目标共振一般应向长波漂移；若仿真却给出相反趋势，就要优先检查材料模型和符号约定。

第五，对强开放候选模，不要只记录 ωr\omega_r 或 QQ，而应同时记录使用的归一化、是否采用 QNM / PML 框架、以及微扰判断使用的是厄米代理还是非厄米左/右模公式。否则后续章节即便复用了同一个符号，也未必在说同一层数学对象。

从 CWT 本征矢量提取远场偏振：方位角、椭圆率与 Stokes 参量

本节给出一组可直接写进代码的公式：如何从 第 10 章、第 11 章 的四波 CWT 本征矢量读出远场偏振态。 在实际器件设计中，除了关心模式的辐射常数和品质因子，还需要直接从 CWT 解中提取远场的偏振特性。设第 jj 个超模的本征矢量为

Vj=[cx,jsx,jcy,jsy,j]T,\bm V_j = \begin{bmatrix} c_{x,j} & s_{x,j} & c_{y,j} & s_{y,j} \end{bmatrix}^{\mathsf T},

其中下标 x,yx,y 分别表示沿 ±x\pm x 和 ±y\pm y 方向传播的基本波对。在最小四波模型中，远场主要由前向波 cx,jc_{x,j} 和 cy,jc_{y,j} 贡献（后向波 sx,j,sy,js_{x,j},s_{y,j} 主要形成驻波反馈）。因此，远场电场的横向分量可近似写为

Eff≈cx,j x^+cy,j y^.\bm E_{\mathrm{ff}} \approx c_{x,j}\,\hat{\bm x} + c_{y,j}\,\hat{\bm y}.

由此可直接定义三个关键偏振参数：

1. 偏振方位角 (azimuth angle)

线偏振主轴相对于 xx轴的夹角为

ψj=12atan2⁡ ⁣(2 Re ⁣[cx,jcy,j∗],  ∣cx,j∣2−∣cy,j∣2),ψj∈(−π2,π2].\psi_j = \frac{1}{2}\operatorname{atan2}\!\left( 2\,\mathrm{Re}\!\left[c_{x,j}c_{y,j}^{\ast}\right],\; |c_{x,j}|^2-|c_{y,j}|^2 \right), \qquad \psi_j \in \left(-\frac{\pi}{2},\frac{\pi}{2}\right].(28.68)

特别地，当相位差接近零时，式 (28.68) 会退化到熟悉的幅值比公式；当 ∣cx,j∣=∣cy,j∣|c_{x,j}|=|c_{y,j}| 且相位差为零时，ψj=45∘\psi_j=45^\circ，对应对角线偏振；当 cy,j=0c_{y,j}=0 时，ψj=0∘\psi_j=0^\circ，对应纯 xx 偏振。

2. 椭圆率 (ellipticity)

若考虑复振幅的相对相位，先定义椭圆角 χj\chi_j：

χj=12arcsin⁡ ⁣(2 Im ⁣[cx,j∗cy,j]∣cx,j∣2+∣cy,j∣2),ej=tan⁡χj,ej∈[−1,1],\chi_j = \frac{1}{2}\arcsin\!\left( \frac{2\,\mathrm{Im}\!\left[c_{x,j}^{\ast}c_{y,j}\right]}{|c_{x,j}|^2 + |c_{y,j}|^2} \right), \qquad e_j = \tan\chi_j, \qquad e_j \in [-1, 1],(28.69)

其中 ej=0e_j=0 对应完全线偏振，∣ej∣=1|e_j|=1 对应圆偏振。

3. 斯托克斯参量 (Stokes parameters)

更一般地，可以用归一化斯托克斯参量完整刻画偏振态：

s0=∣cx,j∣2+∣cy,j∣2,s1=∣cx,j∣2−∣cy,j∣2s0,s2=2Re(cx,j∗cy,j)s0,s3=2Im(cx,j∗cy,j)s0.\begin{aligned}s_0 &= |c_{x,j}|^2 + |c_{y,j}|^2, \notag\\ s_1 &= \frac{|c_{x,j}|^2 - |c_{y,j}|^2}{s_0}, \notag\\ s_2 &= \frac{2\mathrm{Re}(c_{x,j}^* c_{y,j})}{s_0}, \notag\\ s_3 &= \frac{2\mathrm{Im}(c_{x,j}^* c_{y,j})}{s_0}. \end{aligned}(28.70)

这四个参量满足s12+s22+s32=1s_1^2+s_2^2+s_3^2=1（完全偏振光），且在庞加莱球上构成一个单位向量。结合 第 8 章 中的对称性分析，可以预测不同模式在不同对称破缺下的轨迹演化。

群论工具箱：C4vC_{4v} 点群完整参考

本节把 第 8 章 只按需引用的群论事实集中成一份完整参考：生成元、特征标表、直积表、选择规则与 C4v→C2vC_{4v}\to C_{2v} 分支律。它服务于查阅而非线性阅读。

生成元与群元素：C4vC_{4v} 由两个生成元生成： itemize C4C_4: 绕 zz 轴旋转 90∘90^\circ（逆时针方向为正） σv\sigma_v: 关于 xzxz 平面（即 y=0y=0 平面）的镜面对称 itemize 全部 8 个群元素为：{E,C4,C42=C2,C43,σv,x,σv,y,σd,xy,σd,−xy}\{E, C_4, C_4^2=C_2, C_4^3, \sigma_{v,x}, \sigma_{v,y}, \sigma_{d,xy}, \sigma_{d,-xy}\}。其中： itemize EE: 恒等操作 C4,C43C_4, C_4^3: ±90∘\pm 90^\circ 旋转（这两者属于同一类） C2=C42C_2=C_4^2: 180∘180^\circ 旋转 σv,x,σv,y\sigma_{v,x}, \sigma_{v,y}: 垂直镜面（分别关于 yzyz 和 xzxz 平面） σd,xy,σd,−xy\sigma_{d,xy}, \sigma_{d,-xy}: 对角镜面（分别关于 y=xy=x 和 y=−xy=-x 平面） itemize

特征标表回顾与应用：

tabularc|ccccc|c C4vC_{4v} & EE & 2C42C_4 & C2C_2 & 2σv2\sigma_v & 2σd2\sigma_d & 典型基函数

A1A_1 & 1 & 1 & 1 & 1 & 1 & zz, x2+y2x^2+y^2, f(r)f(r)

A2A_2 & 1 & 1 & 1 & -1 & -1 & RzR_z

B1B_1 & 1 & -1 & 1 & 1 & -1 & x2−y2x^2-y^2

B2B_2 & 1 & -1 & 1 & -1 & 1 & xyxy

EE & 2 & 0 & -2 & 0 & 0 & (x,y)(x,y), (Rx,Ry)(R_x,R_y)

tabular

表中每一行的物理含义：当对该 irrep 的基函数施加相应的对称操作时，函数乘以对应的特征标值（+1+1 保持不变，−1-1 改变符号）。例如，A1A_1 表示在所有 8 个操作下都不变——这是"完全对称"的表示；而 EE 是二维表示，其特征标为 0 意味着两个基函数在该操作下互相交换或线性组合。

PCSEL 中的模式分类实操步骤： enumerate 提取近场分布: 从仿真中获得 ∣E(x,y)∣2|\mathbf{E}(x,y)|^2 或特定分量的实部/虚部 应用对称操作: 对场分布依次施加 C4C_4 旋转和四个镜面反射 记录响应: 每次操作后，判断场是 (a) 不变 (+1)、(b) 反相 (-1)、还是 (c) 变成另一个正交模式（二维表示的特征） 查表比对: 将得到的特征标序列与特征标表逐行对比，确定归属

快捷技巧：实际上只需要测试三个关键操作即可唯一确定大多数情况：C4C_4（区分单重态/双重态）、σv\sigma_v（区分 A/BA/B 型）、σd\sigma_d（区分下标 1/2）。 enumerate

基于上述流程，PCSEL 模式可分为两类： itemize 单重态 (A1,A2,B1,B2A_1, A_2, B_1, B_2): 非简并模式，电场具有确定的宇称变换性质。这些模式往往是 BIC 候选（若不与辐射通道 irrep 匹配）。 双重态 (EE): 二重简并模式，对应一对正交偏振态。在无扰动情况下能量严格相等，是典型的亮模候选。 itemize

选择规则的群论表述：考虑任意微扰算符 O^\hat{O}（可以是介电常数扰动、表面形变、电流注入等）诱导的两个模式 ∣ψi⟩|\psi_i\rangle 和 ∣ψf⟩|\psi_f\rangle 之间的耦合矩阵元：

Mfi=⟨ψf∣O^∣ψi⟩.M_{fi} = \langle\psi_f|\hat{O}|\psi_i\rangle.

该矩阵元非零的必要条件是三者的直积包含全对称表示 A1A_1：

Γi⊗ΓO⊗Γf⊃A1,\Gamma_i \otimes \Gamma_O \otimes \Gamma_f \supset A_1,(28.73)

其中 Γ\Gamma 表示所属的 irrep。利用 C4vC_{4v} 的直积表可以快速判断：

tabularc||cccc|c ⊗\otimes & A1A_1 & A2A_2 & B1B_1 & B2B_2 & EE

A1A_1 & A1A_1 & A2A_2 & B1B_1 & B2B_2 & EE

A2A_2 & & A1A_1 & B2B_2 & B1B_1 & EE

B1B_1 & & & A1A_1 & A2A_2 & EE

B2B_2 & & & & A1A_1 & EE

EE & & & & & A1⊕A2⊕B1⊕B2A_1\oplus A_2\oplus B_1\oplus B_2

tabular

实例分析 1：四极矩辐射的选择规则

考虑四极矩辐射算符的分量 Qxy∝xyQ_{xy} \propto xy。从特征标表可知 xyxy 属于 B2B_2 表示。现在判断几组可能的跃迁： itemize A1↔B2A_1 \leftrightarrow B_2 跃迁：A1⊗B2⊗B2=A1⊗A1=A1A_1 \otimes B_2 \otimes B_2 = A_1 \otimes A_1 = A_1 ✓ （对称性允许） A1↔A1A_1 \leftrightarrow A_1 跃迁：A1⊗B2⊗A1=B2⊅A1A_1 \otimes B_2 \otimes A_1 = B_2 \not\supset A_1 ✗ （对称性禁戒） B1↔EB_1 \leftrightarrow E 跃迁：B1⊗B2⊗E=A2⊗E=E⊅A1B_1 \otimes B_2 \otimes E = A_2 \otimes E = E \not\supset A_1 ✗ （禁戒） itemize 这解释了为什么某些对称性保护的 BIC 无法向自由空间辐射——是因为对称性从根本上禁止了这种耦合，而非因为"损耗太小"。

实例分析 2：法向辐射通道的 irrep

在 Γ\Gamma 点，法向辐射对应的是横向均匀平面波，其两支线偏振 (x^,y^)(\hat{x},\hat{y}) 构成 EE 表示的基。因此： itemize 若腔模属于 EE 表示，则 Γmode⊗Γrad=E⊗E=A1⊕A2⊕B1⊕B2⊃A1\Gamma_{\rm mode} \otimes \Gamma_{\rm rad} = E \otimes E = A_1 \oplus A_2 \oplus B_1 \oplus B_2 \supset A_1 ✓ —— 允许辐射（亮模） 若腔模属于 A1A_1 表示，则 A1⊗E=E⊅A1A_1 \otimes E = E \not\supset A_1 ✗ —— 禁止法向辐射（BIC 候选） itemize 因此，高阶带边模式（通常是 EE 表示）可表面发射，而基模（通常是 A1A_1）可能成为暗模。

对称性破缺与不可约表示分解：当结构发生某种变形或引入各向异性时，对称性会降低。设原始对称群为 GG，破缺后的子群为 H⊂GH \subset G。此时 GG 的每个 irrep 按照一定规则"分解"为 HH 的若干个 irrep 的直和。这个过程称为分支律(branching rule)。

对于 PCSEL 常见的 C4v→C2vC_{4v} \to C_{2v} 破缺（例如引入椭圆形孔或矩形晶格），分解规则如下：

A1(C4v)→A1(C2v),A2(C4v)→A2(C2v),B1(C4v)→A1(C2v),B2(C4v)→A2(C2v),E(C4v)→B1(C2v)⊕B2(C2v).\begin{aligned}A_1(C_{4v}) &\to A_1(C_{2v}), \nonumber\\ A_2(C_{4v}) &\to A_2(C_{2v}), \nonumber\\ B_1(C_{4v}) &\to A_1(C_{2v}), \nonumber\\ B_2(C_{4v}) &\to A_2(C_{2v}), \nonumber\\ E(C_{4v}) &\to B_1(C_{2v}) \oplus B_2(C_{2v}). \end{aligned}(28.74)

关键的物理后果： itemize 双重态分裂: EE 表示分解为两个不同的单重态 B1B_1 和 B2B_2。这意味着原本简并的一对模式会发生能量劈裂，劈裂大小正比于对称性破缺强度。

偏振锁定: 分裂后的两个模式分别继承原 EE 表示的两个正交基，对应固定的线偏振方向。这在实验中表现为输出光的偏振方向被"钉扎"在某个晶轴方向。

选择规则放宽: 由于 C2vC_{2v} 的特征标表更简单（只有 4 个一维 irrep），原来在 C4vC_{4v} 中被禁戒的某些跃迁可能在 C2vC_{2v} 中变得允许。这可以用来解释为何轻微的结构不对称会导致 BIC 泄漏出有限的辐射。 itemize

实际应用建议

在设计 PCSEL 时，可以主动利用对称性工程： enumerate 若要实现稳定的单线偏振输出，故意打破 C4vC_{4v} 对称性（如采用椭圆孔或矩形占空比） 若要追求最高品质因子的 BIC，保持完美的 C4vC_{4v} 并确保工作点是 A1/A2/B1/B2A_1/A_2/B_1/B_2 之一 若要实现双波长或双偏振激光器，可以利用接近简并但略微分裂的准-EE 态 enumerate

关键量的首次定义位置索引

为便于读者在正文各章之间快速定位，本节给出全书最重要的派生量与合成量在何处首次引入或定义的参考索引。注意：基础物理量（如 E\E、ω\omega、k\kvec 等）不列，只列出由多个基础量组合而成、或由特定物理机制引入的派生量。

符号/量首次定义章节及位置
Γopt\Gamma_{\mathrm{opt}}第 7 章，式 (7.5)：光学约束因子，定义为模式与有源区的能量重叠。
αrad,  αint\alpha_{\mathrm{rad}},\;\alpha_{\mathrm{int}}第 7 章，式 (7.7)：辐射损耗与内部损耗的等效空间系数。
gthg_{\mathrm{th}}第 7 章，式 (7.12)：阈值增益，αrad+αint\alpha_{\mathrm{rad}}+\alpha_{\mathrm{int}} 与 Γopt\Gamma_{\mathrm{opt}} 之比。
κ1,  κ2\kappa_1,\;\kappa_2第 10 章，式 (10.10)：CWT 中的耦合系数，来自主导几何 Fourier 分量与纵向模重叠。
C1D\bm{C}_{\mathrm{1D}}第 11 章，式 (11.37)：一维反馈耦合矩阵。
Crad\bm{C}_{\mathrm{rad}}第 11 章，式 (11.39)：辐射耦合矩阵，非厄米，来自开放通道。
C2D\bm{C}_{\mathrm{2D}}第 11 章，式 (11.41)：二维高阶耦合矩阵。
venv_{\mathrm{en}}第 7 章，式 (7.10)：能量速度，将空间增益/损耗映射为时间增长率的参数。
n~2=n2+2ik0nα~\tilde n^2=n^2+2\ii k_0 n\tilde\alpha第 11 章，式 (11.5)：含增益/损耗的复折射率。
ξm,n(z)\xi_{m,n}(z)第 11 章，式 (11.7)：光子晶体层的 Fourier 系数。
ϕ0(z)\phi_0(z)第 11 章，式 (11.15)：背景层状波导的基模纵向剖面。
uA1,  uB1,  uE,x/y\bm{u}_{A_1},\;\bm{u}_{B_1},\;\bm{u}_{E,x/y}第 10 章，式 (10.17)、式 (10.21)、式 (10.28)：四波子空间在 C4vC_{4v} 下的 irrep 基底。
G(z,z′)G(z,z')第 11 章，式 (11.22)：辐射波的纵向 Green 函数（近似局域均匀背景）。
Gm,n(z,z′)G_{m,n}(z,z')第 11 章，式 (11.34)：高阶波在光子晶体层内的倏逝 Green 函数。
δ=β−β0\delta=\beta-\beta_0第 11 章，式 (11.17)：偏离 Bragg 条件的传播常数偏移。
αrad,j\alpha_{\mathrm{rad},j}第 11 章，式 (11.50)：第 jj 模的辐射常数，  2 Im(λj)\;2\,\mathrm{Im}(\lambda_j)（亦可写为 −2 Im(ω~j)/vg,j-2\,\mathrm{Im}(\tilde\omega_j)/v_{g,j}）。
ηNz\eta_{N_z}第 11 章，式 (11.4)：纵向模截断能量占比。
Σrad(ω)\mathbf{\Sigma}_{\mathrm{rad}}(\omega)第 11 章，式 (11.53)：辐射通道引入的自能项。
ΓPC\Gamma_{\mathrm{PC}}第 13 章，式 (13.26)：光子晶体层的光学限制因子，用于二维等效介电常数修正。
ε~a,  ε~b\tilde\varepsilon_a,\;\tilde\varepsilon_b第 13 章，式 (13.27)：送入二维 PWEM 的修正等效介电常数。
αrad=2 Im(λ)=−2 Im(ω~)/vg\alpha_{\mathrm{rad}}=2\,\mathrm{Im}(\lambda)=-2\,\mathrm{Im}(\tilde\omega)/v_g第 11 章，式 (11.50)：辐射常数与 CWT 本征值虚部（及复频率虚部）之间的换算关系。
IthI_{\mathrm{th}}第 7 章，式 (7.15)：阈值电流，是器件级电-热-光结果，需经电注入模型换算。
ΔG⋆j\Delta G_{\star j}第 7 章：目标模与第 jj 竞争模之间的模鉴别裕量。
∂ωm/∂T,  ∂γm/∂T\partial\omega_m/\partial T,\;\partial\gamma_m/\partial T第 11 章，式 (11.61)：三维 CWT 输出的温度灵敏度导数，用于热耦合接口。

这条索引的设计目的只有一个：让读者在遇到某个量时，能在 30 秒内找到它在主线上的定义位置，而不是在各章符号之间来回猜测。本索引不追求穷举，只覆盖"跨越两个以上章节使用"的派生量。基础符号的定义应在各章首次使用时已经明确，此处不再重复。

回看本章

数学附录的作用是统一全书的语法，而非堆公式。时间因子决定损耗与增益的符号，傅里叶卷积解释了周期介质中的耦合矩阵，QNM / resonant states 说明开放模为何不是普通束缚态，左/右模与广义归一化说明为什么开放系统不能直接照搬厄米内积，复频率给出 QQ 与光子寿命的精确定义，互易性与一阶微扰则为耦合和灵敏度分析提供了共同骨架。只要这些约定统一，正文中的光学、器件和仿真章节就能共享同一套数学接口；若这些约定混乱，后面的一切“高级模型”都可能只是表面连贯。

练习题

enumerate [基础] 从式 (28.1) 出发，推导在本书时间因子约定下被动介质为何对应正的 ℑ ε\Im\,\varepsilon。

[基础] 从能量衰减定义出发，重新推导式 (28.55)，并说明为什么振幅寿命和能量寿命差一个因子 2。

[进阶] 利用式 (28.13)，解释为什么保留有限个倒格矢后，本征问题会变成矩阵本征值问题。

[进阶] 从式 (28.64) 出发，解释为什么开放系统一阶微扰需要左模，而不能只用右模与自身的复共轭关闭投影。

[综合] 给出一个一阶微扰公式明显失效的 PCSEL 场景，并说明失效来自哪一条前提被破坏。 enumerate

延伸阅读

itemize 下一章将把这些数学语法压缩成器件级建模的符号链，特别是从材料增益到阈值电流的映射。 阅读时应同时核对本章变量在后续模型中的定义和适用范围。

若想进一步打牢开放系统与扰动理论，可对照 [Joannopoulos et al., 2008] [Sakoda, 2005] 中关于周期介质与模式扰动的附录阅读。 阅读时应重点核对开放系统归一化、QNM 与非厄米微扰的标准假设。

若想系统补强开放系统归一化、QNM 与非厄米微扰，可继续阅读 [Leung et al., 1994] [Leung et al., 1994] [Muljarov et al., 2010] [Sauvan et al., 2013] [Kristensen et al., 2015] [Alpeggiani et al., 2017]。 阅读时应重点核对开放系统归一化、QNM 与非厄米微扰的标准假设。 itemize

本附录数学工具在正文中的应用位置
数学工具正文应用具体用途
时间因子约定第 3,7,17 章统一损耗/增益符号
Bloch-Floquet 展开第 4,5,12 章PWEM/RCWA 的平面波基底
QNM 双正交归一化第 11 章非厄米微扰公式
互易定理第 8,13 章远场 - 近场关联
留数定理第 13 章RCWA 散射矩阵极点分析
卷积矩阵构造第 12,13 章周期介质的傅里叶因子化
复频率-QQ关系第 3,7,8b 章BIC/quasi-BIC的QQ值定义

PCSEL KNOWLEDGE BASE

PCSEL 理论、仿真与器件物理导论

正在载入可检索的全书知识站。JavaScript 不可用时，可直接阅读各章独立页面或下载原书 PDF。

从第一章开始 · 下载原书 PDF