10.6mm激光輻照碲鎘汞紅外探測器熱損傷研究

栗興良，牛春暉，馬牧燕，呂 勇,耿 蕊

?

10.6mm激光輻照碲鎘汞紅外探測器熱損傷研究

栗興良，牛春暉，馬牧燕，呂 勇,耿 蕊

（北京信息科技大學 儀器科學與光電工程學院，北京 100192）

研究了HgCdTe紅外探測器的結構以及材料特性，闡述了激光損傷HgCdTe紅外探測器的機理，建立了HgCdTe紅外探測器三維仿真模型，利用有限元分析法，對10.6mm CO2激光輻照HgCdTe探測器的溫度變化情況進行了仿真，并通過參考已有文獻的實驗數據，驗證了模型的準確性。當HgCdTe探測器受到峰值功率密度為5×107W/cm2的單脈沖激光輻照時，HgCdTe晶體的Hg離子開始析出，探測器性能降低，并不可恢復；當激光峰值功率為108W/cm2，探測器HgCdTe晶體開始出現熔融現象，此時激光能量密度為1J/cm2；當激光峰值功率為2×108W/cm2時，銦柱達到熔融溫度，探測器會出現銦柱脫落現象，被徹底損壞。

HgCdTe；紅外探測器；單脈沖激光；激光損傷；有限元分析法

0 引言

碲鎘汞材料的發明，揭開了紅外科技的新篇章，使紅外技術進入了嶄新的時代。碲鎘汞卓越的性能，是其他紅外材料沒法相比的。經過半個多世紀的發展，碲鎘汞探測器經歷了3個階段：第一代陸軍通用組件，第二代高性能紅外焦平面器件和第三代高清、多波段紅外焦平面器件[1]。在碲鎘汞材料發明前，紅外探測主要采用非本征半導體來實現，例如鍺摻雜非本征紅外探測器，這類探測器需要大幅度降低溫度來抑制材料固有的高熱激發噪聲，工作溫度一般低于28K，實際運用難度很大，制造成本也很高。碲鎘汞材料則很大程度上解決了這些難題，碲鎘汞探測器的研發經歷了艱辛的歷程，要把碲、鎘、汞3種元素組合在一起，需要很高的汞蒸汽壓，巨大的汞壓讓研究設備發生了很多次爆炸，研究人員不得不降低溫度，這樣就直接影響了材料的機械性能。直到研究人員認識到用CdTe和HgTe兩種化合物按比例混合來合成碲鎘汞晶體，碲鎘汞的研究道路才平坦了起來。碲鎘汞材料吸收光子的高效性決定了其受到高能激光輻照之后溫度會迅速升高，探測效果受到干擾，甚至永久性損傷。所以研究激光對碲鎘汞探測器熱損傷就非常有必要。碲鎘汞探測器市場價格非常高，決定了大部分情況無法通過大量實驗來研究其損傷特性，計算機仿真可以解決這一難題，那么設計相對準確的仿真模型成為一種需求。

1 碲鎘汞探測器特性和結構

1.1 碲鎘汞材料特性

碲鎘汞化學式Hg1－CdTe有量子效率高、響應速度快的特點，可以通過改變鎘組分（）讓它的光電響應覆蓋1～3mm、3～5mm、8～12mm三個紅外“大氣窗口”，以及18mm以上的甚長波紅外波段，碲鎘汞材料在碲鎘汞探測器中以外延薄膜材料應用為主。

HgCdTe材料在80K時對于10.6mm波長的介電常數實部為3.54，虛部為0.24[2]；295K時，實部為3.54，虛部為0.78；CdTe介電常數實部為2.6實部為0[3]。

HgCdTe材料的導熱系數為：6.394041－0.020149591＋5.954089×10－52－1.134691×10－73＋8.148373×10－114，為熱力學溫度。

1.2 碲鎘汞探測器模型

一代碲鎘汞探測器為光導線列，一維；二代為光伏陣列器件，二維，如240×4，288×4，480×4等；三代則為（二維，光伏）大面陣凝視器件，如640×480，624×512等[4]。光導器件結構相對簡單，為第一代探測器，光伏器件較光導器件功耗小，響應快，但結構復雜。現在只分析其熱特性，故將其結構簡化。探測器分為4層，首先在CdTe層上生長一層HgCdTe晶體作為光敏元，然后焊接銦柱陣列，最下邊是硅基。其中CdTe層厚50mm，HgCdTe光敏層厚10mm，In柱高10mm，直徑20mm，Si基厚50mm，整個模型為11×11個像素元。第三代碲鎘汞紅外探測器工作溫度為77K，需要液氮來冷卻，模型采用將表面溫度設置為恒溫77K的方式，來模擬液氮冷卻場。

圖1 探測器模型結構

碲鎘汞是一種主要由離子鍵結合的三元半導體材料，構成元素汞非常不穩定，汞容易從碲鎘汞材料中析出從而造成材料缺陷、材料不均勻以及器件性能不均勻[5]。實驗表明碲鎘汞材料在350K左右的時候即會產生Hg析出現象，破壞HgCdTe晶體結構，影響器件性能[6]；銦柱熔點為426K，所以當銦柱溫度達到426K時，銦柱就有可能脫落；HgCdTe熔點約為990K，當HgCdTe層溫度達到熔點時，探測器將被徹底破壞。

1.3 激光加熱理論

當激光照射到材料表面時，主要以光熱模型為主[7]，一部分激光能量被反射，另一部分能量被激光作用區材料吸收，在這個過程中，激光束等效為一個具有一定時間和空間分布的熱源。對于材料吸收的功率密度，我們可以根據玻意耳定律表示：

式中：()為距表面處單位體積材料吸收的輻射功率；0為材料表面接收到的激光功率；1－為材料吸收率；為材料的吸收系數，cm－1。

熱量傳遞分為導熱、對流、輻射3部分，它們的機理也是各不相同的，方程也不一樣。總的傳熱方程為：

式中：是材料密度；p是材料比熱容；是絕對溫度；是速度矢量；?為幾何梯度算子；為材料導熱系數；()為激光熱源；等號左邊第一項為熱量積累，第二部分為對流項，等號右邊第一項為熱傳導項。

采用三節點三角形對模型進行網格剖分，三角形單元對復雜幾何輪廓有較強的適應能力，很容易通過增加三角形網格數量來精確的逼近復雜的幾何邊界[8]。

1.4 有限元理論

有限元方法，基礎是變分原理和加權余量法，基本思路是把求解域離散成有限個互不疊加的單元，然后在每個單元內選擇合適的節點作為插值點，把待求的偏微分方程中的因變量改寫成依據及節點上的值的插值函數組成的線性方程組[9]，從而求得所需要的解。有限元方法能夠處理復雜區域，精度可根據需要調節，缺點是計算量大，但高性能計算機已經解決了這個問題。

2 激光輻照碲鎘汞探測器仿真

2.1 波長為10.6mm、功率密度為106 W/cm2激光輻照探測器

不同強度激光輻照碲鎘汞探測器產生的結果有很大的差別，以功率密度為106W/cm2為例。激光脈寬為10ns，仿真時間為200ms，假設在光斑范圍內，激光能量分布均勻。設置探測器上接收到的激光光斑半徑為15mm，功率密度為106W/cm2。經過計算，CdTe層和HgCdTe層吸收的激光能量密度分別如圖2所示，橫坐標為弧長，縱坐標為吸收能量密度。碲鎘汞層和銦柱層溫度分別如圖3所示，橫坐標為時間，縱坐標為相應溫度。

圖2 CdTe層和HgCdTe層對激光的吸收

圖3 碲鎘汞層和銦柱層溫度

我們可以看出，碲鎘汞探測器受到10.6mm激光輻照時，CdTe層對激光能量幾乎沒有吸收，而HgCdTe晶體對激光能量的吸收達到93%。高效的激光能量吸收效率使HgCdTe材料受激光輻照時溫度迅速上升，從如圖3我們可以看出碲鎘汞光敏層吸收激光能量之后，溫度迅速上升，在極短的時間內達到84K左右，然后溫度上升趨勢變緩，在1.5ms左右達到最大值85.5K，然后緩慢下降。

開始溫度迅速上升是因為激光脈沖輻照在探測器上面，探測器吸收了激光能量而發生的能量轉換；然后溫度升高變緩是由于溫度升高之后，熱傳遞和熱輻射作用增強；最后緩慢降低，是因為激光脈沖結束，探測器內部液氮冷卻和熱傳遞作用導致。最后由于液氮制冷HgCdTe光敏層溫度降為77K。

銦柱溫度上升所需要的能量基本上是由于HgCdTe層熱傳遞而來。銦柱溫度在10ms時間內持續上升，上升速度越來越緩慢，在12ms的時候溫度停止升高，最大溫度80K左右，由于制冷作用此后開始降低，最后降低為77K。

圖4 探測器溫度分度圖

模型整體分布溫度如圖4所示，從圖中可以看出，GeTd層溫度較低，距離光斑較遠的銦柱溫度也較低；激光輻照區域的碲鎘汞層溫度較高，銦柱次之，距離光斑越遠，溫度越低，符合熱力學理論。

2.2 不同功率密度的單脈沖激光輻照探測器

以不同強度的單脈沖激光輻照碲鎘汞探測器進行仿真并記錄實驗數據。激光波長為10.6mm，脈寬為10ns。激光功率密度依次為：1×106W/cm2、1×107W/cm2、3×107W/cm2、4×107W/cm2、8×107W/cm2、1×108W/cm2、1×109W/cm2。分別計算銦柱和HgCdTe層的溫度，其中3×107W/cm2和3×107W/cm2的溫度曲線如圖5所示，其余數據經進一步整理成表格形式如表1所示。

圖5 不同功率密度輻照下碲鎘汞層和銦柱最高溫度曲線

表1 不同功率密度輻照下碲鎘汞層和銦柱最高溫度以及激光能量密度

從圖5和表1可以看出來，當激光功率密度為3×107W/cm2時，探測器溫度最高到340K，并未達到碲鎘汞晶體Hg離子析出和晶體熔化溫度，此時探測器只受到軟損傷，停止激光照射，探測器的性能還可以恢復。也就是說，當激光功率密度小于此閾值時，激光無法對碲鎘汞焦平面造成硬損傷。激光功率密度增大到5×107W/cm2時，碲鎘汞層可以溫度升高到500K，此時HgCdTe晶體會有Hg元素析出，探測器性能降低，并且無法恢復。激光功率密度增大到108W/cm2時，HgCdTe溫度超過1000K，達到了HgCdTe材料的融化溫度，此時探測器受到熔融性損壞，無法修復。當激光功率密度增大到2×108W/cm2時，銦柱的溫度達到了450K，此時銦柱達到熔點，出現銦柱脫落現象，探測器徹底損壞。激光脈寬為10ns，可以計算出HgCdTe發生熔融性損壞時，激光能量要大于1J/cm2。

HgCdTe層最高溫度隨功率密度變化如圖6所示。從圖中可以明顯的看出來，激光功率和HgCdTe光敏層最高溫度大致成線性關系，隨著激光功率密度的增大，HgCdTe層可以達到的最高溫度也隨之增大。

圖6 HgCdTe層最高溫度隨功率密度變化

2.3 仿真數據與實驗對比

華中科技大學的蔡虎等人用CO2單脈沖激光器輻照碲鎘汞晶片，當激光峰值功率密度達到1.67×108W/cm2時，HgCdTe晶片表面出現凹凸不平現象，顯示融化過。同時出現網狀龜裂和長的裂縫[10]。仿真的出的損傷功率密度閾值與之相符。

四川大學和中國人民解放軍63891部隊聯合用10.6mm脈沖CO2激光輻照PC型HgCdTe探測器進行破壞閾值試驗研究[11]，當激光脈寬為50ns時，探測器的破壞閾值為12.38J/cm2，此時探測器的輸出電壓和電阻均增大，仿真的能量損傷閾值和實驗基本相符。

3 結語

本文針對波長為10.6mm的激光輻照HgCdTe紅外探測器產生的損傷效應進行了仿真計算，得出激光功率密度在108W/cm2級別時，便能夠對CCD探測器造成永久性損傷，若激光脈寬為10ns，則激光能量密度為1J/cm2，仿真數據量級與實驗基本相符，變化趨勢與實驗一致。本仿真模型中加入了液氮的冷卻，但是未考慮脈沖激光的爆轟波效應，使計算結果略有偏差，有待進一步改進。

[1] 曾戈虹. HgCdTe艱辛的歷程、輝煌的成就: 一、偉大的發明[J]. 紅外技術, 2011, 33(5): 249-251.

ZENG Gehong.The arduous journey and brilliant achievements of HgCdTe: 1. great invention[J]., 2011, 33(5): 249-251

[2] Edward D. Palik. Handbook of Optical Constants of Solids Ⅱ[Z]. San Diego: academic press.1998: 682

[3] Edward D. Palik. Handbook of Optical Constants of Solids Ⅰ[Z]. San Diego: academic press. 1998: 413.

[4] 王憶鋒, 唐利斌. 第三代碲鎘汞器件的研發進展[J]. 光電技術應用, 2009(5): 17-22+66.

WANG Yi-feng, TANG Li-bin.Advances in third-generation HgCdTe devices[J].. 2009(5): 17-22+66.

[5] 史衍麗. 第三代紅外探測器的發展與選擇[J]. 紅外技術, 2013, 35(1): 1-8.

SHI Yan-li.Choice and Development of the third-generation infrared detectors[J].. 2013, 35(1): 1-8.

[6] 蔣志平, 梁天驕, 陸啟生, 等. 激光輻照PC型HgCdTe探測器熱效應的計算[J]. 應用激光, 1995(4): 155-156.

JIANG Zhiping, LIANG Tianjiao, LU Qisheng, et al.Heat effect calculations of PC type HgCdTe detectors when irradiated by Laser[J]., 1995(4): 155-156.

[7] 張震, 江天, 程湘愛, 等. CCD強光串擾效應的串擾線缺口現象及其機制[J]. 強激光與粒子束, 2010, 22: 1505-1510.

ZHANG Zhen, JIANG Tian, CHENG Xiang’ai, et al. Investigation into the dynamical subordinate facula in the video of the CCD irradiated by high repetition pulses laser[J]., 2010, 22: 1505-1510.

[8] 王昂, 郭鋒, 朱志武, 等. 連續激光與單脈沖納秒激光對CMOS的損傷效應[J].強激光與粒子束, 2014, 26(9): 43-47.

WANG Ang, GUO Feng, ZHU Zhiwu, et al. Comparative study of hard CMOS damage irradiated by CW laser and single-lulse ns laser[J]., 2014, 26(9): 43-47.

[9] 王剛, 安琳. COMSOL Multiphysics工程實踐與理論仿真[M]. 北京: 電子工業出版社, 2012: 22-23.

WANG Gang, AN Lin.[M]. Beijing: Publishing house of electronics industry, 2012: 22-23.

[10] 蔡虎. TEA-CO2脈沖強激光對Hg0.8Cd0.2Te晶片的破壞研究[D]. 武漢: 華中科技大學, 2005.

CAI Hu. Investigation of Hihg Intensity Pulsed TEA-CO2Laser-induced Damage on Hg0.8Cd0.2Te Wafer[D]. Wuhan: Huazhong University of Science and Technology, 2005.

[11] 張紅, 薛建國, 成斌, 等. 10.6μm CO2激光對HgCdTe探測器破壞閾值的實驗研究[J]. 光電工程, 2006(5): 41-43+56.

ZHANG Hong, XUE Jianguo, CHENG Bin, et al. Experiment study on the damage thresholds of HgCdTe detector irradiated by 10.6μm CO2laser[J]., 2006(5): 41-43+56.

Research on the Thermal Damage of HgCdTe Infrared Detector under Laser Irradiation of 10.6mm Wavelength

LI Xingliang，NIU Chunhui，MA Muyan，LYU Yong，GENG Rui

(,,100192,)

The structure and material properties of the HgCdTe infrared detector was studied, the mechanism of laser damaging HgCdTe infrared detector was described，and the 3D simulation model of HgCdTe infrared detector was built. The temperature changing of HgCdTe detector under 10.6mm laser irradiation was simulated depending finite element analysis method. Accuracy of the model was verified by contrasting experimental data of the reference literatures. The simulation results show when the peak power density of the single pulse laser is 5×107W/cm2, Hg ion will separate from HgCdTe crystal and the properties of HgCdTe infrared detector is damaged and could not be restored. When the peak power density of the single pulse laser is 108W/cm2, HgCdTe crystal of the detector starts to melt, and the laser fluence is 1J/cm2. When the peak power density of the single pulse laser is 2×108W/cm2, indium column reaches the fusion point and will shed, then the HgCdTe infrared detector is completely damaged.

HgCdTe，infrared detector，single-pulse laser，laser damage，finite element method

TN216

A

1001-8891(2016)01-0006-04

2015-07-20；

2015-09-24.

栗興良（1986-），男，山東臨沂人，碩士，主要從事激光應用與光電檢測方向研究。

北京市教委科研計劃項目（KM201511232005）。