Chasing the neutrino blazar candidates II: SED modeling with hadronic model

Resumo Técnico: Caçando candidatos a blazar de neutrinos II: Modelagem SED com modelo hadrônico

Declaração do Problema
Embora os blazares sejam considerados candidatos promissores para fontes de neutrinos de alta energia, a origem física de sua emissão de neutrinos permanece incerta. Embora o cenário leptônico (espalhamento Compton inverso) seja o quadro padrão para modelar as distribuições de energia espectral (SEDs) de blazares, a detecção de neutrinos de alta energia sugere um papel significativo para hádrons relativísticos. Trabalhos anteriores dos autores (Artigo I) identificaram 199 candidatos a blazares de neutrinos (NBCs) usando uma abordagem de aprendizado com poucos exemplos. Este estudo visa expandir esse trabalho investigando as propriedades físicas desses candidatos por meio de modelagem SED de banda larga dentro de um quadro hadrônico, especificamente para estimar a produção máxima possível de neutrinos e restringir a composição de partículas dos jatos de blazares.

Metodologia
Os autores construíram SEDs de banda larga para 103 NBCs (selecionados dos 199 originais com base na cobertura multiespectral suficiente) e para a fonte confirmada TXS 0506+056. A modelagem empregou um código hadrônico dependente do tempo de uma única zona (AM3) com as seguintes suposições e procedimentos principais:

• Domínio Hadrônico: Para estimar o limite superior da emissão de neutrinos, o modelo assume que fótons de alta energia (HE) e de muito alta energia (VHE) são produzidos principalmente via interações $p\gamma$ e cascatas eletromagnéticas associadas. A contribuição do espalhamento Compton inverso (IC) leptônico é fortemente suprimida.
• Configuração Física: A região de emissão é modelada como um "blob" esférico movendo-se com um fator de Lorentz de bulk $\Gamma$ em um campo magnético $B$. Elétrons e prótons não térmicos são injetados com espectros de lei de potência.
• Processos Modelados: O modelo considera radiação síncrotron, produção de píons foto-pion ($p\gamma \to n\pi^+$ ou $p\pi^0$), produção de pares Bethe-Heitler ($p\gamma \to pe^+e^-$), produção de pares fóton-fóton e o subsequente decaimento de píons e múons para produzir neutrinos. A atenuação pela luz de fundo extragaláctica (EBL) é aplicada a raios gama de alta energia.
• Estratégia de Ajuste: Os autores utilizaram o algoritmo de Otimização por Enxame de Partículas (PSO) via o pacote PySwarms para ajustar nove parâmetros livres (por exemplo, intensidade do campo magnético, tamanho da região, luminosidades de injeção de partículas e índices espectrais) contra dados multiespectrais que abrangem as faixas de rádio a $\gamma$-ray. Dados de rádio abaixo de $10^{11}$ Hz foram excluídos devido às limitações de auto-absorção síncrotron no modelo de uma única zona. Um refinamento de "ajuste visual" foi aplicado para garantir representações SED aceitáveis.
• Análise de Correlação: Uma análise de correlação parcial foi realizada para investigar as relações entre a luminosidade de neutrinos ($L_\nu$) e as luminosidades eletromagnéticas (rádio, óptico, raios X, $\gamma$-ray), bem como as propriedades do motor central (massa do buraco negro $M_{BH}$ e luminosidade do disco de acreção $L_{disk}$), controlando a dependência do desvio para o vermelho.

Principais Contribuições e Resultados

1. Modelagem SED e Restrições de Parâmetros: O estudo restringiu com sucesso nove parâmetros-chave para 103 NBCs. A modelagem revela que a intensidade do campo magnético tende a ser alta (66% das fontes têm $\log B \ge 1.0$) para dissipar eficientemente a energia dos elétrons via radiação síncrotron, suprimindo assim a emissão IC e permitindo que os prótons dominem o pico de alta energia. A potência carregada por prótons relativísticos ($P_p$) domina sobre a potência dos elétrons ($P_e$) e a potência do campo magnético ($P_B$), com uma razão média $\log(P_p/P_B) \simeq 1.67$.
2. Previsão da Faixa de MeV: Uma característica distintiva do modelo é a previsão de emissão síncrotron de prótons proeminente, com pico na faixa de MeV. Para 99 dos 103 NBCs, o pico de síncrotron de prótons cai dentro de 0,1–100 MeV. Em 90 dessas fontes, o fluxo de síncrotron de prótons na faixa de MeV excede o fluxo de produção de pares Bethe-Heitler na faixa de GeV.
3. Correlações com Emissão Eletromagnética:
   • Raios Gama e Raios X: Correlações positivas fortes foram encontradas entre a luminosidade de neutrinos e as luminosidades de $\gamma$-ray ($r=0,78$) e raios X ($r=0,52$). Isso apoia uma origem hadrônica onde prótons de alta energia produzem tanto neutrinos quanto fótons de alta energia.
   • Óptico (Faixa R): Existe uma correlação fraca/moderada ($r=0,30$) entre a luminosidade de neutrinos e a luminosidade óptica na faixa R. Os autores sugerem que isso pode surgir porque fótons síncrotron ópticos do jato servem como fótons semente para a produção de píons foto-pion, ou potencialmente devido a uma conexão com o disco de acreção.
   • Rádio: Nenhuma correlação significativa foi encontrada entre a luminosidade de neutrinos e a luminosidade de rádio em 1,4 GHz.
   • Motor Central: Nenhuma correlação estatisticamente significativa foi encontrada entre a luminosidade de neutrinos e a massa do buraco negro ($M_{BH}$) ou a luminosidade do disco de acreção ($L_{disk}$), embora uma tendência positiva fraca com $L_{disk}$ fosse visualmente aparente.
4. Perspectivas de Detecção: Com base nos fluxos máximos de neutrinos previstos pelo modelo:
   • IceCube: 3 NBCs são potencialmente detectáveis.
   • Instalações Futuras: Até 22 fontes podem ser detectáveis pelo KM3NeT, 45 pelo NEON e 62 pelo TRIDENT.
   • O fluxo de neutrinos somado da amostra não excede os limites atuais de empilhamento do IceCube, sugerindo que o cenário de "máximo neutrino" é fisicamente viável.

Significado e Alegações
O artigo afirma que seus resultados fornecem previsões testáveis para futuras observações multimensageiras e oferecem novos insights sobre a composição dos jatos de blazares. Especificamente:

• A Janela de MeV: Os autores destacam a faixa de MeV como uma janela observacional crítica para distinguir entre cenários leptônicos e hadrônicos, pois o modelo hadrônico prevê um pico de síncrotron de prótons distinto nesta faixa, ausente nos modelos leptônicos padrão.
• Composição do Jato: A modelagem apoia um cenário onde prótons relativísticos dominam o orçamento de energia dos jatos nestes candidatos, consistente com a evidência multimensageira de TXS 0506+056.
• Limites Superiores: Os resultados representam um limite superior na emissão de neutrinos; os fluxos reais podem ser menores, o que explica por que muitos candidatos ainda não foram confirmados como emissores de neutrinos, apesar de estarem acima das curvas de sensibilidade atuais no modelo.
• Observações Futuras: O estudo enfatiza que missões de raios gama de MeV de próxima geração (por exemplo, e-ASTROGAM, AMEGO) e detectores de neutrinos de próxima geração (KM3NeT, NEON, TRIDENT) são essenciais para verificar esses cenários hadrônicos e restringir a composição de partículas dos jatos de blazares.

Os autores permanecem modestos quanto à origem da correlação entre neutrinos e óptico, notando que ela pode decorrer de efeitos de seleção de amostra ou múltiplos mecanismos (jato vs. disco), e reconhecem que a exigência de potências de jato super-Eddington em algumas fontes pode exigir modelos de múltiplas zonas para uma reconciliação completa.