Constraints on Annihilating Dark Matter from Gamma-Ray Background-Galaxy Shape Correlations: Model-independent Null Results and Moderate Template-based Signals

技术摘要：通过伽马射线背景—星系形状相关性对湮灭暗物质的约束

问题与动机
宇宙暗物质（DM）的本质仍然是现代物理学中的主要谜团。弱相互作用大质量粒子（WIMPs）是领先的候选者，它们预言了暗物质湮灭产物与宇宙大尺度结构（LSS）之间存在空间互相关。虽然未解析的河外伽马射线背景（IGRB）预计包含来自暗物质湮灭的分量，但将其从天体物理源（例如活动星系核、恒星形成星系）中分离出来具有挑战性。先前的观测研究结果褒贬不一：早期的研究由于天空覆盖范围有限而报告了零检测；而最近一项利用暗能量调查（DES）第一年数据和 9 年 Fermi-LAT 数据的分析，通过匹配滤波器方法报告了一个信噪比为 8.9 的探测结果。然而，该探测依赖于特定的模型假设。本文重新审视了 IGRB 与星系形状（宇宙剪切）之间的互相关，旨在利用显著更大的数据集和傅里叶空间分析框架，提供更稳健、模型无关的暗物质湮灭约束。

方法论
作者分析了约 $12,000 \text{ deg}^2$ 的共同天空区域，结合了：

1. 伽马射线数据： 14 年的 Fermi 大面积望远镜（LAT）Pass8 观测数据（2008 年 8 月至 2022 年 8 月）。数据覆盖能量范围从 100 MeV 到 1 TeV，分为 40 个对数微分箱，随后被分组为 9 个宏观能量箱（0.63–1000 GeV）。通过从掩模高纬度区域（$|b| \ge 30^\circ$）减去银河系前景色和已解析源（使用 4FGL-DR4 目录）来提取未解析的 IGRB。
2. 星系形状数据： 来自 DES 第 3 年（Y3）和 DECam 全数据处处（DECADE）项目的弱引力透镜剪切星表。这些数据集包含约 1 亿和约 1.7 亿个星系。星表被分为四个层析红移箱（平均红移 $\bar{z} \approx 0.34, 0.53, 0.75, 0.95$）。

分析框架：

• 傅里叶空间估计器： 不同于之前使用实空间相关性的 DES Y3 分析，本研究采用了傅里叶空间估计器来计算 IGRB 强度与 E 模宇宙剪切之间的互功率谱（$C_{\gamma E}$）。这种方法允许对统计不确定性和系统误差进行受控的量化。
• 功率谱测量： 利用伪-$C_\ell$ 算法（NaMaster），作者计算了所有 9 个能量箱与 12 个形状星表（4 个红移箱 $\times$ 3 个调查窗口）组合的互功率谱，总计 108 个功率谱。分析限制在多极矩 $30 < \ell < 512$ 之间，以避免像素化伪影。
• 测试模型：
  • 模型无关模型： 通过将测量的功率谱与零假设（无相关性）进行比较进行的 $\chi^2$ 测试。
  • 暗物质湮灭模型： WIMP 湮灭为 $b\bar{b}$、$\tau^+\tau^-$、$\mu^+\mu^-$ 和 $W^+W^-$ 通道的理论预测。这些模型包括初级 $\gamma$ 射线以及来自逆康普顿（IC）散射上散射背景光子的次级辐射。考虑了子结构提升因子（$b_{sh}$）为 3（LOW）、30（MID）和 100（HIGH）的情况。
  • 现象学模板： 使用来自文献 [11] 的两种模型（幂律和对数抛物线参数化）来测试与先前探测声明的一致性。

关键结果

1. 模型无关测试中的零检测： 主要结果是对互相关信号的零检测。在所有调查区域（DES Y3, DECADE NGC, DECADE SGC）中，每个自由度的信噪比（$S/N$）均与零一致。堆叠的 B 模（指示系统误差）也与零一致，表明系统误差得到了控制。
2. 对湮灭截面的约束： 基于零检测结果，作者推导了在暗物质质量为 7 GeV 至 10 TeV 之间时，暗物质湮灭截面 $\langle \sigma v \rangle$ 的 95% 置信度上限。
   • 对于约 30 的中等子结构提升因子，该分析排除了质量为 2–3 TeV 的类 wino 暗物质的增强湮灭截面。
   • 对于较高的提升因子 $\sim 100$，约束显著增强，排除了 7–40 GeV 暗物质湮灭为 $b\bar{b}$ 或 $\tau^+\tau^-$ 时的典型热残余截面（$\langle \sigma v \rangle = 3 \times 10^{-26} \text{ cm}^3/\text{s}$）。
   • 包含次级 IC 辐射显著加强了这些约束，特别是对于轻子通道和较重的暗物质质量。
3. 衰变暗物质： 研究还推导了衰变情景下暗物质粒子寿命的 2$\sigma$ 下限，根据衰变通道不同，范围在 $\sim 10^{26}$ 至 $10^{27}$ 秒之间。
4. 基于模板的分析： 当应用文献 [11] 中的现象学模板时，数据表现出对幂律模型的偏好，置信水平为 $\sim 3.9\sigma$（$q = 0.62 \pm 0.16$），而对数抛物线模型则被排除。一个假设简单 $\ell^{-1}$ 标度的匹配滤波器分析给出的 $p$ 值约为 0.01（对应于总 $S/N \approx 7.6$）。然而，作者指出，推断出的显著性高度依赖于所假设的信号模型，且幂律拟合在 100 GeV 附近显示出与预期平均强度标度的偏差。

意义与主张
本文声称提供了迄今为止利用 $\gamma$ 射线与星系形状互相关对暗物质湮灭最严格的宇宙学约束，其优势在于利用了最大的天空覆盖范围（$\sim 12,000 \text{ deg}^2$）。其主要意义在于模型无关的零检测结果，这挑战了此前高度依赖特定晕模型模板的探测声明。作者证明，虽然基于模板的分析可以产生 $\sim 3\sigma$–$4\sigma$ 级别的信号，但当在没有先验假设的情况下针对零假设进行测试时，这些信号并不稳健。

研究得出结论，宇宙学互相关提供的约束与局部 $\gamma$ 射线探测（例如矮星系球状星团）相当且具有互补性。结果特别排除了在高提升因子情景下，类 wino 暗物质和热残余暗物质的某些参数空间。作者强调，模型无关的零结果与中等程度的基于模板的信号之间的张力，凸显了未来通过与多个大尺度结构示踪者的联合分析，通过其环境依赖性来分离天体物理学贡献与暗物质贡献的必要性。