Perturber-Driven Dynamics of Supermassive Black Hole Binaries in Galaxy Merger

Resumen Técnico: Dinámica de Binarias de Agujeros Negros Supermasivos Impulsada por Perturbadores

Planteamiento del Problema
La excentricidad orbital de las binarias de agujeros negros masivos (MBHB) en el momento de su formación es un parámetro crítico que determina el fondo estocástico de ondas gravitacionales (GWB) detectable por las matrices de temporización de púlsares (PTA). Simulaciones previas de alta resolución de $N$-cuerpos han demostrado que la excentricidad de formación exhibe una dispersión significativa entre ejecuciones, particularmente para órbitas de fusión casi radiales. Si bien se ha demostrado que esta dispersión disminuye con la resolución como $\sim 1/\sqrt{N}$, lo cual es consistente con el ruido de muestreo de Poisson, no está claro si la subestructura física —específicamente perturbadores masivos como cúmulos globulares, núcleos despojados o agujeros negros de masa intermedia (IMBH)— introduce una estocasticidad astrofísica genuina que supere el umbral de ruido numérico. La cuestión central es si el efecto acumulativo de muchos perturbadores de baja masa (régimen difusivo) o de encuentros masivos poco frecuentes (régimen impulsivo) puede aleatorizar significativamente la excentricidad de formación de las MBHB en galaxias elípticas masivas.

Metodología
Los autores realizaron re-simulaciones de alta resolución de una fusión de galaxias mayor (Merger 6 de la simulación cosmológica IllustrisTNG100-1) utilizando el código de $N$-cuerpos $Griffin$. La configuración de la simulación involucró:

• Condiciones Iniciales: Una fusión casi radial ($e_{orb} \approx 0.987$) entre dos galaxias modeladas como sistemas colisionales (bulbos estelares, halos de materia oscura y MBH centrales) siguiendo perfiles de Hernquist. La galaxia primaria retuvo una pequeña fracción de gas ($f_{gas} \approx 0.19$), pero la simulación trató el sistema como colisional para centrarse en sistemas estelares de larga duración.
• Escenarios de Perturbadores: Se compararon tres suites de simulaciones, cada una con cuatro realizaciones independientes:
  1. Control: Sin perturbadores masivos; el bulbo estelar suave sirve como fondo.
  2. Suite de $10^7 M_\odot$: $3,162$ perturbadores de masa igual ($M_p = 10^7 M_\odot$) redistribuyendo el $10\%$ de la masa del bulbo primario. Esto corresponde a una relación de masa perturbador-binaria $\mu_p \approx 3.2 \times 10^{-3}$.
  3. Suite de $10^8 M_\odot$: $316$ perturbadores de masa igual ($M_p = 10^8 M_\odot$) redistribuyendo la misma fracción de masa total. Esto corresponde a una $\mu_p \approx 3.2 \times 10^{-2}$.
• Criterio de Formación de la Binaria: El tiempo de formación de la binaria ($t_b$) se definió como el primer mínimo local de la excentricidad tras la unión gravitacional permanente. La excentricidad de formación ($e_b$) y su desviación estándar ($\sigma_e$) a través de las realizaciones fueron las métricas principales.
• Diagnósticos Dinámicos: El estudio analizó los residuales en la energía orbital ($E$), el momento angular ($J_z$), el torque ($\tau_z$) y el errante del centro de masa (COM) de la binaria respecto a las ejecuciones de control emparejadas, alineados por $t_b$.

Contribuciones Clave y Resultados
El artículo establece una transición entre los regímenes de dispersión difusiva y casi impulsiva basada en la relación de masa del perturbador respecto a la binaria ($\mu_p$):

1. Umbral de Ruido Numérico: La suite de control sin perturbadores produjo una dispersión de excentricidad de formación de $\sigma_e \approx 0.11$, consistente con el umbral de ruido de Poisson esperado para la resolución adoptada ($N_{bulge} \approx 2 \times 10^6$).
2. Régimen Difusivo ($10^7 M_\odot$): La suite de $10^7 M_\odot$ produjo una $\sigma_e \approx 0.115$, la cual es estadísticamente indistinguible del control. Los residuales dinámicos en energía y momento angular fueron suaves, y no se observaron picos de torque significativos. Esto confirma que para $\mu_p \approx 3.2 \times 10^{-3}$, el sistema permanece en el régimen difusivo donde los encuentros individuales son demasiado débiles para dejar una impronta discernible en la excentricidad de formación.
3. Régimen Casi Impulsivo ($10^8 M_\odot$): La suite de $10^8 M_\odot$ exhibió una dispersión significativamente mayor, $\sigma_e \approx 0.26$ (un factor de $\sim 2.4$ por encima del control). Esta dispersión excesiva fue acompañada de:
   • Saltos discretos en el semieje mayor y la energía orbital.
   • Grandes residuales alineados con eventos en el momento angular y picos de torque cerca de $t_b$.
   • Un errante incrementado del centro de masa de la binaria ($\sim 1.25$–$1.4\times$ mayor que el control).
     Estas firmas indican que los encuentros cercanos individuales en este régimen entregan patadas impulsivas lo suficientemente grandes como para alterar el momento angular y la excentricidad de la binaria antes de que comience la fase de endurecimiento (hardening).
4. Transición Teórica: Los resultados se alinean con la teoría de dispersión binaria-simple, donde la transición está gobernada por $\mu_p$. El caso de $10^8 M_\odot$ se aproxima al límite casi impulsivo ($\Delta e_{single} \sim \mu_p$), mientras que el caso de $10^7 M_\odot$ permanece difusivo.

Significancia y Reivindicaciones
El artículo sostiene que, si bien los perturbadores masivos pueden inducir estocasticidad física en la excentricidad de formación de las MBHB, este efecto depende altamente de la escala de masa de los perturbadores en relación con la binaria.

• Contexto Astrofísico: En las galaxias elípticas masivas (los anfitriones primarios de las MBHB relevantes para las PTA), la población dominante de perturbadores consiste en cúmulos globulares ($10^5$–$10^6 M_\odot$) y núcleos despojados. Para las masas de las binarias relevantes para las PTA ($10^8$–$10^{10} M_\odot$), estos objetos corresponden a $\mu_p \lesssim 10^{-4}$, situándolos firmemente en el régimen difusivo.
• Conclusión: Por consiguiente, es poco probable que la aleatorización de la excentricidad impulsada por perturbadores sea una vía significativa para la mayoría de las fusiones de MBHB relevantes para la detección del GWB. La estocasticidad observada en las simulaciones actuales está dominada por el ruido numérico en lugar de la subestructura física para poblaciones de perturbadores realistas.
• Advertencias: Los autores señalan que su conclusión depende del margen de la relación de masa más que de la significancia estadística de la única mejora de $10^8 M_\odot$ (que es marginal con solo cuatro realizaciones). También reconocen que las fusiones ricas en gas (no modeladas completamente aquí) podrían introducir dinámicas diferentes a través de nubes moleculares gigantes o discos circumbinarios, y que la configuración orbital específica (casi radial) puede influir en el umbral exacto. Sin embargo, la conclusión astrofísica dominante sigue siendo que las poblaciones de perturbadores realistas en galaxias elípticas masivas se encuentran por debajo del umbral casi impulsivo requerido para aleatorizar significativamente la excentricidad de formación.