耦合波导理论

第二章线性电光效应的耦合波理论2001年，She等人提出一种全新的理论，它从麦克斯韦方程出发，考虑二阶非线性极化强度（也就是只考虑线性电光效应），忽略其余高阶极化强度，推出关于线性电光效应的耦合波方程，得到在电场作用下的晶体中光的两个独立电场分量的解析解。这种方法，可运用于研究光在任意一个方向的电场作用下沿任意方向传播的各种线性电光效应的情况，并且不单可以用于研究光的振幅调制，也可以容易去解决光的相位调制问题。另外对于给定的一个晶体（点群），能根据需要利用该理论进行优化设计。这全新的耦合波理论相对折射率椭球理论来说，它的物理图象清晰，得到的结果是解析解，不用再作任何数学变换。我们不单可以方便地进行优化设计，而且也可用于电光调制器等电光器件性能的分析。它的出现拓展电光材料的选择范围和优化调制器的调制方式，从而引起了电光效应研究领域内新一轮的探索。2.1理论推导波在介质中传播时，能够通过介质内的非线性极化而相互作用将导致形形色色的非线性光学现象，如高次谐波、参量转换、受激散射等等。电光效应就是其中的一种非线性光学现象。电（波）与光（波）的互作用，实质上又可以看作是几个处于不同波段的电磁波在非线性介质中的波耦合过程，因此可以象非线性光学那样，通过求解耦合波方程来获得电光作用的有关知识。对于普克尔效应，是入射波为光3）+电波@m）产生一个输出光波（«+«m）的三波耦合过程。对于电光效应，它涉及到的是光与物质的相互作用，光是由麦克斯韦方程或场方程描述，物质体系是由光学布洛方程描述。于是我们采用类似非线性光学方法，首先给出相应的非线性极化强度，把电场所感生的附加极化矢量当成一个微扰量AP，再将它视为新的极化光源引入麦克斯韦波动方程，通过整理最后可得到相应的耦合波方程。线性电光效应耦合波理论就是以麦克斯韦波动方程为基础和出发点推导出来的。我们可以由麦克斯韦方程组和物质方程推导出：TOC\o"1-5"\h\z\o"CurrentDocument"V*V*E()+竺地!"^2P*(t)

c2dt20dt2根据矢量运算规则，VxVxE=V(V・E)-V2E(2-2)这样可得：82「fT-E(t)]d2Pnls(t)/、\o"CurrentDocument"V[V-E(t)]-V2E(t)+[8「刀=f―8^~)(2-3)3为介质的相对介电张量，r0为真空中的磁导率，c为真空中的光速，E(t)为介质中的总电场强度，Pnls(t)为只与电场强度E(t)有关的介质非线性极化强度，暂不考虑旋光效应。当光尚方向传播时，电场强度可分为平行和垂直于1的两个分量，因为此时光波理想化为单色平面电磁波，平行r的分量E〃(t)为零，所以我们只需保留E(f)垂直于传播方向r上的分量E±(t)。在没有自由电荷的均匀介质中和在P<<£oE的情况下，有V•E=0，这样方程(2-3)可变为:-V2-V2E上(t)+1fW(2-4)其中在单色波近似下，外加电场后晶体中总的电场强度可表示为：―1_一E(t)=E(0)+[_E(o)e-s+c.c.](2-5)2E(0)为外加直流电场或频率远小于o的低频电场；c.c.表示电场的复共轭部分；将(2-5)式代入(2-4)式的左边，可得：-V2E.(t)+82[m°)L=-2e&tV2E「o)-2竺周•E(o)]【e-血(2-6)由于电光晶体所产生的线性电光效应比其所产生的二次电光效应强得多，并且在实际应用中常利用立方晶系晶体或均质体来产生二次电光效应，因此由电光品体产生的二次电光效应就显得不重要了。在这里我们只考虑线性电光效应的贡献，而认为由于相位失配其它各二阶非线性效应以及更高阶非线性效应可以被忽略，所以在求解(2-4)式时，把非线性激励项作为一种微扰来处理。所以有：1「，、Pnls(t)=P(2)(t)=2P(2)(o)e-iot+c.c.(2-7)是方程(2-4)式的右边:TOC\o"1-5"\h\z\o"CurrentDocument"一^82P：LSG)=1目①2P(2)(①)e-s(2-8)\o"CurrentDocument"。8t22。^由(2-6)和(2-8)式，则式(2-4)可变为:\o"CurrentDocument"V2E(w)+竺B-E(w)]=―目。①2P(2)(①)(2-9)一般说来，在相位失配的情况下，频率为W的单色平面波在各向异性晶体中传播时没有倍频产生，电场可分为两个相互正交的偏振分量E«w)，E2(w)，设K、K2分别为E1(w)，E2(®)所对应的波矢，因此我们可定义:E(w)=E(w)+E(w)=E(r)e%.r+E(r)eKr(2-10)TOC\o"1-5"\h\z1212如果k=k，E(w)，E(w)分别表示光电场强度的两个相互垂直的分量；如果1212k卫k，E(w)，E(w)分别代表两个折射率不同，在晶体的传播中各自独立的电1212场强度。例如，在各项异性晶体中E1(w)，E2(w)分别表示。光和e光的电场强度。故(2-9)式可变为：w2一，、_,、一一£V2E(w)+一[U.E(w)]=—^w2P(2)(w)(2-11)j1c2j!10j1j=1,2线性电光效应可以与二阶电极化率张量X(2)联系起来，应只包含二阶非线性极化强度，忽略高阶的，其表达式为:(2-12)P(2)(w)=2s0X⑵(w，0):E(w)E(0)=2sX(2)(w，0):E(r)E(0)e.r+2sX⑵(w，0):E(r)E(0)e，k2.0102s0为真空中的介电常数，X⑵(w，0)为二阶极化率张量。(2-12)另外以方面又有:V2[E(r)eikr]=[—k2E(r)+2ikdEi(r)+d2Ei(r)询灯(2-13)drdr2在线性响应和介质无耗的情况下，偏振矢量和场振幅E(r)都是恒定的，与波通过介质时所运行的距离r无关。而在非线性响应的情况下，即使介质是无耗的，偏振矢量和场复振幅也都是r的函数。然而，因为非线性激励项是作为对线性效应

的一种微扰来处理的，因此我们可以认为电场复振幅因子E（r）是1的慢变化函数。于是考虑慢变振幅近似竺也«k竺史，由式（2-4），（2-13）和c2=-^，我们dr2dr80七可以得到：ik。"空也+ikev冬1dr2drdrdr(2-14)①2①2=——z(2)(①，0):E(r)E(0)e"-一*⑵(①，0):E(r)E(0)皿C21C22在忽略离散效应的情况下，我们记:E(r)=E(r)a(2-15)（2-16）E(r)=E(r)bE(0)=E0c我们分别用a、b来点乘方程(2-14(2-15)（2-16）dE1(r)-i①2a侦(2)(3,0):bcE(r)Eei^kr+i-^^a•*⑵"):acE(r)Edrkc220kc210dE2(r)-i①2b•*(2)(①，0):acE(r)Eei心+i①2b•*⑵(必，0):bcE(r)Edrk2c20k2c20其中Ak-k2—匕，对于无损耗介质的，*⑵（。0）是实数且满足全对称性排列，即有：a*⑵（必，0）:bc-b•*（2）（必，0）:ac（2-17）设n设n，，n2分别为两光波E，E2在介质中的折射率，有:（2-18）（2-19）（2-20）k-nk-n—k（2-18）（2-19）（2-20）TOC\o"1-5"\h\z1101c2202c把(2-17)和(2-18)式代入方程组(2-16)，则该方程组变为：dE(r)kk,八、，、i-i0a•*⑵(—，0):bcE(r)EeiAkr+10a•*(2)(—，0):acE(r)Edrn20n10\o"CurrentDocument"dE(r)kk、——-i—0a•*(2)(—，0):bcE(r)EeiAkr+1—b•*⑵(—，0):bcE(r)Edrn210n220又有电光张量兀元/二阶非线性极化率之间的关系:1，、*(2)(①，0)-—_(88)yjkl2jjkkjkl

,、,、,为有效电光系数，其表达式为:eff1eff2eff3r=£(e8)(aybc)eff1jjkkjjkikij^dr=£(88)(ayac)（2-21）eff2jjkkjjklklr=乙(e8)(bybc)eff3jjkkjjklklj,k,l2x(2)(—①，①，0)

jkl（2-21）8jj8kk又得dE(r)drdE(r)drdE(r)

——2—

drk—i—ar2neff1021k—i—arEE2n2eff1kEE(r)e&r—i—^rEE(r)1

ke-iAkr一i—卜rEE2（2-22）1(r)2(r)我们再令d1=d1=2TrfE01AkRd3=肃refE02d2=才rf211/k口d4=肃ref3E0.2（2-23）于是经过（2-22）经过整理最终可得耦合波方程组如下（2-24）匪1(r)=—idE(r)eiMr一idE(r)dr1221北Ir)=—idE(r)e-心一idE(r)dr31422.2线性电光效应耦合波理论的应用（2-24）在电光效应的应用问题上，主要存在两个问题：第一是在确定了直流电场的方向、入射光的传播方向以及电光张量元的情况下，如何得到出射光场的两个相互正交的偏振分量E1（®）和E2@）的表达式?这在上节的基本理论的介绍中已经详细的分析过了。第二是关于调制深度的问题，即我们如何在一个给定的直流电场下达到最大的调制深度？对于这问题，可以应用线性电光效应的耦合波方程组的普遍解来解决。它在电光调制方面的应用，不单可以进行纵向和横向调制，而

且还可以进行斜向调制。在解决问题的过程中可以看到，所有分析的关键是要算出有效电光系数ef1ef苛ef苛同之前的其他理论方法相比较，我们的方法简单明了。不竞可以分析单轴晶体，对双轴晶体的问题也可以容易解决。现在我们按k0和k0两种情)兄进行分析。2.2.1k0时的分析对于强度调制，在以下两种情兄中是k0的。一、光波沿单轴晶体的光轴传播时七。。带入式(2-21)、(2-23)，经过整理得:kn^1(222kd工［(42k)csin2七。。带入式(2-21)、(2-23)，经过整理得:kn^1(222kd工［(42k)csin2ccos2](2-25)1k)ccos22csin2]k6k0,且旅界位1，2,3)，由式(2-25)则可算出d10，此种情兄不适用于耦合波方程组，不需要做进一步讨论。但对于弘不全为零的晶体，我们可以适当调节角，使d10，而d2d40，此时可获得最大的调制深度。(i)6kck0，取/4时有d2d40，这种情兄可以在32,3m，6m2的k对称点群的晶系中实现.(ii6kck0，使k(2kik)cktan2(2-26)26kckk就有d2d40，就可以在42m，4，3，32，3m，6，6m2对称点群的晶体中实现。二、光波在无中心反演对称性的立方晶体中传播时43m和23对称点群的晶体就是属于没有中心反演对称性的立方晶体，它们的非零电光系数之间的关系为(2-27)同理，为了得到最大的振幅调制深度，必需找到d2-d4=0时的d1的最大值。运用方程(2-21),(2-23)和(2-27)我们可以得到d=(2-27)同理，为了得到最大的振幅调制深度，必需找到d2-d4=0时的d1的最大值。运用方程(2-21),(2-23)和(2-27)我们可以得到d=、”0凡r[c(ab+ab)+c(ab+ab)+c(ab+ab)]1263123322311331221(2-28)在a-a=1，b-b=1，c-c=1以及a-b=0的条件下，运用拉格朗日乘法我们可以在d2-d4=0的情况下得到气的最大值。如果直流电场的方向是任意的，就会有无穷多组(a,b,c)使得d1取到最大值；如果运用无限条件c•a=1或者c•b=1，就可得到最小的半波电压。2.2.2Rk。0时的分析当Ak>>d时，i即光波稍偏离光轴传播时，利用积分中值定理得电光效应耦合波方程组近似为竺#=-idE(r)dr21dE(r)—才——=-idE(r)(2-29)方程组的解为E(o)=E(0)ei(k1-d2)rE2(①)=E2(0)ei(k2-d4)r(2-30)这与r=0时所得结果相同。对于相位调制，我们可取E(0)=0，或E(0)=0,reff1则有E(①)=E(0)ei(k1-d2)rE?(①)=0(2-31)R(①)=0E2(®)=E2(0)ei(k2-d4)r(2-32)当Ak»d1时，振幅调制的结果是根据方程（2-29）计算的，由于方程不包含d,d3的项，所以对于电光强度调制，不需要考虑^,d3的值。相位调制的结果是根据（2-31）或（2-32）计算的，同样不包含d，d的项，所以同样不需要考虑d，131从方程（2-21）和（2-23），我们可以得到（2-33）（2-34）d=*0.0气Z（ysin2甲+ycos2甲一2ysin甲cos甲）c221k2（2-33）（2-34）d=*0气£（n4ycos20cos2甲+n4ycos20sin0+n4ysin20-42n01k02ke3k2kn2n2ysin20sin甲一n2n邛sin20cos甲+n2n2ycos20）c如果是双轴品体的情况，设n，n，n分别为晶体的三个主折射率，两个独立的尤yz偏振分量的单