Rates of tidal disruption events from constrained cosmological simulations of the local Universe: population properties and implications for transient surveys

Resumen Técnico: Tasas de Eventos de Disrupción de Marea a partir de Simulaciones Cosmológicas Restringidas del Universo Local

Planteamiento del Problema
Los eventos de disrupción de marea (TDE, por sus siglas en inglés) sirven como sondas críticas para estudiar las poblaciones de agujeros negros supermasivos (SMBH) quiescentes. Sin embargo, existe una discrepancia significativa entre las predicciones teóricas y las mediciones observacionales. Los modelos idealizados de relajación de dos cuerpos predicen tasas de $\gtrsim 10^{-4} \text{ yr}^{-1} \text{ gal}^{-1}$, mientras que las restricciones observacionales de los sondeos ópticos, de rayos X y de infrarrojo medio suelen arrojar tasas un orden de magnitud menores ($\sim 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$). Además, los enfoques analíticos tradicionales suelen depender de extrapolaciones en 2D y suposiciones simplificadas sobre los perfiles de densidad galáctica (por ejemplo, esferas isotérmicas singulares), fallando al no contabilizar la compleja interacción de la formación de estructuras cosmológicas, las historias de fusión y la diversidad morfológica específica (con cúspides o núcleos) de las galaxias anfitrionas. A medida que se aproxima el Legacy Survey of Space and Time (LSST) del Observatorio Vera C. Rubin, se requieren estimaciones de las tasas de TDE espacialmente resueltas y refinadas para interpretar con precisión los próximos datos de transitorios.

Metodología
Los autores utilizan la "Simulación del Entorno Local" (SLOW), una simulación cosmológica restringida de última generación del Universo local, centrándose específicamente en seis regiones de "zoom-in" de alta resolución: los supercúmulos de Coma, Hércules, Shapley, Virgo y Perseo, y el cúmulo de galaxias de Fornax. Estas regiones abarcan extensiones radiales de 55–92 Mpc, cubriendo entornos que van desde escalas de grupos hasta supercúmulos masivos ($M_{500} \approx 3.5 \times 10^{13} - 1.7 \times 10^{15} M_\odot$).

La metodología procede a través de varios pasos rigurosos de filtrado y cálculo:

1. Extracción de Agujeros Negros: Se extraen los SMBH de los volúmenes de simulación, registrando sus masas ($M_\bullet$) y parámetros de espín adimensionales ($\tilde{a}$).
2. Filtrado de la Secuencia Principal: Para asegurar la validez de los cálculos de la dinámica del cono de pérdida (loss-cone), un esquema de filtrado estricto aísla los agujeros negros de la "secuencia principal" dinámicamente estables. Esto excluye sistemas no ligados, dominados por rotación o erosionados por marea mediante cuatro criterios:
   • Corte de Desplazamiento (Offset Cut): Los BH deben estar dentro de 1.5 kpc del centro de su subhalo.
   • Corte de Velocidad: La velocidad relativa debe ser inferior a la dispersión de velocidad estelar local.
   • Corte de Anisotropía Fraccional (FA Cut): El entorno estelar local debe exhibir un $FA < 0.4$, asegurando una esfera aproximadamente isotrópica y soportada por presión, requerida para los formalismos estándar del cono de pérdida.
   • Corte de Relación de Materia Oscura: Se requiere una relación de masa de materia oscura a masa estelar $> 1$ para excluir núcleos desnudos y erosionados.
3. Clasificación Morfológica: Las galaxias anfitrionas se clasifican como "con cúspide" ($g \ge 2$) o "con núcleo" ($g \le 1$) basándose en la pendiente de densidad estelar en 3D medida dentro de 1 kpc. Este radio se elige para evitar artefactos de suavizado numérico y servir como un proxy de la estructura nuclear no resuelta.
4. Corrección de Eficiencia Relativista: El estudio aplica un factor de eficiencia relativista ($\eta$) derivado de Kesden (2012). Este factor cuantifica la fracción del cono de pérdida que resulta en un TDE observable en lugar de una captura directa, contabilizando el espín del agujero negro y el límite de masa de Hills.
5. Cálculo de la Tasa: Las tasas volumétricas de TDE ($\Gamma_{vol}$) se calculan utilizando el marco de Stone & Metzger (2016), que aplica tasas base distintas para entornos de cúspide y de núcleo, ponderadas por la eficiencia relativista $\eta$.

Contribuciones Clave y Resultados

• Tasas Volumétricas y por Agujero Negro: El estudio deriva una tasa volumétrica promedio de TDE de $\approx 600 \text{ Gpc}^{-3} \text{ yr}^{-1}$ a través de los seis entornos. La tasa promedio de TDE por agujero negro se encuentra estable en $\approx 4.5 \times 10^{-5} \text{ yr}^{-1}$, un valor consistente con las restricciones observacionales actuales pero significativamente menor que las predicciones teóricas idealizadas.
• Distribución Espacial vs. Masa Total: El artículo demuestra que los rendimientos de TDE están impulsados por la demografía de los agujeros negros y la concentración espacial, más que por la masa total del cúmulo. Mientras que el presupuesto absoluto de TDE está abrumadoramente dominado por las galaxias satélites con cúspide en halos de cúmulos extendidos, la tasa volumétrica ($\Gamma_{vol}$) está fuertemente influenciada por la densidad numérica de los agujeros negros ($n_\bullet$).
• Depleción del Núcleo y Eficiencia: Las tasas de los núcleos centrales se reducen drástamente debido a la depleción dinámica (fusiones que destruyen subhalos de baja masa) y a las restricciones de captura directa (límite de Hills). En consecuencia, el presupuesto total de TDE no está dominado por las galaxias centrales masivas, sino por la población de satélites de menor masa y con cúspide en las periferias de los cúmulos.
• Dependencia Ambiental:
  • Hércules: A pesar de ser un supercúmulo masivo y en ensamblaje activo con una alta densidad numérica de agujeros negros, exhibe la eficiencia por agujero negro ($\tilde{\Gamma}$) más baja debido al crecimiento de la masa de los agujeros negros por fusiones, lo que empuja a muchos BH hacia el límite de masa de Hills donde $\eta \to \text{0}$.
  • Fornax: Este entorno de baja masa es el más eficiente por agujero negro ($\tilde{\Gamma} \approx 1.07$) porque alberga una fracción desproporcionada de agujeros negros de baja masa y no fusionados que evitan el corte de Hills, manteniendo altas eficiencias del cono de pérdida. Crucialmente, Fornax exhibe la densidad numérica de agujeros negros ($n_\bullet$) más baja y el radio de número medio más grande ($R_{1/2,\bullet} = 31.26$ Mpc) entre los seis entornos, indicando la distribución más dispersa espacialmente de los agujeros negros.
• Discrepancia con Extrapolaciones en 2D: El estudio destaca que las extrapolaciones analíticas tradicionales en 2D sobreestiman las densidades de superficie centrales. Las densidades de tasa de superficie en 2D simuladas son aproximadamente un orden de magnitud menores que los valores de la literatura derivados de proyecciones en 2D, ya que el verdadero presupuesto de TDE se distribuye a través del halo extendido en lugar de concentrarse en el núcleo.

Significado
El artículo afirma cerrar la brecha entre las predicciones teóricas y las mediciones observacionales mediante el uso de un marco cosmológico totalmente en 3D que resuelve los entornos estelares locales. Al ir más allá de las simplificadas extrapolaciones en 2D y los perfiles de densidad idealizados, los autores proporcionan un "presupuesto de TDE" más realista que tiene en cuenta la compleja historia de ensamblaje del Universo local.

Los resultados sugieren que la tasa observada de $\sim 3 \times 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$ está sólidamente respaldada por las simulaciones cosmológicas cuando se contabiliza la depleción del núcleo dinámico y los límites de captura relativista. La tasa derivada de $4.5 \times 10^{-5} \text{ yr}^{-1}$ se presenta como un límite superior teórico robusto para los próximos sondeos de transitorios, reconociendo que la física de sub-malla de la simulación puede sobreproducir ligeramente las galaxias enanas de baja masa que dominan el rendimiento absoluto. El estudio enfatiza que las tasas de TDE no son constantes universales, sino que están fundamentalmente dictadas por las propiedades macroscópicas de la red cósmica, específicamente la concentración espacial de los agujeros negros y la historia de fusión de sus entornos anfitriones.