谭小波,闫欣,易涛,何凯,邵铮铮,周凯凯,高贵龙,汪韬,张军,庄钊文
(1 国防科技大学电子科学学院,长沙410073)
(2 中国科学院西安光学精密机械研究所超快诊断技术重点实验室,西安710119)
(3 中国工程物理研究院激光聚变研究中心,四川绵阳621900)
惯性约束聚变(Inertial Confinement Fusion,ICF)是将高能短脉冲激光注入黑腔对靶丸加热压缩以实现聚变点火的方法,超快诊断技术在ICF 研究中不可或缺[1-2]。ICF 发生在极小的空间(0.1~1 mm),极短的时间(~ns),且部分过程如内爆热斑的持续时间只有100~200 ps,这需要诊断技术具备ps 级甚至百fs 级的超高时间分辨能力以及μm 级的高空间分辨能力[3]。传统电真空超快探测技术通常基于光电转换原理,其性能受到材料响应、空间电荷效应等因素的制约,比如,采用微通道板(Microchannel Plate,MCP)技术可以将分幅相机的时间分辨率提升到~30 ps,但MCP 带来的空间电荷效应限制了其空间分辨率(~50 μm)[4-6];DIXI(DIlation X-ray Imager)利用电子脉冲时间放大技术以获得更短的选通时间,实现了较高时间分辨的成像(~10 ps),但DIXI 中会产生大量次生电子,影响成像的信噪比和空间分辨率(~100 μm)[7-8];条纹相机虽能够达到200 fs 的超高时间分辨率,但其本身不具备二维成像能力,较大地限制了其应用[9-10]。
美国劳伦斯利弗莫尔国家实验室提出了基于半导体光折变效应的超快探测技术,时间响应可达ps~百fs,还具有全光全固态抗强辐射能力等优势[11-12]。LONDON R A 等对GaAs 材料的能量沉积时空特性进行了理论分析和数值计算,研究表明GaAs 的时间响应低至百fs 级,空间分辨优于5 μm,光折变材料可用于高时空分辨超快成像研究[13]。BAKER K L 等采用硒化镉(CdSe)作为响应材料实现了时间分辨率为5 ps 的两分幅X 射线成像,而CdSe 非平衡载流子复合时间较长,需要采用差分方式获得较高时间分辨的图像,其成像的空间分辨率达到30 lp/mm,画幅约为2.5 mm×1.5 mm[12]。
国内,中国工程物理研究院对CdSe 的超快动力学开展了研究,得到其载流子复合时间为470±20 ps[14]。西北核技术研究所将光折变材料应用于Z 箍缩,采用磷化铟(InP)实现了时间分辨率约为1.5 ns、空间分辨率约为140 μm 的成像,并分析了其在5 eV~100 keV 范围内的光谱响应[15]。中科院西安光机所等单位利用砷化镓/铝镓砷(GaAs/AlGaAs)多量子阱结构半导体材料得到了3 ps 高时间分辨的六分幅可见光成像,该材料的非平衡载流子寿命低至2.5 ps,响应时间远小于CdSe,无需进行差分即可实现高时间分辨成像[16-17]。王博等以低温GaAs 为响应材料,实现了l0 ps 时间分辨X 射线探测,但未进行二维成像[18]。钟梓源等对低温生长铝镓砷(Low-Temperature grown AlGaAs,LT-AlGaAs)的光折变效应进行了理论和实验研究,得到LT-AlGaAs 的非平衡载流子的复合时间低至2.08 ps,适用于超快响应芯片的研制[19],但未开展X 射线成像研究。目前为止,尚未出现报道基于光折变超快响应芯片的大画幅X 射线成像实验研究。
因此,本文采用LT-AlGaAs 为材料研制的超快响应芯片进行X 射线成像实验,利用高能量纳秒激光打靶产生X 射线,以此为信号源对LT-AlGaAs 光折变X射线响应芯片的空间性能开展实验研究。实验结果验证了利用LT-AlGaAs 超快芯片进行高空间分辨、大动态范围、大画幅X 射线成像的可行性,为光折变超快X 射线成像系统的研制及应用提供了实验和数据支撑。
半导体的光折变效应是光折变X射线半导体响应芯片的研制基础,当X 射线信号的光子能量大于半导体带隙时,材料将吸收能量并在内部产生非平衡载流子。由于半导体的带收缩、带填充以及自由载流子吸收等效应,这些载流子会在极短的时间尺度(fs~ps)减速并导致材料的介电常数和折射率发生变化[13,18]。折射率变化量Δn与信号强度I在一定条件下成近似线性关系Δn=αI,其中α为常系数[13]。当探针光同时入射到材料内时,其相位将会随着材料折射率的变化而改变[13],相位改变量为
式中,d为半导体材料的厚度,λ为探针光波长。因此,通过材料的光折变效应,信号光的信息被调制成探针光的相位变化量,其大小和信号光的强度成线性关系,通过进一步的相位提取可恢复出目标信号。
本文采用相位光栅编码进行衍射的方法实现相位信息的提取、设计,研制了对应的光折变响应芯片,其主要由光折变材料LT-AlGaAs、光学镀膜(探针光高反膜),金属掩膜(一维光栅,70 lp/mm)组成。其工作原理如图1所示,当X 射线信号入射到芯片正面时,首先被金属光栅调制成一维周期性结构,透过光栅的信号光将到达光折变材料内部,导致材料折射率发生周期性变化,形成一维瞬态相位光栅。另一侧,探针光同步入射到芯片背面,经高反膜反射,探针光的相位将会随材料折射率的改变而发生变化,并被产生的一维相位光栅衍射。
图1 光折变响应芯片的工作原理示意图Fig.1 Principle diagram of the photorefractive chip
探针光经一维相位光栅衍射后,再通过光学4f滤波系统进行相位信息提取,如图2所示,在频谱面上放置滤波器,设置小孔只让±1 级光通过,最终到达像面进行成像。根据相位光栅的衍射理论[20],各个级次的衍射效率
图2 相位光栅编码提取相位原理图Fig.2 Schematic diagram of phase extraction by phase grating coding
因此,滤波让+1 或-1 级通过,最后到达探测器的光强度
式中,Iprobe表示探针光的强度。由此,将相位变化量转为探针光光强的变化,实现了相位信息的提取。
测量实验利用高能量激光脉冲轰击金属靶材产生的X 射线,以此作为信号来激发超快芯片的响应。装置布局如图3所示,高能量纳秒激光器发射出打靶激光脉冲(脉冲宽度2 ns,能量从5~155 J 分段调节),通过抛物面镜(Off-Axis Parabolic Mirror,OAP)汇聚到靶点(Al靶),产生X 射线。X 射线从靶点入射到芯片的正面,通过金属掩膜调制后诱发芯片内半导体折射率的变化。另一侧,由探针光激光器发射出探针光(785 nm),途经反射镜和分光镜后入射到芯片的背面,而后经过时变的折射率调制再反射进入4f系统(由两个焦距为30 cm,口径为5 cm 的透镜组成),滤波器只允许+1 级通过,最后被增强型电荷耦合器件(Intensified CCD,ICCD)接收。实验中,配有一套X 射线二极管探测器(X-Ray Diode,XRD)对X 射线的能量进行实时在线监测。
图3 实验装置布局图Fig.3 Layout of the experimental setups
在芯片正面添加多个“8”字图案的金属掩膜对X 射线信号进行调制,以此作为实验的目标场景,如图4所示。ICCD 的参数设置为增益80,门宽150 ns。实验中通过调节打靶激光脉冲的能量,分析不同X 射线总能量条件下成像系统的空间分辨能力。X 射线成像结果如图5所示,分别为打靶激光能量为5 J、23 J、46 J 和155 J 的实验结果,对应的XRD 测量值换算为X 射线能量并归一化后分别为0.008 3、0.17、0.57 和1。成像的单幅画幅尺寸达到6.7 mm×6.7 mm,在实验所采用的激光能量下,均能较好成像,能够清晰展示目标的细节,随着能量的增强,成像整体的灰度值也随之增加。图像右半部分有数个亮点,这是由芯片加工出现瑕疵带来的强散射光,可以通过进一步的工艺迭代优化进行改善。
图4 芯片表面金属掩膜示意图和ICCD 直接拍摄掩膜的结果Fig.4 Schematic diagram of the mask on chip surface and imaging result of the mask using ICCD
进一步对图5 进行分析,利用数字“8”左下角(图5 右下图中的红实线)的水平强度变化曲线作为边缘扩散函数(Edge Spread Function,ESF),以此来计算成像系统的调制传递函数为了减少误差,对图像进行适当的旋转后,再对其沿竖直方向合理范围内进行平均,计算得到的MTF 如图6所示。
图5 不同打靶激光能量下的成像结果Fig.5 Results under different pump laser energy
图6 不同X 射线能量(归一化)下成像调制传递函数计算结果Fig.6 Calculation results of MTF under different normalized X-ray energy
根据成像质量评价理论,通常将MTF 为0.1 时的空间频率作为系统的最大空间分辨率,由图6 可得该成像系统的最优空间分辨率可达35 lp/mm,优于LLNL 公开报道的结果[12]。受限于实验条件,实验中杂散光控制、系统的稳定性等因素未能进一步优化,成像系统的空间分辨能力在不同X 射线能量下有一定的波动,但均能达到25 lp/mm @MTF=0.1 的空间分辨率。
对不同能量下的成像选定相同的ROI 区域(图5 右下图中红虚框内的区域),对区域内的所有像素的灰度求平均,得到像素平均灰度值与X 射线归一化能量的关系,如图7所示。对实验数据进行线性拟合(图7中红线),拟合的R-square 值(确定系数)为0.97,说明像素平均灰度值和X 射线归一化能量基本满足线性关系。针对图像传感器,动态范围一般定义为信号饱和值与绝对灵敏度阈值之比[21]。受实验打靶发次限制,未能采集到信号饱和值。但在最强X 射线能量下能够清晰成像,也就是低于信号饱和值。另外,绝对灵敏度阈值一般是图像中信噪比为1 时的信号值[21]。由实验结果图5 可知,最弱X 射线能量下的信噪比也是大于1 的。因此,本实验中的X 射线响应芯片的线性动态范围应当优于。
图7 ROI 区域平均灰度值与X 射线归一化能量的关系Fig.7 Relationship between average gray value of ROI and normalized X-ray energy
分析系统的MTF 要考虑以下四方面因素的影响:1)半导体内X 射线能量沉积带来的空间弥散足够小[13],允许成像系统优于5 μm 的分辨率;2)金属光栅的调制会影响成像的空间分辨率,本文采用的芯片表面刻蚀的光栅为70 lp/mm;3)探测器像元的尺寸为12.8 μm,本成像系统的放大率为1,则探测器能表示的最大空间频率约为39 lp/mm;4)光学4f系统的衍射受限效应也会影响空间分辨,当频谱面不加滤波器时,可由透镜焦距及口径计算出光学系统可通过的最大空间频率约为106 lp/mm[22],加滤波器后,由于通光孔的直径需小于衍射零级和一级的间距,其截止作用将大大较低光学系统的空间分辨率。因此,当前系统的MTF 主要受到滤波器上小孔尺寸和探测器像元尺寸的限制,通过增大调制光栅的线对数,使得频谱面零级和±1 级的间距增大,这样可以在屏蔽掉零级光的同时可一定程度增大小孔的尺寸,同时减小探测器像元尺寸或增大成像放大率,系统的空间分辨率会得到进一步提升。
当前系统的MTF 结果表示响应芯片及后端相位提取系统的组合分辨能力,在实际应用中,目标场景要先通过一定放大倍率的成像系统(如针孔成像系统、KB 显微镜等)后成像到芯片上,这样可以进一步提升整系统的物方空间分辨[3]。对于观测相同的目标区域,大画幅可搭配更大的放大率,从而获得更高的物方空间分辨。当前系统中的超快响应芯片实现了6.7 mm×6.7 mm 的大画幅成像,相比于LLNL 研究中约为2.5 mm×1.5 mm 的成像画幅[12],本系统具备更佳的空间分辨能力。
本文实验验证了在120∶1 的X 射线能量动态范围内的大画幅高空间分辨成像的能力,基于LT-AlGaAs材料本身具备的2.08 ps 的超短非平衡载流子复合时间,可综合得出LT-AlGaAs 超快响应芯片可用于大动态范围(≥120∶1)、高时间分辨(~ps)、高空间分辨(~20 μm)的大画幅(≥6.7 mm×6.7 mm)成像。结合先进的超快分幅技术,例如色散分幅[23]、偏振啁啾分幅[16]等技术,可基于LT-AlGaAs 超快响应芯片研制具有高时空分辨的大画幅成像系统,这是下一步研究的重点方向。
基于光折变效应的超快成像技术是一种新型诊断技术,具有高时空分辨、全光全固态抗辐射等优点。本文对LT-AlGaAs 超快响应芯片的空间性能进行研究,验证了成像系统在X 射线谱段大画幅(6.7 mm×6.7 mm)高空间分辨的成像能力,在120∶1 的X 射线入射能量动态范围内,X 射线成像的空间分辨率的最优空间分辨≥35 lp/mm @MTF=0.1。本文的研究结果进一步验证了基于光折变材料的超快X 射线诊断技术的可行性,未来将通过与超快分幅技术相结合,实现多分幅高时空分辨的二维X 射线成像,着力推进ICF 超快诊断的研究进展。