利用数值方法模拟铒镱共掺光纤激光器的特性.doc

利用数值方法模拟铒镱共掺光纤激 光器的特性 课程名称： 科学计算选讲 学 院： 精仪学院 专 业：光电子与光子学技术 姓 名： 张海婷 学 号： 1013202030 一、 铒镱共掺光纤激光器振荡特性的理论基础 1、铒镱共掺系统辐射机理 铒镱共掺光纤掺杂的两种稀土元素 Er和Yb均为镧系金属，其中铒原子序数为 2 2 6 2 6 2 10 6 2 10 6 2 12 68，电子壳层： 1s 2s 2p 3s 3p 4s 3d 4p 5s 4d 5p 6s 4f 。元素镱原子序数为 70， 2 2 6 2 6 2 10 6 2 10 6 2 14 电子壳层排列为： 1s 2s 2p 3s 3p 4s 3d 4p 5s 4d 5p 6s 4f 。其中镱离子不发生受 激辐射而是主要起敏化作用， 铒离子才是产生激光的元素。 铒镱系统的能级如下 图所示： 4 2 I9/2 F Cc r 5/2 4 I11/2 4 I13/2 W56 W65 W21 W12 4 I15/2 2 F7/2 Yb3+ Er3+ 图 1 镱铒准三能级系统受激辐射 图中 Ccr为交叉弛豫系数， W12，W21， W56，W65分别为激光受激吸收、受激 辐射几率，泵浦光受激吸收、受激辐射几率。泵浦光的入射使 Yb的2F7/2 能级上 的电子吸收能量发生受激吸收跃迁到 Yb的2F5/2 能级，通过 Yb3+与 Er3+离子对的 能量传递过程把这些能量传递到 Er3+的 4I15/2 能级，该能级的电子吸收能量跃迁到激发态 4I11/2 ，4I11/2 很不稳定导致多声子弛豫使能量被传递到 4I1,3/2 ， 4I1,3/2 寿命相对较长， 随着泵浦的继续吸 , 收将形成粒子数反转， 在这个准三能级系统中同时还存在着合作上转换过程， Er3+的 4I15/2 到 4I11/2 的受激吸收、由 4I15/2 跃迁到 4I9/2 交叉弛豫过程、铒离子能级 4I9/2 和能级 4I13/2 的多声子弛豫过程等，有的对激光产生没有起到主要影响，有的影响微小甚至可以忽略。 2、铒镱共掺光纤激光器的速率方程组 将图中能级 4I15/2 、 4I13/2 、4I11/2 、4I9/2 、4F7/2 、4F5/2 分别编号为能 级 1、2、3、4、5、6，Ni （i=1 ，2， ,6 ）分别表示各能级上的单位体积的粒子 数（即粒子数密度）， W12为从能级 1 到能级 2 的跃迁几率， W21为从能级 2 到能级 1 的跃迁几率， W56为从能级 5 到能级 6 的跃迁几率， W65 为从能级 6 到能级 5 的跃迁几率。 当掺杂浓度一定时，有 N5+N6=NYb (1-1) N1+N2+N3+N4=NEr (1-2) 同时，对铒离子和镱离子上能级粒子数求偏导数， 存在下面的关系 (A32，A43 为非辐射驰豫速率， C14、Cup、 Ccr 是上转换和交叉弛豫系数 ) ： N1 W12 N1 W31 N1 W21 N 2 N 2 Ccr N1 N 6 Cup N 22 Cup N32 C14 N1 N 4 (1-3) t er N 2 W12 N1 W21 N 2 N 2 A32 N 2 2Cup N 22 2C14 N1 N4 (1-4) t er N 3 W13 N1 A32 N3 A43 N 4 C cr N1 N 6 2C up N 32 (1-5) t N 4 C up N 22 C up N 32 C14 N1 N 4 A43 N 4 (1-6) t N 6 W56 N 5 N 6 W 65N 6 Ccr N1 N 6 (1-7) t yb 当粒子数处于稳态时，上面的偏导数结果均为 0，经过整理后得 N2 W12 N1 W21 N2 N2 Ccr N1N6 Cup N22 0 (1-8) t er N 6 W56 N 5 N6 W 65N6 Ccr N1 N 6 0 t yb (1-9) 由于 4I11/2 和4I15/2 能级粒子数 N3和N4可以被忽略，因此上述 (1-2) 式也可 以简写成： N+N=N （ 1-10 ） 12er 上面（ 1-1 ）、（ 1-8 ）、（ 1-9 ）、（ 1-10 ）四式可联立成方程组，其中能 级间的跃迁几率其中能级间的跃迁几率 Wij 表达式如下（σ12、σ21 是铒离子的吸 收和发射截面面积）： W12 ( z, t) s 12 P( z) s （1-11 ） h s A W21 ( z, t) W56 ( z, t)  21 Ps (z) h s A 56 Pp ( z) h p A  (1-12) (1-13) W65 ( z, t)  p 65 h p A  Pp (