激光器速率方程
第二章激光器的速率方程理论
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
matlab对半导体激光器速率方程进行求解
半导体激光器速率方程求解半导体激光器是一种利用电流驱动产生激光的器件,广泛应用于通信、激光打印、医疗和科学研究等领域。
了解和理解半导体激光器的工作原理对于设计和优化激光器具有重要意义。
在本文中,我们将介绍半导体激光器速率方程以及如何使用MATLAB进行求解。
半导体激光器速率方程半导体激光器速率方程描述了半导体激光器内部发生的光与电子相互作用的过程。
速率方程包含了四个关键项:电子注入项I in,非辐射复合项R nr,自发辐射项R sp和受激辐射项R spon。
方程可以写成如下形式:dN dt =I ineV−Nτnr+Nτsp+Nτspon其中,N表示激光器内的载流子密度,t表示时间,e为元电荷,V为激光器的有效体积,τnr为非辐射寿命,τsp为自发辐射寿命,τspon为受激辐射寿命。
使用MATLAB求解速率方程MATLAB是一种功能强大的数值计算软件,适合用于求解微分方程。
我们可以利用MATLAB的ODE求解器来求解半导体激光器速率方程。
首先,我们需要定义速率方程的函数形式。
在MATLAB中,我们可以定义一个函数文件,例如rate_equation.m,其中包含速率方程的定义:function dNdt = rate_equation(t, N)% 参数定义I_in = 1; % 电子注入项V = 1; % 有效体积tau_nr = 1; % 非辐射寿命tau_sp = 2; % 自发辐射寿命tau_spon = 0.5; % 受激辐射寿命% 速率方程求解dNdt = (I_in / (e * V)) - (N / tau_nr) + (N / tau_sp) + (N / tau_sp on);end然后,我们可以使用MATLAB的ODE求解器对速率方程进行求解。
例如,可以使用ode45求解器来得到方程的数值解。
% 定义时间范围和初始载流子密度tspan = [010]; % 时间范围N0 = 0; % 初始载流子密度% 求解速率方程[t, N] = ode45(@rate_equation, tspan, N0);% 绘制载流子密度随时间变化的曲线plot(t, N);title('半导体激光器载流子密度随时间变化');xlabel('时间');ylabel('载流子密度');运行以上MATLAB代码,我们可以得到半导体激光器载流子密度随时间变化的曲线图。
第二节 典型激光器速率方程
l 为第l模式单程损耗
将W21 21 ( )vNl g2 W12 12 ( )vNl 21 ( )vNl 代入, 得 g1
单模振荡速率方程组
dn3 dt n1W13 n3 ( S32 A31 ) dn g2 2 n S ( n n1 ) 21 ( )vNl n2 ( A21 S 21 ) 3 32 2 g1 dt n n n n 1 2 3 dNl (n2 g 2 n1 ) 21 ( )vNl N l dt g1 Rl
量子效率
E3向E2无辐射跃迁的量子效率 : S32 1= S32 A30 E2向E1跃迁的荧光效率:
2=
A21 A21 S 21
总量子效率: 发射荧光的光子数 F=1 2= 工作物质从光泵吸收的 光子数
吸收截面和发射截面 2.发射截面
同理, 可定义上能级原子的发 射截面 21 ( ) W21 21 ( )vNl
A 21 又W21 g ( ) N l n
A 21 A 21v 2 21 ( ) g ( ) g ( ) 2 n v 8 0
其中 :
N N l 为各模式光子数密度的 总和;
单模振荡速率方程组 2. 四能级系统速率方程组(Nd:YAG,He-Ne激光器)
E3
W03 S30 A30 S32
S21
W A21 W21 12
激光原理第18讲 速率方程、小信号增益系数
第18讲 速率方程、小信号增益系数
18.1 单模振荡速率方程
• 1、三能级系统速率方程
– W13为抽运几率 – A31为自发辐射几率
E2
– S32、S31为热驰豫(无辐射跃迁)几率
S32
– S31<<S32、A31<<S32 – S32远大于S31和A31,基态E1上的粒子被
四能级系 统速率方程 其中忽略了 n3W30,因
为n3太小
dn3
dt
n0W 03
n3 S 32
A30
dn2
dt
n2
g2 g1
n1
21
,
0
v
N
l
n2
A21
S21
n3 S12
d
n
0
dt
n1 S 1 0
n0W03
n3 A30
n0 n1 n2 n3 n
dNl dt
能级的粒子数,τ1、τ2为E1、E2能级的寿命,τ21表示E2
上的粒子由于各种因素跃迁到E1造成的有限寿命;
• 这种表示方式采用激励速率和能级寿命来描述粒子数变化
速率而不涉及具体的激励及跃迁过程,而前面给出的速率
方程则忽略了激光下能级的激励过程。
连续激光器的增益和工作特性
• 增益特性是分析激光器震荡条件、模式竞争、输 出功率和激光放大器净增益系数的基础。
n1S10 n0W03 n3A30
vNl n2
A21 S21
n3S12 W03
A30
S30
n0 n1 n2 n3 n
E0
S32 E2
S21 A21 W21 W12
E1 S10
半导体激光器速率方程
第二章 光注入半导体激光器的速率方程模型2.1 光反馈半导体激光器光反馈或光注入半导体激光器的速率方程是分析和模拟系统特性的理论基础,本节先推导光反馈半导体激光器的电场速率方程―Lang-Kobayashi 方程[29],并分析了振荡条件。
为方便分析,将半导体激光器的参量及各参量的关系分别列入表2-1和表2-2。
表2-1 激光器参量的意义符号 物理量 单位 电量 C 有源区体积 m 3 载流子寿命 ns 光子寿命 ps 限制因子 --- 阈值载流子密度 m -3 透明载流子密度 m -3 增益饱和系数 m 3 线宽增强因子 --- 微分增益 m 3s -1 自发辐射因子 --- 端面强度反射率 ---波长nm表2-2 参量之间的关系Table 2-2 Relationships of parameters2.1.1 图2-1 光反馈Fabry-Perot 谐振腔示意图图2-1为光反馈的示意图,激光谐振腔两端面的反射率分别为1R 、2R ,腔长为L ,外部反射镜的反射率为e R 、距离为/2e L c τ=,τ为激光在外腔中环行一次的时间。
E +、E-分别表示正向、负向传播的时变电场的复振幅。
激光的动态变化行为取决于增益,因此可以将增益作为算子。
激光在腔内环行一次的增益为int 2())r G i kL Γg L α=-+- (2-1)将其变为指数形式,上式可变为int exp(2())r m G i kL Γg L αα=-+-- (2-2)其中/k n c ω=为波数。
实际上,激光器有源区内载流子密度()N t 随时间的变化将导致介质折射率和振荡频率的变化。
因此将波数在无光反馈阈值点(th n ,th ω)展开()()g th th th th th n n n nN N c c c N cωωωωω∂≈+-+-∂ (2-3) 其中,g th nn n ωω∂=+∂为介质的群折射率。
将(2-3)式代入(2-2)中,并将r G 分解成1r G G G ω=,其中:频率无关项1int exp[()]exp((2/)())m th th nG Γg L i L c N N Nααω∂=----∂ (2-4) 频率相关项22exp[()]g th th th n Ln L G ii c cωωωω=--- (2-5) 由于2th th n L c ω是2π的整数倍,并且角频率为ω的单色波电场满足关系式di dtω=,G ω可改写为算子exp()exp()th L LdG i dtωωττ=- (2-6) 由于激光器振荡频率在阈值附近,即th ωω≈,因此对时变复电场()e t 可引入慢变化复电场振幅()|()|exp(())E t E t i Φt =,即()()exp()th e t E t i t ω= (2-7)其中th d dtΦωω=-。
激光器速率方程讲解
n:腔内单位体积中频率处于附近单位频率间隔内 的光波模式数
得到:一个模式内的一个光子引起的受激跃迁 几率等于分配到同一模式上的自发跃迁几率。
f2 W21 al nl ,W12 al nl f1
• 对W21作出近似计算 • 设谱线的总自发辐射跃迁几率为A21,谱线 宽度为,并假设A21均匀分配在所包含 的所有模式上,则分配在一个模式上的自发 辐射跃迁几率为 A21
对表达式进行修正
dn21 ( ) sp n2 A21 ( )d n2 A21 dt dn21 ~( , ) d ( ) st n2W21 ( )d n2 B21 g 0 dt
该积分与辐射场的带宽有关。 1 原子和连续光辐射场的相互作用, 2 原子和准单色光辐射场相互作用,
3 原子和准单色光相互作用
• 由于激光的高度单色性,认为原子和准单色光相互 作用,辐射场的中心频率为 ,带宽为,且 << 。被积函数只在中心频率附近的一个极 ~( , )可 窄范围内才有非零值。在此频率范围内, g 0 以近似看成不变。 表示频率为的准 • 引入函数=(-)
f 2 A21 ~ 12 ( , 0 ) g ( , 0 ) 2 f1 8 0
2
匀加宽物质和非均匀加宽物 质的发射截面分别为
2 A21 ln 2 2 A21 21 2 2 , 21 3 2 2 4 0 H 4 0 D
~( , ) ~( , ) ~( , ) A21 g A21 g A21 g 8 2 0 0 0 W21 Nl N lV nl n 3 n n V n V c ~( , ) A21 g W21 0 al n 为腔内第 l 模内的总光子数 l nl n V
典型激光器单模振荡速率方程
• 量子理论:对光频电磁波和物质原子都作 量子化处理,并将二者作为一个统一的物 理体系加以描述(量子电动力学)。只是 在需要严格地确定激光的相干性和噪声以 及线宽极限等特性时才是必要的 • 速率方程理论:量子理论的简化形式,从 光子(量子化的电磁场)与物质原子的相 互作用出发,忽略了光子的相位特性和光 子数的起伏特性,只能给出激光的强度特 性
中心频率处的发射截面和吸收截面最大。
三、典型激光器单模振荡速率方程
• 三能级系统 • 四能级系统 根据跃迁过程写出速率方程组
思路小结: • 爱因斯坦采用唯象法得到光和物质相互作用 的关系式 • 考虑线型函数后必要的修正:几率按频率的 分布函数 • 原子和准单色光相互作用 • 单模振荡速率方程组(三能级系统和四能级 系统) • 多模振荡速率方程组
书中图4.4.3
小结(非均匀加宽): 1、频率为1 的准单色光的增益系数 、光强为 I1
gi ( 1 , I1 )
gi0 ( 1 ) 1 I1 Is
非均匀加宽工作物质的增益饱和效应的强弱 与频率无关 2、烧孔效应 反转集居数烧孔效应(书中图4.5.1)、强光入 射时弱光的增益系数(图4.5.2) 3、多普勒非均匀加宽驻波腔激光器中,强光在 弱光的增益曲线上对称地烧2个孔(图4.5.3)
1 1 1 H ( ) (下能级为基态) 2 2 2 s nr 1 1 1 H ( ) 2 2 1 1
(下能级不为基态)
在固体工作物质中占主导地位的均匀加宽是晶 格振动引起的加宽,它随温度的升高而增加。
•3、非均匀加宽
特点:原子(分子、离子)体系中每个原子只对谱线 内与它的表观中心频率相应的部分有贡献,因而可以 区分谱线上的某一频率范围是由哪一部分原子发射的 1)气体工作物质的多普勒加宽 由于气体原子的热运动,原子在光传输方向上具有 热运动速度z ,原子在自发辐射和受激辐射跃迁时表 现出来的中心频率不再是0,而是0=0(1+z/c)。 0 称作表观中心频率。 由于气体原子的热运动速度服从麦克斯韦分布,导 致了谱线的非均匀多普勒加宽。其线型函数具有高斯 线型
半导体激光器速率方程
第二章 光注入半导体激光器的速率方程模型2.1 光反馈半导体激光器光反馈或光注入半导体激光器的速率方程是分析和模拟系统特性的理论基础,本节先推导光反馈半导体激光器的电场速率方程―Lang-Kobayashi 方程[29],并分析了振荡条件。
为方便分析,将半导体激光器的参量及各参量的关系分别列入表2-1和表2-2。
表2-1 激光器参量的意义符号 物理量 单位 电量 C 有源区体积 m 3 载流子寿命 ns 光子寿命 ps 限制因子 --- 阈值载流子密度 m -3 透明载流子密度 m -3 增益饱和系数 m 3 线宽增强因子 --- 微分增益 m 3s -1 自发辐射因子 --- 端面强度反射率 ---波长nm表2-2 参量之间的关系Table 2-2 Relationships of parameters2.1.1 图2-1 光反馈Fabry-Perot 谐振腔示意图图2-1为光反馈的示意图,激光谐振腔两端面的反射率分别为1R 、2R ,腔长为L ,外部反射镜的反射率为e R 、距离为/2e L c τ=,τ为激光在外腔中环行一次的时间。
E +、E-分别表示正向、负向传播的时变电场的复振幅。
激光的动态变化行为取决于增益,因此可以将增益作为算子。
激光在腔内环行一次的增益为int 2())r G i kL Γg L α=-+- (2-1)将其变为指数形式,上式可变为int exp(2())r m G i kL Γg L αα=-+-- (2-2)其中/k n c ω=为波数。
实际上,激光器有源区内载流子密度()N t 随时间的变化将导致介质折射率和振荡频率的变化。
因此将波数在无光反馈阈值点(th n ,th ω)展开()()g th th th th th n n n nN N c c c N cωωωωω∂≈+-+-∂ (2-3) 其中,g th nn n ωω∂=+∂为介质的群折射率。
将(2-3)式代入(2-2)中,并将r G 分解成1r G G G ω=,其中:频率无关项1int exp[()]exp((2/)())m th th nG Γg L i L c N N Nααω∂=----∂ (2-4) 频率相关项22exp[()]g th th th n Ln L G ii c cωωωω=--- (2-5) 由于2th th n L c ω是2π的整数倍,并且角频率为ω的单色波电场满足关系式di dtω=,G ω可改写为算子exp()exp()th L LdG i dtωωττ=- (2-6) 由于激光器振荡频率在阈值附近,即th ωω≈,因此对时变复电场()e t 可引入慢变化复电场振幅()|()|exp(())E t E t i Φt =,即()()exp()th e t E t i t ω= (2-7)其中th d dtΦωω=-。
光纤激光器速率方程 matlab程序
光纤激光器速率方程 matlab程序光纤激光器速率方程是描述光纤激光器中光强变化随时间演化的数学模型。
光纤激光器是一种利用光纤作为传输介质的激光器,具有高功率、高效率、窄线宽等优点,在通信、材料加工、医疗等领域有着广泛的应用。
了解光纤激光器速率方程对于优化激光器设计和提高激光器性能具有重要意义。
光纤激光器速率方程描述了光纤激光器中光强的时域演化过程。
光纤激光器中的光强受到多种因素的影响,包括增益、损耗、自发辐射、受激辐射等。
光纤激光器速率方程可以通过对这些因素的定量描述,得到光强随时间的变化规律。
光纤激光器速率方程的基本形式为:dE/dt = (G - α)E - βEL其中,E为激光场强度,t为时间,G为增益系数,α为损耗系数,β为非线性系数,L为光纤长度。
这个方程可以看作是对光强的变化率进行描述,右侧第一项表示增益和损耗对光强的影响,第二项表示非线性效应对光强的影响。
在实际应用中,为了更准确地描述光纤激光器中光强的演化,可以考虑其他因素的影响,如色散、自相位调制等。
这些因素的引入可以使光纤激光器速率方程更加准确地描述光强的变化规律。
为了求解光纤激光器速率方程,可以使用数值方法进行仿真计算。
其中,最常用的方法是有限差分法。
有限差分法将时间和空间离散化,将连续的方程转化为离散的差分方程,然后通过迭代求解差分方程,得到光强随时间的变化。
在MATLAB中,可以通过编写相应的程序来求解光纤激光器速率方程。
首先,需要定义方程中的各个参数和初始条件。
然后,通过差分方法离散化方程,并利用迭代算法求解离散化后的方程。
最后,根据得到的结果,可以绘制光强随时间的变化曲线,以及其他感兴趣的物理量。
光纤激光器速率方程的求解对于优化光纤激光器的设计和改进激光器性能具有重要意义。
通过求解方程,可以获得光纤激光器中光强随时间的变化规律,进而分析和优化激光器的工作状态。
此外,光纤激光器速率方程的求解也为理论研究提供了重要的工具,可以用于研究激光器的非线性效应、自脉冲形成等现象。
- 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