Gamma-Ray Constraints on Heavy Axion-Like-Particle Decays from Fermi-LAT and H.E.S.S. Blazar Spectra

Resumo Técnico: Restrições de Raios Gama sobre Decaimentos de Axions de Massa Pesada a partir dos Espectros de Blazares do Fermi-LAT e H.E.S.S.

Enunciado do Problema
Raios gama de energia muito alta (VHE; $E_\gamma \ge 100$ GeV) de fontes extragalácticas, como blazares, são atenuados durante a propagação pelo espaço intergaláctico devido à produção de pares ($\gamma\gamma \to e^+e^-$) com os fótons do Fundo de Luz Extragaláctica (EBL). Embora os modelos padrão de EBL baseados na evolução de galáxias e levantamentos profundos descrevam com sucesso essa atenuação, potenciais anomalias na transparência do Universo motivaram a exploração de física Além do Modelo Padrão (BSM). Especificamente, este estudo investiga um cenário envolvendo Axions de Massa Pesada (ALPs) com massas $m_a \sim 2,5\text{--}20$ eV. Diferente de ALPs mais leves, onde as oscilações fóton-ALP poderiam aumentar a transparência, ALPs pesadas nesta faixa de massa são hipotetizadas para decair em dois fótons ($a \to \gamma\gamma$) ao longo de escalas de tempo cosmológicas. Se os ALPs constituírem uma fração significativa da matéria escura, seu decaimento injetaria fótons adicionais no fundo cósmico infravermelho/óptico, aumentando, portanto, a intensidade do EBL e aumentando a profundidade óptica de raios gama. Este artigo busca restringir a massa ($m_a$) e o acoplamento de fóton ($g_{a\gamma}$) de tais ALPs pesadas, testando se os espectros observados de blazares são consistentes com um modelo de EBL aumentado pelas contribuições de decaimento de ALPs.

Metodologia
A análise utiliza um conjunto de dados conjunto composto por observações do Fermi Large Area Telescope (Fermi-LAT) e do High Energy Stereoscopic System (H.E.S.S.).

• Seleção de Dados: A amostra consiste em 11 blazares selecionados com base em sua detecção em energias VHE e redshifts onde a absorção pelo EBL é significativa ($\tau > 2$). O conjunto de dados inclui 14 anos de dados do Fermi-LAT (2008–2022) e observações do H.E.S.S.
• Segmentação Temporal: Para lidar com a variabilidade da fonte, os autores empregaram o algoritmo de blocos Bayesianos para segmentar as curvas de luz. A análise espectral foi realizada nesses segmentos para garantir que a evolução espectral intrínseca não enviesasse os resultados. Para fontes variáveis, os dados contemporâneos do Fermi-LAT e H.E.S.S. foram pareados; para segmentos pré-Fermi, apenas os dados do H.E.S.S. foram utilizados.
• Modelagem Espectral: Uma análise conjunta de máxima verossimilhança foi realizada utilizando o pacote Gammapy. Os espectros intrínsecos das fontes foram modelados usando funções de lei de potência ou log-parábola. A análise incluiu um parâmetro de perturbação para contabilizar um deslocamento coerente na escala de energia nos dados do H.E.S.S. ($\sim 10\text{--}15\%$ de incerteza).
• Formalismo de Decaimento de ALP: A contribuição dos ALPs em decaimento para o EBL foi modelada usando o formalismo de [54–56]. A intensidade específica do EBL induzido por ALP depende da massa do ALP, do acoplamento e da densidade ($\Omega_a$). A profundidade óptica total foi calculada como a soma da contribuição padrão do EBL (usando o modelo de Domínguez et al. [9] como base, com comparações com outros modelos) e a atenuação adicional induzida pelos fótons de decaimento de ALP.
• Estrutura Estatística: Os autores construíram uma razão de verossimilhança perfilada, $\lambda(m_a, g_{a\gamma})$, somando as estatísticas-C em todas as fontes e segmentos de tempo. Limites de exclusão foram derivados no nível de confiança de 95% (C.L.) sob a suposição de que os ALPs constituem toda a densidade de matéria escura ($\Omega_a = \Omega_{DM}$).

Principais Contribuições

• Primeira Análise Conjunta: Este trabalho representa a primeira análise de assinaturas de decaimento de ALPs pesadas usando um framework de verossimilhança unificado combinando dados do Fermi-LAT e H.E.S.S., maximizando a robustez estatística através de uma ampla faixa de energia (MeV a TeV).
• Cobertura do Espaço de Parâmetros: O estudo varre a faixa de massa de ALP de $m_a \approx 2,5\text{--}20$ eV. Esta faixa é fisicamente motivada porque os fótons de decaimento nesta janela de massa caem na banda óptica ao infravermelho próximo, o que afeta mais fortemente a atenuação de raios gama TeV via produção de pares.
• Tratamento de Sistemas: A análise incorpora explicitamente as incertezas de escala de energia do H.E.S.S. e testa a robustez dos resultados contra cinco modelos diferentes de EBL (Domínguez, Franceschini, Finke, Saldana-Lopez, e Porras-Bedmar & Meyer).

Resultados

• Limites de Exclusão: A análise combinada produz limites de exclusão de 95% C.L. para o acoplamento fóton-ALP $g_{a\gamma}$. As restrições mais fortes são encontradas em massas mais altas, onde o tempo de vida é mais curto, levando a uma maior emissividade de fótons.
  • Em $m_a \approx 15$ eV, o limite atinge $g_{a\gamma} \sim 7 \times 10^{-12}$ GeV$^{-1}$.
  • Em $m_a \approx 20$ eV, o limite é aproximadamente $g_{a\gamma} \sim 5,4 \times 10^{-12}$ GeV$^{-1}$.
• Consistência com a Hipótese Nula: Os dados não mostram preferência estatística significativa por um componente de ALP sobre o modelo padrão de EBL. A vasta maioria dos segmentos individuais de fontes e o conjunto de dados combinado são consistentes com a hipótesa nula (sem decaimento de ALP).
• Robustez do Modelo: As restrições derivadas são robustas contra incertezas na modelagem do EBL, com variações entre diferentes templates de EBL permanecendo dentro de um fator de aproximadamente 1,5.
• Comparação com Alegações Anteriores: Os resultados desfavorecem o espaço de parâmetros sugerido em trabalhos anteriores (ex: [55, 67]) que propunham que decaimentos de ALP poderiam explicar excessos nos espectros de fundo óptico cósmico ou de EBL. Especificamente, a região sugerida por [55] ($m_a \sim 8\text{--}20$ eV, $g_{a\gamma} \sim (3\text{--}6) \times 10^{-11}$ GeV$^{-1}$) é excluída sob a suposição de $\Omega_a = \Omega_{DM}$.

Significância e Alegações
O artigo afirma fornecer as restrições mais fortes até o momento sobre ALPs pesadas na faixa de $m_a \sim 15\text{--}20$ eV baseadas em observações de raios gama. Ao fechar uma região de espaço de parâmetros anteriormente não restringida, estes resultados oferecem um teste complementar às buscas espectroscópicas por linhas de decaimento estreitas e experimentos de laboratório. Os autores enfatizam que seu método depende da contribuição integrada de decaimento para o EBL e da modificação resultante da profundidade óptica de raios gama, tornando-o independente de suposições sobre fatores-D específicos de fontes ou fundos de linhas estreitas. Embora a análise não encontre evidência de decaimento de ALP pesada, ela estabelece limites superiores rigorosos que restringem a viabilidade de tais partículas como um componente dominante de matéria escura capaz de alterar significativamente a transparência de raios gama do Universo. O trabalho demonstra a eficácia de análises de raios gama multi-instrumento para sondar a física BSM e estabelece uma linha de base para futuras buscas mais sensíveis com instrumentos de próxima geração, como o Cherenkov Telescope Array (CTAO).