Nuclei in high-energy neutrino sources: A multimessenger study of in-source propagation

Resumo Técnico: Núcleos em Fontes de Neutrinos de Alta Energia: Um Estudo Multimensageiro de Propagação Intra-fonte

Enunciado do Problema
A descoberta de neutrinos astrofísicos de alta energia pelo IceCube levantou questões críticas sobre a natureza de suas fontes. Enquanto os modelos padrão assumem que esses neutrinos se originam de interações próton-próton ($pp$) ou próton-fóton ($p\gamma$), a falta de contrapartidas de raios gama TeV–GeV sugere que essas fontes são "escondidas" ou opacas a fótons de alta energia. Estudos anteriores assumiram amplamente uma composição puramente de prótons para as partículas aceleradas. No entanto, a composição química das partículas injetadas (por exemplo, núcleos leves vs. pesados) permanece uma questão em aberto. Se a fonte contiver uma fração significativa de núcleos, a propagação intra-fonte envolve cascatas nucleares complexas (fotodesintegração) e processos eletromagnéticos que alteram os rendimentos resultantes de neutrinos e raios gama. O problema central abordado é como o relaxamento da hipótese de prótons puros e a consideração da composição nuclear afetam a fenomenologia multimensageira, especificamente a razão entre os fluxos de neutrinos e raios gama e a luminosidade necessária da fonte.

Metodologia
Os autores realizam simulações detalhadas de Monte Carlo de cascatas nucleares e eletromagnéticas dentro de uma fonte astrofísica de alta energia, utilizando a galáxia Seyfert NGC 1068 como um estudo de caso representativo. A simulação é dividida em duas etapas:

1. Interações Fotonucleares e Cascata Nuclear: O código rastreia a evolução de núcleos primários (de $Z=1$ a $Z=90$) através de um campo de fótons de fundo. Inclui:
   • Fotodesintegração (PD): Modelada usando seções de choque TALYS/TENDL-2023 para núcleos pesados e parametrizações dedicadas para núcleos leves ($A \le 9$). Isso inclui canais de ruptura exclusiva (ex: $(\gamma, n)$, $(\gamma, p)$, $(\gamma, \alpha)$) e fluorescência de ressonância nuclear.
   • Produção de Fotomésons (PM): Modelada usando o gerador de eventos SOPHIA com uma aproximação de superposição para núcleos.
   • Produção de Pares de Bethe-Heitler (BH): Tratada como um processo de perda de energia contínua para núcleos carregados.
   • Produção de Pares com Captura (PPC): Uma inclusão inovadora de produção de pares onde o elétron é capturado na camada L (PPCL), seguida de desexcitação atômica imediata emitindo fótons de linha.
   • Decaimento Nuclear: Incluído para completude, mas considerado negligenciável para os intervalos de tempo considerados.
2. Cascata Eletromagnética: Fótons secundários, elétrons e pósitrons da etapa (i) são propagados através do campo de radiação. Isso inclui o Espalhamento Compton Inverso (ICS) padrão e a Produção de Pares por Elétrons (EPP), bem como processos raros como produção de pares dupla (DPP) e produção de pares de múons (MPP), embora este último seja observado como negligenciável nos intervalos de energia estudados.

A fonte é modelada como uma esfera homogênea de raio $R$ imersa em um campo de fótons de disco-coroa (disco UV + coroa de raios X). O estudo varia o tamanho da fonte ($R$) e a composição nuclear injetada (prótons, nêutrons e vários núcleos como Carbono, Ferro e Zircônio) para mapear as profundidades ópticas e a partição de energia.

Principais Contribuições

• Estrutura Numérica: O artigo apresenta um novo tratamento numérico de cascatas induzidas por núcleos que estende os processos fotonucleares até isótopos de $Z=90$ e contabiliza explicitamente a produção de pares com captura na camada L, um processo anteriormente negligenciado em tais contextos.
• Análise Sistemática de Composição: Fornece uma análise abrangente de como a composição nuclear injetada altera a razão de rendimento neutrino-raio gama, demonstrando que essa relação é não-monotônica e altamente dependente da profundidade óptica da fonte (densidade de fótons e tamanho).
• Reanálise do COMPTEL: Os autores realizam uma reanálise dedicada dos dados arquivais do COMPTEL para NGC 1068 usando frameworks modernos (GammaLib/ctools), derivando limites atualizados de nível de confiança de 95% na banda de MeV que são aproximadamente um fator de 2 mais rigorosos do que os limites anteriores nos bins de menor energia.

Resultos
O estudo identifica dois regimes distintos baseados na compacidade da fonte (densidade de fótons):

1. Regime Compacto/Denso ($R \lesssim 2 R_S$):

   • A fonte é opticamente espessa para fotodesintegração. Núcleos pesados são rapidamente fragmentados em nucleons livres (prótons e nêutrons) antes de escaparem.
   • Uma fração significativa da energia é transferida para nêutrons livres. Como os nêutrons não sofrem perdas de energia de Bethe-Heitler (BH), a energia disponível para a produção de neutrinos é preservada.
   • Resultado: A injeção de núcleos pesados é mais eficiente na produção de neutrinos por unidade de luminosidade de raios cósmicos do que os prótons. Consequentemente, para corresponder ao fluxo de neutrinos observado pelo IceCube, é necessária uma luminosidade de raios cósmicos total menor para núcleos pesados. A cascata eletromagnética resultante (raios gama de MeV–GeV) é mais fraca (em ~28–35%) comparada ao caso do próton, relaxando assim as restrições das observações de raios gama.

2. Regime Estendido/Rarefeito ($R \gtrsim 5 R_S$):

   • A fotodesintegração é incompleta; uma fração substancial da energia injetada permanece bloqueada em fragmentos nucleares carregados durante a travessia da fonte.
   • Esses fragmentos pesados experimentam perdas de energia BH aumentadas (escalando como $Z^2/A$) em comparação com prótons na mesma energia por núcleon.
   • Resultado: A injeção de núcleos pesados é menos eficiente na produção de neutrinos. Para corresponder ao fluxo do IceCube, é necessária uma luminosidade de raios cósmicos maior (aumento de 20–30%). Além disso, as perdas BH aumentadas e a maior luminosidade levam a uma cascata eletromagnética significativamente mais forte (fator de ~2.0–2.5 maior no fluxo de raios gama), fazendo com que os cenários de núcleos pesados ultrapassem os limites de raios gama existentes de forma mais severa do que os cenários de prótons.

Significância e Alegações
O artigo afirma que a composição nuclear das partículas aceleradas é um parâmetro crítico, embora frequentemente negligenciado, na modelagem de fontes de neutrinos de alta energia. A significância do trabalho reside em demonstrar que a interpretação dos dados multimensageiros (neutrinos vs. raios gama) está intrinsecamente ligada à compacidade e à composição injetada da fonte.

Especificamente, os autores argumentam que:

• A banda de raios gama de MeV é uma ferramenta diagnóstica poderosa para fontes de neutrinos "escondidas", pois a cascata eletromagnética reprocessa fótons de alta energia para esta banda.
• Os dados atuais de raios gama (especificamente os limites do COMPTEL reanalisados) já fornecem restrições leves, mas significativas, sobre a composição e o tamanho da fonte.
• Em fontes compactas, a injeção de núcleos pesados pode explicar o fluxo de neutrinos com menos potência total e menos emissão de raios gama, potencialmente resolvendo a tensão com as não-detecções de raios gama. Por outro lado, em fontes menos compactas, os núcleos pesados são desfavorecidos pelos limites de raios gama.
• Futuras instalações de raios gama de MeV (ex: COSI, AMEGO-X, e-ASTROGAM) serão essenciais para refinar essas sondagens e distinguir entre cenários de injeção de prótons e nucleares na emergente era da astrofísica multimensageira.

Os autores permanecem modestos, observando que seu modelo é um "modelo-brinquedo" simplificado do sistema disco-coroa e que trabalhos futuros precisarão incorporar perdas de sincrotron e geometrias mais complexas (2D/3D) para capturar totalmente a física da fonte.