|

引言
) \/ A4 D5 e. n% J) B1 ?6 n! g2 l硅基光电子技术已成为开发复杂光波导线路的关键技术,在微型化和集成方面具有显著优势。然而,硅绝缘体(SOI)波导中的高折射率对比使其对制造变化非常敏感,这会显著影响器件性能。宽度和厚度的尺寸偏差可能发生在一个芯片上的不同位置、晶圆上的不同芯片之间以及不同晶圆之间。这些变化在实现高产量方面带来了重大挑战,特别是对于干涉和谐振波导线路[1]。; g9 g, v @: _3 m, C* ` `- U
exxq1tr4mk2640591452.png
! O+ r' j# ?6 B; e+ ]. M4 g% l
6 ^" g3 l) K6 O; c18 V% H' i& B c5 S0 g1 j
波导参数提取简介7 q4 ~ _; A% o# O, Y/ d% |# o
传统的波导尺寸提取方法通常涉及多个器件,如低阶马赫-曾德尔干涉仪(MZI)或环形谐振器。这些方法存在局限性,如占用面积增加、需要多次测量、以及假设参数在器件间均匀性等。本文探讨了一种创新方法,使用单个高阶MZI准确提取SOI条形波导尺寸,特别关注在IMEC的iSiPP50G硅基光电子平台中实现的技术[1]。/ J9 s9 t0 p9 J( j& |5 I* i
$ W& v, R. L5 U) R
波导参数提取过程包括两个不同步骤:首先,从线路的测量谱中提取有效折射率(neff)和群折射率(ng);其次,从这些提取的指数确定宽度和厚度。我们将探讨新颖的映射算法和迭代方法如何提高提取精度同时简化过程。& I& Y) g$ ^' G) k, Y. F2 N
( d# A* z, \( o* s% v% \2* I) q- W/ M+ Q- O7 [
波导几何结构和制造变化
0 r2 f; A! }7 P7 z' Y- U首先了解本研究中使用的二氧化硅包覆波导的横截面模型。
6 I' _" Y: t) A, y' ^
4ks4kgpt3pb640591552.png
* u$ Q) `2 l6 c5 X7 _" C
图1:二氧化硅包覆波导几何横截面模型,显示具有宽度w、厚度t和侧壁角α的硅核心。
3 m/ [7 l7 I3 |0 H
3 | Y' ]4 O) K# L5 Z如图1所示,波导在核心底部有宽度(w)和厚度(t),通常呈等腰梯形,底角约为85°,而非完美矩形剖面。对制造的标准波导(设计宽度w = 450 nm和厚度t = 220 nm)的高精度测量显示,实际平均尺寸约为458 nm宽和212 nm厚,最大偏差约为宽度35 nm和厚度3 nm。
/ ?; O) s& V6 i' J+ z这些制造变化直接影响波导的光学特性。波导模式的有效和群折射率对这些几何参数高度敏感,因此准确表征这些参数对可靠的线路设计很重要。
; @; u# _7 A* u2 {- k! b
; Y# A1 l1 A6 X+ L* ?1 a3% ^7 t4 q# U7 V5 v
折射率数据库对模拟的影响
0 j$ V! S L4 ^硅折射率数据库的选择显著影响模拟结果。比较来自不同来源的数据,如H.H. Li(1980年)、D. Franta(2017年)和Lumerical数据库,可以发现明显差异。
. f ~1 W' D! d" `0 O+ y8 L8 C7 C
nt5aalqozlx640591652.png
, F8 _/ P& Y+ @图2:来自不同数据库的硅折射率,显示不同波长下的值变化。* f' r. F `4 q- Z
' X% N% n% P8 O, ]7 L0 {
图2说明了不同数据库之间折射率的变化。特别值得注意的是,在1.532μm波长处,Lumerical数据库中观察到的不连续性导致计算的群折射率跳变0.08。这种不连续性可能引入提取宽度的变化达40 nm,显著影响参数提取精度。对于本研究,选择了Franta数据库,原因是其较新(2017年)相比Li的1980年数据。
9 M$ m) o0 l! g: K4 R. B0 G6 `+ f! S1 @9 F8 ~6 ^: S$ }1 m7 N m
4
- Q- z7 x: Q' t; [波导有效折射率模拟
5 g6 Q2 {, B) x+ L) d. v波导的有效折射率(neff)随宽度和厚度而变化,如模拟所示。
. f: Y2 X9 S0 N L D! B o& [
2jl35f5hurl640591752.png
0 i7 {1 t3 a" F) D; |1 b% M/ X图3:具有不同宽度和厚度的波导的有效折射率,显示neff随两个参数增加而增加。. Y5 c7 M: e# f2 G, x
) I- H# N! [1 \; {# x1 m5 \
图3展示了宽度范围从400 nm到640 nm和厚度从200 nm到224 nm的波导的模拟结果。对于典型波导(w = 450 nm,t = 212 nm),有效折射率对宽度和厚度的敏感度计算为:
% u! ?. O* u3 D% l' U1 ?- v) z( K6 e( a. ~4 t( X& B/ x" N
?neff/?w = 0.0024 nm?1 ?neff/?t = 0.0037 nm?1: Q4 J0 h: B' `3 U' R
3 X, F# |4 c* j& | Z5 R& I$ B
这表明虽然neff对厚度变化比对宽度更敏感,但由于宽度的制造变化(35 nm)比厚度(3 nm)大,宽度偏差对neff总变化的贡献显著(88%)。
& n* M0 j# w# `% m& y6 H) l
0 x% r, Z% A2 h! @5
" ~; _5 r: U- V& r: N( p群折射率计算和敏感度分析7 T; I& W6 I/ _- N8 l8 A" \; [
波导的群折射率(ng)也随宽度和厚度变化,但敏感度模式不同。
& e0 i& g/ p' d/ |6 D% D
isive0o0zee640591853.png
% D- T2 c5 ^) D$ ]# N! [图4:具有不同宽度和厚度的波导的群折射率,显示ng如何随两个参数变化。1 q) z; g* r5 f" ^
% u! y2 G( v- b3 |. B* [) O, A
图4显示了群折射率如何随波导尺寸变化。对于典型波导,敏感度计算为:
+ C' i0 ~- w' z/ |$ d" R; ?
. J) p. S( ~1 `2 M- ]& u% L$ @8 |, R?ng/?w = -0.0020 nm?1 ?ng/?t = 0.0020 nm?1
% E0 ?6 M0 m0 ]% K5 V% n! _ @5 x; E( V. w$ y* X
有趣的是,这些敏感度的绝对值相当,这构成了提出的提取方法的基础。通过准确测量群折射率并使用这些敏感度之间的关系,波导宽度的不确定性可以从最大制造误差(超过30 nm)显著降低到约3 nm。9 y5 v2 ^. V6 z5 v D
5 [7 u: f3 w* ~4 D) y1 {; H68 r& T! i# _0 \8 G
参数提取的新型映射算法
+ X# _, ]+ L$ D- s# e$ R2 X研究人员开发了一种快速收敛迭代算法来解决传统映射方法的局限性。之前使用多项式模型的方法需要大量系数,仍然表现出显著的映射误差。
6 W0 B" m! Z: ?) U3 J0 K9 ]
pzpom5dqbiw640591953.png
1 z7 F! y5 Q8 b) f
图5:不同宽度和固定厚度212 nm的波导的Rng,显示敏感度比率如何在不同波导配置间变化。) g; X. }2 Z3 H7 |. {4 ^
0 c. v3 x4 z6 Q# @图5说明了因子Rng(群折射率对厚度敏感度与对宽度敏感度的比率)对不同波导配置的情况。这个因子有助于量化与初始计算宽度相关的不确定性。对于宽度范围为480 nm至580 nm的C波段氧化物包覆波导,宽度变化可以最小化到约0.9倍厚度差异。# l1 U+ G" o7 L! @: m' h
* T5 U8 J, } C( k: {" t u
该算法通过以下方式工作:使用ng对制造的波导宽度进行初步估计(假设预定厚度)将此宽度估计与neff一起计算厚度迭代细化这些值直至收敛* l) f; j" v0 W: U
[/ol]
2 W% | O. k5 d* ]. V H+ s这种方法在提高准确性的同时显著降低了计算复杂性。映射误差非常小:?wmapping = 0.006 nm和?tmapping = 0.002 nm,比传统映射方法有了显著改进。
* b) _3 W9 i* P" M
6 [4 B, E; d4 }% W- J( g7
) q8 i. O( X u- @- H/ z( i# M- \用于参数提取的MZI设计
' O0 K! v6 p7 F# {- n基于制造偏差分析及其对光学特性的影响,计算了用于参数提取的MZI所需的最小延迟长度。
, }& W5 [- P3 }; I
ntk1ylwfma4640592053.png
! K: s& a! C& g6 d3 @% w/ g" l# i
图6:设计的宽度扫描MZI,显示具有不同臂宽度的折叠MZI结构。1 O# \& Q' ?/ x& J
* t) v& n: Y3 k$ | X图6显示了实施的折叠MZI结构,具有两个作为分配器和组合器的方向耦合器(DC)。两个臂的宽度精心校准,范围从482 nm到546 nm,增量为4 nm。所有MZI的延迟长度设置为100 μm,足以在不引起干涉顺序歧义的情况下提取有效折射率和群折射率。
& u* {6 g. e5 D0 K X5 b; u- K* m( R* {! B9 O8 ^
83 \$ V7 t$ y' g5 U
测量设置和数据分析& s# F2 h6 |8 E- `5 ]
测量设置包括Luna光学矢量分析仪(LUNA OVA)、Polatis开关、六足对准系统和光纤阵列,使能够在单个芯片上对多个线路进行自动测量。
9 [( X P4 B* ^2 y3 O% R5 \, T4 V
row0vi2uhls640592153.png
; K& o4 I+ ]" G, X, K! d, s/ H图7:用于测量的设备,显示LUNA OVA、Polatis开关和光纤阵列设置。
0 I! E- _& F$ `! U, u* h( I/ Y, w6 D% m8 c5 Z
图7说明了测量设置,使能够在波导波长范围1529.30 nm至1571.27 nm内精确捕捉八个不同芯片上所有器件的传输谱。
! a" y% h( {- W k0 L/ E在分析测量谱时,理解MZI传输特性很重要。
3 ?- D' O3 j8 |# o0 T9 p' s
ewy201taknn640592253.png
( H& R/ s( B* L8 C, ^1 ^- R图8:一个设计的MZI的栏端口传输谱,显示原始测量和过滤信号。; F5 H; s8 j+ ], Y" r
4 T. U" F h. ?, d+ T
图8显示了一个MZI的栏端口传输谱,蓝线表示原始测量,红线显示应用低通滤波器后的信号,移除了光路中界面的反射。绿线表示上包络,计算用于消除光纤到芯片耦合效应。6 i* k6 x+ o& G
5g55504b0ob640592353.png
' q& ~1 ]! [5 h. v* e( f图9:一个MZI的归一化谱,显示用于参数提取的干净干涉模式。) l& t9 ^/ H4 t' z6 E
4 b( S* G& b# q/ A5 M4 o图9呈现了作为数学模型拟合目标函数的归一化谱。适当的归一化对准确的参数提取很重要。
0 s, ?% n2 B' [3 n
2 V8 K, i; M E* P4 k! U9
' `9 o9 Q; U6 J5 `' D: E提取结果和分析; z0 H b9 t* R/ i: _$ s4 B ?
将提取算法应用于测量谱后,确定了波导的有效和群折射率。3 B& v# r; n; A
uf0fxyt5vie640592453.png
2 A' G6 r9 Y# O' u& H
图10:所有线路中臂波导的提取有效和群折射率,显示这些参数如何随波导宽度变化。4 a, h7 c! Q6 @$ I- g
0 b9 A: r) d' j3 G+ w, l. h
图10显示了不同波导宽度的提取有效和群折射率。neff随设计的波导宽度增加而增加,而ng随着宽度扩大表现出下降趋势。平滑连续的变化表明在提取过程中没有关于干涉顺序的混淆。' z: K3 {2 w6 A0 o
# t) s" I j* C/ V3 Q* A
使用这些提取的光学特性,确定了波导的几何参数(宽度和厚度)。
5 _9 x- a) k5 P7 K& H9 W& p+ r
rw5csakir3q640592553.png
, {& c. f6 y9 M+ i* S6 }
图11:设计的MZI中臂波导的提取波导宽度偏差和厚度。
6 u: s2 l- ~1 l5 i; i! m% J# f* T4 H- V- |7 I
图11呈现了结果,显示宽度偏差(δw)在-0.90 nm和7.87 nm之间,厚度在210.95 nm和212.70 nm之间。计算的波导宽度最大制造偏差为8.77 nm,显著小于先前报告的35 nm测量值。提取的平均厚度为211.77 nm,与通过椭圆偏振仪测量的212 nm平均厚度紧密一致。
- N: U$ \7 @* V9 m9 ~( b& z: X' g0 ^) b a
10
! k0 s; l6 @ B; }0 e误差分析和讨论8 N( h6 T2 t3 @ F! v
几个因素导致总体提取误差:( L1 E- g. S; X
1. 数据库变化:使用Li的数据库而非Franta的结果,提取的波导宽度约窄4 nm,厚度约厚2.4 nm。& |% R) F. S' D
2. 归一化和拟合误差:由于数据归一化和拟合导致提取的neff和ng的标准偏差分别为1.46×10??和1.59×10??,转化为宽度0.058 nm和厚度0.030 nm的标准偏差。! E0 S2 B+ _! r# B/ h: a2 z" |# ?
3. 波导侧壁角变化:底角1°的变化对应宽度-3.610 nm和厚度-0.037 nm的变化。
7 E, t& u( ^: } M5 M0 ?4 _ M6 L! _0 t8 w9 L: { J' c1 r
总随机误差估计为:
. k4 U" k9 y6 m7 I6 q/ g宽度:±1.89 nm: e, ^& I/ y% {
厚度:±0.06 nm
( y0 i0 i0 d1 _1 a- H8 G
: J- n2 o D: S! G这种方法的一个重要限制是可以准确测量的波导厚度变化范围有限,不会在干涉顺序中产生歧义。对于100 μm延迟长度MZI,最大允许厚度偏差为3 nm。
' G5 p/ L- s" ?3 V1 G) ]% ]# o' C% H9 W: w
117 U4 |( Q. B- b$ J5 l
结论5 F3 S8 o# E9 d3 r
本文探讨了一种创新的映射算法,与高阶MZI一起操作,准确提取波导参数,包括宽度、厚度、有效折射率和群折射率。快速收敛迭代算法使能够使用100 μm延迟长度MZI从栏端口传输谱中提取参数,不会在干涉顺序上产生歧义。/ _3 Z& a, c5 f1 v' W: c' s
7 u% a# T% c/ Q; |2 m: c( m该方法在简化数据处理要求的同时提高了提取精度。提取的波导参数与高精度测量得到的参数紧密一致,证明了这种紧凑设计结合迭代算法对过程监控的有效性。
0 P: E- p$ z7 V* [( w ]8 G! B( p0 a3 c1 s) D
对误差来源的严格分析揭示,绝对提取误差主要源于材料折射率数据库中的不确定性,而提取宽度的相对误差主要归因于侧壁角度的不确定性。这些因素呈现了不能通过线路设计或参数提取算法改进轻易减轻的基本限制,突出了未来研究重点领域。. [7 Z! E$ B& e E
' s% n" }; w% ]2 r
参考文献
2 d( o+ V+ r( d+ @5 m4 X[1] Y. Liu, U. Khan, and W. Bogaerts, "Accurately extracting silicon waveguide dimensions from a single high-order Mach-Zehnder Interferometer," Optics Express, vol. 33, no. 6, pp. 13530-13546, Mar. 2025, doi: 10.1364/OE.558406.4 G2 v) h; ]8 T4 A5 K2 {+ `$ ~
END* H- }$ ^' b" [, x- X! n
" P2 O$ ^3 K0 A/ C软件申请我们欢迎化合物/硅基光电子芯片的研究人员和工程师申请体验免费版PIC Studio软件。无论是研究还是商业应用,PIC Studio都可提升您的工作效能。+ b( A& b8 n) @
点击左下角"阅读原文"马上申请* {) C& W8 l5 w( m, Z# N
( _% m: f) m! }欢迎转载
0 w7 Q, B" Q9 u2 N9 m- a; g
) Q. B# x$ G1 t; w' \! J( V) Z转载请注明出处,请勿修改内容和删除作者信息!
0 Y5 y1 c% C& D$ q$ o
" e1 C# j/ T7 B
+ Y% x, F4 v' S5 U$ i
: X2 @ V: i$ j
x2sr2ytidye640592653.gif
; J' @4 k* Q: t9 r7 E$ [& P2 `: n! J
关注我们* k* K4 x3 H8 o) i
$ I" L! d( o% C* L0 B- M; c V; j6 M3 W* \
datxdwbmr1k640592754.png
% k# H8 t, Z! a! Z% N9 Z5 @
| ; A6 ^) m3 Q2 |, h
wvaqljvlkmo640592854.png
8 \( \! K0 P$ h! s; u* J' F | 5 j: O4 f/ T8 ` [
bxzuqgby2ua640592954.png
" U; |# I2 Z# l! x0 E | ) A( d5 j5 |) Y7 {' `% _
% |! d( Q& |) X
# T) m% |0 s! E% n; K; z
0 C4 v. _ n% B
关于我们:
; z9 Y+ R7 \: }: _/ Y* J. \深圳逍遥科技有限公司(Latitude Design Automation Inc.)是一家专注于半导体芯片设计自动化(EDA)的高科技软件公司。我们自主开发特色工艺芯片设计和仿真软件,提供成熟的设计解决方案如PIC Studio、MEMS Studio和Meta Studio,分别针对光电芯片、微机电系统、超透镜的设计与仿真。我们提供特色工艺的半导体芯片集成电路版图、IP和PDK工程服务,广泛服务于光通讯、光计算、光量子通信和微纳光子器件领域的头部客户。逍遥科技与国内外晶圆代工厂及硅光/MEMS中试线合作,推动特色工艺半导体产业链发展,致力于为客户提供前沿技术与服务。
" U R5 w3 ^( w8 D/ {+ A: F$ Z5 Q2 b# T" Y+ _
http://www.latitudeda.com/ A3 R, }2 V) y- L! ]& ~
(点击上方名片关注我们,发现更多精彩内容) |
|