 上传我的文档
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
毕业论文《铌酸锂晶体光学性能的仿真研究》
下载积分：1000
内容提示：毕业论文《铌酸锂晶体光学性能的仿真研究》
文档格式：DOC|
浏览次数：24|
上传日期： 15:50:21|
文档星级：
该用户还上传了这些文档
毕业论文《铌酸锂晶体光学性能的仿真研究》.DOC
官方公共微信铌酸锂晶体的研究与分析摘要本文基于新型线性电光效应耦合波理论，通过设定超晶格周期极化铌酸锂晶体倒格矢 参数，从而弥补双折射时 o 光和 e 光折射率不同造成的相位失配。计算有效电光系数，推 倒耦合波方程的解析解。并利用 matlab 进行线性仿真，研究温度，波长，外加电场和晶体 占空比变化时对于电光效应中的转换效率的影响
。 仿真的数值结果表明：随着温度与相位匹配时对应的温度的差值的增大，相位失配量 将增加，从而导致转换效率呈峰值逐渐降低的波动形式趋于零，当温度满足相位匹配时转 换效率最高；此外晶体极化周期数量的增加，将使得转换效率的波动更加剧烈，其值也降 低的更快，波动次数也将增多。改变波长的情况基本类似于温度，仅在波动细节上有细微 差距。 而电场对转换效率的影响则是成正比的线性关系。 此外我们所取的占空比约等于 0.25 和 0.75 时将可以使转换效率取到最大值。关键词 铌酸锂；电光效应；耦合波；转换效率I铌酸锂晶体的研究与分析AbstractBased on the new wave coupling theory of linear electro-optic effect. By setting the grating wave vector parameters of periodically poled LiNbO 3 crystal, we compensate for the phase-matched which is caused by different index of refraction of the o-ray and e-ray when birefringence happens.We can calculate the effective electro-optic coefficient of the system . Through the analytical solution of the wave coupling equations, use matlab to do linear simulation, and study on the influence of the conversion efficiency in the electro-optical effect when temperature, wavelength, electric field intensity and crystal duty cycle change. Numerical simulation results show that, with the increasing difference made by temperature which corresponds to the temperature and phase matching, the amount of phase mismatch will increase.it results in an increase of phase mismatch, thus causing the conversion efficiency to assume the fluctuating form tending to zero which the peak val when the temperature satisfies phase match, conversion efficiency is the highest. In addition when the number of the crystal polarization cycle increases,the fluctuation of the conversion efficiency will be more violent,the value will also reduce faster and the number of fluctuations will increase. The situation of the wave length is similar with the temperature.There is just only little discrepancy on the fluctuation in details of the conversion efficiency.But the influence of the electric field to transfer efficiency is the proportional linear relationship. In addition, when the duty cycle is equivalent to about 0.75 and 0.25,the conversion rate can be taken to the maximum.Key wordsLiNbO 3 ; electro- conversion efficiencyII铌酸锂晶体的研究与分析目录摘 要 ....................................................................................................................................I Abstract .................................................................................................................................. II 第 1 章 绪论 ......................................................................................................................... 1 1.1 引 言 ...................................................................................................................... 1 1.2 电光效应的理论发展 ............................................................................................ 1 1.3 研究方向和内容 .................................................................................................... 2 1.4 本章小结 ................................................................................................................ 3 第 2 章 LiNbO 3 晶体电光效应理论 ................................................................................... 4 2.1 电光效应基本椭球理论 ........................................................................................ 4 2.2 LiNbO 3 晶体的电光效应 ...................................................................................... 6 2.3 周期性极化 LiNbO 3 晶体（PPLN）的制备 ....................................................... 8 2.4 线性电光效应耦合波理论 .................................................................................... 9 2.5 本章小结 .............................................................................................................. 11 第 3 章 LiNbO 3 的晶体结构和性质 ................................................................................. 12 3.1 LiNbO 3 晶体结构 ................................................................................................ 12 3.2 LiNbO 3 晶体基本性质 ........................................................................................ 13 3.3 LiNbO 3 晶体特点 ................................................................................................ 13 3.4 PPLN 晶体的应用 ............................................................................................... 14 3.5 LiNbO 3 晶体 Sellmeier 方程 ............................................................................... 14 3.6 本章小结 .............................................................................................................. 15 第 4 章 系统结构和参数设定 ........................................................................................... 16 4.1 相关参数说明 ...................................................................................................... 16 4.2 PPLN 结构参数设定 ........................................................................................... 16 4.3 有效电光系数设定 .............................................................................................. 17 4.4 本章小结 .............................................................................................................. 18 第 5 章 线性仿真与讨论 ................................................................................................... 19 5.1 温度 T 的改变对转换效率? 的影响 .................................................................. 19 5.2 5.3 5.4 5.5 波长 ? 的改变对转换效率? 的影响 ................................................................... 22 外电场 E 的改变对转换效率? 的影响 .............................................................. 24 晶体占空比 D 的改变对转换效率? 的影响 ...................................................... 24 本章小结 .............................................................................................................. 25结 论 ................................................................................................................................. 26 参考文献 ............................................................................................................................... 27 致 谢 ................................................................................................................................. 28III铌酸锂晶体的研究与分析第 1 章 绪论1.1 引 言根据光的电磁理论我们知道，光波是一种电磁波。光波在晶体中的传播性质可以用一 个折射率椭球来描述，当晶体处在一个外加电场中时，晶体的折射率会发生变化，使传播 光波受到影响，折射率的改变正比与外加电场，这就是所谓的电光效应。电光效应可分为 两类，表现为介质折射率同外加电场成线性变化的电光效应称为线性电光效应（或称泡克 耳斯效应），表现为介质折射率同外加电场幅度的平方成比例的电光效应称为二次电光效 应（或称克尔效应）。线性电光效应只发生在没有反演中心的晶体中，由于它比二次电光 效应强很多，所以目前使用的电光调制器主要是基于线性电光效应。 而作为研究电光效应最基础的材料LiNbO 3 晶体，1949年首次发现铌酸锂具有铁电性。 从1965年Ballman等报道利用Czochralshi技术成功地生长LiNbO 3 单晶，以及1968年Larner 等报道了大直径，同成分的铌酸锂晶体生长以来，LiNbO 3 被广泛研究和应用，它在集成光 学和光波导应用中是一个重要的材料。有大的热电、压电、电光和光电常数等特性，使它 成为应用最广泛的电光材料之一。诸如应用于:声波转换器、声波迟缓器、声波过滤器、光 放大调制器、二次谐波器、光束转向器、相连接器、介电波导、存储元件、全息（光） 数 据处理装置等等。铌酸锂晶体已经成为人们研究非线性物理过程的模型晶体，很多分线性 物理过程，都是以它为研究平台展开的。1.2 电光效应的理论发展电光效应理论发展很早，早在十八世纪人们便发现了电光效应。二十世纪六十年代人 们已经利用电光效应进行调制和偏转，并且在光扫描，光存储，光显示等若干领域中有着 广泛的应用。现在电光效应的理论已经发展成熟，并且以此理论为基础得到各个方面的应 用。为了更好利用电光效应，人们不断的提出新理论并运用理论解决电光效应问题。新理 论方法的提出深化了我们对电光效应的认识，推动了电光效应应用的发展。自从1893年电 光效应被发现以来，人们从理论和实验中获得的几种主要的电光效应理论方法。 折射率椭球理论：由于光在晶体中的传播特性可以用折射率椭球完全的描述，所以人 们主要用电场对折射率椭球的影响来描述电光效应，建立在折射率椭球模型上的理论被称 为折射率椭球理论。折射率椭球模型简单直观易理解，所以长期以来人们都倾向于用它来 解决电光效应问题，但是在运用该理论来分析电光效应的过程中，存在着难以绕过的工作: 如何找到合适的坐标变换，从而使加电场后的折射率椭球方程主轴化。这个工作往往比较 复杂，有时甚至是不可能办到的，是折射率椭球理论运用中的难点。即使可以成功的得到 加电场后主轴化的折射率椭球方程，也仅知道加电场后的三个主折射率，难以获得此时晶 体中沿任意一个方向传播的偏振光的信息。所以人们只能用此理论研究电光效应的几种特1铌酸锂晶体的研究与分析殊情况，这就限制了电光效应在实际中的应用。所以，为了突破折射率椭球理论的局限性， 人们要寻找更有效的解决电光效应问题的理论。 特殊耦合波理论：特殊耦合波理论是针对电光效应一些特殊情况的研究，该理论值得 借鉴的地方在于它从麦克斯韦方程和晶体的电光效应出发，导出了入射光沿光轴方向传播 时的耦合波方程组，给出了单轴晶体中两偏振光（o光和e光）的解析解。但是作者给出的 是特殊情况下的耦合波方程组，导致文献给出的最终结果的实用价值有限。但它提出的这 种新的想法给了人们一些重要启示。将电场所感生的附加极化矢量视为一个微扰量，再将 这个微扰量当作新的极化波源引入麦克斯韦方程组中，建立起耦合波方程，通过求解方程 给出电光效应的衍射效率公式。它提出了一个很好的想法，但可惜它不能用来研究入射光 沿任意一个方向入射时的情况，而且还受到入射光方向和初始值等因素的限制，所以很难 用于电光调制器性能的优化。 平面波本征方程的微扰理论：优于以上介绍的两种理论，平面波本征方程的微扰理论 可以给出任意传播方向上的两偏振模式的折射率的改变量。由于在此理论中电光效应表示 的是微扰电场引起的一阶变化，所以在电磁场的波长达到电光晶体尺寸的数量级这个条件 下，可将微扰理论加入到本征矢量方程中来研究电光效应。该理论从电磁场的波动方程出 发，把晶体（包括各向同性的晶体、单轴晶体和双轴晶体中的电光效应当成微扰来处理， 得出了对应的微扰情况下的本征方程。于是，通过解出对应的本征值和本征矢量，可最终 得到任意方向的电场作用下，沿任意方向传播的光波的两种偏振模式的折射率改变量。虽 然这套理论在研究电光效应上有着非常大的进步，但是它无法给出这两个偏振模式在出射 面的场强表达式，而且在使用上也受到电磁场波长的限定，所以不能彻底克服折射率椭球 理论的局限性。 线性电光效应的耦合波理论：从折射率椭球理论到平面波本征方程的微扰理论，前面 所提到的这几套分析电光效应的理论都存在些不足和局限，对晶体上的外电场方向、对入 射光的偏振态和传播方向、对所使用的电光晶体的对称点群等方面，都有一定的限制。而 She等人所提出的线性电光效应的耦合波理论就可以很好地满足以上的要求，该理论从麦 克斯韦方程出发， 考虑到介质的二阶非线性光学效应， 建立了线性电光效应的耦合波理论， 给出了耦合波方程组及其普遍解。此解可以用来描述，在任意方向的外加电场的作用下， 任意偏振态的入射光在任意点群的电光晶体中沿任意方向传播时的情况。我们可以用这套 理论来研究电光调制器的温度特性，以及进行包括降低半波电压、提高消光比、提高调制 度等的调制器优化。本文我们就是从此出发讨论电光效应中转换效率等问题。1.3 研究方向和内容线性电光效应是电光调制器的物理基础。以折射率椭球理论为代表的传统的线性电光 效应理论各有所长，但是在使用时受到诸多限制，我们需要一种更方便的可用来解决线性 电光效应问题的理论。She等人提出的线性电光效应的耦合波理论从麦克斯韦方程组出发， 给出了偏振态不受限的光波在任意方向的外加电场作用下，在任意点群的电光晶体中沿任 意方向传播时出射光光强的表达式。我们的工作内容就是以该理论为基础。 我们由线性电光效应耦合波理论入手，在选定波长和温度的条件下，通过设定 PPLN2铌酸锂晶体的研究与分析晶体倒格矢（即极化周期）参数来弥补双折射情况下产生的 o 光和 e 光的相位失配量，从 而达到相位匹配进行电光调制。设定入射光线和加电场的方向，计算出此时系统有效电光 系数，解析耦合波理论中的微分方程。利用 matlab 进行线性仿真，研究温度，波长，外加 电场强度和晶体占空比变化时对于电光效应中的转换效率的影响。1.4 本章小结本章通过介绍电光效应理论的发展历程，确定了以线性电光效应耦合波理论作为理论 基础的必然性，为后文详细研究基于 LiNbO 3 光波导的电光效应中温度和外电场对于转换 效率的影响奠定基础。3铌酸锂晶体的研究与分析第 2 章 LiNbO 晶体电光效应理论32.1 电光效应基本椭球理论光在晶体中传播时，折射率随传播方向和偏振而异。在绝大多数晶体中，光的各向异 性性质是自然产生的，由晶体内部结构确定。不过人们也发现，通过各种物理效应，这种 特性可以由外部感生出来，电光效应便是其中一种。电光效应是指在直流电场（或低频电 场）的作用下引起材料折射率明显变化的一种现象。也就是说外加电场改变了介质的光学 性质。在某些材料中折射率的变化与所加电场的强度成线性关系，即线性电光效应，亦称 普克尔（Pockels）效应。线性电光效应可认为是入射光场与直流电场混合作用在物质中产 生的二阶非线性极化，由于线性电光效应是用二阶非线性极化率描写的，因此它只能在具 有空间非对称的晶体中发生。在有空间中心对称的材料中，比如液体或玻璃，折射率的变 化与所加电场的平方成正比，这就是二次效应或称克尔（Kerr）电光效应。与线性电光效 应类似，它可用三阶非线性极化来描写。除此之外，还有更高次的电光效应。但一般情况 下，高阶效应要比一次效应弱的多，所以在铌酸锂晶体中，我们只需考虑普克尔一次电光 效应。 光在各向异性介质中的传播特性可以通过求解麦克斯维方程并考虑到极化的各向异 性得出，不过数学过程相当繁复。人们发现，如果用几何图形来表示传播规律则显得十分 方便。为此人们引入了光率体，光率体又称为折射率椭球。电场的作用使晶体折射率椭球 主轴的方向和大小发生了变化。 在各向异性光学晶体中，光电场的电位移矢量D和电场强度E之间的关系写成分量式：D i ? ? 0 ? ? ij Ej j（2-1）或用下式表示：?D X ? D ? Y ?DZ ? ? ? ? XX 　 ? XY 　 ? XZ ? ? ? ? 0 ? YX 　 ? YY 　 ? YZ 　 ? ? ? ? ? ZX 　 ? ZY 　 ? ZZ ? ? ??E X ?? E ?? Y ??E Z ?? ? ? ? ? ?（2-2）无光学吸收损耗晶体的介电张量 ? ij 是一个对称矩阵，只有六个独立的张量元，即 ? xx , ? yy , ? zz , ? xy ? ? yx , ? zx ? ? xz , ? zy ? ? yz 。数学上一个对称矩阵可通过正交变换实现对角化。 物理上表示存在一个新坐标（XYZ），通过（XYZ）坐标系到（XYZ）坐标系的变换 使得（2-2）式具有简明的形式：?D X ? D ? Y ?DZ ? 0 0 ? ? ? XX 　　　　　 ? ? ? ? 0 0 　　 ? YY 　 0 　　 ? ? ? ? 0 　　　 　　 ? ZZ 0 ? ? ??E X ?? E ?? Y ??E Z ?? ? ? ? ? ?（2-3）这一新的坐标系就是晶体折射率主轴系统，晶体的介电张量在该坐标中是一对角矩 阵。 晶体中光电场的能量密度：4铌酸锂晶体的研究与分析U ? 1 2 D?E ? 1 2? 0 ? ? ij E i Eijj（2-4）在上述的主轴系统中，能量密度可写成U ? 1 2? 0 ? [ DXXX 2?DY2? YY?DZ2? ZZ]（2-5）（2-5）式表明，在 D 空间中光电场的等能面是一个椭球面，如图（2-1）。如设X ? ( 1 2 ? 0U1) 2 D X ,Y ? (1 2 ? 0U1) 2 DY , Z ? (1 2 ? 0U1)2 DZ（2-6）图 2-1 晶体折射率椭球则（2-4）式变成X2? XX?Y2? YY?Z2? ZZ?1（2-7）在（X，Y，Z）坐标系中，由（2-6）式决定的曲面为折射率椭球面。在这一主轴坐标 系中折射率椭球方程取最简洁的形式。 光在各向异性晶体中的传播特性可以用折射率折射率椭球来描述：过原点作一与晶体 内任意方向传播光波波矢垂直的平面，该平面与折射率椭球相交的截面是一个椭圆，椭圆 的长短轴分别为光波在晶体内该方向传播时的两个折射率，长短轴的方向为D矢量的偏振 方向。当晶体外加电压时，由于电光效应折射率椭球发生变化。这时椭球方程应取普遍的 形式( 1 n2)1 X2?(1 n2)2Y2?(1 n2)3 Z2? 2(1 n2) 4 YZ2? 2(1 n2) 5 ZX2? 2(1 n2) 6 YX2?1（2-8）将（2.8）式与（2.7）式比较，可知当没有外电场时，( 1 n2)1 ?1 n2 X,(1 n2)1 ?1 n2 Y,(1 n2)1 ?1 n2 Z,(1 n2)4 ? (1 n2)5 ? (1 n2)6 ? 0（2-9）当晶体加上外场时，则 (1 n2) i 量的变化为5铌酸锂晶体的研究与分析?( 1 n2)i ???j ijEj（2-10）采用矩阵的写法，则有? ? 11 　　 ? ? 　　 ? 21 ? ? 31 　　 1 (? 2 ) i ? ? n ? ? 41 　　 ? ? 　　 51 ? ? ? 61 　　 ?? 12 　　 ? 13 ? ? 22 　　 ? 23 ? 32 　　 ? 42 　　 ? 52 　　 ? 62 　　? ??E X ? 33 ? ? ? EY ? 43 ? ? ?E Z ? ? 53 ? ? ? 63 　 ? ?? ? ? ? ?（2-11）式中 ? ij 为线性电光系数， 它给出了 (1 n2) i 随所加电场强度 Ej增加时的变化。E X , E Y , E Z是外加电场在主轴坐标系中的三个分量。2.2 LiNbO 晶体的电光效应3LiNbO 3 晶体为单轴的铁电晶体，在没有外加电场时其标准的折射率椭球方程为X no2 2?Y2 2?Z2 2?1（2-12）none其中Z为光轴，结晶轴XYZ构成折射率主轴坐标系。LiNbO 3 晶体的点群对称群为3m， 其电光张量具有如下形式：? 0 　　　 ? ? 22 　 ? 13 ? ? ? 0 　　　 ? 22 　 ? 13 ? ? ? 0 　　　　0 　 ? 33 ? ?i ? ? ? 0 ? ? 0 　　　 ? 51 　 ? ? 　　　0 　 0 ? 51 ? ? 0 ? ? ? ? 22 　 0 　 　　　 ? ?（2-13）外加电场时，由于电光效应使LiNbO 3 折射率椭球发生的改变由下式给出：? ? 0 　　　 ? 22 　　 ? 13 ? ? ? 0 　　 ? 22 　　 ? 13 ? ??E ? x ? 0 　　　 0 　　 ? 33 ? ? ? 1 （? 2 ） ? ? ??E y ? i n ? 0 　　 ? 51 　　　 0 ? ? ? ? ? 　 　　　　 0 ? ? E z ? 0 51 ? ? ? ? 22 　　0 　 　 0 　 ? ? ? ?（2-14）其中 ? 13 ? 9 . 6 , ? 22 ? 6 . 8 , ? 33 ? 30 . 9 , ? 42 ? 32 . 6 （单位 10 ? 12 m/V）为LiNbO 3 晶体的电光系 数。6铌酸锂晶体的研究与分析晶体的折射率椭球方程则变为以下形式：( 1 n2 o??22E y ? ? 13 E z ) X2?(1 n2 o??22E y ? ? 13 E z )Y2?(1 n2 e? ? 33 E z ) Z2? 2 ? 51 E y YZ ?2 ? 51 E x XZ ? 2 ?22E x XY ? 1（2-15）我们考虑仅施加 Y 向电场时晶体折射率椭球的变化，即在式（2-15）中 E x ? E z ? 0 ， 代入式（2-15）得到( 1 n2 o??22Ey )X2?(1 n2 o??22E y )Y2?1 ne2Z2? 2 ? 51 E y YZ ? 1（2-16）在式（2-16）中仅存在 YZ 交叉项，做如下变换：?X ? X ' ? ' ' ? Y ? Y cos ? ? Z sin ? ? ' ' ? Z ? Y sin ? ? Z cos ?（2-17）将式（2-17）带入式（2-16），令交叉项为 0，则得到新的主轴坐标系下的方程为：( 1 n2 o??22Ey )X'2?(1 n2 o??22E y ? ? 51 E y tan ? )Y'2?(1 n2 e? ? 51 E y tan ? ) Z'2?1（2-18）其中 ? 满足下式：tan 2? ? 2 ? 51 E y 1 no2?1 ne2（2-19）由于 ? 极小，所以就有以下近似式：tan ? ? 2 ? 51 E y 1 no2?1 ne2（2-20）由式（2-18）可知，在仅对晶体施加Y向电场时，晶体将由单轴晶体变为双轴晶体， 且新的主轴 Y ' 和 Z ' 相对原主轴Y和Z绕X轴转动了 ? 角，如图（2-2）所示图 2-2 仅对晶体施 Y 向电场时 LiNbO 3 晶体的折射率椭球变化新主轴坐标系里沿着三个主轴方向上的折射率分别为7铌酸锂晶体的研究与分析1 ? ' 3 ? n X ? n o ? 2 ? 22 E y n o ? 1 1 ? ' 3 3 ? n Y ? n o ? ? 22 E y n o ? ? 51 E y n o tan ? 2 2 ? 1 ? ' 3 ? n Z ? n e ? ? 51 E y n e tan ? 2 ?（2-21）因为 ? 22 比较小，因此在后面的应用中我们只考虑加Y向电场时晶体光轴的偏转，而忽 略晶体折射率大小的改变。2.3 周期性极化 LiNbO 晶体（PPLN）的制备3铁电体具有自发极化特性（spontaneous electric polarization），其电极化强度与电场强 度间的关系上呈现电滞回线。自发极化Ps的存在与否不取决于外加电场，即使没有外加电 场作用，铁电物质中的自发极化亦能产生。但是外加电场的作用能使自发极化方向反转， 即电畴反转。电畴实际上是一些方向不同的自发极化区域，在每一个这样的区域内，铁电 体的永久偶极子沿同一方向排列，故存在固有电偶极矩。 在铁电体内形成周期性电畴结构是目前为止实现准相位匹配最有效的途径，它通过周 期性的反转铁电晶体的晶向，使得有效非线性系数在 ? d eff 和 ? d eff 之间交替变化，从而实 现非线性系数的空间周期调制。 周期极化LiNbO 3 晶体结构中奇数片电畴与偶数片电畴自发 极化矢量相反，因而这些电畴与奇数阶张量相关的物理性质，如倍频系数、电光系数及压 电系数等的符号亦相反，因此，晶体的物理性质也是空间坐标的周期函数。 实验证明外加电场法是制备周期极化铌酸锂最为有效的方法，它可以实现精确的周期 结构和完全贯穿的垂直电畴壁。其方法是，首先在单畴化铌酸锂晶体的一面（+z面或-z面） 淀积或溅射周期结构的金属电极，另一面制作均匀电极。然后施加与晶体自发极化方向相 反方向的外加电场，当外加电场超过晶体的矫顽场时，其自发极化方向便发生反转。利用 微电子工业的光刻技术，使用干涉测量反馈控制（interferometric feedback control），使得 电极周期结构位置误差限制在很小的范围内，能够实现其他方法难以得到的小周期极化结 构。在周期性电场极化的铌酸锂晶体中，除了非线性系数以外，其他如电光系数，弹光系 数等也同样会由于晶体铁电畴的周期性反转结构得到周期性的调制。 早在1962年，Armstrong和Frallken等人就分别提出了使用周期光栅实现相位匹配这一 概念，但真正将此想法付诸实现，制成可用器件却存在很大困难。为此，科学家进行了不 懈的努力，直到九十年代后，利用外加周期电场调制非线性极化率技术的日趋成熟，周期 极化材料的制备才取得突破进展。这里简单介绍一下周期极化LiNbO 3 晶体的制备方法。 首先在双面抛光LiNbO 3 晶体Z轴表面镀一层金属导电膜， 通常使用Ti、 A1和Cr等金属， 膜厚保持在 100 ~ 300nm 左右。然后，利用半导体光刻工艺制备出周期图案的金属条纹； 随后，在金属条纹电极上涂一层厚的绝缘胶，使各金属电极之间保持良好的绝缘隔离。外 加电场通过液体电极加在LiNbO 3 晶体的金属电极上，也可以将外加电场直接加在LiNbO 3 晶体的金属电极上，所有这些都要保证外电场和金属电极有良好欧姆接触。为防止高压对 空气击穿， 极化过程通常都是在高真空或高压绝缘油中完成。 所用外电场为脉冲高压电场，8铌酸锂晶体的研究与分析对LiNbO 3 晶体，脉冲电压要大于23kV/InIn，脉冲周期长短与次数依具体实验条件而定。 当晶体表面运输电荷达到 Q ? 2 PsA 时（其中Ps为LiNbO 3 晶体自发极化强度，A为极化面 积），开始缓慢降低脉冲电压，持续一段时间，保证已经极化反转的畴不会再自行返回， 最后关掉脉冲电压，完成周期极化过程。目前采用上述方法不仅成功制备了极化厚度达 0.5mm、通光长度超过50mm的均匀周期畴结构的LiNbO 3 晶体。图2-3 周期性极化铌酸锂晶体中的电光效应图（2-3）为对周期性极化铌酸锂晶体施加均匀的Y向电场时晶体电光效应的示意图， 如上一节我们所讨论的，当对铌酸锂晶体施加Y向电场时，晶体的折射率椭球将发生偏转， 也就是晶体的光轴将沿+Z轴偏转 ? 角， ? 由式（2-20）给定。对于周期性极化铌酸锂晶体 来说，由于晶体的周期性畴结构，负畴与正畴的光轴偏转角虽然大小相同，但方向相反， 如图上所示。因此，PPLN晶体上施加均匀的Y向电场之后，晶体的光轴也呈现周期性的偏 转，此种结构正如折叠式的Solc滤波器中晶体光轴的交错排列结构。2.4 新型线性电光效应耦合波理论[10]在上节我们介绍了折射率椭球理论及其在调制中的应用，这种传统的方法直观易懂， 但是存在着不可忽略的局限性:在外加电场的作用下， 电光晶体中的折射率椭球将会随之发 生变化，为了使折射率椭球方程在新的坐标系中主轴化，我们需要找到新旧坐标系的线性 变换，而求得新旧坐标系线性变换的过程大多很复杂，有时甚至是不可能办到的，所以折 射率椭球理论仅适用于某些情况，它对电场方向、入射光的偏振态、入射方向都有着比较 高的要求，如果换做双轴晶体情况，折射率椭球理论就更难被运用。2001年，She等人提 出的线性电光效应的耦合波理论，突破了以上的局限性，可被用于拓展电光材料的选择范 围，优化调制器的调制方式，它的出现引起了电光效应研究领域内的新探索。而本文我们 以2006年发表的新型线性电光效应的耦合波理论展开的。下面我就详细的介绍新型线性电 光效应的耦合波理论的基本内容。 线性电光效应可以看成由光波导和外加电场相互的非线性作用。总的外加电场 E 在线 性电光效应的过程可以被表示为 ? ? ? w E ( t ) ? E ( 0 ) ? [ E exp( ? i ? t ) / 2 ? c .c .] （2-22） ? ? w [ E exp( ? i ? t ) / 2 ? c .c .] 是光学中的频率 ? 。 E ( 0 ) 是直流电场或缓慢渐变电场。c.c.表示 复共轭。总的来说，单色波（频率 ? ）在双折射晶体中传播时存在 2 个独立平面电磁波。? ? ? ? ? w w w E ? E 1 ? E 2 ? E 1 ( r ) exp( ik 1 r ) ? E 2 ( r ) exp( ik 2 r )9（2-23）铌酸锂晶体的研究与分析当 k 1 ? k 2 时 E 1w 和 E 2w 表示为 2 个互相垂直的光场分量，当 k 1 ? k 2 时则表示有着不同的 折射率的 2 个独立光场分量，让? ? ? a,b,c? ?（2-24） 为三个单位矢量， A1 ( r ) 和 A 2 ( r ) 为 2 个波的振幅， n 1 和 n 2 表示 o 光和 e 光的折 射率。我们可以假设晶体是在固定的条件下，由于反压电和光弹性效应被抑制，并且二阶 非线性效应很弱（因为相位不匹配）所以只考虑线性电光效应。由麦克斯韦方程组和以线 性电光效应作为扰动，作慢变振幅近似 耦合波方程组表示为：dA 1 ( r ) dr dA 2 ( r ) dr ? ? i ? g ( r ) A 2 ( r ) exp( i ? k ' r ) ? i ? 1 g ( r ) A1 ( r )? E1 (r ) ?? / n 1 A1 ( r ) a?? E 2 (r ) ?? / n 2 A 2 ( r )b?? ? E (0) ? E 0 c（2-25A）? ? i ? g ( r ) A1 ( r ) exp( ? i ? k ' r ) ? i ? 2 g ( r ) A 2 ( r )（2-25B）这里的 g（r）是材料的结构函数， ? k ' ? k 2 ? k 1 ，且? ?k0 2 n1 n 2 reff 1 E 0?1 ?k0 2 n1reef 2 E 0?2?k0 2n2re f 3f E 0（2-26）如果 g（r）是 r 的周期函数（以 ? 为周期） ，由于周期性电光系数的影响，可以写成 类似傅里叶级数：g (r ) ?m ? ???G??mexp( i ? m r )（2-27）谐波光栅波矢量（倒格矢） ? m ? 2 ? m / ? 十分接近 ? ? k ' 。把（2-26）代入（2-25A） 和（2-25B）， 忽略那些由于相位不匹配而对电光效应贡献很小的成分，我们得到dA 1 ( r ) dr dA 2 ( r ) dr ? ? i ? q A 2 (r) exp(i ? kr) - i ? 1q A 1 (r)（2-28A）? ? i ? q A 1 (r) exp(-i ? kr) - i ?*2qA 2 (r)（2-28B）当 ? k ? ? k '? ? m ，且? q ? ?G m ，? q ? ?G ? m ，*? 1q ? ? 1G 0 ，? 2q ? ? 2G 0（2-29）方程组（2-28A）和（2-28B）是准相位匹配的线性电光效应方程组，设 2 个光振幅为 A 1 （0）和 A 2 （0） ，然后解（2-28A）和（2-28B）得到：A1 ( r ) ? ? 1 ( r ) exp[ i ? r + i ? 1 (r)] A 2 ( r ) ? ? 2 ( r ) exp[ i ( ? ? ? k )r + i ? 2 (r)]（2-30） （2-31）10铌酸锂晶体的研究与分析1? 1 ( r ) ? ? A1 ( 0 ) cos ( ? r ) ? [2 2? ?? A1 ( 0 ) ? ? q A 2 ( 0 ) ??2 2 2 ] sin ( ? r ) ? ?（2-32）? 1 ( r ) ? arg[ A1 ( 0 ) cos( ? r ) ? i? A1 ( 0 ) ? ? q A 2 ( 0 ) ?sin( ? r )]（2-33）1? 2 ( r ) ? ? A 2 ( 0 ) cos ( ? r ) ? [2 2? ? ? ?? A 2 ( 0 ) ? ? A1 ( 0 )* q?*? ? 2 2 ] sin ( ? r ) ? ? ?2（2-34）? 2 ( r ) ? arg[ A 2 ( 0 ) cos( ? r ) ? i? ? A 2 ( 0 ) ? ? q A1 ( 0 )?) ? 4? q ? q2 *sin( ? r )]（2-35） （2-36）? ?? ?1 2 (?1 2( ? k ? ? 1q ? ?2q2q? ? 1q ? ? k )? ?1 2( ? k ? ? 1q ? ?2q)（2-37）此解可以用来描述，在任意方向的外加电场的作用下，任意偏振态的入射光在任意点 群的电光晶体中沿任意方向传播时的情况。2.5 本章小结本章主要详细介绍电光效应发生的原理和过程，以及线性电光效应耦合波理论的主要 内容。我们知道在电光效用过程中，电光系数，o 光和 e 光折射率，温度，波长，外加电 场等参数会对电光效应的转换效率产生非常大的影响。后面章节的仿真内容将以线性电光 效应耦合波理论为基础建立模型。11铌酸锂晶体的研究与分析第 3 章 LiNbO 的晶体结构和性质33.1 LiNbO 晶体结构[1]3自1965 年Ballman成功的利用Czochralski提拉法生长出铌酸锂单晶后，铌酸锂晶体得 到了广泛的研究。铌酸锂是目前以知的居里点最高（ 1210 ? C ），自发极化最大（室温时约 ? 2 0 . 7 C / m ）的铁电体，顺电相和铁电相的空间群分别为 R 3 C 和 R 3 C ，其结构如图（3-1） 所示。a） 铁电相 水平线代表氧平面 图3-1 铌酸锂晶体结构示意图b） 顺电相氧八面体以共面的形式叠置起来形成堆垛，公共面与氧八面体三重轴（即极轴）垂直。 许多堆垛再以八面体共棱的形式连接起来形成晶体。在顺电相，Li和Nb分别位于氧平面和 氧八面体中心，无自发极化。在铁电相，Li和Nb都沿c轴发生位移，前者离开了氧八面体 的公共面，后者离开了氧八面体的中心。由于Li和Nb的移动，造成了沿c轴的电偶极矩， 即出现了自发极化。 该结构也可以看成由垂直于极轴且相互等距的氧平面组成。 顺电相时， Nb位于两个氧平面中央，Li位于第三个氧平面内（实际上Li分布于氧平面和氧平面上下各 0.037nm处，其平均位置在氧平面）。铁电相时，Nb和Li都沿+c轴移动。结构分析表明， 室温时，Nb沿+c轴偏离氧八面体中心约0.026nm，Li沿+c轴偏离氧平面0.044nm。下面只介 绍与极化有关的铁电相。铁电相的LiNbO3晶体含有一个三重对称轴，属三角晶系。此外， 它还有一个对称面， 三个成60°角平面相交形成一个三重旋转轴。 这两个对称操作LiNbO 312铌酸锂晶体的研究与分析晶体归类为3m点群（C6v），它也属于空间群。在三角晶系中，可选择两种完全不同的晶 胞:六方晶胞和三角晶胞。 对于惯例的LiNbO 3 的六方晶胞， c轴被定义为晶体的三重旋转轴。 确定c轴方向的标准方法是：在c轴方向压缩晶体，显负电性的面为+c；确定+c轴第二种方 法是冷却晶体，显正电性的为＋c方向。两种方法可从Li、Nb离子与氧八面体的相对运动 进行理解。当受挤压时，Li、Nb离子都向接近于顺电相的方向发生位移，减小了自发极化， +c面的负电荷过剩而使晶面呈负电性。当晶体冷却时，离子的热能降低，弹力把Li、Nb离 子推得远离氧八面体中心及邻近的氧平面， 增强了晶体的自发极化， 使晶体+c面呈正电性。 在1966年精确确定晶体结构之前，人们不知道铌酸锂化学计量中可能存在的偏差。铌 酸锂的晶格参数与精确化学组成的依赖关系是于1968年建立起来的。说明某晶体化学计量 比的很可靠、很精确参数之一是居里温度。通过比较已知化学计量样品的居里温度与待测 铌酸锂样品的居里温度，能够极好地确定样品的化学组成。 根据晶体结构可解释铌酸锂的晶格常数——热膨胀特征。现已发现，温度升高，铌氧 八面体的倾斜度增大，其原因是六方晶格参数a的热膨胀几乎是线性的。在 600 ~ 1000 ? C 温 度范围内六方晶格参数c的收缩，是由于随着Nb离子朝着仲电相位置的移动，八面体的边 长缩短。3.2 LiNbO 晶体基本性质3LiNbO 3 晶体是一种无色或淡黄色的透明晶体，其莫氏硬度为5，和软玻璃相似，它的 努氏显微硬度值为600，在（001）方向硬度值大约高25%。LiNbO 3 晶体能够被普通的金刚 石道具切开，用普通的光学加工技术也能够很好的完成晶体的研磨和抛光。 在 4 ? C ，LiNbO 3 晶体密度为 4644kg/m 3 。其居里温度很高约为 T c ? 1210 ? C ，仅仅比 _ 其熔点低几十度。 在此温度以上晶体属三方晶系 3 m 点群，为顺电相；在居里温度 Tc 以下， 晶体属三方晶系 3m 点群（可用六方晶系来表示） ，为铁电相。由于 LiNbO 3 晶体的居里温 度很高，因而又称为高温铁电体，它具有良好压电性，热释电性，铁电性，电光和非线性 光学性能，又是多功能的晶体材料。LiNbO 3 单晶的介电系数随温度 T 升高而增大，在 0 . 4 ? m ~ 5 ? m 波长范围内，可连续通光。3.3 LiNbO 晶体特点3LiNbO 3 晶体在集成光学和光波导应用中是一个重要的材料， 尤其是近些年来， 稀土掺 杂工程，畴工程和近化学比晶体生长鱼加工技术的完善使得有关于 LiNbO 3 波导的光电子 器件的的功能和性能的研究急剧增加。其具有以下的特点： （1）优良的电光，双折射，非线性光学，声光，光折变，压电，热释电，铁电与光生伏 打效应等物理特性。 （2）机械性能稳定，耐高温，抗腐蚀。 （3）易于生长大尺寸晶体，容易加工，成本低。 （4）实施不同掺杂后能呈现出各种各样的特殊性能，使之在光波导，电光调制器，倍频 转换，全息存储等方面有着广泛应用。13铌酸锂晶体的研究与分析3.4 PPLN 晶体的应用周期性极化 LiNbO 3 （PPLN）材料是技术含量很高的非线性光学频率转换晶体。1998 年为美国国家研究理事会所编:《HarnessingLight-optical seience and Engineering for 21st century》一书中把 PPLN 材料及其应用作为下世纪，非线性频率变换材料的唯一重点研究 对象，建议国家重点资助。它通过倍频、光参量放大和振荡、差频等二阶非线性光学过程， 将来广泛应用于光传输、光存储、光显示和遥感探测等方面。其主要用途有: （l）光存储:通过倍频转换得到的短波长光源，可以用于高密度的光存储，是蓝绿光半导 体激光器的有力竞争者。 （2）光显示：蓝绿光光源作为高纯度三元色可以用于高清晰度显示。 （3）全光通讯：利用差频效应，可以制作出未来全光 DWDM 通讯系统中的关键器件一波 长转换器。与其它类型波长转换器相比，它具有在通讯系统中严格透明的优点。 （4）遥感、探测、生物医学等：利用参量放大和振荡产生可调谐近、中红外光源。应用 于空间分子探测及其它军事方面的应用。另外，小型红外光源在医学、科研方面均有很大 的应用场合。 （5）其它应用：电光调制器、电光偏转器和电光透镜等。3.5 LiNbO 晶体折射率方程3LiNbO 3 晶体在光学上为单轴晶体， 不同于正单轴晶体 n o ? n e ） （ LiTaO 3 的是， LiNbO 3 为负单轴晶体（ n o ? n e ） ，一般条件下，LiNbO 3 在 0.4 ~ 5 ? m 的波长范围内均是无色透明 的，在补偿晶体界面的反射损失时，投射率可达 74%。 LiNbO 3 晶体在氢气中被加热到 670 ~ 870K 后，会由最初的无色透明变为褐色。在 ? ? 0 . 5 ? m 和 ? ? 0 . 8 ~ 1 . 1 ? m 处出现两个新的吸收带，并且在 ? ? 2 . 87 ? m 处形成强的吸收 带。晶体在空气中退火并极化后呈浅黄色。 LiNbO 3 晶体在一些常用激光器的输出波长和几个其他波长处的寻常折射率 n o 和异常 折射率 n e 对温度的依赖关系见表（3-1） 。 通过对实验数据分析可以得到在波长为 400 ~ 4000nm ，计算 LiNbO 3 晶体在不同温度和波 长下的 Sellmeier 方程[1]为：n o ? 4 . 9130 ?21 . 173 ? 1025? 1 . 65 ? 102?2T2 2 2? ? ( 2 . 12 ? 10?? 2 . 7 ? 105?5? 2 . 78 ? 10?8?2（3-1）T )n e ? 4.5567 + 2.605 ? 102-7T20 . 970 ? 102? 2 . 70 ? 102?2 5T2 2 2? ? ( 2 . 01 ? 10? 5 . 4 ? 10 T )? 2 . 24 ? 10?8?2（3-2）式中，T 为绝对温度（K） ? 是以 nm 为单位的波长。 ，表 3-1 LiNbO 3 晶体对不同波长的折射率14铌酸锂晶体的研究与分析? /nm激光 He-Cd Ar Ar Ar Ar Ar Ar Ar Ar Ar He-Ne 红宝石 GaAs Nd He-Ne化学计量比（T=25 C） no ne 2.5 2.0 2.3 2.9 2.3 2.6 2.0 2.1506?同成分熔体（T=24.5 C） no 2.5 2.7 2.9 2.1 2.7 2.6 2.3 2.2225 ne 2.0 2.2 2.1 2.6 2.5 2.9 2.1 2.1519?441.6 457.9 465.8 472.7 476.5 488.0 496.5 501.7 514.5 530.0 632.8 693.4 840.0 0.02.6 2.6 2.3 2.5 2.0 2.0 2.2 2.23203.6 本章小结LiNbO 3 在集成光学和光波导应用中是一个重要的材，为应用最广泛的电光材料之一。 本章我们系统全面的介绍LiNbO 3 晶体结构组成和本征。尤其是最后提到的LiNbO 3 晶体的 Sellmeier方程，将再matlab仿真中多次被用到，是至关重要的一个方程，因此需要重点注 意。15铌酸锂晶体的研究与分析第4章 系统结构和参数设定4.1 相关参数说明转换效率? ：由式（2-30） 和式（2-31）我们已经知道了耦合波方程组的解析解，其 中 A1 ( r ) 和 A 2 ( r ) 分别为 o 光和 e 光的波动振幅。我们假设入射光为 o 光，设定归一化后的 振幅 A1 ( 0 ) ? 1, A 2 ( 0 ) ? 0 。由于 A1 ( r ) 和 A 2 ( r ) 是复数形式，所以转换效率? 采用其模值之比， 即? ?A2 (r ) A1 ( r )（4-1）线性耦合波理论中考虑的情况由于反压电和光弹性效应被抑制，并且二阶非线性效应 很弱（因为相位不匹配） ，所以所以只考虑线性电光效应。这样的转换效率? 理论最大值可 以近似达到 1。但实际运用过程当中最大值仅仅为 40 % ~ 50 % 左右。 占空比系数 D：由图（2-3）我们知道周期极化的晶体中含有若干个极化周期 ? 。而每 一个极化周期 ? 中又由正畴和负畴两部分组成，其长度我们分别用 L ? 和 L ? 表示。而占比 系数 D ? L? /? （4-2）4.2 PPLN 结构参数设定PPLN 的结构函数为? ? 1 　　　（ r 在正极） g (r ) ? ? ? ? 1 　　　（ r 在负极）（4-3）其傅立叶变换为Gm ? 1 [1 ? cos( 2 ? mD ) ? i sin( 2 ? mD )] 　　　　 ? ? ? i? m ? 2 D ? 1 　　　　　　　　　　 　　 ? （ m ? 0） （ m ? 0）（4-4）其中，m 为准相位匹配阶数。经大量研究发现占空比系数 D 的最适值为 0.5 （当 m=1 时） ，0.25 或 0.75（当 m=2 时） ，当 D 取这些值时 G m 有最大值。 我们假设在常温 T=293K 下制作 PPLN 晶体，入射光线我们选择 ? ? 1550nm 波长的 半导体激光器，由寻常光线和异常光线的 Sellmeier 方程式（3-1）和（3-2） ，我们可以计 算出常温下 PPLN 晶体中两种光线的折射率为：n o ? 2 . 2128 , n e ? 2 . 137316铌酸锂晶体的研究与分析我 们 知 道 入 射 光 波 矢 量 k 0 ? 2? / ? 0 ， 所 以 o 光 和 e 光 的 相 位 失 配 量? k ? k1 ? k 2 ? k 0 (n o ? n e ) ? 2? ( n o ? n e )?0我 们 所 构 造 的 PPLN 晶 体 的 倒 格 矢 ? m ? 2 ? m / ? ， 为 了 满 足 相 位 匹 配 条 件 ， 使 ?k ? ? m ? 0 。 我们考虑在准相位阶数 m=2 的情况情况下， 此时假设 PPLN 晶体的占空比系 数 D=0.25，可以计算出 PPLN 晶体的极化周期：? ? m?0 no ? ne ? 2 ? 1550 ? 10?92 . 2128 ? 2 . 1373? 4.1094 ? 10?5m其中负畴长度为 0 . 75 ? ? 3.0821 ? 10 ? 5 m ，正畴长度为 0 . 25 ? ? 1.0273 ? 10 ? 5 m 。我们设定含 有 100 个极化周期的 PPLN 晶体，其长度 r ? 100 ? ? 4.1094 ? 10 ? 3 m 。4.3 有效电光系数设定由前文我们已经知道 ? 12 ? ? 61 ( ? ? ? 22 ), ? 13 ? ? 23 , ? 22 , ? 33 , ? 42 ? ? 51 为 LiNbO 3 晶体的非零 电光系数。 而有效电光系数定义为：? ? eff 1 ? ? ( ? jj ? kk )( a j ? jkl b k c l ) ? j , k ,l ? ? ? eff 2 ? ? ( ? jj ? kk )( a j ? jkl a k c l ) j , k ,l ? ?? ? ? ( ? jj ? kk )( b j ? jkl b k c l ) eff 3 ? j , k ,l ? (2) 2 x jkl ( ? ? , ? , o ) ? ? ? jkl ? ? ? jj ? kk ?（4-5）单下标对双下标的约化为： xx→1，yy→2，zz→3，yz=zy→4，xz=zx→5，xy=yx→6，x→1，y→2，z→3 （4-6） 由式（4-6）我们转换下标得到 LiNbO 3 晶体的非零电光系数有：? ? 112 ? ? ? 113 ? ? ? 222 ?? ? 333 ? ? 232 ? ? ? 121 ? ? ??223 211? ? 3 .4? 8 .6? 3 .4 ? 30 . 8 ? ? 322 ? ? 131 ? ? 311 ? 28（4-7）又知道单轴晶体的独立偏振光场分量的单位矢如下（光不沿光轴方向传播） ，坐标系 Z 轴为晶体光轴所在方向，需要特别指出的是，当入射光线沿着晶体光轴方向传播时，将不17铌酸锂晶体的研究与分析会发生双折射现象。其中入射波矢的角度 ? , ? 由图（4-1）给出：图 4-1 入射波矢角度关系图? ? cos ? cos ? ? ? ? ? o 光： e 光： b ? ? ? cos ? sin ? ? ? sin ? ? ? ? ? ? 我们假设入射 o 光波矢角度 ? ? 0 , ? ? 0 . 5? ，则 a ? ( 0 , ? 1, 0 ), b ? ( 0 , 0 ,1) ，即入射光线沿 ? 轴正向。而所外加的正弦（直流）电场沿 Y 轴正向，即 c ? ( 0 ,1, 0 ) ，此时由于电场方向 ? s i n? ? ? ? ? a ? ? ? c o s? ? ?0 ? ? ?X 垂直于光传播方向，我们可以进行横向电光调制。比起纵向电光调制（电场方向平行于光 传播方向）我们可以发现横向调制更符合实际需要，但由于存在自然双折射而引起的电光 延时，使得对于温度的变化非常敏感。? a 1 ? 0 , a 2 ? ? 1, a 3 ? 0 ? 由上我们知道 ? b 1 ? 0 , b 2 ? 0 , b 3 ? 1 ? c ? 0 , c ? 1, c ? 0 2 3 ? 1由于矢量 a，b，c 三项每一项都只存在一个数值，我们可以得出有效电光系数 ? eff 1 , ? eff 2 , ? eff 3 中 分 别 只 可 能 存 在 ? 232 , ? 222 , ? 332 一 项 ， 这 样 可 以 计 算 出 2 2 4 4 ? eff 1 ? n e n 0 E 0 ( a 2 b 3 ? 232 c 2 ), ? eff 2 ? 2 n 0 E 0 ( a 2 a 2 ? 222 c 2 ), ? eff 3 ? 2 n e E 0 ( b 2 b 2 ? 222 c 2 ) 。 而对于每一次改变入射波矢角度 ? , ? 和外加电场的方向所对应的有效电光系数 ? eff 1 , ? eff 2 , ? eff 3 都将改变而需要重新运算，本文没有将其列入考虑范围。具体的数值计算将 在 matlab 仿真中同步进行。4.4 本章小结本章我们重点构建了 PPLN 系统参数的极化周期 ? ，从而达到相位匹配；并且在我们 所设定的入射光线角度和外加电场方向下给出了所对应的有效非线性系数 ? eff 1 , ? eff 2 , ? eff 3 的 求值方法，从而建立了完整的 LiNbO 3 晶体线性电光效应过程的数学理论模型，为下一章 的 matlab 线性仿真提供了系统模型。18铌酸锂晶体的研究与分析第 5 章 线性仿真与讨论5.1 温度 T 的改变对转换效率? 的影响由图（5-1）我们发现在所测定的温度范围内当 T=293K 时 ? k ? 0 ，以此温度为中心左 右分别出现相位的负，正失配。在数值上是同样是以 T=293K 为中心对称的。而由于温度 改变导致相位失配从而对电光效应中转换效率? 的影响我们将由图（5-2）给出。400020000相位失配量k-2000-4000-6000-8000-280290300310 温 度 T(K)320330340350图 5-1 用 PPLN 弥补后的相位失配量与温度 T 的关系从图（5-2）我们能够发现转换效率? 约在 T ? 292 ~ 293 K 之间时达到峰值，这与我们 所设定的 PPLN 晶体在 T=293K 温度下达到相位匹配存在极其细微的偏差。考虑到我们所 计算的 PPLN 系统结构参数并不是一个精确值，而是近似四舍五入的只保留到小数点后四 位， 这样一来当 T=293K 时， 相位匹配量 ? k 只是近似的等于零。 而真正实现相位匹配 ? k ? 0 的温度值则是我们图中波形出现最大峰值时所对应的温度值。 仿真结果表明随着温度 T 的改变转换效率? 并不是呈单一的曲线上升或下降，而是以 温度 T=293K 为中心向两侧对称性减小。我们能够联系到温度 T 的改变导致了 o 光和 e 光 折射率的改变，其最终结果是导致 ? k 的改变，从而使得相位失配。随着相位失配量 ? k 的 增大，图中即表现为温度上与 T=293K 的差值的增大，转换效率? 逐渐降低。 当温度为绝对零度（273K）时转换效率依然存在，这一结果仅为 matlab 线性仿真的理19铌酸锂晶体的研究与分析论结果，现实中由于无法实验室达到绝对零度而并不存在。1 0.9 0.8 0.7 0.6转换效率0.5 0.4 0.3 0.2 0.1 0 270280290300310 320 温 度 T（ K）330340350图 5-2 温度 T 的改变对转换效率? 的影响而转换效率 ? 不是线性减少而是成波动形式趋于零则是考虑到是由于有效电光系数 2 2 4 2 ? eff 1 , ? eff 2 , ? eff 3 的计算分别和 n o n e , 2 n o , 2 n e 成线性关系而不是和 n o , n e 呈线性关系所造成，我 们可以参照图（5-3）和图（5-4） 。 图（5-3）显示 n o , n e 随着温度 T 增加而相应的增加。但由于 LiNbO 3 为负单轴晶体所 以表现随着温度 T 增加 n o 增长率大于 n e 。2.9 2.139 2.9 2.5 2.13952.138折射率2.8 2.8折射率280 290 300 310 320 温 度 T（ K） 330 340 3502.13752.1372.8 2..136 270280290300310 320 温 度 T（ K）330340350a） 寻常光线 b）异常光线 图 5-3 温度 T 对 o 光和 e 光折射率的影响20铌酸锂晶体的研究与分析而图（5-4）则显示有效电光系数 ? eff 1 , ? eff 2 , ? eff 3 随着温度 T 增加的关系。我们发现并不 是单一的上升和下降，由于上升量和下降量的差值导致了转换效率 ? 不是线性减少而是成 波动形式减少。-6.25 x 10-101.6306x 10-10-8.98x 10-17-6.2551.6306-8.99有 效 非 线 性 系 数 reff1（ m/V）有 效 非 线 性 系 数 reff2（ m/V）-6.26有 效 非 线 性 系 数 reff3（ m/V）280 290 300 310 320 温 度 T（ K） 330 340 3501.6305-9-6.2651.6305-9.01-6.271.6304-9.02-6.2751.6304-9.03-6.28 270280290300310 320 温 度 T（ K）3303403501.-9.04 270280290300310 320 温 度 T（ K）330340350a） ? eff 1b） ? eff2c） ? eff3图 5-4 温度 T 对有效电光系数的影响接下来我们考虑改变 PPLN 晶体极化周期数量 n 后温度 T 的改变对于转换效率? 产生 的影响。1 0.9 0.8 0.7 0.6 1 0.9 0.8 0.7 0.6转换效率0.5 0.4 0.3 0.2 0.1 0 270转换效率280 290 300 310 T(K) 320 330 340 3500.5 0.4 0.3 0.2 0.1 0 270280290300310 T(K)320330340350a） n=50 图 5-5 极化周期数量对转换效率? 的影响b） n=150由图（5-5）能够看出最大转换效率依然发生在 T=293K 的时候，而不同的是图（5-5a） 的波动幅度较小，波动趋势较缓慢；而图（5-5b）的波动幅度较大，波动趋势较剧烈。由 此我们可以的出一个结论：相同的温度 T 的改变量，随着晶体极化周期数量 n 的增加，转 换效率? 的波动将更加剧烈， ? 值也降低的更快，波动次数将增多。21铌酸锂晶体的研究与分析5.2 波长 ? 的改变对转换效率? 的影响由图 （5-6） 我们发现波长 ? 的改变， 类似于前面讨论的温度 T， 大约在我们所设定的 PPLN 结构参数所对应的波长，即 ? ? 1550 nm 时，转换效率? 达到最大值，并且以此波长为中心 向两侧成峰值递减波动形式， 最终趋近于零。 我们还发现靠近理论波长附近， 转换效率? 的 下降趋势极其明显， 这表明从相位匹配变化到相位失配这一过程中，? k 极其细微的变化将 对转换效率? 产生非常巨大的影响。1 0.9 0.8 0.7 0.6转换效率0.5 0.4 0.3 0.2 0.1 0 145015001550 波 长 (nm)16001650图 5-61 0.9 0.8 0.7 0.6波长 ? 对转换效率? 的影响1 0.9 0.8 0.7 0.6转换效率0.5 0.4 0.3 0.2 0.1 0 1450转换效率 波 长 （ nm） 0.5 0.4 0.3 0.2 0.1 0 145015001550 波 长 （ nm）16001650a）n=50b）n=150 图 5-7 极化周期数量对转换效率? 的影响22铌酸锂晶体的研究与分析而改变 PPLN 晶体极化周期数量 n 后波长 ? 对于转换效率? 产生的影响。我们由图 （5-7）给出。对比图（5-5）和图（5-7）我们不难发现晶体极化周期数量 n 的增加，使得 转换效率? 的波动更加剧烈，波动次数将明显增多。 由此我们可以得出结论，波长 ? 和温度 T 的改变，同样是由于改变了 o 光和 e 光的折 射率从而造成相位失配，进而造成转换效率? 的改变。其本质原因都是相位失配量 ? k 值的 改变造成转换效率? 的变化，所以情况基本相同。我们可由图（5-8）看出正负相位失配量 ? k 的值对于转换效率 ? 的改变是对称关系的。1 0.9 0.8 0.7 0.6转换效率0.5 0.4 0.3 0.2 0.1 0 - - -0 相位失配量k2000300040005000图 5-8 相位失配量 ? k 对转换效率? 的影响2.218 2.217 2.216 2.2152.142 2.141 2.14 2.139折射率折射率 波 长 (nm) 2.214 2.213 2.212 2.211 2.21 2.209 14502.138 2.137 2.136 2.135 2.134 145015001550 波 长 (nm)16001650a）寻常光线b）异常光线图 5-9 o 光和 e 光随波长 ? 变化的折射率曲线23铌酸锂晶体的研究与分析至于图（5-2）和图（5-6）波动形式上的细节区别则是考虑到在式（3-1）和式（3-2） 所给出的折射率方程中，波长 ? 和温度 T 的改变对于 n e 和 n o 的改变形式不同而造成，这一 点我们可以比较图（5-3）和图（5-8） ，我们能够发现不同于温度 T， n o , n e 随着波长 ? 增加 反而线性减少。5.3 外电场 E 的改变对转换效率? 的影响我们给定外加电场为正弦电场 E ? 0 . 9 ? 10 5 sin( 0.004 ? t ? 0.0625 ? ) V / m ，其中时间 t 取 ms 为单位，其波形图由图（5-9）给出。1 0.8 0.6 0.4 x 1051 0.9 0.8 0.7 0.6电 场 强 度 E（ V/m）0.2 0 -0.2 -0.4 -0.6 -0.8 -1转换效率0 200 400 600 800
t（ ms） 00 20000.5 0.4 0.3 0.2 0.1 00200400600800 t（ ms）1400160018002000图 5-10 外正交电场波形图 5-11 电场参数 t 对转换效率? 影响我们通过比较图（5-10）和图（5-11） ，不难发现转换效率? 随着时间 t 的改变成波动 形式，波动的周期为外加正弦电场周期的一半，转换效率 ? 的最大值出现在外电场曲线出 现波峰和波谷的时候， 由此我们可以得出转换效率? 与外加正弦电场的幅度值成线性关系， 当幅值增大时转换效率? 也相应的增加，反之则减小。 这可以证实线性电光效应表现为介质折射率同外加电场成线性变化。如果外加电场为 直流电场，那么转换效率? 将正比与直流电场的场强数值大小。5.4 晶体占空比 D 的改变对转换效率? 的影响前面几节的讨论我们是在假定占空比 D=0.25 的情况下进行的。现在我们保持准相位 阶数 m 和 PPLN 晶体的极化周期 ? 不变而改变其占空比 D 来讨论其对转换效率的影响。 由图（5-12）我们发现转换效率? 的曲线是以占空比 D=0.5 为中心而近似对称的两个 半波组成， 当转占空比 D=0.25 和 0.75 时转换效率? 取最大值， 而当 D=0.25 时转换效率? =0， 此时无法发生 o 光和 e 光间的相互转换。 其原因可见分析式（4-4）中结构函数 G m ，由于线性电光效应耦合波方程中占空比 D 仅仅与傅立叶变换后的结构函数 G m 有关系。分析是由于其中的三角函数项最终造成转换 效率? 成两个半波形式出现。24铌酸锂晶体的研究与分析1 0.9 0.8 0.7 0.6转换效率0.5 0.4 0.3 0.2 0.1 000.10.20.30.4 0.5 0.6 晶体站空比D0.70.80.91图 5-12 占空比 D 对转换效率? 的影响5.5 本章小结本章基于线性电光效应耦合波理论，在第四章所设定的系统结构和参数下运用 matlab 进行线性仿真，在一定的误差允许范围内。发现随着与我们所设定的系统温度和波长的差 值越大，转换效率 ? 的值成波动形势逐渐下降，在系统参数 T ? 293 K , ? ? 1550 nm 时，转换 效率最大。 此外 PPLN 晶体极化周期数量 n 后温度 T 对于转换效率? 产生不可忽略的影响。 而外加电场 E 对转换效率? 的影响则是成正比的线性关系。从根本上解释则是由于相位失 配量 ? k 发生的变化导致转换效率 ? 的改变。另外数值结果表明若要取得最大的转换效率 ? ，当准相位阶数 m=2 时，我们所设计的 PPLN 结构应该采用占空比 D 为 0.25 或 0.75。25铌酸锂晶体的研究与分析结论本文我们由浅入深，由理论到仿真，首先学习 LiNbO 3 晶体电光效应的相关理论知识， 着重研究线性电光效应耦合波理论，并以此为基础进行参数设定和线性仿真。我们假设在 常温 T=293K 下制作 PPLN 晶体，入射光线我们选择 ? ? 1550nm 波长的半导体激光器，由 此设定 PPLN 的极化周期 ? 。 另外一方面我们选择入射角度 ? ? 0 , ? ? 0 . 5? 的寻常光线进行 入射，由此我们可以进行有效电光系数 ? eff 1 , ? eff 2 , ? eff 3 的计算。 从 matlab 线性仿真的数值结果我们可以得出结论： 当温度 T 和波长 ? 等于我们所设定 的值时转换效率? 值最大，此时 o 光可以完全转换为 e 光。随着温度 T 和波长 ? 距离我们 所设定的理论值的差距的增大，转换效率 ? 将呈波峰值逐渐降低的波动形式渐渐趋近零。 造成波动的原因则是由于有效电光系数 ? eff 1 , ? eff 2 , ? eff 3 的计算分别和 n o2 n e2 , 2 n o4 , 2 n e2 成线性关 系而不是和 n o , n e 呈线性关系，并且是有增有减。此外随着晶体极化周期数量 n 的增加，转 换效率? 的波动将更加剧烈，? 值也降低的更快， 波动次数将增多； 而电场 E 对转换效率? 的影响则是随着电场 E 幅度值的增减转换效率? 相应的增减。研究表明转换效率? 的改变 实际上是由于参数的改变导致了相位失配量 ? k 产生变化， ? k 值越大转换效率 ? 就越小。 所以我们进行电光调制的过程中，必须调节好系统参数从而使得相位匹配 ? k ? 0 ，这样 o 光与 e 光的转换效率将达到最为理想的状态。此外研究表明占空比约在 D=0.25 和 0.75 时 转换效率? 取最大值，并以 D=0.5 为中心呈两个半波形式，其中 D=0.5 时转换效率最小， 几乎不发生 o 光与 e 光的相互转换。 由于时间有限，本文只初步研究了温度 T，波长 ? ，外加电场 E 和晶体占空比 D 对转 换效率? 的影响，至于电光效应中入射光线的角度和外加电场的方向等参数对转换效率 ? 所产生的影响也是不能忽略的。本文未能更深入讨论，这有待进一步研究。26铌酸锂晶体的研究与分析参考文献1 孔勇发,许京军,张光寅,陆猗.多功能光电材料——铌酸锂晶体[M].科学出版社.2005 2 叶佩弦.非线性光学[M].中国科学技术出版社.1999 3 张一兵.铌酸锂的晶体结构[J].上饶师范学院学报.):52-56 4 张旭，薜冬峰.铌酸锂晶体的结构与性能关系研究[J].信息功能材料.):53-57 5 陈洪云.准位相匹配非线性光学中基于周期极化铌酸锂的电光效应研究[D].江西师范大学 硕士论文. 6 董孝义，高希才.电光学及其应用[J].半导体光电.):189-201 7 李家泽，朱宝亮，魏光辉.晶体光学[M].北京理工大学出版社.1989 8 纪磊.周期极化铌酸锂的制备及应用研究[D].天津大学博士论文.44-50 9 R.T.White,W.P.Bowen,T.Mekinnieete.Effieient tunable ultraviolet generation in periodically poled lithium niobate[J].ELECTRONICSLETTERS.):642-643 10 Guoliang Zheng,Hongcheng Wang,Weilong She.Wave coupling theory of Quasi Phase Matched linear electro-optic effect[J].OPTICS EXPERESS.):27铌酸锂晶体的研究与分析致谢本文是在钟东洲老师的悉心指导和热情关怀下完成的，他严以律己、宽以待人的处事 风范，朴实无华、乎易近人的人格魅力，令人如沐春风，使我明白了许多待人接物和为人 处世的道理。钟东洲老师在许多具体问题上给予了巨大的帮助，与钟老师愉快的讨论使我 对一些问题的朦胧的看法最终转化为可行的方案，让我对这一课题有了更加深入的认识， 将本课题推向细致、深入，感谢钟老师的教导和包容。在此向钟老师表示诚挚的谢意。 快乐，紧张和兴奋的四年大学生活即将结束，回首过往，在我四年的大学学习生涯中， 感谢曾经教导过我的各位老师，是你们给了我丰富的专业知识，让我收益非浅。感谢我同 一届的舍友张文林，詹火炬，康小伟，在跟他们一起的学习生活中，我看到了自己的不足， 积极改进才有了今天的进步。 同样感谢同班的各位同学， 已毕业和未毕业的各位师兄师姐， 现在的师弟师妹，他们平时对我在学习上和生活上的帮助都很多，我就不一一尽述了，在 这里，向你们表达我最诚挚的谢意和最衷心的祝福。28
铌酸锂晶体的研究与分析——提供以文本文档的格式的各类文档免费下载和在线浏览。