Cosmological Zoom-In Simulation of Odd Radio Circles as Merger-Driven Shocks in Galaxy Groups

Problema
Os Círculos de Rádio Ímpares (ORCs) representam uma classe recém-descoberta de objetos de rádio distintos, caracterizados por morfologias de brilho de borda em forma de anel ou casca, frequentemente envolvendo galáxias do tipo inicial em altas latitudes galácticas. Sua origem permanece incerta, pois suas propriedades peculiares — especificamente seus tamanhos enormes (frequentendo excedendo o raio de virial de seus sistemas hospedeiros) e altos poderes de rádio — desafiam os entendimentos atuais de fontes astrofísicas para emissão de rádio difusa. Embora vários cenários de formação tenham sido propostos, incluindo choques de terminação de fluxos galácticos (outflows), fusões de buracos negros supermassivos binários e jatos de AGN, o canal de formação impulsionado por choque via fusões de galáxias majoritárias foi sugerido, mas não totalmente validado no contexto da massa de halo necessária. Simulações anteriores (ex: Dolag et al. 2023) tentando modelar ORCs como choques impulsionados por fusões em grupos de galáxias previram luminosidades de rádio significativamente menores do que as observadas, potencialmente devido ao uso de halos de menor massa ($10^{12} M_\odot$) e tratamentos de pós-processamento de raios cósmicos (CRs) e campos magnéticos.

Metodologia
Para abordar as limitações de estudos anteriores, os autores realizaram uma simulação de "zoom-in" cosmológica de alta resolução de uma fusão de grupo de galáxias massivo usando o código OpenGadget3. As principais características metodológicas incluem:

• Sistema Alvo: Um grupo de galáxias com uma massa de virial final de $M_{200} \sim 1.8 \times 10^{13} M_\odot$, consistente com as massas estimadas das galáxias hospedeiras de ORCs.
• Física: Uma simulação magnetohidrodinâmica (MHD) não radiativa empregando Magnetohidrodinâmica de Partículas Suavizadas (SPMHD) e um solver de gravidade Tree-PM.
• Raios Cósmicos: Um tratamento on-the-fly de elétrons de CR espectrais usando o módulo CRESCENDO (um solver de Fokker-Planck). A injeção de CR é impulsionada exclusivamente pela Aceleração de Choque Difusiva (DSA) baseada em um localizador de choque on-the-fly que determina as propriedades do choque (número de Mach, obliquidade).
• Campos Magnéticos: A simulação resolve equações MHD no regime não ideal. No entanto, reconhecendo que a simulação não pode resolver totalmente o dínamo turbulento necessário para uma amplificação eficiente do campo magnético, os autores também testaram modelos analíticos para a força do campo magnético (equilíbrio de pressão térmica e equilíbrio de pressão turbulenta) para estimar luminosidades de rádio realistas.
• Observáveis: O estudo deriva o brilho de superfície de raios X e a emissão de rádio sincrotrão de forma autoconsistente. Observações simuladas foram geradas em frequências correspondentes ao LOFAR, GMRT e ASKAP, incluindo mapas de índice espectral e frações de polarização (parâmetros de Stokes), contabilizando a rotação de Faraday.

Principais Contribuições

1. Primeira Simulação MHD de Alta Massa de ORCs: Este trabalho apresenta a primeira simulação de zoom-in de uma fusão de grupo de galáxias ($M \sim 10^{13} M_\odot$) com tratamento de CR e campo magnético consistentes e on-the-fly, visando especificamente a formação de ORCs.
2. Validação da Morfologia de Choque de Fusão: A simulação demonstra que uma fusão majoritária entre duas galáxias desencadeia choques fortes no meio circumgaláctico (CGM). Quando visualizados perpendicularmente ao eixo de fusão, esses choques formam uma estrutura de anel gigante em rápida expansão que coincide com a morfologia e escala (centenas de kpc a $\sim 1$ Mpc) dos ORCs observados.
3. Consistência de Raios X: O sistema simulado reproduz as propriedades de raios X de grupos de galáxias observados nesta faixa de massa. Notavelmente, ele explica com sucesso o deslocamento (offset) do ORC Physalis, onde a emissão de raios X está centrada em apenas uma das galáxias em fusão, mostrando que a remoção (stripping) de gás deixa uma galáxia com um núcleo significativamente mais denso que a outra.
4. Análise de Polarização: Pela primeira vez, o artigo fornece uma previsão autoconsistente da fração de polarização para um ORC impulsionado por fusão, comparando-a diretamente com dados observacionais.

Resultos

• Morfologia e Evolução: Pouco após a fusão (aprox. 70 Myr), duas grandes cascas de choque emergem, coalescendo em um anel com diâmetro de 300–400 kpc, que posteriormente se expande para $\sim 1$ Mpc, ultrapassando o raio de virial. A frente de choque é claramente traçada pela densidade de energia de elétrons de CR e quedas no brilho de superfície de raios X.
• Luminosidade de Rádio: Apesar do aumento da massa do halo e do tratamento de CR autoconsistente, a luminosidade de rádio simulada permanece aproximadamente três ordens de magnitude abaixo dos ORCs observados (ex: ORC Physalis). Mesmo aplicando aumentos analíticos na força do campo magnético (assumindo equipartição com a pressão térmica ou turbulenta), a luminosidade prevista falha em atingir os níveis observados.
• Propriedades Espectrais: O espectro de emissão sugere um número de Mach inferido por rádio de $\mathcal{M} \approx 4.8$ na frente de choque, consistente com choques fortes em aglomerados de galáxias. O índice espectral torna-se mais íngreme com o tempo e com a frequência devido a perdas radiativas e diminuição da força do choque.
• Polarização: A simulação prevê uma fração de polarização de 20–40% através do anel. Isso é sistematicamente superior aos valores observados para o ORC 1 e ORC J0356–4216, que tipicamente mostram 10–20% no interior e até $\sim 30\%$ no anel.
• Identificação de Choque: O estudo confirma que o localizador de choque na simulação resolve choques em escalas de $\sim 10$ kpc, enquanto a filtragem Gaussiana de mapas de pressão revela as bordas de pressão circulares nítidas responsáveis pelas características de rádio.

Significância e Alegações
O artigo conclui que fusões majoritárias em grupos de galáxias são um canal de formação viável para a morfologia e a grande extensão dos ORCs. A simulação reproduz com sucesso a estrutura em forma de anel e as características de raios X do ambiente hospedeiro. No entanto, os autores afirmam explicitamente que o mecanismo de aceleração de choque difusiva (DSA) sozinho é insuficiente para explicar as propriedades totais de rádio dos ORCs.

A discrepância na luminosidade de rádio e a sobreprevisão da polarização sugerem que a simulação carece de componentes físicos críticos. Os autores propõem que populações fósseis de raios cósmicos (injetadas por eventos prévios de Núcleos Galácticos Ativos ou feedback estelar) são provavelmente necessárias para explicar o poder de rádio observado. Esses CRs pré-existentes poderiam ser re-acelerados pelo choque de fusão, aumentando o brilho de rádio e aumentando o número de centros de emissão ao longo da linha de visão, reduzendo assim a fração de polarização líquida para coincidir com as observações. O artigo afirma que trabalhos futuros devem incorporar formação estelar, feedback de AGN e tratamento de CR fósseis autoconsistente para resolver plenamente a origem dos ORCs.