光学学报, 2021, 41 (7): 0730003, 网络出版: 2021-04-11   

三级次太阳极紫外窄波段无狭缝成像光谱仪光学设计 下载: 687次

Optical Design of 3-Order Solar Extreme Ultraviolet Narrowband Slitless Imaging Spectrometer
邢阳光 1,2,3李林 1,*彭吉龙 2,**王姗姗 1,***成一诺 4
作者单位
1 北京理工大学光电学院, 北京 100081
2 北京卫星环境工程研究所, 北京 100094
3 可靠性与环境工程技术国家重点实验室, 北京 100094
4 清华大学精密仪器系, 北京 100091
摘要
对太阳进行光谱成像观测是太阳物理和空间天气研究的重要数据来源。在极紫外波段对动态的太阳大气进行光谱成像观测面临着严峻挑战。传统的狭缝式成像光谱仪受瞬时视场的限制,必须通过耗时的推扫才能获取二维日面图像,系统不具有高时间分辨率,无法捕获太阳过渡区域和日冕的快速演化过程。尽管极紫外成像仪能够实现大二维视场和高时间分辨率的观测,但却无法获取光谱分辨率信息。本文基于像差校正的椭球面变线距光栅,提出了一款新型无狭缝成像光谱仪,该系统同时工作在三个光栅衍射级次(m=-1,0,1)上,这样的新颖设计不需要任何元件的机械运动,单次快照便可同时获得三幅视场为20 arcmin×20 arcmin的极紫外窄波段(29.4~31.4 nm)太阳日面图像,该仪器的时间分辨率高且视场大。不发生色散的0级次系统相当于一台极紫外成像仪,可以直接获得高分辨率空间信息(0.6 arcsec);发生色散的±1级次图像携带有空间和光谱的混叠信息,类似于计算层析成像原理。通过数据反演算法,可以从三个级次的图像中提取高分辨率的光谱信息(0.0035 nm)。
Abstract
The observation of solar spectral imaging provides an importance data source for solar physics and space weather research. However, it faces a serious challenge from the dynamic solar atmosphere in the extreme ultraviolet (EUV) band. Only by time-consuming push-broom imaging can the traditional slit-type imaging spectrometer acquire 2D images of the solar disk due to a limited instantaneous field of view (FOV). Moreover, its system is unable to capture the rapid evolution of the solar transition region and corona since it does not have a high time resolution. Although an EUV imager can observe at a large 2D FOV and high time resolution, it cannot get the spectral resolution. For this reason, relying on the aberration-corrected elliptical varied line-space (EVLS) grating, we proposed a new slitless imaging spectrometer in this paper, which can simultaneously work in three diffraction orders (m=-1, 0, or +1). Besides, three EUV narrowband (29.4--31.4 nm) images of the solar disk with a FOV of 20 arcmin×20 arcmin can be obtained in a single snapshot without the mechanical motion of any component, which demonstrates the high time resolution and large FOV of our instrument. Additionally, the zero-order system without spectral dispersion, equivalent to an EUV imager, can directly obtain high-resolution spatial information (0.6 arcsec). In contrast, the images from the ±1-order systems with spectral dispersion carry spatial and spectral overlapping information, which is similar to the imaging principle of computed tomography. High-resolution spectral information (0.0035 nm) can be extracted from the 3-order images by a data inversion algorithm.

1 引言

太阳是唯一一颗可以在所有参数(包括时间、空间和波段)上进行高分辨率观测的恒星[1]。在太阳上层的大气中,具有高度动态性的日冕和过渡区域会辐射出丰富的极紫外(EUV: 10~120 nm)谱线[2]。空间太阳极紫外光谱学观测(包括成像与分光光谱测量)和诊断是研究太阳大气各区域基本物理过程的最重要手段之一,可以研究日冕物质抛射[3](CMEs)、过渡区爆发事件[4](EEs)、磁重联[5]、耀斑[6]等爆发活动及其伴随现象的发生、发展过程和形成机制。此外,太阳活动是空间天气变化的源头,可能会引起日地环境的剧烈变化,对航天、通信、导航等产生重要影响[7],可以通过空间太阳观测来研究太阳活动与地磁暴、电离层扰动等地球空间天气的因果链关系。

太阳大气在空间结构上具有高度的动态性,表现为在极短的时间尺度内发生快速的结构演化,例如典型的EEs的时间尺度为60~90 s,红翼和蓝翼的演化速度为100~400 km/s[8-9],而在一些特殊的太阳活动(CMEs、耀斑爆发等)期间,空间结构的演化会更加剧烈。因此,为了最大化地获得具有高度动态变化的太阳大气在极紫外波段的空间和光谱信息,理想的极紫外观测仪器应该具有高时间分辨、高空间分辨、高光谱分辨和大视场(二维日面)的性能。但是,要在二维日面内实现这样的观测,现阶段还面临着严峻挑战。目前,对太阳活动在极紫外波段的二维光谱成像的观测主要通过狭缝扫描式成像光谱仪来实现,然而这种仪器有其自身的短板,缺乏一些研究太阳大气动态演化所需的重要信息。

狭缝式成像光谱仪通过窄狭缝限制视场,利用凹面光栅和平面探测器,可以同时对多组极紫外谱线和沿狭缝一维空间的消像散成像(即实现一维方向的高空间和高光谱分辨的同时观测),通过狭缝的扫描或者前置望远系统的旋转可以实现对二维日面的高空间和高光谱分辨的不同时观测。典型的狭缝式成像光谱仪有2006年日本发射的Hinode卫星上的极紫外仪器EIS[10]、2013年NASA发射的极紫外探空火箭仪器EUNIS-13[11]、2020年2月由欧洲航天局和NASA联合发射的太阳极区探测器Solar Orbiter上搭载的极紫外仪器SPICE[12-13],以及邢阳光等[14]设计的太阳极紫外三波段成像光谱仪。根据谱线Doppler效应(红移或蓝移),可以得到等离子体运动的视向速度;根据谱线对的强度比值和不同谱线的轮廓宽度等信息,还可以计算等离子体温度和电子密度等。但狭缝扫描式成像光谱仪具有如下缺陷:1)狭缝的存在限制了仪器的视场,很难捕获CMEs、EEs、耀斑等爆发活动的源区,只有狭缝恰好在爆发源区位置时才能获得有效的观测;2)对于太阳活动区域的二维观测依赖于狭缝的扫描,由于曝光积分时间的限制,完成二维日面扫描需要数小时甚至十几个小时的时间,不能得到高时间分辨率的信息,无法实时观测二维日面内太阳活动的快速演化过程;3)需要在考虑光谱仪狭缝扫描速度的情况下重构太阳活动现象的速度图,而通过狭缝扫描获取的二维空间信息和光谱信息在时间上不同步,不同时间扫描得到的图像可能会提供错误的信息,例如可能会将形态在时间上的演化误认为是空间上传播的波。

为了在大视场范围内实现极紫外波段二维日面的高空间分辨、高光谱分辨的同时成像,本文设计了一款面向卫星遥感应用的新型成像光谱仪,该仪器抛弃了传统成像光谱仪用于视场扫描的狭缝,采用具有三个衍射级次(m=-1,0,+1)的无狭缝成像结构。该仪器仅由一个离轴双曲面主镜、一个椭球面变线距(EVLS)次镜光栅和三个级次的平面探测器构成,主镜和次镜表面使用相同的周期性SiC/Mg多层膜,以提高极紫外下的反射效率。

2 三级次无狭缝光谱成像的概念

传统的狭缝式成像光谱仪通过一个冗长的连续扫描来重构具有高空间和高光谱分辨信息的太阳极紫外图像的三维(附加时间轴)数据立方体I(x,y,λ,t),但通过这样的成像和光谱观测获取的I(x,y,λ)不具有同时性,空间信息I (x, y)和时间信息I (t)严重混叠。为了改进传统仪器的观测方式,2001年,Kankelborg等[15]首次提出了多级次无狭缝的概念,并于2006和2015年分别成功研制出多级次太阳极紫外成像光谱仪MOSES。2006年,MOSES随探空火箭首次升空,经过不到5 min的飞行,得到了27幅图像。Fox等[16]首次发现了过渡区爆发事件在He II 30.4 nm处的响应,成功反演了二维日面的速度场分布,这样的观测通过传统仪器是无法做到的。MOSES采用的是Wadsworth结构,去除了限制视场的狭缝,其光学系统由球面光栅主镜、平面次镜和三个平面探测器组成。入射到主镜光栅上的太阳极紫外光线经过光栅的三个衍射级次色散并通过平面镜的光路折叠后分别聚焦成像在三个探测器上,如图1(a)所示。MOSES可以实现对10 arcmin×20 arcmin日面区域的窄波段高分辨率观测,空间分辨率为0.6 arcsec,光谱分辨率为0.0030 nm,其他技术指标如表1所示。单次快照(<10 s)可获得三幅太阳极紫外图像,如图1(b)所示(FOV:filed of view)。

图 1. 三级次无狭缝成像光谱仪的示意图。(a)MOSES的光学布局;(b)单次快照获得的三级次太阳图像

Fig. 1. Schematic of three-order slitless imaging spectrometer. (a) Optical layout of MOSES; (b) solar images with three diffraction orders by single snapshot

下载图片 查看所有图片

表 1. 太阳极紫外无狭缝成像光谱仪的技术指标

Table 1. Specifications of the solar EUV slitless imaging spectrometer

InstrumentWavelength /nm2D FOV /(arcmin×arcmin)Spatialresolution /arcsecSpectralresolution /(10-4 nm)PrimarymirrorSecondarymirrorLength /mm
MOSES-0629.4--31.410×200.630SULS gratingFold flat2400
Our design29.4--31.420×200.635HyperbolaEVLS grating910

查看所有表

尽管多级次无狭缝成像光谱仪的成像特点使其同时具有高时间、空间和光谱分辨率的性能,但MOSES的设计存在以下主要局限性:

1)仪器尺寸的限制。为了获得高空间分辨率,MOSES具有接近2.4 m的轴向长度;为了获得高光谱分辨率以及避免探测器与主镜光栅的干涉,MOSES的色散平面不得不正交于纸面,导致其在色散方向上具有0.56 m的宽度。这样的仪器很难实现基于卫星平台的应用。

2)缺少像差的校正。由于平面次镜不具有校正像差的功能,而球面等线距光栅无法同时实现三个衍射级次的最优像差校正,因此MOSES的±1级和0级具有显著差别的点扩散函数(PSFs)。另外,球面光栅无法实现大视场下的离轴像差校正,仪器的视场受限。

3)缺少谱线“污染”的控制。由于MOSES没有中间像所形成的视场光阑,视场外的极紫外目标谱线以及非目标谱线也衍射到像面上,对视场内的正常谱线成像造成了“污染”,因此不同衍射级次携带的空间和光谱信息存在差异,增大了光谱数据反演的难度[15]

为了解决MOSES在设计上的局限性,本文提出了新的三级次无狭缝成像光谱仪设计方案,如图2所示(λ为衍射波长)。望远主镜采用离轴双曲面镜,次镜为像差校正的椭球面变线距光栅,望远系统形成太阳二维日面的一次像于视场光阑上,视场光阑上的极紫外光线照明椭球面光栅,经光栅的三个级次在纸面方向衍射后,聚焦成像在三个平面探测器上。采用这样的成像分光元件和仪器布局结构,成像光谱仪不仅具有与MOSES系统同等的空间和光谱分辨率,如表1所示,而且获得了2倍的观测视场(20 arcmin×20 arcmin)、小于0.4倍的轴向长度(约 910 mm)和小于0.6倍的宽度(约310 mm)。新的光学设计方案显著压缩了系统的体积,且各个衍射级次具有优良的成像质量,可应用于卫星遥感,实现对二维日面的持续性快照光谱成像,获得前所未有的太阳活动的新信息。

图 2. 具有新光学布局的三级次无狭缝成像光谱仪示意图

Fig. 2. Schematic of three-order slitless imaging spectrometer with new optical layout

下载图片 查看所有图片

3 设计原理

为了解析出仪器的基本光路和像差组成,建立了如图3所示的两个笛卡尔坐标系统,其中:WUL坐标系统(L轴垂直纸面向里)用于前置离轴望远系统的设计;XYZ坐标系统(Z轴垂直纸面向外)用于EVLS光栅像差的校正,X轴为光栅法线且相交于光栅顶点O,Z轴平行于光栅刻线方向,Y轴垂直于光栅刻线方向。口径为D的离轴双曲面镜(离轴量为Δ),会聚平行的极紫外光线于视场光阑处,中心光线过视场光阑的中心和光栅顶点O并与U轴相交于A0点,中心光线A0O经EVLS光栅的三个级次色散后,成像在具有不同倾斜角度的三个平面探测器上。

图 3. 三级次无狭缝成像光谱仪的光线追迹模型

Fig. 3. Ray-tracing model for three-order slitless imaging spectrometer

下载图片 查看所有图片

3.1 前置望远系统的设计

前置望远系统采用离轴双曲面镜。图4为离轴双曲面反射镜在WU面内的几何光学模型,双曲面顶点位于坐标系原点T处,半视场角为ω,位于口径边缘处的两平行极紫外光线E1V1E2V2照明离轴双曲面,光线EV为口径内的任意光线,其在双曲面上的入射角为α。在WU面内,光线EVW轴的夹角为αi,任意点V位于双曲线上,V点处的切线VMW轴的夹角为α,法线VNW轴的夹角为α,反射光线VAU轴于A点,且其与W轴的夹角为γ

图 4. 离轴双曲面反射镜的几何光学模型

Fig. 4. Geometric optics model for off-axis hyperbolic mirror

下载图片 查看所有图片

离轴双曲面在WU面上的方程为

U=W2RT+RT2-(1+k)W2,(1)

式中:RT为双曲面的曲率半径;k为二次项系数,k<-1。图4中的各角度满足如下关系(单位均为度):

αi=90-ωα+α=9090+γ=2α+ω(2)

曲线上任意点V处的切线斜率为

tanα=dUdW=WRT2-(1+k)W2(3)

光线EV的反射光线为VA,其斜率为

Kω=-tanγ=-cot(2α+ω)=-cot2arctanWRT2-(1+k)W2+ω(4)

由(2)式可知:光线在离轴双曲面上的入射角度α=α+ω,具体表达式为

α=arctanWRT2-(1+k)W2+ω,Δ-D2WΔ+D2(5)

对于任意的离轴双曲面镜和入射光线,根据(4)式和(5)式即可追迹反射光线的聚焦位置。当光线EV为中心光线且ω = 0时,光线标记为E0V0,其延长线交W轴于H0点,此时可知:离轴量Δ=|TH0|,V0点的坐标为(Δ,U0)。由(4)式可知V0点处反射光线的斜率为

K0=-cot2arctanΔRT2-(1+k)Δ2(6)

V0点处的反射光线交U轴于A0点,A0点的坐标为(0,U0-ΔK0),则可知图3fT的表达式为

fT=|V0A0|=Δ1+K02(7)

由(6)式和(7)式可得前置离轴望远系统的焦距。

图3中,若|A0O|=rA,|B0O|=rB且整个成像光谱仪系统的焦距为f,系统的横向放大率为β,EVLS光栅在空间成像方向上的焦距为fG,则根据近轴光学原理可知

f=-fTfG/υυ=rA-fG1/rA-1/rB=1/fGβ=rB/rAf=βfT(8)

(8)式表明:成像光谱仪系统的焦距仅与光栅的横向放大率及望远系统的焦距有关。

3.2 EVLS光栅像差的校正

Huber等[17]、Harada等[18]、薛庆生[19]解析了基于轴上点光源照明的凹面光栅的像差,这样的像差校正已被广泛用于狭缝式成像光谱仪中[20],用来获得优良的轴上点空间和光谱成像性能。但是由于离轴光栅像差的存在,这样的理论无法被应用到扩展光源照明下的凹面光栅上。如图3所示,本文设计的无狭缝成像光谱仪通过前置望远系统像面上的视场光阑照明凹面光栅,具有显著的离轴特性。鉴于此,本文解析了扩展光源照明下EVLS光栅的像差,目的是实现轴上点和离轴点的光栅像差的同时校正。

图5为单个EVLS光栅像差的校正原理图。点A0 (xa, ya, 0)为视场光阑的中心,AA0为平行光栅刻线且位于视场光阑上的任意离轴光源,A0点是点A(xa, ya, za)在XOY平面上的投影。由点A发出的发散光束照明EVLS光栅,入射光线AP为发散光束中的任意一条光线,点P(x,y,z)为光栅表面上的任意一点。入射光线APAO经过EVLS光栅衍射后,成像于点B(xb,yb,zb),B0(xb,yb,0)是点B在平面XOY上的投影。可知:轴上的点A0经过光栅顶点衍射后会成像在B0处,光线A0OX轴的夹角i为入射角,|OA0|=rA为光栅的入射臂长;光线OB0X轴的夹角θ为衍射角,|OB0|=rB为光栅的出射臂长。当iθX轴的同一侧时,符号相同,在X轴的异侧时,符号相反。

图 5. 像差校正的EVLS光栅的原理图

Fig. 5. Schematic of aberration-corrected EVLS grating

下载图片 查看所有图片

由于点P(x,y,z)在光栅表面上,因此满足空间椭球面方程

(x-a)2a2+y2b2+z2c2=1,(9)

式中:abc分别为椭球在X方向、Y方向和Z方向上的半轴长度。将(9)式中的x展开成yz的泰勒级数,得到

x=a2y2b2+z2c2+122y2b2+z2c22+123y2b2+z2c23+(10)

对于变线距的椭球面光栅,其刻线密度沿Y轴变化,总刻线数满足多项式

N(y)=y+ab2ζ2y2+ab22ζ3y3+ab23ζ4y4+d0,(11)

式中:N(y)为纵坐标原点O到任意点P的刻线总数;d0为原点O处的刻线间距;ζ2,ζ3,ζ4 为刻线密度的空间变化参数。

费马原理常被用于分析凹面光栅的像差[21],光线APB的光程函数为

F=<AO>+<OB>=<AP>+<PB>+N(y)·(12)

进一步地,将各点坐标代入到(12)式中,将光程函数展开成yz的泰勒级数,可得

F(y,z)=rA+rB+j,kyjzkFjk,Fjk=Cjk+(/d0)Mjk,(13)

式中:jk不同时为零;Mjk是与刻线密度的空间变化有关的参数;Cjk是仅与椭球面的面型参数有关的参数;Fjk的每一项对应EVLS的一个光栅像差,包括基本像差(像散和离轴离焦)和高阶像差(彗差、场曲和球差),对应关系如表2所示。

表 2. Fjk及其对应的像差

Table 2. Fjk terms with their associate aberrations

TermAberration
F10 and F01Basic grating equation
F20Tangential astigmatism
F02Sagittal astigmatism
F11Off-axis defocusing
F30 and F21Coma
F12Field curvature
F40, F22 and F04Spherical aberration

查看所有表

Mjk的表达式为

Mjk=1,j=1k=0ζj(a/b2)j-1,j=2,3,4k=00,otherwise(14)

Cjk仅与椭球面的面型参数有关,这里仅给出了表2中重要像差对应的Cjk项的具体表达式:

C01=-zarA1+za2ra2-1/2-zbrB1+zb2rB2-1/2C10=-1+za2rA2-1/2sini-1+zb2rB2-1/2sinθ,(15)C20=12cos2irA+cos2θrB-a(cosi+cosθ)b2+za22rA2sin2irA-cos2i2rA+acosi2b2+zb22rB2sin2θrB-cos2θ2rB+acosθ2b2C02=121rA+1rB-a(cosi+cosθ)c2-za24rA23rA-acosic2-zb24rB23rB-acosθc2C11=-zarA2sini-zbrB2sinθ,za=fTtanω,(16)C30=sini2rAcos2irA-acosib2+sinθ2rBcos2θrB-acosθb2C12=12sinirA1rA-acosic2+sinθrB1rB-acosθc2(17)

根据费马原理,对位于EVLS光栅表面上的任意一点P,光程函数F(y,z)都存在极值的充分条件是同时满足方程:

F(y,z)/y=0, F(y,z)/z=0(18)

对于任意的jk,Fjk=0,即所有光栅像差的系数都为零,就可以满足(18)式成立,这是获取光谱理想成像的条件。但在实际中,由所有像差系数为零构成的超越方程组没有严格的解析解。为此,本文校正EVLS光栅的主要像差。

F20=0、F02=0和F11=0,根据(13)式、(14)式和(16)式可得到校正像散和离轴离焦的条件为

cos2irA+cos2θrB-a(cosi+cosθ)b2+2maλζ2b2d0=01rA+1rB-a(cosi+cosθ)c2=0sinθrB=sinirA(19)

对于EVLS光栅,当实现像散和离轴离焦的校正后,系统的场曲被自动消除了;通过控制刻线密度变化参量ζ3ζ4,就可以分别实现对EVLS光栅彗差和球差的校正。

3.3 初始解和优化函数

为了获取无狭缝成像光谱仪的初始结构参数,需要明确一些已知的量,包括入射臂长rA、系统轴向长度L、光栅横向放大率β、消像散的波段范围λ12、中心波长λ0=(λ1+λ2)/2、离轴半视场ω、光栅中心的刻线宽度d0。联立(13)~(15)式和(17)~(19)式可得到系统的初始解为

LrA+RT/2βrAi=2arctanΔRT2-(1+k)Δ-arcsinmλ0d0(1+β)a=2βrAβ-1b2=rA×λ1(cosi+cosθ2)-λ2(cosi+cosθ1)λ1(βcos2i+cos2θ2)-λ2(βcos2i+cos2θ1)c2=rA×cosi+cosθ01+βζ2=(cosθ1-cosθ2)[b2(cosθ1+cosθ2)-rA]d02maβrA(λ2-λ1)ζ3=b2siniarAζ2θξ=arcsin8mrA2λξ/d08rA2-[RT2+(3-k)Δ2]tan2ω+sini,ξ=0,1,2,(20)

式中:θ1为边缘波长λ1的衍射角度;θ2为边缘波长λ2的衍射角度;θ0为边缘波长λ0的衍射角度。

根据(20)式获得的初始解,可以实现EVLS光栅在单个衍射级次(m=+1 )上像差的校正。但是,三级次窄波段成像光谱仪系统同时工作在+1、0和-1的衍射级次上,因此,这样的初始解不是系统的最优解。为了解决这个问题,本文在初始解的基础上,利用光学设计软件ZEMAX,通过多重结构和(21)式所示的像差优化函数,来实现对三级次成像光谱仪系统最优化像差的控制,最小化各个衍射级次的光栅像差,以获得系统的最优解。像差优化函数的表达式为

Imin=m{φ(m)τ[τ1F022(λ0)+τ2F202(λ0)+τ3F112(λ0)+τ4F302(λ0)+τ5F122(λ0)+τ6F402(λ0)]},(21)

式中:λ0是窄波段的中心波长;φ(m)是不同衍射级次的权重因子,取相同的值;τ1~τ6是不同像差系数的权重因子,取不同的值,一般高阶像差的权重因子小于低阶像差的权重因子。

4 光线追迹和仪器性能

在太阳大气的低层过渡区域中,Lyman α的单跃迁离子氦产生30.4 nm(He II)的发射谱线,形成谱线的温度接近50000 K,无论在宁静太阳期间还是在活动太阳期间,He II 30.4 nm都具有极高的谱线强度,比邻近的谱线强度至少高出一个数量级,如图6所示[22]。因此,He II 30.4 nm能够很好地被用来观测太阳过渡区域的爆发事件,进而研究该区域物质流动和能量传输的过程。本文选择29.4~31.4 nm的窄波段作为三级次无狭缝成像光谱仪的观测波段,而将中心波长He II 30.4 nm作为光谱数据反演的目标谱线。仪器的总体设计指标参见表1

图 6. SDO/EVE仪器获得的闪耀太阳光谱

Fig. 6. Sun flare spectrum from SDO/EVE

下载图片 查看所有图片

EVLS光栅具有三个面型自由优化变量和三个刻线密度控制参数,这使得其能够在大离轴视场下实现多个衍射级次的光栅像差的同时校正。为了保证ZEMAX光学设计软件所得优化结果的准确性并保证系统具有快的收敛速度,先利用(20)式计算出EVLS光栅的初始结构参数,再在ZEMAX中利用(21)式所示的优化函数对初始结构参数进行全局优化设计。在获得理想的光谱和空间聚焦性能的同时,还要对轻量化、小体积、大口径等多种应用需求进行权衡,以获得综合性能最优的系统。

4.1 仪器设计结果

利用ZEMAX中的多重结构建立EVLS光栅的+1、0和-1衍射级次,通过对三个级次的空间和光谱成像的同时优化,获得具有三个不同倾斜角度的最优探测器位置。最终优化后的光学系统光路布局如图7所示,整个系统包含7个光学元件:硅滤光片、离轴双曲面望远镜、视场光阑、EVLS光栅和三个级次的探测器。观测波段内的太阳极紫外光线照明离轴双曲面镜,在望远镜的焦平面上形成太阳的一次实像,二维的太阳一次实像经过EVLS光栅的三个级次衍射并聚焦形成二次像于探测器上。0级次探测器上的像是不发生色散的太阳二维图像,相当于一台极紫外成像仪,可直接获取空间分辨率信息;±1探测器上的像是发生色散的且空间和光谱信息混叠在一起的太阳二维图像。通过对三个级次的图像进行反演,可以提取目标谱线He II 30.4 nm的光谱信息。为了提高正入射下极紫外的反射率,离轴望远镜和EVLS光栅表面镀相同的周期性SiC/Mg多层膜,膜层设计文件被添加到ZEMAX中,用以对光谱成像质量和仪器效率进行优化。整个系统的技术指标和元件参数如表3所示。

图 7. 太阳极紫外无狭缝成像光谱仪的光路原理图。(a) 2D光路图;(b) 3D模型图

Fig. 7. Optical layout of solar EUV slitless imaging spectrometer. (a) 2D optical layout; (b) 3D model diagram

下载图片 查看所有图片

表 3. 无狭缝成像光谱仪的技术指标和光学元件参数

Table 3. Specifications and optical element parameters for slitless imaging spectrometer

Specification
Spectral range /nm29.4--31.4
FOV /(arcmin×arcmin)20×20
Spectral resolution /(10-4 nm)35
Spatial resolution /arcsec0.6 (440 km)
Line-of-sight velocity /(km·s-1)>34.5
High cadence /s<10
System focal length /mm4500
Pixel size of detectors /μm13 (2048×2048)
Envelope size of optical system /(mm×mm×mm)910×310×60
Telescope design
RT /mm1520
Conic-1.307
Δ /mm80
Spectral imaging system design
Field stop size /(mm×mm)4.5×4.5
1d0 /mm-13600
m+1 order, 0 order, and -1 order
Grating parameterInitialOptimum
rA /mm150150
β6×5.905×
i /(°)5.135.4
a /mm360.000352.906
b /mm305.010300.957
c /mm302.425300.034
ζ2-0.04448-0.04736
ζ3-0.006852-0.007201
Ruling area /(mm×mm)25×25
Three independent detectors design
OrderWavelength /nmTilt angle of detectors /(°)
-1 order29.4--31.414.66
0 order29.4--31.427.95
+1 order29.4--31.432.71

查看所有表

在EVLS光栅最终设计的结果中,光栅刻线区域的尺寸为25 mm×25 mm,其刻线密度在±12.5 mm色散方向上近似线性变化,如图8所示,光栅中心的刻线密度为3600 line/mm。

图 8. 椭球面变线距光栅刻线密度的分布曲线

Fig. 8. Ruling density distribution of EVLS grating

下载图片 查看所有图片

+1级次和-1级次子午像散和弧矢像散的像差系数随波长的变化曲线如图9所示。由图9可知:像散的像差系数均在10-5量级以下,获得了良好的像散校正效果。彗差的像差系数随波长的变化曲线如图10所示,可知:彗差的像差系数均在10-7数量级,取得了良好的彗差校正效果。通过边缘视场的最大网格畸变在不同波长下的具体值来衡量系统光谱畸变的大小,如图11所示。由图11可知-1级次具有0.6515%的最大网格畸变,+1级次具有1.5684%的最大网格畸变。最大网格畸变均远小于3%,说明本研究取得了良好的畸变校正效果。

图 9. 像散的像差系数曲线图。(a) -1级次的子午像散;(b) +1级次的子午像散;(c) -1级次的弧矢像散;(d) +1级次的弧矢像散

Fig. 9. Curve of aberration term of astigmatisms. (a) Tangential astigmatism for -1 order; (b) tangential astigmatism for +1 order; (c) sagittal astigmatism for -1 order; (d) sagittal astigmatism for +1 order

下载图片 查看所有图片

图 10. 彗差的像差系数曲线图。(a) -1级次彗差;(b) +1级次彗差

Fig. 10. Curve of aberration term of coma. (a) Coma for -1 order; (b) coma for +1 order

下载图片 查看所有图片

图 11. 边缘视场下最大网格畸变随波长变化的曲线图。(a) -1级次;(b) +1级次

Fig. 11. Maximum grid distortion change with wavelength under marginal FOV. (a) -1 order; (b) +1 order

下载图片 查看所有图片

4.2 成像质量评价

图12(a)、13(a)分别为-1和+1级次在中心和边缘视场下,空间和光谱方向上弥散斑的均方根(RMS)半径随波长变化的曲线。可以看出,在29.4~31.4 nm的消像散波段和±10 arcmin的离轴视场内,空间和光谱方向上弥散斑的RMS半径均小于6.5 μm,±1衍射级次获得了良好的空间成像和光谱消像散性能。

图12(b)、13(b)分别为-1和+1级次在中心和边缘波长下,空间和光谱方向上弥散斑的RMS半径随离轴视场的变化曲线图。可见,系统在中心和边缘波长下均获得了良好的离轴性能。

图 12. -1级次下弥散斑RMS半径的变化曲线图。(a)不同离轴视场下RMS半径随波长变化的曲线图;(b)不同波长下RMS半径随视场变化的曲线图

Fig. 12. Variation of RMS radius of spot for -1 order. (a) RMS radius of spot versus wavelength in the different off-axis FOVs; (b) RMS radius of spot versus FOV in the different wavelength

下载图片 查看所有图片

图 13. +1级次下弥散斑RMS半径的变化曲线图。(a)不同离轴视场下RMS半径随波长变化的曲线图;(b)不同波长RMS半径随视场变化的曲线图

Fig. 13. Variation of RMS radius of spot for +1 order. (a) RMS radius of spot versus wavelength in the different off-axis FOVs; (b) RMS radius of spot versus FOV in the different wavelength

下载图片 查看所有图片

图14(a)为0级次像面上的点列图,图14(b)为0级次弥散斑RMS半径随视场的变化曲线,可知,全视场范围内弥散斑的RMS半径均小于4.5 μm,0级次成像性能优异。

图 14. 0级次像面上的点列图以及0级次弥散斑RMS半径随视场变化的曲线图。(a)不同视场下0级次像面上的点列图;(b) 0级次弥散斑RMS半径随视场变化的曲线图

Fig. 14. Spot diagrams in 0 order imaging plane under the different FOVs and RMS radius of spot versus FOV for 0 order. (a) Spot diagrams in 0 order imaging plane under the different FOVs; (b) RMS radius of spot versus FOV for 0 order

下载图片 查看所有图片

图15为成像光谱仪在三个衍射级次下中心波长处的调制传递函数(MTF)曲线,可见,MTF在奈奎斯特频率38.46 lp/mm下均大于0.6,MTF接近衍射极限,说明该系统具有良好的成像质量。

图 15. 光学系统在不同衍射级次下的调制传递函数。(a) -1级次;(b) +1级次;(c) 0级次

Fig. 15. MTFs of design system under different diffraction orders. (a) -1 order; (b) +1 order; (c) 0 order

下载图片 查看所有图片

系统的成像质量评价表明,无狭缝成像光谱仪在三个衍射级次上同时获得了优良的成像质量。由于EVLS光栅是非球面面型并具有5.9倍的横向放大率,因此系统不仅具有良好的空间和光谱成像性能,还具有显著减小的体积(相比于MOSES仪器2400 mm×560 mm×80 mm的探空火箭布局,该系统具有910 mm×310 mm×60 mm的紧凑结构,光学系统的包络尺寸减小了84.3%,适合在遥感卫星平台上应用)。

4.3 空间和光谱分辨率

采用圈入弥散斑80%能量的包围圆直径对应的空间分辨率来评价系统的空间分辨性能。三个衍射级次的包围圆能量曲线如图16所示,可见:在全视场下,-1级次、+1级次和0级次在中心波长30.4 nm处的包围圆能量达到80%时,包围圆直径分别为12.19,12,9.4 μm,其直径对应的空间分辨率分别为0.559,0.550,0.431 arcsec;当包围圆直径为13 μm时,圈入能量分数均大于85%,尽管系统的空间分辨率受限于像元尺寸,但仍可以达到0.596 arcsec。

图 16. 用于评估系统空间分辨率的衍射包围圆能量。(a) -1级次;(b) +1级次;(c) 0级次

Fig. 16. Diffraction enclosed energy to evaluate system’s spatial resolution. (a) -1 order; (b) +1 order; (c) 0 order

下载图片 查看所有图片

无狭缝成像光谱仪系统的光谱分辨率取决于EVLS光栅在±1衍射级次的设计,图17为间隔0.0035 nm的三个波长在±1级次像面上的点列图分布,可见:在中心和边缘视场下,±1级次像面上的点列图被明显分开了,系统具有0.0035 nm的光谱分辨率。

图 17. ±1衍射级次的光谱分辨率仿真验证。(a)(b) -1级次;(c)(d) +1级次

Fig. 17. Spectral resolution simulation verification for ±1 diffraction orders. (a)(b) -1 order; (c)(d) +1 order

下载图片 查看所有图片

尽管成像光谱仪具有0.0035 nm的光谱分辨能力,但对二维日面成像时,由于没有窄狭缝对视场的限制,±1级次像面上的空间信息和光谱信息会产生混叠,如图2所示,λ1λ2这两个不同波长的光谱像在空间上或不同视场上产生了混叠。这样的成像方式类似于层析成像原理,通过对三个投影级次的病态反演,可以提取出光谱的分辨率信息,而空间分辨率信息则可以从0级次的成像中直接获取。

三级次无狭缝成像光谱仪不具有宽波段的覆盖能力,这种新型的光谱成像方式,仅仅适用于窄波段,这是它的劣势,同时也是优势。相比于传统狭缝式成像光谱仪具有的宽波段覆盖能力,这无疑是劣势;但是越窄的波段,越有利于光谱数据的反演,进而获取较高的反演精度。针对三级次无狭缝成像光谱仪数据的反演,Fox等[23]提出了傅里叶后投影和SMART算法,Oktem等[24]提出了基于参数估计的CRB算法,这些反演算法具有超过30%的系统误差,最根本的原因是用于光谱数据反演的投影级次太少,即光栅的衍射级次太少,这也是MOSES和本文设计的三级次成像光谱仪最显著的不足之处。为了降低光谱数据反演的系统误差,更多的光栅衍射级次是必要的,比如采用具有5个衍射级次(±1、0、±2)的无狭缝光谱成像结构。

4.4 仪器的传输性能

由于可见光的辐照强度约是极紫外辐照强度的100万倍,因此可见光必须被抑制。选用0.1 μm厚的硅滤光片来抑制可见光,其在30.4 nm处的透过率约为0.65,其透过率曲线如图18所示。

图 18. 0.1 μm厚硅滤光片的透过率曲线

Fig. 18. Transmission of Si filter with 0.1 μm thickness

下载图片 查看所有图片

为了提高仪器在观测波段内正入射下的极紫外反射率,同时抑制图6所示的Fe XV 28.4 nm和Fe XVI 33.5 nm这两个非观测波段内的强极紫外线,周期性SiC/Mg多层膜被应用在离轴反射镜和EVLS光栅上。SiC/Mg多层膜反射率随波长的变化曲线如图19所示,它在30.4 nm处具有0.45的峰值反射率。设计的周期性SiC/Mg多层膜的参数如表4所示,其中γ为正入射角度,d为一个周期的膜厚,τ为一个周期内SiC层厚度的占比, N为膜层的周期数,δSiC-Mg为SiC层在镁层界面上的粗糙度,δMg-SiC为镁层在SiC层界面上的粗糙度。

图 19. SiC/Mg多层膜随波长变化的反射率曲线

Fig. 19. Reflectance versus wavelength for SiC/Mg multilayer

下载图片 查看所有图片

表 4. 周期性SiC/Mg多层膜的参数

Table 4. Periodic SiC/Mg multilayer parameters

γ /(°)d / nmτNδSiC-Mg /nmδMg-SiC /nm
84.615.270.22402.01.0

查看所有表

本文在COSMOL软件中建立了EVLS光栅的面型、刻槽轮廓、刻线分布和SiC/Mg多层膜。图20为仿真计算得到的多层膜EVLS光栅在±1级次和0级次的衍射效率曲线,±1级次衍射效率曲线的峰值在30.4 nm处,为0.13, 0级次衍射效率曲线的峰值在30.8 nm处,为0.11。

图 20. EVLS光栅效率随波长的变化曲线

Fig. 20. Efficiency versus wavelength for EVLS grating

下载图片 查看所有图片

极紫外波段下仪器的有效面积越大,传输效率就越高。无狭缝成像光谱仪在波长λ处的总效率为ETOT(λ),有效面积为AEFF(λ),它们的计算表达式为

ETOT(λ)=Tfilter(λ)×Rmir(λ)×ηgra(λ)×Eccd(λ)AEFF(λ)=Aape×ETOT(λ),(22)

式中:Tfilter(λ)是滤光片的透过率;Rmir(λ)是离轴反射镜膜层的反射率;ηgra(λ)是EVLS光栅的衍射效率,包含光栅刻槽效率和多层膜反射率;Eccd(λ)是探测器的量子效率; Aape是入射孔径的面积,单位为mm2。探测器选用背照式E2V CCD230-42,其在29.4~31.4 nm波段内具有优于0.4的量子效率。成像光谱仪系统的有效面积曲线如图21所示,±1级次的峰值有效面积为57 mm2, 0级次的峰值有效面积为45 mm2

图 21. 无狭缝成像光谱仪有效面积随波长的变化曲线

Fig. 21. Effective area of slitless imaging spectrometer versus wavelength

下载图片 查看所有图片

5 结论

本文提出了新的三级次无狭缝成像光谱仪,该光谱仪采用无狭缝设计,克服了传统狭缝扫描式成像光谱仪在空间和时间上的信息混叠,通过多个衍射级次同时成像和光谱数据反演算法消除了多层膜成像仪在空间和光谱上的信息混叠。整个系统具有新的结构布局,系统的光学像差被解析了。在系统初始解的基础上,通过ZEMAX的优化,新型成像光谱仪获得了紧凑的空间结构、高的空间分辨率、高的光谱分辨率和大的二维观测视场。该仪器使用的是周期性多层膜,因此具有高的传输效率。整个系统不需要任何元件的机械运动,通过单次快照便可以获得三幅携带有谱线信息的二维太阳极紫外图像,通过对三个投影级次的光谱数据进行反演,可以提取出高分辨率光谱信息。尽管三个级次的光谱数据反演具有超过30%的系统误差,但通过增加衍射投影的级次,可以减小光谱数据的反演误差,使仪器具有良好的实用性。这样的仪器可应用于卫星遥感,实现对二维日面的持续性快照光谱成像,获得前所未有的太阳活动的新信息。

参考文献

[1] 甘为群, 黄宇, 颜毅华. 太阳空间探测的过去与未来[J]. 中国科学: 物理学力学天文学, 2012, 42(12): 1274-1281.

    Gan W Q, Huang Y, Yan Y H. The past and future of space solar observations[J]. Scientia Sinica (Physica,Mechanica & Astronomica), 2012, 42(12): 1274-1281.

[2] Dwivedi B N. EUV spectroscopy as a plasma diagnostic[J]. Space Science Reviews, 1993, 65(3): 289-316.

[3] TimothyH. Coronal mass ejections: an introduction [M]. Burton W B, Bertola F, Cassinelli J P, et al. Transl. 1st ed. Beijing: Springer-Verlag Press, 2011: 8- 40.

[4] Rust T, Kankelborg C C. Quiet-Sun explosive events observed in He II λ304 with MOSES-06[J]. The Astrophysical Journal Letters, 2019, 877(1): 59.

[5] PriestE, ForbesT. Magnetic reconnection[M]. Cambridge: Cambridge University Press, 2000: 100- 136.

[6] Milligan R O. Extreme ultra-violet spectroscopy of the lower solar atmosphere during solar flares (invited review)[J]. Solar Physics, 2015, 290(12): 3399-3423.

[7] 汪景琇, 季海生. 空间天气驱动源:太阳风暴研究[J]. 中国科学: 地球科学, 2013, 43: 883-911.

    Wang J X, Ji H S. Recent advances in solar storm studies in China[J]. Science China: Earth Sciences, 2013, 43: 883-911.

[8] Innes D E, Inhester B, Axford W I, et al. Bi-directional plasma jets produced by magnetic reconnection on the Sun[J]. Nature, 1997, 386(6627): 811-813.

[9] Teriaca L, Banerjee D, Falchi A, et al. Transition region small-scale dynamics as seen by SUMER on SOHO[J]. Astronomy & Astrophysics, 2004, 427(3): 1065-1074.

[10] Culhane J L, Harra L K, James A M, et al. The EUV imaging spectrometer for Hinode[J]. Solar Physics, 2007, 243(1): 19-61.

[11] Brosius J W, Daw A N, Rabin D M. Pervasive faint Fe XIX emission from a solar active region observed with EUNIS-13: evidence for nanoflare heating[J]. The Astrophysical Journal Letters, 2014, 790(2): 112-120.

[12] Anderson M, Appourchaux T, Auchère F, et al. The Solar Orbiter SPICE instrument: an extreme UV imaging spectrometer[J]. Astronomy & Astrophysics, 2019, 24(7): 1-26.

[13] Caldwell M E, Morris N, Griffin D K, et al. The VUV instrument SPICE for Solar Orbiter: performance ground testing[J]. Proceedings of SPIE, 2017, 10397: 1039708.

[14] 邢阳光, 李林, 彭吉龙, 等. 具有非罗兰圆结构的太阳极紫外正入射宽波段成像光谱仪光学设计[J]. 光学学报, 2020, 40(23): 2312005.

    Xing Y G, Li L, Peng J L, et al. Optical design of solar extreme ultraviolet normal-incidence broadband imaging spectrometer with non-Rowland circle mounting[J]. Acta Optica Sinica, 2020, 40(23): 2312005.

[15] Kankelborg C C, Thomas R J. Simultaneous imaging and spectroscopy of the solar atmosphere: advantages and challenges of a 3-order slitless spectrograph[[J]. Proceedings of SPIE, 2001, 4498: 16-26.

[16] Fox J L, Kankelborg C C, Thomas R J. A transition region explosive event observed in He II with the MOSES sounding rocket[J]. The Astrophysical Journal Letters, 2010, 719(2): 1132-1143.

[17] Huber M C E, Tondello G. Stigmatic performance of an EUV spectrograph with a single toroidal grating[J]. Applied Optics, 1979, 18(23): 3948-3953.

[18] Harada T, Takahashi K, Sakuma H, et al. Optimum design of a grazing-incidence flat-field spectrograph with a spherical varied-line-space grating[J]. Applied Optics, 1999, 38(13): 2743-2748.

[19] 薛庆生. 星载宽波段远紫外高光谱成像仪光学系统设计[J]. 光学学报, 2013, 33(3): 0322001.

    Xue Q S. Optical system design of a spaceborne broadband far ultraviolet hyperspectral imager[J]. Acta Optica Sinica, 2013, 33(3): 0322001.

[20] 于磊, 林冠宇, 张子辉. 宽谱段改进型Wadsworth光栅型成像光谱仪光学设计[J]. 光学学报, 2016, 36(6): 0622001.

    Yu L, Lin G Y, Zhang Z H. Optical design of broadband advanced Wadsworth grating imaging spectrometer[J]. Acta Optica Sinica, 2016, 36(6): 0622001.

[21] Beutler H G. The theory of the concave grating[J]. Journal of the Optical Society of America, 1945, 35(5): 311-350.

[22] Del Zanna G, Mason H E. Solar UV and X-ray spectral diagnostics[J]. Living Reviews in Solar Physics, 2018, 15(1): 1-278.

[23] Fox J L, Kankelborg C C, Metcalf T R. Data inversion for the multi-order solar extreme-ultraviolet spectrograph[J]. Proceedings of SPIE, 2003, 5157: 124-132.

[24] Oktem FS, KamalabadiF, Davila JM. Cramer-Rao bounds and instrument optimization for slitless spectroscopy[C]∥2013 IEEE International Conference on Acoustics, Speech and Signal Processing. [S.l.:s.n.], 2013: 2169- 2173.

邢阳光, 李林, 彭吉龙, 王姗姗, 成一诺. 三级次太阳极紫外窄波段无狭缝成像光谱仪光学设计[J]. 光学学报, 2021, 41(7): 0730003. Yangguang Xing, Lin Li, Jilong Peng, Shanshan Wang, Yinuo Cheng. Optical Design of 3-Order Solar Extreme Ultraviolet Narrowband Slitless Imaging Spectrometer[J]. Acta Optica Sinica, 2021, 41(7): 0730003.

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

相关论文

加载中...

关于本站 Cookie 的使用提示

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