Modelling Galactic neutrino emission: contributions from massive star clusters and interstellar cosmic rays

Resumo Técnico: Modelagem da Emissão de Neutrinos Galácticos de Aglomerados Estelares Massivos e Raios Cósmicos Interestelares

Definição do Problema
A origem dos raios cósmicos (RC) mais energéticos do Universo, particularmente aqueles que atingem energias de PeV (PeVatrons), permanece incerta. Embora os remanescentes de supernovas (SNRs) sejam tradicionalmente considerados os principais aceleradores, as observações de raios gama sugerem que eles têm dificuldade em atingir energias de PeV, necessitando de classes de fontes alternativas. Aglomerados Estelares Massivos Jovens (YMSCs) foram propostos como candidatos viáveis a PeVatrons devido ao aumento da turbulência magnética gerada pelo choque de ventos estelares colidentes e múltiplos eventos de SNR. Concomitantemente, o Observatório IceCube detectou um fluxo de neutrinos de alta energia do Plano Galáctico. No entanto, as análises atuais dependem de modelos que consideram apenas a emissão difusa pura de neutrinos gerada pela interação do mar de RC com o gás interestelar. Esses modelos não resolvem explicitamente a contribuição de aceleradores hadrônicos não resolvidos, como a população de YMSCs, o que pode mascarar as assinaturas específicas dessas fontes.

Metodologia
Os autores apresentam um novo arcabouço teórico para computar a emissão de neutrinos esperada da Via Láctea, integrando dois componentes distintos:

1. População Sintética de YMSCs: Para superar a limitação de observar YMSCs apenas em nossa vizinhança solar local ($\lesssim 2$ kpc), os autores geraram 100 realizações independentes da população galáctica de YMSCs. Esta população sintética amostra massas de aglomerados ($10^3 - 6,3 \times 10^4 M_\odot$), idades ($<30$ Myr) e distribuições espaciais (seguindo braços espirais, a Barra Galáctica e o Espúrio Local) baseadas em restrições observacionais.
2. Aceleração e Transporte de Partículas: Para cada aglomerado sintético, o modelo simula a população estelar (Sequência Principal, Wolf-Rayet e progenitores de Supernova) para calcular a luminosidade mecânica coletiva de ventos e supernovas. Os autores modelam as bolhas de vento sopradas resultantes, calculando os perfis de densidade de gás e de campo magnético.
   • Mecanismos de Aceleração: A aceleração de partículas é modelada em dois locais: o choque de terminação de vento coletivo (WTS) e os choques de remanescentes de supernova (SNR).
   • Cenários de Difusão: A propagação de RCs dentro das bolhas de aglomerados é computada sob três suposições de turbulência distintas: difusão de Kolmogorov, Kraichnan e Bohm.
   • Produção de Neutrinos: Os fluxos de neutrinos são calculados a partir das colisões hadrônicas de prótons acelerados com o gás alvo (principalmente na casca varrida densa das bolhas).
3. Componente Galáctico Difuso: O modelo também computa a emissão difusa de neutrinos decorrente da interação do mar de RC galáctico com o gás interestelar. Isso utiliza modelos de propagação de RC (incorporando o endurecimento espectral no interior da Galáxia) e mapas de gás (GALPROP e um modelo 3D alternativo).
4. Geração de Templates: O fluxo total previsto é a soma da contribuição de YMSC e do componente difuso induzido por RC. Eles são sintetizados em templates espaciais e espectrais para comparação com os dados do IceCube.

Principais Contribuições

• Novos Templates de Neutrinos: O artigo fornece os primeiros templates de neutrinos abrangentes que incorporam explicitamente a emissividade da população de YMSC junto ao componente difuso padrão induzido por RC.
• Modelagem de População Sintética: Oferece uma estimativa estatisticamente robusta da produção cumulativa de neutrinos de toda a população de YMSC da Galáxia, em vez de depender apenas dos poucos aglomerados próximos e resolvidos.
• Análise de Múltiplos Componentes: O arcabouço distingue entre a produção "in-situ" de neutrinos dentro dos aglomerados e a emissão difusa de fundo, permitindo a separação das contribuições das fontes.
• Sensibilidade de Parâmetros: O estudo explora a dependência do fluxo previsto nos regimes de difusão de partículas (Kolmogorov, Kraichnan, Bohm) e no endurecimento espectral de RC no interior da Galáxia.

Resultados

• Detectabilidade de Fontes Individuais: O fluxo previsto de YMSCs individuais é considerado marginalmente detectável mesmo por detectores de escala quilométrica (ex: IceCube, KM3NeT). Consequentemente, a emissão cumulativa da população de YMSCs deve aparecer como um componente difuso não resolvido.
• Magnitude do Fluxo: A normalização dos modelos é consistente com o melhor ajuste (best-fit) do IceCube para os templates galácticos existentes. Os autores descobrem que a contribuição não resolvida da emissão de aglomerados pode não ser desprezível.
• Dependência da Difusão:
  • No cenário de difusão de Bohm (confinamento mais eficiente), os YMSCs dominam o fluxo de neutrinos abaixo de $\sim 50$ TeV, potencialmente superestimando o fluxo em comparação com os dados do IceCube nos cenários de Kraichnan e Kolmogorov, caso nenhum endurecimento espectral seja assumido.
  • Nos cenários de Kolmogorov e Kraichnan, o componente difuso induzido por RC geralmente domina ou fornece contribuições comparáveis à componente de YMSC, particularmente em energias mais altas.
• Distribuição Espacial: A contribuição de YMSC está concentrada na faixa de longitude galáctica interna ($|l| \le 30^\circ$, a Crista Galáctica), enquanto a emissão induzida por RC domina em longitudes maiores.
• Endurecimento Espectral: Incluir o endurecimento espectral do mar de RC no interior da Galáxia reduz a contribuição relativa de YMSCs nos cenários de Kolmogorov e Kraichnan, trazendo o espectro total previsto para uma melhor concordância com as observações do IceCube a $\sim 100$ TeV.
• Consistência: As previsões do modelo para a contraparte de raios gama das colisões hadrônicas permanecem consistentes com as observações atuais de raios gama difusos na faixa de TeV-PeV.

Significância e Alegações
Os autores alegam que seu trabalho fornece uma evolução necessária na modelagem teórica da emissão de neutrinos galácticos. Ao ir além dos templates puramente difusos para incluir a emissividade específica da população de YMSC, eles oferecem uma ferramenta para "desmembrar" os dois componentes da emissão de neutrinos galácticos. O artigo postula que a aplicação desses novos templates aos futuros conjuntos de dados de neutrinos permitirá à comunidade restringir o papel dos aglomerados estelares na aceleração extrema de RC e na produção de neutrinos. Os autores enfatizam que seus resultados sugerem que a contribuição não resolvida dos aglomerados é provavelmente não desprezível e que o ajuste conjunto das contribuições das fontes e do componente difuso é crucial para identificar aceleradores galácticos e resolver a origem dos raios cósmicos. O trabalho não pretende ter identificado definitivamente os YMSCs como a fonte do sinal do IceCube, mas sim fornece o arcabouço teórico necessário para testar essa hipótese contra dados observacionais.