总览与电磁基础

周期介质中的 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)

对全空间积分并取实部，得到时间平均功率流

Pout=∮∂V12Re⁡(E×H∗)⋅ds=−∫Vω2Im⁡(εr)ε0∣E∣2 dV.P_{\text{out}} = \oint_{\partial V}\frac{1}{2}\operatorname{Re}(\E\times\Hf^*)\cdot\dd\bm{s} = -\int_V\frac{\omega}{2}\operatorname{Im}(\varepsilon_r)\varepsilon_0|\E|^2\,\dd V.(3.9)

若边界为理想导体（n×E=0\bm{n}\times\E=0）且介质无损（Im⁡εr=0\operatorname{Im}\varepsilon_r=0），则两侧均为零，能量既不能流出也不能耗散，只能永远囚禁在腔内；而 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.10)

物理含义：方程式 (3.10)定义了开放光学系统的本征问题框架。算符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：它足以解释物理上为什么得到实频率，但不是函数分析层面的完整证明。

对两个试探场 E1,E2\E_1,\E_2，构造加权内积

⟨E1,E2⟩ε=∫Vεr(r)E1∗⋅E2 dV.\langle \E_1,\E_2\rangle_\varepsilon=\int_V \varepsilon_r(\rvec)\E_1^\ast\cdot\E_2\,\dd V.

在上述约束子空间内，下面只抓住对“为什么封闭无损时频率为实数”关键的一步：旋度项在该加权内积下可通过分部积分变成对称形式。利用向量恒等式 ∇⋅(A×B)=B⋅(∇×A)−A⋅(∇×B)\nabla\cdot(\bm{A}\times\bm{B})=\bm{B}\cdot(\nabla\times\bm{A})-\bm{A}\cdot(\nabla\times\bm{B}) 对旋度做分部积分，可以把作用在 E2\E_2 上的旋度“转移”到 E1\E_1 上：

∫VE1∗⋅(∇×∇×E2) dV=∫V(∇×E1)∗⋅(∇×E2) dV−∮∂V(E1∗×∇×E2)⋅n^ dS.\int_V \E_1^\ast\cdot(\nabla\times\nabla\times \E_2)\,\dd V = \int_V (\nabla\times\E_1)^\ast\cdot(\nabla\times \E_2)\,\dd V -\oint_{\partial V}(\E_1^\ast\times\nabla\times\E_2)\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 取复值（表示面内衰减），表面积分不再严格抵消，算符也不再严格自伴。对调 E1\E_1 和 E2\E_2 后结果不变，于是

⟨E1,L^E2⟩ε=⟨L^E1,E2⟩ε,\langle \E_1,\hat{\mathcal L}\E_2\rangle_\varepsilon=\langle \hat{\mathcal L}\E_1,\E_2\rangle_\varepsilon,

这也是自伴性的定义。自伴算符的本征值必为实数。

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

开放 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.11)

其中 ω~=ω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.12)

这说明一个反直觉但标准的事实：若把开放共振模按纯出射边界条件延拓到无穷远，它在空间上通常不是衰减而是增长的。 这并不表示它“物理上发散了”，而是说它对应的是有限时间窗口内的辐射拖尾（因果响应），并不要求在 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.13)

一个常用的等价描述是 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.14)

表示无穷远处只有向外传播分量，没有向内入射分量。这与束缚态要求的 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.15)

进一步，品质因子定义为

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

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

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

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

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

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

开放边界条件的数值实现方式对比。(a)辐射边界条件：在计算域边界施加出射波条件，要求波只能向外传播；(b)完全匹配层（PML）：通过在边界处引入复坐标延拓使外向波指数衰减，理想情况下无反射；(c) 周期性边界条件：利用Bloch 定理将无限周期结构约化到单个原胞；(d) 边界条件选择指南：根据问题类型（本征模、散射、周期结构）选择合适的边界条件组合。注意 PML并非物理实体，而是将无限空间折叠进有限计算域的数学技巧。
开放边界条件的数值实现方式对比。(a)辐射边界条件：在计算域边界施加出射波条件，要求波只能向外传播；(b)完全匹配层（PML）：通过在边界处引入复坐标延拓使外向波指数衰减，理想情况下无反射；(c) 周期性边界条件：利用Bloch 定理将无限周期结构约化到单个原胞；(d) 边界条件选择指南：根据问题类型（本征模、散射、周期结构）选择合适的边界条件组合。注意 PML并非物理实体，而是将无限空间折叠进有限计算域的数学技巧。图 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.16) 成立。

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

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

延伸阅读

PCSEL KNOWLEDGE BASE

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

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

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