光学学报, 2018, 38 (9): 0934001, 网络出版: 2019-05-09   

基于Angel型龙虾眼X射线透镜的数值模拟 下载: 1065次

Numerical Simulation of Angel Lobster Eye X-ray Micro Pore Optic Lens
作者单位
1 北方夜视技术股份有限公司南京分公司, 江苏 南京 211106
2 中国科学院国家天文台, 北京 100012
摘要
为了研究Angel型龙虾眼X射线透镜的聚焦成像特性,基于X射线全反射原理和旋转坐标系方法,建立X射线在方孔内壁的数值模型。通过求解X射线在方孔内壁的交点,得到所有光线的传输路径。为了验证模型的准确性,对透镜的焦距和传输效率分别进行实验测试和模拟。在365 mm焦距处,参与反射的X射线被会聚为十字线,中心焦斑光强最大,与模拟所得结果相符。在能量4.5 keV下,透镜镀金属Ir膜前、后的传输效率分别为1.23%和9.18%,模拟结果的传输效率分别为1.44%和10.14%。结果表明:构建的数值模型是合理的,可为龙虾眼透镜的研制提供理论基础。
Abstract
In order to study the focusing and imaging performance of Angel lobster eye X-ray micro pore optic lens, the numerical model of the X-ray in the square channel is established based on X-ray reflection principle and rotating coordinate system. By solving the intersection, the transmission path of all the rays is obtained. In order to verify the accuracy of the model, the focal length and transmission efficiency of the lens are both simulated and tested. At the focal length of 365 mm, the reflected X-rays are focused into cross line, and the center focal spot intensity is the largest, which is consistent with the simulation results. As for the transmission efficiency of the lens before and after coating Ir, the simulation results are 1.23% and 9.18%, respectively, at the energy of 4.5 keV. Accordingly, the experimental results are 1.44% and 10.14%, respectively, which is consistent with the above data. The results show that the numerical model is reasonable and provides the theoretical and simulation basis for the fabrication of the lens.

1 引言

1979年,Angel根据龙虾眼结构,提出基于正方形阵列结构来研制龙虾眼型(LE)X射线天文望远镜[1]。LE系统具有视野大、焦距短、体积小、质量轻、灵敏度高、装调简单、聚焦效率高等优点,符合未来卫星载荷X射线天文观测发展的需求,在X射线聚焦成像领域拥有巨大的应用前景[2-7]。受限于当时的制作工艺和加工水平,直到20世纪90年代,随着一种新型微孔光学技术(MPO)的日渐成熟,Angel型LE系统才得以实现和推广[8-9]。MPO透镜由数百万个指向球心的方孔通道构成,通过模拟X射线在方孔中的传输规律来改进MPO的制造工艺,正成为LE光学系统的一个重要研究方向[10]

关于X射线在龙虾眼微孔内壁的模拟有几种典型的方法。Chapman等[11]提出了几何光学数值模型,但未考虑透射光线的影响。Tich y'[12]和Peele等[13]提出了迭代数值模型,但忽略了X射线反射率的影响。Putkunz等[14]和Zhao等[15]基于光线追踪软件进行成像模拟,但构建模型复杂、耗时太长。本文根据X射线全反射原理和旋转坐标系方法[16],建立X射线在方孔通道传输的数值模型,通过直接求解方程组得到所有光线在方孔内部的传输路径,大大提高了运算效率,缩短了模拟时间,并通过实验验证了构建数值模型的准确性。

2 MPO透镜数值模型

2.1 X射线全反射

X射线是短波长、高能量的电磁波,当以一定角度入射至光滑表面时存在全反射临界角θc。当X射线入射角度小于θc时,会在物质表面全反射。临界角θc通常与入射X射线能量和材料相关[17],可表示为

θc=20ρ/E,(1)

式中:ρ为材料的密度(g/cm3);E为入射X射线的能量(keV);θc为掠入射临界角(mrad)。可以看出,X射线临界角与入射光子能量成反比,与材料密度成正比。为了提高掠入射临界角和X射线反射率,通常在X射线微孔内壁镀上金属反射膜[18]

在实际情况中,方孔内壁表面也存在一定的粗糙度,会引起X射线的漫散射,导致反射率以指数形式下降。当掠入射角θ<θc时,Nevot-Croce(NC)模型的计算过程更为准确[19],在理想反射率公式右侧乘以一个Nevot-Croce因子,X射线在表面有粗糙度的反射率R可表示为

R=R0exp-(4πσsinθ)2n1n2λ2,(2)

式中:R0为理想光滑表面的反射率,n1为X射线在空气中的折射率;n2为X射线在材料中的折射率;σ为表面粗糙度均方根参数;θ为X射线掠入射角;λ为X射线波长。图1给出了在粗糙度σ=1 nm情况下,不同能量下X射线入射角与相对反射率的关系,可以看出,X射线反射率与光线能量和入射角成反比。

图 1. 不同能量下X射线入射角与相对反射率的关系

Fig. 1. Relationship between X-ray incident angle and relative reflectivity under different energies

下载图片 查看所有图片

2.2 数值模型

MPO透镜是一种基于掠入射原理对X射线进行会聚的光学器件,开口面积Aopen约60%,方孔分布为轴对称性结构,如图2所示。对第一象限内的方孔结构[图2(b)]进行建模,再将其他象限的方孔旋转至第一象限进行计算,即可得到MPO镜片所有方孔的分布情况。由于每个方孔通道中心都指向球心,其最外层曲线方程可构建为

f(z)=k0z+b0,(3)

则在第i行、第j列个方孔位置处的上壁方程为

ku(i,j)z+bu(i,j)=y,(4)

式中:ku(i,j)=j+an+0.5k0;bu(i,j)=j+an+0.5b0;系数a与通道的开口面积有关,可表征为a= Aopen2

图 2. (a) MPO透镜截面图;(b)第一象限方形通道

Fig. 2. (a) Cross-section of MPO lens; (b) area of interest of square channel

下载图片 查看所有图片

同理,可以将方形微通道下侧、左侧和右侧内壁方程分别简化为

kdi,jz+bdi,j=y,5

kl(i,j)z+bl(i,j)=x,(6)

kr(i,j)z+br(i,j)=x,(7)

式中:kd(i,j)= j-an-0.5k0;bd(i,j)= j-an-0.5b0;kl(i,j)= i-an-0.5k0;bl(i,j)= i-an-0.5b0;kr(i,j)= i+an+0.5k0;br(i,j)= i+an+0.5b0

假设光源的位置为ps(xs,ys,zs),随机生成透镜(i,j)的入口坐标为pin(xin,yin,zin),则入射光子的方向向量u

u=(ux,uy,uz)=(xin-xsd0,yin-ysd0,zin-zsd0),(8)

式中d0= (xin-xs)2+(yin-ys)2+(zin-zs)2

入射光子所在的直线方程可以表示为

xin-xsux=yin-ysuy=zin-zsuz(9)

通过联立方程组,可以得到X射线与方孔上壁的交点zu

zu=ys-uxuzzs-bu(i,j)ku(i,j)-uyuz(10)

同理,可以求得方形微通道下侧、左侧和右侧交点坐标分别为

zd=ys-uxuzzs-bdi,jkdi,j-uyuz,11

zl=xs-uyuzzs-bl(i,j)kl(i,j)-uxuz,(12)

zr=xs-uyuzzs-br(i,j)kr(i,j)-uxuz(13)

在上述z坐标值中z=zmin,认为对应的z坐标值为入射X射线光子与方形微通道内壁的实际交点,则入射法向量n方程为

n=(nx,ny,nz)=mxmx2+my2+mz2,mymx2+my2+mz2,mzmx2+my2+mz2(14)

对于方形微通道上侧内壁,外法向量在X轴上的分量为0、Y轴上的分量为正。因此,任意z位置的mux,muy,muz可以分别表示为

mux=0;muy=1;muz=ku(i,j)(15)

同理,方形微通道下侧、左侧和右侧法线方程分别为

mdx=0;mdy=-1;mdz=kdi,j,16

mlx=-1;mly=0;mlz=kl(i,j),(17)

mrx=1;mry=0;mrz=kr(i,j)(18)

若满足全反射条件θθc,则入射X射线会发生全反射,且作为下一次反射的入射X射线。反射X射线的方向向量r可以表示为

r=(ux-2nx×u×n,uy-2ny×u×n,uz-2nz×u×n)(19)

r作为入射光线,迭代运算,可得到光线在MPO通道内部的传输路径,实现光线追踪。假设接探测器的焦平面距离MPO镜片的距离为f,出射X射线光子最后一次全反射与方形微通道内壁的交点为pf(xf,yf,zf),则在焦平面上的交点可以表示为(xf,yf):

xf=xs+rxrzf-zf,20

yf=ys+ryrz(f-zf)(21)

MPO透镜具有中心对称的结构,因此可将其余3个象限内的方形微通道均旋转至第一象限内进行计算。由于入射X射线光源的空间位置是随意设置的,因此在旋转其他象限内方形微通的同时,还须将X射线光源旋转相同的角度。以第二象限内的一点p0(x0,y0)为例,经过坐标旋转之后的位置p'(x',y')可表示为

x'=qcosφ,22

y'=qsinφ,(23)

式中q= x02+y02,φ=tan-1( y0x0)- π2

假设MPO球面曲率半径为Rm,X射线点光源物距为S,则焦距f可以根据龙虾眼物像公式[10]可得:

1f-1S=2Rm(24)

如果龙虾眼MPO应用于天文X射线领域,物距默认为S=¥,则平行X射线聚焦成像在焦距f=Rm/2探测器焦平面上。

表 1. MPO透镜参数

Table 1. Parameters of MPO

ParameterExperiment /Simulation
Dimension /(mm×mm)40×40
Channel dimension /μm20
Wall thickness of channel /μm6
Length of channel /mm1.25
Radius of MPO /mm750
Distance of MPO /mm13300

查看所有表

3 实验装置

由于透镜几何参数的限制,常规光学检测手段会产生干涉和衍射现象而无法使用。目前,检测透镜质量最有效的方法是使用准平行X射线光束。为了研究MPO透镜的微孔统计特性,在国家天文台X射线成像实验室,对镀金属反射膜Ir前、后的MPO透镜分别进行X射线聚焦成像测试。该X射线束流系统主要由50 kV的高压电源、50 W的Ti靶X射线管(能量为4.5 keV,微焦斑为50 μm)、六维精密装置台、大型真空腔体(长为13300 mm,真空度为10-3Pa)、MPO透镜(参数见表1),以及互补金属氧化物半导体(CMOS)探测器(图像分辨率为2048 pixel×2040 pixel,成像像素为11 mm)组成,如图3所示。

首先,通过六维精密装置台和激光发生器将X射线光源、MPO透镜及CMOS探测器调整为一个共轴光学系统;然后,经真空腔体约束的准平行X射线束照射到MPO透镜的凸面,进而会聚成十字图像;最终,被放置在焦距处的CMOS探测器接收。测试均在电压为8 kV、电流为200 μA的条件下进行,曝光时间为512 ms。

图 3. MPO透镜的测试系统。(a)测试平台示意图;(b)测试平台实物图;(c)镀膜MPO透镜

Fig. 3. Measuring system for MPO lens. (a) Schematic of experimental setup; (b) photograph of system; (c) physical map of coating MPO lens

下载图片 查看所有图片

4 结果与讨论

4.1 焦距

MPO透镜的模拟参数见表1,方孔数目约2.6×107,光线约为109根,方孔内壁粗糙度为1 nm,光源能量为4.5 keV。MPO透镜相对于CMOS探测器的距离与成像质量直接关联,未镀Ir膜的MPO透镜在不同位置处的模拟与实验成像结果如图4所示。可以看出,X射线聚焦成像结果与模拟图像变化趋势一致,少数光线透射形成背景噪声,大部分参与反射的光线被MPO透镜聚焦为十字线。在焦距为365 mm处,聚焦图像最清晰,焦斑亮度最大。偏离焦距位置处,焦斑光强迅速衰减,十字臂由2条变成了4条,中心亮斑分裂成4个,像质恶化十分明显,与理论模拟结果相吻合。模拟运算时间约15 min。该数值模型大幅度提高了运算效率,节约了模拟时间。

图 4. 不同位置处的MPO聚焦成像结果。(a)模拟;(b)实验

Fig. 4. Focused imaging results at different positions. (a) Simulation; (b) experiment

下载图片 查看所有图片

图5为未镀Ir膜的MPO透镜在焦距处的光强立体分布图,可以看出,实验与模拟结果相吻合,都服从高斯分布。MPO聚焦成像的十字线长度约22 mm,从中心位置到镜片边缘由强逐渐变弱。主要来自于两方面:一方面是X射线在微孔通道内壁反射时,需满足掠入射全反射条件,距离坐标中心的距离越远,全反射角度越大,X射线反射效率越低;另一方面是由于通道角的限制,阻挡比例与中心位置呈线性相关,距离越远,阻挡比例越大,参与反射的X射线由全部变为部分。

图 5. MPO透镜聚焦光斑分布立体图。(a)模拟;(b)实验

Fig. 5. Space diagram of focal spot distribution. (a) Simulation; (b) experiment

下载图片 查看所有图片

4.2 传输效率

为了提高MPO透镜的传输效率,采用原子沉积技术(ALD)[20]在方孔内壁表面沉积一层20 nm的金属反射膜Ir。反射膜的覆盖完整性与一致性对MPO的各项性能都有影响,使用扫描电子显微镜(SEM)对沉积膜层进行测试,在同一条通道随机抽选3个位置处的膜厚进行测试,图6是方孔通道SEM剖面图。测试厚度分别为20.4,21.8,20.6 nm,测试膜厚差别已经与测试设备的测试精度相当。SEM测试结果表明,采用ALD技术在微孔内壁镀制金属反射膜的完整性和一致性较好。

为了验证模拟的准确性,对镀膜前、后透镜的传输效率进行测试,实验结果如图7所示。透镜镀膜前、后的传输效率分别为1.23%和9.18%,模拟结果的传输效率分别为1.44%和10.14%。可以看出,模

拟结果与实验结果基本吻合。实验结果表明MPO透镜镀Ir膜后,掠入射临界角增大,X射线反射率提高,能够会聚到焦斑中心的X射线光子数目增多,焦斑的强度增大,传输效率提高。

图 6. 方孔内壁Ir膜扫描电子显微镜图。(a) MPO方孔通道剖面图;(b)所选区域的放大图

Fig. 6. SEM images of MPO square channel with coating Ir. (a) Cross section of MPO square channel; (b) enlarged view of selected area

下载图片 查看所有图片

图 7. X射线聚焦成像结果。(a)未镀Ir膜;(b)镀Ir膜

Fig. 7. Focused imaging results of X-ray. (a) Glass; (b) coating Ir

下载图片 查看所有图片

4.3 角分辨率

MPO透镜应用于X射线聚焦成像,需要考虑透镜的空间分辨能力,通常采用角分辨率来表征MPO透镜的空间成像分辨能力。理想情况下MPO透镜的半高宽(FWHM)为0.02 mm,成像角分辨率为11″。图8为MPO透镜中心焦斑光强分布曲线,未镀Ir膜的MPO曲线最强峰右侧有不弱的伴峰,可能是因为MPO透镜测试区域存在一定比例的方孔通道指向偏差,导致焦斑高斯曲线右侧存在伴峰,镀膜后中心焦斑光强增益提高约7倍。经计算得到镀Ir膜前、后的FWHM分别为0.44 mm和0.64 mm,对应的角分辨率为4.14'和6.02',与理论值存在偏差。这是因为MPO透镜在制作过程中存在很多工艺缺陷[21-22],导致方孔结构发生畸变。镀膜后透镜角分辨率变差,可能是因为MPO透镜经过镀膜高温过程(约300 °C),球面面型发生畸变,微孔统计特性变差,导致十字线弥散、角分辨率增大。

图 8. MPO透镜中心焦斑光强分布曲线

Fig. 8. X-ray intensity distribution curve of focal spot

下载图片 查看所有图片

5 结论

基于X射线全反射原理和旋转坐标系的方法,建立了X射线在方孔内壁的传输模型,并利用国家天文台X射线束流成像装置,在能量为4.5 keV下,对镀Ir膜前、后的MPO透镜进行了聚焦成像测试。实验结果表明,MPO透镜的焦距和X射线传输效率与模拟结果相吻合,表明构建的数值模型合理。镀Ir膜有助于提高掠入射临界角和X射线反射率,中心焦斑增益增大约7倍,传输效率由1.23%提升至9.18%。但镀膜工艺也可能会引起球面面型产生畸变,角分辨率由4.14'增大为6.02',降低成像质量。该数值模型可推广应用于X射线在任何形状的微孔传输特性模拟,有助于指导MPO透镜的研制。

参考文献

[1] Angel J R P. Lobster eyes as X-ray telescopes[J]. The Astrophysical Journal, 1979, 233: 364-373.

[2] Ables J G. Fourier transform photography: a new method for X-ray astronomy[J]. Publications of the Astronomical Society of Australia, 1968, 1(4): 172-173.

[3] 王飞翔, 邓彪, 王玉丹, 等. 基于毛细管分光的同步辐射X射线立体成像[J]. 光学学报, 2016, 36(8): 0834004.

    Wang F X, Deng B, Wang Y D, et al. Synchrotron radiation X-ray stereo imaging based on capillary beam splitting[J]. Acta Optica Sinica, 2016, 36(8): 0834004.

[4] Petre R, Serlemitsos P J. Conical imaging mirrors for high-speed X-ray telescopes[J]. Applied Optics, 1985, 24(12): 1833-1837.

[5] 李磊, 王林元, 蔡爱龙, 等. 基于等值线拟合的X射线双能CT快速投影分解算法[J]. 光学学报, 2016, 36(8): 0834001.

    Li L, Wang L Y, Cai A L, et al. Projection decomposition algorithm for X-ray dual-energy computed tomography based on isotransmission line fitting[J]. Acta Optica Sinica, 2016, 36(8): 0834001.

[6] Gondoin P, Aschenbach B, Brauninger H W, et al. X-ray performance of a qualification model of an XMM mirror module[J]. Proceedings of SPIE, 1996, 2808: 390-401.

[7] Priedhorsky W C, Peele A G, Nugent K A. An X-ray all-sky monitor with extraordinary sensitivity[J]. Monthly Notices of the Royal Astronomical Society, 1996, 279(3): 733-750.

[8] Brunton A N, Fraser G W, Lees J E, et al. Metrology and modeling of microchannel plate x-ray optics[J]. Applied Optics, 1997, 36(22): 5461-5470.

[9] Peele A G. Investigation of etched silicon wafers for lobster-eye optics[J]. Review of Scientific Instruments, 1999, 70(2): 1268-1273.

[10] Beijersbergen M, Kraft S, Bavdaz M, et al. Development of X-ray pore optics: novel high-resolution silicon millipore optics for XEUS and ultra-low mass glass micropore optics for imaging and timing[J]. Proceedings of SPIE, 2004, 5539: 104-115.

[11] Chapman H N, Nugent K A, Wilkins S W. X-ray focusing using square channel-capillary arrays[J]. Review of Scientific Instruments, 1991, 62(6): 1542-1561.

[12] Hudec R, Nemcová s. Effective algorithm for ray-tracing simulations of lobster eye and similar reflective optical systems[J]. Experimental Astronomy, 2016, 41(3): 377-392.

[13] Peele A G, Lyngsjo H, Crocker R M, et al. Modeling of the Lobster-ISS X-ray telescope in orbit[J]. Proceedings of SPIE, 2004, 5488: 232-241.

[14] Putkunz C T, Peele A G. Detailed simulation of a Lobster-eye telescope[J]. Optics Express, 2009, 17(16): 14156-14165.

[15] Zhao D H, Zhang C, Yuan W M, et al. Ray tracing simulations for the wide-field X-ray telescope of the einstein probe mission based on Geant4 and XRTG4[J]. Proceedings of SPIE, 2014, 9144: 91444E.

[16] 彭诗棋, 刘志国, 孙天希, 等. 基于多毛细管X光透镜的数值模拟[J]. 光学学报, 2015, 35(2): 0234001.

    Peng S Q, Liu Z G, Sun T X, et al. Numerical simulation of polycapillary X-ray lens[J]. Acta Optica Sinica, 2015, 35(2): 0234001.

[17] Pantell R H, Chung P. Transmission of X-rays through curved waveguides[J]. IEEE Journal of Quantum Electronics, 1978, 14(9): 694-697.

[18] Mitsuighi I, Ezoe Y, Takgi U, et al. Novel ultra-lightweight and high-resolution MEMS X-ray optics[J]. Proceedings of SPIE, 2009, 7360: 73600C.

[19] Spiller E A. Soft X-ray optics[J]. Optics & Photonics News, 1993, 7(2): 60-64.

[20] Windover D, Armstrong N, Cline J P, et al. Characterization of atomic layer deposition using X-Ray reflectometry[J]. Fresenius Journal of Analytical Chemistry, 2005, 788(1/2/3): 144-146.

[21] 李旭, 何飞, 李达, 等. 方孔微通道板结构缺陷对成像质量的影响[J]. 光学学报, 2012, 32(10): 1011001.

    Li X, He F, Li D, et al. Influence of structural defect of square-hole micro channel plate on imaging quality[J]. Acta Optica Sinica, 2012, 32(10): 1011001.

[22] Willingale R, Pearson J F, Martindale A, et al. Aberrations in square pore micro-channel optics used for X-ray lobster eye telescopes[J]. Proceedings of SPIE, 2016, 9905: 99051Y.

金戈, 黎龙辉, 张臣, 杨雪, 张诗鲲, 张智勇, 毛汉祺, 徐昭, 张振, 姜博文. 基于Angel型龙虾眼X射线透镜的数值模拟[J]. 光学学报, 2018, 38(9): 0934001. Ge Jin, Longhui Li, Chen Zhang, Xue Yang, Shikun Zhang, Zhiyong Zhang, Hanqi Mao, Zhao Xu, Zhen Zhang, Bowen Jiang. Numerical Simulation of Angel Lobster Eye X-ray Micro Pore Optic Lens[J]. Acta Optica Sinica, 2018, 38(9): 0934001.

本文已被 4 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!