A Constraint on Dark Matter Self-Interaction from Combined Strong Lensing and Stellar Kinematics in MACS J0138-2155

Resumen Técnico: Restricción a la Auto-Interacción de la Materia Oscura en MACS J0138-2155

1. El Problema Científico

La naturaleza de la materia oscura (DM) sigue siendo uno de los mayores enigmas de la física. El modelo estándar de Materia Oscura Fría (CDM) asume que las partículas de DM son no colisionales. Si bien este modelo explica con éxito la estructura a gran escala del universo, enfrenta tensiones en escalas más pequeñas (como galaxias enanas), donde las simulaciones de CDM predicen perfiles de densidad "cuspidales" (muy concentrados en el centro) que a veces no coinciden con las observaciones.

Una alternativa prometedora es la Materia Oscura Auto-Interactuante (SIDM), donde las partículas de DM pueden dispersarse entre sí mediante una nueva fuerza. Esto puede transformar los perfiles de densidad cuspidales en "núcleos" (cores) más planos o, en escalas de tiempo largas, provocar un colapso gravotermal. Sin embargo, las restricciones existentes sobre la sección transversal de interacción ($\sigma/m$) varían significativamente según la escala de masa (desde galaxias enanas hasta cúmulos de galaxias), sugiriendo que la interacción podría depender de la velocidad.

El objetivo de este trabajo es proporcionar una de las restricciones más precisas a la fecha sobre la SIDM en la escala de cúmulos de galaxias, utilizando el sistema MACS J0138-2155.

2. Metodología

Los autores desarrollan un modelo de masa autoconsistente que combina dos sondas gravitacionales complementarias para restringir el perfil de densidad de la materia oscura en un rango de radios de $\sim 1$ a $\sim 100$ kpc.

• Sistema Objetivo: MACS J0138-2155, un cúmulo de galaxias relajado y no en fusión, que actúa como una lente gravitacional fuerte. Es único por albergar dos supernovas fuertemente lentejadas: SN Requiem y SN Encore.
• Datos Utilizados:
  • Imágenes de Alta Resolución: Datos de HST (programas 14496 y 15663) y JWST (NIRCam, programa GO-6549) para mapear las posiciones de las imágenes lentejadas.
  • Espectroscopía Integral de Campo (IFU): Datos del instrumento MUSE en el VLT, que permiten medir la cinemática estelar resuelta espacialmente de la galaxia central (BCG).
• Enfoque de Modelado:
  • Perfil SIDM: Utilizan el formalismo de Jeans isotérmico para modelar el efecto de la auto-interacción. En lugar de ajustar un perfil paramétrico y luego relacionarlo con SIDM, calculan el perfil de densidad predicho por SIDM en cada paso del modelo. Esto implica un núcleo isotérmico dentro de un radio $r_1$ y un perfil NFW fuera de él.
  • Lente Gravitacional Fuerte: Utilizan el código PyAutoLens para modelar las posiciones de las imágenes múltiples de la fuente (galaxia en reposo a $z=1.95$) y las supernovas. Incluyen halos de materia oscura para el cúmulo, galaxias miembros (escaladas mediante la relación Faber-Jackson) y cizallamiento externo.
  • Cinemática Estelar: Utilizan el paquete JAM (Jeans Anisotropic MGE) para modelar la dispersión de velocidades estelares de la BCG. La luz de la galaxia se modela con una Expansión Multi-Gaussiana (MGE) derivada de imágenes de HST.
  • Geometría 3D: A diferencia de estudios anteriores que asumen simetría elíptica en el plano del cielo, este modelo permite un halo prolato tridimensional, proyectando correctamente la densidad 3D a la densidad superficial 2D observada.
• Análisis Estadístico: Se emplea un muestreador MCMC (emcee) para obtener distribuciones posteriores de los parámetros, incluyendo la masa del halo ($M_{200m}$), la concentración ($c_{200m}$), la relación masa-luz estelar ($\Upsilon_*$) y la sección transversal de SIDM ($\sigma/m$).

3. Contribuciones Clave

1. Integración de Sondas Complementarias: Es el primer estudio que combina de manera autoconsistente lentes gravitacionales fuertes y cinemática estelar resuelta espacialmente en un solo sistema de cúmulo para restringir la SIDM. Esto permite romper la degeneración masa-hoja y restringir la densidad en el interior del halo donde la SIDM tiene mayor efecto.
2. Modelado de SIDM "Desde Adentro": Implementan un método donde el perfil de densidad modificado por SIDM se calcula directamente a partir de las ecuaciones de Jeans en cada iteración, en lugar de usar aproximaciones paramétricas post-hoc.
3. Tratamiento de Geometría 3D: Introducen un modelo de halo prolato tridimensional, abordando el sesgo de concentración y triaxialidad inherente a los cúmulos seleccionados por lentes.
4. Predicción de Supernovas: Proporcionan predicciones actualizadas sobre los tiempos de retardo y magnificaciones de las supernovas lentejadas, cruciales para la cosmografía de tiempo de retardo.

4. Resultados Principales

El análisis se realizó en tres configuraciones: solo lente, solo cinemática y modelo conjunto.

• Restricción a SIDM:
  • La combinación de datos de lente y cinemática produce un límite superior de confianza del 95% sobre la sección transversal de interacción:
    $$\sigma/m < 0.613 \, \text{cm}^2/\text{g}$$
    para una velocidad de interacción promedio $\langle v_{pair} \rangle < 2090 \, \text{km/s}$.
  • Este resultado es competitivo con los límites anteriores de cúmulos, pero derivado de un análisis de sistema único mucho más detallado.
• Consistencia de Modelos:
  • Los resultados individuales de lente y cinemática son altamente consistentes entre sí.
  • El modelo conjunto revela una bimodalidad en los parámetros de $\Upsilon_*$ (relación masa-luz) y $\sigma/m$. Existe una degeneración entre una alta masa estelar con un perfil de DM cuspide y una baja masa estelar con un perfil de DM con núcleo (SIDM). Sin embargo, la inclusión de cinemática suprime la preferencia por valores altos de $\sigma/m$ observada en el modelo solo de lente.
• Propiedades del Cúmulo:
  • Se mide una concentración del halo ($c_{200m} \approx 10$) que es aproximadamente $2\sigma$ más alta que la predicha por la relación masa-concentración estándar de $\Lambda$CDM. Esto se atribuye al sesgo de selección de lentes (halos más concentrados y elongados a lo largo de la línea de visión son más propensos a lensar) y al estado relajado del cúmulo.
  • Se confirma la presencia de un halo de materia oscura asociado a una galaxia roja brillante cercana al cúmulo, lo que mejora el ajuste del modelo.
• Predicciones de Supernovas:
  • El modelo predice que la siguiente aparición de SN Requiem ocurrirá entre enero de 2027 y noviembre de 2028 (con un retraso de ~12 años respecto a la primera aparición), mucho antes de lo predicho por modelos anteriores (que sugerían mediados de la década de 2030).
  • Para SN Encore, se predice una reaparición alrededor de 2033-2034.

5. Significado e Impacto

• Avance en la Física de la Materia Oscura: Este estudio demuestra que los cúmulos de galaxias pueden imponer restricciones estrictas a la SIDM, complementando las restricciones de galaxias enanas. La consistencia entre las sondas valida el uso de la cinemática estelar para refinar modelos de lentes.
• Cosmografía de Tiempo de Retardo: La predicción de la reaparición de SN Requiem en 2027-2028 es crucial. Observar este evento permitirá medir la constante de Hubble ($H_0$) con una precisión sin precedentes desde un solo sistema, ayudando a resolver la tensión actual en la cosmología.
• Metodología Futura: El enfoque desarrollado aquí (combinación de lentes + cinemática + modelado SIDM directo) establece un nuevo estándar para futuros estudios de física de materia oscura y modelado de lentes en cúmulos, especialmente con la llegada de datos más profundos de JWST y MUSE.

En conclusión, el trabajo no solo refina los límites de la auto-interacción de la materia oscura, sino que también transforma a MACS J0138 en un laboratorio de precisión para la cosmología, ofreciendo una ventana temporal clara para la observación de eventos cósmicos repetidos.