Impact of subhalo dynamical friction heating on the formation of the first structures in the universe

Título: Impacto del calentamiento por fricción dinámica de subhalos en la formación de las primeras estructuras del universo

1. Planteamiento del Problema

La historia térmica de los bariones en halos de materia oscura regula fundamentalmente la formación estelar. Si bien se conocen mecanismos de retroalimentación como la fotoionización, las supernovas y la retroalimentación de AGN, el papel del calentamiento gravitacional mediante la fricción dinámica (DF) de subhalos orbitantes ha recibido menos atención, especialmente en el universo temprano.

• Contexto: En el modelo $\Lambda$CDM, las estructuras crecen jerárquicamente. Los subhalos orbitan dentro de sus halos anfitriones, perturbando el gas circundante y creando una estela de fricción dinámica que disipa energía cinética orbital en energía térmica.
• La Incógnita: Se desconocía la magnitud de este calentamiento en halos de baja masa (minihalos) a altos corrimientos al rojo ($z \sim 15$) y su capacidad para competir con el enfriamiento molecular ($H_2$), lo cual es crucial para determinar si se forman estrellas de Población III (Pop III) o si el gas colapsa directamente para formar semillas de agujeros negros (DCBH).

2. Metodología

Los autores desarrollaron un modelo analítico para cuantificar las tasas de calentamiento por fricción dinámica, validado y calibrado con datos de simulaciones cosmológicas.

• Fuente de Datos: Utilizaron la simulación hidrodinámica TNG50 (proyecto IllustrisTNG) para derivar la Función de Masa de Subhalos (SHMF) evolucionada a través de un amplio rango de corrimientos al rojo ($z=15$ a $z=0$).
• Modelo Analítico:
  • Se basaron en la teoría de perturbación lineal dependiente del tiempo de Ostriker (1999) para calcular la fuerza de arrastre y la tasa de calentamiento.
  • Integraron la SHMF sobre la masa de los subhalos para obtener la tasa de calentamiento total en un halo anfitrión.
  • Consideraron correcciones por no linealidad (cuando la perturbación gravitatoria es fuerte) y dependencias del número de Mach, utilizando datos de TNG50 para estimar factores de corrección promedio.
• Análisis de Dispersión: Reconociendo que el número de subhalos varía significativamente entre halos de la misma masa, realizaron un muestreo de Monte Carlo (500 realizaciones) asumiendo fluctuaciones de Poisson para cuantificar la dispersión en las tasas de calentamiento.
• Comparación Térmica:
  • Bajo $z$: Compararon el calentamiento por DF con las tasas de enfriamiento en halos masivos ($10^{10} - 10^{15} M_\odot$) utilizando la biblioteca de enfriamiento Grackle y perfiles de densidad NFW.
  • Alto $z$ ($z \sim 15$): Compararon el calentamiento por DF con el enfriamiento molecular ($H_2$) en minihalos ($10^6 - 10^8 M_\odot$), evaluando cómo afecta a la fracción crítica de $H_2$ necesaria para el colapso.

3. Contribuciones Clave

• Derivación de la SHMF Evolucionada: Proporcionan una parametrización analítica de la SHMF evolucionada desde $z=15$ hasta $z=0$, mostrando cómo la pendiente de la función cambia con el tiempo cósmico debido a la pérdida de masa por marea.
• Cuantificación de la Dispersión: Demuestran que la tasa de calentamiento por DF no es un valor único para una masa de halo dada, sino que varía en dos órdenes de magnitud debido a la dispersión en el número de subhalos (fluctuaciones de Poisson).
• Umbral de Resolución: Establecen que para capturar el calentamiento por DF en simulaciones, es necesario resolver subhalos con una relación de masa subhalo/halo anfitrión de $\psi \gtrsim 0.05$.
• Modelo de Calentamiento para Semillas de AGN: Presentan un mecanismo físico continuo y natural (orbitas de subhalos) que puede suprimir la formación estelar en minihalos, complementando o actuando junto con el fondo de radiación Lyman-Werner (LW).

4. Resultados Principales

A. Universo Temprano ($z \sim 15$) y Formación de Pop III:

• En minihalos libres de metales, las tasas de calentamiento por DF pueden ser comparables a las tasas de enfriamiento de $H_2$.
• El calentamiento por DF aumenta la fracción crítica de $H_2$ ($f_{crit}$) necesaria para que ocurra el colapso gravitacional y la formación de estrellas Pop III hasta en un orden de magnitud.
• Esto hace que la formación de estrellas Pop III sea más difícil, ya que se requiere una mayor abundancia de $H_2$ para superar el calentamiento.
• Incluso con una fracción media de $H_2$ de $\sim 10^{-5}$ (típica en la fase de progenitor de minihalos), el calentamiento por DF puede suprimir la formación estelar si la dispersión en el número de subhalos es alta.

B. Formación de Semillas de Agujeros Negros (DCBH):

• El calentamiento por DF favorece la formación de semillas de agujeros negros de colapso directo (DCBH) en halos de enfriamiento atómico.
• Al suprimir el enfriamiento molecular, el gas puede evitar la fragmentación en estrellas Pop III y colapsar isotérmicamente a través del enfriamiento atómico, formando estrellas supermasivas o quasi-estrellas que colapsan en agujeros negros.
• Este mecanismo es efectivo incluso cuando el fondo de radiación LW es débil ($J_{LW} < J_{crit}$), siempre que la dispersión en el número de subhalos sea suficiente para mantener el calentamiento por encima del enfriamiento.

C. Universo Local y Halos Masivos ($z \approx 0$):

• En halos masivos ($M > 10^{13} M_\odot$), el calentamiento por DF es un mecanismo importante para el apagado (quenching) de galaxias, compitiendo con la retroalimentación de AGN.
• La tasa de calentamiento media es aproximadamente 1-2 órdenes de magnitud menor que el calentamiento gravitacional total estimado en estudios previos (que ignoraban la pérdida de masa por marea), pero sigue siendo significativa.
• La gran dispersión en el número de subhalos implica que el calentamiento por DF puede ser insuficiente para apagar galaxias en halos con pocos subhalos, incluso si la masa total es alta.

5. Significado e Implicaciones

• Revisión de Modelos de Semillas de AGN: El estudio sugiere que el calentamiento por fricción dinámica es un factor crítico que debe incluirse en los modelos semi-analíticos de formación de agujeros negros. Ofrece una vía alternativa a los escenarios que requieren ráfagas extremas de radiación LW o tasas de acreción inusualmente altas.
• Requisitos de Simulación: Los resultados imponen restricciones estrictas a la resolución necesaria en simulaciones cosmológicas futuras. Si no se resuelven subhalos con $\psi \gtrsim 0.05$, se subestimará significativamente el calentamiento del gas y se malinterpretarán las condiciones para la formación de las primeras estrellas y agujeros negros.
• Física del Universo Temprano: Proporciona un mecanismo físico robusto y continuo (no dependiente de eventos transitorios raros) que explica la supresión de la formación estelar en ciertos minihalos, facilitando la formación de los agujeros negros supermasivos observados en el universo temprano ($z \approx 6-7$).

En resumen, el trabajo demuestra que la fricción dinámica de subhalos es un mecanismo de calentamiento gravitacional subestimado que juega un papel dual: suprime la formación de las primeras estrellas (Pop III) en minihalos y promueve la formación de semillas de agujeros negros masivos, con una variabilidad significativa dictada por la estadística de subhalos en cada halo anfitrión.