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

Resumo Técnico: Restrições Conjuntas de Multiperíodo Fermi-LAT e LHAASO sobre Axions-Like Particles de Mrk 421

Problema
Partículas do tipo axion (ALPs) são extensões bem motivadas do Modelo Padrão que podem induzir oscilações fóton-ALP em campos magnéticos externos, potencialmente alterando os espectros de raios gama de fontes distantes. Embora estudos anteriores tenham restringido parâmetros de ALPs, regiões significativas do espaço de parâmetros permanecem não restringidas. Um grande desafio nessas buscas é o tratamento de incertezas sistemáticas, particularmente relacionadas à Radiação de Fundo Extragaláctica (EBL), que atenua os raios gama. Abordagens tradicionais frequentemente tratam a absorção pela EBL e as oscilações fóton-ALP como fatores independentes e multiplicativos, negligenciando sua interação não trivial. Além disso, análises de múltiplos épocas frequentemente assumem que períodos observacionais são estatisticamente independentes, ignorando correlações residuais fracas decorrentes de variabilidade compartilhada da fonte, sistemáticas instrumentais e suposições de modelo. Essas simplificações podem levar a restrições enviesadas ou intervalos de confiança artificialmente estreitos. Adicionalmente, o blazar Mrk 421 exibe uma correlação fraca entre suas emissões de GeV (Fermi-LAT) e TeV (LHAASO), sugerindo origens físicas distintas que complicam a interpretação de dados de banda larga para buscas de ALPs.

Metodologia
Os autores apresentam um framework abrangente para restringir parâmetros de ALPs ($m_a$ e $g_{a\gamma}$) usando cinco épocas observacionais simultâneas de Mrk 421 de março de 2021 a março de 2024, combinando dados do Fermi-LAT e LHAASO.

1. Formalismo de Propagação: Diferente das aproximações fatorizadas, o estudo emprega um framework de propagação de fótons autoconsistente onde a absorção pela EBL é incorporada diretamente no hamiltoniano de mistura fóton-ALP via um termo absorção ($i/2\lambda_\gamma$). A intensidade da EBL é tratada com um fator de normalização ($K_{EBL}$) como um parâmetro de incerteza (nuisance parameter) com uma priori Gaussiana ($\mu=1.0, \sigma=0.2$) para contabilizar incertezas.
2. Modelagem da Fonte: Dois modelos distintos de emissão de jato são testados para abordar a fraca correlação GeV-TeV em Mrk 421:
   • Um modelo híbrido de duas zonas apresentando uma zona interna compacta (emissão flaring VHE) e uma zona externa estendida (emissão quiescente).
   • Um modelo hadrônico de zona única, fisicamente plausível dada a fraca correlação e amplitudes de variabilidade distintas.
     Ambos os modelos utilizam componentes de campo magnético helicoidais e emaranhados e propagam fótons através da fonte, do meio intergaláctico (IGMF) e do Campo Magnético Galáctico (GMF).
3. Framework Estatístico: Os autores desenvolvem uma análise de verossimilhança de perfil conjunta (joint profile-likelihood). Para abordar a limitação de tratar épocas como independentes, eles introduzem um formalismo de cópula Gaussiana com um coeficiente de correlação conservador ($\rho = 0.03$). Isso modela correlações inter-épocas fracas através de distribuições de excesso $\Delta\chi^2$ transformadas. A análise otimiza parâmetros de incerteza (parâmetros espectrais intrínsecos e $K_{EBL}$) globalmente sob a hipótese de ausência de ALPs antes de realizar o perfil da massa e do acoplamento da ALP.
4. Processamento de Dados: As probabilidades de sobrevivência dos fótons são pré-computadas em uma grade quadridimensional (energia, parâmetros de ALP, normalização da EBL) usando o código gammaALPs e interpoladas durante a avaliação da verossimilhança. Modelos espectrais intrínsecos são selecionados por época usando o Critério de Informação de Akaike corrigido para amostras finitas (AICc).

Contribuições Principais

• Modelo de Propagação Unificado: O artigo implementa um tratamento autoconsistente de mistura fóton-ALP e absorção pela EBL dentro da equação de propagação, evitando os vieses da fatoração multiplicativa.
• Análise de Múltiplas Épocas Correlacionadas: A introdução de uma cópula Gaussiana para lidar com correlações fracas entre períodos observacionais fornece um método estatisticamente consistente para combinar conjuntos de dados de múltiplas épocas, evitando o estreitamento artificial de intervalos de confiança.
• Restrições de Banda Larga: O estudo aproveita a cobertura única de banda larga do Fermi-LAT (GeV) e LHAASO (multi-TeV a PeV) para sondar parâmetros de ALPs em uma ampla faixa de massa (0.1–500 neV).
• Robustez de Modelagem: Ao testar cenários de duas zonas híbridos e de zona única hadrônica, o trabalho contabiliza a física de emissão complexa de Mrk 421, garantindo que as restrições não sejam excessivamente dependentes de uma única suposição de emissão.

Resultados
Usando o framework de verossimilhança de perfil conjunta, os autores derivam limites superiores de nível de confiança (CL) de 95% para o acoplamento fóton-ALP ($g_{a\gamma}$):

• Regime de Baixa Massa ($m_a \lesssim 1$ neV): As restrições são altamente rigorosas, com um limite inferior de $5.0 \times 10^{-14}$ GeV$^{-1}$ e um limite superior de $8.33 \times 10^{-14}$ GeV$^{-1}$. Esta região é sensível às restrições do CAST.
• Massa Intermediária a Alta ($0.1 - 500$ neV):
  • Para o modelo híbrido de duas zonas, o limite superior conservador é $g_{a\gamma} < 6.50 \times 10^{-12}$ GeV$^{-1}$.
  • Para o modelo hadrônico de zona única, o limite superior conservador é $g_{a\gamma} < 7.30 \times 10^{-12}$ GeV$^{-1}$.
• Comportamento: Os contornos de exclusão exibem uma característica de "salto-platô" (jump-plateau). Uma transição abrupta ocorre em torno de $m_a \approx 1$ neV, onde o limite superior aumenta em aproximadamente uma ordem de magnitude, seguido por um relaxamento gradual em massas mais altas. O limite inferior permanece estável em toda a faixa de massa.
• Robustez: Testes com coeficientes de correlação alternativos ($\rho = 0, 0.15, 0.3$) mostram que as variações nas restrições permanecem em menos de uma ordem de magnitude, confirmando a estabilidade dos resultados.

Significância
O artigo afirma que estes resultados beneficiam-se significativamente da cobertura de banda larga de energia muito alta do LHAASO e do monitoramento de longo prazo de Mrk 421. O framework desenvolvido oferece um pipeline robusto e escalável adequado para futuras buscas de ALPs multi-mensageiro com observatórios de raios gama de próxima geração. Ao modelar explicitamente as correlações fracas e tratar a absorção pela EBL de forma autoconsistente, o estudo fornece restrições mais confiáveis e conservadoras sobre o espaço de parâmetros de ALPs, abordando vieses sistemáticos presentes em abordagens fatorizadas anteriores. O trabalho demonstra que a fraca correlação GeV-TeV em Mrk 421 necessita de modelos de emissão complexos (dois modelos de duas zonas ou hadrônicos) para restrições precisas de ALPs, e que a metodologia estatística proposta lida efetivamente com as incertezas associadas.