摘要
本文研究了氮化硅微腔在2 μm波段的光学频率梳产生。通过几何设计对氮化硅波导进行色散调控,选取合适的总线波导尺寸,从热吸收理论出发,讨论氮化硅微腔不同调制频率下的热折射噪声。以非线性薛定谔方程为基础模型,研究不同色散作用下的腔体演变过程。数值结果表明,氮化硅在2 μm波段能够更清晰地观察系统的滞后状态转变,即系统向稳定域转变过程的弛豫振荡现象,同时在高阶色散的作用下,腔体能够更快地向稳态孤子转变,这为研究呼吸孤子提供了方案。
光学频率梳(Optical Frequency Comb, OFC)作为离散、等距的相关频率分量,其光谱类似于“梳尺”。在激光的锁模机制被发掘到频率梳集成化发展的70多年里,OFC已经成为研究热门领域之一。2 μm波段为大气低损耗波段,与传统信息带宽相比,采用C+L波段光传输的方式可以进一步拓宽通信带宽。利用空心芯光纤实现2 μm短距损耗通信已经被很好地验
本文针对Lugiato-Lefever Equation(LLE)进行数值建模,研究2 μm波段下Si3N4孤子光梳的锁定生成。通过有限元建模分析不同截面大小的波导参数,从而对色散进行精细调控,进一步讨论了微腔在2 μm波段下的热响应。分析泵浦光源与微腔演变的稳态条件,以2 μm为中心泵浦波长的情况下探讨单孤子光梳的锁定生成。模拟结果表明,在2 μm波段Si3N4微腔向稳态转变过程中,光谱存在周期性的展开与收缩,最后达到稳定孤子的产生条件,在考虑高阶色散的作用下,系统会减少振荡过程,更快地向稳态转变。
以可调谐窄线宽光源作为泵浦源驱动片上微环结构,该结构由直波导和高Q值微环谐振腔组成。经端面耦合后输入光通过耦合间隙进入高Q值谐振腔发生非线性作用,在色散与非线性以及增益与损耗相互作用下达到OFC产生的稳态条

图1 Si3N4波导模型以及频域场变化下有效模场面积
Fig. 1 Si3N4 waveguide model and effective mode field area under frequency domain field variations
有效模场面积:
. | (1) |
波导色散:
, | (2) |
其中,E为光场振幅,λ为波长,neff代表波导有效折射率,c为真空光速。
设计波导高度为0.76 μm,通过改变波导宽度,即宽度由1.1~2.0 μm线性递增,获得的色散曲线如

图2 波导色散特性曲线:(a) 波导高度设定为0.76 μm时,不同波导宽度条件下的波导色散特性曲线;(b) 波导宽度设定为2 μm时,不同波导高度条件下的波导色散特性曲线
Fig. 2 Waveguide dispersion curves: (a) waveguide dispersion curves for different waveguide widths with the waveguide height set to 0.76 μm; (b) waveguide dispersion curves for different waveguide heights with the waveguide width set to 2 μm
增加波导厚度或者波导宽度,都能获得更宽的负色散区域,有利于宽带OFC的产生。但是实际工艺加工过程中,厚氮化硅波导会带来更高的应力,因此实现更厚的波导对工艺提出了更高的要求,针对厚氮化硅波导的制备可以采用循环沉积-退
根据FEM对Si3N4微腔进一步进行研究(总线波导宽度为2.0 μm,波导高度为0.76 μm,半径为224 μm),以微腔共振模式频率ω0进行泰勒展开,方程如下所示:
, | (3) |
其中,中心泵浦波长ω0为2 μm,ωμ代表相对于ω0的参考频率,其中D1/2π为微腔的自由光谱范围(free spectral range, FSR),D2、D3分别代表了二阶色散、三阶色散值。求解得到微腔三阶色散以及集成色散如

图3 氮化硅微腔在中心泵浦频率为150 THz下的色散特性曲线
Fig. 3 Dispersion characteristic curve of a silicon nitride microcavity with a radius of 224 μm at a central pumping frequency of 150 THz
光脉冲在克尔非线性波导中传输遵循广义非线性薛定谔方程(nonlinear Schrödinger equation, NLSE)。通过对麦克斯韦方程组的推导解析,可以发现波导支持有限的模式进行传输,并且每个光模式对应特定的传播常数以及电磁场分)以及电场包络(A)以标量进行改写,以常数的形式对传播过程中色散、非线性、损耗等参数进行标
, | (4) |
其中,为波导单位长度损耗,
表示n阶色散。
NLSE为包含多项的高阶非线性微分方程,针对该方程模拟可以采用分步傅里叶变换(split-step Fourier)、四阶-龙格-库塔(fourth-order Runge-Kutta
, | (5) |
其中,TR是腔内渡越时间,Ain为外部泵浦场,t与测量环腔往返次数m相关,表示慢时间,κ表示环腔功率耦合系数,重新表示为环腔总损耗,为泵浦相对于环腔共振失谐位置,L是环腔往返长度。LLE可以用无量纲方程来建立与NLSE之间的联
. | (6) |
忽略非线性效应与色散等因素,冷腔共振角频率所对应的极值被规定在=0,在非线性作用下,可以观察到极值点,即共振峰值转移到
=X。取X=12,曲线整体呈现为克尔非线性倾斜,并且在黑色与红色曲线之间出现多值化稳态解,如

图4 取X=12时稳态解曲线,黑色曲线对应下分支,红色曲线对应上分支,虚线代表不稳定解
Fig. 4 Take the steady-state solution curve at X=12, the black curve corresponds to the lower branch and the red curve corresponds to the upper branch, dashed lines represent unstable solutions
在孤子存在区域进行参数锁定,从而快速获取稳态解是目前获得单孤子频梳方式之一。LLE很好地解释了OFC的产生以及相关的锁模性质,在LLE中引入拉曼效应从而精确地预测孤子光梳中的拉曼自频移现
连续的窄线宽光源驱动Si3N4微环谐振腔,光源波长由谐振腔共振波长的“蓝失谐”区域向“红失谐”区域过渡时,将会产生一系列丰富的微腔动力学变化。在环境达到稳态条件时,腔内将表现为多孤子态,即时域上表现为多个脉冲。选用2 μm×0.76 μm(波导宽度×波导高度)作为Si3N4总线波导尺寸。Si3N4具有正的热光系数(通常TO=2.45×1
, | (7) |
. | (8) |
取总体腔损耗κ/2π=4×1

图5 半径为224 μm的氮化硅微腔在中心泵浦频率为150 THz下的热响应曲线
Fig. 5 Thermal response curve of a silicon nitride microcavity with a radius of 224 μm at a central pumping frequency of 150 THz
为建立腔体与孤子之间的联系,通过改变泵浦的失谐量来获取腔内状态,通过合适的步长设计从而准确地获得腔体的反馈。方程不考虑拉曼、自陡峭、多光子吸收和自由载流子等效应,仅对色散、损耗、耦合以及非线性项进行研究。泵浦功率Pin参数为1.2 W,在经过共振波长后腔内的能量会产生一个骤降的过程,进一步控制失谐的位置以及速度可以在低功率的下分支获取孤子光频梳的稳态解。在考虑相关热参量的影响,利用LLE可以得到腔体的能量演变如

图6 (a) 二阶色散以及考虑二阶、三阶色散作用下的腔谱演变;(b) 不同色散作用下呼吸振荡边界向稳定孤子转变过程
Fig. 6 (a) The second-order dispersion and the evolution of the cavity spectrum considering the effects of the second- and third-order dispersions; (b) the transition process of respiratory oscillation boundary to stable soliton under different dispersions
在数值模拟过程,我们忽略了模式交叉带来的影响,在呼吸的失谐边带附近,光谱会周期性地展宽与收缩,如

图7 氮化硅在2 μm波段下呼吸孤子的失谐周期振荡以及最终孤子光梳的形成
Fig. 7 Detuned-period oscillations of breathing solitons in silicon nitride in the 2 μm band and the formation of final soliton optical combs
为了2 μm波段的光学频率梳产生,我们通过CMOS工艺对氮化硅微腔进行了制备与测试。利用780 nm厚度的氮化硅晶圆,进行器件的制备。首先利用匀胶机在表面涂覆一层负性光刻胶,考虑后续刻蚀工艺,刻蚀气体与光刻胶之间的选择比(1:1.5),因此负性光刻胶的胶厚为500 nm。通过电子束曝光(electron beam lithography, EBL)进行图案定制,经过显影后,我们采用CHF3与氧气(O2)混合气体进行刻蚀(比例为1∶5),为了减少氢键带来的吸收损耗,最后在退火炉中1 200 ℃的环境进行退火3 h得到最终的器件,器件以及波导截面SEM图如

图8 SEM图与透射谱:(a) 100 GHz氮化硅微腔SEM图;(b) 波导截面SEM图;(c) 氮化硅微腔2 μm透射谱
Fig. 8 SEM images and transmission spectrum: (a) SEM image of 100 GHz silicon nitride microcavity; (b) SEM image of waveguide cross section; (c) 2 μm transmission spectrum of silicon nitride microcavity
本文基于Si3N4参数模拟了2 μm波段光频梳产生过程并分析了孤子特性。在经过严格的色散调控后以2 μm×0.76 μm厚氮化硅波导作为2 μm孤子光频梳产生的模型,并且讨论了材料在该波段下的热吸收响应。在确定LLE相关参量后,根据该方程进行建模仿真,研究单孤子光梳的锁定生成方式,为实现2 μm为中心泵浦所产生的孤子光频梳奠定了基础。
模拟结果表明2 μm波段下,能更加容易地观察到系统的弛豫转变过程,这为研究呼吸孤子提供了一种新的思路,高阶色散的作用也可以进一步促进微腔混沌态的过渡。基于2 μm中心泵浦波长进行单孤子的生成可以进一步拓宽OFC的光学带宽,这可以为更大容量的通信提供解决方案。总之,本文从微腔色散调控以及热吸收响应出发,较为系统地研究了2 μm波段氮化硅微腔的动力学过程,同时对相关器件进行了制备与表征,可以提供相关的参考,以及推动进一步的研究工作。
References
Bradley T D, Jasion G T, Hayes J R, et al. Antiresonant hollow core fibre with 0.65 dB/km attenuation across the C and L telecommunication bands[C]//45th European Conference on Optical Communication (ECOC 2019). IET, 2019: 1-4. [百度学术]
Pfeifle J, Brasch V, Lauermann M, et al. Coherent terabit communications with microresonator Kerr frequency combs[J]. Nature Photonics, 2014, 8(5): 375-380. [百度学术]
Niu R, Li M, Wan S, et al. kHz-precision wavemeter based on reconfigurable microsoliton[J]. Nature Communications, 2023, 14(1): 169. [百度学术]
Reimer C, Kues M, Roztocki P, et al. Generation of multiphoton entangled quantum states by means of integrated frequency combs[J]. Science, 2016, 351(6278): 1176-1180. [百度学术]
Suh M G, Vahala K J. Soliton microcomb range measurement[J]. Science, 2018, 359(6378): 884-887. [百度学术]
Karpov M, Pfeiffer M H P, Liu J, et al. Photonic chip-based soliton frequency combs covering the biological imaging window[J]. Nature communications, 2018, 9(1): 1146. [百度学术]
Diddams S A, Vahala K, Udem T. Optical frequency combs: Coherently uniting the electromagnetic spectrum[J]. Science, 2020, 369(6501): eaay3676. [百度学术]
Wang Jun, Zhang Xiaoyan, et al. Research Progress of Two-Dimensional Nonlinear Optical Materials and Devices[J]. Chinese Journal of Lasers, 2017, 44(7): 0703004.王俊, 张晓艳, 张赛锋, 等. 二维非线性光学材料与器件研究进展[J]. 中国激光, 2017, 44(7): 0703004. [百度学术]
Wang X, Liu Y J, Zhang Z M, et al. Research progress in 2 μm waveband on-chip photonic integrated devices (Invited)[J]. Infrared and Laser Engineering, 2022, 51(3): 104-115. [百度学术]
王希,刘英杰,张子萌,等. 2 μm波段片上光子集成器件的研究进展(特邀)[J]. 红外与激光工程, 2022, 51(03): 104-115. [百度学术]
Gaeta A L, Lipson M, Kippenberg T J. Photonic-chip-based frequency combs[J]. Nature Photonics, 2019, 13(3): 158-169. [百度学术]
Herr T, Brasch V, Jost J D, et al. Temporal solitons in optical microresonators[J]. Nature Photonics, 2014, 8(2): 145-152. [百度学术]
Brasch V, Geiselmann M, Herr T, et al. Photonic chip–based optical frequency comb using soliton Cherenkov radiation[J]. Science, 2016, 351(6271): 357-360. [百度学术]
Guo H, Karpov M, Lucas E, et al. Universal dynamics and deterministic switching of dissipative Kerr solitons in optical microresonators[J]. Nature Physics, 2017, 13(1): 94-102. [百度学术]
Yao S, Bao C, Wang P, et al. Generation of stable and breathing flat-top solitons via Raman assisted four wave mixing in microresonators[J]. Physical Review A, 2020, 101(2): 023833. [百度学术]
Tikan A, Riemensberger J, Komagata K, et al. Emergent nonlinear phenomena in a driven dissipative photonic dimer[J]. Nature Physics, 2021, 17(5): 604-610. [百度学术]
Lu J, Puzyrev D N, Pankratov V V, et al. Two-colour dissipative solitons and breathers in microresonator second-harmonic generation[J]. Nature Communications, 2023, 14(1): 2798. [百度学术]
Kippenberg T J, Gaeta A L, Lipson M, et al. Dissipative Kerr solitons in optical microresonators[J]. Science, 2018, 361(6402): eaan8083. [百度学术]
Luke K, Okawachi Y, Lamont M R E, et al. Broadband mid-infrared frequency comb generation in a Si3N4 microresonator[J]. Optics Letters, 2015, 40(21): 4823-4826. [百度学术]
Pfeiffer M H P, Kordts A, Brasch V, et al. Photonic Damascene process for integrated high-Q microresonator based nonlinear photonics[J]. Optica, 2016, 3(1): 20-25. [百度学术]
Agrawal G P. Fiber-optic communication systems[M]. John Wiley & Sons, 2012. [百度学术]
Tzolov V P, Fontaine M, Godbout N, et al. Nonlinear self-phase-modulation effects: a vectorial first-order perturbation approach[J]. Optics Letters, 1995, 20(5): 456-458. [百度学术]
Hult J. A fourth-order Runge–Kutta in the interaction picture method for simulating supercontinuum generation in optical fibers[J]. Journal of Lightwave Technology, 2007, 25(12): 3770-3775. [百度学术]
Lugiato L A, Lefever R. Spatial dissipative structures in passive optical systems[J]. Physical Review Letters, 1987, 58(21): 2209. [百度学术]
Coen S, Randle H G, Sylvestre T, et al. Modeling of octave-spanning Kerr frequency combs using a generalized mean-field Lugiato–Lefever model[J]. Optics Letters, 2013, 38(1): 37-39. [百度学术]
Leo F, Coen S, Kockaert P, et al. Temporal cavity solitons in one-dimensional Kerr media as bits in an all-optical buffer[J]. Nature Photonics, 2010, 4(7): 471-476. [百度学术]
Pasquazi A, Peccianti M, Razzari L, et al. Micro-combs: A novel generation of optical sources[J]. Physics Reports, 2018, 729: 1-81. [百度学术]
Coen S, Erkintalo M. Universal scaling laws of Kerr frequency combs[J]. Optics Letters, 2013, 38(11): 1790-1792. [百度学术]
Karpov M, Guo H, Kordts A, et al. Raman self-frequency shift of dissipative Kerr solitons in an optical microresonator[J]. Physical Review Letters, 2016, 116(10): 103902. [百度学术]
Guo H, Lucas E, Pfeiffer M H P, et al. Intermode breather solitons in optical microresonators[J]. Physical Review X, 2017, 7(4): 041055. [百度学术]
Wang W, Lu Z, Zhang W, et al. Robust soliton crystals in a thermally controlled microresonator[J]. Optics Letters, 2018, 43(9): 2002-2005. [百度学术]
Xue X, Xuan Y, Liu Y, et al. Mode-locked dark pulse Kerr combs in normal-dispersion microresonators[J]. Nature Photonics, 2015, 9(9): 594-600. [百度学术]
Helgason Ó B, Arteaga-Sierra F R, Ye Z, et al. Dissipative solitons in photonic molecules[J]. Nature Photonics, 2021, 15(4): 305-310. [百度学术]
Arbabi A, Goddard L L. Measurements of the refractive indices and thermo-optic coefficients of Si 3 N 4 and SiO x using microring resonances[J]. Optics letters, 2013, 38(19): 3878-3881. [百度学术]
Huang G, Lucas E, Liu J, et al. Thermorefractive noise in silicon-nitride microresonators[J]. Physical Review A, 2019, 99(6): 061801. [百度学术]
Guo H, Lucas E, Pfeiffer M H P, et al. Intermode breather solitons in optical microresonators[J]. Physical Review X, 2017, 7(4): 041055. [百度学术]
Yi X, Yang Q F, Zhang X, et al. Single-mode dispersive waves and soliton microcomb dynamics[J]. Nature communications, 2017, 8(1): 14869. [百度学术]