高斯光束研究
高斯光束通过非线性介质的自聚焦现象 The Self-focusing phenomenon of Gaussian beams passing through a nonlinear medium
姓 名: 王睿 学 号: [1**********] 班 级: 2005级4班 指导教师: 刘永明 职 称: 教授
2009年5月完成
目 录
摘要……………………………………………………………………… 3 Abstract ………………………………………………………………… 3 一、引言………………………………………………………………… 4 二、预备知识…………………………………………………………… 4
(一)波动方程 …………………………………………………… 4 (二)赫姆霍茨方程……………………………………………… 6 (三)高斯光束表达式推导………………………………………… 7 三、问题研究…………………………………………………………… 10
(一)光束在介质上表面折射透射……………………………… 11 (二)光束在介质中的传播……………………………………… 14 四、模型建立及分析…………………………………………………… 17 参考文献 ………………………………………………………………… 24 附录程序……………………………………………………………………25 致谢 …………………………………………………………………… 27
高斯光束通过非线性介质的自聚焦现象
摘要:随着信息技术和纳米技术的迅速发展,要求光信息存储器件中的最小信息位尺寸、大规模集成电路和微电子技术中的光刻线宽和光学显微镜的分辨率等均能达到纳米量级(
本文主要研究高斯激光束通过非线性均匀绝缘介质后光强的改变。由电磁场基本原理,推导出高斯光束是缓变振幅条件下波动方程的近似解,研究其在介质突变面处的反射透射。重点研究高斯激光束在非线性介质中的传播问题,这一过程中有自聚焦现象。研究过程主要采用数值计算方法用差分方程代替偏微分方程研究问题的数值解。比较光强的变化。
关键词:高斯光束,非线性,自聚焦,差分方程
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
一、引言
随着信息技术和纳米技术的迅速发展,要求光信息存储器件中的最小信息位尺寸、大规模集成电路和微电子技术中的光刻线宽和光学显微镜的分辨率等均能达到纳米量级(
实验中我们常常采用高斯光束作为光源进行问题研究。高斯光束是波动方程在缓变振幅下的一个特解,非线性介质的折射率随光强的变化而变化,因而高斯光束通过非线性介质发生自聚焦和衍射现象,从而改变能量分布。本文主要研究光强的变化,通过具体数值建立数学模型,采用差分方程代替偏微分方程以求得问题的数值解,研究光束通过非线性介质后能量的变化。
二、预备知识 (一)波动方程
波动理论认为,光是一定频率范围内的电磁波,其运动规律可用Maxwell 方程组来描述:
⎧∂B ∇⨯E =- ⎪ ∂t ⎪ ⎪⎪∇⋅D =ρ
⎨
⎪∇⨯H =J +∂D (1-1)⎪∂t
⎪⎪⎩∇⋅B =0
∂ ∂ ∂
+j +k 其中,∇≡i ∂x ∂y ∂z
上式中E 为电场强度,D 为电位移,H 为磁场强度,B 为磁感应强度,一般情
况下他们都是矢量且为时间空间坐标的函数,还满足物质方程:
⎧D =ε0E +P ⎪ ⎪
⎨B =μ0(H +M ) ⎪
(1-2)⎪⎩J =σE
式中P 为电极化强度,J 为电流密度,ρ为自由电荷密度,σ为电导率,M 为磁化强度。
ε0=8.854⨯10-14AS /Vcm 真空中的介电常数μ0=1.257⨯10VS /Acm 真空中的磁导率
在线性极化情况下
-8
P =χε0E
=0 j =0中,由(1-1)的第
式中χ为介质的线性极化率。
在非磁,各向同性均匀介质中,M =0 ,在区域ρ
二式、(1-2)中第一式,有∇⋅E =0,将(1-2)第二式代入(1-1)第一式,等式两边取旋度,有
∂
∇⨯∇⨯E =-μ0∇⨯H
∂t (1-3)
()
由(1-1)第三式、(1-2)第一、三式可得
∂E ∂P ∇⨯H =σE +ε0+
∂t ∂t (1-4)
(∇⨯E )=∇将(1-4)代入(1-3),由∇⨯
(∇⋅E )-∆E 可得
22 ∂E ∂E ∂P
∇∇⋅E -∆E =-μ0σ-μ0ε02-μ02
∂t ∂t ∂t (1-5)
()
因为∇⋅E =0,(1-5)整理后可得
22 ∂E ∂E ∂P
∆E -μ0σ-μ0ε02-μ02 =0
∂t ∂t ∂t (1-6)
对于无损介质(等效于σ=0)有
2 1∂E 1∂P
∆E -22-22 =0
c ∂t ε0c ∂t (1-7)
2
式中c 为真空中的光速:
c 2
=
1
μ0ε0
(1-6)、(1-7)为线性光学的基本方程。
(二)赫姆霍茨方程
激光光学中常用复数E (公式中用E 代替方便输入)表示电场强度:
E =1 2
(E + E *
)
E (x , y , z , t )= E
(x , y , z )e i ωt
介质的电极化强度也可以用复数表达式ϕ:
P =12(ϕ +ϕ
*)
ϕ
(x , y , z , t )=ε
0χE (x , y , z , t )
ϕ
(x , y , z )=ε
0χE (x , y , z )
式中带“*”量为共轭量。
利用(2-1)—(2-5)式可将(1-7)式化为
(1-8)
(2-1)
(2-2)
(2-3)
(2-4)
(2-5)
22 ∆E (x , y , z )+ηk E (x , y , z )=0
为复折射率
式中η
(2-6)
=η (2-7)
在标量场假设下,(2-6)式成为
2k 2E x , y , z =0∆E (x , y , z )+η()
在真空中,η
=1,于是有
(2-8)
∆E (x , y , z )+k 2E (x , y , z )=0
(2-8)、(2-9)式都称为赫姆霍茨方程。
(三)高斯光束表达式推导
由前面分析可知稳态传输电磁场满足赫姆霍茨方程
(2-9)
∆E (x , y , z )+k 2E (x , y , z )=0
式中E (x , y , z )与电场强度的复表式E (x , y , z , t )间有关系:
(3-1)
E (x , y , z , t )=E (x , y , z )e i ωt
(3-2)
由数理方程基本知识可知,平面波和球面波都是(3-1)式的特解。高斯光束则不同,它不是(3-1)式的精确解,而是在缓变振幅近似下的一个特解。设
E (r , z )=A (r , z )e -ikz
在SV A (缓变振幅)近似下有
(3-3)
∂A
kA ∂z
∂2A ∂A
k 2∂z ∂z (3-4)
利用(3-4)式可将(3-1)式在柱坐标
(r , ϕ, z )下写为
∂2A 1∂A 1∂2A ∂A
++2-2ik =0 22∂r r ∂r r ∂ϕ∂z (3-5)
在旋转对称情况下
A 与ϕ
无关,(3-5)式简化为如下的抛物方程
∂2A 1∂A ∂A
+-2ik =0 2∂r r ∂r ∂z (3-6)
为了求得(3-6)式的一个特解,可设在z =0处有一振幅为
r 2
2
w 0
A (r , z )z =0=A (r ,0)=A 0e
的高斯光束,然后求在任意处
-
(3-7)
z 的A (r , z )。式中A 0为振幅常数,如果只考虑相对
1
的r 值,称为e
值,则可由归一化条件求出。ω0定义为z =0处场振幅减小到最大值腰斑(或光腰,束腰),它是高斯光束光斑半径的最小值。
设试探解为
-f 2(z )
r 2
2
ω0
A (r , z )=A 0f 1(z )e
式中f 1(z )、f 2(z )为待定函数,满足
(3-8)
f 1(0)=f 2(0)=1
将(3-8)式微分后代入(3-6),整理得到
(3-9)
⎡2f 22(z )f 1(z )ikf 2'(z )f 1(z )⎤2⎡2f 1(z )f 2(z )⎤
'++ikf 1(z )⎥=0⎢⎥r -⎢422
ω0ω0ω0⎣⎦⎣⎦ (3-10)
由(3-10)对任意r 成立条件得到下面两个关系式
⎧2f 22(z )
+ikf 2'(z )=0⎪2
⎪ω0⎨
⎪2f 1(z )f 2(z )+ikf 'z =0
1()2⎪ω0⎩ (3-11)
d
式中 ' 表示。微分方程组(3-11)在边界条件为(3-9)式时的解为
dz
f 1(z )=f 2(z )=
11-i
2z 2k ω0
=
11-i
z Z 0
(3-12)
式中
2
12πω0
Z 0=k ω0=
2λ
(3-13)
称为共交参数。
于是我们证明了,形如
r 2-
2
ω0
A (r , z )=
A 01-i
Z 0
e
1-i
Z 0
(3-14)
的高斯光束是赫姆霍茨方程(3-1)在SV A 近似下的一个特解。其物理意义为:如果在
z =0处有一形如(3-7)式的高斯光束,则它将以(3-14)式非均匀高斯球面波的形式在空间传播。
(3-14)式可改写为
A 0ω0-ω2(z )
A (r , z )=e ⋅e
ωz 式中
r 2
⎡kr 2⎤-i ⎢-ψ⎥⎢2R z ⎦⎥⎣
(3-15)
ω(
z )=ω0
高斯光束的光斑半径
⎛z Z 0⎫
R (z )=z 0 +⎪——高斯光束的等相面曲率半径
⎝Z 0z ⎭ -1z ψ=tan ——高斯光束的相位因子Z 0
(3-16)
利用(3-16)式可将E (r , z )改写为
A 0ω0-ω2(z )
E (r , z )=e ⋅e ωz 相位部分
振幅部分
r 2
⎧⎫⎤⎪⎡r 2⎪-i ⎨k ⎢+z ⎥-ψ-ωt ⎬⎪⎢⎪⎣2R z ⎥⎦⎩⎭
r 2
⎧⎤⎫⎪⎡r 2⎪-i ⎨k ⎢+z ⎥-ψ⎬
⎢2R z ⎦⎥⎪⎪⎩⎣⎭
(3-17)
A 0ω0-ω2(z )
E (r , z , t )=e ⋅e ωz (3-18) 相位部分
振幅部分
三、问题研究:
下面研究高斯光束在非线性薄膜介质中的折射系数透射系数的计算问题。计算中物理量取常用单位。
如上图所示,为了使问题简单化,我们假设高斯光束垂直入射介质,非线性薄膜
介质绝缘,面积无穷大,厚度为h 。折射率η表达式为:
η(r , I )=η0+ηI I (4-1)
式中I 为入射光光强。在波动光学中,光强为振幅的平方。η0、ηI 为常数。
实际上η值还与z 有关,但由于待研究的非线性介质薄膜厚度极小,简化问题,我们默认非线性薄膜介质垂直方向η值不随z 的改变而改变。
我们将问题研究分为三个:
1、高斯光束在非线性薄膜上表面发生的反射透射 2、进入介质后光束的传播
3、光束在薄膜下表面再次发生反射透射
问题1和3道理相同,我们暂时先不考虑问题3。主要研究问题1,问题2,而关键在于问题2光波在非线性介质中的传播。
(一)光束在介质上表面反射透射
光波是波长很短的电磁波,因此光的反射、折射现象就是电磁波在不同介质分界面上的反射、折射。任何波动在两个不同介质分层面上的反射、折射都属边值问题,因此电磁波在两种不同介质分界面上的反射与折射是由电磁场E 和B 在分界面上所满足的边值关系确定的。
边值关系:
⎧n ⋅D 2-D 1=σ⎪
⎪n ⨯(E 2-E 1)=0⎨ (4-2) n ⋅B -B =0⎪(21)
⎪n ⨯(H -H )=π
21⎩
()
式中n 表示突变面(即分界面)的上的单位外法线。
如图,在z =0处平面为分界面的两种不同的均匀各向同性介质中,由(3-16)可知
R (0)→∞,ψ=0入射波i ,反射波r ,折射波d 的电场强度可表示为
⎧E i (r , t )=E 0⋅e i ωt ⎪i ω't
'E r , t =E ⋅e ⎨r ()0
(4-3)
⎪E (r , t )=E ''⋅e i ω''t
0⎩d
式中E 0=A 0e
-r 2
2
ω0
'、E 0''未知。 ,E 0
因为在绝缘介质分层面上电荷面密度σ=0,电流面密度π=0,因此在z =0分界面上只需满足如下边值关系
⎧⎪⎡⎣E d -(E i +E r )⎤⎦t =0
⎨ (4-4) H -H +H =0⎡⎤()i r ⎦t ⎪⎩⎣d
式中的下标t 代表切向分量。综上,要满足(4-4)条件,首先要求入射波,反射波,折射波的相位在z =0平面上任何一点,任何时刻都要相同,即
[i ωt =i ω't =i ω''t ]z =0 (4-5)
由此必有ω=ω'=ω'',说明反射波,折射波与入射波频率相等。
接下来继续分析反射波,折射波的振幅。 由于z =0,入射的高斯光束表达式为
-r 2
2
ω0
E i (r , t )=A 0e
由此可得薄膜折射率为
⋅e i ωt (4-6)
η=η0+ηI ⋅I 0e
对r 进行泰勒展开后取前两项,有
-
2r 2
2ω0
η≈η0+ηI ⋅I 0-
2ηI ⋅I 0
ω02
r 2
η≈n 0(1+cr 2)
其中
n 0=η0+ηI ⋅I 0
c =-2ηI ⋅I 0
(ηη
0+I ⋅I 0) ω20
正入射情况下,有透射波d :
E t )=2
d (r , η+1
E i (r , t )
(二)进入介质后光束的传播
当折射率η为上述定义表示时,赫姆霍茨方程有
∆E +k 2n 2
0(1+cr 2) 2E =0
其中
(1+cr 2) 2=1+2cr 2+c 2r 4≈1+2cr 2=1-Γ2r 2Γ2=-2c
整理后,有
∆E +k 2n 220(1-Γr 2) E =0
(若考虑吸收,吸收系数定义为
a (r , z ) =a 0(1+β(z ) r 2)
设复折射率η'=η-i
a
2k
,可将方程改为 ∆E +k 2
n '20(1-Γ'2r 2) E =0
n 'a 0
=-2c +i a 其中0=n 0-i 2k
,Γ'0
n k (β-c ) )
E (r , z ) =A (r , z ) e
在SAV 条件下,可得
-ikn 0z
∂2A 1∂A ∂A 222
+-2ikn -2k n 0ΓA =0 02∂r r ∂r ∂z
设该方程解具有高斯函数形式
A 0
A (r , z ) =e
s (z )
带入方程,对任意r 成立,满足微分方程组
ikn r 2-
2q (z )
⎧dq 22
=1+Γq ⎪⎪dz
⎨1ds 1⎪=⎪⎩s dz q
根据边界条件
q (z ) z =0=s (z ) z =0=q 1=iZ 0
有
1⎧
q 1cos Γz +sin Γz ⎪q =⎪
-q 1Γsin Γz +cos Γz ⎨⎪1
⎪s =sin Γz +q 1cos Γz ⎩Γ
式E 均用此式表示
根据上式,我们得到了高斯光束在薄膜中的传播方程,下面讨论过程中电场表达
1-η
通过折射率我们可知薄膜下表面反射系数为R 1=
1+η
2η
,透射系数为T 1=
1+η
当第一次下表面反射波返回上表面的时候,有上表面反射系数R 2=R 1, 透射系数
T 2=T 1,同时传播距离为2d ,由此分析我们可知,第一次下表面透射波
22η
E 1(r , d ) =⋅⋅E (r , d )
η+1η+1
第二次下表面透射波:
22
E 2(r , d ) =⋅R 1T 1E (r ,3d )
η+1
第N 次下表面透射波:
22n
E n (r , d ) =⋅R 1T 1E (r ,(2n +1) d )
η+1
薄膜透射波有
E d (r , d ) =lim ∑E i (r , d )
n →∞
i =1
n
实验中,由于折射率η的影响,几次后透射波会小到忽略不计,可以通过数值计算来求出合理的答案
需要带入实际参数来计算结果