高斯光束研究.doc

高斯光束通过非线性介质的自聚焦现象华东师范大学理工学院数学系2009届学士学位论文高斯光束通过非线性介质的自聚焦现象The Self-focusing phenomenon of Gaussian beams passing through a nonlinear medium姓 名： 王睿 学 号： 10051520118 班 级： 2005级4班 指导教师： 刘永明 职 称: 教授 2009年5月完成目 录摘要 3Abstract 3一、引言 4二、预备知识 4（一）波动方程 4（二）赫姆霍茨方程6（三）高斯光束表达式推导7 三、问题研究 10（一）光束在介质上表面折射透射11（二）光束在介质中的传播14 四、模型建立及分析 17 参考文献 24附录程序25致谢 27高斯光束通过非线性介质的自聚焦现象摘要：随着信息技术和纳米技术的迅速发展，要求光信息存储器件中的最小信息位尺寸、大规模集成电路和微电子技术中的光刻线宽和光学显微镜的分辨率等均能达到纳米量级（100nm），而由于光衍射本身的限制，无法达到实际需求。非线性薄膜材料的研究，通过选择非线性强的光学薄膜材料，调节激光能量和控制薄膜厚度及结构，在非线性薄膜结构的出射面能使光斑尺寸进一步下降，实现纳米光斑。该光斑通过近场耦合作用在信息存储薄膜或光刻薄膜上，从而实现纳米信息存储、纳米光刻或纳米成像。本文主要研究高斯激光束通过非线性均匀绝缘介质后光强的改变。由电磁场基本原理，推导出高斯光束是缓变振幅条件下波动方程的近似解，研究其在介质突变面处的反射透射。重点研究高斯激光束在非线性介质中的传播问题，这一过程中有自聚焦现象。研究过程主要采用数值计算方法用差分方程代替偏微分方程研究问题的数值解。比较光强的变化。关键词：高斯光束，非线性，自聚焦，差分方程Abstract：With the development of information and nano technology, it is desirely requried that some aspects can reach as exactly small as the scale of nano, including the the Size of the minial bit of information in the Pieces of optical information storage, Lithography linewidth in Large-scale integrated circuits and Microelectronic technology, Resolution optical microscope. However, due to the natural restrictions of Ray diffraction, it can be hardly satisfied as it practical required. While the research of nonlinear thin film materials can realize the obtaination of the nano facula, which uses a nonlinear thin film materials as well as the modulation of Laser energy and the thickness along with the structure of a film control thus the exit surface of the nonlinear thin film can furhter decrease the spot size. The facula can have effect on information storge or lighography film through the Near-filed coupling so that it can be applied into Nano-information storage, nano-lithography or nano-imaging. In this article, we majorly have research on the change of the light power when Gaussian beams passes through the nonlinear insulation thin film。Based on the principle of Electromagnetic field，we conclude that Gaussian beams is the approximate result of Wave equation under the condition of slowly varying amplitude，we also focus on its reflex transmission at the point of Medium surface mutations. While the most important aspect that we have research on is the issue of transpotation of Gaussian beams in a nonlinear medium，in which there is a phenomena of the self-focusing. We use Numerical method as a primary method, in which we use Difference equation instead of Partial Differential Equations to calculate the numerical solution of the issue and consequently compare the change of the light intensity. Key words: Gaussian beams, nonlinear, Self-focusing, Differential Equations一、引言随着信息技术和纳米技术的迅速发展，要求光信息存储器件中的最小信息位尺寸、大规模集成电路和微电子技术中的光刻线宽和光学显微镜的分辨率等均能达到纳米量级（100nm），而由于光衍射本身的限制，无法达到实际需求。而通过非线性薄膜材料的研究，通过选择非线性强的光学薄膜材料，调节激光能量和控制薄膜厚度及结构，在非线性薄膜结构的出射面能使光斑尺寸进一步下降，实现纳米光斑。该光斑通过近场耦合作用在信息存储薄膜或光刻薄膜上，从而实现纳米信息存储、纳米光刻或纳米成像。实验中我们常常采用高斯光束作为光源进行问题研究。高斯光束是波动方程在缓变振幅下的一个特解，非线性介质的折射率随光强的变化而变化，因而高斯光束通过非线性介质发生自聚焦和衍射现象，从而改变能量分布。本文主要研究光强的变化，通过具体数值建立数学模型，采用差分方程代替偏微分方程以求得问题的数值解，研究光束通过非线性介质后能量的变化。二、预备知识（一）波动方程波动理论认为，光是一定频率范围内的电磁波，其运动规律可用Maxwell方程组来描述： (1-1)其中，上式中为电场强度，为电位移，为磁场强度，为磁感应强度，一般情况下他们都是矢量且为时间空间坐标的函数，还满足物质方程： (1-2)式中为电极化强度，为电流密度，为自由电荷密度，为电导率，为磁化强度。在线性极化情况下式中为介质的线性极化率。在非磁，各向同性均匀介质中， ，在区域中，由(1-1)的第二式、(1-2)中第一式，有，将(1-2)第二式代入(1-1)第一式，等式两边取旋度，有 (1-3)由(1-1)第三式、(1-2)第一、三式可得 (1-4)将(1-4)代入(1-3)，由可得 (1-5)因为，(1-5)整理后可得 (1-6)对于无损介质（等效于）有 (1-7)式中为真空中的光速： (1-8)(1-6)、(1-7)为线性光学的基本方程。（二）赫姆霍茨方程激光光学中常用复数E（公式中用代替方便输入）表示电场强度： (2-1) (2-2)介质的电极化强度也可以用复数表达式： (2-3) (2-4) (2-5)式中带“”量为共轭量。利用(2-1)(2-5)式可将(1-7)式化为 (2-6)式中为复折射率 (2-7)在标量场假设下，(2-6)式成为 (2-8)在真空中，于是有 (2-9)(2-8)、(2-9)式都称为赫姆霍茨方程。（三）高斯光束表达式推导由前面分析可知稳态传输电磁场满足赫姆霍茨方程 (3-1)式中与电场强度的复表式间有关系： (3-2)由数理方程基本知识可知，平面波和球面波都是（3-1）式的特解。高斯光束则不同，它不是（3-1）式的精确解，而是在缓变振幅近似下的一个特解。设 (3-3)在SVA（缓变振幅）近似下有 (3-4)利用(3-4)式可将(3-1)式在柱坐标下写为 (3-5)在旋转对称情况下与无关，(3-5)式简化为如下的抛物方程 (3-6)为了求得(3-6)式的一个特解，可设在处有一振幅为 (3-7)的高斯光束，然后求在任意处的。式中为振幅常数，如果只考虑相对值，则可由归一化条件求出。定义为处场振幅减小到最大值的值，称为腰斑（或光腰，束腰），它是高斯光束光斑半径的最小值。设试探解为 (3-8)式中、为待定函数，满足 (3-9)将(3-8)式微分后代入(3-6)，整理得到 (3-10) 由(3-10)对任意成立条件得到下面两个关系式 (3-11)式中 表示。微分方程组(3-11)在边界条件为(3-9)式时的解为 (3-12)式中 (3-13)称为共交参数。于是我们证明了，形如 (3-14)的高斯光束是赫姆霍茨方程(3-1)在SVA近似下的一个特解。其物理意义为：如果在处有一形如(3-7)式的高斯光束，则它将以(3-14)式非均匀高斯球面波的形式在空间传播。(3-14)式可改写为 (3-15)式中 (3-16)利用(3-16)式可将改写为 (3-17) (3-18)三、问题研究：下面研究高斯光束在非线性薄膜介质中的折射系数透射系数的计算问题。计算中物理量取常用单位。如上图所示，为了使问题简单化，我们假设高斯光束垂直入射介质，非线性薄膜介质绝缘，面积无穷大，厚度为h。折射率表达式为： (4-1)式中为入射光光强。在波动光学中，光强为振幅的平方。、为常数。实际上值还与有关，但由于待研究的非线性介质薄膜厚度极小，简化问题，我们默认非线性薄膜介质垂直方向值不随的改变而改变。我们将问题研究分为三个：1、高斯光束在非线性薄膜上表面发生的反射透射2、进入介质后光束的传播3、光束在薄膜下表面再次发生反射透射问题1和3道理相同，我们暂时先不考虑问题3。主要研究问题1，问题2，而关键在于问题2光波在非线性介质中的传播。（一）光束在介质上表面反射透射光波是波长很短的电磁波，因此光的反射、折射现象就是电磁波在不同介质分界面上的反射、折射。任何波动在两个不同介质分层面上的反射、折射都属边值问题，因此电磁波在两种不同介质分界面上的反射与折射是由电磁场和在分界面上所满足的边值关系确定的。边值关系： (4-2)式中表示突变面（即分界面）的上的单位外法线。如图，在处平面为分界面的两种不同的均匀各向同性介质中，由(3-16)可知，入射波i，反射波r，折射波d的电场强度可表示为 (4-3)式中，、未知。因为在绝缘介质分层面上电荷面密度，电流面密度，因此在分界面上只需满足如下边值关系 (4-4)式中的下标t代表切向分量。综上，要满足(4-4)条件，首先要求入射波，反射波，折射波的相位在平面上任何一点，任何时刻都要相同，即 (4-5)由此必有，说明反射波，折射波与入射波频率相等。接下来继续分析反射波，折射波的振幅。由于，入射的高斯光束表达式为 (4-6)由此可得薄膜折射率为对r进行泰勒展开后取前两项，有即其中 正入射情况下，有透射波d： （二）进入介质后光束的传播当折射率为上述定义表示时，赫姆霍茨方程有其中整理后，有（若考虑吸收，吸收系数定义为设复折射率，可将方程改