总览与电磁基础

Bloch 定理与二维光子晶体基础

读完本章,你应能

读前准备

第 3 章第 2 章 中的频域 Maxwell 方程、坐标约定和倒格矢记号。

阅读主线

本章所处模型层级:L1

本章先从“周期性如何约束解”讲起,再给出 Bloch 定理的证明思路,随后介绍二维光子晶体几何和布里渊区,最后把这些理想周期概念放回具有有限厚度的 PCSEL slab 结构中,解释 TE/TM 为什么经常只是近似语言。

为什么周期性如此重要

上一章已经说明,PCSEL 的被动光学问题来自开放系统中的 Maxwell 方程。现在我们加入新的结构信息:介电函数在平面内是周期的。设存在一组晶格平移 R\bm{R},使得

ε(ρ+R,z)=ε(ρ,z).\varepsilon(\rho+\bm{R},z)=\varepsilon(\rho,z).(4.1)

这条周期性条件直接改变了解空间的结构。它意味着:虽然结构在平面上无限延展,但解并不需要在整个无限空间里“毫无组织地”寻找。周期性会把解空间分解成一族由波矢 k\kvec 标记的子空间。Bloch 定理正是这种分解的数学表达。

对初学者来说,Bloch 定理最容易被误记成一句形式化定义

Ek(r)=eikρuk(r),\E_{\kvec}(\rvec)=\ee^{\ii\kvec\cdot\rho}\bm{u}_{\kvec}(\rvec),

然后就往下继续算。这样当然可以做题,但难以说明这一形式为何与倒格矢、高对称点和能带图相联系。下面从平移算符说明其来源。

Bloch 定理的来源:平移算符与共同本征态

考虑晶格平移算符 T^R\hat T_{\bm{R}},定义为

T^RE(ρ,z)=E(ρ+R,z).\hat T_{\bm{R}}\E(\rho,z)=\E(\rho+\bm{R},z).

若介质满足 式 (4.1),则 Maxwell 算符在平面内与所有 T^R\hat T_{\bm{R}} 对易。这个对易关系最好显式写一次。设频域算子记为

L^=×μr1(ρ,z)×ω2c2εr(ρ,z),\hat{\mathcal L} = \nabla\times \mu_r^{-1}(\rho,z)\nabla\times -\frac{\omega^2}{c^2}\varepsilon_r(\rho,z),

则对任意场 E\E

T^RL^E(ρ,z)=L^(ρ+R,z)E(ρ+R,z).\hat T_{\bm R}\hat{\mathcal L}\E(\rho,z) = \hat{\mathcal L}(\rho+\bm R,z)\E(\rho+\bm R,z).

一旦结构满足

εr(ρ+R,z)=εr(ρ,z),μr(ρ+R,z)=μr(ρ,z),\varepsilon_r(\rho+\bm R,z)=\varepsilon_r(\rho,z), \qquad \mu_r(\rho+\bm R,z)=\mu_r(\rho,z),

就有 L^(ρ+R,z)=L^(ρ,z)\hat{\mathcal L}(\rho+\bm R,z)=\hat{\mathcal L}(\rho,z),从而

T^RL^E=L^T^RE,[T^R,L^]=0.\hat T_{\bm R}\hat{\mathcal L}\E = \hat{\mathcal L}\hat T_{\bm R}\E, \qquad\Longrightarrow\qquad [\hat T_{\bm R},\hat{\mathcal L}]=0.

对易意味着可以寻找同时是 Maxwell 算符和平移算符本征态的场。于是存在一组复数 λ(R)\lambda(\bm{R}) 使得

T^RE(ρ,z)=λ(R)E(ρ,z).\hat T_{\bm{R}}\E(\rho,z)=\lambda(\bm{R})\E(\rho,z).

现在利用群性质:

T^R1T^R2=T^R1+R2.\hat T_{\bm{R}_1}\hat T_{\bm{R}_2}=\hat T_{\bm{R}_1+\bm{R}_2}.

对同一本征态作用后得到

λ(R1)λ(R2)=λ(R1+R2).\lambda(\bm{R}_1)\lambda(\bm{R}_2)=\lambda(\bm{R}_1+\bm{R}_2).

满足这种可乘性质的连续表示只能写成

λ(R)=eikR,\lambda(\bm{R})=\ee^{\ii\kvec\cdot\bm{R}},

其中 k\kvec 是某个二维向量。一个简要论证如下:若 λ(R)\lambda(\bm{R}) 连续且满足 λ(R1)λ(R2)=λ(R1+R2)\lambda(\bm{R}_1)\lambda(\bm{R}_2)=\lambda(\bm{R}_1+\bm{R}_2),则可写 lnλ(R)=cR\ln\lambda(\bm{R})=\bm{c}\cdot\bm{R},其中 c\bm{c} 为常向量。无损周期系统要求 λ=1|\lambda|=1,因此 c\bm{c} 必为纯虚数,记作 ik\ii\kvec,于是得到指数形式。若允许有损或倏逝 Bloch 解,则 k\kvec 可取复数(“复 Bloch 波矢”),对应沿晶格方向的衰减或放大。

于是便有

Ek(ρ+R,z)=eikREk(ρ,z),\E_{\kvec}(\rho+\bm{R},z)=\ee^{\ii\kvec\cdot\bm{R}}\E_{\kvec}(\rho,z),(4.2)

再把快速相位因子提出,就得到Bloch 形式

Ek(ρ,z)=eikρuk(ρ,z),uk(ρ+R,z)=uk(ρ,z).\E_{\kvec}(\rho,z)=\ee^{\ii\kvec\cdot\rho}\bm{u}_{\kvec}(\rho,z), \qquad \bm{u}_{\kvec}(\rho+\bm{R},z)=\bm{u}_{\kvec}(\rho,z).(4.3)
一维 Bloch定理直观演示。(a)周期势场V(x)(黑色)与 Bloch波函数|_k(x)|^2(彩色);(b) Bloch波函数的分解:包络函数|u_k(x)|^2(红色虚线,具有与晶格相同的周期性)乘以平面波因子(蓝色振荡);(c) 电子系统与光子系统的类比:尽管控制方程不同(Schrödinger 方程vs Maxwell 方程),但由于两者都是线性波动方程且具有周期系数,Bloch 定理的数学形式完全相同。这使得我们可以将固体
一维 Bloch定理直观演示。(a)周期势场V(x)V(x)(黑色)与 Bloch波函数ψk(x)2|\psi_k(x)|^2(彩色);(b) Bloch波函数的分解:包络函数uk(x)2|u_k(x)|^2(红色虚线,具有与晶格相同的周期性)乘以平面波因子(蓝色振荡);(c) 电子系统与光子系统的类比:尽管控制方程不同(Schrödinger 方程vs Maxwell 方程),但由于两者都是线性波动方程且具有周期系数,Bloch 定理的数学形式完全相同。这使得我们可以将固体物理中的能带理论直接迁移到光子晶体。图 4.1 · 原始 LaTeX/TikZ 图源独立渲染

图 4.1通过一维示例直观展示了Bloch 定理的物理图像。

Bloch 定理的本质是利用平移对称性把无限系统的解空间按表示分类,而非经验性地假设“平面波乘周期函数”。

这里还隐含了一个对初学者很重要的前提:本章讨论的是无损、无限、严格周期的理想背景,因此我们把 k\kvec 视为实向量,并要求

λ(R)=1.|\lambda(\bm{R})|=1.

这个条件的物理意思是:把模式平移一个晶格周期后,只允许它获得相位,而不允许幅值沿晶格方向系统性增长或衰减。若问题改成有限周期链、有损介质或倏逝 Bloch 解,则复 Bloch 波矢同样可能出现;但那已经超出了本章“无限周期模式分类”的任务边界。

MathematicalFoundationBoxBloch 定理的严格数学证明 设周期势场满足 V(r+R)=V(r)V(\mathbf{r}+\mathbf{R})=V(\mathbf{r}),其中 R\mathbf{R} 为任意晶格矢量。定义平移算符 TRT_\mathbf{R} 使得 TRf(r)=f(r+R)T_\mathbf{R}f(\mathbf{r}) = f(\mathbf{r}+\mathbf{R})

第一步:证明 [H,TR]=0[H, T_\mathbf{R}] = 0。动能项 22m2-\frac{\hbar^2}{2m}\nabla^2 是空间导数算符,与平移操作对易,因此只需验证势能项:

TR[V(r)ψ(r)]=V(r+R)ψ(r+R)=V(r)TRψ(r),\begin{aligned}T_\mathbf{R}[V(\mathbf{r})\psi(\mathbf{r})] &= V(\mathbf{r}+\mathbf{R})\psi(\mathbf{r}+\mathbf{R}) \\ &= V(\mathbf{r})T_\mathbf{R}\psi(\mathbf{r}),\end{aligned}

HTR=TRHHT_\mathbf{R} = T_\mathbf{R}H

第二步:由量子力学基本原理,可对易的算符拥有共同本征态。设 ψ\psiHHTRT_\mathbf{R} 的共同本征态:

TRψ(r)=c(R)ψ(r).T_\mathbf{R}\psi(\mathbf{r}) = c(\mathbf{R})\psi(\mathbf{r}).

第三步:确定 c(R)c(\mathbf{R}) 的形式。由于平移不改变波函数归一化,要求 c(R)2=1|c(\mathbf{R})|^2=1,故可写 c(R)=eiθ(R)c(\mathbf{R})=e^{i\theta(\mathbf{R})}。又因平移操作的群性质 TR1TR2=TR1+R2T_{\mathbf{R}_1}T_{\mathbf{R}_2}=T_{\mathbf{R}_1+\mathbf{R}_2},可得:

c(R1)c(R2)=c(R1+R2).c(\mathbf{R}_1)c(\mathbf{R}_2) = c(\mathbf{R}_1+\mathbf{R}_2).

该函数方程的唯一连续解为 c(R)=eikRc(\mathbf{R}) = e^{i\mathbf{k}\cdot\mathbf{R}},其中 k\mathbf{k} 为与 R\mathbf{R} 无关的常矢量。

第四步:构造 Bloch 函数。令 ψ(r)=eikru(r)\psi(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}}u(\mathbf{r}),代入本征方程:

TRψ(r)=ψ(r+R)=eik(r+R)u(r+R).T_\mathbf{R}\psi(\mathbf{r}) = \psi(\mathbf{r}+\mathbf{R}) = e^{i\mathbf{k}\cdot(\mathbf{r}+\mathbf{R})}u(\mathbf{r}+\mathbf{R}).

另一方面,由本征值关系:

TRψ(r)=eikRψ(r)=eikReikru(r).T_\mathbf{R}\psi(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{R}}\psi(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{R}}e^{i\mathbf{k}\cdot\mathbf{r}}u(\mathbf{r}).

比较两式立得 u(r+R)=u(r)u(\mathbf{r}+\mathbf{R})=u(\mathbf{r}),即 uu 具有晶格周期性。证毕。

关键注记:上述证明中假设了 λ(R)=1|\lambda(\mathbf{R})|=1,这对应于无损耗的保守系统。若介质存在增益或损耗,或者考虑有限结构的倏逝模,则 Bloch 波矢可以是复数:k=k+ik\mathbf{k} = \mathbf{k}' + i\mathbf{k}''。此时 c(R)=ekR1|c(\mathbf{R})| = e^{-\mathbf{k}''\cdot\mathbf{R}} \neq 1,但仍满足乘法群同态关系。这也是激光腔中有源区建模时需要推广的理论框架。 MathematicalFoundationBox

Bloch 定理讨论的是无限周期结构中的模式分类。它告诉我们如何组织候选模,但并不自动告诉我们有限器件中哪个模会激射,也不直接给出阈值电流。

晶格平移与倒格矢记号

若二维晶格的原胞基矢为 a1,a2\bm{a}_1,\bm{a}_2,则倒格矢基矢 b1,b2\bm{b}_1,\bm{b}_2 定义为

aibj=2πδij.\bm{a}_i\cdot\bm{b}_j=2\pi\delta_{ij}.(4.4)

这里写 2π2\pi 而不是 11,是为了保证任何晶格平移对应的相位因子 exp(iGR)\exp(i\gvec\cdot\bm{R}) 对所有整数 mi,nim_i,n_i 都严格等于 1,从而 Bloch 条件与傅里叶展开可以直接用整数标记。部分固体物理教材采用 aibj=δij\bm{a}_i\cdot\bm{b}_j=\delta_{ij} 的约定,两种写法等价但符号不同,查阅文献时需注意对照。 任何晶格平移都可写为

R=m1a1+m2a2,m1,m2Z,\bm{R}=m_1\bm{a}_1+m_2\bm{a}_2,\qquad m_1,m_2\in\mathbb{Z},

任何倒格矢都可写为

G=n1b1+n2b2,n1,n2Z.\gvec=n_1\bm{b}_1+n_2\bm{b}_2,\qquad n_1,n_2\in\mathbb{Z}.

之所以在一开始就引入倒格矢,是因为 PCSEL 的许多耦合图像本质上都是“相差一个倒格矢的平面波分量之间发生散射”。如果读者不熟悉倒空间记号,后面讨论 Bloch 波、布里渊区、耦合波理论(CWT)和 RCWA 时会不断中断。

对于方格晶格,若 a1=ax^\bm{a}_1=a\hat{\bm{x}}a2=ay^\bm{a}_2=a\hat{\bm{y}},则由 式 (4.4) 立刻得到 b1=(2π/a)x^\bm{b}_1=(2\pi/a)\hat{\bm{x}}b2=(2π/a)y^\bm{b}_2=(2\pi/a)\hat{\bm{y}}。这个最简单的例子会在后文多次作为参考系出现。

为什么波矢只需取第一布里渊区

Bloch 波矢并非任意二维向量都互不相同。若 G\gvec 是倒格矢,则

ei(k+G)R  =  eikReiGR  =  eikR,\ee^{\ii(\kvec+\gvec)\cdot\bm{R}} \;=\; \ee^{\ii\kvec\cdot\bm{R}}\ee^{\ii\gvec\cdot\bm{R}} \;=\; \ee^{\ii\kvec\cdot\bm{R}},

因为 GR\gvec\cdot\bm{R} 总是 2π2\pi 的整数倍。于是 k\kveck+G\kvec+\gvec 对应同一组 Bloch 相位条件。因此,我们只需要在倒空间中选取一个基本区域来代表所有不等价的 k\kvec,这个区域就是第一布里渊区(first Brillouin zone)。

这一步的物理意义在于:它说明在周期结构中,物理上不等价的是波矢模掉倒格矢后的等价类,而非绝对波矢。后面我们画能带、讨论 Γ\Gamma 点和高对称路径时,都是在这个意义下工作的。

约化区方案与扩展区方案

在固体物理和光子晶体文献中,能带图有两种常见的画法约定:

description [简约区方案 (Reduced Zone Scheme)] 将所有能带都折叠到第一布里渊区内显示。这是最常用的方式,因为它紧凑且突出了布里渊区边界处的能隙 opening。当你看到一条连续的色散曲线在边界处似乎"断开"并跳到另一支能带时,那就是折叠效应的结果。

[扩展区方案 (Extended Zone Scheme)] 允许波矢超出第一布里渊区,在不同区域内画出相应的能带分支。这种方式更能反映自由粒子色散关系的连续性,但在比较不同材料的能带时不如简约区方便,因为你需要时刻记住哪些区域是等价的。 description

本书默认采用简约区方案,除非特别说明。理解这两种方案的区别对于正确解读文献中的能带图至关重要。

初学者在读图时还会遇到另一个常见困惑:为什么有时文献把所有波矢都折回第一布里渊区,有时又把色散画成跨越多个布里渊区的连续曲线?

这对应两种完全等价的记账方式。

阅读时可采用以下顺序:先用约化区方案理解高对称点和模式分类,再用扩展区方案建立“散射把哪几个平面波连到一起”的直觉。后者对耦合波理论尤其重要。

区折叠、Umklapp 过程与 Γ\Gamma 点候选模

要把“画图约定”升级成“物理机制”,关键是把扩展区中的波矢显式映射回第一布里渊区。对任意扩展区波矢 kext\kvec_{\mathrm{ext}},总能找到某个倒格矢 G\gvec 使

kred=kextG,kredBZ1.\kvec_{\mathrm{red}}=\kvec_{\mathrm{ext}}-\gvec,\qquad \kvec_{\mathrm{red}}\in \mathrm{BZ}_1.(4.5)

式 (4.5) 就是区折叠(zone folding)的数学表达:不同扩展区分支在折回后可能落到同一个约化区点,从而在该点发生简并或近简并重组。

这一过程的物理图像是 Umklapp(折返)散射:晶格周期性允许“动量守恒到一个倒格矢”。最小写法是

kout=kin+G,\kvec_{\mathrm{out}}=\kvec_{\mathrm{in}}+\gvec,(4.6)

其中 G\gvec 可由晶格提供(或吸收)。因此在折叠语境下,看似“不同方向/不同区”的平面波分量,实际上可能在约化区内属于同一耦合家族。

这也是 Γ\Gamma 点模式形成的核心入口:当扩展区中的若干波分量经过 式 (4.5) 同时映射到 Γ\Gamma(即 kred=0\kvec_{\mathrm{red}}=\bm{0})附近时,它们会在周期耦合下重组为一组区中心候选模。对 PCSEL 来说,这组模式之所以关键,是因为区中心模式更容易与法向辐射通道建立相干叠加关系。

对常见方格晶格 PCSEL,可把二阶 Bragg 工作点直观理解为“主导导模波数与一阶倒格矢同量级”,即

βslabG=2πa.\beta_{\mathrm{slab}}\sim |\gvec|=\frac{2\pi}{a}.(4.7)

在该条件附近,±x/±y\pm x/\pm y 方向的主导平面波分量经区折叠与周期耦合汇聚到 Γ\Gamma 邻域,形成后续第5章讨论的区中心模家族与模式选择问题。

二维晶格、倒晶格与布里渊区

图 4.3 展示了二维实空间晶格和倒空间晶格的概念示意。图 4.2进一步给出了正方晶格的完整实例,包括实空间原胞构造(左)和对应的第一布里渊区及高对称路径(右)。

正方晶格的实空间原胞与对应的第一布里渊区示意图。
正方晶格的实空间原胞与对应的第一布里渊区示意图。(a)左图给出晶格基矢、散射体位置以及Wigner-Seitz原胞;(b)右图给出倒空间中的第一布里渊区、高对称点 Γ\GammaXXMM 及常用能带路径 ΓXMΓ\Gamma\to X\to M\to\Gamma图 4.2 · 原始 LaTeX/TikZ 图源独立渲染

表 表 4.1总结了PCSEL中最常见的三种二维晶格类型及其基本性质。

表 4.1 常见二维晶格类型及其性质
晶格类型正方三角蜂窝
基矢a1=ax^\bm{a}_1=a\hat{x}a1=ax^\bm{a}_1=a\hat{x}a1=ax^\bm{a}_1=a\hat{x}
a2=ay^\bm{a}_2=a\hat{y}a2=a2(x^+3y^)\bm{a}_2=\frac{a}{2}(\hat{x}+\sqrt{3}\hat{y})a2=a2(x^+3y^)\bm{a}_2=\frac{a}{2}(\hat{x}+\sqrt{3}\hat{y})
倒格矢正方形六角形六角形
第一布里渊区正方形正六边形正六边形
高对称点Γ,X,M\Gamma, X, MΓ,K,M\Gamma, K, MΓ,K,M\Gamma, K, M
点群C4vC_{4v}C6vC_{6v}C6vC_{6v}
典型应用大多数PCSEL三角晶格激光器光子石墨烯

实空间基矢告诉我们孔阵在物理空间如何重复,倒空间基矢则告诉我们哪些傅里叶分量最先参与散射与耦合。

更正式地,二维倒格基矢由

aibj=2πδij\bm a_i\cdot \bm b_j = 2\pi\delta_{ij}(4.8)

定义。对任意二维基矢 a1,a2\bm a_1,\bm a_2,可写成

b1=2πz^×a2a1×a2,b2=2πa1×z^a1×a2.\bm b_1 = 2\pi\frac{\hat{\bm z}\times \bm a_2}{|\bm a_1\times\bm a_2|}, \qquad \bm b_2 = 2\pi\frac{\bm a_1\times \hat{\bm z}}{|\bm a_1\times\bm a_2|}.(4.9)

式 (4.9) 是后面处理三角晶格时应反复调用的工具,而不是只对正方晶格记一个“旋转九十度再乘 2π/a2\pi/a”的特殊结果。

二维实空间晶格与倒空间晶格示意。
二维实空间晶格与倒空间晶格示意。图 4.3 · 原始 LaTeX/TikZ 图源独立渲染

第一布里渊区的构造最好也不要只记最终图形。它的定义是:在倒空间中取所有比其他倒格点更接近原点 Γ\Gamma 的点所组成的区域,也就是倒空间里的 Wigner--Seitz 原胞。实际作图时,最稳妥的步骤是先画出离原点最近的一圈倒格点,再对连接原点与这些倒格点的线段作垂直平分线,最后取围住原点的最小闭区域。

方格晶格第一布里渊区的构造示意。蓝色区域是围住原点的 Wigner--Seitz 原胞。
方格晶格第一布里渊区的构造示意。蓝色区域是围住原点的 Wigner--Seitz 原胞。图 4.4 · 原始 LaTeX/TikZ 图源独立渲染

图 4.5 则给出了标准方格晶格高对称点标注。中心点 Γ\Gamma 对应 k=0\kvec=\bm{0},而边中心 XX 与角点 MM 分别代表两类不同的特殊相位条件。后面讨论 PCSEL 时,Γ\Gamma 点之所以格外重要,是因为它对应法向辐射、区中心模以及多组等价平面波分量之间的特殊耦合结构。

方格晶格第一布里渊区与标准高对称点标注。常用路径为 Γ X MΓ。
方格晶格第一布里渊区与标准高对称点标注。常用路径为 ΓXMΓ\Gamma\rightarrow X\rightarrow M\rightarrow\Gamma图 4.5 · 原始 LaTeX/TikZ 图源独立渲染

从平面波展开理解能带为什么会出现

如果把上一章的 Maxwell 方程和本章的 Bloch 形式合并起来,可以把周期介质中的模式进一步写成

Ek(ρ,z)=Gck+G(z)ei(k+G)ρ.\E_{\kvec}(\rho,z) = \sum_{\gvec}\bm{c}_{\kvec+\gvec}(z)\, \ee^{\ii(\kvec+\gvec)\cdot\rho}.

这个展开式说明:在周期结构中,一个 Bloch 模并非单一平面波,而是许多相差倒格矢的平面波分量的相干叠加。若没有周期调制,这些平面波彼此互不耦合;一旦周期调制存在,它们就通过介电函数的傅里叶分量彼此混合,于是原本简单的自由色散被“切开”和“重排”,这就是能带结构出现的根本原因。

这一步容易被误解为单纯的“频率对波矢作图”。实际上,能带是平面波在周期介质中相互耦合后的本征值图谱。理解这种耦合图像后,“带隙”“带折叠”“带边”和“Gamma 点简并”都可回到“哪些平面波被耦合在一起”这一问题上。

对一维周期结构,最基本的带隙可以看成两个满足布拉格条件的反向平面波发生耦合后形成的分裂。二维 PCSEL 虽然复杂得多,但逻辑是一致的:只是参与耦合的波分量不再只有一对,而可能是一组等价方向的波。

二维光子晶体的常见晶格:方格与三角晶格

PCSEL 文献中最常见的二维晶格是方格晶格和三角晶格。选择哪一种,不是审美问题,而会直接影响:

方格晶格的优势是几何和对称性直观,很多最小耦合波图像都容易从这里建立;三角晶格则往往在六重对称条件下展现不同的模式简并与辐射特征。对初学者而言,先在方格晶格上建立直觉,再迁移到三角晶格,是比较稳妥的学习路径。

三角晶格速览:后面为什么会频繁出现 KK 点与 C6vC_{6v}

三角晶格的一个常用基矢选择是

a1=ax^,a2=a(12x^+32y^).\bm a_1 = a\hat{\bm x}, \qquad \bm a_2 = a\left(\frac{1}{2}\hat{\bm x}+\frac{\sqrt{3}}{2}\hat{\bm y}\right).(4.10)

将它代入 式 (4.9),得到

b1=2πa(1,13),b2=2πa(0,23).\bm b_1 = \frac{2\pi}{a}\left(1,-\frac{1}{\sqrt{3}}\right), \qquad \bm b_2 = \frac{2\pi}{a}\left(0,\frac{2}{\sqrt{3}}\right).(4.11)

此时第一布里渊区是正六边形,最常见的高对称点是区中心 Γ\Gamma、边中心 MM 和顶点 KK。在标准笛卡尔坐标系下,这些点的坐标通常取为

Γ=(0,0),M=(πa,π3a),K=(4π3a,0).\Gamma = (0,0), \qquad M = \left(\frac{\pi}{a},\frac{\pi}{\sqrt{3}a}\right), \qquad K = \left(\frac{4\pi}{3a},0\right).(4.12)

注意 MM 点是第一布里渊区六边形边的中点,而 KK 点是顶点。有时也会见到用倒格矢基底的约化坐标表示,例如 K=(1/3,1/3)K=(1/3,1/3),这取决于具体约定。这一节的任务是给后文一个最小几何坐标系:当 第 8 章第 9 章 讨论 C6vC_{6v} 小群、双重简并和六重对称散射时,读者至少知道这些标签对应的倒空间位置是什么,而非立刻把三角晶格全部算透。

三角晶格第一布里渊区示意。最常用的高对称路径是 Γ M KΓ。
三角晶格第一布里渊区示意。最常用的高对称路径是 ΓMKΓ\Gamma\rightarrow M\rightarrow K\rightarrow\Gamma图 4.6 · 原始 LaTeX/TikZ 图源独立渲染

初读高对称路径时应看什么

很多书或软件会直接给出一条高对称路径,比如方格晶格中的 ΓXMΓ\Gamma\rightarrow X\rightarrow M\rightarrow\Gamma。常见误区是把注意力全部放在“哪条曲线高、哪条曲线低”,却忽略了更重要的问题:

  1. 哪些点是高对称点,为什么只看这些点;

  2. 某条带在这些点附近是否变平,意味着哪些耦合增强;

  3. 哪些带在高对称点重合,意味着可能存在简并与后续可控分裂。

因此,第一次看高对称路径,不要急着评价“哪条带更好”,而要先学会读出“哪里在发生模式重组”。

若在方格晶格中只保留最主要的四个平面波分量,许多 Γ\Gamma 点耦合现象就能得到很清楚的图像。这也是为什么很多教材和早期解析模型喜欢从方格晶格开始,而不是一上来就讨论所有可能晶格。

把二维理想图像放回真实 PCSEL:slab 结构会带来什么变化

到目前为止,我们讨论的 Bloch 理论主要是“平面内无限周期”的语言。但真实 PCSEL 是具有有限厚度的 slab 结构,而非严格二维材料。也就是说,场不仅在平面内受周期调制,还在 zz 方向受外延层栈约束。这一事实会带来两个关键后果。

第一,实际模式是全矢量三维模式。即使在某些参数区间内人们习惯说“TE-like”或“TM-like”,那也通常意味着“主要电场分量或主要磁场分量具有某种偏向”,而不是二维教科书中严格分离的纯 TE/TM 解。

第二,slab 的存在允许导模、泄露模和辐射模共存。某个在二维理想模型中看似被很好分类的 Bloch 模,放到有限厚度 slab 中以后,可能会通过垂直辐射通道与自由空间耦合,从而成为 PCSEL 需要追踪的候选模。

二维 Bloch 理论的作用是告诉我们“平面内有哪些组织良好的模式家族”;slab 结构则提醒我们“这些家族在真实三维器件里会怎样与纵向约束和辐射通道重新混合”。PCSEL 正是在这两层之间工作。

TE/TM 分类为什么常常只是近似

在严格二维、无厚度、标量化处理的模型中,TE 与 TM 可以作为清晰的极化分类。但在真实 PCSEL slab 中,下列因素都会破坏这一理想图像:

因此,后文若使用“TE-like”或“TM-like”这样的词,应理解为工程上的近似标签,而不是严格数学分类。本书采用一个可操作的判别标准:若模式电场能量中 EzE_z 分量占比 <5%<5\%,称为 TE-like;反之若 HzH_z 分量占比 <5%<5\%,称为 TM-like。这里的“like”提醒读者:只要 slab 上下不对称,EzE_zHzH_z 都不可能严格为零,真实模式是混合极化的。这个提醒很重要,因为许多关于偏振和远场的误判,都源自把近似 TE/TM 标签当成了真实矢量模的完整描述。

一个对初学者更安全的学习顺序

如果读者刚接触 PCSEL,不必一上来就放弃所有 TE/TM 直觉。更安全的做法是分三步:

  1. 先用二维标量或 TE-like 图像理解带结构、带隙和高对称点;

  2. 再用全矢量 slab 模型检查哪些结论仍然成立,哪些开始偏移;

  3. 最后在有限尺寸和开放边界模型中验证偏振与远场。

这种“先建立直觉,再主动拆掉过强近似”的路径,比一开始就把所有三维矢量细节压到新手身上更有效。

CommonPitfallBox把二维 TE/TM 当成三维 slab 的严格本征极化 常见误读:在二维模型里标了 TE/TM,就直接把它当成真实器件中的“纯 TE/纯 TM 模”。

正确做法:把二维结果仅作为 TE-like/TM-like 的初筛标签,再用全矢量三维模型(含开放边界与有限尺寸)复核偏振、远场和阈值排序。

典型后果: itemize 把混合极化模误判为“纯偏振模” 误判法向辐射允许性与偏振选择规则 误把二维排序当成器件阈值排序 itemize CommonPitfallBox

Bloch 理论能告诉我们什么,又不能告诉我们什么

到这一章结尾,必须把 Bloch 理论的能力边界说清楚。

它能告诉我们: itemize 周期结构中模式如何按 k\kvec 和对称性分类; 第一布里渊区内哪些高对称点值得重点关注; 某些简并与选择规则为何出现。 itemize

它不能单独告诉我们: itemize 有限器件边缘效应会不会改变模式排序; 有源区增益分布是否偏向目标模; 电注入和温升下哪一个模最终先起振并保持稳定。 itemize

因此,Bloch 理论是模式分类的起点,而不是 PCSEL 器件结论的终点。

通向下一章:为什么 Gamma 点成为重点

本章已经说明:二维周期结构中的模式可以用 Bloch 波矢分类,而 Γ\Gamma 点是布里渊区中心的特殊位置。下一章将进一步回答:为什么在 PCSEL 中,Γ\Gamma 点附近的模式经常成为重点关注对象?简并、对称性、辐射耦合和有限尺寸修正又会怎样共同决定候选激射模?

回看本章

Bloch 定理把无限周期结构的求解组织成了按波矢分类的问题,并由此引出倒格矢、布里渊区和高对称点等核心概念。对 PCSEL 而言,这一步的意义在于建立候选模的“目录系统”。但这还只是目录,目录之后还要讨论哪些模式与辐射、增益和器件工作条件相匹配。

练习题

enumerate [基础] 对方格晶格和三角晶格,分别写出一组可能的实空间基矢与倒格矢基矢。

[基础] 从平移算符的角度,重述一次 Bloch 定理的证明思路。

[基础] 证明 k\kveck+G\kvec+\gvec 对应相同的 Bloch 相位条件。

[进阶] 解释为什么 Γ\Gamma 点的重要性不仅仅来自“群速度小”,还来自模式对称性和辐射通道。

[综合] 给出一个 slab 结构中 TE/TM 近似可能失效的具体场景。 enumerate

延伸阅读

itemize [Joannopoulos et al., 2008] [Sakoda, 2005] 中关于 Bloch 定理与二维光子晶体的基础章节。

读完下一章后再回看本章,通常会更容易理解 Γ\Gamma 点为何成为 PCSEL 模选择的核心位置。 itemize

PCSEL KNOWLEDGE BASE

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

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

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