High-energy Multi-messenger Emission from Galaxy Clusters in the Local Universe

Resumo Técnico: Emissão Multi-mensageira de Alta Energia de Aglomerados de Galáxias no Universo Local

Enunciado do Problema
A origem dos fundos difusos de neutrinos e raios gama de alta energia permanece um grande problema não resolvido na astrofísica. Embora núcleos galácticos ativos (AGN) e galáxias com formação estelar explosiva (starburst) sejam contribuintes conhecidos, o papel específico dos aglomerados de galáxias como fontes desses mensageiros cósmicos não está totalmente restringido. Aglomerados de galáxias, caracterizados por gás quente do meio intra-aglomerado (ICM) ($\sim 10^8$ K) e campos magnéticos fortes ($\sim 1$ $\mu$G), são teoricamente capazes de acelerar raios cósmicos (CR) via choques de fusão, choques de acreção e turbulência. Esses prótons acelerados podem interagir com o ICM para produzir raios gama e neutrinos secundários via processos hadrônicos. No entanto, estudos anteriores frequentemente recorreram a modelos simplificados, esfericamente simétricos, ou assumiram injeção de CR exclusivamente a partir dos núcleos centrais dos aglomerados, potencialmente distorcendo a distribuição espacial e o fluxo total da emissão difusa. Além disso, distinguir a emissão difusa do aglomerado da emissão de fontes pontuais (por exemplo, de AGN centrais como NGC 1275 ou M87) permanece desafiador devido à resolução angular limitada nas observações atuais.

Metodologia
Este trabalho emprega um quadro numérico rigoroso combinando simulações magnetohidrodinâmicas (MHD) de alta resolução com métodos de transporte Monte Carlo para modelar a emissão multi-mensageira de aglomerados de galáxias no Universo local ($z \lesssim 0.1$).

1. Simulações MHD (SLOW): O estudo utiliza a simulação "Simulating the LOcal Web" (SLOW), que incorpora condições iniciais constrangidas por observações baseadas no catálogo CosmicFLow. Diferentemente de modelos anteriores, o SLOW fornece uma descrição 3D autoconsistente e de alta resolução das distribuições de densidade de gás, temperatura e campo magnético dentro dos aglomerados (especificamente Virgo, Perseu e Coma), capturando morfologias não esféricas e subestruturas decorrentes de fusões.
2. Injeção de Raios Cósmicos: Adota-se um modelo de injeção fisicamente motivado onde a injeção de CR segue a distribuição espacial da densidade de gás do aglomerado. Esta abordagem contabiliza a aceleração de CR em regiões de alta densidade associadas a turbulência e choques em todo o volume do aglomerado, em vez de restringir a injeção ao núcleo do aglomerado. O orçamento total de energia de CR é parametrizado pela razão $X_{CR} = E_{CR}/E_{cluster} \simeq 0,01$, consistente com os limites superiores observacionais.
3. Transporte de CR e Produção Secundária: A propagação de CRs (assumidos como prótons) dentro do ICM e do meio intergaláctico (IGM) é simulada usando o código Monte Carlo CRPropa. O código resolve a equação de transporte de Parker, contabilizando advecção, difusão espacial, difusão de momento e perdas de energia. Canais de interação chave incluem colisões próton-próton ($pp$) (dominantes abaixo de $10^{17}$ eV), produção de fotopíons e cascateamento eletromagnético (espalhamento Compton inverso sobre a CMB e EBL, perdas por síncrotron). Raios gama e neutrinos secundários são gerados de forma autoconsistente a partir dessas interações.

Contribuições Principais

• Modelagem Ambiental Realista: A integração de simulações MHD do SLOW permite a modelagem da propagação de CR em ambientes de aglomerados realistas e não esféricos, superando as suposições simplificadas de simetria esférica de estudos analíticos ou semi-analíticos anteriores.
• Esquema de Injeção Estendido: A adoção de um modelo de injeção de CR espacialmente estendido, traçando o perfil de densidade de gás, representa uma melhoria metodológica significativa sobre trabalhos anteriores que assumiam injeção apenas central. Isso reduz o tempo de confinamento efetivo de CR e as taxas de interação no núcleo, levando a um perfil de emissão mais difuso.
• Estimativa Abrangente de Fluxo: O estudo calcula tanto os fluxos individuais de aglomerados (Virgo, Perseu, Coma) quanto o fluxo difuso cumulativo de toda a população de aglomerados massivos ($M \gtrsim 5 \times 10^{13} M_\odot$) dentro de uma distância comóvel de 500 Mpc.

Resultados

• Aglomerados Individuais: Os fluxos de raios gama difusos previstos para os aglomerados Virgo, Perseu e Coma estão bem abaixo dos limites superiores atuais relatados pela colaboração LHAASO. Os fluxos previstos também são inferiores à emissão de raios gama observada pelo MAGIC a partir da fonte central NGC 1275 em Perseu, consistente com a expectativa de que o MAGIC detecta uma fonte pontual enquanto o modelo prevê emissão difusa.
• Fluxo Cumulativo: O fluxo integrado de raios gama de todos os aglomerados dentro de 500 Mpc é significativamente menor que o fundo difuso de raios gama medido pelo Fermi-LAT. O estudo observa que, em energias abaixo de $\sim 100$ GeV, outras fontes (blazares, AGN desalinhados, galáxias starburst) dominam o fundo.
• Emissão de Neutrinos: O fluxo de neutrinos de todos os sabores previsto da população local de aglomerados está bem abaixo do fluxo difuso observado pelo IceCube, mas é comparável aos limites superiores existentes do IceCube especificamente para aglomerados de galáxias ($M > 10^{14} M_\odot$). O fluxo de neutrinos previsto do aglomerado Virgo é encontrado como comparável à sensibilidade projetada do detector KM3NeT/ARCA.
• Características Espectrais: Os índices espectrais dos fluxos previstos ($\alpha \approx 2,0 - 2,5$) são consistentes com as expectativas para aceleração por choques de fusão e acreção. Os resultados indicam que aglomerados de menor massa contribuem eficientemente para a emissão de raios gama de $\sim 10$ GeV, enquanto sistemas massivos dominam o fluxo de neutrinos de alta energia ($\gtrsim 10$ TeV) devido às suas capacidades superiores de confinamento de CR.

Significado e Afirmações
O artigo afirma que seus resultados são consistentes com as restrições observacionais atuais, particularmente a não detecção de emissão difusa significativa pelo LHAASO e os limites superiores estabelecidos pelo IceCube. Os autores postulam que a tensão entre previsões de modelos anteriores (que frequentemente assumiam fluxos mais altos) e os limites superiores do IceCube é substancialmente reduzida pela adoção do modelo de injeção espacialmente estendido e pelo uso de simulações SLOW de alta resolução.

O estudo conclui que, embora os aglomerados de galáxias sejam improváveis de dominar o fundo total difuso de neutrinos de alta energia, eles podem representar uma fração não desprezível dele. O trabalho fornece um quadro teórico que motiva buscas dedicadas por neutrinos difusos do aglomerado Virgo usando KM3NeT e IceCube-Gen2, e por raios gama difusos usando CTA e LHAASO. Os autores enfatizam que suas descobertas se alinham com as perspectivas de futuros observatórios multi-mensageiros, oferecendo um teste oportuno dos processos de aceleração hadrônica de CR em estruturas de grande escala.