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

Problema
Los Círculos de Radio Impares (ORCs, por sus siglas en inglés) representan una clase de objetos de radio distintos recientemente descubiertos, caracterizados por morfologías de borde brillante en forma de anillo o cáscara, que a menudo rodean a galaxias tempranas centrales en altas latitudes galácticas. Su origen sigue siendo incierto, ya que sus propiedades peculiares —específicamente sus enormes tamaños (que a menudo superan el radio virial de sus sistemas anfitriones) y su alta potencia de radio— desafían las comprensiones actuales de las fuentes astrofísicas para la emisión de radio difusa. Si bien se han propuesto varios escenarios de formación, incluyendo choques de terminación de flujos galácticos, fusiones de agujeros negros supermasivos gigantes y chorros de AGN, el canal de formación impulsado por choques mediante fusiones de galaxias mayores ha sido sugerido pero no validado completamente en el contexto de las masas de halo requeridas. Simulaciones previas (por ejemplo, Dolag et al. 2023) que intentaban modelar los ORCs como choques impulsados por fusiones en grupos de galaxias predijeron luminosidades de radio significativamente menores a las observadas, potencialmente debido al uso de halos de menor masa ($10^{12} M_\odot$) y tratamientos de post-procesamiento de rayos cósmicos (CRs) y campos magnéticos.

Metodología
Para abordar las limitaciones de estudios previos, los autores realizaron una simulación de "zoom-in" cosmológica de alta resolución de una fusión de un grupo de galaxias masivo utilizando el código OpenGadget3. Las características metodológicas clave incluyen:

• Sistema Objetivo: Un grupo de galaxias con una masa virial final de $M_{200} \sim 1.8 \times 10^{13} M_\odot$, consistente con las masas estimadas de las galaxias anfitrionas de los ORCs.
• Física: Una simulación magnetohidrodinámica (MHD) no radiativa empleando Magnetohidrodinámica de Partículas Suavizadas (SPMHD) y un resolvedor de gravedad Tree-PM.
• Rayos Cósmicos: Un tratamiento "on-the-fly" de electrones de CR espectrales utilizando el módulo CRESCENDO (un resolvedor de Fokker-Planck). La inyección de CR es impulsada únicamente por la Aceleración de Choque Difusiva (DSA) basada en un buscador de choques "on-the-fly" que determina las propiedades del choque (número de Mach, oblicuidad).
• Campos Magnéticos: La simulación resuelve las ecuaciones MHD en el régimen no ideal. Sin embargo, reconociendo que la simulación no puede resolver completamente el dínamo turbulento requerido para una amplificación eficiente del campo magnético, los autores también probaron modelos analíticos para la fuerza del campo magnético (equilibrio de presión térmica y equilibrio de presión turbulenta) para estimar luminosidades de radio realistas.
• Observables: El estudio deriva la brillantez superficial de rayos X y la emisión de radio sincrotrón de manera autosuficiente. Se generaron observaciones simuladas a frecuencias correspondientes a LOFAR, GMRT y ASKAP, incluyendo mapas de índice espectral y fracciones de polarización (parámetros de Stokes), contabilizando la rotación de Faraday.

Contribuciones Clave

1. Primera Simulación MHD de Alta Masa de ORCs: Este trabajo presenta la primera simulación de "zoom-in" de una fusión de grupo de galaxias ($M \sim 10^{13} M_\odot$) con un tratamiento autosuficiente y "on-the-fly" de CR y campos magnéticos, dirigida específicamente a la formación de ORCs.
2. Validación de la Morfología de Choque de Fusión: La simulación demuestra que una fusión mayor entre dos galaxias desencadena choques fuertes en el medio circumgaláctico (CGM). Cuando se observa perpendicularmente al eje de fusión, estos choques forman una estructura de anillo gigante que se expande rápidamente, la cual coincide con la morfología y escala (cientos de kpc a $\sim 1$ Mpc) de los ORCs observados.
3. Consistencia en Rayos X: El sistema simulado reproduce las propiedades de rayos X de los grupos de galaxias observados en este rango de masa. Notablemente, explica con éxito el desfase observado en el ORC Physalis, donde la emisión de rayos X se centra en solo una de las galaxias en fusión, al mostrar que el barrido de gas deja a una galaxia con un núcleo significativamente más denso que la otra.
4. Análisis de Polarización: Por primera vez, el artículo proporciona una predicción autosuficiente de la fracción de polarización para un ORC impulsado por fusión, comparándola directamente con datos observacionales.

Resultos

• Morfología y Evolución: Poco después de la fusión (aprox. 70 Myr), emergen dos grandes cáscaras de choque, que se fusionan en un anillo con un diámetro de 300–400 kpc, el cual posteriormente se expande a $\sim 1$ Mpc, superando el radio virial. El frente de choque es trazado claramente por la densidad de energía de los electrones de CR y las caídas en la brillantez superficial de rayos X.
• Luminosidad de Radio: A pesar del aumento de la masa del halo y el tratamiento autosuficiente de CR, la luminosidad de radio simulada permanece aproximadamente tres órdenes de magnitud por debajo de los ORCs observados (por ejemplo, ORC Physalis). Incluso aplicando aumentos analíticos a la fuerza del campo magnético (asumiendo equipartición con la presión térmica o turbulenta), la luminosidad predicha no logra alcanzar los niveles observados.
• Propiedades Espectrales: El espectro de emisión sugiere un número de Mach inferido por radio de $\mathcal{M} \approx 4.8$ en el frente de choque, consistente con choques fuertes en cúmulos de galaxias. El índice espectral se vuelve más pronunciado con el tiempo y con la frecuencia debido a las pérdidas radiativas y la disminución de la fuerza del choque.
• Polarización: La simulación predice una fracción de polarización del 20–40% a través del anillo. Esto es sistemáticamente superior a los valores observados para ORC 1 y ORC J0356–4216, que típicamente muestran un 10–20% en el interior y hasta un $\sim 30\%$ en el anillo.
• Identificación de Choques: El estudio confirma que el buscador de choques en la simulación resuelve choques en escalas de $\sim 10$ kpc, mientras que el filtrado Gaussiano de mapas de presión revela los bordes de presión circulares y nítidos responsables de las características de radio.

Significancia y Reivindicaciones
El artículo concluye que las fusiones mayores en grupos de galaxias son un canal de formación factible para la morfología y la gran extensión de los ORCs. La simulación reproduce con éxito la estructura de tipo anillo y las características de rayos X del entorno anfitrión. Sin embargo, los autores declaran explícitamente que el mecanismo de aceleración de choque difusiva (DSA) por sí solo es insuficiente para explicar las propiedades de radio completas de los ORCs.

La discrepancia en la luminosidad de radio y la sobrepredicción de la polarización sugieren que la simulación carece de componentes físicos críticos. Los autores proponen que las poblaciones de rayos cósmicos fósiles (inyectados por eventos previos de núcleos galácticos activos o retroalimentación estelar) son probablemente necesarias para explicar la potencia de radio observada. Estos CRs preexistentes podrían ser re-acelerados por el choque de la fusión, aumentando la brillantez de radio e incrementando el número de centros de emisión a lo largo de la línea de visión, reduciendo así la fracción de polarización neta para coincidir con las observaciones. El artículo afirma que el trabajo futuro debe incorporar la formación estelar, la retroalimentación de AGN y el tratamiento autosuficiente de CR fósiles para resolver plenamente el origen de los ORCs.