总览与电磁基础

周期介质中的 Maxwell 方程

想象一只“关不上门”的腔：顶部完全封闭，光出不来；顶部一旦打开，能量又必然泄漏。PCSEL 必须在这两者之间工作，所以它的模式从一开始就不是普通封闭腔里的实频本征态，而是带辐射寿命的开放共振。本章要建立的 Maxwell 语言，正是为了把这扇门怎样打开、打开后频率和损耗怎样一起出现说清楚。

读完本章，你应能

读前准备

大学电动力学中的 Maxwell 方程、边界条件与向量恒等式。

阅读主线

本章所处模型层级：L1

本章先从 Maxwell 方程写起，再转入频域形式，随后讨论在何种条件下它成为本征值问题。接着我们解释为什么 PCSEL 从根本上是开放系统，并补上开放问题关键的数学直觉：共振模为什么不再属于普通平方可积函数空间、QNM / resonant states 应如何作最小理解。最后再把这一章和后面的 Bloch 理论、阈值分析接起来。

从时域方程出发：我们到底在求解什么

本章先明确所求解的电磁问题：PCSEL 中的光场由什么方程决定？这些方程中的哪些部分来自材料本身，哪些部分来自边界，哪些部分来自我们选择的近似？

在宏观电磁学层面，时域 Maxwell 方程写为

∇⋅D=ρf,∇⋅B=0,∇×E=−∂B∂t,∇×H=Jf+∂D∂t,\begin{aligned}\nabla\cdot\D &= \rho_f, \\ \nabla\cdot\B &= 0, \\ \nabla\times\E &= -\frac{\partial \B}{\partial t}, \\ \nabla\times\Hf &= \J_f + \frac{\partial \D}{\partial t},\end{aligned}

其中 ρf\rho_f 与 Jf\J_f 是自由电荷和自由电流。在本书讨论的腔模、本征模和被动共振问题中，我们通常取无外加源情形，即 ρf=0\rho_f=0、Jf=0\J_f=0。这并非说器件里没有载流子，而是说在电磁子问题中，载流子的作用被压缩进本构关系，例如介电函数、增益和损耗之中。

为了把 Maxwell 方程变成可计算问题，还需要本构关系。最基本的写法是

D=ε0εr(r,ω)E,B=μ0μr(r,ω)H.\D=\varepsilon_0\varepsilon_r(\rvec,\omega)\E,\qquad \B=\mu_0\mu_r(\rvec,\omega)\Hf.

在大多数半导体光学问题中，磁响应不显著，后文通常取 μr≈1\mu_r\approx 1。这一近似在近红外到可见光频段（约 400 nm--2 μ\mum）的非磁性半导体（GaAs、InP、GaN 等）中成立，因为轨道磁响应极弱，对应磁化率通常满足 χm=μr−1∼10−5\chi_m=\mu_r-1\sim 10^{-5}。对 THz 频段或人工磁响应超材料，此近似不成立；本书讨论的材料体系均满足 μr≈1\mu_r\approx 1。主要复杂性集中在 εr\varepsilon_r：它既包含空间分布，又可能包含频率色散、材料损耗、载流子注入和温度依赖。也就是说，PCSEL 的全部复杂性，最终都会不同方式写进 εr(r,ω)\varepsilon_r(\rvec,\omega) 之中。

转到频域：为什么教材几乎总在频域讨论模式

若系统在单频附近工作，或者我们要研究本征模与共振，最自然的做法是采用 第 2 章 中约定的时间因子 e−iωt\ee^{-\ii\omega t}，把场写成

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].

代回 Maxwell 方程时，关键步骤如下。以 Faraday 定律为例：

∇×E(r,t)=−∂B∂t=−∂∂t[μ0μrH(r)e−iωt]=iωμ0μrH(r)e−iωt.\nabla\times \E(\rvec,t)=-\frac{\partial \B}{\partial t} =-\frac{\partial}{\partial t}\bigl[\mu_0\mu_r\Hf(\rvec)\ee^{-\ii\omega t}\bigr] =\ii\omega\mu_0\mu_r\Hf(\rvec)\ee^{-\ii\omega t}.

消去公共因子 e−iωt\ee^{-\ii\omega t} 后即得频域 Faraday 方程。注意正号的来源：时间导数给出 −(−iω)=+iω-(-\ii\omega)=+\ii\omega，两个负号相消。若改用 e+iωt\ee^{+\ii\omega t} 约定，则 ∂t→+iω\partial_t\to +\ii\omega，所有含 iω\ii\omega 的符号均翻转——不同文献中频域 Maxwell 方程正负号不同，即源于时间因子约定不同。

对 Amp\`ere 定律也应把符号完整走一遍，而不只是“照抄结果”。从

∇×H(r,t)=∂D∂t=∂∂t[ε0εrE(r)e−iωt]=−iωε0εrE(r)e−iωt\nabla\times \Hf(\rvec,t)=\frac{\partial \D}{\partial t} =\frac{\partial}{\partial t}\bigl[\varepsilon_0\varepsilon_r\E(\rvec)\ee^{-\ii\omega t}\bigr] =-\ii\omega\varepsilon_0\varepsilon_r\E(\rvec)\ee^{-\ii\omega t}

可见，在 e−iωt\ee^{-\ii\omega t} 约定下，Amp\`ere 定律右端会直接给出负号。消去公共因子 e−iωt\ee^{-\ii\omega t} 后，频域方程写成

∇×E=iωμ0μrH,∇×H=−iωε0εrE.\begin{aligned}\nabla\times \E &= \ii\omega\mu_0\mu_r \Hf, \\ \nabla\times \Hf &= -\ii\omega\varepsilon_0\varepsilon_r \E. \end{aligned}(3.5)、(3.6)

对 式 (3.5) 再取一次旋度：

∇×(∇×E)=∇×(iωμ0μrH)=iωμ0μr∇×H.\nabla\times(\nabla\times \E)=\nabla\times(\ii\omega\mu_0\mu_r \Hf) =\ii\omega\mu_0\mu_r\nabla\times\Hf.

若 μr\mu_r 空间变化，则需写成 ∇×[μr−1(∇×E)]=iωμ0(∇×H)\nabla\times[\mu_r^{-1}(\nabla\times\E)]=\ii\omega\mu_0(\nabla\times\Hf)。再用 式 (3.6) 消去 ∇×H\nabla\times\Hf，得到电场主方程

∇×[μr−1(r)∇×E]−k02εr(r)E=0,k0=ωc.\nabla\times\left[\mu_r^{-1}(\rvec)\nabla\times\E\right]-k_0^2\varepsilon_r(\rvec)\E=0, \qquad k_0=\frac{\omega}{c}.(3.7)

这是后文常用的频域 Maxwell 方程。形式上它已经非常像一个本征值问题，但若要把它当作本征值问题，还需要说明两个问题：边界条件是什么？介质是否损耗、是否开放？

坡印廷定理与能量流向

回答这两个问题的最直接工具是能量流。利用矢量恒等式 ∇⋅(E×H∗)=H∗⋅(∇×E)−E⋅(∇×H∗)\nabla\cdot(\E\times\Hf^*) = \Hf^*\cdot(\nabla\times\E) - \E\cdot(\nabla\times\Hf^*)，并将 式 (3.5)、式 (3.6) 代入，经过直接计算可得

∇⋅(12E×H∗)=iω2(μ0μr∣H∣2−ε0εr∗∣E∣2).\nabla\cdot\left(\frac{1}{2}\E\times\Hf^*\right) = \frac{\ii\omega}{2}\left(\mu_0\mu_r|\Hf|^2 - \varepsilon_0\varepsilon_r^*|\E|^2\right).(3.8)

式 (3.8) 是无外加源、实频率下的关系。对被动介质，若定义

Pabs=∫Vω2ε0Im⁡(εr)∣E∣2 dV(Pabs≥0),P_{\mathrm{abs}} = \int_V\frac{\omega}{2}\varepsilon_0\operatorname{Im}(\varepsilon_r)|\E|^2\,\dd V \quad (P_{\mathrm{abs}}\ge 0),(3.9)

那么直接取实部会得到 Pout=−PabsP_{\mathrm{out}}=-P_{\mathrm{abs}}。这个负号不是说“被动腔一边吸收、一边凭空向外辐射”，而是在提醒我们：无源、被动、实频率的开放边界问题没有非平凡稳态本征解。实频散射问题必须加入外部电流或入射端口，此时功率守恒写成

Psrc=Pout+Pabs.P_{\mathrm{src}}=P_{\mathrm{out}}+P_{\mathrm{abs}}.(3.10)

无源开放腔则应使用复频率，回到时域后写成

dUdt=−Pout−Pabs,\frac{\dd U}{\dd t}=-P_{\mathrm{out}}-P_{\mathrm{abs}},(3.11)

其中 UU 是腔内储能。若边界为理想导体且介质无损，则两项均为零，允许实频率封闭腔本征模；PCSEL 允许光穿过开放边界，因而无源本征问题通常必须用复频率描述。判断“封闭”或“开放”要看边界上的能量通量和本征问题的频率条件，不能只看几何形状。

算符观点：为什么有时本征频率是实数，有时必须是复数

为了建立后文需要的概念，我们把 式 (3.7) 记成算符形式

L^(ω)E=0,L^(ω)=∇×μr−1∇×−k02εr(r,ω).\hat{\mathcal L}(\omega)\E = 0, \qquad \hat{\mathcal L}(\omega)=\nabla\times\mu_r^{-1}\nabla\times-k_0^2\varepsilon_r(\rvec,\omega).(3.12)

物理含义：方程式 (3.12)定义了开放光学系统的本征问题框架。算符L^(ω)\hat{\mathcal L}(\omega)的第一项∇×μr−1∇×\nabla\times\mu_r^{-1}\nabla\times描述磁场的旋度能量，第二项k02εrk_0^2\varepsilon_r描述电场在介质中的存储能力。求解det⁡[L^(ω)]=0\det[\hat{\mathcal L}(\omega)]=0即得到系统的本征频率和本征模式。

如果材料无损、无色散，边界又是理想封闭的，那么这个算符可以在适当函数空间上写成自伴形式。此时本征频率是实数，不同模式之间具有良好的正交关系。这是很多电磁教材首先讲授的情形，因为它数学上干净，也有助于建立能量守恒图像。最小类比是一自由度 RLC 共振器：理想无损时本征频率为实数，串入电阻后频率变成“振荡频率 + 衰减率”的复数——厄米对应“无电阻”，非厄米对应“有电阻”。

但 PCSEL 不是这样的系统。它要向自由空间辐射输出；只要允许辐射，系统就不再是封闭腔，模式就不再是严格束缚态，而会变成准束缚态（quasi-bound state）或共振态。在后面的 第 9 章 中（完整内容见第 8b 章，位于 Ch08 与 Ch09 之间）我们会进一步看到，某些 Γ\Gamma 点模甚至会以 BIC / quasi-BIC 的形式组织辐射通道。此时本征频率通常写成

ω~=ω0−iγ,\tilde{\omega}=\omega_0-\ii\gamma,

其中 γ\gamma 不再是“数值误差”，而是模式能量逸出或耗散的速率编码。因此，PCSEL 的被动模分析从一开始就必须带着“非厄米开放问题”的意识进行。图 3.1 展示了这些模式在复频率平面上的位置和物理含义。 ComparisonBox封闭谐振腔与开放 PCSEL 系统的对比

tabularl l l

性质 & 封闭腔 (金属壁) & 开放腔 (PCSEL)

本征频率 & 实数 ωn\omega_n & 复数 ωn−iγn\omega_n-i\gamma_n

模式正交性 & 标准内积正交 & 需要双正交基

边界条件 & n×E=0\bm{n}\times\bm{E}=0 & Sommerfeld 辐射条件

QQ 因子来源 & 导体损耗/介质损耗 & 辐射损耗 + 介质损耗

完备性 & 离散谱完备 & 离散谱 + 连续谱

tabular

ComparisonBox

复频率平面示意。被动开放系统的共振模位于下半平面（γ>0），距实轴越近表示 Q 越高。连续谱束缚态理想地处于实轴上（γ=0），准连续谱束缚态则刚刚离开实轴。
复频率平面示意。被动开放系统的共振模位于下半平面（γ>0\gamma>0），距实轴越近表示 QQ 越高。连续谱束缚态理想地处于实轴上（γ=0\gamma=0），准连续谱束缚态则刚刚离开实轴。图 3.1 · 原始 LaTeX/TikZ 图源独立渲染

对 PCSEL 常见近红外工作点（ω0∼1015 rad/s\omega_0\sim10^{15}\ \mathrm{rad/s}），由

Q=ω02γ⇒γω0=12QQ=\frac{\omega_0}{2\gamma}\quad\Rightarrow\quad \frac{\gamma}{\omega_0}=\frac{1}{2Q}

可得：Q=103Q=10^3 时 γ/ω0≈5×10−4\gamma/\omega_0\approx5\times10^{-4}，Q=105Q=10^5 时 γ/ω0≈5×10−6\gamma/\omega_0\approx5\times10^{-6}。因此“离实轴多远算高品质因子”可以直接用这个比例估算，而不必只靠定性描述。

若要把“无损、无色散、封闭边界下 Maxwell 问题可写成自伴形式”说得更正式，较自然的写法不是直接把 式 (3.7) 当普通矩阵本征值问题，而是先把它限制在

X={E∈H(curl;V) | ∇⋅ ⁣(εE)=0, 满足封闭边界条件}\mathcal X= \left\{ \E\in H(\mathrm{curl};V) \ \middle|\ \nabla\cdot\!\bigl(\varepsilon\E\bigr)=0,\ \text{满足封闭边界条件} \right\}

这一约束子空间上，再把它改写成广义本征值问题或等价的加权自伴算子问题。下面的命题与证明只给出教学化 proof sketch：它足以解释物理上为什么得到实频率，但不是函数分析层面的完整证明。

更稳妥的写法是把 Maxwell 方程视为广义本征值问题

AE=λBE,A=∇×μr−1∇×,B=εr,λ=(ωc)2.A\E=\lambda B\E, \qquad A=\nabla\times\mu_r^{-1}\nabla\times, \quad B=\varepsilon_r, \quad \lambda=\left(\frac{\omega}{c}\right)^2.(3.13)

这里 AA 与 BB 的 Hermitian 性先在普通 L2L^2 内积 ⟨E1,E2⟩=∫VE1∗⋅E2 dV\langle\E_1,\E_2\rangle=\int_V\E_1^*\cdot\E_2\,\dd V 下判断。利用向量恒等式对旋度做分部积分，得到

∫VE1∗⋅(∇×μr−1∇×E2) dV=∫V(∇×E1)∗⋅μr−1(∇×E2) dV−∮∂V[E1∗×μr−1(∇×E2)]⋅n^ dS.\int_V \E_1^\ast\cdot\bigl(\nabla\times\mu_r^{-1}\nabla\times \E_2\bigr)\,\dd V = \int_V (\nabla\times\E_1)^\ast\cdot\mu_r^{-1}(\nabla\times \E_2)\,\dd V -\oint_{\partial V}\bigl[\E_1^\ast\times\mu_r^{-1}(\nabla\times\E_2)\bigr]\cdot\hat{\bm n}\,\dd S.

若边界为理想导体（PEC），则在边界上 n^×E=0\hat{\bm n}\times\E=\bm{0}，表面积分项为零；若为周期边界，相对两侧的法向相反且场值满足 Bloch 条件，表面项成对抵消。更具体地，在两侧边界 r\rvec 与 r+R\rvec+\bm{R} 上有

E(r+R)=E(r)eik⋅R,\E(\rvec+\bm{R})=\E(\rvec)\ee^{\ii\kvec\cdot\bm{R}},

表面积分的相位因子与法向翻转带来的负号相抵，因而相消。这一步只在 k\kvec 为实向量时成立；若 k\kvec 取复值（表示面内衰减），表面项一般不再抵消。对无损介质，μr\mu_r 与 εr\varepsilon_r 为正实 Hermitian 张量，因此 AA 与 BB 都是 Hermitian，且 BB 正定。于是

λ=⟨E,AE⟩⟨E,BE⟩∈R≥0,\lambda = \frac{\langle\E,A\E\rangle}{\langle\E,B\E\rangle} \in\mathbb{R}_{\ge0},(3.14)

从而 ω\omega 可取实数。等价地，也可以用 B−1/2AB−1/2B^{-1/2}AB^{-1/2} 把它变成普通 Hermitian 本征问题。需要注意：这里证明的是广义本征对 (A,B)(A,B) 的 Hermitian 结构，并没有证明旋度算符在 ε\varepsilon 加权内积下直接自伴。

但对开放辐射边界，∂V\partial V 延伸到无穷远，出射波携带能量使表面积分不为零且不对称，算符不再自伴。本征频率必须取复数 ω~=ω0−iγ\tilde{\omega}=\omega_0-\ii\gamma，其中 γ>0\gamma>0 是场振幅衰减率，储能衰减率为 2γ2\gamma。

开放 Maxwell 问题的函数空间直觉：为什么它不再是普通闭腔本征值问题

这里需要把一个常被轻描淡写、但对 PCSEL 极关键的数学事实讲清楚。封闭腔本征模通常属于普通的平方可积函数空间，因为场在腔外被边界截断或指数衰减，能量积分是有限的。开放系统中的共振模则不同。若某个模满足出射波条件，其远场在三维直觉上可写成

E(r)∼eik~rrf(r^),k~=ω~c,\E(\rvec)\sim \frac{\ee^{\ii \tilde{k} r}}{r}\bm{f}(\hat{\rvec}), \qquad \tilde{k}=\frac{\tilde{\omega}}{c},(3.15)

其中 ω~=ωr−iγ\tilde{\omega}=\omega_r-\ii\gamma，且 γ>0\gamma>0。把指数项拆成振荡部分与增长部分后可得

eik~r=eiωrr/c e+γr/c.\ee^{\ii \tilde{k} r} = \ee^{\ii \omega_r r/c}\,\ee^{+\gamma r/c}.(3.16)

这说明一个反直觉但标准的事实：若把开放共振模按纯出射边界条件延拓到无穷远，它在空间上通常不是衰减而是增长的。 这并不表示它“物理上发散了”，而是说它对应的是有限时间窗口内的辐射拖尾（因果响应），并不要求在 t→+∞t\to+\infty 的真实场中保持该增长；在严格的散射/Green 张量解析延拓意义下，它是共振极点态，而不是普通 Hilbert 空间里的束缚态。

这一步必须讲明白，因为它直接解释了为什么开放系统中的归一化、正交性和微扰理论不能直接照搬封闭厄米问题：问题的根源不只是“损耗让本征值变成复数”，而是函数空间本身已经换了。

QNM / resonant states 的最小定义

在本书语境中，一个准正规模（quasinormal mode, QNM）或共振态（resonant state），最小可以定义为：在无外加源条件下满足 Maxwell 方程，并同时满足纯出射边界条件的非平凡解，

L^(ω~n)EnR=0,仅允许 outgoing-wave solution.\hat{\mathcal L}(\tilde{\omega}_n)\E_n^{\mathrm{R}}=0, \qquad \text{仅允许 outgoing-wave solution}.(3.17)

一个常用的等价描述是 Sommerfeld 辐射条件。以一维标量场为例：

lim⁡r→∞r(∂E∂r−ik0E)=0,\lim_{r\to\infty} r\left(\frac{\partial E}{\partial r}-\ii k_0 E\right)=0,(3.18)

表示无穷远处只有向外传播分量，没有向内入射分量。这与束缚态要求的 E→0E\to 0（r→∞r\to\infty）完全不同；正因为 ω~\tilde{\omega} 的虚部存在， QNM 在空间上可指数增长，导致其 L2L^2 范数发散。数值计算中采用 PML， 正是将这一出射条件以可计算的方式实现。这里上标 R\mathrm{R} 表示右模（right mode）， 是为了和后面会出现的伴随模区分。对应的复频率 ω~n\tilde{\omega}_n 可等价地理解为 Green 张量或散射矩阵的极点位置；也就是说，开放系统的模 不是先于散射存在的独立对象，而正是散射响应在复频率平面上的结构化奇点 [Fan & Joannopoulos, 2002] [Sauvan et al., 2013] [Alpeggiani et al., 2017]。

这个定义已经足够支撑本书后面的物理主线。它告诉我们：

围绕这一极点态身份重新建立的归一化、伴随模与双正交框架，见 第 28 章 中“伴随模、双正交与广义归一化”一节，尤其可对照 式 (28.48) 及其前后说明。

为什么 PCSEL 作为开放系统

到这里，我们可以回答一个看似简单、实际上极其关键的问题：为什么本书反复强调 PCSEL 是开放系统？

原因很直接。PCSEL 的目标之一就是法向出光。既然要出光，就必须允许电磁能量离开器件进入自由空间。能量能够离开，模式就不可能是理想封闭腔中的永不衰减本征态。于是：

因此，即使只进行能带分析，也需要保留开放性的影响范围。一个在封闭近似下看起来很好的模式，到了开放系统中可能因为辐射耦合太强而损耗很大；反过来，一个能有效法向出射的模式也一定不是完全束缚的。PCSEL 的设计任务始终围绕“控制目标模式的辐射通道和泄漏速率”，而不追求“完全抑制泄漏”。

复频率、品质因子与光子寿命

设某个被动共振模的复频率为

ω~=ω0−iγ.\tilde{\omega}=\omega_0-\ii\gamma.

若把时间依赖写成 ℜ ⁣[E(r)e−iω~t]\Re\!\left[\E(\rvec)\ee^{-\ii\tilde{\omega} t}\right]，则模场振幅满足

∣E(t)∣∝e−γt.|\E(t)|\propto \ee^{-\gamma t}.

于是，振幅本身的衰减时间尺度是 τamp=1/γ\tau_{\mathrm{amp}}=1/\gamma。但器件分析里更常用的是储能或光子数寿命。由于电磁能量与场强平方成正比，

U(t)∝∣E(t)∣2∝e−2γt,U(t)\propto |\E(t)|^2\propto \ee^{-2\gamma t},

因此损耗功率满足

Ploss=−dUdt=2γU.P_{\mathrm{loss}}=-\frac{\dd U}{\dd t}=2\gamma U.

这就解释了为什么品质因子的分母是 2γ2\gamma 而不是 γ\gamma：QQ 比较的是单位振荡周期内储能衰减得有多慢，所以它天然对应能量寿命，而不是振幅寿命。于是可定义光子寿命

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

进一步，品质因子定义为

Q=ω02γ=ω0τp.Q=\frac{\omega_0}{2\gamma}=\omega_0\tau_p.(3.20)

这两个式子在后文会不断出现，因为它们把“复频率”这种看似抽象的本征值，转换成了器件工程师更熟悉的时间尺度和损耗尺度。

但 QQ 只是被动冷腔信息，式 (3.20) 是后续阈值分析的输入之一，而不是终点2。

开放边界的实现：数学边界条件与数值边界条件不是一回事

在数学上，开放系统的理想边界即 式 (3.18) 的出射波条件（outgoing-wave condition）：远处只允许向外跑的波，不允许从无穷远无源返回的入射波。对矢量 Maxwell 问题，严格形式会更复杂，但物理含义完全相同。边界条件本身，就是开放系统定义的一部分。

数值计算无法把空间做到无穷大，于是需要用辐射边界条件、端口条件或完全匹配层（perfectly matched layer, PML）来近似它。

三类常用边界机制及其组合规则。(a) 出射辐射边界直接施加外向波条件；(b) PML 通过复坐标延拓衰减外向波；(c) Bloch 边界把无限周期结构约化为单胞。底部选择条强调：周期性与开放性属于两个不同方向上的问题，常需组合使用，并分别做相位与边界收敛检查。
三类常用边界机制及其组合规则。(a) 出射辐射边界直接施加外向波条件；(b) PML 通过复坐标延拓衰减外向波；(c) Bloch 边界把无限周期结构约化为单胞。底部选择条强调：周期性与开放性属于两个不同方向上的问题，常需组合使用，并分别做相位与边界收敛检查。图 3.2 · 原始 LaTeX/TikZ 图源独立渲染

图 3.2对比了几种常用的开放边界实现方式及其适用场景。

PML 的核心思想值得在这里做一个最小介绍。它在计算域边界处引入一层特殊介质，通过对坐标做复数延拓

x→x+iσ(x)ω,x \to x + \frac{\ii\sigma(x)}{\omega},

使出射波在该层中指数衰减，但由于变换是连续的、在物理--PML 界面上梯度为零，波在进入该层时不产生反射。需要注意的是，σ/ω\sigma/\omega 显式依赖频率；对QNM 这类复频率 ω~\tilde{\omega} 问题，取用哪一个 ω\omega 会影响虚部精度，因此常需做收敛性检查。

若想看到“为什么会衰减”的最小代数，可取一维右行波 E(x)=eikxE(x)=\ee^{\ii kx}。代入复坐标延拓后：

E→exp⁡ ⁣[ik ⁣(x+iσxω)]=eikx e−kσx/ω.E \to \exp\!\left[\ii k\!\left(x+\frac{\ii\sigma x}{\omega}\right)\right] = \ee^{\ii kx}\,\ee^{-k\sigma x/\omega}.

后面的指数因子 e−kσx/ω\ee^{-k\sigma x/\omega} 就是 PML 区域中的吸收包络；它直接展示了“复坐标 →\rightarrow 指数衰减”的物理效果。

一旦开放边界实现不当，就会发生两类严重误判：

CommonPitfallBox开放边界 vs 封闭边界 错误理解：将 PCSEL 当成封闭谐振腔，使用理想导体（PEC）边界条件。封闭边界会人为产生反射驻波，导致 QQ 值被严重高估（可能差 3--5 个数量级）、模式频率偏移、远场方向图完全错误。

判别方法：如果仿真得到 Q>106Q>10^6 但没有使用 BIC 机制，很可能是边界条件错误。反过来，“复频率的虚部很小，说明数值做得很准”这种说法也没有物理意义：虚部大小首先反映的是模式本身的损耗与开放性，只有在收敛和边界条件正确的前提下，才谈得上数值精度。 CommonPitfallBox

这一章和后文的关系：为什么先讲 Maxwell，再讲 Bloch

到这里，读者应该能看出：本章建立的是问题类型，而不是最终解法。我们已经知道 PCSEL 的被动光学子问题来自开放系统中的 Maxwell 方程，但还没有利用二维周期性。周期性一旦被引入，问题就可以借助 Bloch 定理进行分类和简化，这也是下一章的主题。

从逻辑上说，本章回答的是“PCSEL 的电磁问题是什么”；下一章回答的是“在周期性存在时，这个问题如何被组织和分类”；再下一章才会进一步问“在这些被组织好的模式中，哪些最有可能成为 PCSEL 关心的候选模”。

回看本章

本章把 PCSEL 的电磁基础明确成了一个开放、通常非厄米的 Maxwell 本征问题。频域 Maxwell 方程给出模式条件，复频率给出损耗时间尺度，开放边界决定模式是否能正确反映真实面发射结构，而 QNM / resonant states 则给出了开放共振模最小且严格的数学身份。后续所有能带、Γ\Gamma 点、耦合波和阈值分析，都建立在这一本章奠定的框架之上；其中更严格的归一化、伴随模与非厄米微扰公式将在 第 12 章、第 28 章 中继续展开。

练习题

  1. [基础] 从频域 Maxwell 方程出发，推导磁场形式的主方程，并说明它与电场形式在数值上可能各有什么优缺点。

  2. [基础] 解释为什么“面发射”这一目标本身就迫使 PCSEL 成为开放系统。

  3. [进阶] 在本书的时间因子约定下，证明 式 (3.20) 成立。

  4. [进阶] 利用式 (3.15)、式 (3.16)，解释为什么开放共振模通常不属于普通平方可积函数空间。

  5. [综合] 试举出一个边界条件实现错误会如何误导 PCSEL 设计判断的例子。

延伸阅读

PCSEL KNOWLEDGE BASE

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

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

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