1. 引言
随着光通信技术的发展,光通信波长开始向近红外甚至中红外波段扩展,2 μm工作波长在通信、传感等领域受到了与日俱增的关注 [1] [2] 。2 μm波段应用的蓬勃发展带来了与日俱增的器件需求,要实现不同光模块之间的互连,功耗成为关键,面对未来光器件的需求,绿色低功耗特性成为了众多研究者的关注点。
硅基单片集成自从被提出就受到了举世瞩目的关注,目前硅基集成光子学器件主要利用典型的绝缘体上硅(SOI, Silicon-on-Insulator)实现,其波导层为硅材料,硅材料折射率比较高,与包层和衬底的折射率对比度大,有利于实现紧凑结构,但是波导侧壁粗糙给它带来了较大的辐射损耗,并且硅材料本身存在双光子吸收和自由载流子吸收效应,以上均表明硅材料不能适应目前低功耗器件的需求 [3] 。相比于硅波导,氮化硅材料与包层和衬底的折射率对比度适中,器件尺寸略大,侧壁粗糙引起的辐射损耗较小,并且氮化硅材料中无双光子吸收,因此氮化硅波导损耗更低,更适应目前硅基集成器件低功耗的发展趋势和需求 [4] 。
虽然氮化硅波导具有低功耗的优势,但是氮化硅波导器件尺寸较大,波导中光场强度较低,且氮化硅材料的非线性系数(n2 = 2.4 × 10−19 m2/W)比硅的非线性系数(n2 = 4.5 × 10−18 m2/W)低很多,导致其非线性效应弱,而非线性效应是实现波长转换、光频梳、光逻辑等的基础,因而氮化硅波导在非线性领域的应用受到了限制 [5] [6] 。
为了提高氮化硅波导的非线性效应,可利用氮化硅槽波导,相比于条波导和脊波导,氮化硅槽波导可将光场限制在低折射率的槽中,若在槽中填充具有高非线性系数的材料,可有效提高氮化硅槽波导的非线性效应。2021年,法国巴黎第十一大学的研究小组通过结构优化,在氮化硅槽波导中实现了30%以上的功率限制因子,但是未对其进行非线性效应的研究 [7] 。
本文采用非线性效应的全矢量模型和槽波导的有限元模场求解法,系统地研究了2 μm工作波长下氮化硅槽波导结构参数对其非线性系数和功率限制因子的影响。主要研究了氮化硅条宽度、槽宽度和波导高度对功率限制因子、有效模式面积和非线性系数的影响,然后,通过结构设计与优化,得到最大的非线性系数和最大的功率约束因子。
2. 方法
氮化硅槽波导结构如图1所示,衬底为二氧化硅材料,两个氮化硅条互相靠近从而构成槽波导,最后在整个波导上方包括槽中填充具有高非线性系数的材料。此处氮化硅槽波导取为对称结构,即槽两侧的氮化硅条的宽度相等,记为W,槽的宽度为Wslot,氮化硅槽波导的高度为H。在我们的模型中,包层材料选为DDMEBT,它在近红外波段具有较高的三阶非线性克尔系数和较低的材料吸收系数。氮化硅、二氧化硅和包层材料DDMEBT的折射率分别取为1.96、1.44和1.8 [8] [9] ,DDMEBT的非线性折射率系数在2 μm工作波长下为0.89 × 10−17 m2/W [10] ,比氮化硅和二氧化硅的非线性折射率系数分别高2和3个数量级。
![](//html.hanspub.org/file/5-2960278x7_hanspub.png?20230316174148717)
Figure 1. Diagram of the silicon nitride slot waveguide
图1. 氮化硅槽波导示意图
对于图1中所示的氮化硅槽波导结构而言,槽的方向与光传输平面垂直,被称为垂直型槽波导。由于包层和衬底的材料不同,因此,波导中存在的模式非严格的横电模式和横磁模式,而是准横电(Quasi Transverse-Electric, Quasi-TE)模式和准横磁(Quasi Transverse-Magnetic, Quasi-TM)模式,以下分别简称为TE和TM。其中,TM模式的电场强度沿图1中y方向,在垂直于氮化硅和上空气包层交界面的方向上是不连续的,所以TM模式的光场并不会集中在槽区域中。而TE模式的电场强度沿图1中x方向,在垂直于氮化硅和槽交界面的方向上是不连续的,但是其电位移矢量是连续的,由物质方程可得
(1)
其中E1和E2分别为交界面处氮化硅中的电场强度和空气槽中的电场强度,n1和n2分别为氮化硅材料和空气的折射率。由上式可知,由于氮化硅材料的折射率大于空气的折射率,所以槽中的电场强度要高于氮化硅中的电场强度,所以氮化硅槽波导可将光场主要集中在槽中,槽中的功率占比被称为槽波导的功率限制因子,一般用Γ来表示。
波导中三阶非线性相互作用的强度由非线性参数γ描述,γ取决于波导的几何形状和非线性材料的非线性折射率系数。为了准确地评价非线性,本文采用全矢量模型,将波导的非线性系数表示为 [11]
(2)
其中λ为工作波长,Aeff为有效模式面积,
为平均非线性折射率。Aeff和
的表达式为 [11]
(3)
(4)
其中E和H分别为传播模的电场和磁场,ez为指向正z方向的单位矢量,也是导模的传播方向,ɛ0和μ0分别为真空的介电常数和磁导率,n(x, y)和n2(x, y)分别为线性和非线性折射率。由式(3)可知,n2可以看作是非线性折射率在一个非均匀截面上对场分布加权的平均。
根据公式(2)和(4)可知,想要得到高非线性,需要减小有效模场面积或提高非线性材料中的功率,或两者同时进行。需要注意的是,要充分发挥槽波导的优势,必须保证槽区域的高功率限制。
3. 结果与讨论
3.1. 功率限制因子
为了对模式进行有效地限制和导引,波导的有效折射率应大于衬底包层/衬底的折射率,以减少光模式的泄露。氮化硅槽波导的初始参数分别为:波导高度为800 nm,氮化硅条宽度为300 nm,槽宽度为50 nm。其中每个结构参数的取值范围参考了现有加工工艺水平下氮化硅波导器件的尺寸。
围绕着氮化硅槽波导的初始结构参数,采用了控制变量法逐次研究了氮化硅槽波导的氮化硅条宽度、槽宽度和波导高度对槽中光功率限制因子的影响规律。
首先研究了氮化硅条宽度对氮化硅槽波导功率限制因子的影响,其变化情况如图2中黑色正方形标识线所示,此时槽宽度保持为50 nm,波导高度保持为800 nm,氮化硅条宽度从700 nm增加到1500 nm,步长为50 nm。可以看出,当控制其他结构参数不变时,随着氮化硅条宽度增加,氮化硅槽波导TE模式的功率限制因子几乎为线性减小,这是因为随着氮化硅条宽度的增加,波导的有效折射率明显增加,波导对光的约束能力增强,大部分的光被限制在氮化硅条中,因而泄漏到衬底和包层的功率以及限制在槽中的功率均减少。
![](//html.hanspub.org/file/5-2960278x14_hanspub.png?20230316174148717)
Figure 2. Variation of the power confinement factor with the slot width of the silicon nitride slot waveguide
图2. 氮化硅槽波导的功率限制因子随槽宽度的变化曲线
随后研究了槽宽度对槽中功率限制因子的影响,波导高度保持为800 nm,槽宽度从50 nm增加到250 nm,步长为50 nm,其变化情况如图2所示,此时有效折射率小于二氧化硅材料折射率的数据未显示。可以看出,当控制其他结构参数不变时,随着槽宽度从50 nm增加到250 nm时,槽内的功率限制因子增加,槽宽为250 nm时得到最大功率限制因子为7.9%。这是因为随着槽宽度的增加,槽区面积增加,因而限制在槽中的功率增加,槽内的功率限制因子增加。
最后研究了氮化硅槽波导高度对其功率限制因子的影响,其变化情况如图3所示,此时槽宽度保持为250 nm,氮化硅槽波导高度从800 nm增加到1500 nm,步长为100 nm。可见,当控制其他结构参数不变时,随着氮化硅波导高度的增加,氮化硅槽波导的功率限制因子增加。当波导高度从700 nm增加到1500 nm时,槽中的功率限制因子从7.9%增加到16.4%。这也是因为随着波导高度的增加,槽区面积增加,因而限制在槽中的功率增加,槽内的功率限制因子增加。
![](//html.hanspub.org/file/5-2960278x15_hanspub.png?20230316174148717)
Figure 3. Variation of the power confinement factor with the height of the silicon nitride slot waveguide
图3. 氮化硅槽波导的功率限制因子随波导高度的变化曲线
3.2. 有效模式面积
围绕着氮化硅槽波导的初始结构参数,采用了控制变量法逐次研究了氮化硅槽波导的氮化硅条宽度、槽宽度和波导高度对有效模场面积的影响规律。
首先研究了氮化硅条宽度对氮化硅槽波导有效模式面积的影响,其变化情况如图4中黑色正方形标识线所示,此时槽宽度保持为50 nm,波导高度保持为800 nm,氮化硅条宽度从700 nm增加到1500 nm,步长为50 nm。可以看出,当控制其他结构参数不变时,随着氮化硅条宽度增加,氮化硅槽波导的有效模式面积增加,这是因为随着氮化硅条宽度的增加,波导的有效折射率明显增加,波导对光的约束能力增强,大部分的光被限制在氮化硅条中,因为模式的有效模式面积增加。
随后研究了槽宽度对有效模式面积的影响,波导高度保持为800 nm,槽宽度从50 nm增加到250 nm,步长为50 nm,其变化情况如图4所示,可以看出,当控制其他结构参数不变时,随着槽宽度从50 nm增加到250 nm时,有效模式面积增加,槽宽为50 nm时可得到最小有效模式面积为2.96 μm2。这是因为光场主要现在槽区内,槽区的面积为波导高度乘以槽宽度,随着槽宽度的增加,有效模式面积必然增加。
最后研究了氮化硅槽波导高度对其有效模式面积的影响,其变化情况如图5所示,此时槽宽度保持为50 nm,氮化硅槽波导高度从800 nm增加到1500 nm,步长为100 nm。由图可见,当氮化硅波导高度从800 nm增加到1000 nm时,氮化硅槽波导的最小有效模式面积逐渐减小,均在氮化硅条宽度临界截止位置出现,当氮化硅波导高度从1000 nm增加到1500 nm时,最小有效模式面积几乎不变,对应的氮化硅条宽度也几乎未变,可得到最小有效模式面积为2.57 μm2,此时氮化硅条宽度为600 nm。这是因为当波导高度小于1000 nm时,随着波导高度的增加,氮化硅条宽度允许出现的最小值减小,因而最小有效模式面积减小,当波导高度大于1000 nm时,虽然波导高度继续增加,但是最小模式面积出现时氮化硅条宽度几乎未变,即模式分布未随波导高度而变化,因此最小有效模式面积几乎未变。
![](//html.hanspub.org/file/5-2960278x16_hanspub.png?20230316174148717)
Figure 4. Variation of the effective mode area with the slot width of the silicon nitride slot waveguide
图4. 氮化硅槽波导的有效模式面积随槽宽度的变化曲线
![](//html.hanspub.org/file/5-2960278x17_hanspub.png?20230316174148717)
Figure 5. Variation of the effective mode area with the height of the silicon nitride slot waveguide
图5. 氮化硅槽波导的有效模式面积随波导高度的变化曲线
3.3. 非线性克尔系数
在前两小节的研究基础上,围绕着氮化硅槽波导的初始结构参数,采用了控制变量法逐次研究了氮化硅槽波导的氮化硅条宽度、槽宽度和波导高度对非线性系数的影响规律。
首先研究了氮化硅条宽度对氮化硅槽波导非线性系数的影响,其变化情况如图6中黑色正方形标识线所示,此时槽宽度保持为50 nm,波导高度保持为800 nm,氮化硅条宽度从700 nm增加到1500 nm,步长为50 nm。可以看出,当控制其他结构参数不变时,随着氮化硅条宽度增加,氮化硅槽波导的非线性系数减小,这是因为随着氮化硅条宽度的增加,波导的有效模式面积增加,导致非线性系数减小。
随后研究了槽宽度对非线性系数的影响,波导高度保持为800 nm,槽宽度从50 nm增加到250 nm,步长为50 nm,其变化情况如图6所示,可以看出,当控制其他结构参数不变时,随着槽宽度从50 nm增加到250 nm时,非线性系数增加,但是槽宽为200 nm的氮化硅波导晚于槽宽为250 nm的氮化硅波导截止,因此可得到最大非线性系数,为3.05 /W/m。
![](//html.hanspub.org/file/5-2960278x18_hanspub.png?20230316174148717)
Figure 6. Variation of the nonlinear coefficient with the slot width of the silicon nitride slot waveguide
图6. 氮化硅槽波导的非线性系数随槽宽度的变化曲线
最后研究了氮化硅槽波导高度对其非线性系数的影响,其变化情况如图7所示,此时槽宽度保持为250 nm,氮化硅槽波导高度从800 nm增加到1500 nm,步长为100 nm。由图可见,虽然随着波导高度的增加,相同氮化硅条宽度对应的非线性系数减小,但是由于波导高度越高,波导的截止氮化硅条宽度越小,所以会得到更大的非线性系数,因此,当波导高度为1400 nm时,氮化硅条宽度为450 nm时,得到了最大非线性系数,为4.36/W/m。
![](//html.hanspub.org/file/5-2960278x19_hanspub.png?20230316174148717)
Figure 7. Variation of the nonlinear coefficient with the height of the silicon nitride slot waveguide
图7. 氮化硅槽波导的非线性系数随波导高度的变化曲线
4. 结论
本文主要研究了2 μm工作波长下氮化硅槽波导在不同结构参数下的功率限制因子,有效模式面积和非线性系数。仿真结果表明,氮化硅槽波导的功率限制因子随着氮化硅条宽的减小、随着槽宽的增加、随着波导高度的增加而增加。氮化硅槽波导的有效模式面积随着氮化硅条宽度的减小和槽宽的减小而减小,但是当波导高度高于1000 nm后,最小有效模式面积和所对应的氮化硅条宽度几乎不变。氮化硅槽波导的非线性系数随着氮化硅条宽度的减小、随着槽宽的增加、随着波导高度的减小而增加,但是由于高波导对应的截止氮化硅条宽度更小,所以波导高度更高时可以得到更大的非线性系数。本文在槽宽为250 nm,氮化硅条宽度为450 nm,波导高度为1500 nm时,可得到最大的功率限制因子和最小的非线性系数,分别为16.4%和4.36/W/m,此值与硅槽波导相比仍存在一定差距,但是本文的工作为2 μm工作波长下氮化硅槽波导的研究提供了理论基础,对基于氮化硅波导的非线性器件的研制与开发具有重要的参考价值。