基于CSST多色测光的时序测光精度分析.pdf
《基于CSST多色测光的时序测光精度分析.pdf》由会员分享,可在线阅读,更多相关《基于CSST多色测光的时序测光精度分析.pdf(12页珍藏版)》请在咨信网上搜索。
1、第6 4卷第4期2023年7 月doi:10.15940/ki.0001-5245.2023.04.008天文学报ACTA ASTRONOMICA SINICAVol.64 No.4Jul.,2023基于CSST多色测光的时序测光精度分析王婉好 1,2,3(1南京大学天文与空间科学学院南京2 10 0 2 3)(2平湖实验室(北京无线电测量研究所)平湖3140 0 0)(3现代天文与天体物理教育部重点实验室南京2 10 0 2 3)摘要中国巡天空间望远镜(Chinese Survey Space Telescope,CSST)是中国的首个大型空间光学望远镜,将对包括系外行星探测在内的诸多科学目
2、标开展研究,有望取得前沿科学进展.时序测光精度是CSST重要的性能指标,受到物理噪声和仪器噪声的影响,需要通过数值模拟对其分析和评估.模拟基于目前公布的CSST主要技术参数,建立了时序的恒星信号和噪声模型,以CSST的i波段为例,分析凝视观测模式下的测光精度.通过数值仿真,展示了孔径测光中各项噪声的贡献,特别是由指向抖动和像素响应不均匀性导致的抖动噪声,模拟结果还给出了测光孔径的推荐范围.为了获得更高的信噪比,可以减小仪器抖动振幅和像素不均匀性,或者采用参考星较差测光的方式.结果为CSST后续的时序测光精度与不同指标参数的相关性分析、系外行星探测能力评估以及测光数据处理提供了模拟数据的支撑.关
3、键词望远镜,行星和卫星:探测,技术:测光,方法:数值中图分类号:P141;文献标识码:A1引言在过去的二十多年中,系外行星领域取得了巨大的进展借助Kepler、T ESS(T r a n s i t i n g Ex-oplanet Survey Satellite)2等多项空间任务,截至2023年5月16 日,已有5419 颗系外行星被证认1.凌星法是目前运用最广泛的系外行星探测方法,可以给出行星轨道和行星-恒星半径比等信息基于Gaia(Global Astrometric Interferometer for As-trophysics)3、LA M O ST (La r g e Sk y
4、 A r e a M u l t i-Object Fiber Spectroscopic Telescope)4等望远镜提供的恒星基本参数,我们可以得到行星更准确的半径、质量等信息随着系外行星样本的扩大,系2022-06-08收到原稿,2 0 2 2-0 6-2 0 收到修改稿*国家自然科学基金项目(1197 30 2 8、119330 0 1),中国载人航天工程巡天空间望远镜专项科学研究项目(CMS-CSST-2021-B12、C M S-C SST-2 0 2 1-B0 9),十四五民用航天基金项目(D050105)资助1https:/exoplanetarchive.ipac.calt
5、ech.edu/刘慧根 1,2,3外行星的大气特征刻画也成为研究的热点.对系外行星大气的研究目前主要来自凌星行星的透射光谱观测,即行星在经过其主星前方的凌星过程中,部分恒星光透过行星大气,表现出与波长相关的吸收特征通过测量凌星时的透射光谱,可以推断行星大气特性,如行星的化学成分、温度-压强结构、云霾存在等.无论系外行星探测还是行星大气刻画都需很高的测光精度.中国巡天空间望远镜(Chinese Survey SpaceTelescope,CSST)的科学目标之一是对系外行星开展多波段测光和光谱观测研究,获得系外行星系统的大气光谱,为未来寻找和确认系外生命特征信号周济林1,345-164卷奠定重要
6、基础.CSST的主光机是一个2 m口径的离轴三反消像散(three-mirror anastigmat,T M A)光学系统,后端仪器包括了多色成像和无缝光谱巡天模块、多通道成像仪、积分视场光谱仪、系外行星成像星冕仪、高灵敏度太赫兹模块等5.本文将基于CSST的多色成像测光观测,对CSST的测光精度进行模拟分析测光的精度受到多种因素影响,例如恒星表面的黑子和热斑等导致恒星表面不均匀引起的光度变化6 .更多地,测光精度受到多种物理和仪器噪声的干扰,例如天光背景、仪器指向抖动、像素不均匀性等噪声的大小会直接体现在后续测光或光谱探测的信噪比中,进而影响行星的探测以及行星大气的光谱反演7 ,因此对噪声
7、进行准确的模拟并尽可能地对其进行校正是十分必要的.此前,有一些工作针对PLATO(PLAnetary Transits and Os-cillations of stars)、T ESS和JWST(James WebbSpace Telescope)等科学仪器的恒星光变曲线进行了模拟8-10 在这里我们主要基于CSST的基本参数,对恒星信号和噪声建模,从而计算测光信噪比.该模拟方法对光谱信噪比计算同样适用,可以应用在光谱信噪比估计中。在本工作中,我们模拟了CSST的恒星信号和各项噪声,特别讨论了仪器抖动和像素不均匀性对测光观测精度的影响,为后续CSST的参数设计和数据处理模式以及CCD选型做准
8、备.第2 节描述了信号和噪声模拟的计算模型;第3节介绍了本文使用的参数和测试结果;第4节总结了本文工作并对未来可以继续进行的研究展开讨论.2沙测光精度计算方法本节主要描述CSST测光精度模拟的计算框架和算法模拟分为3个部分首先,通过恒星能谱和视星等(AB星等)计算恒星的辐射流量,并通过CSST仿真数据得到点扩散函数(point spreadfunction,PSF),考虑像素响应函数(pixel responsefunction,PR F)和仪器抖动(jitter)的影响,通过卷积获得模拟的恒星信号.随后,对泊松噪声、天光背景、暗电流、读出噪声、稳定性误差等噪声项进行建模和计算.最后,计算孔径
9、测光信噪比(signal天文学报to noise ratio,SNR).2.1(CCD上的恒星信号模拟2.1.1恒星辐射流量在望远镜入瞳的位置处,恒星波长入处辐射的流量密度Fobs()(单位为ergs-1.cm-2.nm-1)可由该波段的AB星等ms计算得到,即:Fobs(A)=10-0.4 ma Fo,其中,Fo为该波段的零星等流量.Fobs()或可由恒星模型流量密度谱Fs(入)计算得到,即:Fobs(A)=Fs(A)Rd2其中,R.为恒星半径、d为恒星距望远镜的距离.设恒星光子入射到有效面积为Atel的望远镜上,经过模块效率为(含光学元件和探测器)的调制,则CCD(c h a r g e-
10、c o u p l e d d e v i c e)单位时间内接收的电子数(单位为e-s-1)为:n),s=Fobs()Ateln/其中,h为普朗克常数、c为光速、袭为观测波段中心波长入处的光子能量、入为观测波段的带宽.2.1.2PSF模型恒星经过望远镜光路,在焦平面上形成的星象轮廓可用其PSF表征,PSF与视场位置和波长都相关.CSST项目组发布的主巡天模块Cycle 5数据产品提供了PSF静态仿真数据,在光学设计模型的基础上,考虑了镜面加工、装调、CCD不平整度、重力变化和热变形等误差因素,给出了视场不同位置处各波段的静态PSF数据.本文为了研究不同抖动参数的影响,采用了静态PSF数据,作
11、为后续模拟动态PSF(加入像素响应不均匀性和抖动的影响)的基础。CSST提供的PSF静态仿真数据,在焦平面的每个探测器上,对均匀分布的30 30 个视场位置进行了采样,以保证PSF的插值精度另外,在多色成像的每个波段内,提供了4个典型波长的采样,展示PSF对不同波长的依赖关系.在给定波长的情况下,可以根据不同视场的PSF采样,通过插值的方法得到视场内任意位置4期(1)(2)hc入,(3)45-264卷处的PSF.在CSST的PSF仿真中,采用的是IDW(inverse-distance weighting)插值1 方法,PSF场连续变化,任一位置处的PSF由其周围的4个相邻点的PSF采样插值给
12、出,插值权重与距离的平方成反比.图1显示了i波段中一个典型波长(7 0 6.1nm)在位置(0,0)处插值得到的PSF.每个PSF给出了2 56 x王婉好等:基于CSST多色测光的时序测光精度分析4期256像素空间上的能量分布,采样的像素大小为5m.图中的横纵坐标表示PSF在c-y像素平面上的位置坐标.为了提高计算效率,我们没有对不同波长的PSF进行插值累加,在观测波段较窄的情况下,可以选择典型波长近似代替整个波段的PSF(详见第3节).640.0620.041axId/A0.0200.00-210050-1005000-5050100 x/pixel图1PSF插值仿真示意图.左:全空间2 5
13、6 2 56 范围内的PSF分布图.右:12 12 范围内的PSF分布图,覆盖了目标源91.9%的光子能量。横纵坐Fig.1 Simulated PSF distribution using interpolation.Left:PSF distribution in an area of 256 256 pixels.Right:PSFdistribution in an area of 12 x 12 pixels,covering 91.9%photos energy of the target star.The positional coordinates of the PSFon t
14、he-y pixel plane are represented by the horizontal and vertical axes,labeled as and y,respectively.2.1.3像素响应函数(PRF)模型恒星流量在CCD上的数值受到PRF的影响.PRF在像素内部不同位置及不同像素之间存在差异,分别称为像素内不均匀性和像素间不均匀性.因此,目标恒星在计算流量和PSF后,还需与PRF做卷积,获得探测器上最终读数.如果望远镜指向稳定,我们可以通过测量平场来改正像素间响应的不均匀性.但结合下文2.1.4节所述的仪器抖动的影响,像素不均匀响应将与抖动误差相耦合,带来测光噪声
15、。CSST探测器性能目前尚未确定,本文采用Kepler探测器在相应波段的像素模型12 ,将每个像素划分为4545的亚像素采样,计算探测器不同像素位置处的响应情况图2 左图显示了7 0 0 nm波长处的像素内响应函数(intra-pixel response func-tion,I PR F),右图在左图响应函数的基础上给不同像素增加了1%的随机差异.0.070.060.050.0490.030.02-40.01-60.00-6-4-2024x/pixel标a、y 表示PSF在c-y像素平面上的位置坐标.2.1.4仪器抖动影响CSST在轨观测期间,由于微振动环境和精密稳像残差导致仪器会在一定频率
16、上产生小幅度的抖动,使得星象在CCD上的位置发生偏移,从而使目标源对应的PSF和PRF都发生变化,对恒星信号产生影响.为了分析抖动引起的恒星信号的时序变化,本文采用数值模拟来分析抖动叠加PRF带来的测光误差.首先通过二维随机游走原理生成抖动序列.从平面上起始坐标为(o,yo)的一点出发,每经过t的时间随机游走一步设每一步c、9方向上的分量i、y;均符合均值为0,方差为的高斯分布.由中心极限定理可知,n(n 1)步之后,即t=nt 时刻,、方向上的总分量均符合高斯分布.设at、y t 分别为t时刻抖动序列的、y分量,则at=Z(r i)、t =(A y i)均 N(0,no).位移R=VAat)
17、+(A y t)2 则服从瑞利分布,概率645-364卷密度函数为:天文学报图3显示了由CSST抖动指标(0.0 5 30 0 s)生成的30 0 s内的一组随机抖动序列.我们取时间步长t=0.05s,则单次抖动的、y 方向分量符合2nR2eno4期均值为0,标准差为oo0.0 0 7 p i x e l 的高斯分布.01axId/K10Fig.2Simulated PRF distribution.Left:Intra-pixel effect.Right:Inter-pixel effect in an area of 12 x 12 pixels.1.00Axt0.75Ayt0.500.
18、250.00-0.25-0.50-0.75-1.000图3由CSST抖动指标(0.0 5 30 0 s)生成的30 0 s内的一组随机抖Fig.3 A sample random array in 300 s for jitter simulationbased on the jitter technical parameters(0.05300 s)of结合上述仪器抖动的影响,将恒星在CCD上的位置随时间的变化描述为:t=Co+t,yt=yo+y t.012345678910111210X/pixel图2 像素响应函数分布图.左:像素内效应.右:12 12 范围内像素间效应。50100动序列
19、CSST1.00.80.60.40.20.0456789 101112x/pixel在每个时刻t,目标源的点扩散函数为(t,yt)位置处对应的PSF,且由于星象落在不同的像素上,对应的PRF也不同经过长的曝光时间tobs,最终CCD对该目标源的读数表现为积分效应,像素(c,y)处接收到的总的电子数(单位为e-)为:tobsN入,s=n入,s*PSF(ac,y)*PRF(c,y)dt,(5)Jo其中,*表示卷积,n入,s为恒星辐射对应的电子流量,150200t/s250300(4)由(3)式给出,点扩散函数PSF(c,y)和像素响应函数PRF(c,y)分别由2.1.2 节和2.1.3节给出.2.
20、2噪声计算2.2.1散粒噪声通过测量望远镜处的光强给出的是接收光子的平均数量,望远镜任一时刻接收的光子数满足泊松分布,因此散粒噪声也称之为泊松噪声(Poissonnoise).泊松噪声在粒子数足够大时趋向于正态分布,其标准差等于光子数的平方根.对于特定的观测系统和波长,探测器产生的电子数正比于其接收到的光子数,因此有:45-464卷其中,N,s由(5)式给出.2.2.2天光背景噪声影响CSST测光观测的主要天光背景来源于黄道光(Zodiacal light).Pascale等13 通过下式模拟黄道光的反射与发射分量的总和(单位为Wm-2.m-1.sr-1):Izodi()=3.5 10-14
21、B(5500K)+3.58 10-8 Bx(270K),其中,B(5500K)和B(270K)分别为550 0 K和270K温度下的普朗克黑体辐射强度,是与黄纬有关的参数,Sarkar等10 提出用如下多项式对其进行拟合:=0.22968868c7+1.12162927c61.72338015c5+1.13119022c40.95684987c3+0.2199208c2 0.05989941C+2.57035947,其中,=lg(lecl+1),lecl为黄纬(单位为).我们取的平均值为1.3,则(7)-(8)式计算出的黄道光与HST14给出的黄道光平均值2 2.7 magarcsec-2基本
22、一致.CCD单个像素、单位时间因黄道光辐射产生的电子数(单位为e-s-1pixel-1)为:hcnZodi=IzodiAtelpixn/入(9)其中,pix为像素立体角.经过时间tobs,对npix个像素进行采样,黄道光噪声可表示为:Obg=Vnzodinpixtobs:2.2.3暗电流噪声探测器主要的热噪声来源是暗电流(dark cur-rent),指探测器没有接收光子时仍会产生的微小电流.暗电流大小与温度有关,并随积分时间累积,同样服从泊松分布.王婉好等:基于CSST多色测光的时序测光精度分析OPs=VNA,s,(6)(7)(8)(10)4期设单位时间单个像素上暗电流平均大小为ndc(单位
23、为e-s-1pixel-1),经过曝光时间tobs,对npix个像素进行采样,暗电流噪声可表示为:Odc=Vndenpixtobs.2.2.4读出噪声CCD读出信号时将像素上的电荷转换为数值会产生读出噪声(readout noise).各个像素上的读出噪声不均匀,通常用电子数的均方根(root meansquare,R M S)来表征,设为nro(单位为e-RMS/pixel),则对npix个像素进行采样,读出nexp次,读出噪声可表示为:Oro=Vnpixnexpnro.2.2.5#抖动噪声根据2.1.4节所述仪器抖动的影响,抖动使得恒星在探测器上的位置发生偏移,而不同位置对应了不同的PSF
24、和PRF,因此会给测光带来额外的噪声,即抖动噪声.我们通过生成多组随机的抖动序列,对同样的目标源进行重复测量,给出CCD像素接收到的总电子数的序列N入,sArray,计算其均方差(standard deviation,ST D)作为tobs曝光时间、npix个像素内的抖动噪声:jt=STD(N,sArray),其中,N,sArray由(5)式多次计算给出.2.3信噪比计算模型根据2.1节和2.2 节所述,目标源信号由(5)式给出,噪声由(6)-(13)式给出.总噪声为:Otot=VoPs+0g+0%+0%+0元,其中,OPs、O b g、d e、O r o 和ojit分别为泊松噪声、天光背景噪
25、声、暗电流噪声、读出噪声和抖动噪声.信噪比计算为:NA.=SNR=Otot(11)(12)(13)(14)N入.s(15)45-564卷其中,各参数的含义与2.1和2.2 节中相同.上述计算得到的恒星测光信噪比可以用来判断在该模式下能否探测到系外行星类太阳恒星周围热木星的典型凌星深度约为1%量级,Ke-pler通常将判定行星候选体的凌星信号信噪比限制为SNR7.115,也有工作探讨了SNR6 的合理性16 这里我们选取了SNR6作为凌星信号探测的阈值.将热木星凌星深度的典型值和凌星信号信噪比的阈值相除,我们认为当恒星观测的噪声水平,即恒星测光信噪比6 0 0 时,有能力探测到系外类木行星.3子
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 CSST 多色 测光 时序 精度 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。