The Density Profile of Dynamical Halos

Resumen Técnico: El Perfil de Densidad de los Halos Dinámicos

Planteamiento del Problema
Las propiedades fundamentales de los halos de materia oscura, particularmente sus perfiles de densidad, se caracterizan tradicionalmente mediante ajustes analíticos (por ejemplo, NFW, Einasto) a la distribución de masa total dentro de un radio de referencia como $R_{200m}$. Sin embargo, este enfoque confunde dos poblaciones dinámicas distintas: las partículas que están gravitacionalmente ligadas y en órbita dentro del halo, y las partículas que están en proceso de caída (infalling) pero que aún no están ligadas. Propuestas recientes sugieren definir un "halo dinámico" estrictamente como la colección de partículas en órbita. Si bien esta definición evita mezclar la estructura interna con la acreción en curso, las propiedades estructurales (media y dispersión) de estos perfiles en órbita no habían sido caracterizadas completamente. Específicamente, no estaba claro qué grados de libertad estructurales permanecen a una masa en órbita fija, la magnitud de la dispersión intrínseca, y si las correlaciones con la historia de formación son propiedades genuinas del componente en órbita o artefactos de la mezcla de poblaciones.

Metodología
Los autores utilizaron una simulación de Materia Oscura Fría (CDM) ($1024^3$ partículas en una caja de $1 \, h^{-1}$ Gpc) ejecutada con GADGET-2. Generaron un catálogo de halos dinámicos utilizando el buscador de halos Rockstar como punto de partida, y luego aplicaron el algoritmo de García et al. para etiquetar las partículas como "en órbita" (orbiting) o "en caída" (infalling) basándose en su tiempo de acreción ($a_{acc}$) y su velocidad radial media en el tiempo.

• Definición de Órbita: Las partículas con tiempos de acreción tempranos y velocidades radiales medias cercanas a cero se etiquetan como en órbita; aquellas con tiempos de acreción tardíos y velocidades negativas se etiquetan como en caída.
• Agrupación por Masa: Los halos se agruparon en nueve estrechos intervalos de masa ($\log M_{orb} \in [13.40, 15.00]$).
• Modelo Paramétrico: Los perfiles de densidad en órbita se ajustaron utilizando el modelo de Salazar et al., definido como:
  $$\rho_{orb}(x|M_{orb}) = A \left(\frac{x}{\epsilon}\right)^{-\alpha(x|M_{orb})} \exp\left(-\frac{x^2}{2}\right)$$
  donde $x \equiv r/r_h$, $r_h$ es el radio del halo (escala de truncamiento exponencial), y la pendiente $\alpha(x)$ transiciona de 0 a un valor asintótico $\alpha_\infty$. El parámetro $\epsilon$ se fijó en 0.037.
• Estrategia de Ajuste: Se ajustaron halos individuales para determinar los mejores valores de $r_h$ y $\alpha_\infty$. Posteriormente, los autores investigaron la correlación entre estos parámetros y la dependencia con la historia de acreción, caracterizada por una variable de tiempo de formación relativa $a_{RF}$.

Contribuciones Clave y Resultados

1. Un Único Grado de Libertad Efectivo:
   El estudio demuestra que, a una masa en órbita fija, los dos parámetros libres del perfil ($r_h$ y $\alpha_\infty$) están fuertemente correlacionados. Los halos más extendidos (mayor $r_h$) poseen pendientes internas más pronunciadas (mayor $\alpha_\infty$). Esta correlación permite que el perfil de un halo dinámico individual sea caracterizado por un único parámetro libre: el radio del halo $r_h$. La relación viene dada por:
   $$\Delta\alpha_\infty = 0.05 + \ln R$$
   donde $R = r_h / r_{h,st}(M_{orb})$ y $\Delta\alpha_\infty$ es la desviación de la pendiente promediada (stacked).

2. Dispersión en el Radio del Halo:

   • A una masa en órbita fija, la distribución de los radios de los halos es aproximadamente log-normal. La dispersión intrínseca en $\ln r_h$ es $\sigma_{\ln r_h | M_{orb}} \approx 0.16$ (16%).
   • Los autores investigaron la dependencia de esta dispersión con la historia de acreción. Definieron un tiempo de formación relativa $a_{RF}$ basado en el percentil 60 de la distribución del tiempo de acreción de las partículas.
   • Existe una correlación entre $r_h$ y $a_{RF}$: los halos de formación tardía son más compactos (menor $r_h$).
   • Considerar esta correlación reduce la dispersión en el radio del halo a una masa y tiempo de formación fijos a $\sigma_{\ln r_h | M_{orb}, a_{RF}} \approx 0.11$ (11%). Los autores señalan que esta reducción es marginal, lo que sugiere que el tiempo de formación relativa por sí solo es insuficiente para explicar completamente la dispersión.

3. Comparación con la Tasa de Acreción:
   Los autores compararon la correlación entre el radio del halo y el tiempo de formación relativa frente a la correlación entre el radio del halo y la tasa de acreción por unidad de masa. Encontraron que las tendencias son cualitativamente similares, con la correlación de la tasa de acreción mostrando restricciones ligeramente más estrictas pero colas no gaussianas más pronunciadas. Advierten que no pudieron estimar de manera autoconsistente las tasas de acreción para los halos dinámicos en este análisis específico.

4. Desempeño del Modelo:
   Los ajustes paramétricos a los perfiles en órbita producen residuos comparables a los ajustes NFW aplicados a los perfiles de masa total (residuo mediano $\approx -0.02$, dispersión $\sim 0.15$). Esto valida el uso del modelo de Salazar et al. para el componente en órbita.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera caracterización de la dependencia de la historia de acreción de los perfiles en órbita para halos definidos dinámicamente. Al aislar el componente en órbita, los autores argumentan que las correlaciones medidas (por ejemplo, entre estructura y tiempo de formación) son propiedades genuinas del sistema ligado y no artefactos definicionales causados por la mezcla de material en caída.

La significancia principal reside en la capacidad de "leer directamente" la masa y el tiempo de formación de un halo a partir de sus propiedades estructurales ($r_h$ y $\alpha_\infty$). A diferencia de los halos de sobredensidad esférica estándar, donde el radio está definicionalmente ligado a la masa encerrada, los halos dinámicos exponen cómo su extensión y estructura surgen de la dinámica orbital. Esta calibración sirve como una línea base necesaria para el modelado directo (forward modeling), permitiendo descomposiciones dinámicas más precisas de los perfiles de densidad de los halos y mejorando la interpretación de las distribuciones de masa inferidas de datos de lensing y agrupamiento (clustering). Los autores enfatizan que, si bien su modelo no reemplaza los ajustes empíricos para las distribuciones de materia total, es esencial para comprender la estructura dinámica subyacente.