1、j sinisinj2)j 为相位厚度, 有Njd j cos j3)如公式( 2)所示,j 的表示为一个 2 2 的矩阵形式,其中每个矩阵元都没有任何实际物理意义,它只是一个计算结果,其推导过程将在第二部分给出。2.传输矩阵法在了解了传输矩阵的基础上,下面将介绍传输矩阵法的定义:传输矩阵法是将磁场在实空间的格点位置展开, 将麦克斯韦方程组化成传输矩阵形式,变成本征值求解问题。从其定义可以看出,传输矩阵法的实质就是将麦克斯韦方程转化为传输矩 阵,也就是传输矩阵法的建模过程, 具体如下: 利用麦克斯韦方程组求解两个紧 邻层面上的电场和磁场, 从而可以得到传输矩阵, 然后将单层结论推广到整个介 质
2、空间,由此即可计算出整个多层介质的透射系数和反射系数。传输矩阵法的特点:矩阵元少( 4 个),运算量小,速度快;关键:求解矩 阵元;适用介质:多层周期性交替排列介质。二、 传输矩阵的基础理论 薄膜光学理论1. 麦克斯韦方程组麦克斯韦方程组由四个场量: D、E、B、 H,两个源量: J、 以及反映它们 之间关系的方程组成。而且由媒质方程中的参数 、 、 反映介质对电磁场的 影响。方程组的实质是描述电磁场的传播, 即:一个变化的磁场引起邻近区域的 电场变化,而此电场的变化又引起邻近磁场的变化, 如此进行下去, 便可抽象出 电磁场的传播。如图 3 所示。HE图 3 电磁场传播的模拟图将媒质方程带入麦
3、克斯韦方程组,并对方程组求解可得以下两个重要结论:cv (4n)jk式(4)中,N即为介质的 光学导纳 ,单位为西门子。特别说明:光波段时,约等于 1,N 数值上等于折射率。自由空间导纳0.002653772E c22E2t(5)E1 E0H1 H0(E8)0H0,通电场和磁式(8)中,上标为 +的代表入射波, 得:- 表示反射波。又由导纳定义式( 4)可21(v62) 2t式( 5)为电场的波动方程,与经典波导方程( 6)相比可得常把光速 c 和电磁波在介质中速度之比定义为折射率,即得折射率公式:n (7)2. 边界条件及反射折射电磁波在介质交界处满足切向分量连续的边界条件。 垂直入射时,N
4、0(k E0 )N 0 ( k (9E0 )H1N1(k E1 )(10)将式( 9)、( 10)代入8)中,整理可得反射系数定义式:r E0 N0E0 N0N1场均与入射面垂直,则它们的切向分量既是本身。根据边界条件可得:(11)r 为反射系数, R 为反射率。 透射系数原理相同,在此不再推导上面讨论的是垂直入射的情况, 斜入射时情况类似, 只是用修正导纳 0 、 1 代替( 11)中的 N0、N1。其实,无论电磁波入射情况如何, 电磁波只有两种情况: 一种是电场 E 平行 入射面即 TM波( P分量),此时电场的切向分量 Etg Ecos ( 为入射角), 而磁场的切向分量是其本身,因此由
5、( 4)式可得:H Htg N(k Etg ) N(k Ecos ) (k E) (12)cos将(12)式与( 4)式对比可得到P 分量的修正导纳,同理可得TE波(S分量)的修正导纳:(13)N cos可得一般情况下的反射、透射系数表达式:01(14) 0 1介质的传光特性可以由反射、透射系数所表征, 而由以上讨论可知, 这两个参数与导纳紧紧联系。因此,求解介质的传光特性就可以转换为求解导纳问题, 这也是传输矩阵法所解决的核心问题之一。 其实,传输矩阵法就是通过求得介质 的导纳,从而得到介质的反射透射系数。3.传输矩阵这一部分将应用薄膜光学理论详细推导介质的传输矩阵, 以及如何求得介质 导纳
6、,根据第一部分传输矩阵的介绍可以知道,它其实是每层特征矩阵的乘积, 所以,这一部分的推导就从单层薄膜的特殊矩阵入手, 进而推广到整个介质空间 推导出介质的传输矩阵。面就详细介绍单层薄膜的特殊矩阵。电磁波通过厚度为 d1 的单层薄膜过程如图 4 所示。E11N0E11 N1E12E2N2图 4 电磁波通过单层薄膜图 5 单层薄膜等效为介质面的示意图薄膜是存在一定厚度的,电磁波从 E0透过薄膜变为 E2 的过程,与简单的穿 过介质面相比多了个 E1的中间变换,如果可以将 E0和 E2 通过导纳直接联系起来,那么薄膜就可以等效为一个介质面(如图 5 所示),前面所介绍的反射透射 公式便可用。因此,我
7、们第一步完成从薄膜到介质面的等效推导。 令薄膜导纳(介 质面 1 和介质面 2 的组合导纳)为 Y ,则可得到薄膜的透射反射系数:20(15)由式( 15)可知,求得 Y 便可求得 r、t。由导纳定义并对薄膜的第一介质面应用边界连续条件可得:(16)H0 Y(k E0 )kE0 k(17)11 H 111(kE11)图 4 中的 E11 、E11 表示刚刚穿过介质面一的瞬时状态。 E12 、 E12 表示即将穿过介质面二的瞬时状态。 这两个瞬时状态的唯一不同只是因为薄膜厚度引入的相位因子,即有:E11eN 1d1 cos 1(18)将式( 18)代入式( 17)中可得式( 19),并将其转为矩
8、阵形式( 20):(k(k E12)eiE12 ) 1eE12)ei11e i 1E12)(19)ei11ei1i11e1k(20)同理,薄膜的第二介质面有如下关系式:E 2 E12k E2 k(21)H2H12E2)212(k1 H 2111(22)(23)式(20)、(23)分别表示介质面一、二两侧空间电磁场之间的联系,若将 式(23)代入式( 20)中相乘,则所得到的结果就表示整个薄膜两侧空间电磁场 之间的联系,即:k E0e i 1 2k E21e i 1 1cos 1i 1 sin(24)从式( 24)中得到了第一层的特征矩阵:M1(25)H 0 Y(kE0)(26)H 2 2(k考
9、虑到导纳定义有如式( 26)的关系,则可对式(24)进一步化简:E0)1Yi sin(k E2)(27)B 令C为为膜系的特征方程,则有关系式:Ci1 sin(28)对比式( 24)等号左边的形式,由导纳定义可得整个单层薄膜的组合导纳:(29)从而由式( 15)可求得单层薄膜的反射、透射系数至此完成了第一步,即从薄膜到介质面的等效推导。将将单层得到的结论推广到整个介质空间可得:cos j sin j j(30)i j sin j cos jj Njdj cos j(31)ACBD(32)B CM(z)(33)(34)YCB(35)式( 30)为介质第 j 层的特征矩阵,需要注意的是特征矩阵的行
10、列式值为 1。 由式( 32)即可得到整个介质的传输矩阵。至此,完成了多层介质传输矩阵的建模过程。 值得一提的是,在讨论单层薄膜时,得到单层薄膜的反射率后,若对薄膜的 光学厚度 H(H=nd,n为薄膜折射率, d为薄膜实际厚度 )求导,可得如图 6 的结 果。从结果中我们可以看出, 在厚度为 时,反射率根据折射率的不同可达到最4图 6 反射率与光学厚度的关系大或最小值。三、 传输矩阵法的应用举例传输矩阵法的典型应用是对多层周期性交替排列介质的分析, 具有这样结构的器件实例有:光子晶体、光栅、量子阱结构、 DBR结构器件等。具体应用过程请参见文献 传输矩阵法分析一维光子晶体的传光特性 四、 小结
11、(1)传输矩阵法概念:将麦克斯韦方程组转换为传输矩阵的形式,应用传输矩 阵分析的计算方法。(2)传输矩阵:形式为每层特征矩阵的乘积。(3)典型应用:( 4)解决问题:传光特性( R、 T)、场强度( E、H)。注意:(3)、(4)共同决定传输矩阵法对所研究问题的适用性。( 5)重要结论:导纳 N、折射率定义 n ,光波段下,导纳无意义,它就是折射率。(6)传输矩阵的推导(薄膜光学理论)是繁琐的,但实际应用中可忽略推导, 直接应用结论式( 30)( 35)。(7)用传输矩阵法求解问题过程:1 )应用已有结论式( 30 )( 35)建立介质模型并求解:2)建立实际问题的模型。3)模型整合。(8) 额外的结论:薄膜厚度选为 的原因。参考文献1唐晋发,郑权 . 应用薄膜光学 . 上海科学技术出版社 .1984: 1-51.2贾习坤 . 基于传偷矩阵法对垂直腔半导体光放大器小信号增益特性的研究 . 西南交通大学. 2002: 6-133匡萃方,张志峰.传输矩阵法分析一维光子晶体的传光特性 . 激光杂志. 2003, (24)4: 38-39.
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1