《NaI(Tl)闪烁晶体原理》由会员分享,可在线阅读,更多相关《NaI(Tl)闪烁晶体原理(9页珍藏版)》请在金锄头文库上搜索。
1、.附录一 NaI(Tl)闪烁晶体闪烁体按其化学性质可分为两类:一类是无机晶体闪烁体,通常是含有少量杂质(称为激活剂)的无机盐晶体,如碘化钠(铊激活)单晶体、即NaI(Tl),碘化铯(铊激活)单晶体、即CsI(Tl),硫化锌(银激活)、即ZnS(Ag)等;另一类是有机闪烁体,它们都是苯环碳氢化合物。闪烁体的发光机制比较复杂,在此对无机晶体闪烁体的发光机制作一些简要的定性介绍。无机晶体闪烁体属离子型晶体,原子(离子)之间结合得比较紧密相互之间影响比较大,晶格中原子电子能级加宽成为一系列连续的能带。其中最低能量状态已为电子所填满,故称为满带;价电子都处于稍高的能量状态,这种能带称为“价带。假设价带未
2、填满,那么在外电场作用下将有净电流产生;假设价带已填满,那么必须有电子被激发到更高的能带导带上去,才能产生电流,此时价带上有一空穴,导带上有一电子,即产生了一个自由电子空穴对。价带与导带之间的空隙中不存在电子能级,称为禁带;禁带有一宽度Eg,它和晶体的导电性质密切相关,导体在0.1eV左右,半导体在0.632.5eV之间,无机闪烁体为绝缘透明物质,Eg3eV,NaI为7.0eV。也存在另一种情况:在闪烁晶体中产生的电子空穴对仍束缚着,称为“激子,它们在晶格中一起运动,在外电场中无净电流产生,其能带在导带之下,称为“激带。自由的导带电子和价带空穴可以复合成激子,激子也可以吸收热运动能量变成自由电
3、子空穴对。当核辐射进入闪烁体时,既可产生自由电子空穴对,也可以产生激子。而后电子从导带或激带跃迁到价带,退激过程中放出光子;也存在着竞争过程非辐射跃迁,即通过放热(晶格振动)退激。有一点需要指出,纯的NaI晶体不是有效的闪烁体。一是因为相应禁带宽度的光子能量在紫外光围,不是可见光;二是退激发出的光子尚未逸出晶体就会被晶体自身吸收。为了解决这一问题,在纯晶体中掺入少量杂质原子(如Tl),称为“激活剂,它们成为发光中心,形成一套激发能级,能量比导带低,而基态却比价带高,这样跃迁产生的光子能量就比禁带宽度Eg小,那么它就不可能再使价带上的电子激发到导带上去,从而避免自吸收。碘化钠闪烁晶体能吸收外来射
4、线能量使原子、分子电离和激发,退激时发射出荧光光子。NaI(Tl)晶体的密度较大(=3.67g/cm3),而且高原子序数的碘占重量的85%,所以对射线的探测效率特别高,同时相对发光效率大;它的发射光谱最强波长为415nm左右,能与光电倍增管的光谱响应较好匹配。此外,晶体的透明性也很好,测量射线时能量分辨率也是闪烁体中较好的一种。一个需要指出的问题是:在闪烁体的选取上要注意闪烁体对所测的粒子要有较大的阻止本领,以使入射粒子(特别是能量较大的粒子)在闪烁体中能损耗较多的能量而退激产生光子。原先使用的国产NaI(Tl)晶体尺寸为f205mm,这一厚度对定标时测高能(E1MeV)时的效率不够高,而且对
5、高能粒子的计数率也比较低;本装置的闪烁探测器采用的尺寸为f2020mm的NaI(Tl)晶体可以说是一大改进,一方面可以提高探测高能部分的效率,另一方面也提高了实验中高能粒子的计数率。NaI(Tl)晶体的缺点是容易潮解,吸收空气中的水分就会变质失效;因此我们采用了200mm的铝来密封;这就需要对穿过Al膜时的能量损失进行修正。在实验中我们发现,对于不同的粒子能量的损失不尽相同;所以在实际的实验和数据处理中进行了能量损失的合理修正。附录二 光电倍增管光电倍增管是一种常用的灵敏度很高的光探测器,它由光阴极、电子光学输入系统、倍增系统及阳极组成,并且通过高压电源及一组串联的电阻分压器在阴极打拿极(又称
6、“倍增极)阳极之间建立一个电位分布。光辐射照射到阴极时,由于光电效应,阴极发射电子,把微弱的光输入转换成光电子;这些光电子受到各电极间电场的加速和聚焦,光电子在电子光学输入系统的电场作用下到达第一倍增极,产生二次电子,由于二次发射系数大于1,电子数得到倍增。以后,电子再经倍增系统逐级倍增,阳极收集倍增后的电子流并输出光电流信号,在负载电阻上以电压信号的形式输出。K光阴极;F聚焦极;D1D10打拿极;A阳极。根据打拿极的几何形状和排列方式,光电倍增管分为聚焦型(环状、直线)和非聚焦型(百叶窗式、盒栅式)。本装置采用GDB44F型百叶窗式光电倍增管,其优点为脉冲幅度分辨率较好,适用闪烁能谱测量。它
7、的主要指标应该包括以下几方面:光电转换特性、电子倍增特性、噪声或暗电流、时间特性等;在此主要介绍光电转换特性和电子倍增特性。1) 光电转换特性光阴极的光谱响应和灵敏度光阴极是接收光子并放出光电子的电极,一般是在真空中把阴极材料蒸发在光学窗的表面上,形成半透明的端窗阴极;光阴极材料的品种有数十种,但最常用的只是五、六种,如锑铯化合物等。一般光电倍增管光阴极前的光学窗有两种:硼玻璃窗或石英窗,前者适用于可见光,后者可透过紫外光。光阴极受到光照射后发射光电子的几率是波长的函数,称为光谱响应。在长波端的响应极限主要由光阴极材料的性质决定,而短波端的响应主要受入射窗材料对光的吸收所限制。了解光电倍增管的
8、光谱响应特性有利于正确选择不同管子使之与闪烁体的发射光谱相匹配。在实际应用中,光电转换特性通常使用另一个宏观定义,即一定通量F的白光照射阴极所能获得的光电子流ik称为光阴极光照灵敏度: 1其中ik单位为微安;F为光通量,单位为“流明(lm)。2) 电子倍增特性光电倍增管的放大倍数及阳极灵敏度 光电倍增管的放大倍数增益M由于打拿极的倍增作用,从光阴极发射出来的电子不断被倍增,最后可在阳极上得到大量电子。从光阴极射出,到达第一打拿极的一个电子,经过多次倍增后在阳极得到的电子数,称为光电倍增管电流放大倍数增益。在理想情况下一般可写成: 2式中d是平均二次发射系数,n为打拿极的级数。二次发射系数d是极
9、间电压的函数,可用经验公式表示: 3其中VD为打拿极之间的电压,a、b为经验常数。如果打拿极电子传递效率为g,那么增益M比较实际的表达式可写成: 4对设计良好的聚焦型管子g约等于1,对非聚焦型管子g1。 阳极光照灵敏度S放大倍数是光电倍增管的重要参数之一,但往往有些技术说明书不直接给出它的数值,而是在给出光阴极光照灵敏度Sk的同时,给出光电倍增管的“阳极光照灵敏度Sa,它们之间的关系是: 5其中Sa的单位为A/lm,gc为第一打拿极对光电子的收集效率。阳极光照灵敏度的物理意义是:当一个流明的光通量照在光阴极上时,在光电倍增管阳极上输出的电流阳极电流ia的数值。当入射光通量F增大时,阳极电流ia
10、在相当宽的围是线性增大的;但F太大时,就出现偏离线性。原因之一是打拿极发射二次电子疲劳,使放大倍数减小;其二是最后几级打拿极和阳极上有空间电荷堆积;也有可能是分压电阻选择不当,使最后几级打拿极以及阳极之间的电压降低,放大系数减小,这一问题可以通过调整分压电阻来解决。阳极光照灵敏度Sa和总电压的关系由式(3)、(4)、(5)可知:,故,两个量的对数成线性关系;因而随着电流增加到某一数值会出现非线性,logSa增加变得缓慢;一般说来,加在光电倍增管上的高压在1000V之线性还是比较理想的。需要指出的是:闪烁探测器的线性问题是由多个因素共同作用的结果,不仅光电倍增管是个重要因素,闪烁晶体本身也存在能
11、量线性问题。因此在实际应用中,必须考虑多方面的因素,比如各部件的匹配等,而常用的解决方法那么是调整光电倍增管的工作参数。光电倍增管的管脚插入底座。底座是由分压器与射极跟随器组成。1分压器光电倍增管中各电极的电位由外加电阻分压器抽头供给。本实验使用正高压电路,阴极接地,阳极处于高电位,输出端使用耐高压电容隔开。所加电压应根据说明书或不同用途以及管子的性能进行考虑;建议用户在使用本实验装置时采用我们的推荐值。2射极跟随器射极跟随器具有电流放大作用(放大倍数一般为几十一百以上),但其电压放大倍数恒小于1而接近于1,且输出电压和输入电压同相,因此具有电压跟随的特点,频率响应较好。附录三 多项式拟合及能
12、谱的光滑处理和寻峰一 多项式拟合1概述在实际的物理实验中,我们经常测出某一未知形式函数在假设干点处的值i=1,2n,所要做的是根据这些数值点推导的近似表达式,这就是所谓的曲线拟合问题。在介绍曲线拟合之前首先要谈一下函数插值。所谓函数插值是根据插值原理建立一个次数不高于n的插值多项式作为函数的近似表达式,基本要拟合曲线通过所有数据点。而曲线拟合不要求通过所有数据点,只要求得出的近似函数能反映数据的基本关系,因此曲线拟合比函数插值得到的结果更能反映客观实际。在某种意义上来说曲线拟合还更有实用价值,因为实际问题中所提供的数据点往往很多,用函数插值势必得出次数很高的插值多项式,导致计算上的很多麻烦,而
13、且高次多项式的振荡问题也是一大弊端。因此实际工作中更多采用了曲线拟合的处理方法。那么对给出的数据点作拟合曲线时,怎样才算“拟合得最好?或者说采取的原那么是什么?一般希望使各观测数据与拟合曲线的方差和最小,就能使拟合曲线更接近真实函数,即最小二乘原理。2算法用一个m次多项式来拟合n个观测数据点,其中m远小于n。误差平方和:,要求F的极小值那么必须分别对a0、a1、am求偏导数并使之为零,因此得到:; 31即:,即系数a0、a1、am应满足如下方程组: 32此方程组中有m+1个方程,通常称为法方程;其系数矩阵为一个对称矩阵,并且是正定的,可以求出唯一解a0、a1、am而后代入m次多项式即可得拟合结
14、果。3多项式拟合阶数的选取在对实验数据点进行多项式拟合的过程中,拟合阶数的选取是一个至关重要的问题;如果拟合阶数选取过低那么无法找准y与x的关系,假设选取过高那么高次多项式的振荡问题在所难免。在讨论这一问题之前先介绍几个概念。假设用来拟合数据点,可以根据上面介绍的方法算出系数、代入拟合多项式,那么拟合值,单个的残差为: 33残差平方和为: 34可以用来估计观测值方差: 35因此拟合阶数K的选取是否合适将直接关系到观测数据的方差估计。K选取过小那么由于未找准y与x的关系而引起很大误差;K选取过大那么由于R趋于稳定,而式35根号中分母变小同样会使方差估计偏大。原那么上选取对应于R随K值开始变化不大
15、时的K值。什么叫“R随K值开始变化不大呢?这需要一个定量的判据,下面介绍一个统计性判据。按K阶拟合可以得到相应的残差平方和RK,按K1阶拟合可以得到RK1。根据R的统计性质可以证明,参见复旦、清华、北大合编“原子核物理实验方法第十八章的附录四。假如由这两种阶数多项式拟合得到的方差在某一显著水平上没有明显差异,那么不必用K阶而考虑用K1阶来拟合;反之假设方差有较大差异,那么考虑用K阶来拟合。要检验两种分布的方差是否一致,可以用F分布来检验。假设,令,由分布加减法可以得到: 36而和互相独立,因此由F分布定义:,即 37也就是说FK服从第一自由度为1,第二自由度为nk1的F分布有关F分布可以参考概率论和数理统计的相关章节。具体方法为:用最小二乘法对观测数据分别作K阶和K1阶的多项式拟合,算出相应的RK,RK1和SK;按显著水平或查F分布表,找出相对于自由度为1和nk1的F分布的临界值,比较和。当说明由K、K1阶两种拟合结果一致,因此不必用K;当,说明两者有显著差异,有必