Accretion history dependence of the halo depletion radius

Aquí presento un resumen técnico detallado del artículo "Dependencia de la historia de acreción en el radio de agotamiento del halo" (Accretion history dependence of the halo depletion radius) de Zhou y Han, basado en la versión de borrador proporcionada.

1. Planteamiento del Problema

La definición física de los límites externos de los halos de materia oscura ha sido un desafío persistente en la cosmología. La definición tradicional basada en la sobredensidad esférica (radio virial, $R_{vir}$) es arbitraria y sufre de "pseudo-evolución". Aunque el radio de salpicadura ($R_{sp}$) se propuso como un límite dinámico asociado al primer apocentro de la materia recién acrecida, existe otro límite dinámico complementario: el radio de agotamiento ($R_{id}$).

El radio de agotamiento se define como el radio donde la evolución temporal de la densidad local cambia de signo (de crecimiento a disminución), marcando la transición entre la región de crecimiento interno y la región de agotamiento externo. Aunque se ha observado que $R_{id} \approx 2 R_{vir}$, su dependencia cuantitativa con la historia de acreción, la masa del halo, el corrimiento al rojo ($z$) y otras propiedades secundarias no había sido caracterizada sistemáticamente. El objetivo de este trabajo es investigar cómo la historia de acreción moldea $R_{id}$ y cómo se compara con $R_{sp}$.

2. Metodología

Los autores utilizaron una simulación cosmológica de cuerpo N a gran escala de la suite CosmicGrowth:

• Simulación: $3072^3$ partículas de materia oscura en una caja de $600 \, h^{-1}\text{Mpc}$, con una masa de partícula de $5.54 \times 10^8 \, h^{-1}M_\odot$. Cosmología $\Lambda$CDM plana ($\Omega_m = 0.268$).
• Catálogo de Halos: Se identificaron halos usando el algoritmo Friends-of-Friends (FoF) y se construyeron árboles de fusión con el código HBT+.
• Selección de Muestra: Se aplicaron criterios estrictos para asegurar resolución y pureza:
  • Masa virial mínima: $\log_{10}[M_{vir}/(h^{-1}M_\odot)] \ge 12$.
  • Historia evolutiva definida: El subhalo central debe permanecer identificado como tal durante al menos un tiempo dinámico ($t_{dyn}$).
  • Aislamiento: Se excluyeron halos dentro de $2.5 R_{vir}$ de un halo más masivo.
• Definición de $R_{id}$: Se derivó directamente de la ecuación de continuidad. El radio de agotamiento corresponde al mínimo local en el perfil de tasa de flujo de masa (MFR, Mass Flow Rate), definido como $MFR(r) = 4\pi r^2 \rho(r) v_r(r)$.
• Tasa de Acreción de Masa (MAR): Se utilizó la tasa de crecimiento logarítmico de la masa virial sobre un tiempo dinámico ($\Gamma_{dyn}$).
• Análisis: Se apilaron perfiles de densidad, velocidad radial y MFR para halos con propiedades similares, normalizados por cantidades viriales ($R_{vir}, V_{vir}, M_{vir}$). Se compararon halos con y sin fusiones mayores recientes.

3. Contribuciones Clave y Resultados

A. Dependencia de la Tasa de Acreción ($\Gamma_{dyn}$)

• Existe una fuerte correlación entre $R_{id}/R_{vir}$ y $\Gamma_{dyn}$.
• Tendencia: A medida que aumenta la tasa de acreción reciente ($\Gamma_{dyn}$), el radio de agotamiento normalizado disminuye. Los halos que acrecen rápidamente tienen una capa de crecimiento más compacta.
• Relación Funcional: La dependencia se ajusta bien a una forma exponencial:
  $$\frac{R_{id}}{R_{vir}} = A + B e^{-\Gamma_{dyn}/C}$$
  Donde para $z=0$, la relación es $\frac{R_{id}}{R_{vir}} = 0.79 + 1.79 e^{-\Gamma_{dyn}/3}$.

B. Dependencia de la Masa y el Corrimiento al Rojo ($z$)

• Masa: La dependencia de la masa del halo es débil. A un $\Gamma_{dyn}$ fijo, la posición de $R_{id}$ es casi invariante con la masa, aunque los halos más masivos muestran perfiles de flujo ligeramente más profundos.
• Corrimiento al Rojo: A diferencia del radio de salpicadura ($R_{sp}$), que muestra una evolución débil con $z$ cuando se escala con $R_{vir}$, el radio de agotamiento ($R_{id}$) exhibe una evolución sistemática y clara con el corrimiento al rojo.
  • A mayor $z$, la relación $R_{id}/R_{vir}$ vs. $\Gamma_{dyn}$ es más pronunciada (pendiente más suave y valores más altos para $\Gamma_{dyn} \gtrsim 1$).
  • Esta evolución es impulsada principalmente por eventos dinámicos (fusiones).

C. El Papel de las Fusiones Mayores

• Se identificó que las fusiones mayores recientes (acreción de un progenitor con relación de masa > 1:9 en el último tiempo dinámico) son el motor principal de la evolución con el corrimiento al rojo.
• Al excluir halos con fusiones mayores:
  • La dependencia de $z$ en la relación $R_{id}/R_{vir} - \Gamma_{dyn}$ se suprime significativamente.
  • La relación se vuelve universal y se ajusta a un modelo de colapso esférico auto-similar donde $\Omega_m \approx 1$.
  • En este régimen "suave", $R_{id}$ y $R_{sp}$ difieren por una fracción constante de $R_{vir}$: $R_{id} \approx 0.54 R_{vir} + R_{sp}$.

D. Dependencia de Propiedades Secundarias (Dicotomía)

El estudio revela una dicotomía en los impulsores de $R_{id}$:

1. Acreción Lenta ($\Gamma_{dyn} \lesssim 1$): La tasa de acreción es el factor dominante. La dependencia de otras propiedades es mínima.
2. Acreción Rápida ($\Gamma_{dyn} \gtrsim 1$): La tasa de acreción no es suficiente. Propiedades secundarias como la forma, el espín, la concentración y, crucialmente, el tiempo de formación del subhalo central ($a_{1/2}$) juegan un papel significativo.
   • Halos con formación temprana ($a_{1/2}$ pequeño) pero acreción total rápida muestran un $R_{id}$ más pequeño. Esto indica un crecimiento dominado por fusiones de satélites masivos en lugar de acreción suave, lo que altera el perfil de flujo de masa (supresión del flujo saliente).

E. Explicación de la Universalidad $R_{id} \approx 2 R_{vir}$

El trabajo explica por qué estudios anteriores reportaron un ratio casi constante de $\approx 2$:

• No es una constante física fundamental.
• Es un efecto estadístico resultante de dos tendencias compensatorias:
  1. La tasa de acreción mediana ($\Gamma_{dyn}$) aumenta con el corrimiento al rojo (lo que tendería a reducir $R_{id}/R_{vir}$).
  2. La relación $R_{id}/R_{vir} - \Gamma_{dyn}$ se desplaza hacia arriba con el corrimiento al rojo (lo que tendería a aumentar $R_{id}/R_{vir}$).
• La compensación de estas dos tendencias en una muestra de halos trazada desde $z=0$ produce un ratio aparentemente constante a lo largo del tiempo cósmico.

4. Significado e Implicaciones

• Nuevo Diagnóstico Físico: El radio de agotamiento se establece como una sonda física sensible no solo a la tasa de acreción reciente, sino también al modo detallado de acreción (suave vs. fusiones, anisotrópico).
• Complementariedad: Mientras que $R_{sp}$ y $R_{id}$ están fuertemente correlacionados, $R_{id}$ es más sensible a la historia de fusiones y a la evolución cosmológica de fondo, ofreciendo información complementaria a $R_{sp}$.
• Desafíos Teóricos: Los resultados desafían los modelos de colapso esférico puro, especialmente en regímenes de acreción rápida donde la asimetría y las perturbaciones de las fusiones son dominantes.
• Aplicabilidad Observacional: Dado que $R_{id}$ ha sido detectado mediante lentes gravitacionales débiles y cinemática de satélites, estos resultados permiten utilizarlo para inferir la historia de acreción y el estado dinámico de los halos en observaciones reales, más allá de la simple estimación de la tasa de acreción.

En resumen, el artículo demuestra que el radio de agotamiento es un indicador dinámico rico que codifica la historia de crecimiento de los halos, siendo especialmente sensible a la naturaleza de las fusiones mayores en el universo local.