A MINOT-based Study of Gamma-ray emission from SPT-CL J2012-5649/Abell 3667

Resumo Técnico: Um Estudo Baseado em MINOT da Emissão de Raios Gama de SPT-CL J2012-5649/Abell 3667

Declaração do Problema
Os aglomerados de galáxias são as maiores estruturas virializadas do universo, servindo como reservatórios de raios cósmicos (RC) de alta energia acelerados por choques impulsionados por fusões. Embora as buscas iniciais do Fermi-LAT por emissão de raios gama de aglomerados tenham produzido resultados nulos, estudos recentes relataram detecções tentativas, incluindo um sinal de $\sim6.1\sigma$ do aglomerado em fusão SPT-CL J2012-5649 (Abell 3667, daqui em diante A3667). No entanto, a origem dessa emissão permanece incerta devido à presença de galáxias de rádio próximas dentro da função de espalhamento pontual do Fermi-LAT e a uma tensão significativa entre o índice espectral observado ($\Gamma \approx -3.6$) e as previsões hadrônicas teóricas ($\Gamma \approx -2.4$). Além disso, as contribuições relativas dos processos hadrônicos (interações próton-próton) versus processos leptônicos (espalhamento Compton inverso de elétrons de raios cósmicos) exigem quantificação rigorosa para interpretar a física não térmica do meio intra-aglomerado (ICM).

Metodologia
Os autores empregam a estrutura MINOT (Modelagem e Interpretação de Emissão Não Térmica de Aglomerados de Galáxias do Tipo Via Láctea) para construir um modelo físico autoconsistente de A3667. O estudo utiliza as seguintes entradas e procedimentos:

• Dados Termodinâmicos: Perfis de densidade e temperatura de elétrons são derivados de observações de raios X do Chandra de arquivo (catálogo ACCEPT), cobrindo o núcleo do aglomerado até $0.22 R_{500}$.
• Extrapolação de Perfis: Os perfis termodinâmicos observados são extrapolados para um raio de truncamento ($R_{trunc} = 5000$ kpc $\approx 3.7 R_{500}$) usando um modelo $\beta$ para densidade de elétrons e um modelo de lei de potência para temperatura.
• Configuração de Física Não Térmica:
  • Campo Magnético: Modelado como escalando com a densidade de elétrons ($B \propto n_e^{\eta_B}$) com $\eta_B = 0.5$ e um campo de referência de $5 \mu$G em 100 kpc.
  • Raios Cósmicos: As distribuições espaciais de prótons ($CR_p$) e elétrons ($CRe$) são escaladas em relação à densidade térmica ($n_{CR} \propto n_e^{\eta_{CR}}$). A linha de base assume escalamento isodens ($\eta_{CR}=1$).
  • Espectros: $CR_p$ seguem uma lei de potência ($dN/dp \propto p^{-2.4}$). $CRe$ primários seguem um espectro de injeção contínua com uma quebra de resfriamento em 5 GeV. $CRe$ secundários são computados de forma autoconsistente via interações $pp$ usando o modelo Pythia8.
• Cálculo de Emissão: A estrutura calcula fluxos de raios gama de dois canais:
  1. Hadrônico: Decaimento de píons neutros ($\pi^0 \to \gamma\gamma$) de colisões $pp$.
  2. Leptônico: Espalhamento Compton Inverso (IC) de $CRe$ sobre a Radiação Cósmica de Fundo (CMB).
• Integração: Os fluxos são integrados sobre volumes esféricos de 1 kpc a $R_{trunc}$ e comparados com as observações do Fermi-LAT (1–300 GeV) e limites superiores do DAMPE, COMPTEL e INTEGRAL.

Principais Resultados

1. Previsão de Fluxo Hadrônico: O fluxo de raios gama hadrônico previsto na faixa de 1–300 GeV é de $2.82 \times 10^{-11} \text{ cm}^{-2} \text{ s}^{-1}$ dentro de $R_{500}$. Quando integrado até o raio de truncamento ($3.7 R_{500}$), o fluxo sobe para $1.15 \times 10^{-10} \text{ cm}^{-2} \text{ s}^{-1}$. Este valor está em concordância de ordem de grandeza com o fluxo relatado pelo Fermi-LAT de $1.3 \times 10^{-10} \text{ cm}^{-2} \text{ s}^{-1}$.
2. Distribuição Espacial: Aproximadamente 76% do fluxo hadrônico total previsto origina-se de regiões além de $R_{500}$, uma consequência do perfil de densidade externa raso ($\beta \approx 0.33$) que sustenta uma coluna de gás significativa em grandes raios.
3. Subdominância Leptônica: A contribuição IC de elétrons primários e secundários é subdominante ao componente hadrônico por um fator de $\sim20$ dentro de $R_{500}$ na faixa de 1–300 GeV. Mesmo quando integrado até $R_{trunc}$, o fluxo IC ($1.49 \times 10^{-11} \text{ cm}^{-2} \text{ s}^{-1}$) permanece uma ordem de grandeza abaixo do fluxo observado do Fermi-LAT e da previsão hadrônica.
4. Tensão Espectral: Embora a normalização do fluxo em $R_{trunc}$ corresponda às observações, o índice espectral apresenta uma discrepância significativa. O espectro observado é íngreme ($\Gamma = -3.61 \pm 0.33$), enquanto o modelo hadrônico prevê um espectro mais duro ($\Gamma \approx -2.4$ a $-2.6$ para $\alpha_p = 2.4$). Mesmo suavizando o espectro de prótons para $\alpha_p = 2.8$, não é possível reproduzir a inclinação observada.

Significado e Conclusões
O artigo conclui que A3667 é um alvo promissor, mas complexo, para estudos de emissão não térmica. O principal significado do trabalho reside em demonstrar que:

• O nível de fluxo de raios gama observado pode ser reproduzido por um modelo hadrônico padrão se a integração se estender às bordas do aglomerado ($R_{trunc}$), destacando a importância de modelar a emissão além de $R_{500}$.
• O canal IC é negligenciável para explicar o sinal do Fermi-LAT em A3667, confirmando que qualquer sinal detectável é de origem hadrônica sob suposições padrão de RC.
• Existe uma tensão entre a forma espectral observada e as previsões hadrônicas. Os autores afirmam que um modelo ICM hadrônico puro com prótons de RC isotérmicos padrão não pode explicar totalmente a inclinação espectral observada. Eles sugerem que essa discrepância pode indicar contaminação por fontes de fundo ou que o espectro de prótons de RC foi modificado por processos de reaceleração e transporte não capturados por um modelo simples de lei de potência.

Os autores enfatizam que seus resultados são estimativas de ordem de grandeza sujeitas a incertezas na distribuição de RC, estrutura do campo magnético e na suposição de simetria esférica, que não captura a morfologia assimétrica dos choques de fusão e frentes frias de A3667. Trabalhos futuros são sugeridos para incorporar dados de rádio com resolução espacial e modelagem 3D para refinar as previsões de brilho superficial.