JWST lensed quasar dark matter survey IV: Stringent warm dark matter constraints from the joint reconstruction of extended lensed arcs and quasar flux ratios

Resumen Técnico: JWST Lensing Survey IV – Restricciones de Materia Oscura Cálida mediante Reconstrucción Conjunta

Planteamiento del Problema
La naturaleza de la materia oscura (DM) sigue siendo desconocida, y el paradigma de la Materia Oscura Fría (CDM) predice una jerarquía sin escala de halos hasta masas muy pequeñas. Teorías alternativas, como la Materia Oscura Cálida (WDM), predicen una supresión de la estructura por debajo de una escala de streaming libre específica, caracterizada por la masa de modo medio ($m_{hm}$). El lente gravitacional fuerte de cuásares con imágenes cuádruples ofrece una sonda única de esta subestructura al medir perturbaciones en las posiciones de las imágenes y las razones de flujo. Sin embargo, los análisis previos que dependían únicamente de las posiciones de fuentes puntuales y las razones de flujo se han visto limitados por las degeneraciones entre el perfil de masa del deflector principal (macromodelo) y la población de subestructura. Las incertidumbres en el macromodelo (por ejemplo, elipsidad, cizalla externa y estructura angular) pueden imitar o enmascarar las firmas de los subhalos de materia oscura, lo que conduce a restricciones más débiles sobre las propiedades de la DM. Además, el microlente estelar ha contaminado históricamente las mediciones de la razón de flujo para fuentes compactas, aunque esto se ha mitigado mediante las recientes observaciones de JWST de regiones de polvo cálido extendidas.

Metodología
Este artículo presenta un marco de modelado hacia adelante que reconstruye simultáneamente arcos lanteados extendidos, posiciones de imágenes y razones de flujo para una muestra de 28 cuásares con imágenes cuádruples. El análisis utiliza un marco de inferencia bayesiana para restringir los hiperparámetros que describen la población de subestructura de DM.

• Integración de Datos: El estudio combina las razones de flujo de JWST MIRI (provenientes de la región de polvo cálido del cuásar, evitando el microlente estelar) con datos astrométricos e imágenes de arcos lanteados extendidos de JWST (NIRCam y MIRI) y HST.
• Modelado Conjunto: A diferencia de trabajos anteriores que trataban los arcos lanteados y las razones de flujo por separado, este análisis emplea pesos de muestreo de importancia ($w_{img}$) para propagar las restricciones de los arcos extendidos hacia los parámetros del macromodelo. Esto rompe la covarianza entre el macromodelo y la población de subestructura, estrechando las restricciones en las razones de flujo predichas por el modelo.
• Modelado Físico:
  • Subestructura: El marco genera poblaciones completas de subhalos, halos en la línea de visión y cúmulos globulares utilizando el software pyHalo. Incorpora un modelo de evolución de marea (basado en Du et al.) que predice los perfiles de densidad de los subhalos basándose en las propiedades de la caída (infall), contabilizando la pérdida de masa y los cambios de concentración.
  • Efectos de WDM: El modelo modifica la función de masa de halo de CDM y la relación concentración-masa para dar cuenta de la supresión por streaming libre en WDM, parametrizada por la masa de modo medio $m_{hm}$.
  • Macromodelo: El deflector principal se modela utilizando un perfil de ley de potencia elíptica (EPL) aumentado con perturbaciones multipolares elípticas (m=1, 3, 4) para dar cuenta de la estructura angular (por ejemplo, boxiness/diskiness) que de otro modo podría hacerse pasar por subestructura.
• Inferencia: La función de verosimilitud se evalúa generando millones de realizaciones (más de 174 millones para la muestra completa) de modelos de lente. El análisis margina sobre los parámetros de molestia (macromodelo, luz de la fuente, realizaciones de subestructura) para inferir las distribuciones posteriores de los hiperparámetros, específicamente la masa de modo medio $m_{hm}$ y la amplitud de la función de masa de subhalos $\Sigma_{sub}$.

Contribuciones Clave

1. Ruptura de Degeneraciones: El principal avance metodológico es el modelado simultáneo de arcos lanteados extendidos junto con las razones de flujo y las posiciones de las imágenes. Este enfoque reduce significativamente la incertidumbre en el macromodelo, rompiendo así la covarianza entre la abundancia de subhalos ($\Sigma_{sub}$) y la escala de streaming libre ($m_{hm}$).
2. Mejora en el Modelado de Subhalos: El artículo implementa un modelo de evolución de marea más riguroso físicamente que predice las propiedades de los subhalos en el momento de la caída, en lugar de en el redshift del lente. Esto tiene en cuenta la dependencia de la erosión por marea de la concentración de la caída, lo cual es particularmente relevante en WDM donde los halos son menos concentrados.
3. Fuentes de Perturbación Exhaustivas: El análisis incluye explícitamente cúmulos globulares y perturbaciones multipolares elípticas como fuentes de ruido en la razón de flujo, asegurando que las restricciones sobre WDM no estén sesgadas por estructuras de pequeña escala no modeladas.
4. Tamaño de la Muestra: Este trabajo analiza la mayor muestra de cuásares con imágenes cuádruplas hasta la fecha para la inferencia de DM (28 sistemas), aprovechando todas las capacidades del sondeo de lentes de JWST.

Resultos

• Restricciones de Materia Oscura Cálida:
  • Asumiendo un prior log-uniformo en la amplitud de la función de masa de subhalos, la inclusión de arcos lanteados mejora el límite de la masa de modo medio en 0.6 dex en comparación con el uso exclusivo de posiciones de imágenes y razones de flujo.
  • Al aplicar priors informativos sobre $\Sigma_{sub}$ derivados del modelo semianalítico galacticus, el análisis infiere $m_{hm} < 10^{7.4} M_\odot$ (factor de Bayes de 10:1).
  • Al aplicar priors basados en simulaciones N-cuerpo (Symphony suite), la restricción se ajusta a $m_{hm} < 10^{7.2} M_\odot$.
  • Estos límites corresponden a límites inferiores en la masa de una partícula de reliquia térmica de WDM de 6.5 keV y 7.4 keV, respectivamente.
  • Los límites de exclusión al 95% son $m_{hm} < 10^{6.8} M_\odot$ y $m_{hm} < 10^{6.7} M_\odot$, que corresponden a masas de partículas $> 9.6$ keV y $> 10.2$ keV.
• Abundancia de Subhalos de CDM:
  • Asumiendo el límite de CDM ($m_{hm} < 10^5 M_\odot$), el estudio mide la densidad de masa superficial proyectada de los subhalos como $1.7^{+2.6}_{-1.2} \times 10^7 M_\odot \text{kpc}^{-2}$ con una confianza del 95%.
  • Esto corresponde a una fracción de masa proyectada de $3.2^{+5.1}_{-2.3}\%$.
  • Esta medición es estadísticamente consistente con las predicciones de galacticus, pero muestra una ligera tensión ($\sim 1.8\sigma$) con simulaciones recientes de N-cuerpo, que predicen una abundancia de subhalos menor.

Significancia
El artículo afirma que estos resultados representan algunas de las restricciones más fuertes sobre WDM hasta la fecha y la medición más precisa de la abundancia de subhalos alrededor de lentes fuertes. La significancia radica en la demostración de que el modelado conjunto de arcos lanteados extendidos y razones de flujo aísla eficazmente la señal de la subestructura de materia oscura de las incertidumbres en el perfil de masa del deflector principal. Los autores señalan que sus restricciones superan los límites previos de lentes fuertes y son competitivas o más fuertes que otras sondas de pequeña escala como el bosque Lyman-$\alpha$ y el conteo de galaxias enanas. Este trabajo establece un marco robusto para futuras inferencias de DM utilizando las grandes muestras de lentes que se esperan de los próximos sondeos (Euclid, Rubin, Roman), destacando el potencial transformador de combinar la imagen de alta resolución basada en el espacio con el modelado hacia adelante avanzado.