On gamma rays as predictors of UHECR flux in AGNs

Resumo Técnico: Sobre os raios $\gamma$ como preditores do fluxo de UHECR em AGNs

Definição do Problema
Núcleos Galácticos Ativos (AGNs) são candidatos primários para as fontes de Raios Cósmicos de Ultra-Alta Energia (UHECRs). No entanto, modelos fenomenológicos recentes que utilizam o fluxo observado de raios $\gamma$ de muito alta energia (VHE) como um proxy direto para a luminosidade de UHECR ($L_{CR}$) falham em reproduzir os dados experimentais do Observatório Pierre Auger. Especificamente, esses modelos têm dificuldade em ajustar simultaneamente o espectro de energia, a composição de massa e o mapa de direção de chegada. Um ponto crítico de falha é a superestimação da contribuição do blazar Markarian 421 (Mkn 421), o que cria um hotspot previsto e uma direção de dipolo inconsistente com as observações. Além disso, esses modelos frequentemente produzem um dipolo de direção "espúrio" e um ajuste pobre à forma espectral, particularmente quando misturam diferentes classes de AGN (por exemplo, galáxias de rádio e blazares) na mesma análise.

Metodologia
Os autores propõem que a discrepância decorre de uma suposição implícita em trabalhos anteriores: a de que a emissão de UHECR é direcionada (beamed) da mesma maneira que a emissão de raios $\gamma$ de jatos relativísticos. O artigo argumenta que, enquanto os raios $\gamma$ estão sujeitos a um forte aumento Doppler ($D^4$) devido ao movimento relativístico de "blobs" dentro do jato, os UHECRs provavelmente sofrem isotropização (decolimação) antes de escaparem do ambiente da fonte.

O estudo procede através das seguintes etapas:

1. Revisão Teórica: Os autores revisam os mecanismos de aceleração em jatos de AGN (choques, reconexão, aceleração de cisalhamento) e a propagação de partículas através dos lobos do jato. Eles argumentam que os campos magnéticos nos lobos e no coco da fonte espalham os UíHECRs, isotropizando efetivamente sua emissão mesmo que se originem em fluxos relativísticos.
2. Correção de Proxy: Em vez de usar a luminosidade de raios $\gamma$ observada ($L^{obs}_{\gamma}$), os autores propõem o uso da luminosidade intrínseca de raios $\gamma$ ($L^{int}_{\gamma}$) ou da potência radiativa do jato como o proxy para $L_{CR}$. Isso envolve corrigir o fluxo observado pelo fator Doppler ($D$) elevado à potência $q$, onde $L^{int}_{\gamma} = D^{-q}L^{obs}_{\gamma}$. Dois cenários são testados: $q=4$ (assumindo que $L_{CR} \propto L^{int}_{\gamma}$) e $q=2$ (assumindo que $L_{CR} \propto$ potência do jato radiativo).
3. Simulação e Ajuste: Utilizando o framework CRPropa3, os autores realizam simulações de propagação 1-D para cinco núcleos primários (H, He, N, Si, Fe) através de uma grade de energias e distâncias. Eles comparam o fluxo simulado contra os dados do Pierre Auger (espectro e momentos de $X_{max}$) e mapas de direção de chegada.
4. Análise Comparativa: Três casos são comparados:
   • Caso (a): Usando $L^{obs}_{\gamma}$ como o proxy (abordagem padrão).
   • Caso (b): Usando $D^{-2}L^{obs}_{\gamma}$ como o proxy.
   • Caso (c): Usando $D^{-4}L^{obs}_{\gamma}$ como o proxy.

Principais Contribuições e Resultados

• Supressão da Dominância de Blazares: A correção reduz significativamente o peso de blazars altamente direcionados (como Mkn 421) em relação às galáxias de rádio (como Cen A, M87 e Fornax A). No modelo padrão, Mkn 421 domina o fluxo; nos modelos corrigidos, sua contribuição é suprimida de 2 a 5 ordens de magnitude, tornando-o negligenciável.
• Melhoria no Ajuste de Espectro e Composição: O ajuste combinado dos dados de espectro de energia e composição melhora significativamente quando se utiliza os proxies de luminosidade intrínseca. O qui-quadrado reduzido ($\chi^2/NDF$) diminui de 4,6 (usando $L^{obs}_{\gamma}$) para 3,1 (usando $D^{-2}L^{obs}_{\gamma}$). Esta melhoria é atribuída a uma contribuição relativa mais realista de fontes próximas versus fontes distantes, evitando as distorções espectrais causadas pela sobrecontribuição de fontes distantes e altamente direcionadas.
• Resolução da Tensão do Dipolo: A tensão entre as direções de dipolo modeladas e observadas é reduzida. Para o campo de visão de Auger em energias $E > 8$ EeV, a discrepância cai de 5,9$\sigma$ para 3,5$\sigma$. Para $E > 32$ EeV, a tensão reduz de 2,1$\sigma$ para 1,1$\sigma$.
• Realinhamento de Hotspots: Os mapas de direção de chegada previstos mudam de serem dominados por Mkn 421 para mostrar hotspots ao redor de Cen A, Fornax A, M87 e blazars distantes específicos (Mkn 180, 1ES2344+514). Esses novos hotspots alinham-se mais proximamente com os dois hotspots relatados pelo Observatório Pierre Auger e o um relatado pelo Telescope Array.

Significância e Alegações
O artigo afirma oferecer uma solução para as contradições encontradas em modelos anteriores de fontes de UHECR baseados em AGN ao desacoplar as suposições de direcionamento (beaming) de raios $\gamma$ e UHECRs. Os autores afirmam que:

1. A associação entre raios $\gamma$ e UHECRs é fraca, mas não nula; embora os raios $\gamma$ rastreiem regiões de aceleração de partículas, os UHECRs são provavelmente isotropizados por campos magnéticos nos lobos da fonte antes de escaparem.
2. Usar o fluxo de raios $\gamma$ observado como um proxy assume implicitamente que os UHECRs são direcionados, uma suposição que os autores argumentam ser fisicamente improvável para UHECRs escapando de ambientes magnéticos complexos.
3. Corrigir para a luminosidade intrínseca ou potência do jato radiativo ($D^{-2}L^{obs}_{\gamma}$ ou $D^{-4}L^{obs}_{\gamma}$) fornece um proxy mais robusto que reconcilia os dados de espectro, composição e direção de chegada do Observatório Pierre Auger.

Os autores concluem que, embora seu modelo seja um ajuste combinado simplificado (não contabilizando distribuições completas de $X_{max}$ ou deflexões detalhadas de campos magnéticos), os resultados fortalecem a hipótese de AGNs como fontes de UHECR e fornecem um framework mais consistente para interpretar dados de múltiplos mensageiros de diferentes classes de AGN. Eles observam que trabalhos futuros devem abordar a decolimação dependente da energia e os efeitos específicos de campos magnéticos galácticos e extragalácticos regulares.