激光器速率方程
第二章激光器的速率方程理论

dn WDn n dt
均匀加宽
W
c3 8
2
A21 g ( )
c3 1 W A21 2 2 2 (0 ) 2 (V / 2) 2
第 个模的光场与第 个原子作用的受激辐射速率为
2 c3 A21 W sin 2 k z 2 2 2 ( ) 2 2
速地转移到激光上能级E2,其跃迁几率用S32表示。 3. 处在E2能级的粒子,能通过自发辐射、非辐射跃迁和 受激辐射,跃迁到激光下能级E1,其跃迁几率分别用A21, S21和W21表示。
4. 处在E1能级的粒子,能通过受激吸收到达E2,或非辐
射跃迁到E0,其跃迁几率分别用 W12和S10表示。 各能级上的粒子数密度N0,N1,N2,N3如果变化? 光子数密度n如果变化?
在脉冲开始建立的时间内,光子数和反转粒子数为
n(t ) ni e t
( WDi )
D(t ) Di {1
Wni
[1 e t ]}
2.5 均匀加宽的激光器的多模振荡
纵模
q qc / 2 L
c q q 1 q 2L
谱线线型 激光器阈值 增益饱和 均匀加宽 多模振荡
n 0, dD 0 dt
( a b ) 1 D0 [(a b ) N0 ( a b )]/ 2 2
引入
1 || ( a b ) 2
dD || ( D D0 ) 2WDn dt
dd dt
dd dt
|| (d d0 ) 2Wd n
当激光器在阈值之上不太高时,激光光子数不太大,有
dn (G0 )n Cn 2 dt
激光原理-4.2 典型激光器的速率方程

f2 f1
B21
在辐射场 的作用下的总受激跃迁几率 W21 中,
分配在频率 处单位频带内的受激跃迁几率为:
W21 B21 B21 g% ,0
太原理工大学物理与光电工程学院
4、公式的修正
dn21 dt
sp
n2 A21
d
n2 A21g% , 0
d n2 A21
' '
'
此时有: '
激光器中的情形即 是如此!
0 '
太原理工大学物理与光电工程学院
g%
',
0
'
d
'
中的被积函数
只在辐射场中心频率 附近很窄范围内才不为零。
g% ',0 g% ,0
' ' 且:
因 ' 很小,则有:
'
0
' '
'
d
' 1
太原理工大学物理与光电工程学院
物理意义: 由于谱线加宽, 外来光的频率 并不 一定要精确等于原子发光的中心频率 0 才能产生 受激跃迁,而是主要在=0 附近的一个频率范围内 都能产生受激辐射。当ν偏离中心频率ν0时,跃 迁几率急剧下降。
太原理工大学物理与光电工程学院
6、受激辐射、受激吸收几率的其它表达形式
W21
B21 g% ,0
—吸收截面
中心频率处发射截面和吸收截面最大!
太原理工大学物理与光电工程学院
均匀加宽工作物质中心频率发射截面
21 0
A21 2
4
2
2 0
H
非均匀加宽工作物质中心频率发射截面
激光原理(4)-速率方程

= n2 ∫
+∞ −∞
A21 g (ν ,ν 0 )dν
= n2 A21
谱线加宽对单位时间内自发辐射跃迁的原子数没有影响 3. n2 数 个原子中单位时间内发生受激辐射跃迁的原子总
+∞ dn21 ( ) st = ∫ n2W21 (ν )dν −∞ dt
2. 线宽
1 = = g (ν 1 ,ν 0 ) g (ν 2 ,ν 0 ) g (ν 0 ,ν 0 ) 2
∆ν = ν 2 − ν 1
线宽:光谱线的宽度
FWHM = Full width at half maximum 半幅线宽
NJUPT
谱线加宽的机理
自然加宽(Natural broadening) 这种谱线加宽是不可避免的 (1) 经典理论 处于激发态的发光粒子,在自发辐射的发光过程中, 辐射功率不断衰减,导致光谱线有一定宽度。 经典电子理论:原子是一个正电中心和 一个负电中心组成的偶极子 γt − 2 i 2πν 0 t 0
③ 全量子理论。本质上是量子电动力学体系,其特点是,将激光场看成是 遵循量子化规律的光子群的集合,将与激光场发生作用的工作物质看成是遵循 量子力学规律的微观粒子的集合,在此基础上进而将两者看成是一个统一的体 系而加以量子理论处理。这种理论体系的主要优点,是它能对涉及到激光与物 质相互作用过程中出现的各种现象与效应,给出严格而又全面的物理描述;其 不足之处,是这种理论的数学处理过程过于繁杂而不便求解。基于全量子理论, 在一定前提下还可派生出一些往往是十分简洁有用的专门理论。如在忽略量子 化激光场的位相特性(或光子数目起伏)的前提下,可简化为速率方程理论, 能非常方便地用它来描述激光的产生、振荡与放大等过程中的粒子数输运和激 光功率方面的动态特性。
电子科大龙格库塔法解半导体激光器速率方程

姓名:李清学号:************ 学院:电子科学技术研究院龙格库塔法解半导体激光器速率方程1、光强与载流子随时间变化曲线图1图22、分析半导体激光器工作原理从图1中我们可以看出激光器工作开始时反转粒子数不断增加,当超过阈值后发生激光的激射。
同时,观察图2我们还可以发现,当发生激射后,反转粒子数还在不断增加,激光光强不断增加。
由于激光的产生是以消耗反转粒子数为代价的,因此载流子数开始减少,小于阈值后便不会继续产生激光。
接着反转粒子数被不断激励,数目增加,超过阈值后又发生激光激射,这就是半导体激光器的工作原理。
3、使用稳态分析推导阈值电流的大小在稳态时,增益等于损耗,也就是G=r,同时电场和载流子数均不随时间变化,将这些带入第二个方程,即可解得结果如下,与之电流为0.058417067A。
4、源程序:t0=0;h=1e-12;tn=1e-8;n=(tn-t0)/h+1;E=zeros(1,n);N=zeros(1,n);E(1)=0.1;N(1)=1e8;t=t0:h:tn;for i=1:n-1E1=f1(N(i),E(i));E2=f1(N(i)+h/2,E(i)+E1*h/2);E3=f1(N(i)+h/2,E(i)+E2*h/2);E4=f1(N(i)+h,E(i)+E3*h);E(i+1)=E(i)+(E1+2*E2+2*E3+E4)*h/6;N1=f2(E(i),N(i));N2=f2(E(i)+h/2,N(i)+N1*h/2);N3=f2(E(i)+h/2,N(i)+N2*h/2);N4=f2(E(i)+h,N(i)+N3*h);N(i+1)=N(i)+(N1+2*N2+2*N3+N4)*h/6;endNn=N(n);In=abs(E(n)).*abs(E(n));N0=1.5e8;g=3.6e3;r=252e9;Nx=r/g+N0;q=1.6e-19;re=1.66e9;I0=re*Nx*q+r*q*E(n);subplot(211)plot(t,abs(E))title('电场强度曲线')xlabel('t')ylabel('E')subplot(212)plot(t,N)title('载流子数变化曲线')xlabel('t')ylabel('N')figure(2),plot(t,abs(E).*abs(E)) title('光强变化曲线')xlabel('t')ylabel('光强')function f1=f1(N,E)a=3;g=3.6e3;N0=1.5e8;G=g*(N-N0);r=252e9;f1=0.5*(1+1i*a)*(G-r)*E;function f2=f2(E,N)I=90e-3;q=1.6e-19;re=1.66e9;g=3.6e3;N0=1.5e8;G=g*(N-N0);f2=I/q-re*N-G*(abs(E))^2;。
matlab对半导体激光器速率方程进行求解

matlab对半导体激光器速率方程进行求解文章标题:深入探讨Matlab对半导体激光器速率方程的求解1. 简介半导体激光器作为一种重要的光电器件,在通信、医疗、材料加工等领域具有广泛的应用。
而速率方程是描述半导体激光器内部过程的重要数学模型,通过对速率方程的求解,可以更好地理解半导体激光器的工作原理和特性。
在本文中,我将结合Matlab软件,就如何利用Matlab对半导体激光器速率方程进行求解展开深入探讨。
2. 半导体激光器速率方程简介半导体激光器速率方程是描述半导体激光器内部电子和光子之间相互作用的重要数学模型。
其一般形式如下:\[ \frac{dn}{dt} = G - \frac{n}{\tau_{n}} - \frac{nI}{I_{s}} \]\[ \frac{dI}{dt} = \frac{e\eta V}{q}n - (\alpha+g)nI \]其中,n为激子数密度,I为激光光强,G为外界注入的光子数密度,τn为激子寿命,I s为饱和光子密度,η为激子与光子的电荷,V为激光器波导体积,q为电子电荷量,e为元电荷,α为损耗系数,g为增益系数。
3. Matlab在对半导体激光器速率方程求解中的应用Matlab作为一种强大的科学计算软件,提供了丰富的数学建模和仿真工具,非常适合用于对半导体激光器速率方程的求解。
利用Matlab,可以通过编写相应的数学模型和算法,实现对速率方程的数值求解。
Matlab提供了丰富的绘图和数据分析功能,可以对求解结果进行直观展示和分析。
4. 在Matlab中编写半导体激光器速率方程求解程序在使用Matlab对半导体激光器速率方程进行求解时,首先需要编写相应的数学模型和算法。
可以利用Matlab的ODE求解器对速率方程进行数值求解。
还可以结合Matlab的优化工具,对速率方程的参数进行拟合和优化,得到更准确的结果。
在编写程序时,应注意处理数值求解的收敛性和稳定性,避免出现数值不稳定或发散的情况。
激光器速率方程

思路小结
• 爱因斯坦采用唯象法得到光和物质相互作用 的关系式 • 考虑线型函数后必要的修正:几率按频率的 分布函数 • 原子和准单色光相互作用 • 单模振荡速率方程组(三能级系统和四能级 系统) • 多模振荡速率方程组
R1、R2为单位体积中,在单位时间内激励至E1 、E2能级的粒子数(激励速率);τ1、τ2为E1 、E2能级的寿命; τ21为E2能级由于至E1能级 跃迁造成的有限寿命。
• In Figure we show the ground state 0 as well as the two laser levels 2 and 1 of a four-level laser system. The density of atoms pumped per unit time into level 2 is taken as R2, and that pumped into 1 is R1. Pumping into 1 is, of course, undesirable since it leads to a reduction of the inversion. In many practical situations it cannot be avoided. The actual decay lifetime of atoms in level 2 at the absence of any radiation field is taken as t2. This decay rate has a contribution tspont which is due to spontaneous (photon emitting)2→1 transition as well as to additional non-radiative relaxation from 2 to 1. The lifetime of atoms in level 1 is t1.
激光原理-4.2 典型激光器的速率方程

太原理工大学物理与光电工程学院
1、三能级系统的能级跃迁特点和跃迁示意图
W13
A31
S31
E3 泵浦上能级
S32(热弛豫)
E2
激光上能级 (亚稳态)
A21
S21
W21
W12
S31, A31 S32; S31 A31
E1(激泵光浦下下能能级级)
S21 A21
太原理工大学物理与光电工程学院
F
太原理工大学物理与光电工程学院
思考:分别求洛仑兹线型和高斯线型下简
化线型函数对应的等效谱宽 。
21 , 0
A21 2
8
h
2 0
g% ,0
21 l , 0 Nl 21 Nl 21 N
l
l
太原理工大学物理与光电工程学院
根据简化模型, 四能级多模速率方程
dn3 dt
n0W03
' '
'
此时有: '
激光器中的情形即 是如此!
0 '
太原理工大学物理与光电工程学院
g%
',
0
'
d
'
中的被积函数
只在辐射场中心频率 附近很窄范围内才不为零。
g% ',0 g% ,0
' ' 且:
因 ' 很小,则有:
'
0
' '
'
d
' 1
太原理工大学物理与光电工程学院
n2
f2 f1
n1
21
,
0
Nl
R1,R2为单位体积中,在单位时间内激励至 E1,E2能级的粒子数;τ1,τ2为E1, E2能级的寿命; τ21为E2能级由于至E1能级的跃迁造成的有限寿命。
第四章典型激光器 的速率方程

dDn Dn Dn 21 n l ,n 0 )vNl n0W03 dt 2
2
1 A21 S21
dDn Dn ) Dn 21 n l ,n 0 vN n0W03 dt 2 n0W03 2 n0 w03 2 Dn 0 Dn 21 n 1 ,n 0 ) 2 In1 In 1 1 21 n 1 ,n 0 )vN 2 1 1 hn 0 I s n 1 ) dDn0 Dn0 1 I s n1 ) n0W03 In1 Nh n1v n 0 , n 0 , 3 1 dt 2
S10
dn0 n1 S10 n0W03 n3 A30 dt
dNl Nl Nl f2 n2W21 n1W12 (n2 n1 ) 21 n ,n 0 )vN l dt Rl f1 Rl
n0 n1 n2 n3 n
忽略n3W30 , n2A21 ?
• 小信号增益曲线的形状完全取决于谱线线型函数 均匀加宽介质
中心频率处小 信号增益系数
2 ) D n 2 0 0 H n ) g H n 0 ) gH 2 n n 0 ) Dn H 2)2
g n 0 ) Dn 21
0 H 0
v A21 Dn 2 4 2n 0 Dn H
ln 2
g n)
Dn n )
• 增益线宽~ (自发辐射)荧光线宽DnF 氦氖 Nd:YAG 钕玻璃 若丹明 6G GaAlAs (0.85mm) InGaAsP (1.55mm)
荧光线宽(s-1) 1.5×109 1.95×1011 7.5×1012 5×1012~3×1013 1013 1012~1013
1 2 n
N l N l 1 ,N l 2 N l n
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
dn3 = n1W13 − n3 ( S32 + A31 ) dt dn2 f2 = −(n2 − n1 )σ 21 (ν ,ν 0 )υN l − n2 ( S21 + A21 ) + n3S32 dt f1 n1 + n2 + n3 = n dN l f2 Nl = (n2 − n1 )σ 21 (ν ,ν 0 )υN l − dt f1 τ Rl
W21 = σ 21 (ν ,ν 0 )υN l W12 = σ 12 (ν ,ν 0 )υN l
υ为工作物质中的光速
4 发射截面和吸收截面
• σ21(ν,ν0)和σ12(ν,ν0)分别称为发射截面和吸收截面, 它们具有面积的量纲 A21υ 2 ~ 中心频率处的发射截面与吸 g (ν ,ν 0 ) σ 21 (ν ,ν 0 ) = 2 8πν 0 收截面最大。当ν=ν0时,均
dt = ( n2 − f1 n1 )σ 21 (ν ,ν 0 )υN l −
τ Rl
• 1)假设各个模式的衍射损耗比腔内工作物质 的损耗及反射镜损耗小很多,因而可以认为 各个模式的损耗是相同的。 ~ ~ (ν ,ν ) g ′(ν ,ν 0 ) 代 • 2)将线型函数 g 0 用一矩形谱线 替,并使矩形谱线的高度与谱线轮廓中心点 的高度相等,矩形谱线所包含的面积与原有 ~ ~ 谱线包含的面积相等。即 g ′(ν ,ν 0 ) = g (ν 0 ,ν 0 )
1 δν = ~ g (ν 0 ,ν 0 )
• 对洛仑兹线型与高斯线型,等效线宽分别为
δν = π
2 ∆ν F , δν = 1 π ∆ν F 2 ln 2
• 按照以上简化模型,四能级多模振荡的速率 方程可写为(见下页),式中N为各模式光 子数密度的总和;σ21为中心频率处的发射截 面;η1为E3能级向E2能级无辐射跃迁的量子 效率; η2为E2能级向E1能级跃迁的荧光效率
dN f2 N = (n2 − n1 )σ 21υN − dt f1 τR dn3 S32 = n0W03 − n3 dt η1 dn2 f2 A21 = −(n2 − n1 )σ 21υN − n2 + n3 S32 dt f1 η2 dn0 = n1S10 − n0W03 dt n0 + n1 + n2 + n3 = n
∫
+∞
−∞
ρν′ dν ′ = ∫ ρδ (ν ′ −ν )dν ′ =ρ
−∞
+∞
单色光辐射场的总 能量密度,Jm-3
+∞ dn21 ~ ~ ( ) st = n2 B21 ∫ g (ν ′,ν 0 ) ρδ (ν ′ − ν )dν ′ = n2 B21 g (ν ,ν 0 ) ρ −∞ dt
• 在频率为ν的单色辐射场的作用下,受激跃Leabharlann 几率为• 四能级系统速率方程组
dn3 = n0W03 − n3 ( S32 + A30 ) dt dn2 f2 = −(n2 − n1 )σ 21 (ν ,ν 0 )υN l − n2 ( S 21 + A21 ) + n3 S32 dt f1 dn0 = n1S10 − n0W03 + n3 A30 dt n0 + n1 + n2 + n3 = n dN l Nl f2 = (n2 − n1 )σ 21 (ν ,ν 0 )υN l − dt f1 τ Rl
1 爱因斯坦采用唯象法得到光和物质 相互作用的关系式
dn21 ( ) sp = A21n2 dt dn21 ( ) st = W21n2 , W21 = B21 ρν dt dn12 ( ) st = W12 n1 , W12 = B12 ρν dt A21 8πhν 3 = , B12 f1 = B21 f 2 3 B21 c
dn2 f2 = − ∑ ( n2 − n1 )σ 21 (ν l ,ν 0 )υN l − n2 ( S21 + A21 ) + n3S32 dt f1 l
•
~ g (ν l ,ν 0 ) 值不同,必须建 由于每个模式的频率、损耗、
立m个光子数密度速率方程,其中第l个模的光子数 密度速率方程为 dN l f2 Nl
• 四能级系统另外一种粒子数密度速率方程
dn2 n2 f2 = R2 − − ( n2 − n1 )σ 21 (ν ,ν 0 )υN l dt f1 τ2 dn1 n n f = R1 − 1 + 2 + ( n2 − 2 n1 )σ 21 (ν ,ν 0 )υN l dt τ 1 τ 21 f1
S32 η1 = S32 + A30
η2 =
A21 A21 + S 21
7 总量子效率
• η F = η1η2 为总量子效率,它的意义:由光泵 抽运到E3能级的粒子,只有一部分通过无辐 射跃迁到达激光上能级E2,另一部分通过其 它途径返回基态。而到达E2能级的粒子,也 只有一部分通过自发辐射跃迁到达E1能级并 发射荧光,其余粒子通过无辐射跃迁而跃 迁到E1能级。
f 2 A21υ σ 12 (ν ,ν 0 ) = 2 f1 8πν 0
2
匀加宽物质和非均匀加宽物 ~ g (ν ,ν 0 ) 质的发射截面分别为
ln 2υ 2 A21 υ 2 A21 , σ 21 = 3 2 2 σ 21 = 2 2 4π ν 0 ∆ν H 4π ν 0 ∆ν D
~ ~ ~ A21 g (ν ,ν 0 ) A21 g (ν ,ν 0 ) A21 g (ν ,ν 0 ) 8πν 2 W21 = Nl = N lV = nl nν = 3 nν nν V nν V c ~ W21 A21 g (ν ,ν 0 ) = al = nl为腔内第l模内的总光子数 nl nν V
~ ~ W21 = B21 g (ν ,ν 0 ) ρ ,W12 = B12 g (ν ,ν 0 ) ρ
由于谱线加宽,和原子相互作用的单色光的频率ν并 不一定要精确等于原子发光的中心频率ν0才能产生 受激跃迁,而是在ν= ν0附近一个频率范围内都能产 生受激跃迁。 激光器内ρ与第l模内的光子 数密度Nl的关系为ρ= Nlhν
发射荧光的光子数 ηF = 工作物质从光泵吸收的光子数
思路小结
• 爱因斯坦采用唯象法得到光和物质相互作用 的关系式 • 考虑线型函数后必要的修正:几率按频率的 分布函数 • 原子和准单色光相互作用 • 单模振荡速率方程组(三能级系统和四能级 系统) • 多模振荡速率方程组
• 采用激励速率和能级寿命来描述粒子数变化 速率而不涉及具体的激励及跃迁过程 • 前面的速率方程忽略了激光下能级的激励过 程,对大部分激光工作物质来说,这一忽略 是允许的。 • 根据所研究工作物质的激励与跃迁过程选择 或建立适用的速率方程
6 多模振荡速率方程
• 如果激光器中有m个振荡模,其中第l个模的频 率、光子数密度、光子寿命分别为νl、Nl及τRl 。 则E2能级的粒子数密度速率方程为
nν:腔内单位体积中频率处于ν附近单位频率间隔内 ν 的光波模式数
得到:一个模式内的一个光子引起的受激跃迁 几率等于分配到同一模式上的自发跃迁几率。
f2 W21 = al nl , W12 = al nl f1
• 对W21作出近似计算 • 设谱线的总自发辐射跃迁几率为A21,谱线 宽度为∆ν,并假设A21均匀分配在∆ν所包含 的所有模式上,则分配在一个模式上的自发 辐射跃迁几率为 A21
对表达式进行修正
+∞ dn21 ( ) sp = ∫ n2 A21 (ν )dν = n2 A21 −∞ dt +∞ +∞ dn21 ~ ( ) st = ∫ n2W21 (ν )dν = n2 B21 ∫ g (ν ,ν 0 ) ρν dν −∞ −∞ dt
该积分与辐射场ρν的带宽∆ν′有关。 1 原子和连续光辐射场的相互作用,∆ν ′ >> ∆ν 2 原子和准单色光辐射场相互作用,∆ν ′ << ∆ν
−∞
+∞
A21(ν)表示在总自发跃迁几率A21中,分配在频率ν处单 位频率内的自发跃迁几率; W21(ν)表示在总受激跃迁几 率W21中,分配在频率ν处单位频率内的受激跃迁几率
c3 c3 A21 (ν ) B21 = A21 = ~ 8πhν 3 8πhν 3 g (ν ,ν 0 )
c3 ~ B21 (ν ) = B21 g (ν ,ν 0 ) = A (ν ) 3 21 8πhν ~ W21 (ν ) = B21 (ν ) ρν = B21 g (ν ,ν 0 ) ρν
2 考虑线型函数后必要的修正
• 线型函数可以理解为几率按频率的分布函数 ~ ~ P(ν ) = Pg (ν ,ν ) = n hν A g (ν ,ν ) = n hν A (ν )
0 2 0 21 0 2 0 21
~ A21 (ν ) = A21 g (ν ,ν 0 )
∫
+∞
−∞
~ A21 (ν )dν = ∫ A21 g (ν ,ν 0 )dν = A21
al =
A21nl A21 W21 = = Nl nνV∆ν nν ∆ν
nν V∆ν
f2 f 2 A21nl f 2 A21 N l W12 = W21 = = f1 f1 nνV∆ν f1 nν ∆ν
5 单模振荡速率方程组
• 三能级系统速率方程组:各能级集居数随时间变 化的方程和激光器腔内的光子数密度随时间变化 的规律
n为单位体积 工作物质内的 总粒子数,第 l个模式的光 子寿命为τRl, 工作物质长度 l等于腔长L。
dn3 = n1W13 − n3 ( S32 + A31 ) dt dn2 = n1W12 − n2W21 − n2 ( S 21 + A21 ) + n3 S32 dt n1 + n2 + n3 = n dN l Nl = n2W21 − n1W12 − dt τ Rl