[论文解读] The Denoised, Deconvolved, and Decomposed Fermi $γ$-ray sky - An application of the D$^3$PO algorithm

Marco Selig, V. Vacca+2 more|Oct 16, 2014

Astrophysics and Cosmic Phenomena参考文献 48被引用 30

一句话总结

该论文将D³PO贝叶斯推断算法应用于6.5年的费米LAT数据（0.6–307.2 GeV），在不依赖光谱或空间模板的情况下，对弥漫和点状伽马射线天空进行去噪、去卷积和分解。该方法揭示了两组弥漫发射结构——由宇宙射线与冷稠密ISM相互作用产生的强子成分（软，冷ISM）和与费米气泡相关的轻子成分（硬，费米气泡）——其大尺度角功率谱指数为2.47±0.02，并生成了一个包含3,106个源候选者的星表，其中包括在具有射电晕的星系团方向检测到的伽马射线发射。

ABSTRACT

We analyze the 6.5yr all-sky data from the Fermi LAT restricted to gamma-ray photons with energies between 0.6-307.2GeV. Raw count maps show a superposition of diffuse and point-like emission structures and are subject to shot noise and instrumental artifacts. Using the D3PO inference algorithm, we model the observed photon counts as the sum of a diffuse and a point-like photon flux, convolved with the instrumental beam and subject to Poissonian shot noise. D3PO performs a Bayesian inference in this setting without the use of spatial or spectral templates;i.e., it removes the shot noise, deconvolves the instrumental response, and yields estimates for the two flux components separately. The non-parametric reconstruction uncovers the morphology of the diffuse photon flux up to several hundred GeV. We present an all-sky spectral index map for the diffuse component. We show that the diffuse gamma-ray flux can be described phenomenologically by only two distinct components: a soft component, presumably dominated by hadronic processes, tracing the dense, cold interstellar medium and a hard component, presumably dominated by leptonic interactions, following the hot and dilute medium and outflows such as the Fermi bubbles. A comparison of the soft component with the Galactic dust emission indicates that the dust-to-soft-gamma ratio in the interstellar medium decreases with latitude. The spectrally hard component exists in a thick Galactic disk and tends to flow out of the Galaxy at some locations. Furthermore, we find the angular power spectrum of the diffuse flux to roughly follow a power law with an index of 2.47 on large scales, independent of energy. Our first catalog of source candidates includes 3106 candidates of which we associate 1381(1897) with known sources from the 2nd(3rd) Fermi catalog. We observe gamma-ray emission in the direction of a few galaxy clusters hosting radio halos.

研究动机与目标

• 在不依赖预定义光谱或空间模板的前提下，从费米LAT数据中重建弥漫和点状伽马射线天空。
• 解析全天弥漫伽马射线发射在形态和光谱上的复杂性。
• 利用非参数贝叶斯推断方法，分离弥漫通量中的强子与轻子贡献。
• 生成包含3,106个候选源的源星表，并与已知的费米源星表进行验证。
• 研究弥漫伽马射线发射的角功率谱及其能量无关性。

提出的方法

• D³PO算法对费米LAT光子计数图进行贝叶斯推断，将观测计数建模为弥漫和点状通量分量与仪器响应函数卷积及泊松噪声的组合。
• 该方法使用非参数、基于稀疏性的先验，实现去噪、去卷积和分解，而无需假设参数化模板。
• 该算法同时估计弥漫和点源通量分量，及其不确定性，以及弥漫发射的角功率谱。
• 重建在HEALPix球谐框架下使用NIFTy软件包完成。
• 该方法通过伪彩色可视化和角功率谱分析，实现光谱指数制图及形态学分解，将发射分为软（强子）和硬（轻子）成分。
• 源候选者从去卷积后的点源图中识别，并与第二版和第三版费米LAT源星表进行交叉匹配。

实验结果

研究问题

• RQ1 能否在不依赖预定义光谱或空间模板的前提下重建弥漫伽马射线天空？
• RQ2 弥漫伽马射线发射在全天范围内的形态和光谱结构是什么？
• RQ3 弥漫发射是否与强子和轻子过程的两组分模型一致？
• RQ4 弥漫伽马射线通量的角功率谱是什么？其是否随能量变化？
• RQ5 是否存在新的伽马射线源候选者？它们能否与已知天体（如星系团）关联？

主要发现

• 弥漫伽马射线天空最适宜由两个独立成分描述：一个软的强子成分，对应于冷而致密的星际介质；一个硬的轻子成分，与费米气泡及喷流相关。
• 弥漫通量的角功率谱在所有能量 band 中均遵循幂律，指数为2.47±0.02，且与能量无关。
• 星际介质中尘埃与软伽马射线的比值随银纬升高而降低，表明强子发射效率具有空间变化性。
• D³PO算法成功重建了高达数百GeV的弥漫通量形态，实现了首次全天弥漫伽马射线背景的光谱指数图。
• 首个D³PO费米源星表包含3,106个候选源，其中1,381个（1,897个）与第二版（第三版）费米LAT星表中的源匹配。
• 在若干已知具有射电晕的星系团方向检测到伽马射线发射，提示其可能与相对论性电子的逆康普顿发射有关。