TDEs on FIRE: Illuminating the Cosmic Evolution of Tidal Disruption Rates

Resumen Técnico: TDEs en FIRE

Planteamiento del Problema
Aunque los eventos de disrupción de marea (TDE, por sus siglas en inglés) han sido caracterizados extensamente en el universo local, su prevalencia y evolución cósmica en altos desplazamientos al rojo ($z > 1$) permanecen mayormente sin constreñir, tanto observacional como teóricamente. Las muestras observacionales actuales están limitadas por el flujo y son abrumadoramente locales, con solo un puñado de candidatos con chorro (jetted) detectados en $z \sim 1$. Las predicciones teóricas enfrentan incertidumbres significativas debido a la falta de perfiles estelares resueltos en galaxias de alto desplazamiento al rojo y a las pobres restricciones de la función de masa de agujeros negros (BHMF) para agujeros negros de baja masa, que dominan la población de TDE. En consecuencia, no existe una predicción robusta para la evolución en el desplazamiento al rojo de las tasas de TDE (TDR) ni de su dependencia con las propiedades de la galaxia anfitriona a través del tiempo cósmico.

Metodología
Este estudio utiliza las simulaciones de zoom cosmológico FIRE-2 (Feedback In Realistic Environments), específicamente la suite Massive Halo (simulaciones A1, A2, A4 y A8), para computar las tasas intrínsecas de TDE en el rango de desplazamiento al rojo $z = 1–10$. La metodología integra tres componentes clave:

1. Marco de Simulación: Las simulaciones FIRE-2 emplean el código GIZMO con alta resolución espacial (alcanzando escalas de parsecs) y un modelado exhaustivo de la formación estelar, evolución y retroalimentación estelar. La suite Massive Halo incluye un modelo sofisticado para la formación y el crecimiento de agujeros negros masivos (MBH), sembrándolos con $M_{seed} = 1.4 \times 10^4 M_\odot$ cuando los halos de materia oscura exceden un umbral de masa específico. Los agujeros negros son tratados como partículas sumidero (sink particles) sin colisiones con masa dinámica artificial para asegurar que permanezcan en el centro de la galaxia.
2. Análisis del Perfil Estelar: Los autores extraen perfiles de densidad estelar radial dentro de 1 kpc de los centros galácticos (en lugar del centro del BH, para considerar posibles desplazamientos por fricción dinámica). Estos perfiles se ajustan con una función de ley de potencia doble. Los sistemas con pendientes internas no físicas (más pronunciadas que las pendientes externas) o ruido numérico son filtrados. La muestra se restringe a galaxias con $M_{gal} > 10^7 M_\odot$ y $M_{BH} < 10^8 M_\odot$ (el umbral por encima del cual las estrellas son tragadas enteras en lugar de ser disrumpidas por marea).
3. Cálculo de la Tasa: Las tasas de TDE se calculan utilizando la teoría del cono de pérdida (loss cone theory) y el código de acceso público phaseflow. El código resuelve la función de distribución estelar $f(E)$ mediante la inversión de Eddington a partir de los perfiles de densidad ajustados. La tasa de consumo total se deriva integrando la dispersión de estrellas hacia el cono de pérdida debido a la relajación de dos cuerpos, asumiendo una población estelar monocromática similar al Sol ($m_\star = 1 M_\odot, r_\star = 1 R_\odot$).

Contribuciones Clave y Resultados

• Evolución del TDR con el Desplazamiento al Rojo: El estudio proporciona la primera predicción de la evolución del TDR a través de $z = 1–10$. El TDR promedio aumenta desde el universo temprano, alcanza un máximo de aproximadamente $4 \times 10^{-4} \text{ yr}^{-1}$ cerca de $z \sim 2.5$, y declina a $\sim 10^{-5} \text{ yr}^{-1}$ en $z = 1$. Se proporciona un ajuste de ley de potencia para $z \ge 2.5$: $\log_{10}(\text{TDR}) \approx z^{-0.2} - 4$.
• Correlación con la Formación Estelar: Se encuentra una fuerte correlación entre el TDR y la tasa de formación estelar (SFR) de la galaxia anfitriona. El TDR sigue la historia de la SFR cósmica, alcanzando su máximo en $z \sim 2.5$. La relación de escala se deriva como $\text{TDR} \propto \text{SFR}^{0.69}$.
• Dependencia de las Masas ($M_{BH}$ y $M_{gal}$):
  • Masa del Agujero Negro: El TDR aumenta con $M_{BH}$ en el régimen de agujeros negros de masa intermedia (IMBH), alcanza un máximo en $M_{BH} \sim 1.8 \times 10^5 M_\odot$, y luego declina para los agujeros negros supermasivos (SMBH). Esta tendencia de "giro" (turnover) coincide cualitativamente con los hallazgos en el universo local, a pesar de las diferencias en la normalización y las pendientes.
  • Masa de la Galaxia: De manera similar, el TDR aumenta con $M_{gal}$ para galaxias pequeñas, alcanza un máximo en $M_{gal} \sim 5 \times 10^9 M_\odot$, y declina para sistemas más masivos.
• Impulsores Estructurales Estelares: El TDR está fuertemente impulsado por las propiedades estelares centrales. Densidades estelares internas ($\rho_0$) más altas y pendientes de densidad interna ($\gamma$) más pronunciadas se correlacionan con tasas más altas. El estudio destaca que los TDR más altos ocurren en sistemas con $\rho_0$ alta, $\gamma$ pronunciada y radios de influencia ($r_{inf}$) de los agujeros negros pequeños.
• Galaxias Satélite: Las galaxias satélite (anfitrionas no primarias) exhiben TDR comparables a las galaxias primarias. Su contribución fraccional al TDR total aumenta significativamente en altos desplazamientos al rojo, lo que sugiere que son cruciales para sondear los IMBH y el ensamblaje temprano de galaxias. La mayoría de los TDE de satélites en la muestra caen dentro de los límites de resolución angular de instalaciones de próxima generación como el Telescopio Espacial Nancy Grace Roman, el Observatorio Vera Rubin y ULTRASAT.

Significancia y Reivindicaciones
El artículo sostiene que las simulaciones cosmológicas, específicamente FIRE-2, ofrecen una vía prometedora para constreñir la evolución cósmica de las tasas de TDE. El trabajo demuestra que las tendencias cualitativas del TDR con respecto a la masa del agujero negro y de la galaxia, observadas en el universo local, persisten desde el alto desplazamiento al rojo hasta el día de hoy, lo que sugiere relaciones de escala consistentes entre el agujero negro y la galaxia a través del tiempo cósmico.

Los autores reconocen modestamente las limitaciones, incluyendo la ausencia de retroalimentación de AGN en la suite Massive Halo (lo que puede aumentar artificialmente las SFRs y las densidades centrales en $z \lesssim 2$), la incapacidad de resolver completamente los cúmulos estelares nucleares debido a la resolución finita, y la restricción de la simulación a $z=1$ (lo que impide establecer restricciones en $z=0$). A pesar de estas salvedades, el estudio establece un marco para comparar las predicciones teóricas con las observaciones venideras de los sondeos de área amplia y alta cadencia, con el objetivo de validar modelos de crecimiento de agujeros negros, dinámica estelar y coevolución de galaxias.