Joint Multi-Period Fermi-LAT and LHAASO Constraints on Axion-Like Particles from Mrk 421 Using Profile Likelihood with Gaussian Copula Correlation

技术摘要：结合 Fermi-LAT 与 LHAASO 多周期观测对 Mrk 421 中类轴子粒子（ALPs）的联合约束

问题陈述
类轴子粒子（ALPs）是极具理论吸引力的标准模型扩展，它们可以在外部磁场中诱导光子-轴子振荡，从而可能改变来自遥远源的高能伽马射线能谱。尽管先前的研究已经对 ALP 参数进行了约束，但仍有显著的参数空间未被覆盖。此类搜索中的一个主要挑战是处理系统不确定性，特别是与消减伽马射线的河外背景光（EBL）相关的部分。传统的处理方法通常将 EBL 吸收和光子-ALP 振荡视为独立的乘法因子，忽略了两者之间非平凡的相互作用。此外，多历元分析经常假设观测周期在统计上是独立的，忽略了由共同的源变率、仪器系统误差和模型假设所导致的微弱残余相关性。这些简化处理可能会导致偏差的约束或人为收紧置信区间。此外，耀变体 Mrk 421 的 GeV（Fermi-LAT）与 TeV（LHAASO）辐射之间存在微弱的相关性，表明其物理起源可能不同，这增加了对宽带数据进行 ALP 搜索时的解释难度。

方法论
作者提出了一个综合框架，利用 2021 年 3 月至 2024 年 3 月期间来自 Fermi-LAT 和 LHAASO 的五个同步观测历元的数据，来约束 Mrk 421 的 ALP 参数（$m_a$ 和 $g_{a\gamma}$）。

1. 传播形式化： 不同于因子化近似，本研究采用了一个自洽的光子传播框架，其中 EBL 吸收通过吸收项（$i/2\lambda_\gamma$）直接纳入光子-ALP 混合哈密顿量中。EBL 强度通过一个归一化因子（$K_{EBL}$）进行处理，该因子作为一个具有高斯先验（$\mu=1.0, \sigma=0.2$）的琐碎参数（nuisance parameter）来应对不确定性。
2. 源建模： 为了解决 Mrk 421 中微弱的 GeV-TeV 相关性问题，研究测试了两种不同的喷流发射模型：
   • 一种是双区混合模型，包含一个紧凑的内区（爆发式 VHE 发射）和一个扩展的外区（静息发射）。
   • 另一种是单区强子模型，考虑到微弱的相关性和截然不同的变率幅度，该模型在物理上是合理的。
     两者均利用螺旋和缠绕磁场分量，并使光子穿过波源、星际介质（IGMF）和银河系磁场（GMF）。
3. 统计框架： 作者开发了一种联合剖面似然分析（joint profile-likelihood analysis）。为了解决将历元视为独立的局限性，他们引入了高斯 Copula 形式化，并使用了一个保守的相关系数（$\rho = 0.03$）。该方法通过变换后的 $\Delta\chi^2$ 剩余分布来模拟历元间的微弱相关性。分析在无 ALP 假设下，全局优化琐碎参数（固有光谱参数和 $K_{EBL}$），进而对 ALP 质量和耦合进行剖面分析。
4. 数据处理： 光子生存概率使用 gammaALPs 代码在四维网格（能量、ALP 参数、EBL 归一化）上预先计算，并在似然评估期间进行插值。固有光谱模型根据经修正后的赤池信息准则（AICc）按历元进行选择。

核心贡献

• 统一传播模型： 本文在传播方程中实现了光子-ALP 混合与 EBL 吸收的自洽处理，避免了乘法因子化带来的偏差。
• 相关多历元分析： 通过引入高斯 Copula 来处理观测周期之间的微弱相关性，为结合多历元数据集提供了一种统计一致的方法，防止了置信区间的人为收紧。
• 宽带约束： 研究利用 Fermi-LAT（GeV）和 LHAASO（多 TeV 至 PeV）独特的宽带覆盖范围，探测了广泛质量范围内的 ALP 参数（0.1–500 neV）。
• 稳健建模： 通过测试双区和单区强子情景，本工作考虑了 Mrk 421 复杂的发射物理，确保约束不会过度依赖于单一的发射假设。

结果
利用联合剖面似然框架，作者得出了光子-ALP 耦合（$g_{a\gamma}$）的 95% 置信水平（CL）上限：

• 低质量区域（$m_a \lesssim 1$ neV）： 约束非常严格，下限为 $5.0 \times 10^{-14}$ GeV$^{-1}$，上限为 $8.33 \times 10^{-14}$ GeV$^{-1}$。该区域对 CAST 约束具有敏感性。
• 中等及高质量区域（$0.1 - 500$ neV）：
  • 对于双区混合模型，保守上限为 $g_{a\gamma} < 6.50 \times 10^{-12}$ GeV$^{-1}$。
  • 对于单区强子模型，保守上限为 $g_{a\gamma} < 7.30 \times 10^{-12}$ GeV$^{-1}$。
• 行为特征： 排除轮廓呈现出特征性的“跳跃-平台”（jump-plateau）特征。在 $m_a \approx 1$ neV 附近发生剧烈转变，此时上限增加约一个数量级，随后在高质量端逐渐放宽。下限在整个质量范围内保持稳定。
• 稳健性： 使用替代相关系数（$\rho = 0, 0.15, 0.3$）进行的测试表明，约束的变化保持在不到一个数量级的范围内，证实了结果的稳定性。

意义
本文声称，这些结果显著受益于 LHAASO 的宽带极高能覆盖范围以及对 Mrk 421 的长期监测。所开发的框架提供了一个稳健且可扩展的流水线，适用于未来基于下一代伽马射线天文台的多信使 ALP 搜索。通过显式建模微弱相关性并自洽地处理 EBL 吸收，本研究提供了更可靠且更保守的 ALP 参数空间约束，解决了以往因子化方法中存在的系统偏差。这项工作证明，Mrk 421 中微弱的 GeV-TeV 相关性要求使用复杂的发射模型（双区或强子模型）来进行准确的 ALP 约束，并且所提出的统计方法论能有效处理相关的各种不确定性。