Gravothermal Collapse: Robust Against Baryonic Feedback

Resumen Técnico: Colapso Gravotérmico en Halos de SIDM bajo Retroalimentación Bariónica

Planteamiento del Problema
El modelo estándar de Materia Oscura Fría (CDM) predice cúspides de densidad central pronunciadas en los halos de materia oscura, lo que a menudo entra en conflicto con los núcleos poco densos observados en galaxias enanas. Aunque los mecanismos de retroalimentación bariónica (por ejemplo, flujos de salida impulsados por supernovas) pueden inyectar energía para transformar cúspides en núcleos, observaciones recientes revelan una amplia diversidad en las distribuciones de materia oscura —que van desde núcleos poco densos hasta cúspides extremadamente densas— que el modelo CDM tiene dificultades para reproducir. Además, el lente gravitacional fuerte y las observaciones de corrientes estelares han identificado perturbadores compactos con densidades significativamente superiores a las predicciones típicas de CDM. En el marco de la Materia Oscura Autointeractuante (SIDM), las autointeracciones de la materia oscura impulsan la evolución gravotérmica, produciendo naturalmente tanto núcleos de baja densidad como estados colapsados de alta densidad. Sin embargo, estudios previos se han basado mayoritariamente en simulaciones de solo SIDM o potenciales bariónicos estáticos; dado que los potenciales bariónicos realistas moldeados por la retroalimentación son dinámicos y no adiabáticos, sigue sin estar claro si la predicción definitoria de SIDM —el colapso gravotérmico— sigue siendo robusta frente a la inyección de energía y las fluctuaciones de potencial causadas por una fuerte retroalimentación bariónica.

Metodología
Los autores realizan simulaciones controladas de $N$-cuerpos utilizando el código GADGET-2 con un módulo de SIDM para poner a prueba el colapso gravotérmico bajo una fuerte retroalimentación bariónica.

• Condiciones Iniciales: Se simulan dos halos de materia oscura con masa $M_{200} = 2 \times 10^9 M_\odot$: un halo de alta concentración ($c_{200} = 42$, que representa una fluctuación de $3\sigma$) y un halo de concentración mediana ($c_{200} = 15$, la mediana cosmológica). Ambos siguen perfiles de Navarro–Frenk–White (NFW).
• Parámetros de SIDM: La sección eficaz de autointeracción por unidad de masa se fija en $\sigma/m = 50 \text{ cm}^2/\text{g}$, consistente con modelos que explican los perturbadores de lente fuerte de alta densidad y la galaxia enana difusa Crater II.
• Modelo de Retroalimentación: Se implementa un potencial oscilante semianalítico para imitar la retroalimentación bariónica. Se añade un potencial de Plummer central con una masa dependiente del tiempo $M_{ext}(t)$ al cálculo de la fuerza. El modelo incluye "Retroalimentación Fuerte" (donde $M_{ext}$ oscila entre valores positivos y negativos, creando aceleración repulsiva) y "Retroalimentación Débil" (que oscila entre cero y valores positivos). El periodo de oscilación se establece en $P = 0.2$ Gyr, con una amplitud de masa máxima $M_{max} = 10^7 M_\odot$ y un radio de escala $a = 0.1$ kpc.
• Escenarios: El estudio compara la evolución de solo SIDM frente a casos con retroalimentación continua fuerte/débil, y retroalimentación episódica (activa durante intervalos específicos, por ejemplo, 0–4 Gyr u 8–10 Gyr) seguida ya sea de su eliminación o de su reemplazo por un potencial estático.

Resultados Clave
Las simulaciones revelan que la resiliencia del colapso gravotérmico depende críticamente de la competencia entre la escala de tiempo de termalización de SIDM ($t_{th}$) y el periodo de oscilación de la retroalimentación ($P$).

1. Halos de Alta Concentración ($c_{200} = 42$):

   • La escala de tiempo de termalización es corta ($t_{th} \sim 0.01\text{--}0.04$ Gyr), significativamente más corta que el periodo de retroalimentación ($P = 0.2$ Gyr).
   • En consecuencia, la dinámica de SIDM domina sobre la retroalimentación bariónica.
   • Bajo retroalimentación fuerte, el colapso solo se retrasa levemente (por $\sim 1$ Gyr) y nunca se detiene. La densidad central alcanza eventualmente niveles consistentes con las simulaciones de solo SIDM y las restricciones observacionales para perturbadores densos (por ejemplo, JVAS B1938+666).
   • La retroalimentación débil tiene un efecto insignificante o acelera ligeramente el colapso debido a la contracción durante las fases de entrada de flujo.

2. Halos de Concentración Mediana ($c_{200} = 15$):

   • La escala de tiempo de termalización es más larga ($t_{th} > P$ durante la mayor parte de la evolución), lo que hace al halo susceptible al calentamiento impulsado por la retroalimentación.
   • Retroalimentación Fuerte: Durante la fase de expansión del núcleo, la retroalimentación fuerte puede generar rápidamente núcleos grandes (hasta $\sim 0.8$ kpc) y reducir la densidad central a $\sim 25\%$ de su valor inicial. Sin embargo, una vez que la retroalimentación cesa, el colapso gravotérmico se reanuda, aunque a un ritmo más lento.
   • Retroalimentación Episódica: El perfil de densidad final es altamente sensible a la historia de la retroalimentación.
     • Si la retroalimentación ocurre temprano (0–4 Gyr) y se elimina, el halo colapsa más tarde que en el caso de solo SIDM.
     • Si el potencial oscilante es reemplazado por un potencial estático después de la retroalimentación, el colapso se acelera.
     • Si la retroalimentación ocurre tarde (8–10 Gyr) durante la fase de colapso, la detiene temporalmente pero no revierte el proceso.
   • Esta sensibilidad produce una amplia diversidad en las densidades centrales finales, que van desde núcleos poco densos hasta cúspides densas, dependiendo de la línea de tiempo específica de la retroalimentación.

Significancia y Reivindicaciones
El artículo establece que el colapso gravotérmico es una predicción robusta de SIDM, resiliente incluso ante una retroalimentación bariónica fuerte y oscilante.

• Para Halos de Alta Concentración: La corta escala de tiempo de termalización asegura que el colapso proceda independientemente de la intensidad de la retroalimentación. Esto refuerza la interpretación de los perturbadores de lente fuerte densos y los objetos compactos como halos de SIDM en colapso de núcleo. También apoya los escenarios de SIDM para la formación de semillas de agujeros negros supermasivos en el Universo temprano, incluso en entornos con formación estelar intermitente.
• Para Halos de Concentración Mediana: La interacción entre la retroalimentación y las autointeracciones puede explicar la diversidad observada en las distribuciones de materia oscura (desde núcleos poco densos hasta cúspides densas) dentro de halos de masa similar. El estado final no está determinado únicamente por la masa del halo o la concentración inicial, sino que es modulado por la historia episódica de la retroalimentación bariónica.
• Conclusión: Los resultados sugieren que el colapso gravotérmico sigue siendo una "huella distintiva" (smoking-gun) de SIDM. La diversidad predicha en los perfiles de densidad, impulsada por la combinación de autointeracciones y la historia de la retroalimentación, ofrece un marco comprobable para futuras observaciones de curvas de rotación de galaxias y la naturaleza de la materia oscura.