An Enhanced Isothermal Jeans Approach to Constraining Dark Matter Self-Interactions from Galactic Kinematics

Resumen Técnico: Un enfoque de Jeans isotérmico mejorado para restringir las autointeracciones de la materia oscura a partir de la cinemática galáctica

Planteamiento del problema
El paradigma cosmológico estándar $\Lambda$CDM enfrenta desafíos en las escalas subgalácticas, particularmente en lo que respecta a la diversidad estructural de las galaxias enanas. Si bien la retroalimentación bariónica puede explicar algunas discrepancias, la Materia Oscura Autointeractuante (SIDM, por sus siglas en inglés) ofrece un mecanismo alternativo donde las partículas de materia oscura experimentan frecuentes auto-dispersiones, impulsando una evolución de dos etapas: la formación inicial de un núcleo seguida de un potencial colapso de núcleo gravotérmico. Los intentos previos de restringir la sección eficaz de auto-dispersión de la SIDM utilizando curvas de rotación galáctica se han visto limitados por los costes computacionales (dependiendo de simulaciones de fluidos 1D o de $N$-cuerpos completas) o por la incapacidad de los modelos semianalíticos para capturar la fase de colapso del núcleo. Además, las implementaciones anteriores del modelo de Jeans isotérmico excluían sistemas con centros dominados por bariones debido a problemas de robustez numérica.

Metodología
Los autores presentan un modelo semianalítico mejorado basado en el enfoque de Jeans isotérmico, construido sobre los marcos establecidos por Jiang et al. (2023) y Yang et al. (2024). La metodología involucra tres mejoras técnicas clave:

1. Secciones eficaces dependientes de la velocidad: El modelo incorpora la dispersión de Rutherford elástica mediante un potencial de Yukawa, reemplazando la suposición de sección eficaz constante por una sección eficaz efectiva ($\sigma_{\text{eff}}$) que depende de la velocidad circular máxima del halo ($V_{\text{max}}$). Esto se parametriza mediante una normalización de baja velocidad ($\sigma_0$) y una escala de velocidad ($\omega$).
2. Tratamiento empírico del colapso del núcleo: Para abordar la ruptura del supuesto isotérmico durante el colapso gravotérmico, los autores adoptan un recurso empírico. Utilizan la solución de alta densidad de las ecuaciones de Jeans, que típicamente se descarta, y la reflejan alrededor del tiempo de fusión ($t_{\text{merge}}$) donde las soluciones de baja y alta densidad convergen. Esto permite al modelo aproximar los perfiles de densidad post-colapso.
3. Robustez numérica mejorada: Los autores refinan el procedimiento de "costura" (empalme del núcleo isotérmico con el halo exterior de tipo CDM) mediante la implementación de umbrales de convergencia estrictos ($\delta^2 = 10^{-5}$) y la partición de los intervalos de búsqueda. Esto permite la identificación fiable tanto de las soluciones de baja densidad (crecimiento del núcleo) como de las de alta densidad (colapso del núcleo), incluso para galaxias con centros dominados por bariones.

El modelo se aplica a una muestra de 68 galaxias de la base de datos SPARC. El proceso de ajuste implica:

• Modelar las distribuciones bariónicas (gas y estrellas) utilizando un perfil de Hernquist, con relaciones masa-luminidad estelar derivadas de ajustes de múltiples perfiles previos.
• Utilizar métodos de Monte Carlo por cadenas de Markov (MCMC) para inferir la masa virial del halo ($M_{\text{vir}}$), la concentración ($c$) y el producto degenerado $p \equiv \sigma_{\text{eff}} t_{\text{age}}$.
• Convertir $p$ en $\sigma_{\text{eff}}$ estimando las edades de los halos mediante la teoría de Press–Schechter extendida.
• Realizar un análisis conjunto para restringir los parámetros globales $\sigma_0$ y $\omega$ minimizando el desacuerdo entre las predicciones del modelo y las curvas de rotación observadas.

Contribuciones clave y resultados

• Expansión de la muestra: La mejora en la estabilidad numérica permite la inclusión de galaxias con centros dominados por bariones, expandiendo significamente el tamaño de la muestra en comparación con estudios previos que utilizan métodos semianalíticos similares.
• Bimodalidad y estados evolutivos: Aproximadamente 1/6 de la muestra (11 galaxias) exhibe posteriors bimodales, donde tanto una solución de crecimiento de núcleo como una solución de colapso de núcleo proporcionan ajustes igualmente buenos. Esta degeneración surge porque los halos en diferentes etapas evolutivas pueden poseer densidades centrales casi idénticas. Las galaxias restantes se encuentran predominantemente en la fase de crecimiento del núcleo o cerca de la etapa de núcleo máximo.
• Restricciones dependientes de la velocidad: Las restricciones conjuntas a través de la muestra revelan una clara dependencia de la velocidad en el espacio de parámetros permitido. Los resultados forman una degeneración en forma de "L" en el plano $(\sigma_0, \omega)$, lo que indica que dos regímenes distintos son viables:
  • Secciones eficaces bajas, casi constantes: $\sigma_0 \sim 2 \, \text{cm}^2/\text{g}$ con $\omega \gtrsim 500 \, \text{km/s}$.
  • Modelos fuertemente dependientes de la velocidad: $\sigma_0 \sim 100 \, \text{cm}^2/\text{g}$ con $\omega \sim 60 \, \text{km/s}$.
  • Adoptando la interpretación de crecimiento del núcleo, los valores de mejor ajuste son $\sigma_0 \simeq 5 \, \text{cm}^2/\text{g}$ y $\omega \simeq 250 \, \text{km/s}$.
• Comparación con CDM: Al compararse con modelos de Materia Oscura Fría (CDM), el modelo SIDM supera a los perfiles simples de contracción adiabática y rivaliza con los perfiles de retroalimentación empírica basados en bariones (específicamente el modelo de Di Cintio et al. 2014) en el ajuste de las curvas de rotación. Sin embargo, a diferencia de los modelos de retroalimentación bariónica, las secciones eficaces de SIDM inferidas no muestran una correlación significativa con la relación masa estelar-halo, lo que sugiere que la SIDM proporciona un mecanismo distinto para la diversidad estructural.
• Relaciones de escala: El estudio encuentra que los modelos de masa de SIDM resultan en concentraciones de halo sistemáticamente más bajas en comparación con las relaciones estándar de CDM (por ejemplo, Dutton & Macciò 2014). Además, aunque algunas galaxias en los ajustes de SIDM exhiben eficiencias de formación estelar altas (acercándose a la unidad), esta tendencia también se observa en estudios previos basados en CDM, lo que indica que la SIDM no resuelve de forma única las tensiones en la relación masa estelar-halo.

Significancia
El artículo sostiene que el paradigma de la SIDM sigue siendo un marco convincente para explicar la diversidad de la estructura a pequeña escala. Las restricciones derivadas son notablemente consistentes con sondas independientes, incluyendo la cinemática estelar en enanas ultra débiles, el lente fuerte y las funciones de masa de agujeros negros supermasivos en el universo temprano (específicamente los "Little Red Dots"). Los autores enfatizan que, si bien la degeneración en forma de L impide un modelo de sección eficaz único que se ajuste a toda la muestra de galaxias, los resultados resaltan una región específica del espacio de parámetros que satisface simultáneamente las restricciones ortogonales de la cinemática de galaxias locales y la demografía de agujeros negros de alto desplazamiento al rojo. El estudio concluye que los futuros datos cinemáticos de baja masa, particularmente para enanas ultra débiles ($V_{\text{max}} \lesssim 30 \, \text{km/s}$), serán cruciales para romper las degeneraciones actuales en los modelos de sección eficaz dependiente de la velocidad.