集成光学2课件

上节课内容

in1n2n3xzdO

k1’i

i波导的模方程：传播常数：=n1k0sini有效折射率neff:neff=

/k0=n1sini导波存在条件：k2<<k1n2<neff<n1为上下两个波导界面的全反射相移1上节课内容in1n2n3xzdOk1’ii波导1第二章平面介质光波导平面介质光波导概述平板光波导的分析方法射线光学法

波动方程法条形光波导的分析方法马卡梯里法等效折射率法数值方法2第二章平面介质光波导平面介质光波导概述22麦克斯韦方程从麦克斯韦方程，建立光波在介质波导中的电磁场分布方程（波动方程），结合边界条件导出传播模式的特征方程，进而讨论介质波导中光传播的特性。时谐电磁场的麦克斯韦方程组

3麦克斯韦方程从麦克斯韦方程，建立光波在介质波导中的电磁场分布3波动方程将矢量各分量展开，得：

并且考虑到y方向是均匀的，即

设波沿着z方向传播，则沿z方向场的变化可用一个传输因子exp(iβz)来表示TE模(横电模)TM模(横磁模)

4波动方程将矢量各分量展开，得：并且考虑到y方向是均匀的，即4波动方程（TE模）

xn1n3n2yz包层薄膜层衬底5波动方程（TE模）xn1n3n2yz包层薄膜层衬底55波动方程的解上式为波动方程，也叫做Helmholtz方程，他的通解可表示为：其中，通常个成为横向波矢。a1,a2,j为待定系数。6波动方程的解上式为波动方程，也叫做Helmholtz方程，他6波动方程的解（场分布）根据物理意义可以预见在导波层内是驻波解，可用余弦函数表示，而在覆盖层、衬底层中是倏逝波，应是衰减解，用指数函数表示。导模存在条件：kx、a3、a2均应为实数，故须满足与射线法结果一致7波动方程的解（场分布）根据物理意义可以预见在导波层内是驻波解7边界条件边界条件为：边界处切向Ey分量连续，切向分量Hz也连续，由知连续(1)x=-a处，8边界条件边界条件为：边界处切向Ey分量连续，切向分量Hz也连8边界条件(2)x=a处，9边界条件(2)x=a处，99特征方程（本征值方程）TE模关于的函数10特征方程（本征值方程）TE模关于的函数1010TE模的特征方程：与射线法得到结果的一致：特征方程（本征值方程）关于的函数关于q的函数TM模的特征方程：11TE模的特征方程：与射线法得到结果的一致：特征方程（本征值方11特征方程（本征值方程）传播常数横向波矢12特征方程（本征值方程）传播常数横向波矢1212特征方程在这里，我们求解对称波导的情况，n1=n2(1)(2)考虑到：13特征方程在这里，我们求解对称波导的情况，n1=n2(1)(213图解法求解特征方程模式的阶数：mm越大，kx越大，b越小。

in1n2n3xzdO

k1’i

i14图解法求解特征方程模式的阶数：min1n2n3xzdO14图解法求解特征方程模式数量向下取整15图解法求解特征方程模式数量向下取整1515图解法求解特征方程0阶模总是存在1阶模存在条件:2阶模存在条件:16图解法求解特征方程0阶模总是存在1阶模存在条件:2阶模存在条16图解法求解特征方程单模条件17图解法求解特征方程单模条件1717截止波长对于对称平板波导，TE0和TM0的截止波长均为无限长截止波长18截止波长对于对称平板波导，TE0和TM0的截止波长均为无限长18特征方程与色散曲线特征方程中有4个参数（n1,n2,a,），改变任何一个结构参数都要对方程重新求解，不利于应用。为此作归一化处理。传播常数范围：归一化传播常数：波导参数V：19特征方程与色散曲线特征方程中有4个参数（n1,n2,a,）19（m=0,1,2,---）波导参数一旦确定，对应模的数量就确定；m=0模的传播常数最大，随着m的增大，传播常数减小；特征方程表示的是TE波（S波），习惯用模的阶数作为偏振光的下标，如TE0模，如TE1模等。特征方程20（m=0,1,2,---）波导参数一旦确定，对应模的数量20当给定V值后，波导中可能存在的导波模数量就确定了；b=0称为截止，从图中可见，TE1模的截止V值等于/2，如果V值小于/2，则只有一个模，称其为单模区；V值大于TE1模的截止V值，称其为多模区；TM波的导波参数V与TE波稍有不同，如果相对折射率不到1%，则同阶模的V值可认为相同。21当给定V值后，波导中可能存在的导波模数量就确定了；2121色散曲线均匀空间平面波的色散曲线：三层对称平板波导色散曲线：（一条直线）22色散曲线均匀空间平面波的色散曲线：三层对称平板波导色散曲线：22平板波导模式分布将特征方程的解，代入上式，并确定各个系数，求得Ey。而后根据右式，确定其余场分量。23平板波导模式分布将特征方程的解，代入上式，并确定各个系数，求23平板波导模式分布24平板波导模式分布2424场在覆盖层和衬底中是按指数函数衰减的，衰减的快慢分别由衰减系数a2和a3确定。平板波导模式分布a2和a3的值大，则场衰减快，穿透深度1/a2和1/a3就浅，说明场主要束缚在导波层中。反之，a2和a3的值小，则场衰减慢，穿透深度就深，说明波导束缚场的能力差。

a2和a3的大小与覆盖层、衬底的折射率有关，同时还与模序数m密切相关。由模式本征方程可以导出，m越大，则β越小，a2和a3也越小。这表明高阶模的电磁场可延伸到导波层外的距离较远。25场在覆盖层和衬底中是按指数函数衰减的，衰减的快慢分别由衰减系25求传播常数的顺序波导参数：a,n1,n2V参数色散曲线归一化传播常数b传播常数电场分布群速度色散波长：26求传播常数的顺序波导参数：a,n1,n2V参数色散曲线归一化26导模和辐射模qiCladdingCoreqiSubstrateCladdingCoreSubstrate辐射模的特点：传播常数连续。沿传播方向有损耗。与导模一起组成一个完备的正交函数集。导模的特点：传播常数取分立的值。理论上没有损耗。各个导模正交。27当s=m,为1，其余为0导模和辐射模qiCladdingCoreqiSubstrat27转移矩阵方法28麦克斯韦方程：均匀介质波动方程：转移矩阵方法28麦克斯韦方程：均匀介质波动方程：28各区域内都满足波动方程：区域与区域的边界满足边界条件：所有各层介质中的z方向传播满足:

E(x,y,z)=Et(x,y)*exp(ibz)转移矩阵方法29模式的基本性质：n1n2n3xzdO

电场E和磁场H的切向分量连续各区域内都满足波动方程：转移矩阵方法29模式的基本性质：n129转移矩阵方法30n1n2n3xzdO

波动方程的解(TE)：上式为一元二阶齐次微分方程，通解可写成：取特解：转移矩阵方法30n1n2n3xzdO波动方程的解30转移矩阵方法31n1n2n3xzdO

考虑n1区域：上式可写成矩阵形式：转移矩阵方法31n1n2n3xzdO考虑n1区域31转移矩阵方法32n1n2n3xzdO

M1是一个2*2的矩阵，称为n1区域的转移矩阵，它描述了n1区域中在两个分界面处的场的关系。转移矩阵方法32n1n2n3xzdOM1是一个232转移矩阵方法33……nN+1xzd1OnNnN+2

考虑多层平板波导：n1dN-1dN多层平板波导也可以用一个矩阵来表示最下面一个界面的场与最上面一个界面的场的关系。转移矩阵方法33……nN+1xzd1OnNnN+2考虑多33转移矩阵方法34……nN+1xzd1OnNnN+2

如何得到特征方程？n1dN-1dN最上方和最下方两个区域的场可表示为：将上式代入转移矩阵方程转移矩阵方法34……nN+1xzd1OnNnN+2如何得34转移矩阵方法35……nN+1xzd1OnNnN+2

得到：n1dN-1dN两边同时乘以上式为多层平板波导特证方程，用于求解b。转移矩阵方法35……nN+1xzd1OnNnN+2得到：35特征模的展开任意电场分布的光波入射如何转变成特征模？处理方法：将任意电场分布展开，分解成不同特征模的电磁场分布。数学上用正交函数展开，如傅立叶级数等，称之为特征模展开；各导波模以相应阶数模的传播常数传播；随着光的传播，不同模之间的相位差将发生变化，导致导波模叠加以后的电磁场分布也随着传播过程而变化，光束像蛇一样反复蠕动前进。特征模的展开任意电场分布的光波入射如何转变成特征模？36弯曲波导弯曲波导-光路变换的重要单元弯曲波导弯曲波导-光路变换的重要单元37弯曲波导模场分布弯曲波导的导模场分布偏向拐弯的外侧弯曲波导模场分布弯曲波导的导模场分布偏向拐弯的外侧38弯曲波导一般来说，

弯曲损耗↓：弯曲半径↑。波长↓。折射率差↑。角速度相同-----越往外的光场传播速度越快有一部分光场有可能跟不上-----形成弯曲损耗。弯曲波导一般来说，

弯曲损耗↓：角速度相同----39几种折射率差的SiO2波导*单模光纤：2a=8.9μm，Δ=0.27%(加了折射率匹配油)。*对于=1.55μm，90度弯曲波导的损耗小于0.1dB。40几种折射率差的SiO2波导*单模光纤：2a=8.9μm，Δ=40光波导中的各种损失在单模波导中导波模只有基模，其余展开分量全部转变成耦合损失，所以为减小耦合损失，应尽量使入射光束的形状与波导基模的形状相同。光波导中的各种损失在单模波导中导波模只有基模，其余展开分量全41渐变折射率波导在扩散性波导中，折射率分布多为渐变式：对称型渐变折射率波导42渐变折射率波导在扩散性波导中，折射率分布多为渐变式：对称型渐42渐变折射率波导中光的传播2axtpS0xzX光纤向折射率高的方向偏折渐变折射率波导中光的传播2axtpS0xzX光纤向折射率43渐变折射率波导的应用渐变折射率波导中的光纤传输呈周期性聚焦：若长度取1/4周期，则输出为平行光：中间可插入波片、偏振片、滤波片，方便进行光信号处理。P渐变折射率波导的应用渐变折射率波导中的光纤传输呈周期性聚焦：442.1平面介质光波导2.1.1平面介质光波导概述2.1.2平板光波导的分析方法射线光学法波动方程法2.1.3条形光波导的分析方法马卡梯里法等效折射率法数值方法452.1平面介质光波导2.1.1平面介质光波导概述4545条形光波导的分析方法光在平板波导中传播时，在无约束的方向上发散。为了避免这种情况，在集成光学中通常采用条形波导。和平板波导相比，条形波导的分析要复杂得多。通常采用近似的方法对此进行分析。

马卡梯里法

等效折射率法

数值方法

（有限差分，有限元）46条形光波导的分析方法光在平板波导中传播时，在无约束的方向上发46马卡梯里法(Marcaili)

由于导波模的大部分能量集中在波导芯层内传输，而在波层中的能量很少，图中阴影部分的能量就更少。马卡梯里(Marcaili)方法的近似处理是忽略图中的四个阴影区域，只考虑图中五个区域。矩形光波导的近似分析严格分析，必须将空间区域分为9块，每块中有6个场矢量的分量，同时必须在12个界面满足边界条件。Marcaili近似分析以Ex、Hy为主的模

以Ex、Hy为主的模47马卡梯里法(Marcaili)由于导波模的大部分能量集中在47从麦克斯韦尔方程出发，假设Hx分量为0，将所有场分量用Hy表示：模Hy满足波动方程：48从麦克斯韦尔方程出发，假设Hx分量为0，将所有场分量用Hy表48假设场分布E(x,y)可以表示成如下分离变量的形式将上式分裂成如下两个方程，即：分立变量法其中：对以上两个方程分开求解，然后在合并。49假设场分布E(x,y)可以表示成如下分离变量的形式将上式49

波动方程的解

其中为z方向传播常数kx,ky为x,y方向传播常数由分离变量法得到上述波动方程的解：50

50特征方程根据处和连续的条件，可以得到根据处和连续的条件，可以得到用数值方法可获得(kx,ky),之后再求出模斑分布可以用类似的方法得到51特征方程根据处和连续的条件，可以得到根据处和连续的条件51矩形波导的模场分布

(a)m=1,n=1(b)m=2,n=1;(c)m=1,n=2;Ex场分布：52矩形波导的模场分布(a)m=1,n522.1平面介质光波导2.1.1平面介质光波导概述2.1.2平板光波导的分析方法射线光学法波动方程法2.1.3条形光波导的分析方法

马卡梯里法

等效折射率法

数值方法532.1平面介质光波导2.1.1平面介质光波导概述5353等效折射率法(1)将三维波导分解成平板波导的组合(如图所示)；(2)对每个平板波导，沿着y方向根据平板波导本征方程求解neff(x)(图示neff2，neff3)；(3)根据步骤(2)中得到的neff2，neff3，得到沿着x方向的折射率分布(二维平板波导)；(4)对x方向的平板波导再解一次平板波导本征方程，可以得到传播常数。在整个求解过程中，需要注意的是，若求解TE模Ex，则在求neff(x)时，应取TE模，而在求时，应取TM模；求解TM模Ey，则在求neff(x)时，应取TM模，而在求时，应取TE模；54等效折射率法(1)将三维波导分解成平板波导的组合(如图所54等效折射率法等效折射率方法对于脊形波导特别适用：55等效折射率法等效折射率方法对于脊形波导特别适用：5555马卡梯里法VS等效折射率法56三维矩形光波导，n1>n2=n3=n4=n5：b-v关系：马卡梯里法和等效折射率法，均为近似方法，只在特定情况下比较准确。马卡梯里法VS等效折射率法56三维矩形光波导，n1>n2562.1平面介质光波导2.1.1平面介质光波导概述2.1.2平板光波导的分析方法射线光学法波动方程法2.1.3条形光波导的分析方法

马卡梯里法等效折射率法

数值方法572.1平面介质光波导2.1.1平面介质光波导概述5757数值解法解析解比较困难，近似会存在比较大的误差；数值方法被应用到光波导模式解中，他们可以应用于任何形式的波导结构和折射率分布。线方法(theMethodofLines,MoL)矩方法(theMethodofMoments,MoM)有限差分方法(theFiniteDifferenceMethod,FD)边界元方法(theBoundaryElementMethod,BEM)有限元方法(theFiniteElementMethod,FEM)58数值