J-O理论的计算讲解.doc

J-O理论的计算 在应用J-O理论模型来计算各种光学参数的过程中,参数的单位均采用高斯制(CGS)的表达方式,又称为厘米 ·克 ·秒制.以下以掺Er3+的样品为例,详细说明了各参数的计算方法,此方法可以应用到其他稀土掺杂玻璃的材料上,具有很强的通用性. 根据实验所测得的吸收光谱上的吸收峰列出其所对应的吸收跃迁过程如表1 .吸收过程 计算实验振子强度fexp以及理论振子强度fmd 根据表1可以看出,吸收光谱中共有11个吸收峰,对于每个吸收峰而言,都需要计算振子强度fexp和fmd.计算fexp的关键在于计算各个吸收峰的平均波长以及光密度的积分.在计算平均波长时应考虑各波长点上吸收强度的贡献.而对于离散数值计算来说,积分的计算就是求和的运算.由此根据式(1)不难算出实验振子强度.根据式(2)和表1可以看出在计算fmd时,公式中的J为吸收跃迁基态能级的总角动量量子数,无论计算哪个吸收峰的fmd, J的取值都是15/2.对标量ΔJ进行检测,根据不同的检测结果,计算出不同的矩阵元〈γJ‖L+ 2 S‖γ′J′〉.如此便可以无一遗漏地计算出存在磁偶极跃迁的理论振子强度fmd. 计算强度参数Ωt( t = 2 , 4 , 6)以及理论振子强度fed 由于强度参数Ωt在整个J-O理论模型中处于核心的地位,因此这三个参数的计算也变得十分重要.值得说明的是,本文将多次以表达式f(i)的形式来表示第i个跃迁过程中物理量f的值,目的是为了与C或MATLAB等程序语言中的数组概念联系起来，方便理解和阐述.可以通过建立数组U与各吸收峰所对应的张量算符U( t)的约化矩阵元值,这里采用文献[ 7 ]中的数据.根据最小二乘法的定义可以知道,拟合得出的Ω2, Ω4和Ω6需使目标函数F= 取到最小值,其中的Sed(i)可以根据公式fed( i) = fcal( i) - fmd( i)计算得到,这里假设fexp(i) = fcal(i) .因此可以通过多元函数求极值的方法,联立方程组解得参数Ω2,Ω4和Ω6.随后重新运这三个参数计算Sed( i)和Fed( i) , i = 1 ,2 ,…,11 .从而方便地利用公式(8)算出拟合误差. 计算自发辐射几率Arad[γJ ,γ′J′] 从式(9)可以看出,计算Arad的关键在于计算自发辐射跃迁中电偶极和磁偶极的跃迁振子强度Sed和Smd.由拟合得到的三个强度参数Ω2,Ω4和Ω6,结合文献[ 11 ]中的Er3 +在各能级间自发辐射跃迁的约化矩阵元〈γJ‖U( t)‖γ′J′〉2的值,根据式(3)便可以计算出自发辐射跃迁振子强度Sed.在表2中列出了Er3 +可能发生自发辐射跃迁的能级对. 从表中可以看出共存在21对自发辐射跃迁能级对,每个辐射跃迁过程均对应了一个Sed, Smd和Arad.在计算磁偶极辐射跃迁振子Smd时,不同于吸收跃迁过程,其初始能级的J是不同的.具体数值详见表2 . 计算荧光分支比β和辐射寿命时间τrad 由表2可知共存在21个辐射跃迁过程.对应的每个跃迁过程都需要计算荧光分支比β.值得注意的是,跃迁过程的初态能级只有9个.由β的定义可知,要计算第u个辐射跃迁过程的荧光分支比β( u)时,分子为相应的自发辐射跃迁几率Arad( u) ,而分母为所有发自于此初态能级i的自发辐射跃迁几率之和Asum( i) =∑uArad( u) ,其中i , u与Asum( i)的关系见表3 　　例如,在计算β(2)和β(3)时,分子分别为辐射跃迁几率Arad(2)和Arad(3) ,对照表3的第3行和第4行的数据可知,分母均为Asum(2) = ∑3u = 2Arad( u) .而由荧光寿命的定义可知,只需将表3中的Asum( i)分别取倒数便可以计算出所需要的9个自发辐射初态的荧光寿命时间τrad.结合计算出来的参数以及测量得到的吸收和荧光光谱,根据式(10)和(11) ,也不难算出吸收和自发辐射跃迁截面.