5颗导航用X射线脉冲星计时分析_闫林丽.pdf
《5颗导航用X射线脉冲星计时分析_闫林丽.pdf》由会员分享,可在线阅读,更多相关《5颗导航用X射线脉冲星计时分析_闫林丽.pdf(16页珍藏版)》请在咨信网上搜索。
1、526588-1航空学报Acta Aeronautica et Astronautica SinicaFeb.15 2023 Vol.44 No.3ISSN 1000-6893 CN 11-1929/V5颗导航用 X射线脉冲星计时分析闫林丽1,葛明玉2,*,庹攸隶2,周庆勇3,4,叶文韬2,郑世界2,韩大炜21安徽建筑大学 数理学院,合肥2306012中国科学院高能物理研究所 粒子天体物理重点实验室,北京1000493西安测绘研究所,西安7100544地理信息工程国家重点实验室,西安710054摘 要:脉冲星计时特性分析是开展脉冲星导航和脉冲星时研究的基础,为其提供了精确的脉冲星计时模型参数和
2、辐射特征。本文在系统总结脉冲星计时分析方法的基础上,选择了当前导航中应用最广泛的 5颗脉冲星,并利用国内外在轨 X 射线观测卫星的最新观测数据对其分析,包括:“慧眼”硬 X 射线调制解调望远镜(HXMT)针对 Crab脉冲星和 PSR B1509-58的观测,中子星内部结构探测器(NICER)针对 PSRs J1821-2452A、J1939+2134和 J0030+0451的观测。通过对“慧眼”HXMT 和 NICER 在 20172021年高精度计时观测数据的分析,获得了较长时间段内脉冲星最新状态的自转特性和物理信息,一方面给定了脉冲星最新的、覆盖时间长、自转参数精度较高的 X射线星历,证
3、明在 X射线波段也可独立给出较高精度的星历;另一方面建立了它们清晰显著的积分脉冲轮廓,其中在选定的能量段范围内 3颗毫秒脉冲星的轮廓是当前最精确的,可为脉冲星导航研究提供最新的标准模板。本文仅对脉冲星的 X射线数据进行了分析,未来利用国内外多个望远镜,开展多波段联合计时分析,将是脉冲星计时研究的重要方向。关键词:脉冲星;导航;X射线;计时分析;星历中图分类号:V11;P124.4;P145.6 文献标识码:A 文章编号:1000-6893(2023)03-526588-16脉冲星是一种高度磁化的、旋转的中子星或白矮星,具有非常稳定的周期,在地球上接收到此类天体的射电信号是脉冲的形式,因此被命名
4、为脉冲星1-2。自 1967 年首次被观测证认3,50多年的时间内已经发现 3 000多颗脉冲星。首颗脉冲星是在射电波段被发现的,随着空间天文的发展,对脉冲星的观测已经推广到红外、光学、X射线、射线波段。利用脉冲星精确的周期性可以进行引力波探测4-5、监测和修正原子时的稳定度6-7、脉冲星导航8-9等研究。脉冲星是宇宙空间天然的时钟,它们的自转建立在恒星级质量天体的物理过程基础上,能够提供一种基于遥远自然天体并持续数百年至数十亿年的时间频率,不易被干扰并且不受地面系统控制,可以作为独立的时间基准10。实现脉冲星时建立及服务的前提是充分了解脉冲星的时间特征,并长期监测脉冲星的时间特征变化。对脉冲
5、星时间特征的研究称为计时分析,最终目标http:/ 引用格式:闫林丽,葛明玉,庹攸隶,等.5 颗导航用 X 射线脉冲星计时分析 J.航空学报,2023,44(3):526588.YAN L L,GE M Y,TUO Y L,et al.Timing analysis for five navigation X-ray pulsars J.Acta Aeronautica et Astronautica Sinica,2023,44(3):526588(in Chinese).doi:10.7527/S1000-6893.2022.26588收稿日期:2021-10-29;退修日期:2021-1
6、1-22;录用日期:2022-02-13;网络出版时间:2022-03-1510:02网络出版地址:https:/ X射线甚至 射线波段,脉冲星的空间观测具有几乎不受星际介质的影响,不受到人为干扰的优势,但是大多数脉冲星的高能辐射信号强度弱,并且由于空间载荷重量的限制,有效探测面积比较小,需要较长时间观测才能获得比较准确的计时参数。X 射线、射线与射电波段联合观测,可充分利用各自优势,以实现高效的计时参数获取和探索不同观测资源的优化配置。脉冲星的工程应用研究之一 X 射线脉冲星导航为一种新型天文导航方式,旨在通过分析处理脉冲星的 X 射线辐射信号,获取航天器的位置、姿态、时间等完整的导航信息,
7、具有自主性强、抗干扰能力强、可靠性高等优点11。X 射线脉冲星的空间位置、辐射特性及自转周期等信息是 X 射线脉冲星导航的基本输入参数,是实现脉冲星导航的前提条件,需要精确测定这些参数。由于大多数脉冲星的 X 射线辐射流量随能量呈幂律分布,通常选择辐射较强的软 X射线波段作为观测窗口,选择能够同时产生射电和 X 射线信号的脉冲星作为导航脉冲星,通过地面射电天文台和空间 X 射线卫星联合观测,以获取脉冲星的空间位置、自转周期、X 射线脉冲波形和相位等信息,这是当前的最佳脉冲星计时方案。1999 年,ARGOS 卫星携带的非常规恒星定位 实 验 仪 器(Unconventional Stellar
8、 Aspect)进行了首次 X 射线空间导航尝试12。2017 年 11月,NASA 在其官网发布了“空间站 X 射线计时与导航技术”试验项目的结果。该项目是国际上首次开展基于毫秒脉冲星的空间自主导航试验,利用搭载在国际空间站上的中子星内部结构探测 器(Neutron star Interior Composition Explorer,NICER),对 4 颗毫秒脉冲星进行了为期2 d 的观测,获得了 78 组观测数据,并利用 8 h 的实验数据,实现了空间站位置 16 km 误差内的导航,其中数据较好一段的可优于 5 km 的导航精度13。可见,X 射线脉冲星导航具有完全自主导航的潜力,并
9、期望其在深空探测领域的应用将逐步变得成熟,为太阳系内乃至太阳外的探测器提供自主的导航定位服务。中国科学院高能物理研究所和空间应用中心的研究人员利用中国载人航天天宫二号空间实验室上的 射线暴偏振探测器(-ray Burst Polarimeter,POLAR)对蟹状星云脉冲星进行了 1 个月的探测,完成了中国首次脉冲星导航空间实验,实现了对天宫二号的定轨,定轨方法采用脉冲轮廓显著性与卫星轨道的关联分析,实现的定轨精度如下:轨道位置精度在X、Y、Z方向分别为 11.5、13.1和 9.1 km,对应的轨道速度精度分别为 0.3、0.31、0.01 km/s 14。2019年,该研究团队利用“慧眼”
10、硬 X 射线调制望远 镜(Hard X-ray Modulation Telescope,HXMT)采用改进的定轨方法,进一步验证了利用脉冲星定轨的可行性,精度可达 10 km9。还需提到的是,XPNAV-1 也进行导航实验,利用Crab 脉冲星的观测,其轨道的确定精度为 38.4 km15。西安测绘研究所的研究人员对脉冲星时空基准建设也开展了系列研究16-19。目前脉冲星计时分析方法的介绍集中在射电波段,而在 X 射线波段研究相对较少。X 射线与射电波段的观测手段不同,记录的数据信息不同,因此计时过程存在一定的差异。在射电波段,地面望远镜接收脉冲星的辐射信号,首先需要消色散,即射电脉冲通过星
11、际介质后会产生延迟现象,然后再进行计时分析。而在 X 射线波段,则没有色散延迟。但在分析数据时,首先需要先进行太阳系质心修正,再累积脉冲轮廓进行计时分析。因此,每个波段的观测,进行完各自的修正之后,后续的计时分析过程对于射电与 X射线观测无差异。本文通过分析处理“慧眼”HXMT 和 NICER对 5 颗脉冲星的观测数据,说明脉冲星 X 射线计时分析的步骤和方法,并给出它们计时模型,以便于应用于脉冲星导航时间和空间精度的分析研究工作中。航空学报526588-31卫星介绍及脉冲星的观测信息1.1“慧眼”HXMT卫星2017-06-15,中国第一颗空间 X 射线天文卫星“慧眼”HXMT 在酒泉卫星发
12、射中心成功发射,运行于高度 550 km、倾角 43的近地圆轨道上。“慧眼”HXMT 卫星携带了高能 X 射线望远镜(High Energy X-ray telescope,HE)、中能 X 射线望远镜(Medium Energy X-ray telescope,ME)和低 能 X 射 线 望 远 镜(Low Energy X-ray telescope,LE)3 种科学载荷和空间环境监测器。慧眼 HXMT 卫星对脉冲星采取定点观测,其优点包括覆盖能段宽、有效面积大、时间分辨率高、探测死时间很小、对强源观测没有光子堆积效应,因此该卫星具有独特的研究天体多波段快速光变的能力。“慧眼”HXMT 卫
13、星的基本指标如表 1所示。1.2NICER2017-06-03,NICER 被搭载在国际空间站上,旨在测量来自脉冲星的 X 射线脉冲轮廓,以便更好地研究中子星状态方程。与 X 射线计时探测器(Rossi X-ray Timing Explorer,RXTE)相比,NICER 在软 X 能段上工作,并有更高的能量分辨率、计时分辨率和灵敏度。NICER主要性能有:大 有 效 面 积:约 1 900 cm2 1.5 keV;能量范围:0.2 keVE12.0 keV;绝对计时精度:300 ns;能谱分辨率:0.58 keV时6 E/E100 s 的总观测次数。NICER 也对Crab脉冲星进行了大量
14、观测,更详细 X 射线计时分析工作见已公开的论文25。2X射线计时分析流程2.1数据筛选和预处理“慧眼”HXMT 和 NICER观测得到的数据传输至地面系统后,经过解算和原始筛选,观测信息被保存成 FITS 文件以供工程或科学分析。科学工作者将获得的观测数据进一步筛选以尽量减少干扰源的影响,例如剔除卫星指向偏离较大、被地球遮挡、经过大西洋异常区、空间环境污染比较严重的时间段内的光子。数据筛选后,光子到达探测器的时间、光子能量以及卫星的轨道信息是进行 X 射线计时分析主要采用的。在利用这些信息进行计时分析之前,需要对光子到达探测器的时间进行修正。FITS 数据文件记录的光子时间是 X 射线光子到
15、达探测器的时间,由于周年视差、Roamer延迟、引力时延等效应2,在进行计时分析之前需计算这些效应引起的时间偏离量,并进行时间修正,即将光子到达探测器的时间转换至光子到达太阳系质心位置的时间。表 35颗脉冲星的观测信息Table 3Basic observational information of five pulsars名称CrabPSR B1509-58PSR B1509-58PSR J0030+0451PSR J1821-2452APSR J1939+2134观测器HXMTHXMTNICERNICERNICERNICER观测时间范围2017-08-272021-01-142017-07
16、-062018-09-132017-06-292017-08-112017-07-242020-12-312017-06-252020-12-312017-06-242020-12-31观测次数3404813495358500曝光时间/(103s)604(HE)769(ME)94(HE)145(ME)163 1529891 786表 25颗脉冲星的基本参数Table 2Basic parameters of five pulsars名称CrabPSR B1509-58PSR J0030+0451PSR J1821-2452APSR J1939+2134赤经05:34:31.94015:13:5
17、5.81100:30:27.42818:24:32.00819:39:38.561赤纬+22:00:52.20059:08:09.600+04:51:39.71124:52:10.834+21:34:59.126自转周期P0/ms29.947151.8404.8653.0541.558周期导数 P14.21010131.52910121.01710201.61910181.0511019距离/kpc2.0004.4000.3255.5003.500表面磁场/G3.7910121.5410132.251082.251094.09108自转能损率/(ergs s1)4.510381.710373.
18、510332.210361.11036特征年龄/a1.261031.571037.581092.991072.35108注:1 kpc=3.081019 m;1 G=104 T;1 ergs s1=107 J s1。航空学报526588-5需要注意的是,即便进行了数据筛选,筛选后的数据除了来自脉冲星辐射的光子外,还包括视场内脉冲星星云的辐射、空间弥散本底、探测器的噪声等“无关光子”形成的背景噪声。尽管好时间段选择可以有效降低背景噪声的影响,但并不能完全剔除背景光子,这也将成为计时分析的主要干扰。当脉冲星的脉冲辐射远强于背景时,“无关光子”影响可以忽略;当脉冲星本身的辐射比较弱时,应采取一定措施
19、消除干扰或者增加曝光时间,利用信噪比高的观测,提高计时结果的准确性。其中,选取辐射光子流量比较强的能量段也是排除干扰、提高脉冲信噪比的有效方法之一。在本文中,使用“慧眼”HXMT 的数据处理软件 HXMTDAS v2.0426对慧眼数据进行数据预处理。“慧眼”HXMT 的数据筛选过程,包含对温度增益的修正、事例的分级和重建,以及好时间段(Good Time Interval,GTI)的判断预筛选过程。使用软件中的批处理指令 hpipeline27,以及推荐的 GTI筛选条件28,最后产生经过质心修正以及筛选的事例文件以供后续计时分析使用。使 用 HXMTDAS 中 的 指 令 hxbary 对
20、“慧 眼”HXMT 的观测数据进行太阳系质心修正。对于“慧眼”HXMT 的载荷 HE,筛选了 27250 keV的观测能段,对于 ME 选择了 1035 keV 的观测数据。对于 LE 的观测,由于电子学读出有约850 s 的延迟,并会使脉冲轮廓的形状产生形变29-30,因此在本文中未使用 LE的数据。针对 NICER 观测数据的筛选可通过天文数据 处 理 软 件 HEASOFT v.6.27 中 包 含 的NICERRDAS v.5 软 件 包 完 成。与 上 述 慧 眼HXMT 的数据筛选过程类似,NICER 观测数据的筛选先后使用指令 nicerl2、niprefilter2、nimak
21、etime 以 及 nicermergeclean 实 现,最 后 用 指 令barycorr对光子到达时间进行太阳系质心修正即可进行后续的计时分析。详细的数据筛选指令可 参 考 HEASOFT 网 站 数 据 处 理 指 导 手 册 及NICER相关工作31-32。2.2X射线计时分析步骤如图 1 所示,X 射线计时分析分为以下几个步骤:周期搜寻,搜索短时间观测内脉冲星的自转频率,并折叠该时间段内的脉冲轮廓,即光子在相位范围 0,1 内的统计分布曲线;到达时间计算,计算累积脉冲轮廓的到达时间(Time of Arrival,TOA);TOA 序列生成,对不同时间段的观测重复前两步,形成脉冲星
22、 TOA 序列文件;计时分析,利用 Tempo233或其它计时分析工具拟合计时模型参数,更新脉冲星星历。在接下来的 2.2.12.2.3节中,详细描述 X 射线计时分析的步骤。2.2.1自转频率搜寻与轮廓折叠脉冲星自转的过程中向外发射 X 射线光子,到达 X 射线探测器后被记录下来,这些光子来自脉冲星不同的自转周期内,但可以根据自转频率计算出每个光子在一个周期内所处的相位,统计各个相位区间光子的分布,即可得到脉冲星在一个周期内的光子数分布脉冲轮廓曲线。设t0为参考时间零点,ti为第 i个光子到达太阳系质心的时间,f0、f1、f2分别为脉冲星在参考时间零点对应的自转频率及其一阶、二阶导数,第 i
23、个光子相对于t0的绝对相位为i=f0(ti-t0)+f1(ti-t0)22!+f2(ti-t0)33!(1)该光子在一个周期内的相位位置只需要将i对 1图 1X射线计时分析流程Fig.1X-ray timing analysis processes航空学报526588-6取余即可,因此得到脉冲轮廓的关键参数就是脉冲星的自转频率及其各阶导数。NICER 对 3 颗毫秒脉冲星的单次观测时间约几千秒,可不考虑式(1)第 2项以后各部分引起的相位偏差,在进行自转频率搜寻时只考虑f0即可。自转频率搜寻的方法是:对给定的一段光子到达时间序列,选定合适的频率搜索范围和步长,计算各个自转频率对应的统计量,根据
24、统计量与自转频率的对应关系曲线,确定最佳的自转频率参数。不同的统计量选取方式,对应不同的自转频率搜寻方法,常用的方法有 Pearson 2检验34-35、Z2m检验36和 H检验37。例如,对于第 1种方法,假设相位 0,1 被平均分为 m 个子相位区间,这一序列的总光子数为 N,第 i个子相位区间的光子数为ni个,则对应的 Pearson 2量为2=i=1m(ni-N/m)2N/m(2)当2量最大时,对应的自转频率为脉冲星产生的这一序列光子的最佳自转频率。对不同时间段的光子序列进行自转频率搜寻,可以观察到脉冲星的自转频率随时间演化的大致趋势。图 2列举出毫秒脉冲星 PSR J1939+213
25、4某单次观测自转频率搜寻结果及折叠出的脉冲轮廓,这次观测的周期搜寻结果比较理想,即自转频率与统计量的变化曲线在最佳自转频率附近比较对称,且折叠出的累计脉冲轮廓信噪比较高。然而,由于本底光子的影响,仍有很多观测不能够搜索出有效的脉冲信号,那么数据处理过程中需剔除这部分观测。2.2.2TOA的计算利用 Tempo2软件进行计时分析的对象就是脉冲到达时间,可由单次观测的脉冲轮廓计算得到。将观测轮廓与标准轮廓进行交叉相关,拟合出观测轮廓与标准轮廓之间的时间间隔,再加上标准轮廓的到达时间,即为观测轮廓的 TOA。当实际观测的起始时间为t0,以该时刻为参考的积分脉冲轮廓的最高峰对应相位为0,搜索到的最佳自
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 导航 射线 脉冲 计时 分析 闫林丽
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。