高斯光束研究(可编辑修改word版)

高斯光束通过非线性介质的自聚焦现象摘要：随着信息技术和纳米技术的迅速发展，要求光信息存储器件中的最小信息位尺寸、大规模集成电路和微电子技术中的光刻线宽和光学显微镜的分辨率等均能达到纳米量级（<100nm），而由于光衍射本身的限制，无法达到实际需求。

非线性薄膜材料的研究，通过选择非线性强的光学薄膜材料，调节激光能量和控制薄膜厚度及结构，在非线性薄膜结构的出射面能使光斑尺寸进一步下降，实现纳米光斑。

该光斑通过近场耦合作用在信息存储薄膜或光刻薄膜上，从而实现纳米信息存储、纳米光刻或纳米成像。

本文主要研究高斯激光束通过非线性均匀绝缘介质后光强的改变。

由电磁场基本原理，推导出高斯光束是缓变振幅条件下波动方程的近似解，研究其在介质突变面处的反射透射。

重点研究高斯激光束在非线性介质中的传播问题，这一过程中有自聚焦现象。

研究过程主要采用数值计算方法用差分方程代替偏微分方程研究问题的数值解。

比较光强的变化。

关键词：高斯光束，非线性，自聚焦，差分方程一、引言随着信息技术和纳米技术的迅速发展，要求光信息存储器件中的最小信息位尺寸、大规模集成电路和微电子技术中的光刻线宽和光学显微镜的分辨率等均能达到纳米量级（<100nm ），而由于光衍射本身的限制，无法达到实际需求。

而通过非线性薄膜材料的研究，通过选择非线性强的光学薄膜材料，调节激光能量和控制薄膜厚度及结构，在非线性薄膜结构的出射面能使光斑尺寸进一步下降，实现纳米光斑。

该光斑通过近场耦合作用在信息存储薄膜或光刻薄膜上，从而实现纳米信息存储、纳米光刻或纳米成像。

实验中我们常常采用高斯光束作为光源进行问题研究。

高斯光束是波动方程在缓变振幅下的一个特解，非线性介质的折射率随光强的变化而变化，因而高斯光束通过非线性介质发生自聚焦和衍射现象，从而改变能量分布。

本文主要研究光强的变化， 通过具体数值建立数学模型，采用差分方程代替偏微分方程以求得问题的数值解，研究光束通过非线性介质后能量的变化。

二、预备知识 （一）波动方程波动理论认为，光是一定频率范围内的电磁波，其运动规律可用 Maxwell 方程组来描述：⎧ ∂⎪∇⨯ E = - B ⎪∂t ⎪∇ ⋅ D = ⎨∂ ⎪∇⨯ H = J + D (1-1)⎪ ∂t ⎪ ∇ ⋅ ⎩⎪B = 0其中， 上式中 为电场强度， 为电位移， 为磁场强度，为磁感应强度，一般情况下E D H B他们都是矢量且为时间空间坐标的函数，还满足物质方程：⎪D = 0 E + PM 0⎧⎨B = 0 (H + M )⎪⎪⎩J =E(1-2)式中 P 为电极化强度， J 为电流密度， 为自由电荷密度，为电导率， M 为磁化强度。

= 8.854 ⨯10-14AS / Vcm 真空中的介电常数 = 1.257 ⨯10-8VS / Acm 真空中的磁导率在线性极化情况下式中为介质的线性极化率。

P = 0 E在非磁，各向同性均匀介质中，= 0 ，在区域= 0j = 0 中，由(1-1)的第二式、(1-2)中第一式，有∇ ⋅ E ＝0 ，将(1-2)第二式代入(1-1)第一式，等式两边取旋度， 有∇⨯ ∇⨯∂(E )= -0 ∂t∇⨯ H(1-3)由(1-1)第三式、(1-2)第一、三式可得∂ ∂ ∇⨯ H =E +0E + P ∂t ∂t(1-4)将(1-4)代入(1-3)，由∇ ⨯（∇ ⨯ E ）=∇ (∇ ⋅ E ) - ∆E 可得∂∂2∂2∇(∇ ⋅ E )- ∆E = -E- E - P因为∇ ⋅ E ＝０，(1-5)整理后可得∂t0 0 ∂t 2 0 ∂t 2(1-5)∂∂2∂2 ∆E- E -E- P =0 0∂t对于无损介质（等效于= 0 ）有0 0 ∂t 2 0 ∂t 2(1-6)E - E 0P2式中c 为真空中的光速：∆ 1 ∂2 - c 2 ∂t 2 1 c 2 ∂2∂t 2=0(1-7)c 2 =(1-6)、(1-7)为线性光学的基本方程。

100(1-8)（二）赫姆霍茨方程激光光学中常用复数 E （公式中用 E 代替方便输入）表示电场强度： E =1*E + E 2(2-1)E ( x , y , z , t ) = E ( x , y , z ) e 介质的电极化强度也可以用复数表达式：it(2-2)1*P =(+ )(2-3)( x , y , z , t ) =()E x , y , z , t(2-4)( x , y , z ) =( )E x , y , z式中带“ * ”量为共轭量。

利用(2-1)—(2-5)式可将(1-7)式化为(2-5)∆ 2 2式中 为复折射率E ( x , y , z ) +k E ( x , y , z ) = 0(2-6)在标量场假设下，(2-6)式成为= 1+(2-7)()∆E ( x , y , z ) + 2k 2E ( x , y , z ) = 0在真空中， = 1，于是有∆E ( x , y , z ) + k 2E ( x , y , z ) = 0(2-8)、(2-9)式都称为赫姆霍茨方程。

（三）高斯光束表达式推导由前面分析可知稳态传输电磁场满足赫姆霍茨方程∆E ( x , y , z ) + k 2E ( x , y , z ) = 0式中 E ( x , y , z ) 与电场强度的复表式 E ( x , y , z , t ) 间有关系：E ( x , y , z , t ) = E ( x , y , z )e i t(2-8)(2-9)(3-1)(3-2)由数理方程基本知识可知，平面波和球面波都是（3-1）式的特解。

高斯光束则不同，它不是（3-1）式的精确解，而是在缓变振幅近似下的一个特解。

设E (r , z ) = 在 SVA （缓变振幅）近似下有∂AA (r , z )e -ikz(3-3)∂z kA ∂2 A ∂z 2k ∂A ∂z(3-4)利用(3-4)式可将(3-1)式在柱坐标(r ,, z ) 下写为∂2 A + 1 ∂A + 1∂2 A - ∂A =∂r 2 r ∂r r 2 ∂2 2ik ∂z0 (3-5)在旋转对称情况下 A 与无关，(3-5)式简化为如下的抛物方程∂2 A + 1 ∂A - ∂A =∂r 2 r ∂r 2ik ∂z0 (3-6)为了求得(3-6)式的一个特解，可设在 z = 0 处有一振幅为w 2 r ⎢⎥ ⎢ ⎩2A (r , z )z =0 = A (r , 0) = r 2 A 0e 0(3-7)的高斯光束，然后求在任意处 z 的 A (r , z ) 。

式中 A 0 为振幅常数，如果只考虑相对值， 则可由归一化条件求出。

定义为 z = 0 处场振幅减小到最大值 1的 r 值，称为腰斑e（或光腰，束腰），它是高斯光束光斑半径的最小值。

设试探解为2- f 2 ( z ) 2A (r , z ) = A 0 f 1 ( z ) e(3-8)式中 f 1 ( z ) 、 f 2 ( z ) 为待定函数，满足f 1 (0) = 将(3-8)式微分后代入(3-6)，整理得到f 2 (0) = 1(3-9)⎡ 2 f 2 ( z ) f ( z ) ikf '( z ) f ( z ) ⎤ ⎡ 2 f ( z ) f ( z ) ' ⎤2 1 + 2 1 4 2 r 2 - 1 2 2 + ikf 1 ( z )⎥ = 0 ⎣ 0 0 ⎦ ⎣ 0 ⎦ (3-由(3-10)对任意 r 成立条件得到下面两个关系式 10)⎧ 2 f 2 ( z ) ' ⎪ 2+ ikf 2 ⎪ 0 ⎨( z ) = 0 ⎪ 2 f 1 ( z ) f 2 ( z ) + ikf ' ( z ) = 0⎪ 2 1 0d(3-11)式中 ' 表示 dz。

微分方程组(3-11)在边界条件为(3-9)式时的解为f ( z ) = f ( z ) = 1 = 11 21- i 2z k 21- i z Z0 0(3-12)式中-( )Z Z称为共交参数。

于是我们证明了，形如Z 0 = 1 k222=r 22 -0z(3-13)A (r , z ) = A 0 1- i z1-i eZ 0(3-14)的高斯光束是赫姆霍茨方程(3-1)在 SVA 近似下的一个特解。

其物理意义为：如果在z = 0 处有一形如(3-7)式的高斯光束，则它将以(3-14)式非均匀高斯球面波的形式在空间传播。

(3-14)式可改写为- r 2 -⎡ kr 2 ⎤A2( z )i ⎢ 2 R ( z ) -⎥A (r , z ) = 0 0 ez 式中⋅ e⎢⎣ ⎥⎦(3-15)⎛ z ⎫2( z ) = 01+ ⎪ ⎝ 0 ⎭ — — 高斯光束的光斑半径R ( z ) = z ⎛ z + Z 0 ⎫ — —高斯光束的等相面曲率半径Z z ⎪ ⎝ 0 = tan -1 z — — Z 0⎭高斯光束的相位因子(3-16)利用(3-16)式可将 E (r , z ) 改写为r 2-⎧⎪ ⎡ r 2⎤ ⎫⎪ A - 2i ⎨k ⎢ 2 R ( z ) + z ⎥-⎬ E (r , z ) = 0 0 e ( z ) ⋅ e ⎩⎪ ⎢⎣ ⎦⎥⎪⎭ ( z )相位部分振幅部分(3-17)r 2- ⎧⎪ ⎡ r 2 ⎤ ⎫⎪A - 2 i ⎨k ⎢ 2 R ( z ) + z ⎥--t ⎬E (r , z , t ) = 0 0 e ( z )⋅ e⎩⎪ ⎢⎣ ⎦⎥ ⎪⎭( z )(3-18)三、问题研究：相位部分振幅部分下面研究高斯光束在非线性薄膜介质中的折射系数透射系数的计算问题。

计算中物理量取常用单位。

如上图所示，为了使问题简单化，我们假设高斯光束垂直入射介质，非线性薄膜 介质绝缘，面积无穷大，厚度为 h 。

折射率表达式为：(r , I ) =0 +I I式中 I 为入射光光强。

在波动光学中，光强为振幅的平方。

0 、I 为常数。

(4-1)实际上值还与 z 有关，但由于待研究的非线性介质薄膜厚度极小，简化问题，我们默认非线性薄膜介质垂直方向值不随 z 的改变而改变。

我们将问题研究分为两个：⎧⎩E ⎪ 1、高斯光束在非线性薄膜上表面发生的反射透射 2、进入介质后光束的传播（一）光束在介质上表面反射透射光波是波长很短的电磁波，因此光的反射、折射现象就是电磁波在不同介质分界面上的反射、折射。

任何波动在两个不同介质分层面上的反射、折射都属边值问题， 因此电磁波在两种不同介质分界面上的反射与折射是由电磁场 E 和 B 在分界面上所满足的边值关系确定的。

边值关系：⎪ n ⋅ (D 2 - D 1)=⎪n ⨯( E 2 - E 1 ) = 0 ⎨n ⋅ ( B - B ) = 0 (4-2)⎪ 2 1 ⎪n ⨯( H 2 - H 1 ) = 式中n 表示突变面（即分界面）的上的单位外法线。