New insights on low-mass dark matter subhalo tidal tracks via numerical simulations

Resumen Técnico: Nuevas perspectivas sobre las trayectorias de marea de subhalos de materia oscura de baja masa mediante simulaciones numéricas

Planteamiento del problema
El artículo aborda la evolución de los subhalos de materia oscura (DM) de baja masa ($\lesssim 10^6 M_\odot$) mientras orbitan dentro de una galaxia anfitriona, centrándose específicamente en el entorno de la Vía Láctea (MW). Aunque estudios previos han establecido que los subhalos pierden masa debido al despojo de marea (tidal stripping), existe un debate continuo sobre la resiliencia de los subhalos con cúspides de densidad interna (específicamente "cúspides rápidas" o prompt cusps con pendientes $\gamma \geq 1.5$) y la evolución precisa de sus parámetros estructurales. Una brecha crítica en la literatura es la falta de un análisis detallado sobre cómo evolucionan estos parámetros en los pericentros orbitales (donde las fuerzas de marea son más fuertes) en comparación con los apocentros, particularmente cuando el halo anfitrión incluye componentes bariónicos realistas (disco y bulbo) y un potencial gravitatorio que evoluciona en el tiempo. Además, la relación entre la pérdida de masa y la "trayectoria de marea" resultante de los parámetros estructurales ($V_{max}$, su radio $r_{max}$ y la concentración) requiere una resolución más alta y una exploración del espacio de parámetros más amplia de lo que estaba disponible anteriormente.

Metodología
Los autores emplean el código DASH (un integrador de $N$-cuerpos optimizado para el despojo de marea), utilizando una versión mejorada introducida por Aguirre-Santaella et al. (2023). Las simulaciones rastrean un único subhalo de DM orbitando dentro de un potencial definido analíticamente y en evolución que representa la Vía Láctea, el cual incluye un halo de DM, un disco estelar y un bulbo.

• Propiedades del subhalo: Los subhalos son completamente oscuros (sin bariones) con una masa inicial de $10^6 M_\odot$. Se inicializan con un perfil de densidad NFW generalizado, variando la pendiente interna ($\gamma$) para incluir perfiles NFW estándar ($\gamma=1$) y cúspides rápidas ($\gamma=1.5$).
• Espacio de parámetros: El estudio varía la concentración inicial del subhalo ($c$), el redshift de acreción ($z_{acc}$), la energía orbital ($x_c$), la circularidad ($\eta$) y el ángulo de inclinación ($\theta$).
• Resolución: Las simulaciones utilizan hasta $N=2^{25}$ partículas para asegurar resultados robustos incluso cuando los subhalos pierden una cantidad significativa de masa (hasta fracciones de masa ligada $f_b < 10^{-3.5}$).
• Análisis: Los autores extraen $V_{max}$ y $r_{max}$ en cada instantánea (snapshot), distinguiendo entre apocentros y pericentros orbitales. Derivan una "concentración de velocidad" ($c_V$) basada en estas cantidades, la cual evita la dependencia de un radio de virial que puede dejar de existir tras el despojo.

Contribuciones clave

1. Análisis de pericentros: A diferencia de trabajos previos que analizaron principalmente las trayectorias de marea en los apocentros orbitales (donde los subhalos están más cerca del equilibrio), este estudio deriva y compara explícitamente las trayectorias de marea en los pericentros, donde el subhalo experimenta el máximo estrés de marea.
2. Evolución de cúspides rápidas: El artículo investiga la evolución de marea de subhalos con cúspides rápidas ($\gamma=1.5$), una población predicha como altamente resiliente pero que a menudo carece de resolución suficiente en las simulaciones cosmológicas estándar.
3. Potencial bariónico en evolución temporal: Las simulaciones incorporan un potencial anfitrión que evoluciona con el tiempo e incluye componentes bariónicos, ofreciendo una representación más realista del entorno de la Vía Láctea en comparación con los potenciales estáticos o de solo DM utilizados en la literatura previa.
4. Trayectorias de marea de alta resolución: Al lograr una resolución de partículas sin precedentes, los autores mapean la evolución estructural de los subhalos hasta alcanzar fracciones de masa ligada extremadamente bajas, proporcionando nuevas funciones de ajuste para las trayectorias de marea.

Resultos

• Evolución estructural: Tanto $V_{max}$ como $r_{max}$ disminuyen con el tiempo, pero $r_{max}$ se contrae más significativamente que $V_{max}$. Este encogimiento diferencial conduce a un aumento continuo en la concentración del subhalo ($c_V$).
• Aumento de la concentración: La concentración de velocidad aumenta aproximadamente dos órdenes de magnitud desde la acreción hasta el presente día. Esto es aproximadamente un orden de magnitud superior al aumento observado para los halos de campo. El aumento es más pronunciado para las trayectorias de marea de pericentro en comparación con las de apocentro.
• Trayectorias de marea ($V_{max}$ vs. $f_b$ y $r_{max}$ vs. $f_b$):
  • Para perfiles NFW estándar, los autores confirman la existencia de una trayectoria de marea donde la evolución estructural depende primordialmente de la fracción de masa ligada ($f_b$).
  • Pericentros vs. Apocentros: Las trayectorias de pericentro muestran una curvatura más fuerte para pérdidas de masa modestas y alcanzan los mismos valores que las trayectorias de apocentro solo tras un agotamiento significativo ($\log_{10} f_b < -3$).
  • Cúspides rápidas: Los subhalos con cúspides rápidas son más resilientes; en $\log_{10} f_b = -2.5$, $V_{max}$ cae a $\sim 40\%$ de su valor inicial para las cúspides rápidas, comparado con un $\sim 20\%$ para los perfiles NFW.
• Motores de la dispersión: La dispersión en las trayectorias de marea es impulsada por la relación pericentro-apocentro, la distancia precisa del pericentro, el redshift de acreción y la circularidad orbital. El redshift de acreción se identifica como el principal motor de la dispersión en la evolución de la concentración de velocidad.
• Límites de resolución: Los autores señalan que para $\log_{10}(r_{max}/r_{max,i}) \lesssim -1.5$, los resultados pueden verse afectados por la relajación de dos cuerpos artificial, limitando la fiabilidad de la trayectoria de marea en los regímenes de despojo más extremos.

Significancia
El artículo sostiene que estos hallazgos son cruciales para interpretar correctamente los datos de las poblaciones de satélites de galaxias, las búsquedas de subhalos mediante lentes gravitacionales o corrientes estelares, y las búsquedas indirectas de materia oscura (por ejemplo, rayos gamma de aniquilación).

• Detección indirecta: El hallazgo de que las concentraciones de los subhalos aumentan significativamente (en $\sim 2$ órdenes de magnitud) implica que las señales de aniquilación de estos subhalos podrían ser sustancialmente más altas de lo estimado previamente, ya que la fuerza de la señal escala con la densidad al cuadrado.
• Tasas de supervivencia: La identificación de las cúspides rápidas como altamente resilientes sugiere que una población mayor de subhalos oscuros de baja masa podría sobrevivir en las regiones internas de la Vía Láctea de lo que se pensaba anteriormente, particularmente aquellos en órbitas excéntricas que pasan cerca del centro galáctico.
• Precisión del modelado: El estudio destaca que ignorar los componentes bariónicos y la evolución temporal del potencial anfitrión, así como descuidar la dinámica de los pericentros, conduce a una comprensión incompleta de la evolución estructural de los subhalos. Las trayectorias de marea derivadas proporcionan funciones de ajuste actualizadas para modelos semianalíticos y futuras simulaciones.