第2章介质的非线性极化.docx
《第2章介质的非线性极化.docx》由会员分享,可在线阅读,更多相关《第2章介质的非线性极化.docx(19页珍藏版)》请在冰豆网上搜索。
第2章介质的非线性极化
第2章介质的非线性极化
本章主要问题:
光在介质中传播的波动方程有哪些不同形式?
介质极化率如何定义,有那些对称性质?
极化率实部和虚部有何物理意义,其间有何关系?
2.1非线性介质的波方程2.1.1非线性介质的麦克斯韦方程
光波在非性线介质中传播时也服从麦克斯韦方程:
„dD
'HJ
dt
'D二、
'B=0
D=0EP
B=%(HM)
J-E
式中E、D――电场强度、电感应强度
、B――磁场强度、磁感应强度
、M——电极化强度、磁极化强度
、」o——真空介电系数、真空磁导率
――电导率(代表介质的吸收损耗)
J――电流密度,
——自由电荷密度
iJ0
t
在非线性介质中,P可以展开为E的幕级数:
X(n)是n阶电极化率,它是个n+1阶张量。
将式(2.1.10)代入(2.1.5)可得
这里
(2.1.12)
是介质的线性介电系数;其中X11是线性极化率。
在各向异性介质中X11和&二者
都是复数二阶张量。
般非线性介质是绝缘体(J=0,?
=0)和非磁性材料(M=0),
则非线性介质的麦克斯韦方程组可表为:
'E=
(2.1.13)
ct
„cD
'HE
ct
(2.1.14)
D=&EPnl
(2.1.15)
2.1.2各向异性介质的时域波方程
将(2.1.13)的两边进行I运算,再将式(2.1.14)代入,并用式(2.1.15)得到
这就是描述光在各向异性非线性介质中传播的时阈波动方程。
该方程比线性波动方程仅多了右边的一项。
相当存在一个次波源。
第二项与介质的吸收损耗有关,若介质为无损耗的,即厂-0,再利用c=1/.'幕;0,式(2.1.16)表为
这是光在无损耗各向异性非线性介质中传播的时阈波动方程。
为解方程求得场强
E,必须首先求出非线性极化强度Pnl0
2.1.3各向异性非线性介质的频域波方程
将E(r,t)和Pnl(r,t)展开成i=1,2,3…个单色平面波的组合(傅里叶展开)
(2.1.18)
(2.1.19)
E(r,t)「Ei(二'Eei("iJ)
ii
Pnl(r,t)八Pnl(ki「i)八PiNLei(g心
i
(2.1.20)
(2.1.19)代入式(2.1.仃),消去两边的求和号和
2
八C)2£]E(k,J
SgC
这是各向异性非线性介质的单色平面波的波方程。
2.1.4各向同性非线性介质频域波方程
在方程(2.1.20)中,利用''E八CE)-'2E,考虑各向同性介质,有
'E=0;再用关系式k。
=d:
/c,k二kon和n二,;/;。
,则得
i2E(k,)k2E(k,)二Pnl(k,)(2.1.21)
这是各向同性非线性介质的单色平面波的波方程。
它是一个非齐次二阶微分方
程,难于求解,一般都要做近似简化处理,慢变振幅近似是一种常用的方法。
现在考虑一个沿z方向传播的稳态单色平面波,振幅随z变化,但不随时间
变化。
电场强度和非线性极化强度分别表为:
E(z,J=E(z)e匕t)
PNL(z,J=PNL(z)尹'」)
式中k和k'分别是原光波和极化波的波失。
将式(2.1.22)代入(2.1.21),其中式
(2.1.21)左边第一项为
-2
\、2E(z「)=(2i2kk2)E(z)ei(kz「)
_r.厶r
czcz
因此式(2.1.21)表为
此为在各向同性介质中z向传播的单色平面波的波方程。
假设在波长量级的距离内光波振幅的变化非常慢,满足以下条件:
并假设PNL(z「)随z的变化可以忽略不计,则式(2.1.23)中略去第一项写成
EE⑵ik0pNL(z◎)e-l(kz-Cj)二ik0pNL(z)ei(k'」)z
.z
2;°k'2;°k
式中:
k二k-k。
这样,在慢变近似条件下,各向同性非线性介质中z向传播的
单色波的频域波方程被简化为简单的一阶微分方程,便于求解。
式(2.1.25)描述
在稳态和在慢变近似条件下的各向同性非线性介质中沿z向传播的单色光波的
频域波方程。
若存在介质对光电场的吸收,根据式(2.1.16),式(2.1.25)应改写为
式中「-是介质的吸收系数
2.1.5各向同性非线性介质时域波方程
考虑各向同性介质'E二0及n二,(2.1.17)式变为
此为各向同性非线性介质中的时阈波方程。
设时域下的波场为单色平面波
(2.1.28)
E(z,t)=A(z,t)ei(kz-7
式(2.1.27)中的各项为
-2
V2E(z,t)=(鼻+i2k三-k2)A(z,t)ei(kMfl)
cz&
-2-2-
:
2
-2Pnl(z,t)三-JPnl(z,t)
-t
假设波的振幅随空间和时间皆缓慢变化,满足以下慢变近似条件:
则在(2.1.27)中略去场振幅的二阶时间导数和二阶空间导数,得到一阶波方程:
(2.1.30)
:
A^.1£A^=JPNL(z,t)e」z)
:
zv:
t2;0cn
这是在慢变近似条件下各向同性非线性介质中单色波的时域波方程。
若光波是一个宽脉冲,在(2.1.30)式中v二c/n是光波的相速度;
若光波是一个短脉冲,在(2.1.30)式中v=d•/dk是波包的群速度
2.2非线性极化率2.2.1极化强度的频域表达式
考虑电极化强度P与电场强度E之间的因果关系。
在时刻t,介质感应的电极化强度dP
(1)(t)是由在此之前时刻b=t-dt!
的电场强度E(tj在dt!
时间内的作用所确定,二者呈正比关系,
dP(t)二;。
x
(1)(t-切E(tjdti(221)
考虑E(ti)在t之前所有时间电场强度E(ti)对P⑴(t)的贡献,则有
P⑴(t)二"X⑴(t—ti)E(ti)dti(2.2.2)
实际上,当tit时,E(ti)对P⑴(t)没有贡献,⑴(t—ti)=O。
再取E(ti)和P(i)(t)的傅里叶变换
E(切=.EC)eJtid-.
(2.2.3)
P(i)(t)P()e」td■
将式(2.2.3)代入式(2.2.2),得到频域的表达式
P⑴匕)二;0x(i)CJE()(2.2.4)
式中
x(i)C'H:
:
X(i)(^ti)ei叫ti(2.2.5)
在非线性情况下,P可以展开为E的幕级数,极化强度在频域中表达为
P()=P⑴()P
(2)()P(3)()(2.2.6)
其中
P⑴(•)二;0X⑴C;•)E(■)
(2)
(2)
PC)=OX(;"'2):
E(-i)E(匕)
(2.2.7)
P⑶(■)二X⑶C;'i「'2JJE(「JE
(2)E(匕)
IIOH
(2.2.8)
『(■)=.;Xn)(t-ti,t-t2,…,t-tn)ei[ZFi、Z)]dtidt2…dtn
式中
F面给出各阶电极化强度的直角坐标分量表达式。
介质中的场由n个不同频
率的分量(包含着这些频率的谐波、和频波、差频波等)组成
(2.2.9)
E(t)八E(n)eqnt
n
式中ECn)是复数振幅,n可正可负。
并规定
频率为「的极化强度分量为
(2.2.11)
(2.2.12)
(■)八;o[(;■)E:
(‘)Ct
Pf(J勺凳(;★,2)E(i)E-(JoP
pl3(®M:
$率扁牡;用冷,3巳G(1E)02E)3()(2.2.13)
otPV
式中二i—2—3JII;川1二x,y,z。
2.2.1极化率的对称性
下面指出电极化率张量的对称特性,它反应了介质结构的对称性和电极化强度的实数性。
1.频率置换对称性
可以证明电极化率张量具有以下固有的置换对称性
(i)
(1)/、
x(')=X(—J
X?
*(8;耳,旳2)=孙(斜;虫2,国)=X?
(国2;-创3
(2.2.14)
IIIIH
澀几*(mggjlgn)=Xhngi’IlZn®)
二IHHI二Xhn丄('nJ',-'l」llL'nJ
若外场频率远离介质的共振频率,介质被认为是无色散的和无耗的,则存在着完
全的置换称特性,即式(2.2.14)中的*号可以取消。
有81个独立元。
如果介质具有中心对称结构,即在坐标反演变换{x,y,z}r{_x,—y,—z}时,P和E都要变成反方向。
由(2211)-(2213式可见,P⑴和P⑶的表示式不变,但P
(2)式左边变号,据对称性要求⑺必须等于零,该式才能成立。
也就是说,具
有中心对称介质的偶阶极化率为零。
若只考虑到三阶非线性效应,对于具有中心对称性的介质,没有二阶非线性效应,只有三阶非线性效应。
2.2.3简并因子
(1)若电场强度和电极化强度分别表示为
E(r,t)八E(n)eJntc.c.,(2.2.16)
n
P(r,t)='P(n)eJntc.c.。
(2.2.17)
n
考虑到极化率的对称性,频率为••的n阶极化强度分量表示如下,它是由n个波场所引起,其中有m个相同频率。
P丁(「)二D;0X叫…(;,1「2「3,…/'n)E-()EE(n)(2.2.18)
式中的系数D被称为简并因子,对于式(2.2.16)和(2.2.17)的情况,可以证明:
n!
D=-
m!
E(r,t)=1^[EC'n)e4ntE*^.)eint],
P(r力十.P(=十
对这种情况,极化强度分量式(2.2.18)也成立,但是简并因子变成
D才卫
Im!
几种常见的非线性光学效应的极化率表达式及其相应的两种简并因子列于
非线性过程
阶
极化率
Dm21』(n!
/m!
)
D=n!
/m!
线性色散
1
X(1g)
1
1
线性吸收
1
x⑴(恥)
1
1
电光效应
2
x⑴(g,0)
1
2
二次谐波
2
x⑵(2们;®,®)
1/2
1
和频效应
2
X(①3®Q2)
1
2
差频效应
2
X(国2,国3,1)
1
2
三次谐波
3
(3)
X(38;叱,8)
1/4
1
四波混频
3
X(国4;斜,矶离3)
3/2
6
简并四波混频
3
⑶/\
x(wn
3/4
3
简并四波混频相位共轭
3
(3)/\
x(国c;国1,—灼2,灼p)
3/2
6
光克尔效应(自作用)
3
(3)/、
x2、0,—0、0\
3/4
3
光克尔效应(互作用)
3
x⑶(国阿,七2)
3/2
6
自聚焦
3
(3)
分(时;时,_豹,豹)
3/4
3
饱和吸收
3
x⑶
3/4
3
双光子吸收
3
(3)/\
x(国1;叫卞2冋)
3/2
6
拉曼散射(斯托克斯)
3
(3)/、
»(%;备七1盟S)
3/2
6
拉曼散射(反斯托克斯)
3
x(%;◎,◎,Y*s)
3/4
3
注:
表中极化率括号中的分号之后为入射场频率,分号之前为生成场频率。
2.3Kramers-Kronig色散关系
2.3.1极化率实部与虚部的关系
必需指出,电极化率(巧是一个复数,若表达为
(.)=;;:
(.)「(.),(2.3.1)
其实部和虚部之间有如下关系
1旳7'(⑷')
"(•)Pd',(233)
富*/屮:
'-■■
式中P表示后面的积分为柯西主值积分。
这是著名的Kramers-Kronig色散关系,
简称K-K关系。
由K-K关系可见,只要知道极化率的实部和虚部中任何一个的光谱就可通过此关系求出另外一个。
根据,"(')是「’的奇函数,而'(■')是:
的偶函数。
K—K
关系可以写成如下另一种形式:
'(•)---P"(2')2'd'',(2.3.4)
兀©蛍、一尬
2・3・2极化率实部和虚部的物理意义
1.线性折射率和吸收系数与极化率的关系我们考察一束频率为■■的单色平面波在各向同性介质中沿z方向的传播所
产生的线性极化。
设光电场强度表示为
E(z,,)二E(z)ei(kz"c.c.,(2.3.6)
式中,k是非线性介质的复数波矢,其实部表示波的相位变化(介质的色散),
虚部表示波的振幅的变化(介质的吸收),即
k=k、ik"十n0i-,(2.3.7)
2
式中,k-二-是真空中的波矢;n-和:
-分别表示介质的线性折射率和线性吸
c
收系数。
由电感强度的定义,考虑远离共振情况下的线性极化效应,则有
D二;-EP⑴二;-E-
(1)E=(;-;-
(1))E=E,(2.3.8)
式中,;-为真空的介电系数;;二,-「-⑴为介质的复线性介电系数;⑴为介
质的复线性极化率,可以分为实部和虚部两部分,利用关系;、二;-(1*⑴、),则
;可表为
;=;0•;0⑴'•i;0⑴''=;'•i;0⑴''=;'(1•⑴'')。
(2.3.9)
z'
利用线性折射率n0,式(2.3.9)改为
2艾⑴”
;:
=n2;。
(1i—^)。
(2.3.10)
n。
再利用复线性折射率n=二厂;0和真空光速c=1/暫;0,将介质的复波矢表为
k=—n-■■、"0i。
(2.3.11)
c'
将式(2.3.10)代入(2.3.11),得到
1
(■/⑴''军
k^Jv~1+i—(2.3.12)
In0丿
式(2.3.12)中的根号中第二项的模远小于1,可将'二展成泰勒级数,近似取前
两项得
kU&门0(1i-
⑴''k
)=k°n0i:
⑴"
2n°2n°
(2.3.13)
将式(2.3.13)对比(2.3.7),利用
;'「(「⑴')得到
n°"1
(1)(
11
•■)']2L1-⑴(T,
2
(2.3.14)
:
.厂k°⑴()“_⑴(.丁。
(2.3.15)
n°cn°
可见介质的线性折射率和线性吸收系数分别与一阶极化率的虚部和实部成正比。
2.非线性折射率和吸收系数与极化率的关系
假设介质具有三阶非线性,入射激光是如式(2.3.6)的单色平面波,可以用
以下慢变近似非线性波方程(2.1.25),求解光场E(z)。
这里设厶k二k'-k=0。
购1PNL(z)二如PNL(z)。
(2.3.16)
.'z2;0cn02;0k
非线性极化强度表为
设e⑶C‘)=3⑶(J,⑶二匕)⑶一i(•)⑶'',
22
PNL(Z)=%[兀⑶⑥)’E(z)知於3)⑻)“E⑵]E⑵,
e3)(.JE(z)2ie3)(yE(z)2]E(z)
:
z2cn0
利用I
12
ocn°E(z),式(2.3.19)改为
2
-:
E(z)
:
z
i‘:
,■2⑶2
e()'Ii-
2cn°[I;oCno
e3)c)'=iko^^^Ii
IL2pen。
ocrb
e3)^)"I'e(z)
/(『IE(z)
22
ocno
k_ke3)c)'
kNL一ko
2
0cno
22;ocno
式(2.3.20)变为
-:
E(z)
.:
z
ikNLE(z)。
解得
E(z)=E(o)eikNLz。
由式(2.3.6)
E(乙J
=E(z)ei(kz")
二E(0)JkTz*
其中
△ot
ao
=konoi-=ni
22
e3)()'■e3)()"
-Ii—I
=k°noio2
2Socno
22
oCno
(2.3.18)
(2.3.19)
(2.3.20)
(2.3.21)
(2.3.22)
(2.3.23)
(2.3.24)
(2.3.25)
(2.3.26)
(2.3.27)
由式(2.3.26)和(2.3.27)可见,对三阶非线性介质,其非线性折射率由三阶极化率的实部决定,并与光强成正比;而非线性吸收系数由三阶极化率的虚部决定,也与光强成正比
2.3.3非线性折射率与非线性吸收系数的关系
将式(2326)和(2327)分别代入式(232)可以得到
•:
n(•)-P^-^d'。
(2.3.28)
2兀』_°%©'-国)
据式(234),厶n(J和厶:
(■)的关系还可表为
n(•)=£p^4d・'。
(2.3.29)
假若测得某非线性介质的线性吸收光谱和在强光作用下的吸收光谱,从两光
谱的差值算得■:
:
C-),就可以由式(2.3.29)算出该介质的非线性折射率谱
:
n()。
而且可由式(2.3.29)确定介质在某频率下的三阶极化率⑶'()o
单位制说明
本书采用国际单位制,即
两种单位制的关系为:
MKS/SI单位制,有许多文献米用cgs/esi单位制,
MKS/SI单位制
p(n)=;0(%・
cgs/esi单位制
p(n)”(n)En
线性极化率
⑴(SI)=4二⑴(esu)
二阶极化率
(2)(SI)4二4
(2)(esu)
3"04
三阶极化率
⑶(SI)=8⑶(esu)
')9“08
n阶极化率
(n)(SI)//n4(n)(esu),c=3108
(10c)
2.时间反演对称性
根据电极化强度的实数性可以证明
XU仰厝1厲2,…声n)=盘;人(七;卞1,七2,…,七n)(2.2.15)
3.空间结构对称性
由于介质结构的对称性,当笛卡儿坐标的指标被置换时,紀Lln保持不变,
使非线性极化率张量的独立矩阵元的总数大大减少:
⑵只有27个独立元;2(3)只