谷光子晶体拓扑激光器：从能带计算到光泵浦实测

二维谷光子晶体拓扑激光器这个题目我前前后后跟了快两年。二维光子晶体提供能带设计平台谷自由度给模式加了一个可选择的“取向”拓扑则相当于给激光模式买了一份保险。整套思路一句话可以讲清楚先把光子晶体的能带结构设计成具有谷依赖的拓扑边界态再在这个边界态上引入增益介质让它变成受拓扑保护的单模激光器。它针对的老问题是传统微型激光器对侧壁粗糙度、加工误差和弯曲损耗都很敏感而拓扑边界态对很多局部扰动不敏感模式可以更稳。这篇文章主要面向正在做光子晶体或拓扑光子学研究的同学、以及想把“拓扑”这个概念真正做成器件的工程师我会把从能带计算、器件仿真到光泵浦实测的整个流程和踩坑记录都尽量摊开来讲。1. 二维谷光子晶体拓扑激光器是个什么“物种”1.1 题干拆解谷、光子晶体、拓扑激光各扮演什么角色先把题目拆成三块二维谷光子晶体指的是在二维平面内周期排列的电介质结构比如三角晶格空气孔、蜂窝晶格介质柱通过改变晶格对称性能带会在布里渊区的 K 和 K 点形成能带极值这两个等价的极值点就是“谷”。谷自由度最早在电子体系里被研究得很深比如过渡金属硫族化合物中的能谷电子学后来光子学把它接了过来用光场模拟电子局域在能谷附近的色散和手性行为。拓扑在这里不是“拓扑优化”那个意思而是起源于凝聚态物理的拓扑不变量概念。用最形象的话解释拓扑保护的光模式就像把一段轨道放在地形图里模式在界面上的传播是由体态的拓扑性质决定的即使界面里有小缺陷、拐弯或者局部折射率涨落模式依然不容易被散射到其他通道。拓扑激光器就是把增益介质与这样的边界模式结合让受激辐射发生在一个对扰动不敏感的模式上。这和普通的回音壁微腔激光、法布里-珀罗腔激光不一样。微环激光靠全内反射侧壁粗糙度会带来很大的散射损耗FP腔激光对端面反射率极其敏感稍微有一点刻蚀歪斜阈值就会明显升高。拓扑边界态激光最主要的卖点是光场被“锁”在两种拓扑结构的分界面上传播方向和谷量子数绑定背向散射被抑制所以理论上可以提高激光器的鲁棒性同时保持单模工作。1.2 它到底想解决什么真实问题很多人觉得做拓扑激光器是为了发文章实际上工程上确实有痛点。微型激光器的传统方案里高 Q 值往往意味着精细的特征尺寸和严格的工艺公差。在硅或者三五族半导体上刻空气孔孔的位置偏差几纳米就会改变模式频率和损耗批量一致性很难保证。另一个痛点是弯曲和结构不连续传统波导沟道在拐角处会有辐射损耗必须做弯曲半径优化能绕过一些障碍的紧凑光子回路并不容易设计。谷光子晶体拓扑边界态的出现让这些问题有了新的解决路径。因为模式不是靠全内反射约束而是靠体带的拓扑带隙约束界面上一段 120 度甚至 60 度的急转弯可以保持高透射模式不会轻易“掉出”波导。这样的特性对于集成光路里的光源很有价值激光器不再需要严格放在笔直的腔体里它可以长成不规则的环、折线甚至阵列结构只要拓扑畴的界面保持不变模式就能稳定存在。最后一点从科研角度看这个方向也是验证体边对应原理的很好平台。用光子晶体把拓扑物理翻译成可观测的光谱、近场分布、激射阈值和线宽让抽象的概念变得可以直接测量。1.3 谁适合读这篇总结如果你是做光子晶体仿真的新手前两章能帮你建立从能带到边界态的基本直觉如果你已经在做被动拓扑波导、想往激光方向走第三章和第四章可以直接参考我的建模和测试流程如果你已经推进到器件级第五章的坑可能正是你正在头疼的问题。当然这也不是一篇教科书我不会把每个方程都推一遍但会讲清楚每步操作背后的判断逻辑和使用场景。2. 原理层面三个支柱缺一不可2.1 能带里的“谷”和普通带边有什么区别做一维光子晶体的人都很熟悉带边在布里渊区边界色散变平群速度为零态密度升高这是很多带边激光器的基本思路。二维谷光子晶体的带边更精致它工作在布里渊区角上的 K 点和 K 点。这两个点是六角晶格对称性下的不等价点它们之间存在时间反演对称联系。如果只了解平面波展开的能带图最初看到 K 点和 K 点的模场分布时可能会感觉奇怪K 点在实空间里表现为两种不同手性的旋转场K 点则互为镜像。所以谷自由度天然带有“手性”或“涡旋取向”的含义。就像一个环形山K 点像左旋的台阶K 点像右旋的台阶下山的方向是相反的。在普通各向同性介质里这两个谷是完全简并的一旦我们通过结构形变打破某些镜面对称性K 和 K 会劈裂带隙打开同时这个带隙带有谷锁定的特征。实际设计中打破对称的方式很多。最常见的是把三角晶格空气孔从圆形改成三角形或椭圆并且让这些孔在晶胞内旋转一个角度另一种方式是保持孔形不变但把每个晶胞里的两个子晶格孔半径做得一大一小或者让孔的位置向某一个方向偏移。无论哪种方式本质都是把一个原本具有 C3v 对称性的结构降至 C3 对称从而把退化打开成带隙。2.2 谷陈数拓扑“保险”的具体来源要理解拓扑保护得引入谷陈数作为体带的不变量。在光子晶体里如果把本征模的 Bloch 波矢在 K 点附近绕一圈会积累一个几何相位也就是贝里相位。这个相位对局域参数并不敏感它对远处的结构变化也不敏感只取决于能带在动量空间里的拓扑结构。贝里曲率在 K 谷和 K 谷附近的积分分别对应谷陈数 C_K 和 C_K由于时间反演对称性两者互为相反数通常为 ±1/2 这种半整数。这里的“±”并不是随便给的。它由晶格形变的方向决定比如空气孔旋转的正负、子晶格偏移的方向。两个相邻区域如果 A 区域的谷陈数是 1/2B 区域的谷陈数是 -1/2那么在它们的界面上体带拓扑不变量发生跳变。体边对应原理表明能隙中一定会存在界面模式而且这个模式出现的能量位置、传播方向都和谷陈数跳变一一对应。我当年第一次看到这个对应关系时也很怀疑半整数确实不太像普通的 TKNN 整数陈数。后来在数值上把投影能带画出来看到界面模式被稳稳地夹在体带隙中间才真正接受了这个图像。这个模式还有一个特性在某一个谷附近界面态只能沿着一个方向传播它的群速度锁定在某个方向如果打破界面局域结构它不会轻易背散射到反向这就是我们常说的谷锁定单向传输。用生活类比来说明拓扑保护好比把模式放在了一个很宽的“传送带”上这个传送带的转动方向由两侧晶格的谷极值锁定。小石子落到传送带上会被带着走而不会跳回入口。这就是为什么拓扑激光器对局部扰动不敏感。2.3 边界态怎么变成可以出光的激光器拓扑边界态本身是无源的被动模式要变成激光器需要把增益介质放到模式场分布最强的区域。一般是在晶格空气孔以外的半导体膜里生长多量子阱或者敷设量子点让边界态电场与增益材料有足够大的空间重叠。泵浦光被增益介质吸收后产生电子-空穴对受激辐射放大该频率附近的光场。从模式角度看拓扑边界态可以形成两种腔体。第一种是开放波导构型边界态沿着一条直线或折线传播两端用反射镜或者晶格失配区作为腔镜形成边界态激光第二种是闭合构型把拓扑畴的界面设计成一个环边界态绕成一圈形成一个不需要端镜的环形拓扑腔。第二种更容易实现连续单模激射因为能形成稳定驻波或行波模式不需要外部额外反馈。在设计时要重点看两个指标一是边界态在目标频率处是否有足够的局域场增强二是它和增益谱的匹配程度。如果边界态色散太陡即使频谱重叠模式体积也可能过大如果色散太平虽然态密度高但容易先激励起其他杂散模式。所以理想情况是让边界态尽量落在低损耗窗口同时和增益谱中心接近这样阈值会比较低。3. 结构与参数设计先把能算的算清楚3.1 选晶格构型从能带到拓扑畴的一步我做仿真时最常用的是三角晶格空气孔结构因为工艺上在半导体膜里刻圆孔或椭圆孔非常成熟。另一种是蜂窝晶格介质柱常见于后期想把增益介质放在柱子里、实现更精确控制的情况。两种结构各有取舍我总结了一张常用对比表。设计维度三角晶格空气孔蜂窝晶格介质柱工艺难度低孔阵刻蚀即可高柱阵悬空或部分刻蚀较难增益介质放置增益在膜内与孔壁区域重叠可以放在柱内模式重叠更集中谷带隙实现圆孔偏移、椭圆孔定向旋转、孔半径交替子晶格柱半径交替界面态宽度场集中在界面两侧 2-3 个晶胞场约束通常更强但损耗也高实验报道成熟度最常见文献例子多相对少更多用于理论验证我自己最后选了三角晶格椭圆孔旋转的方案。原因很简单电子束曝光写椭圆孔的保真度高而且椭圆旋转角的正负直接对应谷陈数的正负版图上只要把相邻区域的椭圆方向改成相反界面就自然形成了不用额外做复杂补偿。参数初值可以这样设工作波段在 1550 nm 附近时晶格常数 a 取 600 nm 到 650 nm膜厚 220 nm 左右空气孔椭圆长轴半径约为 0.25a短轴半径约为 0.18a旋转角度通常在 ±10 度到 ±30 度之间。这组参数下K 点带隙宽度大概有 20 到 40 nm 可调。需要强调这只是经验初值一定要在仿真里扫一遍因为每个人的材料和工艺条件不同。3.2 拓扑畴和域壁的构造要点两个拓扑相之间要有一个分界面这个界面就是边界态所在的“轨道”。在版图设计时通常会取一个掉头矩形或者平行四边形区域左半区所有晶胞的椭圆孔朝一个方向偏右半区朝另一个方向偏两侧的谷陈数相反。界面方向也很关键。理论分析表明沿 zigzag 方向也就是连接两个 K 点之间方向的界面更容易出现无间隙的边界态而且背散射抑制效果更好。如果沿 armchair 方向投影能带里边界态往往会与体模式重叠模式变为部分泄漏拓扑保护效果下降。我第一次做的时候因为没注意这一点把界面角度随便转了 30 度最后仿真发现模式损耗非常大后来才意识到是沿 armchair 方向投影成带边散射了。实际版图里我会把界面方向固定为 zigzag每个转角处再仔细检查跨界面的晶胞排列是否完整。即使做 60 度或 120 度折角只要每个分界面都保持同样的相对取向边界态就可以连续绕过去。转角处的晶格缺陷会带来一点散射但对谷锁定模式来说这种散射通常还在可控范围内。3.3 从被动结构到有源腔的建模套路我把这套流程在 COMSOL 和多物理场工具里反复跑了很多遍固定了一个工作流第一步无限周期能带计算。用平面波展开或者有限元的 Floquet 边界条件扫第一布里渊区的 K-Γ-M 路径得到体带能带图找到带隙范围。这一步主要确认带隙中心和宽度是否和设计目标接近。第二步超胞边界态计算。取一个包含两个拓扑畴的超胞保持平行于界面的方向为周期方向垂直界面方向用足够大的宽度并加吸收边界。算出来的本征模里能明显看到两条由界面态组成的色散带它们横跨体带隙而且场分布局域在界面附近。这一步要提取边界态群速度 v_g 和模式的有效折射率 n_eff。第三步有限结构透射仿真。用 FDTD 或者频域有限元做一截真实的拓扑波导加增益之前的被动结构加一个宽谱光源和一个监视器看透射谱里有没有对应边界态的共振峰。这个共振峰就提示这个结构可以支持边界态模式。第四步加增益材料模型。最简单的做法是在有源区叠加一个复折射率虚部为负的材料模拟净增益。更贴近实际的是 FDTD 里用集总速率方程模型给一个泵浦强度参数扫描不同泵浦强度看模式的增长是否越过阈值。在这一点上我建议不要直接用非常高增益来“硬激”否则会掩盖其他模式竞争的问题得不到可靠的阈值曲线。3.4 仿真参数里的几个隐蔽坑第一膜厚的影响。二维近似只适合做定性实际介质膜是有限厚度模场的垂直分布会改变有效折射率也会改变带隙位置。算完二维能带后一定要用三维仿真修正一次。第二空气孔椭圆旋转后网格各向异性会影响本征频率。如果使用三角形网格且网格不细旋转后的结构会带来虚假频移。最好在孔边界加一样细化的边界层并且收敛性测试。第三FDTD 里最大的坑是边界模式很容易被激励成“超模”也就是界面两侧的两个模式同时存在拍频干涉导致透射谱里出现额外抖动。处理办法是减小激发源的时间和空间相干性以及延长仿真时间做足够长的时域响应。4. 从仿真到实测一条相对完整的验证路径4.1 规划极限与测试矩阵光泵浦的拓扑激光器验证路径不能只做一个器件。我的建议是一次版图里放三组结构一组是只有直线边界态的波导腔用来测量最基本的激射特性一组是带一个 90 度或 120 度折角的边界态用来验证弯曲后能不能保持低损耗还有一组是真实环形拓扑腔用来测 Q 值和激射阈值。每组结构里再留几个尺寸渐变点比如晶格常数的偏移量从 -2% 到 2%椭圆旋转角从 ±15 度到 ±30 度。不要小看这张矩阵它能帮你快速筛出工艺窗。你没有必要把所有结构都当成激光器来测试可以先做无源透射对比把最优拓扑参数找出来再做有源器件这能省下很多测激光的时间。测试顺序也很重要。建议先做无源宽带透射或空间反射谱确认边界态透射峰位置和带隙范围再做上升 Twophoton 激发荧光或者边发射谱找到增益谱范围最后才做泵浦强度扫描看是否出现激射阈值。4.2 关键表征透射谱、L-L 曲线和线宽无源表征阶段我会用宽带超连续光源耦合进端面通过一根锥形光纤或者自由空间透镜收集透射谱。透过率峰的位置如果在模拟带隙内说明边界态确实存在。注意透射谱的峰宽不代表损耗需要用短脉冲或者自相关的方式测群延迟来推断群折射率。有源阶段最重要的是 L-L 曲线输出光强-泵浦功率关系。阈值以下时输出随泵浦缓慢增长阈值以上斜率会有一个明显突变。同时还要监测发射光谱的半高全宽通常阈值前后会从几纳米收窄到 0.2 纳米以内。如果你测到阈值以上还有多模纵模或者边模要检查模式竞争和空间烧孔问题。对于边界态激光而言空间分布的测量更关键。我会用近红外相机或者扫描近场显微镜把器件表面的发光分布拍下来。如果发光集中在界面处、并且在弯曲后依然沿着界面走这基本可以证明激射模式确实是拓扑边界态而不是某处的体带模式。4.3 拓扑保护效果的判据设置很多人问“怎么证明拓扑保护”我的做法是用同一个拓扑波导腔做几组对照在界面中间加一个刻蚀留下来的小突刺或者故意放一个尺寸微小的障碍物然后再测 L-L 曲线。如果普通波导激光器在这个突起下阈值会大幅上升而拓扑边界态激光器阈值变化很小那就可以说明它对这类扰动有鲁棒性。不过要注意这个鲁棒性不是无限大的。因为拓扑保护的是“传播方向”和“谷量子数”它不保护吸收损耗也不保护材料辐射效率。如果你加的障碍物已经大到改变了局部晶格的拓扑相甚至破坏了界面两侧的畴结构那边界态照样会消失。所以做对比实验时障碍物尺寸最好控制在 0.2 个晶格常数左右不要一上来就弄个高反射金属柱那样的结果没有指导意义。我还有一个很实用的判据把泵浦光斑位置从界面中心往左右偏移观察激射强度变化。如果泵浦光斑偏离开界面 5 微米后阈值急剧上升说明激射确实来自边界态而不是体区。这是最简单的空间選択性验证不需要复杂近场扫描。5. 实操中躲不开的坑与排查思路5.1 能带和边界态的识别错误仿真时最常见的误判是把投影能带里的连续体模式错当成边界态。判断边界态要同时看三个特征第一色散曲线落在体带隙内不与任何体模式重合第二场分布呈指数衰减地局域在界面两侧第三在超胞计算中边界态场在两个不同取向畴之间对称性会变化。如果只看到频带在带隙里但场分布是扩展到整个超胞的那大概率是有限尺寸下的体模。实验上也有类似问题。如果发光光斑很大很可能激射的是带边模式而不是边界态。一定要用远场/近场成像去确认空间轮廓不要只凭一条窄线宽就下结论。5.2 增益介质引起的频率偏移有源器件的模式频率和无源仿真经常对不上。原因是载流子注入会改变材料折射率增益谱峰值也会随载流子密度移动。无源仿真里算出来的边界态共振在 1555 nm加了增益材料后实际激射可能跑到 1563 nm 甚至更高。这不一定是设计错了而是折射率变化导致的频率红移。处理办法是在仿真初期就给有源区留出频率余量把带隙中心设计在增益峰值稍短波一点点的地方让注入后红移正好落在增益峰值附近。这个偏移量不好精确预测但一般留 10-20 nm 的窗口是合理的。如果做成可调泵浦功率你还会发现激射波长随泵浦强度有轻微漂移这通常是热和载流子效应的综合结果。5.3 工艺误差和刻蚀形貌空气孔刻蚀的垂直度对谷带隙影响非常大。如果孔侧壁倾斜等效介电常数梯度减弱谷陈数带来的带隙会越变越窄。椭圆孔的旋转角偏差比如本来 25 度刻成 22 度带隙中心也会移动。如果条件允许建议每次流片后在扫描电镜下认真测量孔的尺寸和旋向。不要认为相邻两个椭圆孔的旋向按照 GDS 文件就够了电子束曝光低剂量部分往往会出现“椭圆变圆”的趋势。我的经验是把椭圆孔长轴短轴比从 1.4:1 设计到 1.2:1制备后旋转角仍然保真但长轴短轴比低于 1.15 时谷带隙基本会被工艺涨落抹掉。5.4 泵浦热效应与模式竞争光泵浦时大部分泵浦光会转化成热特别是在高重复频率或连续泵浦下。器件温度升高会带来两件事一是增益谱红移二是折射率上升导致模式频率漂移。两者一旦错开激光性能会劣化甚至停振。最好使用低占空比的脉冲泵浦比如 1% 以下的占空比控制平均热功率。如果一定要连续泵浦结构需要做热衬底工艺或者把器件悬空以减少衬底散热路径。模式竞争是另一个容易翻车的地方。拓扑边界态虽然是带隙中的一支模式但体带边也有高态密度泵浦过强时有可能出现边模激射。要抑制边模可以把体能带边缘设计得离边界态稍远一点或者在界面上增加一段渐变损耗区让体模式更容易被漏掉。这个细节往往决定了单模成品率。5.5 仿真规模控制与收敛三维有限元做全尺寸拓扑腔非常耗资源。我通常先用二维模型做参数扫描筛选再对最有希望的几个结构做三维验证。三维验证时不必把整块衬底都仿真只取膜厚和上下包层包层用完美匹配层吸收。还要注意对称性结构有镜面对称时可以在频域里对称/反对称分开算计算量能降到 1/4。但前提是拓扑界面也是对称的如果加了折角或缺陷就不能用对称性简化。6. 后续能落地的方向与我的个人体会6.1 从单个激光器到拓扑激光器阵列单个拓扑激光器只是第一步。边界态天然沿着拓扑畴界面走因此可以把许多个界面并行排列形成阵列。如果相邻界面之间的距离足够远边界态互不泄漏如果设计得近它们之间会有耦合通过调整泵浦分布可以实现相位锁定。相位锁定的阵列激光器在相干合成、自由空间光通信上很有价值。谷光子晶体的平台优势在于模式锁定关系可以直接由拓扑结构决定不需要额外加复杂的表面光栅。在这个方向上我已经看到一些组在做一维阵列输出光场带涡旋相位。通过设计界面的周向闭合顺序边界态可以积累一个 Berry 相位让输出光携带轨道角动量。这个相位不是受工艺误差影响的而是拓扑决定的所以稳定性很好。6.2 传感、非线性、量子光学接口界面态对界面周围介质折射率变化很敏感很适合做微流控生物传感。拓扑激光器作为有源传感器有额外优势它本身是光源加传感腔合一可以监测阈值变化和线宽变化比单纯无源传感器多了模式频率和强度两个读出维度。非线性方面拓扑边界态把光场限制在很小的截面里同时具有慢光增强非线性系数会显著提升。再加上拓扑保护带来的模式稳定性可以去做可用的片上倍频、差频或者光频梳。量子光学方向把量子点嵌入增益区拓扑腔可以自发辐射增强同时通过谷锁定模式隔离环境噪声有可能做出更稳定的单光子源或纠缠源。6.3 一点很私人的总结这个方向最让我觉得有意思的不是“拓扑”听起来高级而是它逼着我把理论、工艺、电光测量全部打通。我踩过最大的坑是一开始太迷信拓扑保护疯狂在仿真里加各种缺陷来展示鲁棒性却忽略了增益谱和边界态频率的匹配。最后器件怎么泵浦都只能出很弱的自发辐射后来才意识到是材料增益峰值差了几十纳米。这次教训让我改了工作习惯先做增益谱测试再做拓扑结构设计很有用。如果你正准备上手这个方向我给的建议很简单先把一个被动无源拓扑波导的透射谱测透再去做激光。不要一上来就期待拓扑保护能补材料缺陷激光器的底线永远是损耗、增益和模式体积的平衡。拓扑能让你在结构上更自由但在物理上不免费。这句话是我这两年最真实的体会。