目录图表
1 热场设计

1 热场设计

2 实验

2 实验

2.1 样品制备和测试

2.1 样品制备和测试

2.2 双加热器热场模拟

2.2 双加热器热场模拟

3 结果与讨论

3 结果与讨论

3.1 数值模拟结果

3.1 数值模拟结果

3.2 双加热器热场对单晶成晶影响分析

3.2 双加热器热场对单晶成晶影响分析

3.3 电阻率均匀性分析

3.3 电阻率均匀性分析

3.4 光学均匀性测试分析

3.4 光学均匀性测试分析

4 结论

4 结论

参考文献:

参考文献:

Φ420 mm超大直径红外锗单晶研制

张 泽敏1王 宇2汪 洲1田 莉1陈 亚新1杨 业1韩 子洋1韩 智君1马 远飞2
1. 北京空间机电研究所质量技术处,北京 100094
2. 有研国晶辉新材料有限公司,河北 廊坊 065001
摘要: 近年来,为了获取更大的观测视场和更高的灵敏度及分辨率,空间红外光学遥感载荷对红外镜头的直径和均匀性均提出了越来越高的要求。超大直径红外(直径≥420 mm)锗单晶作为大视场、高质量天基红外光电成像系统中的首选材料之一,实际制备过程中仍面临着成晶困难、光学和电阻率均匀性差等问题。为此,本文提出了一种新型双加热器结构,结合CGSim模拟软件,建立了大尺度组合式热场理论模型;重点研究分析了底加热器功率对热场温度分布、单晶成晶及电阻率均匀性的影响规律;最终成功实现了Φ420 mm超大直径锗单晶制备。测试结果表明:采用该热场结构获得的Φ420 mm锗单晶光学均匀性达到6.52×10−5,电阻率均匀性达到19%,均优于使用指标要求。研究结果可为同类型大直径单晶材料制备提供较好的参考借鉴。
关键词: 超大直径;红外锗单晶;电阻率;均匀性;双加热器热场
基金: 河北省科技厅重大科技成果转化专项(22281101Z)
收稿日期: 2024-01-17

在大气窗口3~5 μm中波红外波段和8~14 μm长波红外谱段被广泛用于空间对地军事侦察。空间红外光学遥感载荷[1]对成像视场、成像质量以及能量集中度均有较高要求,光学设计时前端采用大口径零膨胀石英构成的反射系统,后端红外光路还需使用硅、锗(germanium,Ge)等透镜构成二次校正成像系统(图1),实现光学镜头冷光阑匹配和复消色差。其中,锗单晶在红外光学透射材料中具有高折射率(n≈4.0)、低色散(阿贝数υ≈1200)(表1)、各向同性、无本征散射和双折射效应,具有良好的力学性能、理化稳定性和抗辐照性能等特性,有利于简化系统、降低系统像差,是大视场高质量红外成像系统的首选材料之一[2]。

图1 红外光学载荷镜头示意图
Fig.1 Schematic diagram of infrared optical loading lens
表1 常用红外透镜材料性能对比
Table 1 Comparison of properties of common infrared lens materials
Materials Ge Si Zinc selenide Chalcogenide glass
Transmission band/μm 1.8~23 1.1~9 0.5~20 0.9~14
Refractive index ~4.0 ~3.4 ~2.4 ~2.6
υ ~1200 ~240 ~55 ~175
Refractive index temperature coefficient/K−1 3.96×10−4 1.6×10−4 6×10−5 ~4.4×10−6
Optical homogeneity 1×10−5 1×10−5 1×10−6 1×10−4
Fracture strength/MPa 90 120 60 ~18
Stability High High Medium Medium
Relative costs Medium Low Low Medium

红外锗单晶制备通常采用直拉法[3](CZ)、垂直梯度凝固法[4](VGF)和水平区熔法等。垂直梯度凝固法、水平区熔法一般只能用于6英寸[5](152.4 mm)以下的锗单晶[6]。直拉法是制备大直径红外锗单晶较为可行的技术途径。目前全球技术领先的红外锗单晶企业主要有比利时的优美科公司(Umicore)、美国的Ⅱ-Ⅵ公司和德国的Photonic Sense等。比利时优美科公司(Umicore)公开报道的红外锗材料最大直径达到Φ600 mm,光学均匀性优于5×10−5。美国Ⅱ-Ⅵ公司能够生产直径Φ450 mm以上的红外锗材料,光学均匀性优于2×10−5。德国Photonic Sense是生产高品质锗红外光学镜片的世界级厂商,拥有世界上最大的晶体生长炉,依靠先进的数控设备,实现了各种复杂锗产品的生产,目前已知生产的最大红外锗材料的直径可达Ф520 mm,纯度5N(99.999%)以上,光学均匀性优于5×10−5[7]。目前,公开资料表明国内从事大直径红外锗材料研制和生产的单位均具备Φ350 mm左右红外级锗单晶的制备能力,但距离国际先进水平仍有较大差距。Φ350 mm左右的红外锗单晶制备工艺已比较稳定[8],而Φ350 mm以上的锗单晶制备成功率不高[9]。表2为大直径高均匀性锗单晶研制技术难点分析。

表2 大直径高均匀性锗单晶研制技术难点
Table 2 Different main heater and bottom heater power control group
Technical content Technical difficulty
Growth technology of super large diameter germanium single crystal High temperature convection intensity,heat field instability;large amount of feed leads to complex thermal convection and crystal defects;easy to form twin defects;it is difficult to ensure a smooth solid-liquid interface
Optical uniformity control technology Control of impurity distribution and uniformity of doping;high dislocation density
Residual stress relief technology Thermal stress in growth process is high;residual stress of ingot is large,and it is easy to crack in processing
Material testing technology High precision refractive index test at low temperature;large size uniformity test

近年来,随着空间红外光学遥感载荷的快速发展,为了获取更大的观测视场和更高的灵敏度及分辨率[10],对锗单晶的直径和电阻率均匀性[11]均提出了更高的要求。其中,锗单晶直径要求不小于Φ420 mm,全口径均匀性不低于1×10−4。目前国内受工艺装备及工艺限制,大直径锗单晶制备成功率不高,其原因主要在于锗单晶直径越大[12],对热场温度控制要求苛刻,成晶难度大,并且固液界面形状难以控制,从而导致电阻率分布不均匀,影响光学均匀性。

因此,为了满足空间红外光学遥感载荷需要的大直径、高均匀性锗单晶[13]需求,本文采用直拉法[14],在研制首台国产大尺寸单晶炉基础上,结合CGSim模拟软件,系统研究分析了新型双加热器热场[15]条件下大直径红外锗单晶成晶[16]规律;并对双加热器热场条件下Φ420 mm锗单晶制备工艺进行了验证与分析。

1 热场设计

在大尺寸单晶炉中设计了40英寸(Φ1000 mm)热场,采用底加热器为主加热器的双加热器组合式热场结构[17],如图2所示,加热器材质为三高石墨,主加热器结构为筒状结构围绕在坩埚四周,底加热器形状为圆盘结构,位于坩埚底部。加热控制的电源系统可以实现主加热器和底加热器功率的单独调节,主加热器和底加热器电源的功率调节范围分别为0~150 kW和0~100 kW。同时,优化设计高效换热的水冷屏,设置水冷屏在大直径单晶炉所处的合理位置,水冷屏与加热器形成合理单晶提拉所需的热场。

图2 40英寸双加热器组合式热场示意图
Fig.2 40-inch dual-heater combined heat field diagram

锗单晶制备过程中,加热器使坩埚内的区熔锗处于熔液状态,坩锅以合适的速度转动保证锗熔液温场及参杂的均匀性,同时籽晶以合适的速度向上提升,在水冷屏冷却作用下,快速实现锗单晶的成晶。

2 实验

2.1 样品制备和测试

在同一台大尺寸单晶炉中,使用同一套双加热器热场系统,选用4组不同底加热器功率作为变量,调节主加热器功率实现晶体正常拉晶,其中底加热器功率为0 kW时,可近似为单加热器系统作为对比。投入实验原料为同一批次的6N(纯度99.9999%)区熔锗锭,投料量同为350 kg,并且保证投等量的锑作为掺杂剂,晶转和埚转同为6∶2,坩埚位置和放肩角度相同,等径时提拉速度同为9 mm·h−1,等径长度为50 mm。当单晶达到预设长度后,手动快速提升晶升将单晶提断,待炉内冷却后,出炉单晶并观察单晶界面形状,测量晶体界面的突出度,以及截面电阻率[18]分布情况,统计单晶成晶率见表3。

表3 不同主加热器和底加热器功率对照组
Table 3 Different main heater and bottom heater power control group
Bottom heater power/kW Feed quantity/kg Single crystal diameter/mm Length of equal diameter/mm Main heater power/kW
0 350 420 50 105
10 350 420 50 94
20 350 420 50 82
30 350 420 50 69

2.2 双加热器热场模拟

实验使用CGsim软件对双加热器热场系统中Φ420 mm锗单晶的生长情况进行模拟。CGsim软件[19]是俄罗斯的一款专业晶体生长模拟器,该模拟工具可以通过建模计算,分析热场中的温度梯度分布、晶体生长过程中的传热、界面形状以及单晶应力分布情况[20]等。本实验使用此软件将4种不同底加热器功率下生长条件进行导入,并通过模拟计算,分析得出Φ420 mm锗单晶生长过程中相应的温度梯度分布情况。

3 结果与讨论

3.1 数值模拟结果

图3为使用CGsim软件在不同底加热功率,相同拉速等工艺参数下,拉制Φ420 mm锗单晶,等径50 mm长时的温度分布和界面形状模拟情况。可以看出,底加热器功率为0 kW时,单晶界面凸出度最大,结晶前沿中心位置和边缘位置的温度差最大。随着底加热器功率的增加,界面趋向于平坦,当底加热器功率为20 kW时,晶体结晶前沿界面趋近平坦,接近较为理想的固液界面形状[21],此时结晶前沿中心位置和边缘位置的温度差为0.4 K,径向温度梯度较小;底加热器功率升至30 kW时,结晶前沿中心位置和边缘位置的温度差更小为0.3 K,但是晶体结晶前沿界面形状成为凹凸有序的双曲面,这不利于单晶成晶和电阻率均匀性控制。因此,通过模拟可分析出,随着底加热器功率的增加,熔体中心温度会逐渐升高[22],熔体径向温度梯度较小,然而当底加热器功率过大时,会使得熔体径向温度梯度过小,造成固液界面形状向不利的方向发展。不同底加热功率条件下熔体径向温度梯度变化情况见图4。可知,当底加热功率为0~20 kW时,熔体径向温度呈现从坩埚中心向边缘逐渐升高的趋势,此时主加热器起主要作用,而当底加热器功率达到30 kW时,熔体中径向温度呈现从中心向边缘逐渐降低的趋势,此时底加热器起主要加热作用。

图3 不同底加热功率条件下拉制Φ420 mm锗单晶,等径50 mm长时的温度分布和界面形状模拟情况
Fig.3 Temperature distribution and interface shape simulation of drawing equal diameter 50 mm length Φ420 mm germanium single crystals under different bottom heating power conditions
(a)0 kW;(b)10 kW;(c)20 kW;(d)30 kW
图4 不同底加热功率条件下熔体径向温度梯度变化情况
Fig.4 Change of radial temperature gradient of melt under different heating powers

3.2 双加热器热场对单晶成晶影响分析

在4种底加热器功率条件下,分别开展3次Φ420 mm锗单晶拉制,其他参数保持相同,进行定角度放肩以及定拉速等径。统计Φ420 mm锗单晶的成晶次数以及拉制过程中发生晶变的最早直径位置的平均值,如表4所示,从而得到双加热器系统下,底加热器功率的调节对单晶成晶的影响情况。可以看出,只有在底加热器功率为20 kW时,成功拉制出未发生晶变的Φ420 mm锗单晶,并且随着底加热器功率的增加,单晶发生晶变的最早直径呈现出先升后降的趋势,如图5所示,底加热器功率过大或者过小都会直接影响Φ420 mm锗单晶的成晶情况。

表4 不同底加热器功率条件下单晶成晶情况
Table 4 Single crystal crystallization under different bottom heater power conditions
Bottom heater power/kW Average diameter of the earliest crystal transformation/mm Success times of germanium single crystal Crystallization rate of germanium single crystal
0 354 0 0
10 384 0 0
20 402 2 66%
30 332 0 0
图5 单晶最早晶变直径与底加热器功率变化关系
Fig.5 Relationship between diameter of the earliest crystal change and power of bottom heater

在相同热场条件下,放肩直径变化越趋近于线性,单晶发生晶变的概率越低。因为当单晶放肩过程中,直径增长速度稳定时,表明热场中的环境越趋于稳定,晶体生长过程中产生的结晶潜热能够顺利释放,并且结晶前沿过冷度不会出现较大变化。这就要求单晶放肩和等径过程中的拉速变化需要与热场相适配。在底加热器20 kW功率条件下,对晶体放肩和等径过程中,拉速与直径的对应关系进行了记录分析,发现单晶生长过程中,拉速范围为0.1~0.25 mm·min−1,单晶生长形态较好,放肩较为均匀,两者具体对应关系见表5,晶体拉速-直径变化曲线如图6所示。

表5 单晶拉速与直径对应关系
Table 5 Corresponding relationship between single crystal pulling speed and diameter
Points 1 2 3 4 5 6 7 8 9 10
Diameter/mm 80 120 160 210 250 280 300 340 380 420
Pulling speed/(mm·min−1) 0.23 0.22 0.21 0.19 0.17 0.16 0.14 0.11 0.1 0.13
图6 单晶拉速-直径变化曲线
Fig.6 Single crystal drawing speed diameter variation curve

3.3 电阻率均匀性分析

单晶的电阻率均匀性,与单晶直拉过程中的固液界面形状相关。对4种底加热器功率条件下,生长的4种不同界面的Φ420 mm锗晶体进行了截面切割,并且依据《GBT 26074—2010锗单晶电阻率直流四探针测量方法标准》中电阻率测量取点方法进行取点测量,为减少误差,取点过程中选取该位置附近的单晶晶粒进行测量,依据测试结果,计算单晶径向电阻率不均匀性见表6。可以看出,4种底加热器功率条件下的拉晶界面情况于模拟结果基本对应,其中底加热器功率为20 kW时,晶体的界面突出度最小,同时电阻率不均匀性也最低。

表6 电阻率均匀性和单晶界面突出度统计
Table 6 Resistivity uniformity and single crystal interface protrusion statistics
Bottom heater power/kW The lowest value of resistivity/(Ω·cm) Highest resistivity/(Ω·cm) Heterogeneity of resistivity/% Interface protrusion of single crystal/mm
0 23 36 56.5 42
10 22 32 45.5 28
20 26 31 19 13
30 24 30 25.0 Hyperboloids

3.4 光学均匀性测试分析

将底加热器功率为20 kW制备出的Φ420 mm锗单晶进行光学均匀性。单晶(未镀膜状态)按《GJB 2019A—2015 红外光学用锗规范》进行了光学均匀性测试,具体数据见表7,通过测试结果可以看出,锗单晶的光学均匀性较好。

表7 光学均匀性测试结果
Table 7 Optical uniformity test results
Thickness/mm Diameter/mm Uniformity Note
45 Φ420 6.52×10−5 Average

4 结论

通过设计并使用双加热器热场系统,在底加热器功率20 kW,投料量350 kg条件下,成功制备出径向电阻率均匀性≤19%的Φ420 mm锗单晶,并且通过使用数值模拟对双加热器系统中热场温度分布和底加热器功率对晶体固液界面形状的影响进行了分析,得出以下结论:

1)底加热器对熔体径向温度梯度分布能够产生较大影响,在合适范围内,底加热功率越高,径向温度梯度越小,越有利于Φ420 mm锗单晶的成晶。

2)底加热器功率影响熔体径向温度会间接影响到晶体的固液界面形状,在一定范围内,随着底加热器功率的增加,固液界面突出度减小。

3)大直径锗单晶的界面突出度直接影响晶体的径向电阻率分布情况,通过双加热器结构可以调节晶体的固液界面形状,从而实现调节径向电阻率。

4)该工艺方法制备的Φ420 mm锗单晶,光学均匀性达到6.52×10−5,电阻率均匀性达到19%,均优于使用指标要求。

通过Φ420 mm大直径红外锗单晶的研制技术研究,突破了拉晶所需高精度双温场控制技术瓶颈,有望打破国外锗单晶制备技术的垄断、加速国产化材料制备和应用,为空间工程化应用奠定了基础,同时为同类型的大直径单晶材料制备提供参考借鉴。

参考文献:

[1] 

(左建坤, 潘美妍, 段辉高, 贾红辉, 胡跃强. 新型红外隐身结构材料研究综述 [J].光电工程, 2023, 50(5): 7.)

[2] 

(黄宇飞, 徐嘉, 白绍竣, 高冀, 李昂, 高洪涛. 卫星红外遥感术及其在防灾救灾中的应用研究 [J]. 航天返回与遥感, 2020, 41(5): 118.)

[3] 

(董汝昆, 吴绍华, 王柯, 尹国良, 史娜娜, 姚杨, 郭晨宇. 锗单晶材料的发展现状 [J].红外技术, 2021, 43(5): 510.)

[4] 

(马奎, 洛桑才仁, 陈超, 张万益, 王丰翔. 全球锗资源分布、供需及消费趋势研究 [J].矿产保护与利用, 2019, 39(5): 16.)

[5] 

(廊坊传媒网. 有研国晶辉: 争做红外材料产业“领头羊”[EB/OL]. 2021-11-19 [2023]. http://www.lfcmw.com/shms/content/2021-11/19/content_921068.html.)

[6] 

(李宝学, 赵燕, 罗玉萍, 孙吉单. 大直径锗单晶制备工艺研究及其指标分析 [J].云南冶金, 2016, 45(6): 50.)

[7] 

(陈安萌, 梁德印, 聂云松, 李昊谦, 韩波. “资源一号”02E星红外相机探测器成像验证方法研究 [J].航天返回与遥感, 2021, 42(4): 37.)

[8] 

(葛桓宇, 肖正航, 王跃. 基于热变形补偿的双材料低温镜头支撑结构研究 [J].航天返回与遥感, 2022, 43(3) : 69.)

[9] 

(马绍芳, 裴玉兰, 冯德伸, 解红霞. 大直径红外光学锗单晶 [J].广东微量元素, 1998, 5(2): 64.)

[10] 

(潘朝猛, 康丽珠, 罗敏, 陶亮, 陈树刚, 陈波, 白忠宏, 崔海, 徐参军, 赵劲松. 空间红外相机技术发展现状及应用 [J].红外技术, 2022, 44(11): 1186.)

[11] 

(张红勇, 姚婕. 锗单晶直拉生长法的自动等径控制方法[P]. 中国, CN 102758250A, 2012.)

[12] 

(高欢欢, 黎建明, 冯德伸, 王霈文, 李葳. <100>p型4英寸无位错锗单晶的研制 [J].稀有金属, 2018, 42(3) : 285.)

[13] 

(丁国强, 苏小平, 张峰燚, 黎建明, 冯德伸. VGF法生长6英寸Ge单晶中热应力的数值模拟优化 [J].稀有金属, 2011, 35(2): 244.)

[14] 

(肖祥江, 惠峰, 董汝昆, 吕春富. 用VGF法生长6英寸锗单晶中籽晶熔接工艺研究 [J].半导体技术, 2017, 42(11): 860.)

[15] 

(陈燕. 直拉法生长锗单晶的水平放肩工艺实践 [J]. 云南冶金, 2002, 31(6): 36.)

[16] 

(原洛渭, 高利强. 双段式加热锗单晶生长设备的结构设计 [J].电子工业专用设备, 2008, 157(2): 11.)

[17] 

(冯德伸, 苏小平, 闵振东, 尹士平. Φ300 mm红外锗单晶生长及性能测试 [J].稀有金属, 2006, 30(5): 694.)

[18] 

(陈晨, 赵堃, 韩焕鹏. 6英寸低位错锗单晶生长热场设计 [J].人工晶体学报, 2021, 50(6): 979.)

[19] 

(闫洪. Φ6.5 mm×8.9 mm半导体锗单晶圆片电阻率的测定 [J].电子质量, 2007, (10): 50.)

[20] 

(单志远, 秦瑞锋, 宁永铎, 李洋, 张果虎. 水平超导磁场直拉法硅单晶熔体过热现象 [J].稀有金属, 2023, 47(7): 1050.)

[21] 

(左建龙, 冯德伸, 李楠. 锗单晶中位错密度的影响因素 [J].稀有金属, 2010, 34(5): 726.)

[22] 

(穆宏赫, 王鹏飞, 施宇峰, 张中晗, 武安华, 苏良碧. 热交换坩埚下降法制备大直径氟化铈晶体的热场设计与优化 [J].无机材料学报, 2023, 38(3): 288.)

顶部