Detectable subhalo impacts in Milky Way streams

Resumen Técnico: Impactos de Subhalos Detectables en Corrientes de la Vía Láctea

Planteamiento del Problema
Los subhalos de materia oscura (DM) por debajo de la escala de $10^8 M_\odot$ están pobremente restringidos por las observaciones actuales, aunque su abundancia y estructura interna ofrecen pruebas críticas para las teorías de DM que van más allá del paradigma estándar de la Materia Oscura Fría (CDM), incluyendo la DM Cálida (WDM), la DM Difusa (FDM) y la DM Autointeractuante (SIDM). Si bien las corrientes estelares son sondas sensibles para subhalos de baja masa ($10^5 - 10^6 M_\odot$), estudios previos se han centrado principalmente en unas pocas corrientes bien conocidas (p. ej., GD-1, Palomar 5) o han dependido de modelos analíticos simplificados. Ha faltado una evaluación sistemática del número esperado de impactos fuertes individualmente detectables a través del catálogo creciente de corrientes de la Vía Láctea, teniendo en cuenta las incertidumbres de modelado y las capacidades observacionales de los próximos sondeos (LSST y Via). Además, definir un impacto "detectable" es complicado debido a que los encuentros de subhalos distantes producen perturbaciones de gran ángulo que son degeneradas con las incertidumbres en el modelado de la corriente suave (p. ej., historial de pérdida de masa y el potencial galáctico).

Metodología
Los autores desarrollaron un flujo de trabajo estadístico para pronosticar el número de impactos de subhalos detectables en 14 corrientes prometedoras de la Vía Láctea, seleccionadas de un catálogo de 131 candidatos (Bonaca & Price-Whelan 2025).

1. Selección y Modelado de Corrientes:

   • Las corrientes fueron seleccionadas basándose en su longitud ($>20^\circ$), masa del progenitor conocida y la capacidad de ser bien modeladas por un algoritmo de pulverización de partículas (particle spray) (Fardal et al. 2015).
   • Las edades de las corrientes ($T_{\text{form}}$) se determinaron ajustando las longitudes de las corrientes simuladas a los datos observados, variando la masa del progenitor y las tasas de pérdida de masa ($\gamma$).
   • Los criterios de selección requirieron que el ancho de la corriente simulada coincidiera con los anchos observados dentro de un factor de 2, filtrando corrientes con morfologías complejas o de orígenes inciertos.

2. Modelo Proxy Polinómico para Corrientes Suaves:

   • Para abordar la divergencia de los impactos detectables con el parámetro de impacto ($b$) y las incertididades en el modelado de la corriente suave, los autores introdujeron un modelo proxy polinómico.
   • Los observables de la corriente suave (densidad $\rho$, ángulo transversal $\phi_2$ y velocidad radial $v_r$) se ajustan con un polinomio de tercer orden.
   • Un impacto detectable se define como una desviación estadísticamente significativa de este ajuste polinómico óptimo, en lugar de de una simulación suave específica. Este enfoque absorbe las variaciones a gran escala causadas por encuentros distantes o incertidumbres de modelado, asegurando que el número de impactos detectables converja.

3. Flujo de Trabajo de Simulación de Impactos:

   • Muestreo: Los encuentros de subhalos se muestrearon durante la vida de la corriente mediante un enfoque analítico (Erkal et al. 2016) con una función de masa de subhalos CDM y un perfil radial.
   • Filtrado de Tres Capas: Para gestionar el coste computacional manteniendo la precisión, se aplicó un proceso de selección de tres niveles:
     1. Corte Analítico: Una estimación conservadora del tamaño de la señal basada en órbitas circulares.
     2. Integración de Órbita (OI): Una simulación de coste moderado que tiene en cuenta las órbitas elípticas y el potencial real de la Vía Láctea.
     3. Simulación de Alta Estadística: Simulaciones completas de pulverización de partículas para los impactos que pasaron los cortes previos, con el fin de generar perfiles de señal precisos.
   • Criterio de Detectabilidad: Un impacto se considera detectable si la relación señal-ruido al cuadrado ($SNR^2$), calculada a partir de los residuos del ajuste polinómico a través de $\rho$, $\phi_2$ y $v_r$, supera un umbral de 15.51 (correspondiente a un nivel de confianza del 95% para 8 grados de libertad).

4. Escenarios Observacionales:

   • Los pronósticos se generaron para dos escenarios: datos actuales (Gaia + DESI) y datos de un futuro cercano (LSST 10 años + sondeo espectroscópico Via).

Resultados Clave

• Tasas de Impacto Esperadas: Bajo el escenario de futuro cercano LSST+Via, solo 5 de las 14 corrientes estudiadas tienen un número esperado de impactos detectables ($N_{\text{det}}$) mayor a 0.2:
  • Jet: $5.15^{+1.10}_{-0.95}$ (el objetivo más prometedor).
  • Orphan-Chenab: $1.40^{+0.62}_{-0.47}$.
  • ATLAS-Aliqa Uma: $1.25^{+0.60}_{-0.44}$.
  • GD-1: $0.55^{+0.43}_{-0.28}$.
  • Palomar 5: $0.40^{+0.39}_{-0.23}$.
  • Con los datos actuales de Gaia+DESI, solo tres corrientes muestran detecciones esperadas no nulas, con tasas significativamente menores (0.1–0.45).
• Convergencia mediante Sustracción Polinómica: Sin la sustracción del proxy polinómico, el número de impactos detectables crece con el parámetro de impacto $b$ y la masa del subhalo $M$, fallando en la convergencia. El método de sustracción logra regular este crecimiento, concentrando los impactos detectables dentro de $b < 5r_s$ (donde $r_s$ es el radio de escala del subhalo).
• Contribución de los Observables: El $SNR^2$ total se comparte entre los observables. En promedio, la densidad ($\rho$) contribuye con ~49%, el ángulo transversal ($\phi_2$) con ~36%, y la velocidad radial ($v_r$) con ~16%. Incluir $\phi_2$ y $v_r$ es crucial para recuperar las tasas de detección que de otro modo se subestimarían en análisis basados únicamente en la densidad.
• Dependencia de los Modelos de DM:
  • WDM y FDM: En los límites de masa actuales ($m_{\text{WDM}} = 5.9$ keV, $m_{\text{FDM}} = 1.4 \times 10^{-20}$ eV), el número de impactos se ve suprimido por un factor de $\sim 4$ en comparación con CDM.
  • SIDM: Un escenario extremo de SIDM con halos de colapso de núcleo ($M < 10^8 M_\odot$) resulta en un aumento de orden $O(1)$ en el número de impactos debido a densidades centrales más altas.
• Incertidumbres Sistemáticas: Las predicciones dependen fuertemente de la población de subhalos asumida. Las variaciones en la distribución radial de los subhalos (p. ej., basadas en mediciones de satélites de la Vía Láctea frente a simulaciones FIRE-2) pueden alterar las tasas de impacto predichas por un factor de pocos puntos.

Significación y Reivindicaciones
El artículo afirma proporcionar el primer pronóstico sistemático de impactos de subhalos individualmente detectables a través de una amplia muestra de corrientes de la Vía Láctea, yendo más allá de los estudios de caso de una sola corriente. Al introducir un modelo proxy polinómico, los autores resuelven el problema teórico de los conteos de impactos divergentes a grandes parámetros de impacto, ofreciendo una definición robusta y basada en datos de la detectabilidad.

Los autores enfatizan que sus resultados representan "estimaciones del mejor caso" bajo supuestos idealizados (potencial de la Vía Láctea estático, función de masa CDM analítica). Declaran explícitamente que distinguir entre modelos de DM (CDM vs. WDM/FDM/SIDM) basándose únicamente en el número de impactos será desafiante debido a las fluctuaciones de Poisson y las incertidumbres sistemáticas en la población de subhalos. En su lugar, el trabajo identifica corrientes específicas (notablemente Jet) como objetivos de alta prioridad para futuros sondeos para restringir las propiedades de la DM y potencialmente reconstruir las propiedades de los perturbadores individuales (masa, órbita, tiempo de impacto). El estudio sirve como una guía para la asignación de recursos observacionales y destaca la necesidad de mejorar el modelado de las corrientes para aprovechar plenamente las corrientes estelares en la física de la materia oscura.