Thermal Stress Structural Optimization of InSb Infrared Focal Plane Array Detector
-
摘要: 在液氮冲击实验中,锑化铟红外焦平面阵列探测器中各层材料之间线膨胀系数的不同将导致热失配产生,过大的热失配应力将引起锑化铟芯片断裂失效。为了降低热失配对锑化铟芯片的影响,基于弹性多层体系热应力计算理论,借鉴平衡复合物结构设计方法,优化平衡复合物结构上表面的热应变,使得平衡复合物结构中硅读出电路上表面的热应变尽可能接近锑化铟芯片下表面的热应变,从而大幅降低锑化铟芯片中的热应力。考虑器件加工工艺成熟度,经一系列计算表明:当硅读出电路的厚度取25 μm时,平衡复合物结构中硅读出电路上表面的热应变与InSb芯片下表面的热应变最为接近,此时锑化铟芯片中的拉应力最小。锑化铟芯片中拉应力的大幅降低,将为消减液氮冲击中锑化铟芯片的碎裂几率提供可以信赖的结构设计方案和实现途径。
-
关键词:
- 红外焦平面阵列探测器 /
- 热失配 /
- 热应变 /
- 平衡复合物结构
Abstract: In liquid nitrogen shock experiments, a thermal mismatch occurs owing to the difference in the linear expansion coefficients of the layered In antimonide infrared focal plane array detector, and the excessive thermal mismatch stress fractures the InSb chip. Based on the calculation theory of the thermal stress suitable for the elastic multilayer system, employing the design method of a balanced composite structure is considered to be effective in reducing the impact of the thermal mismatch on the InSb chip. Accordingly, we optimize the thermal strain on the upper surface of the balanced composite structure. In other words, the optimization involved making the thermal strain on the upper surface of the Readout circuit in the balanced composite structure as close as possible to the thermal strain on the lower surface of the InSb chip. Consequently, the reduced thermal mismatch reduces the thermal stress in the InSb chip. Considering the maturity of the device processing technology, the thickness of the readout circuit is set at 25 μm, which is the thinnest sheet of the readout circuit fabricated in our lab using the chemical mechanical polishing method. For the defined thickness (25 μm) of the readout circuit, the calculation results indicate that the thermal strain on the upper surface of the readout circuit is the closest to the thermal strain on the lower surface of the InSb chip. When these two structures are glued together by the underfill, the tensile stress accumulated in the InSb chip is the smallest. The significant reduction in the tensile stress in the InSb chip provides a reliable structural design scheme and an implementation approach to reduce the fragmentation probability of the InSb chip in the liquid nitrogen impact. -
0. 引言
红外焦平面探测器广泛应用于航空航天、红外遥感、国防、气象、环境、医学及科学仪器等领域[1-4]。在发展较为成熟的红外焦平面探测器中,相较于碲镉汞焦平面探测器(HgCdTe IRFPAs)而言,锑化铟焦平面阵探测器(InSb IRFPAs)的均匀性尤为突出,在3~5 μm波段更大面积的InSb光伏探测器阵列也得到更迅猛的发展,并在红外探测领域占据重要地位。例如,2004年雷神视觉系统公司Raytheon Vision Systems(RVS)生产的InSb探测器阵列规模就已经达到了2048×2048[5]。红外焦平面探测器的结构具有典型的层状特征,由光敏元芯片和硅读出电路通过铟柱阵列经倒装焊工艺连接形成[6-8]。为了抑制背景噪声,提高信噪比,灵敏度高的红外焦平面阵列探测器通常工作于液氮温度(77 K)环境。当探测器由装配温度(300 K)快速降到液氮温度时,面阵列探测器结构中相邻材料间线膨胀系数的不同,将会引起热失配,从而在光敏元芯片中产生热应力,过大的热应力将引起光敏元芯片碎裂[9-11],制约着探测器的成品率。这种失效现象在高性能红外焦平面探测器的生产研制中尤为突出。因此如何降低光敏元芯片中的热应力,提高器件的结构可靠性就成为红外焦平面探测器研制中的一个关键问题。
针对液氮冲击下红外焦平面探测器的碎裂失效问题,相关学者对此进行了深入研究。如中国科学院上海技术物理研究所的陈星[12]等以碲镉汞红外焦平面探测器为研究对象,以探测器结构中的热应力、低温形变关系为基础,对各层材料的结构尺寸、材料参数进行了相应的优化设计和实验验证。赵玉红[13]等以InSb红外焦平面探测器为研究对象,将具有低膨胀系数的混合材料置于探测器中底充胶与硅读出电路之间,降低了探测器芯片上的热应力,实现了探测器结构优化的目的。这些研究大多聚焦于光敏元芯片中张应力的减小,以特定结构、特定参数来阐述热应力减小的效果,在细节及深度方面阐述得不够。为此,本文中我们将从探测器整体结构出发,着眼于降低光敏元芯片下表面与硅读出电路上表面之间的热失配,从而从根本上降低光敏元芯片中的张应力,消除液氮冲击中的光敏元芯片的碎裂几率,提高成品率。
1. 探测器结构及弹性多层体系热应力计算理论
InSb IRFPAs是由InSb芯片和硅读出电路经铟柱阵列连接形成的,InSb芯片位于最上层,用来接收红外信号并将其转化为电信号,随后电信号通过铟柱阵列传输给硅读出电路(Si ROIC),紧挨着InSb芯片的材料是底充胶(Underfill),呈网状结构环绕着铟柱阵列,起到机械支撑作用,并隔绝电信号传输通道,将上述3层结构用低温胶粘结在蓝宝石(Sapphire)基片上,方便电信号传输至探测器组件外部。此外,在蓝宝石基片下表面粘结有柯伐(Kovar)合金,方便制冷器冷却InSb IRFPAs。具体结构如图 1所示。
弹性多层体系热应力计算理论是由C. H. Hsueh提出的,理论模型如图 2所示。即在高温下,n层厚度为ti的薄膜依次生长到厚度为ts的基板上,下标i表示膜层标号,范围为1~n,第一层直接沉积在基板上,详见图 2。
坐标系是这样建立的,衬底与第一层薄膜之间的界面位于z=0处,衬底的下自由表面位于z=-ts处,最外(上)层的自由表面位于z=hn处,第i层与第i+1层的界面位于hi,则hi与ti之间的关系表示为:
$$ h_{i}=\sum\limits_{j=1}^{i} t_{j}, (i=1 \sim n) $$ (1) 当多层弹性系统从生长温度冷却到室温时,其降温范围用ΔT表示,假定各层薄膜之间没有相互束缚,降温后,衬底和薄膜中的自由热应变分别为εs=αsΔT、εi=αsΔT。
多层弹性系统中衬底和薄膜沿z轴方向的平面内正应力σ与系统的总应变ε之间的关系为:
$$ \sigma_{\mathrm{s}}=\frac{E_{\mathrm{s}}}{1-v_{\mathrm{s}}}\left(\varepsilon_{\mathrm{sT}}-\alpha_{\mathrm{s}} \Delta T\right), t_{s} \leq z \leq 0 $$ (2) $$ \sigma_{i}=\frac{E_{i}}{1-v_{i}}\left(\varepsilon_{i \mathrm{~T}}-\alpha_{i} \Delta T\right), 1 \leq i \leq n $$ (3) 式中:Es、Ei分别为衬底和第i薄膜的杨氏模量;νs、νi分别表示衬底和第i薄膜的泊松比;εsT,εiT分别表示衬底和第i薄膜在系统中的表现应变,T为total。
多层弹性系统中总应变ε可以分解为均匀应变分量c和弯曲应变分量,弯曲应变分量用(z-tb)/r来计算。显然,某一位置处的弯曲应变分量与该位置距离弯曲轴的距离(z-tb)成正比,与系统的曲率半径r成反比,弯曲轴tb定义为系统的弯曲应变分量为零时的轴线位置。由此,系统的总应变量ε可表述为:
$$ \varepsilon=c+\frac{z-t_{\mathrm{b}}}{r}, -t_{\mathrm{s}} \leq z \leq h_{n} $$ (4) 当弹性多层体系不发生弯曲变形时,弹性多层体系的总应变量为均匀应变分量c,由均匀应变分量引起的和应力为零;当弹性多层系统仅发生弯曲变形时,系统中的弯曲应变分量产生的合应力为零;相对于弯轴位置tb,弹性多层体系的弯矩之和为零,由上述3个边界条件求解可得到c,tb,1/r三个系统参数,用等式(5)~(7)描述:
$$ c= \frac{\left(\frac{E_{\mathrm{s}}}{1-v_{\mathrm{s}}} \alpha_{\mathrm{s}} t_{\mathrm{s}}+\sum\limits_{i=1}^{n} \frac{E_{i}}{1-v_{i}} \alpha_{i} t_{i}\right) \Delta T}{\frac{E_{\mathrm{s}}}{1-v_{\mathrm{s}}} t_{\mathrm{s}}+\sum\limits_{i=1}^{n} \frac{E_{i}}{1-v_{i}} t_{i}} $$ (5) $$ t_{\mathrm{b}}= \frac{-\frac{E_{\mathrm{s}}}{1-v_{\mathrm{s}}} t_{\mathrm{s}}^{2}+\sum\limits_{i=1}^{n} \frac{E_{i}}{1-v_{i}} t_{i}\left(2 h_{i-1}+t_{i}\right)}{2\left(\frac{E_{\mathrm{s}}}{1-v_{\mathrm{s}}} t_{\mathrm{s}}+\sum\limits_{i=1}^{n} \frac{E_{i}}{1-v_{i}} t_{i}\right)} $$ (6) $$ \frac{1}{r}=\frac{3\left[\frac{E_{\mathrm{s}}}{1-v_{\mathrm{s}}}\left(c-\alpha_{\mathrm{s}} \Delta T\right) t_{\mathrm{s}}^{2}\right]-\sum\limits_{i=1}^{n} \frac{E_{i}}{1-v_{\mathrm{i}}} t_{i}\left(c-\alpha_{i} \Delta T\right)\left(2 h_{i-1}+t_{i}\right)}{\frac{E_{\mathrm{s}}}{1-v_{\mathrm{s}}} t_{\mathrm{s}}^{2}\left(2 t_{\mathrm{s}}+3 t_{\mathrm{b}}\right)+\sum\limits_{i=1}^{n} \frac{E_{i}}{1-v_{\mathrm{i}}} t_{i}\left[6 h_{i-1}^{2}+6 h_{i-1} t_{i}+2 t_{i}^{2}-3 t_{\mathrm{b}}\left(2 h_{i-1}+t_{i}\right)\right]} $$ (7) 根据弹性多层热应力计算理论,在系统降温范围ΔT一定时,根据公式(3)可知,第i层薄膜的面内正应力仅取决于该层膜在系统中的总应变,系统的总应变用公式(4)描述,仅与c,tb,1/r三个参数有关,这3个参数的大小仅取决于材料的杨氏模量、泊松比、材料的线膨胀系数和厚度。当材料确定后,唯一能调整的就是材料的厚度了。故此,我们认为弹性多层系统第i层薄膜中的面内正应力的大小可通过层状材料厚度的调整而减小。对于InSb IRFPAs而言,为了降低InSb芯片受到的张应力,使硅读出电路上表面的热应变尽可能接近InSb芯片下表面的热应变,用于减小二者之间的热失配,从而实现探测器结构的优化。
2. 平衡复合物结构及热应变
由弹性多层系统热应力计算理论可知,在材料选取一定的前提下,平衡复合物结构的热应变只与其厚度有关,故此,为了使硅读出电路上表面的热应变尽可能接近InSb芯片下表面的热应变,需要对底充胶下面的多层材料的厚度进行优化调整。工程上,常用的多层材料的选取与排布如图 3所示,该结构统称为平衡复合物结构[14]。
在平衡复合物结构中,从上到下依次为硅读出电路、蓝宝石基片和柯伐合金。这3层材料中,硅读出电路的线膨胀系数与InSb芯片的线膨胀系数的差距最大,调整效果应该最明显。
3. 探测器的模型建立及参数选择
为了解决InSb面阵探测器的热失配问题,降低InSb芯片中的张应力,基于弹性多层体系热应力计算理论,借鉴平衡复合物结构的设计思路,优化探测器结构,使硅读出电路上表面的热应变尽可能接近InSb芯片下表面的热应变,从而达到降低InSb芯片中张应力的目的。根据弹性多层体系热应力计算理论[15-16]可知,在材料选取一定的前提下,平衡复合物结构的热应变只与其厚度有关,故探测器层状结构中的材料布局和厚度选择就成了实现结构优化的主导因素。
InSb焦平面探测器通常由5种层状材料构造而成,最上层到最下层依次是InSb芯片(10 μm)、底充胶(10 μm)、硅读出电路(25~340 μm)、蓝宝石基片(500 μm)和柯伐合金(635 μm)[17]。在InSb焦平面探测器结构中,将InSb芯片、固化的底层胶、硅读出电路、蓝宝石基片和柯伐合金都视为线弹性材料,参数如表 1所示。InSb面阵探测器的加工过程在300 K完成,在这个温度下,可以假设InSb焦平面探测器内部不存在残余应力。在整个InSb焦平面探测器的降温过程为瞬态降温,可以假设InSb IRFPAs内部的温度是均匀的。当温度从300 K快速降低到77 K,由于各材料线膨胀系数不同,导致热应力在InSb焦平面探测器内部产生。
表 1 InSb IRFPAs结构的材料参数及其具体尺寸Table 1. Parameters of InSb IRFPAsMaterials Elastic modulus (E/GPa) Poison's ratio (μ) Temperatures (T/K) Height/μm Coefficient of linear expansion(10-6K-1) InSb chip 409 0.35 77–300 10 5.1 Underfill 0.0002/α 0.30 77–300 10 23.82 Silicon ROIC 163 0.28 77–300 25–340 2.62 Sapphire 335 0.25 77–300 500 4.8 Kovar 142 0.32 77–300 635 5.5 Note: α=22.46×10-6+5.04×10-8×(T-273), where T is temperature in Kelvin 4. 结果与讨论
为了能够降低InSb芯片中的张应力,调整复合材料组件中硅读出电路的厚度,使硅读出电路上表面的热应变更为接近InSb芯片的自由应变,并以此探究结构优化对InSb芯片中张应力的影响;将上述参数代入复合材料组件模型中进行计算,可得硅读出电路厚度取不同值时,复合材料组件中硅读出电路上表面的热应变变化趋势,如图 4。显然,当硅读出电路的厚度为25 μm时,InSb IRFPA中硅读出电路上表面的热应变与InSb芯片的热应变最为接近。
在液氮冲击中,InSb芯片的自由应变为-1137.3 PPM(part per million,PPM),当硅读出电路从340 μm逐渐减薄到25 μm时,复合材料组件中硅读出电路上表面的热应变由-695 PPM逐渐变化到-1004 PPM,不断接近InSb芯片的自由应变。
根据前述公式(3)~(7)可知:如果将硅片减薄后的BCS结构与InSb芯片连接在一起,当硅读出电路较薄时,复合材料组件中硅读出电路上表面的热应变更接近InSb芯片的自由热应变,预期将InSb芯片与复合材料组件结合后,InSb芯片与复合材料组件之间的热失配更小,从而能够显著降低InSb芯片中的张应力。
为此,我们将硅读出电路厚度取不同值时,整个探测器中的热应力计算出来,一同绘制在图 5中。显然当硅读出电路的厚度为25 μm时,InSb芯片受到的张应力最小。
如图 5,当硅读出电路从340 μm变薄到25 μm,InSb芯片的张应力从263 MPa下降到78 MPa,由弹性多层热应力计算理论可知,InSb芯片中的热应变随着硅读出电路厚度的减小而减小,InSb芯片中的张应力亦随之有所降低。硅读出电路中的应力(负值)为压应力,InSb芯片中的应力(正值)为张应力,都随着硅读出电路厚度的变化而变化,且硅读出电路与InSb芯片之间存在相对的受力关系,由于硅读出电路的线膨胀系数与InSb芯片的线膨胀系数的差距最大,间接控制了InSb芯片的受力程度,故调整硅读出电路的厚度对改善InSb芯片的应力水平有很大的帮助。
在液氮冲击下,由于蓝宝石的线膨胀系数大于硅读出电路的线膨胀系数,导致二者形成的结构呈现中间凸起的状态,以硅读出电路为25 μm为例,由计算可知,二者形成的结构的曲率1/r为0.1399 m-1;这时在蓝宝石基片下粘结上线膨胀系数比蓝宝石大的柯伐材料,计算可得三者所形成结构的曲率为0.2365 m-1,促使结构中心更加的凸起,进一步压缩硅读出电路,在InSb芯片中产生压应力,抵消芯片受到的部分张应力,从而提高探测器的可靠性。
T. Kanno等[18]人为减少红外焦平面探测器中的热失配,采用了硅读出电路减薄的方法,以256×256 HgCdTe红外焦平面探测器为例,为了确定最优厚度的硅读出电路,实验测量了在热冲击下减薄后的硅读出电路以及蓝宝石基片沿对角线的收缩长度,并绘制了对角线收缩长度与硅读出电路厚度的关系图。实验表明,当硅读出电路的厚度小于40 μm时,收缩逐渐减少,直到接近蓝宝石基片的收缩,并且当硅读出电路的厚度减薄到25 μm时,硅读出电路的热收缩程度最为接近HgCdTe芯片的热收缩程度,降低了硅读出电路与HgCdTe芯片的热失配,使得HgCdTe芯片中的张应力有所降低。该实验与本文的计算结果基本一致,因此,调整硅读出电路的厚度是降低硅读出电路与芯片之间的热失配,实现探测器结构优化的一种有效的方法。
5. 结论
为了降低InSb焦平面探测器在液氮冲击中产生的热应力,基于弹性多层热应力计算理论,借鉴平衡复合物结构设计方法,对InSb面阵列探测器进行结构优化,使得平衡复合物结构中硅读出电路上表面的热应变尽可能接近锑化铟芯片下表面的热应变,从而达到降低锑化铟芯片中张应力的目的。结果表明当硅读出电路厚度为25 μm时,复合材料组件中硅读出电路上表面的热应变最为接近InSb芯片的自由应变,并且在该厚度时,InSb芯片受到的张应力最小,提高了探测器的可靠性。
-
表 1 InSb IRFPAs结构的材料参数及其具体尺寸
Table 1 Parameters of InSb IRFPAs
Materials Elastic modulus (E/GPa) Poison's ratio (μ) Temperatures (T/K) Height/μm Coefficient of linear expansion(10-6K-1) InSb chip 409 0.35 77–300 10 5.1 Underfill 0.0002/α 0.30 77–300 10 23.82 Silicon ROIC 163 0.28 77–300 25–340 2.62 Sapphire 335 0.25 77–300 500 4.8 Kovar 142 0.32 77–300 635 5.5 Note: α=22.46×10-6+5.04×10-8×(T-273), where T is temperature in Kelvin -
[1] 雷亚贵, 王戎瑞, 陈苗海. 国外非制冷红外焦平面阵列探测器进展[J]. 激光与红外, 2007(9): 801-805. DOI: 10.3969/j.issn.1001-5078.2007.09.001 LEI Yagui, WANG Rongrui, CHEN Miaohai. Development of foreign uncooled IRFPA detectors[J]. Laser & Infrared, 2007(9): 801-805. DOI: 10.3969/j.issn.1001-5078.2007.09.001
[2] Bhan R K, Dhar V. Recent infrared detector technologies, applications, trends and development of HgCdTe based cooled infrared focal plane arrays and their characterization[J]. Opto-Electronics Review, 2019, 27(2): 174-193. DOI: 10.1016/j.opelre.2019.04.004
[3] Philip Klipstein, Daniel Aronov, Michael ben Ezra, et al. Recent progress in InSb based quantum detectors in Israel[J]. Infrared Physics & Technology, 2013, 59: 172-181. http://www.onacademic.com/detail/journal_1000035936383110_429d.html
[4] Rogalski A. Recent progress in infrared detector technologies[J]. Infrared Physics & Technology, 2011, 54(3): 136-154. http://www.researchgate.net/profile/Antoni_Rogalski/publication/241112884_Recent_progress_in_HgCdTe_infrared_detector_technology/links/5524e12e0cf22e181e73b04e.pdf
[5] Hoffman A W, Corrales E, Love P J, et al. 2K×2K InSb for astronomy[J]. Proceedings of SPIE - The International Society for Optical Engineering, 2004, 5499: 59-67.
[6] HA H D, TANG Z, LU H D, et al. Design and simulation of a novel flip-chip structure for THz detector package[J]. Microsystem Technologies, 2016, 23(7): 1-8. DOI: 10.1007/s00542-016-3021-1
[7] MENG Q D, ZHANG X L, LV Y Q, et al. Function reconsideration of indium bump in InSb IRFPAs[J]. Optical and Quantum Electronics, 2019, 51(9): 304.1-304.13. DOI: 10.1007/s11082-019-2021-7
[8] 陈星, 华桦, 何凯, 等. 红外焦平面探测器封装结构热应力分析[J]. 激光与红外, 2014(6): 59-62. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201406013.htm CHEN Xing, HUA Hua, HE Kai, et al. Thermal stress analysis of IRFPA packaging assembly[J]. Laser & Infrared, 2014(6): 59-62. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201406013.htm
[9] 王雯, 张小雷, 吕衍秋, 等. Si基InSb红外焦平面阵列探测器的研究[J]. 红外与激光工程, 2014(5): 1359-1363. DOI: 10.3969/j.issn.1007-2276.2014.05.001 WANG Wen, ZHANG Xiaolei, LV Yanqiu, et al. InSb infrared focal plane arrays detector based on Si wafer[J]. Infrared and Laser Engineering, 2014(5): 1359-1363. DOI: 10.3969/j.issn.1007-2276.2014.05.001
[10] ZHANG X L, MENG Q D, ZHANG L W, et al. Modeling and deformation analyzing of InSb focal plane arrays detector under thermal shock[J]. Infrared Physics & Technology, 2013, 63: 28-34.
[11] MENG Q D, ZHANG X L, LV Y Q, et al. Local delamination of InSb IRFPAs in liquid nitrogen shock tests[J]. Infrared Physics & Technology, 2017, 86: 207-211.
[12] 陈星. 碲镉汞红外焦平面探测器可靠性相关技术研究[D]. 上海: 中国科学院研究生院(上海技术物理研究所), 2014. CHEN Xing. Research on the Related Technology of HgCdTe Infrared Focal Plane Detector Reliability[D]. Shanghai: Shanghai Institute of Technology Physics, Chinese Academy of Sciences, 2014.
[13] 赵玉红. InSb红外焦平面探测器结构优化研究[D]. 洛阳: 河南科技大学, 2016. ZHAO Y H. Optimization of Structure of InSb Infrared Focal Plane Array[D]. Luoyang: Henan University of Science and Technology, 2016.
[14] Kuzmenko N K. BCS universal ratios in finite systems[J]. Physica C: Superconductivity and its Applications, 2020, 576: 1353709. http://www.sciencedirect.com/science/article/pii/S0921453420300332
[15] Hsueh, C H. Thermal stresses in elastic multilayer systems[J]. Thin Solid Films, 2002, 418: 182-188. http://www.onacademic.com/detail/journal_1000034583558710_41a3.html
[16] Hsueh, C H, Luttrell C R, Lee S, et al. Interfacial peeling moments and shear forces at free edges of multilayers subjected to thermal stresses[J]. J. Am. Chem. Soc. , 2006, 89: 1632-1638. http://www.researchgate.net/profile/Chun-Hway_Hsueh/publication/229870518_Interfacial_Peeling_Moments_and_Shear_Forces_at_Free_Edges_of_Multilayers_Subjected_to_Thermal_Stresses/links/00b7d523d4a135cc8a000000
[17] ZHANG X L, MENG Q D, LV Y Q, et al. Assessment of InSb infrared detector arrays assembly procedure employing ANSYS[J]. Optical and Quantum Electronics, 2019, 51(4): 1-8. http://www.onacademic.com/detail/journal_1000042271909399_7d60.html
[18] Kanno T, Wada H, Nagashima M, et al. 256×256 element HgCdTe hybrid IRFPA for 8- to 10-μm band[C]//Proceedings of SPIE-The International Society for Optical Engineering, 1995, 2552: doi: 10.1117/12.218237.
-
期刊类型引用(1)
1. 方志浩,付志凯,王冠,张磊. 超大面阵碲镉汞探测器低应力设计及有限元分析. 红外. 2024(02): 9-17 . 百度学术
其他类型引用(1)