Measurement of Substructure from the Kinematics of the GD-1 Stellar Stream

Resumen Técnico: Medición de la Subestructura a partir de la Cinemática de la Corriente Estelar GD-1

Planteamiento del Problema
En el paradigma de la Materia Oscura Fría (CDM), se predice que las galaxias albergan una vasta población de subhalos de materia oscura, que van desde escalas galácticas hasta escalas de la masa de la Tierra. Sin embargo, los subhalos por debajo del umbral de formación de galaxias ($\approx 10^8 M_\odot$) se espera que carezcan de materia bariónica, lo que los hace invisibles para la observación directa. En consecuencia, probar el paradigma CDM en escalas pequeñas requiere trazadores gravitacionales indirectos. Las corrientes estelares, que son estructuras cinemáticamente frías formadas por la disrupción de marea de cúmulos globulares, son sensibles a las perturbaciones de subhalos de baja masa. Mientras que estudios previos han utilizado las brechas y espuelas morfológicas en corrientes como GD-1 para inferir interacciones con subhalos, este trabajo aborda la necesidad de restringir la dispersión de velocidad intrínseca de la corriente. Una dispersión de velocidad medida significativamente mayor de la predicha por potenciales suaves de la Vía Láctea indicaría un calentamiento impulsado por la subestructura de materia oscura.

Metodología
Los autores analizan la corriente estelar GD-1 utilizando una combinación de datos observacionales y modelado de teoría de perturbaciones:

1. Recolección de Datos y Medición de la Dispersión de Velocidad:

   • El estudio utiliza 160 estrellas miembros únicas identificadas a partir del catálogo de GD-1 basado en datos (Starkman et al. 2025) con una probabilidad de membresía $>75\%$.
   • Las velocidades radiales ($v_r$) se compilan de cuatro sondeos espectroscópicos: DESI (entrega de datos temprana), SDSS DR9, LAMOST DR8 y MMT Hectochelle.
   • Para controlar la contaminación por estrellas binarias, los autores utilizan observaciones repetidas de 34 estrellas en diferentes instrumentos.
   • Se emplea un modelo bayesiano jerárquico para separar los errores observacionales, el "jitter" de velocidad inducido por binarias ($\sigma_{binary}$) y la dispersión de velocidad intrínseca de la corriente ($\sigma_{\Delta v_r}$). La corriente se divide en tres segmentos a lo largo de su longitud ($\phi_1 \in [-60, -40], [-40, -20], [-20, 0]$ grados) para medir la dispersión local.

2. Modelado Dinámico:

   • Los autores emplean el marco streamsculptor, que utiliza la teoría de perturbaciones para modelar el efecto acumulativo de cientos de impactos de subhalos en la corriente.
   • Los modelos de corrientes no perturbadas se generan mediante un enfoque de "particle-spray" para masas de progenitor de $2 \times 10^4 M_\odot$ y $10^5 M_\odot$, y edades dinámicas de 3.5, 5 y 8 Gyr, dentro de un potencial estático de la Vía Láctea (MilkyWayPotential2022).
   • Los subhalos se modelan como perfiles de Hernquist. El estudio explora dos escenarios para romper las degeneraciones entre la densidad numérica de subhalos y su estructura interna:
     • Modelo I: Asume un escenario de Materia Oscura Cálida (WDM) con una supresión en la función de masa de subhalos (SHMF) por debajo de una masa de modo medio ($M_{hm}$) y una relación masa-tamaño constante ($r_s/r_{s,cdm}$).
     • Modelo II: Asume una SHMF estándar de CDM ($M_{hm}=0$) pero permite que la relación masa-tamaño varíe como una ley de potencia rota con la masa, probando la concentración dependiente de la masa.

3. Inferencia:

   • Los autores utilizan Computación Bayesiana Aproximada (ABC) para comparar las simulaciones de la corriente contra la dispersión de velocidad observada.
   • Las simulaciones se generan para varios parámetros de la población de subhalos (normalización, $M_{hm}$, concentración). Las muestras se aceptan si la distancia euclidiana entre la dispersión de velocidad simulada y la observada (y las restricciones de ancho de la corriente) cae dentro de una tolerancia $\epsilon$.

Contribuciones Clave y Resultados

• Dispersión de Velocidad Intrínseca: El estudio mide la dispersión de velocidad radial intrínseca de GD-1 como aproximadamente $2\text{--}5$ km s$^{-1}$ a lo largo de su longitud. La región con la mayor dispersión ($\phi_1 \in [-40, -20]$ grados) muestra un valor de $4.8^{+1.2}_{-0.7}$ km s$^{-1}$. Esto representa una desviación de $4\sigma$ respecto a los modelos de corrientes no perturbadas formadas en un potencial suave, que predicen dispersiones de solo $\approx 0.7\text{--}1.1$ km s$^{-1}$.
• Restricciones de la Población de Subhalos: La dispersión observada puede explicarse por una población de subhalos de materia oscura. La fracción de masa de la Vía Láctea inferida en subhalos es $f_{sub} = 0.05^{+0.08}_{-0.03}$ (68% de confianza), consistente con las expectativas de CDM ($\approx 7.5\%$).
• Compacidad de los Subhalos: Un hallazgo principal es la preferencia por subhalos que son más compactos que las predicciones estándar de CDM.
  • En el Modelo I, si la fracción de masa de subhalos es del $10\%$, el radio de escala se restringe a $r_s/r_{s,cdm} = 0.3^{+0.3}_{-0.2}$ (68% de confianza).
  • En el Modelo II, hay una preferencia por subhalos cada vez más compactos a bajas masas ($M < 10^7 M_\odot$). Para subhalos con $M_{200} = 10^6 M_\odot$, la concentración $c_{-2}$ se restringe a ser $>190$ con un 68% de confianza, significativamente mayor que la predicción de CDM de $\approx 43$.
• Escenarios Alternativos: Los datos pueden explicarse ya sea por numerosos impactos con subhalos de baja masa o por un único impacto con un subhalo muy compacto ($M \gtrsim 10^8 M_\odot$). Los autores señalan que un único impacto con un subhalo compacto es suficiente para generar la dispersión observada, lo cual es consistente con las interpretaciones previas de la característica de la "espuela" de GD-1.

Significancia y Reivindicaciones
El artículo afirma proporcionar una nueva restricción independiente sobre la subestructura de materia oscura utilizando la cinemática de las corrientes estelares, complementando las restricciones previas derivadas de la morfología de la densidad de la corriente. Los resultados sugieren una posible desviación de la relación masa-tamaño estándar de CDM a bajas masas de subhalos, donde se espera que los subhalos estén dominados puramente por materia oscura. Específicamente, la preferencia por concentraciones más altas (o radios de escala más pequeños) que las predicciones de CDM se alinea con teorías de materia oscura alternativas, como la Materia Oscura Autointeractuante (SIDM) o la Materia Oscura Atómica, que predicen densidades centrales aumentadas.

Los autores enfatizan que sus restricciones son empíricas y agnósticas a modelos específicos de física de partículas, sirviendo como un punto de referencia para probar diversas teorías de materia oscura. Reconocen que existen degeneraciones entre el número de impactos de subhalos y su concentración (la "degeneración tasa-concentración"), que podrían resolverse en trabajos futuros mediante el modelado conjunto de la densidad y la cinemática de la corriente. El estudio concluye que, si bien la fracción de masa de los subhalos es consistente con CDM, la estructura interna de estos subhalos puede requerir física más allá del paradigma estándar de CDM, particularmente a bajas masas.