Detection of a dark matter subhalo in the strongly lensed system PJ011646

Resumen Técnico: Detección de un Subhalo de Materia Oscura en el Sistema Fuertemente Lenteado PJ011646

Planteamiento del Problema
La naturaleza de la materia oscura sigue siendo un problema fundamental sin resolver en cosmología. Mientras que el modelo de Materia Oscura Fría (CDM) predice una formación jerárquica de halos autoenlazados que se extienden hasta masas muy bajas, los modelos de Materia Oscura Cálida (WDM) predicen una supresión de halos de baja masa por debajo de una escala de corte característica. Distinguir entre estos modelos requiere mediciones estadísticas de la población de subhalos, particularmente en el régimen de masas donde las predicciones divergen ($M_{200} \sim 10^7 - 10^9 M_\odot$). Un desafío principal en la detección de estas subestructuras mediante lentes gravitacionales fuertes es la "degeneración" entre las señales genuinas de subhalos y la complejidad angular no modelada en la distribución de masa suave de la galaxia lente. Los modelos simplificados de ley de potencia elíptica (PL) a menudo fallan en capturar desviaciones en galaxias de tipo temprano (ETG), lo que conduce a detecciones espurias de subhalos o al enmascaramiento de los reales. Además, las detecciones anteriores han estado limitadas por la resolución angular de los telescopios ópticos (por ejemplo, HST) y la contaminación por la luz de la galaxia lente.

Metodología
Los autores presentan un análisis de lente gravitacional fuerte del sistema PJ011646 (lente $z=0.555$, fuente $z=2.125$) utilizando observaciones de alta resolución ($\sim 0.1$ segundos de arco) de la continuidad del polvo obtenidas con el Atacama Large Millimeter/submillimeter Array (ALMA) en la Banda 7. El análisis utiliza el conjunto de software de código abierto PyAutoLens, PyAutoFit y PyAutoGalaxy.

La metodología procede en tres etapas principales:

1. Modelado Macro: La distribución de masa de la lente se modela utilizando un perfil de ley de potencia elíptica (EPL) con cizallamiento externo. Para abordar las limitaciones de los modelos elípticos simples, los autores incorporan multipolos de Fourier de orden superior (órdenes $m=3$ y $m=4$, y adicionalmente $m=1$) para capturar desviaciones angulares en la densidad de masa. La convergencia de estos multipolos se parametriza para evitar degeneraciones asociadas con fases periódicas.
2. Reconstrucción de la Fuente: La brillo superficial de la fuente se reconstruye utilizando un método de inversión bayesiana pixelada (Suyu et al. 2006) con pixelación de malla de Voronoi e interpolación de vecinos naturales. Esto se realiza directamente en el plano $uv$ para aprovechar los datos interferométricos completos.
3. Búsqueda de Subhalos: Se realiza una búsqueda basada en una cuadrícula en el plano de la imagen para detectar perturbaciones. Un perfil esférico Navarro-Frenk-White (NFW) se añade como perturbador en cada celda de la cuadrícula, y se calcula la diferencia de log-verosimilitud bayesiana ($\Delta \ln Z$) con respecto al modelo macro de mejor ajuste. Se adopta un umbral de $\Delta \ln Z > 10$ para identificar candidatos robustos. Los autores también prueban un perfil NFW generalizado (gNFW) con una pendiente interna libre ($\gamma$) para verificar desviaciones del perfil CDM estándar.

Resultados Clave

• Mejora del Modelo Macro: La inclusión de multipolos de tercer y cuarto orden ($m=3, 4$) mejora significativamente el ajuste a los datos, produciendo un aumento en la evidencia bayesiana de $\Delta \ln Z = 52.1$ en relación con un modelo PL puro. Las amplitudes de estos multipolos son $\simeq 1.5\%$ de la convergencia, consistente con observaciones de ETG cercanas. Un modelo que incluye un multipolo de primer orden ($m=1$) proporciona una mejora adicional marginal ($\Delta \ln Z = 6.5$), sugiriendo un posible compañero físico pero sin impulsar la señal principal.
• Detección de Subhalos: Utilizando el modelo macro PL+MP de mejor ajuste, la búsqueda en cuadrícula identifica dos regiones en el plano de la imagen con $\Delta \ln Z > 10$. Sin embargo, un análisis de seguimiento revela que estas dos regiones mapean a la misma ubicación en el plano de la fuente, indicando una degeneración. Un modelo de dos subhalos no mejora la evidencia sobre un modelo de un solo subhalo, y la posterior para el segundo candidato produce concentraciones no físicas ($c_{200} > 1000$). En consecuencia, los autores concluyen que solo un subhalo es físicamente real, ubicado a la derecha del arco norte.
• Propiedades del Subhalo: El subhalo detectado tiene una masa $M_{200} = 2.78^{+0.43}_{-0.66} \times 10^{10} M_\odot$ y una concentración $c_{200} = 30^{+5}_{-7}$. Esta concentración es consistente con las expectativas para un subhalo NFW despojado por mareas dentro de un halo padre. La significancia de la detección es de $\sim 5.8\sigma$ en relación con el modelo macro PL+MP.
• Restricciones Robustas de Masa: Debido a la degeneración entre el radio de escala NFW y la normalización, la masa total $M_{200}$ no está fuertemente restringida. En su lugar, los autores definen un "radio robusto" ($R_{sub} \approx 2$ kpc) donde la incertidumbre fraccional sobre la masa encerrada proyectada se minimiza. Dentro de este radio, la masa encerrada es $M_{sub} = 3.57^{+0.16}_{-0.14} \times 10^9 M_\odot$.
• Límites de Sensibilidad: En las celdas de la cuadrícula donde no se detecta ningún subhalo ($\Delta \ln Z < 10$), los autores derivan límites superiores sobre las masas de subhalos detectables. En las regiones más sensibles (a lo largo de los arcos lenteados brillantes), la masa mínima detectable con una significancia de $3\sigma$ es $M_{200} \approx 8 \times 10^8 M_\odot$, asumiendo la relación masa-concentración de Ludlow et al. (2016).

Significado y Afirmaciones
El artículo afirma que este trabajo demuestra la capacidad de la imagenación de continuidad de ALMA con resolución subarcosegundo para sondear la subestructura de materia oscura en un régimen de masas ($M_{200} \sim 10^9 M_\odot$) donde los modelos CDM y WDM predicen abundancias diferentes. Al modelar explícitamente la complejidad angular de la galaxia lente mediante multipolos, los autores mitigan las incertidumbres sistemáticas que históricamente han plagado las búsquedas de subhalos. La detección de un subhalo con propiedades consistentes con las predicciones del CDM, junto con la derivación de límites superiores estrictos en regiones sin detección, proporciona una prueba robusta de la formación de estructura a pequeña escala. Los autores señalan que, aunque la resolución angular es comparable a la del HST, la alta magnificación del sistema y la ausencia de contaminación por luz de primer plano en los datos de ALMA permiten sondear masas más bajas de las alcanzables previamente con datos ópticos. El estudio se posiciona como un paso hacia búsquedas sistemáticas y estadísticamente significativas de subestructura de materia oscura utilizando la muestra creciente de galaxias formadoras de estrellas polvorientas (DSFG) fuertemente lenteadas observadas con ALMA.