唐庆文

(南昌大学物理系，江西 南昌　330031)

近几年来，天体的高能伽玛射线波段巡天成为天文学研究的重要手段，究其原因在于伽玛射线可以反应各类高能物理过程，如宇宙线粒子加速和辐射过程、宇宙暗物质作用过程等[1]。研究天体源的伽玛射线能谱特征，可以证论某些高能物理过程是否主导，如中性π介子衰变可以形成峰值流量在67.5 MeV的鼓包能谱特征[2]。一方面，伽玛射线巡天可以有目标地观测已在其他波段观测过的天体源，如哈勃望远镜(Hubble Space Telescope)的光学/红外巡天目录，斯隆数字巡天(Sloan Digital Sky Survey, SDSS)目录以及澳大利亚帕克斯天文台64 m射电望远镜(Parkes Radio Telescope)巡天目录[3]，观测的数据可进一步限制目标源的物理、化学特性。另一方面，高能伽玛射线探测器会被一些强烈的天体伽玛暂现源触发，例如，伽玛射线暴在GeV能区的存在时标通常大于keV能区的时标(通常不超过100 s)，在暴触发后的数千秒也可能被观测到[4-5]，因此探测器可应用于探究该类天体伽玛射线暂现源的辐射过程。随着我国天文学观测的发展，特别是暗物质卫星(DArk Matter Particle Explorer, DAMPE，悟空号)、硬X射线调制望远镜(Hard X-ray Modulation Telescope, Insight-HXMT，慧眼号)的升空[6-7]，它们的科学目标集中在搜寻一些未知源的多波段辐射并进而解释一些极高能的物理过程，这些过程一般在实验室里很难实现。因此，快速定位未知伽玛射线源是实时跟进未知源的有效手段。

由于在高能伽玛波段探测器的角分辨率较低，研究掌握精准定位伽玛射线源的算法，对于伽玛射线巡天观测起着很重要的作用。实验研究采用国际成熟的费米卫星携带的大面积望远镜(Large Area Telescope, LAT)的伽玛射线巡天观测为实验对象，选取伽玛射线暴为暂现源样本，采用(1)不同时间、(2)不同中心位置计算伽玛射线暂现源的检验统计(Test Statistic, TS)位置图，最终得到最优化的定位。其次将得到的位置(含误差)与其他波段的后随观测(Follow-up Observations)比照，调整实验设定参数，从而得到实验结果与实现研究目标。

1　费米卫星巡天观测及实验样本选择

费米伽玛射线太空望远镜(原名GLAST，现命名为Fermi Gamma-Ray Space Telescope)是多国合作研发的，包括美国、法国、意大利、日本、瑞典等国家。它的科学目标主要针对活动星系核、脉冲星、伽玛射线暴、宇宙线等。搭载的望远镜视野达到全天区的20%，于2008年8月4日开始科学观测，并于一年后实时公布观测数据。LAT伽玛波段的观测能量范围从30 MeV到大于300 GeV，单个入射的天体伽玛光子被各层探测装置逐级记录，Fermi-LAT研究组系统地分析了这些光子事例并进行分类，如最常用的Source和Clean类型光子。官方统计从2008年8月4日到2011年8月4日，LAT共有超过1.8 × 1011个事件，最终有约1.44 × 108个通过更细致的伽玛光子选择标准并发布。图1(a)是Fermi-LAT 5年观测的大于1 GeV伽玛射线分布图。其中，银盘区域是最亮的伽玛射线探测区域，该区域之外存在一些亮点天体源，例如大麦哲伦星云(LMC)、蟹状星云(Crab)。

实验样本为LAT观测的122个伽玛射线暴，时间范围为2008年8月4日至2016年10月30日。图1(b)是GRB080916009(数字为Fermi触发暴的命名方式，前6位为年月日，后3位区分当天观测到的不同伽玛射线暴)的100 MeV～200 GeV光子计数率的时间变化图。从图中可以看出，在触发(Trigger Time，T0)后100 s以内，光子计数率呈现脉冲式分布，而T0之前及T0+100 s后，光子计数率比较平稳。实验样本的每个伽玛暴，选择T0到T0+1 000 s的100 MeV～100 GeV的光子事例。样本选择标准为在触发后1 000 s内有高能光子被探测到，即能量大于100 MeV的光子。通过初步筛选发现,16个伽玛暴，或者只有10～100 MeV的低能伽玛光子探测，或者在较晚期(T0+1 000 s之后)才有大于100 MeV光子探测，因此排除了这些伽玛暴样本，共收集到106个高能伽玛探测的伽玛暴。

图1(a) Fermi-LAT 5年观测到的大于1 GeV伽玛射线分布图,中心较亮的区域是银盘区域； (b) GRB080916009的高能光子(100 MeV～200 GeV)计数率随时间分布图

Fig.1(a) Galactic Gamma-ray distribution observed by the Fermi-LAT in 5-yr mission, the bright part of which is galactic plane; (b) Count rates of high-energy photons as a function of time since GBM trigger (100MeV-200GeV)

2　Fermi科学软件介绍和实验算法设计

实验是基于Fermi官方科学软件的最大似然拟合分析。首先，计算伽玛暴高能伽玛探测的检验统计值(TS,TS= σ2，当达到5 σ 为99.999%置信度)随不同起止时间和不同中心位置选择的分布[8]。单次数据处理流程包括：(1)光子事件类型、中心位置及半径(12°)、观测时间、观测能段的选择(gtselect工具)；(2)飞船位置相关选择(gtmktime工具);(3)弥散伽玛背景计算(gtdiffrsp工具)；(4)光子曝光计算(gtexpmap工具)；(5)最大似然拟合(gtlike工具)；(6)TS位置图计算(gttsmap工具)。算法设计就两个条件做了阈值设定，即：

(1)时间间隔选择。对单个伽玛射线暴，所有的时间起点(T1)都是伽玛暴的触发时刻(T1=T0)，而时间终点(T2)依次累加(δt)，如0 s

(2)中心位置。该步骤中有两个位置。(a)初始位置(P0)，以Fermi-LAT的初步结果为初始位置(P0)；(b)以单次数据处理流程得到TS位置图的最大TS值的位置(P1)，再次执行最大似然拟合,得到本次数据处理流程的最大TS值。

3　实验结果及讨论

图2两个伽玛暴的r、R95、RTS,max、RN,pred随不同T2的变化图
Fig.2Four observables as a function of selected end Time (T2)

位的P3的95%等高误差线以内，而GRB120709883的X射线波段的后随探测位置则稍微偏离了P3的95%等高线内部，但此时R95或r都小于1°，考虑LAT高能伽玛射线在1 GeV处95%的角分辨率(稍大于1°)，P3是有效的定位。

探测类型GRB数目/个探测类型GRB数目/个Strong (强)90Low (低)1Middle (中等)14None (无)1

4　实验总结

本文介绍了高能伽玛射线暂现源的位置探测的一种算法，主要通过探测置信度随时间变化，考察了不同截止时间对探测置信度的影响。研究结果表明，截止时间通常位于T0到T0+1 000 s之间，这通常可能由于飞船运行状态、LAT中心轴位置变化或者本征的伽玛射线光子分布造成的。结果表明实验方案可得到99.1%的伽玛暴样本的明显探测和84.9%的样本的强伽玛探测。综上所述，实验方案可快速地对观测数据进行伽玛射线源的位置确定，该位置可通过伽玛射线暂现源定位网络(GCN)、行星际网络(IPN)等发布信息给其它空间或地面望远镜、探测器，进行有效的跟踪观测，从而对未知天体伽玛射线源进行多波段研究。

实验算法可作为实验室开放性的学科实验，对计算机性能要求适中，适合具备天文学科或高能物理等硕士点以上的院校进行教学和科学研究之用。在各个天文机构的服务器中，教师和学生可将本算法作为一个代码包应用于其他天体或天体系统的精确定位[9-10]。