1. 带隙是什么光在周期结构里为何会被“拒之门外”先说点我自己的经历。第一次接触二维光子晶体带隙仿真的时候我满脑子都是“不就是算个能带嘛和固体物理里的电子能带差不多”结果第一次跑程序就栽了跟头——结构画对了材料参数填对了可出来的能带图里干干净净一个带隙都没有。当时的导师看了眼我的输入文件问了一句“你的介电常数对比度是多少”我才发现自己随手填了一个接近1的值等于让光在一个几乎没有折射率差别的“假周期结构”里传播当然不会有什么带隙。后来我才真正理解光子晶体带隙的本质是周期性介电函数对电磁波模式的重新分配。在均匀介质里光可以以任意频率传播色散关系是一条直线。但当你把介质做成周期排布并且周期的尺度接近波长量级情况就完全不同了。结构不再只是“让光折射”的普通材料它变成了一张频率的筛子——某些频率范围内的光无论从哪个方向入射都无法在结构内部找到允许传播的模式这就是我们说的完全光子带隙也叫完全禁带。对应到能带图上横轴是波矢在布里渊区里的路径纵轴是归一化频率通常是 (\omega a/2\pi c)也就是 (a/\lambda)。带隙就是从某个能带的最低点到下一个能带最高点之间的空白区域。这条空白的“宽度”直接决定了器件的性能上限带隙越宽能容纳的缺陷态频率调节范围越大波导转弯的损耗越小腔体的品质因子也能做得越高。这里必须区分一个概念不完全带隙和完全带隙。有些结构只在特定偏振下存在带隙比如只对TE模式电场垂直于周期平面的模式有禁带对TM模式没有。只有一个偏振方向带隙的结构做波导、做分束器一样能玩出花来但如果你想要那种“不管什么偏振的光都传不进去”的全带隙就得让结构在TE和TM两个偏振下同时打开带隙。而二维光子晶体最经典的魅力恰恰在这里三角晶格空气孔结构在合适的空气填充率下两个偏振都能开带隙这在很多三维结构里是可遇不可求的。我在仿真里最常用的一句话是“带隙不是算出来的是设计出来的。”这句话有两层意思。第一层带隙的宽度和位置取决于结构的对称性、晶格常数、空气孔的半径、背景材料的介电常数这些都是你可以在仿真前通过理论估算预判的。第二层仿真只是把那道频率筛子精确地画给你看帮你验证设计而不是替你发明结构。所以搞清楚带隙的物理来源比急着把仿真软件跑起来重要得多。2. 由布洛赫定理到布里渊区仿真前必须建立的三层认知2.1 从电子带隙到光子带隙同一个数学骨架如果你学过固体物理看到光子晶体能带图的时候一定会觉得眼熟。没错它的数学骨架和电子能带理论几乎完全一致。电子在周期势场中的波函数满足布洛赫定理可以写成一个平面波乘以一个周期函数光子在一张周期介电函数“地形图”里传播电场和磁场的分布同样满足布洛赫形式。你只需要把薛定谔方程换成麦克斯韦方程组把周期势垒换成周期介电常数剩余的一切——倒格子、布里渊区、能带折叠——全部照搬。这种类比极其有用。比如你可以问一个问题为什么光子晶体只能在一部分结构里实现完全带隙而电子晶体几乎都有带隙答案藏在组份的对比度里。电子遇到的是原子势垒深度和周期往往处在同一个量级光子遇到的则是介电常数差。当两种介质的介电常数相差很小比如1.5和2.0周期结构的“散射强度”太弱不同波矢的模式之间几乎没有耦合带隙就不会打开。所以仿真之前第一件事就是确认材料体系的介电常数对比度足够大一般至少要背景和孔的介电常数比值超过5才有可能在二维结构里看到像样的完全带隙。2.2 倒格子、布里渊区、高对称点能带图上那条横轴的意义很多新手拿到能带图只看有没有带隙却不知道横轴是哪里来的。实际上那条横轴是布里渊区的边界路径。二维正方晶格的第一布里渊区是一个正方形高对称点是 (\Gamma)、(X)、(M)三角晶格的第一布里渊区是一个六边形高对称点是 (\Gamma)、(M)、(K)。能带图的横轴一般是沿着 (\Gamma \to X \to M \to \Gamma) 或者 (\Gamma \to M \to K \to \Gamma) 走一圈把所有方向上可能存在的模式都包含进去。为什么要沿高对称点的路径扫描因为周期性结构色散关系的高点和低点往往出现在布里渊区的高对称点上或其附近。完全带隙判据是在整条路径上第一个能带的最高频率仍然低于第二个能带的最低频率。你扫描的路径越完整越能确认这个带隙在“任意传播方向”都存在。如果只扫一个方向而漏了其他方向就可能把方向带隙误判成完全带隙后续做波导时才会发现光在里面沿着某些方向根本传不动。我自己的习惯是第一轮计算先把高对称点路径全部扫一遍确认存在带隙后再看能带边缘对应的模场分布。因为能带的带边valence band maximum 和 conduction band minimum对应的模式决定了带隙两端的物理特性比如慢光效应、负折射的起始频率都跟带边模式强相关。2.3 晶格选择与空气填充率带隙拓扑的直觉来源二维光子晶体最常见的两种构型空气孔在介质背景里介质柱在空气背景里。 空气孔结构背景是高折射率材料孔里是空气适合做W1波导、微腔因为光主要在背景介质里传播缺陷容易刻成波导和腔体介质柱结构背景是空气柱子里是高折射率材料更容易形成远离互连的孤立介质波导适合做拓扑边界态等研究然而在工艺上柱子容易倒。决定带隙的“第一参数”不是晶格常数 (a)而是空气填充率 (f)。对于三角晶格空气孔结构填充率定义为圆孔面积与单胞面积之比(f \frac{2\pi}{\sqrt{3}}(r/a)^2)。当 (r/a) 在0.3到0.45之间TM模式会打开一个较宽的带隙当 (r/a) 更小比如0.25以下TE模式带隙可能占主导。想要TE和TM同时打开完全带隙通常需要把 (r/a) 调到0.35-0.4附近并且背景介电常数不低于12比如硅、砷化镓都行。这些数值不是我拍脑袋定的而是二维光子晶体研究里反复验证过的经验窗口。我的建议是仿真前用简单的公式估算一下填充率然后把 (r/a) 在0.25到0.45之间扫一遍。这个参数扫描是各光子晶体仿真课题组的“入学第一课”因为带隙宽度对 (r/a) 的变化极其敏感差0.02的半径就能让带隙从无到有再从有变无。3. 选对计算工具平面波展开法与有限时域差分法的真实分工3.1 平面波展开法能带图的首选但不是万能钥匙平面波展开法Plane Wave ExpansionPWE是我做大面积参数扫描时最依赖的工具。它的思路很直白既然结构是周期的介电函数 (\varepsilon(\mathbf{r})) 就能展开成傅里叶级数电场/磁场模式也展开成一组平面波的叠加。把麦克斯韦方程组变换成倒格子空间里的矩阵本征值问题求解本征频率就能得到所有能带。PWE的优点非常明显不需要离散化整个实空间结构直接利用周期性速度极快几秒钟就能算一条带隙边界。所以在扫描填充率、扫描晶格类型、扫描折射率这几类任务上PWE几乎不可替代。但PWE有个致命短板——它要求结构严格周期性处理不了缺陷。你可以在超胞近似下用它算带缺陷的态密度但精确度受超胞大小限制一旦缺陷模式与超胞边界杂散耦合结果就会变得不可靠。3.2 FDTD缺陷模式与透射谱的真正落脚点时域有限差分法FDTD则是另一套哲学。它直接把麦克斯韦方程组在时空网格上离散差分用时间步进的方式模拟电磁波在结构里的真实演化。它的优点是不受周期性假设限制可以随便定义缺陷、波导拐弯、耦合器还能加吸收边界条件来模拟开放空间。带隙仿真里FDTD最常用的场景是计算带隙结构的透射谱验证PWE预言的概率带隙在实空间里真实存在计算非周期缺陷结构的本征模式比如微腔的谐振频率和Q值模拟光脉冲在光子晶体波导里的传播过程。FDTD的代价是速度和内存。三维仿真需要非常密的网格才能解析纳米级孔壁每一步时间推进都牵动全部网格点的电磁场更新跑一次宽谱脉冲仿真动辄需要几十分钟甚至更久。所以对于“找带隙”这个目标用FDTD效率极低但一旦你的带隙结构开始变成具体器件FDTD又绕不开。我的经验是把PWE当筛选器把FDTD当验证器两者结合可以省掉无数时间。还有一类工具是基于有限元法FEM的商软仿真在几何灵活性上最强能处理任意形状的缺陷。但求解本征值问题时有限元对网格质量的要求比PWE和FDTD都苛刻网格剖分不当容易产生伪模。如果你只是想快速出带隙图我的优先级是PWE优先FDTD验证FEM做复杂几何的兜底。3.3 参数归一化所有单位制的换算口诀这是二维光子晶体仿真里最容易忽略的环节。因为晶体结构完全由几何比例和介电常数决定没有固定的绝对尺寸仿真里所有长度、频率都是归一化的。晶格常数 (a) 是一把“尺子”给定结构比例改变 (a) 相当于整体等比缩放带隙在归一化频率轴上的位置不变但绝对频率会整体移动。换算关系非常直接[ \lambda_{\text{带隙}} a \times \frac{c}{\omega_{\text{归一化}}} ]如果仿真得到归一化带隙范围是 (0.25 a/\lambda 0.35)而你想让它落在可见光波段比如 (\lambda 600) nm附近那么反推 (a) 应取150-210nm如果想落在通信波段1550nm那么 (a) 大约要取390-540nm。很多新手把结构画好了却忘了归一化频率和绝对波长的换算导致最后实验测不到想要的透射峰。这个换算值得在仿真前就写在草稿纸上。4. 第一个实战案例三角晶格空气孔结构的全流程仿真4.1 结构定义与初始参数为了把流程讲完整我用一个具体案例从头到尾走一遍所有参数都可以直接复现。基片选用折射率 (n3.46) 的硅介电常数约12.1空气孔按三角晶格排列背景介质是硅圆孔里是空气折射率1.0。晶格常数记为 (a)圆孔半径取 (r0.35a)。这个比例下空气填充率 (f\frac{2\pi}{\sqrt{3}}(0.35)^2 \approx 0.44)是典型的完全带隙工作区间。仿真目标分两步先用平面波展开法拿到TE与TM能带图确认完全带隙位置再用FDTD构造一个厚度有限的平板结构计算透射谱验证带隙在实空间里的存在性。实际操作中很多人一上来就用FDTD算透射谱结果结构大、网格密、跑半天最后得到一条看不懂的透射曲线。正确顺序永远是先PWE后FDTD先用秒级计算确认频段再花时间做精细验证。4.2 平面波展开法的计算流程我用一台普通办公电脑在开源的MPB框架下就能完成计算。流程大概是四步定义晶格常数和晶格类型。三角晶格选择hexagonal lattice并指定基矢长度。定义几何体单胞里放一个半径为 (0.35a) 的圆柱材料设为空气背景指定为硅。程序会自动把结构在倒空间展开成傅里叶系数。设置k-path路径(\Gamma \to M \to K \to \Gamma)每个段取20到30个k点采样。分别计算TE和TM偏振。TE模式下电场沿z方向垂直于周期平面的方向TM模式下磁场沿z方向。很多人第一步会忘记程序里TE与TM的定义和习惯相反一定要先确定你关心的偏振方向。跑完输出的能带图里如果TE和TM都能看到一个空白区带隙恭喜你结构选对了。三角晶格、 (r/a0.35)、硅背景这个组合下TE带隙大约落在 (a/\lambda 0.26 \sim 0.36)TM带隙大约落在 (a/\lambda 0.22 \sim 0.40)二者重叠区域就是完全带隙的候选范围。4.3 FDTD的透射谱验证PWE给出的是无限大周期性结构的本征模式而真实器件尺寸有限。为了验证带隙的真实性我把结构切成有限大小在x方向排列15排空气孔y方向用周期边界条件入射光沿x方向传播从结构一侧注入宽谱脉冲另一侧记录透射场。这种做法相当于把光子晶体当一维光栅用观察它是否在带隙频段“堵住”光。透射谱出来之后你会清晰看到在某些归一化频率区间透射率跌到背景值的千分之一以下这就把PWE的带隙在透射层面验干净了。还需要顺手做一步用FDTD的“本征模求解”功能比如在商软里用Modesolver或开源框架的Harminv工具去提取带隙中的缺陷态频率为之后的波导设计做准备。FDTD网格我习惯设为最小特征尺寸的1/10。本案例圆孔边缘到背景都要用网格刻画所以空间步长取 (a/20) 以下总网格数量控制在千万量级。如果你机器的内存不够可以把周期方向和传播方向都用对称边界条件切成一半只算四分之一结构速度提升四倍还不会损失精度。4.4 带隙边界定得越准后续器件越好做仿真结束并不意味着工作结束。我会额外做一件事追踪带隙带边的模场分布。等频率接近带隙上边界的模式其电场往往强烈局域在低介电区域这决定了它在缺陷结构里的束缚能力接近带隙下边界的模式则局域在高介电区域。理解带边模式的空间分布能帮你预测把缺陷刻在哪种位置会形成束缚态——如果你希望微腔模式局域在空气孔中心那就应该从带隙上边界的模式出发设计缺陷。5. 从能带到缺陷态光子晶体器件设计的三条务实路径5.1 缺失孔与线缺陷波导与微腔的本质带隙是“光的禁地”而缺陷态则是禁地里的孤岛。你从结构中拿掉一个空气孔周期被打破带隙里就会涌现出一个局域的本征模式频率落在带隙重叠区间内。模式被周围光子晶体“全反射”困住无法逃逸就形成了一个微腔。基于这个原理所有实际光子晶体器件——波导、微腔、分束器、滤波器——本质上都是带隙里缺陷态的重新排列组合。最经典的是W1波导拿掉一排空气孔带隙里就出现一条波导模式色散线。这条色散线的斜率和位置可以通过调整波导的宽度、旁边的孔大小来调节。如果你需要慢光群速度极低的波导就调整缺陷区域的几何参数让色散线在目标频率处变得非常平坦。这里面的设计自由度比普通介质波导大得多但前提永远是你的完全带隙足够宽给缺陷模式留出足够的调节空间。5.2 TM/TE特异性设计的互补场景完全带隙大于一切不一定。很多优秀器件只需要单一偏振。比如TE偏振的W1波导在硅基平台上能实现极低的传播损耗因为TE模式的带隙在三角晶格空气孔结构里特别稳定设计余量大。TM偏振的优势是模式不容易受加工粗糙度影响对垂直入射光更友好常用于自由空间光入射/出射的滤波器设计。我的建议是先明确器件最终要跟哪一种光源和探测器配合。如果后续要跟片上激光器、探测器对接TE模式通常是首选如果要做自由空间光传感器或垂直耦合器TM模式更省事。别为了“完全带隙”的虚荣值牺牲偏振匹配的实际需求。5.3 参数扫描的技巧别浪费算力在无效区间参数扫描是光子晶体仿真的常规操作但很容易变成无脑循环。我在跑扫描前会先做两道手算估算一下填充率对带隙中心频率的影响方向。孔变大会让有效介电常数变小不利于低频带隙的支撑带隙整体向高频偏移孔变小则相反。确认目标频率与晶格常数的换算避免扫出来的带隙落在完全用不到的波长区间。扫描维度通常包括空气孔直径、背景折射率、晶格常数比、缺陷区几何宽度。每一轮扫描我都建议只动一个变量保持其他参数完全不变否则一维扫描就会变成多维模糊搜索根本没法定位关键参数。如果要同时扫两个维度可以在算完整个二维矩阵后做等高线图找带隙宽度的“高原区”外围边界那是最稳健的设计工作点——尽量远离参数悬崖这样加工误差对带隙的影响能降到最低。5.4 从二维到准三维有限厚度的教训完全二维仿真在数学上是严格的但真实芯片上的光子晶体膜都有厚度通常在100到500nm量级。当你把二维结构“挤压”成一块薄膜里排列空气孔光会被限制在平面内吗答案是可能漏射到空气中。这引出准三维仿真里的一个重要参数薄膜厚度。太薄的膜限制能力弱模式辐射损耗严重太厚的膜又会激发出高阶垂直模式破坏单模工作条件。我自己踩过这个坑二维仿真预言带隙里有个高Q微腔但在实际加工测试中怎么也测不到。后来准三维仿真发现模式在垂直于平板的方向上没有束缚好能量从上下表面漏掉了。从那以后凡是要拿去加工的结构我一定会先做一次3D仿真或至少用有效折射率法把厚度折算进去再决定是否交付流片。6. 盲区与误判带隙仿真最容易翻车的四个细节6.1 网格分辨率与傅里叶阶数的骗局PWE里有一组隐藏参数平面波的数量。平面波截断数太小能带边缘会被人为拉低带隙宽度容易被高估。FDTD里也有网格密度的“骗局”网格太粗带隙两端的边缘频率会向中间靠拢带隙看起来比真实情况窄。这两种误差方向完全相反所以你在论文里对比两种方法的结果时不能指望它们严丝合缝地一致而应该确认它们给出的带隙范围存在足够的重叠。我的经验是先用粗网格跑一遍再加密一倍网格重跑如果带隙边界漂移小于1%说明计算收敛了如果漂移大必须继续加密。PWE平面波数量的判断同理。任何“第一次跑就出漂亮带隙”的结果都要警觉先做收敛性测试再说。6.2 模式遗漏与近简并误判PWE输出的能带图会按频率从低到高排列如果想确定带隙严格办法是比较第N条和第N1条能带在整个k路径上的极值。但有时两条能带非常接近几乎要交叉你在能带图上肉眼看起来像有一个小带隙实际却是假简并。出现这种情况要回到模场分布确认两条能带的模式对称性——只有本征模式对称性不允许耦合时色散线才会出现真正的“反交叉”。近简并区往往是即将打开带隙却因参数不足而没有打开的先兆所以有时几何参数微调一点误判为带隙的缝隙就会彻底关闭。6.3 边界条件对缺陷模拟的干扰用超胞方法计算带缺陷的状态时边界效应是个隐蔽问题。超胞小了缺陷模式和相邻超胞之间的耦合会让频率漂移超胞大了计算成本骤增。我的做法是试一个超胞尺寸把缺陷态频率记为 (f_1)把超胞尺寸扩大一倍再记一个频率 (f_2)。如果 (|f_1-f_2|/f_1 0.1%)说明超胞足够大否则继续扩大。这个过程虽然费时间但它是获得可靠缺陷态频率的唯一办法。6.4 材料损耗被忽略的后果硅在通信波段的吸收损耗很低但有些光子晶体设计用的是聚合物或掺杂半导体材料存在明显吸收。带隙仿真通常假设材料无损耗这会给Q值算出一个虚高的理论值。如果后续要跟实验对照必须把复介电常数加进去看看模式频率是不是真的还在带隙内。有时候频率虚部的增加会让模式“渗出”带隙出现辐射损耗这个效应在二维仿真中完全看不到只有在添加了材料损耗的三维仿真中才会暴露。7. 从仿真到流片带隙设计走向实验的几个现实问题7.1 工艺误差的名义值与真实值仿真时参数都是精确的——半径0.35a就是0.35a。但实际加工里光刻和刻蚀的直径会有一到几纳米的偏差。你可能费尽心思把工作点设计在带隙边缘结果加工后半径偏移2%模式直接从带隙边缘掉到连续态里Q值一落千丈。所以在评估设计鲁棒性时至少要做一次半径偏移±2%的仿真看带隙宽度和缺陷模式频率是否还能接受。7.2 垂直侧壁的倾斜与粗糙度二维光子晶体仿真模型中的空气孔侧壁是垂直、光滑的。但实际工艺里反应离子刻蚀会带来侧壁倾角侧壁粗糙度也会造成散射损耗。侧壁倾斜会把圆孔变成锥形有效填充率随深度变化带隙位置会发生漂移。如果你手里有SEM照片里的侧壁角度数据强烈建议把它折算到仿真模型里哪怕只是一个锥形孔的近似。如果目标器件Q值需要达到百万以上侧壁粗糙度几乎是决定性瓶颈。这个问题在二维仿真里看不见只能通过设计宽容度来对冲尽量把工作频率落在带隙中央而非边缘模式场分布尽量避开孔壁或者用空气桥结构减少与衬底的泄漏。7.3 封装与测量的匹配问题仿真里的理想情况是光完全在光子晶体平板内传播测量时则需要光栅耦合器或者锥形波导把外部激光导入结构。耦合器的模式匹配效率很容易让人头疼——光子晶体波导的模式与普通条波导的模式varying很大失配会产生大量反射。常规的做法是在仿真中加入渐变过渡区让条波导的模场逐渐过渡到光子晶体波导的模场这个过程也需要用三维FDTD仔细优化。我在实际项目里最强的体会是二维带隙仿真回答的是“这个结构理论性能能有多好”而三维带隙仿真加上工艺误差模型回答的才是“这个设计在工厂里能做出什么样的成品”。从前者到后者之间的差距就是光子晶体科研和工程的差距。8. 踩坑实录我在这类项目中花费时间最长的排查过程有一次做一个正方晶格介质柱结构的能带仿真用PWE算完没有带隙我以为是几何参数没选对把填充率从0.2一直扫到0.55整个下午都在循环计算没有任何一个参数出现带隙。问题是后来我发现结构里介电常数的定义出现了问题——背景介质被我设成了空气介质柱用了硅但这只是第一种错误更隐蔽的是程序里对于正方晶格介质柱默认的电场偏振定义和我想要的方向相差了90度。也就是说我一直在算两种偏振都设反了的情况等于把TE偏振当TM偏振算了自然扫不出正确的带隙。这件事给我的教训非常深刻仿真软件里的“TE”和“TM”在不同工具里定义并不统一有的以电场方向为准有的以磁场方向为准有的干脆按z方向与结构平面的关系定义。拿到任何新工具的默认配置第一时间跑一个已知解析解的结构比如均匀介质或真空的能带图是一条标准直线验证符号约定然后再开始正式计算。第二个大坑是超胞边界的杂散态。我做W1波导的缺陷模式时PWE算出的缺陷态频率看起来正常但模场分布却出现了一个不合理的“卫星峰”。排查发现超胞尺寸取得刚好是缺陷周围第3排孔的位置而杂散模式恰好和超胞边界周期产生了布拉格反射在带隙中形成了寄生态。换一个更大的超胞后寄生态彻底消失。所以当你看到缺陷模场不是老老实实局域在缺陷附近而是出现周期性分布时十有八九是超胞边界效应而不是什么有趣的物理新现象。第三个坑发生在FDTD透射谱的归一化上。直接对透射场做傅里叶变换得到的透射率包含了光源本身的频谱形状如果不先算一条无结构的参考透射谱再归一化得到的“透射谷”可能是光源频谱的凹陷。这个错误做光学仿真的人经常犯检查方法很简单换一个完全不同中心频率的光源重跑如果透射谱峰值位置不变说明归一化正确如果整体频谱跟着光源走就说明结果被光源污染了。9. 给初学者的入门梯度从复现文献到独立设计的路径如果你刚接触二维光子晶体带隙仿真我不建议你直接照着一篇论文的参数去“复现”因为没有哪本教材能告诉你论文里那些隐藏的默认设置。我的建议是先跑通三件套用PWE计算均匀介质比如每处折射率都是3.46的能带图你应该得到一条直线色散这是一切周期结构计算正确性的“地基”用一个最简单的二维结构比如正方晶格介质柱阵列算带隙然后把文献里的曲线叠上来对比在此基础上加一个缺陷孔看带隙中是否出现缺陷态体验“禁地里挤出一根针”的乐趣。这三个练习全部通过后再进到参数扫描和优化阶段。优化时建议画一个“带隙宽度 vs (r/a)”的曲线你会看到带隙存在一个最大值。找到这个极值点你就掌握了该结构体系的设计重心。之后再去考虑引入椭圆孔、异质结构等更多自由度。很多初学者在这个过程中会掉进“参数越多越高级”的诱惑里一口气把两三种结构调整参数都打开。实际上一个稳定、可解释、可复现的设计永远比一个“看起来很大”却对误差极其敏感的带隙更值钱。做光子和做其他工程一样少即是多。10. 最后分享一个实用的小技巧快速估算仿真规模每次开始仿真前我心里都会做一个简单估算目标结构尺寸若是 (10a \times 10a)FDTD网格取 (a/20)那么网格点数为200乘200每个网格点存若干个电磁场分量三维再乘上厚度方向的网格层数。这样粗略一算就能判断当前机器能不能扛住要不要用并行计算或缩小仿真区域。这个估算还有一个用处帮你决定哪些结构必须用PWE而非FDTD。如果你要扫描几十组参数的三维结构用FDTD逐个跑任务量可能超过几千小时而PWE在同样工作量下只需几分钟。先PWE后FDTD不只是顺序问题还是算力管理问题。二维光子晶体带隙仿真最有意思的地方在于你要在完全禁带的“无路可走”里设计出一条条精巧的“小路”。每一次对参数的小调整都可能打开或关闭那一小片光的频率空间。这种近乎雕琢的乐趣是均匀介质波导完全无法提供的。仿真只是手段带隙只是工具真正让我觉得奇妙的是我们能通过周期结构的微小塑形让原本“必定能传播”的光老老实实按我们规划的路径走。而这恰恰就是光子晶体这门手艺最让人着迷的地方。