留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

离子能量分析器测量特性的仿真研究

翟红雨 程健 陈银华 陆伟

王凯, 吕学阳, 吴锟霖, 等. 不同顺序中子/γ辐照对双极器件电流增益的影响[J]. 强激光与粒子束, 2020, 32: 044001. doi: 10.11884/HPLPB202032.190333
引用本文: 翟红雨, 程健, 陈银华, 等. 离子能量分析器测量特性的仿真研究[J]. 强激光与粒子束, 2020, 32: 084002. doi: 10.11884/HPLPB202032.190459
Wang Kai, Lü Xueyang, Wu Kunlin, et al. Effects of different sequential neutron/gamma irradiation on current gain of bipolar devices[J]. High Power Laser and Particle Beams, 2020, 32: 044001. doi: 10.11884/HPLPB202032.190333
Citation: Zhai Hongyu, Cheng Jian, Chen Yinhua, et al. Simulation study on measurement characteristics of ion energy analyzer[J]. High Power Laser and Particle Beams, 2020, 32: 084002. doi: 10.11884/HPLPB202032.190459

离子能量分析器测量特性的仿真研究

doi: 10.11884/HPLPB202032.190459
基金项目: 国家自然科学基金项目(11575182)
详细信息
    作者简介:

    翟红雨(1996—),男,硕士研究生,从事嵌入式测量系统的理论与仿真分析;hyzhai@mail.ustc.edu.cn

    通讯作者:

    程 健(1964—),男,副教授,主要从事等离子体测量设备系统的研制与应用;chengj@ustc.edu.cn

  • 中图分类号: O53

Simulation study on measurement characteristics of ion energy analyzer

  • 摘要: 针对空间等离子体及其模拟环境、空间原子氧及其模拟环境对离子能谱测量的需要,利用仿真软件COMSOL,对离子能量分析器的低能离子测量特性进行了仿真研究。介绍了离子能量分析器的工作原理,对离子能谱测量过程进行了公式推导。通过对三种待选仪器设计方案进行离子透过率仿真分析,确定了一种较优的仪器设计方案。多种离子温度下的误差分析结果也表明,该设计方案能够较为准确地测量离子能量分布。分析了电场畸变、等离子鞘层、栅网对齐方式和离子温度对测量结果的影响,根据仿真结果对一些仿真实验现象做出了合理的解释。
  • 在中子和γ射线同时存在的混合辐射下的双极器件同时表现出位移损伤和电离损伤。一般假设这两种效应在辐射损伤评估中是相互独立的,它们对器件的综合影响是单个中子位移效应和总剂量电离效应引起的退化的简单总和。然而,近年来的研究[1-6]指出,这两种效应相互之间具有非线性影响,并且它们的组合效果不等于它们的简单相加。

    Gorelick等人[2] 研究了12种先暴露于中子然后暴露于γ射线的线性集成电路的辐射损伤。实验结果表明,半导体器件在中子照射后表现出对电离效应有一定程度的增强或抑制,总剂量电离效应和位移损伤无法解耦。Barnaby等人[5-6]研究了在质子与γ射线以及质子与中子同时存在辐射环境下双极运算放大器的辐射损伤效应,他们的实验数据显示,这些运算放大器及其输入晶体管对位移效应和电离效应引入的缺陷的综合影响的电响应是非线性的。李兴冀等人[7-8]对由质子和电子辐照导致的PNP和NPN双极晶体管的协同辐照效应进行了实验研究,其实验数据表明,电子导致的电离损伤效应可能会对质子诱导产生的位移损伤产生退火作用或增强效果。王晨辉等人[9]利用半导体器件仿真软件TCAD,通过改变少数载流子寿命、在氧化层中添加陷阱电荷以及增大表面复合速率等措施,对中子和γ射线混合辐射引起的6种横向PNP双极晶体管的电离/位移协同效应进行了数值模拟,结果表明,横向PNP双极晶体管的电离/位移协同效应不是总剂量电离效应和位移效应的简单总和,并且总剂量电离效应可以增强中子位移损伤,导致更大的增益退化。他们认为,氧化物层中的正电荷和γ辐照诱导的Si/SiO2界面陷阱可以增强中子辐照诱导的体缺陷中载流子的复合过程,是导致横向PNP双极晶体管产生电离/位移协同效应的主要原因。Song Yu等人[10]发现,在样品被中子预辐照后,辐射损伤表现出对γ射线剂量的“√”样关系,即损伤对于小的伽马射线剂量异常降低,然后对于大的伽马射线剂量几乎线性地增加。对实验数据分析的总损伤小于中子和γ损伤简单相加之和,即,表现为负协同效应。基于实验数据推导出两种负协同效应与γ射线剂量之间的关系:第一种负协同效应与γ射线剂量呈线性关系,其斜率与初始位移损伤呈二次关系,可归因于硅中中子辐射引起的缺陷的愈合;第二种负协同效应随γ射线剂量指数衰减,其振幅表现出较强的低剂量率增强(ELDRS)效应,可归因于硅中γ射线产生的正电荷对Si/SiO2界面附近中子诱导缺陷的钝化作用。以上研究均表明,电离/位移协同效应并不是电离损伤和位移效应的简单相加之和,电离损伤对已存在的位移损伤具有一定的增强或抑制作用。在以往的研究中,对于双极晶体管电离/位移协同效应的研究主要通过质子或重离子辐照同时产生电离效应和位移效应,无法对单独的电离效应或位移效应进行分析。王晨辉和Song Yu等人的研究中虽然使用中子和γ射线作为辐照源,但是研究对象主要为PNP型晶体管,且前者以模拟为主要研究手段,缺少对NPN型双极晶体管的实验研究,后者主要研究了先中子后γ辐照实验,未研究先γ后中子辐照实验。

    本文采用中子源和γ源对NPN及PNP晶体管进行不同顺序的辐照实验,分析了中子/γ的不同顺序辐照对不同类型双极晶体管造成的电流增益变化与电学特征参数之间的关系,包括先中子后γ辐照条件下和先γ后中子辐照条件下双极晶体管共射极电流增益随集电极电流的变化关系等。

    本文选用的实验样品为A组(NPN型)和B组(PNP型)两种类型的双极晶体管,为尽可能减小器件个体差异的影响,每种类型实验样品都选自同一批次。实验分别从两种类型晶体管中选取六只,其中A组NPN型双极晶体管编号为2#,4#,5#,6#,8#,10#,B组PNP型双极晶体管编号为3#,4#,5#,6#,8#,9#。三只作为先中子后γ辐照实验样品,另外三只作为先γ后中子辐照实验样品,即A组2#,4#,5#样品和B组3#,4#,5#样品进行先中子后γ辐照实验,A组6#,8#,10#样品和B组6#,8#,9#样品进行先γ后中子辐照实验。实验流程如图1所示。

    图  1  实验流程图
    Figure  1.  Experimental flow chart

    辐照实验依托中国工程物理研究院核物理与化学研究所CFBR-II快中子脉冲堆和钴源装置开展,设定中子辐照累积注量为1.2×1012 cm-2,γ辐照累积注量为612.6 Gy(Si)。两轮实验中中子注量率和γ剂量率相同。辐照前后所有器件效应参数均在BC3193半导体分立器件测试系统上测试,测试环境温度保持在(25±0.5)℃,环境湿度保持在55%~65%范围内。

    对A组-NPN型双极晶体管而言,当集电极−发射极电压VCE偏置一定时,先中子后γ辐照或者先γ后中子辐照导致的双极晶体管电流增益hFE的退化很接近。随集电极电流IC增大,电流增益hFE逐渐增大。如图2所示,在集电极电流IC比较小时电流增益退化最大可达到约59.9%,在所测电流范围内,hFE退化最小约22.8%。可以看出,当集电极-发射极电压VCE一定时,晶体管的增益退化程度与注入电流相关,电流越小,增益退化程度越大。

    图  2  A组-NPN型双极晶体管电流增益hFE随集电极电流IC的曲线
    Figure  2.  Curves of current gain hFE vs collector current IC of A-NPN BJTs

    对B组-PNP型双极晶体管而言,当集电极-发射极电压VCE偏置一定时,先中子后γ辐照相较先γ后中子辐照引起的双极晶体管电流增益hFE退化程度不同。如图3所示,可以明显看出,先中子后γ辐照造成的电流增益退化程度比先γ后中子辐照造成的电流增益退化程度大。先中子后γ的增益退化程度在24.4%~47.0%范围内变化,而先γ后中子的电流增益退化程度在10.7%~31.4%范围内变化,两者相差近1倍。随集电极电流IC增大,电流增益hFE的退化程度逐渐减小。先中子后γ辐照情况,在集电极电流IC比较小时电流增益退化最大约47.0%,随IC逐渐增大,退化程度逐渐减小,在所测电流范围内,电流增益退化最小约24.4%;先γ后中子辐照情况,电流增益退化最大约31.4%,最小约10.7%。

    图  3  B组-PNP型双极晶体管电流增益hFE随集电极电流IC的变化曲线
    Figure  3.  Curves of current gain hFE vs collector current IC of B-PNP BJTs

    根据实验结果可知,当集电极−发射极电压偏置一定时,不论是A组-NPN型双极晶体管还是B组-PNP型双极晶体管,辐照后的电流增益变化趋势与辐照前相同,在集电极电流极低的情况下,电流增益下降比较严重,随集电极电流增加,电流增益逐渐增大。

    图4为NPN型晶体管工作与正向有源区时电子流和空穴流示意图,在集电极电流较小情况下,发射区注入到基区的少数载流子很少,而发射极-基极的体内存在着工艺条件引起的复合中心,导致进入基区的少数载流子大部分被复合,进一步降低了基区内少数载流子数量,发射极注入到基区的电子电流减小,通过基区的少数载流子扩散电流减小,发射效率γ降低,最终导致电流增益下降严重[11]。此外,基区表面复合电流增加和发射结空间电荷区复合电流增加也是引起电流增益下降的原因[12]。随着集电极电流进一步增大,发射极注入到基区的少数载流子浓度和空穴浓度增大,基区电子输运速度增大,基极电流线性增加(如图5所示),发射效率增大,因此,电流增益hFE随集电极电流IC增大而逐步增大。

    图  4  NPN型晶体管工作于正向有源区时电子流和空穴流示意图
    Figure  4.  Diagram of electron flow and hole flow when NPN transistor is working in positive active region
    图  5  A组-NPN型和B组-PNP双极晶体管基极电流IB随集电极电流IC的曲线
    Figure  5.  Curves of base current IB vs. collector current IC of A-NPN BJTs and B-PNP BJTs
    图  6  先γ后中子入射辐照效应示意图(以PNP型双极晶体管为例)
    Figure  6.  Schematic diagram of irradiation effect of neutron first followed by γ incidence (taking PNP bipolar transistor as an example)

    根据图23可知,对于双极晶体管,先γ后中子辐照引起的电流增益退化低于先中子后γ辐照引起的增益退化,该现象在A组NPN型双极晶体管中表现的不太明显,但在B组PNP型双极晶体管中十分明显,接近二倍关系。

    基于以上实验数据,我们进行了如下定性分析。如图6所示,先γ辐照时,γ射线在硅体中产生电子-空穴对和少量的空位和间隙原子,在氧化层产生净的正电荷。产生的电子一部分会输运到Si-SiO2界面与氧化层正电荷复合,一部分会与空穴复合。产生的空位与杂质原子、空位相结合形成稳定的缺陷中心,如V-P、V-O、V-V等[13-14]。后续中子辐照产生大量的空位和间隙原子以及点缺陷甚至缺陷簇。中子辐照产生的空位会与γ辐照产生的空位结合形成双空位(V-V)中心,与杂质原子、间隙原子以及γ辐照产生的缺陷中心(V-P、V-O、V-V等)相结合,降低γ辐照产生的缺陷浓度,从而表现为对γ辐照的“退火”作用。如图7所示,先中子辐照时,中子辐照在硅体中产生的空位会形成双空位以及与杂质原子结合后的稳定的缺陷中心,γ辐照产生的空位会与这些缺陷中心相结合增大缺陷密度,同时还会与中子辐照产生的点缺陷或缺陷簇相结合,扩大缺陷的损伤能力,从而表现为后γ辐照对先中子辐照产生的损伤的“诱导”作用。

    图  7  先中子后γ入射辐照效应示意图(以PNP型双极晶体管为例)
    Figure  7.  Schematic diagram of irradiation effect of γ first followed by neutron incidence ( taking PNP bipolar transistor as an example)

    针对两种类型的晶体管损伤差异的不同,可归结为两点原因:(1)A组样品为NPN型双极晶体管,B组样品为PNP型双极晶体管,两种类型的样品结构不同;(2)A组为高频小功率NPN管,B组为小功率PNP管,两种类型的样品截止频率不同。双极晶体管的截止频率fT是指共射极电流增益为1时的频率,通常可表示为[15]

    fT=12πτec
    (1)

    式中:τec为发射极到集电极的时间,由四部分组成

    τec=τe+τb+τd+τc
    (2)

    式中:τe为发射结结电容充电时间;τb为基区渡越时间;τd为集电结耗尽区渡越时间;τC为集电结电容充电时间。其中τbτd为主要影响因素,可表示为

    τb=x2B2Dn
    (3)
    τd=xdcvs
    (4)

    式中:xB为基区宽度;Dn为电子扩散系数;xdC为基极−集电极的耗尽区宽度;vs为电子饱和速度。由此可知,要使截止频率较高,需要有较小的基区宽度和基极−集电极耗尽区宽度。因此A组晶体管的基区宽度小于B组的基区宽度,而γ辐照会造成PN结耗尽层向P侧扩散,A组晶体管的基区宽度将会进一步降低,导致产生的空位浓度降低,对不同顺序辐照损伤差异的影响减小。因此不同顺序中子/γ辐照后,A组增益退化的差异小于B组的差异。

    本文以CFBR-II快中子反应堆和钴源装置为辐照装置,以A型NPN双极晶体管和B型PNP双极晶体管为实验对象,通过研究不同顺序的中子和γ射线辐照引起的双极晶体管增益退化,达到对双极晶体管电离/位移协同效应的研究目的,对双极晶体管在强核辐射场或空间航天器应用中的抗辐射加固以及硬度评估等方面提供了参考。通过分析实验结果可以得出以下结论:

    (1)当集电极电流很小的情况下,电流增益退化比较严重,随集电极电流增大,电流增益逐渐增大。主要原因在于:集电极电流很小时,发射极注入到基区的少数载流子浓度低,基区体复合电流较大,基区表面复合电流和发射结耗尽区复合电流相对增大,引起发射效率γ降低,最终导致集电极电流较小的情况下电流增益退化严重;随集电极电流逐渐增大,基区内少数载流子浓度增加,基区电子输运速度增大,基极电流线性增加,发射效率增大,最终导致电流增益hFE随集电极电流IC增大而逐步增大。

    (2)先中子后γ辐照引起的电流增益退化程度比先γ后中子辐照引起的电流增益退化程度更深,此现象在A组型NPN双极晶体管中表现得不太明显,但是在B组型PNP双极晶体管中,先中子后γ辐照引起的电流增益退化程度与先γ后中子辐照引起的电流增益退化程度约为二倍关系。初步分析原因为:针对实验选用的两类双极晶体管,γ射线辐照引起的点缺陷对中子预辐照产生的点缺陷和缺陷簇具有“诱导加剧”作用,中子辐照引起的点缺陷和缺陷簇对γ预辐照引起的点缺陷和缺陷簇具有“退火淹没”作用。

    (3)两种类型晶体管辐照结果的差异主要由A组NPN型双极晶体管截止频率较高导致,截止频率高,要求在制造时基区宽度要小于非高频晶体管,而γ辐照会造成PN结耗尽层向P侧扩散,A组晶体管的基区宽度将会进一步降低,导致产生的空位浓度降低,对不同顺序辐照损伤差异的响应减小。

  • 图  1  离子能量分析器典型结构

    Figure  1.  Typical structure of ion energy analyzer (IEA)

    图  2  理想情况下离子能量分析器测量曲线

    Figure  2.  Ideal IEA measurement curve

    图  3  模型简化过程

    Figure  3.  Model simplification processes

    图  4  三种设计方案的离子透过率曲线

    Figure  4.  Ion transmission curve of three schemes

    图  5  方案A和方案B的离子透过率曲线

    Figure  5.  Ion transmission curve for scheme A and scheme B

    图  6  离子能量分析器内部电场图

    Figure  6.  Internal electric field of IEA

    图  7  不同入射能量下的离子运动轨迹

    Figure  7.  Ion trajectories at different initial energies

    图  8  栅网对齐的两种情况

    Figure  8.  Two cases of grid alignment

    图  9  两种对齐方式下的离子透过率曲线

    Figure  9.  Ion transmission curve in two alignments

    图  10  不同仪器电位下的离子透过率曲线

    Figure  10.  Ion transmission curve at different instrument potentials

    图  11  不同仪器电位下的I-V曲线

    Figure  11.  I-V curve at different instrument potentials

    图  12  不同温度下的I-V曲线和离子能量分布

    Figure  12.  I-V curve and ion energy distribution at different temperatures

    表  1  不同温度下离子速度拟合结果

    Table  1.   Ion velocity fitting results at different temperatures

    temperature/Kaverage velocity/(m·s−1relative error/%εRMSE/eV
    5007593−0.090.05
    100076260.340.11
    150076570.750.07
    200078112.770.73
    下载: 导出CSV
  • [1] Bourdeau R E, Whipple E C Jr, Donley J L, et al. Experimental evidence for the presence of helium ions based on Explorer VIII satellite data[J]. Journal of Geophysical Research, 1962, 67(2): 467-475. doi: 10.1029/JZ067i002p00467
    [2] Sarkar S, Gwal A K, Parrot M. Ionospheric variations observed by the DEMETER satellite in the mid-latitude region during strong earthquakes[J]. Journal of Atmospheric and Solar-Terrestrial Physics, 2007, 69(13): 1524-40. doi: 10.1016/j.jastp.2007.06.006
    [3] 陳信良. 中華衛星一號電離層電漿電動效應儀 (IPEI) 資料之頻譜分析[D]. 台灣: 國立中央大學, 2009: 4-7.

    Chen Hsinliang. Spectrum analysis of the data of China Satellite No. 1 Ionospheric Plasma and Electrodynamics Instrument (IPEI)[D]. Taiwan: National Central University, 2009: 4-7
    [4] Sgrigna V, Buzzi A, Conti L, et al. The ESPERIA satellite project for detecting seismo-associated effects in the topside ionosphere: First instrumental tests in space[J]. Earth Planets Space, 2008, 60(5): 463-475. doi: 10.1186/BF03352813
    [5] Satir M, Sik F, Turkoz E, et al. Design of the retarding potential analyzer to be used with BURFIT-80 ion thruster and validation using PIC-DSMC code[C]/7th International Conference on Recent Advances in Space Technologies. 2015.
    [6] Goldan P D, Yadlowsky E J, Whipple E C. Errors in ion and electron-temperature measurements due to grid plane potential nonuniformities in retarding potential analyzers[J]. Journal of Geophysical Research, 1973, 78(16): 2907-2916. doi: 10.1029/JA078i016p02907
    [7] Troy B E, Maier E J. Effect of grid transparency and finite collector size on determining ion temperature and density by retarding potential analyzer[J]. Journal of Geophysical Research−Space Physics, 1975, 80(16): 2236-2240. doi: 10.1029/JA080i016p02236
    [8] Knudsen W C. Finite grid radius and thickness effects on retarding potential analyzer measured suprathermal electron-density and temperature[J]. Journal of Geophysical Research−Space Physics, 1992, 97(A9): 13767-13775. doi: 10.1029/92JA00642
    [9] Chao C K, Su S Y. Charged particle motion inside the retarding potential analyzer[J]. Physics of Plasmas, 2000, 7(1): 101-107. doi: 10.1063/1.873817
    [10] Chao C K, Su S Y, Yeh H C. Grid effects on the derived ion temperature and ram velocity from the simulated results of the retarding potential analyzer data[J]. Advances in Space Research, 2003, 32(11): 2361-2366. doi: 10.1016/S0273-1177(03)90566-7
    [11] Klenzing J H, Earle G D, Heelis R A, et al. A statistical analysis of systematic errors in temperature and ram velocity estimates from satellite-borne retarding potential analyzers[J]. Physics of Plasmas, 2009, 16: 052901. doi: 10.1063/1.3125311
    [12] 冯宇波, 王世金, 孔令高. 电离层星载阻滞势分析器的误差仿真分析[J]. 地球物理学报, 2010, 53(11):2535-2543. (FengYubo, Wang Shijin, Kong Linggao. Simulation analysis of errors from ionosphere satellite-borne retarding potential analyzer[J]. Chinese Journal of Geophysics−Chinese Edition, 2010, 53(11): 2535-2543
    [13] Fisher L E, Lynch K A, Fernandes P A, et al. Including sheath effects in the interpretation of planar retarding potential analyzer’s low-energy ion data[J]. Review of Scientific Instruments, 2016, 87: 043504. doi: 10.1063/1.4944416
    [14] Pfaff F, Borovsky E, Young T. Measurements of thermal ion drift velocity and temperature using planar sensors[J]. Geophysical Monograph, 1998, 102: 61-71.
    [15] COMSOL Inc. Particle tracing module updates[EB/OL]. https://www.comsol.com/release/5.3/particle-tracing-module.
    [16] Bilitza D, Altadill D, Reinisch B, et al. The international reference ionosphere: model update 2016[C]//EGU General Assembly Conference Abstracts. 2016: 18.
    [17] The Community Coordinated Modeling Center. International Reference Ionosphere−IRI (2016) with IGRF-13 coefficients [EB/OL]. https://ccmc.gsfc.nasa.gov/modelweb/models/iri2016_vitmo.php.
    [18] Marrese C M, Majumdar N, Haas J, et al. Development of a single-orifice retarding potential analyzer for Hall thruster plume characterization[C]//25th International Electric Propulsion Conference. 1997: 397-404.
    [19] Zheng Xiangzhi, Zhang Aibing, Guan Yibing, et al. Research on retarding potential analyzer aboard China seismo-electromagnetic satellite[J]. Acta PhysicaSinica, 2017, 66: 079401.
    [20] 熊年禄, 唐存琛, 李行健. 电离层物理概论[M]. 武汉: 武汉大学出版社, 1999: 53-55.

    Xiong Nianlu, Tang Chunshen, Li Xingjian. Introduction to ionospheric physics[M]. Wuhan: Wuhan University Press, 1999: 53-55
    [21] Lieberman M A, Lichtenberg A J. Principles of plasma discharges and materials processing[M]. New Jersey: John Wiley & Sons, 2005.
  • 期刊类型引用(2)

    1. 韩星,王永琴,曾娅秋,刘宇,粟嘉伟,林珑君. 双极晶体管空间辐射效应的研究进展. 环境技术. 2024(07): 188-194 . 百度学术
    2. 曾超,许献国,钟乐. 抗辐射电子学研究综述. 太赫兹科学与电子信息学报. 2023(04): 452-471 . 百度学术

    其他类型引用(3)

  • 加载中
图(12) / 表(1)
计量
  • 文章访问数:  1302
  • HTML全文浏览量:  396
  • PDF下载量:  78
  • 被引次数: 5
出版历程
  • 收稿日期:  2019-12-09
  • 修回日期:  2020-06-08
  • 刊出日期:  2020-08-13

目录

/

返回文章
返回