ICM-SHOX III. The case of MACS J0018.5+1626, a radio relic that looks like a radio halo?

Resumo Técnico: ICM-SHOX III. O caso de MACS J0018.5+1626, um relicado de rádio que se parece com um halo de rádio?

Definição do Problema
A emissão de rádio difusa em aglomerados de galáxias, tipicamente classificada como halos de rádio (morfologia central, arredondada) ou relicados de rádio (periféricos, alongados, impulsionados por choques), fornece informações críticas sobre o orçamento de energia não térmica e a dinâmica de fusão do meio intracluster (ICM). No entanto, a interpretação de tal emissão é frequentemente complicada por efeitos de projeção e pela falta de restrições multi-sonda sobre os parâmetros da fusão (massa, parâmetro de impacto, ângulo de visão, estágio da fusão). O aglomerado de galáxias MACS J0018.5+1626 apresenta um desafio específico: observações do LOFAR a 144 MHz revelam uma fonte de rádio difusa, centralmente localizada, com uma morfologia típica de um halo de rádio. Por outro modo, o pipeline Improved Constraints on Mergers with SZ, Hydrodynamical simulations, Optical, and X-ray (ICM-SHOX) sugere que o aglomerado está passando por uma fusão binária próxima à passagem do pericentro, visualizada quase ao longo do eixo de fusão. Esta configuração teoricamente produz dois choques impulsionados pela fusão que, se visualizados de frente (face-on), poderiam mimetizar um halo de rádio central. O problema central abordado é se a aceleração de choque difusiva (DSA) padrão atuando sobre elétrons térmicos em tal configuração de choque face-on pode reproduzir a potência e a morfologia de rádio observadas em MACS J0018.5+1626, ou se mecanismos alternativos (por exemplo, re-aceleração turbulenta ou elétrons fósseis) são necessários.

Metodologia
Os autores utilizam um arcabouço numérico híbrido combinando simulações magnetohidrodinâmicas (MHD) 3D com um resolvedor de Fokker–Planck para modelar a aceleração de partículas e a emissão de sincrotron.

• Simulações MHD: Usando o código de malha móvel AREPO, os autores simulam fusões binárias idealizadas de aglomerados de galáxias. As condições iniciais são restringidas pelos resultados do pipeline ICM-SHOX para o MACS J0018.5+1626, especificamente definindo o redshift ($z \approx 0.5456$), massas dos componentes e perfis de gás. As simulações exploram variações nas velocidades relativas iniciais ($v_i = 2400$ e $3000$ km s$^{-1}$), parâmetros de impacto ($b = 100$ e $250$ kpc) e beta de plasma inicial ($\beta_p = 100$ e $200$).
• Rastreamento de Partículas: Partículas traçadoras são advectadas com o fluxo do fluido para rastrear a evolução lagrangiana do gás. Apenas traçadores que encontram células chocadas (número de Mach $M > 1.3$) são considerados para a aceleração de partículas.
• Resolvedor de Fokker–Planck: Um resolvedor de Fokker–Planck modela a evolução da distribuição de momento dos elétrons de raios cósmicos (CRe). O resolvedor contabiliza a aceleração de choque difusiva (DSA), expansão adiabática e perdas radiativas (sincrotron, Compton inverso, Coulomb, Bremsstrahlung).
• Modelos de Aceleração de Choque: O estudo testa dois modelos de eficiência de aceleração de choque: uma eficiência constante ($\eta_0$) e uma eficiência dependente do número de Mach baseada em Kang et al. (2007). A eficiência de injeção de elétrons é tratada como um parâmetro livre ($\eta_e$) escalonado pela eficiência de prótons, explorando valores de $10^{-3}$ a $10^{-1}$.
• Emissão de Sincrotron: O resolvedor computa a emissividade de sincrotron para os traçadores, gerando mapas 3D projetados ao longo de várias linhas de visada (LOS). Estes mapas são convoluídos com os tamanhos de feixe do LOFAR para comparação direta com observações a 144 MHz.
• Comparação Observacional: Mapas simulados de rádio e raios-X são comparados contra dados do LOFAR (144 MHz) e Chandra (0.5–7 keV). As métricas incluem potência de rádio total, perfis de brilho de superfície radial, índices espectrais e coeficientes de correlação ponto-a-ponto entre rádio e raios-X.

Resultados Principais

1. Configuração de Fusão e Propriedades de Choque: As simulações confirmam que uma fusão visualizada quase ao longo do eixo (frente/ face-on) produz dois choques que se sobrepõem significativamente no centro do aglomerado. A fusão resultante gera choques com números de Mach médios típicos de $M_s \sim 2–3$ (com $\sigma_M \sim 0.5–1.5$).
2. Reprodução da Potência de Rádio: A DSA padrão pode reproduzir a potência de rádio observada do LOFAR ($\sim 37$ mJy a 144 MHz) apenas sob condições específicas:
   • Eficiência: É necessária uma eficiência de aceleração de choque de elétrons de $\eta_e \sim 10^{-3}–10^{-2}$. Isto corresponde ao limite otimista das estimativas teóricas, mas é consistente com alguns resultados de simulações de Particle-in-Cell (PIC) para choques fracos.
   • Campo Magnético: Um beta de plasma inicial de $\beta_p = 100$ é suficiente para gerar a amplificação do campo magnético necessária ($\sim 0.8–1$ $\mu$G nas regiões de choque) para corresponder às observações. Modelos com $\beta_p = 200$ falham em atingir a potência observada.
   • Época: O melhor ajuste ocorre pouco depois da passagem do pericentro (aproximadamente 40–200 Myr pós-pericentro), coincidindo com o pico de dissipação de energia do choque.
3. Morfologia e Ângulo de Visada: A projeção face-on dos dois choques sobrepostos produz naturalmente uma morfologia de rádio circular e centralizada, que se assemelha fortemente a um halo de rádio. O estudo descobre que ângulos de visão $\lesssim 30^\circ$ (para $v_i=2400$ km s$^{-1}$) ou $\lesssim 36^\circ$ (para $v_i=3000$ km s$^{-1}$) fornecem o melhor ajuste aos perfis radiais e à distribuição de brilho de superfície. Ângulos de visão maiores resultam em estruturas tipo relicado, não sobrepostas e inconsistentes com as observações.
4. Propriedades Espectrais e de Correlação:
   • O índice espectral integrado ($\alpha_{int}$) varia de $\sim -1.3$ a $-1.8$, consistentes com valores típicos de halos e relicados de rádio.
   • A correlação ponto-a-ponto entre o brilho de superfície de rádio e raios-X é encontrada como variável. Dependendo do estágio da fusão, do ângulo de visão e da resolução, a correlação pode ser superlinear, sublinear ou linear. Isto desafia a noção de que correlações sublineares são uma assinatura única de halos de rádio (origem turbulenta) versus relicados (origem por choque).

Significância e Alegações
O artigo afirma que o aparente "halo de rádio" em MACS J0018.5+1626 é provavelmente o resultado da superposição de dois relicados de rádio impulsionados por fusão visualizados quase de frente (face-on), em vez de um halo turbulento centralizado. Esta interpretação é sustentada pelas restrições multi-comprimento de onda do pipeline ICM-SHOX, que favorecem uma geometria de fusão que se alinha com o cenário de choque face-on.

O estudo demonstra que a DSA padrão, quando aplicada a uma configuração de choque face-on com eficiências de aceleração otimistas mas fisicamente plausíveis ($\eta_e \sim 10^{-3}–10^{-2}$) e forças de campo magnético, pode reproduzir a potência e a morfologia de rádio do aglomerado sem invocar re-aceleração turbulenta ou populações de elétrons fósseis. No entanto, os autores observam que as eficiências necessárias estão no lado otimista, sugerindo que mecanismos adicionais (como re-aceleração turbulenta ou elétrons fósseis) ainda podem desempenhar um papel.

O trabalho destaca a importância crítica do ângulo de visão e dos efeitos de projeção na classificação da emissão de rádio difusa. Estabelece também o pipeline ICM-SHOX como uma ferramenta poderosa para ligar observáveis térmicos (raios-X, SZ) e não térmicos (rádio) para restringir o estado físico de aglomerados em fusão. Os autores concluem que, embora o seu modelo explique com sucesso as observações sob a DSA padrão, observações de rádio de alta frequência e dados de polarização são necessários para distinguir ainda mais as origens por choque e turbulenta, e para restringir melhor a eficiência de aceleração de elétrons no ICM.