McFACTS IV: Electromagnetic Counterparts to AGN Disk Embedded Binary Black Hole Mergers

Résumé technique : McFACTS IV : Contreforts électromagnétiques des fusions de binaires de trous noirs enchâssées dans les disques d'AGN

Énoncé du problème
Les disques d'accrétion des Noyaux Actifs Galactiques (AGN) constituent un canal théorique majeur pour la formation et la fusion de binaires de trous noirs (BBH), une population de plus en plus détectée par le réseau LIGO-Virgo-KAGRA (LVK). Bien que les fusions de BBH dans les disques d'AGN soient uniques parmi les canaux de formation par leur potentiel à produire des contreforts électromagnétiques (EM) via des interactions avec le gaz environnant, l'absence de modèles statistiques prédictifs entrave actuellement l'identification de ces contreforts dans les données archivistiques et les relevés temporels en temps réel (par exemple, ZTF, LSST). De plus, la fraction d'événements LVK provenant du canal AGN ($f_{BBH,AGN}$) demeure incertaine, et distinguer les véritables contreforts EM de la variabilité intrinsèque des AGN est difficile en raison d'une localisation par ondes gravitationnelles (GW) médiocre et de la nature brillante et variable des AGN.

Méthodologie
Ce travail utilise une version mise à jour du code Monte Carlo For AGN Channel Testing and Simulation (McFACTS v.0.4.0), un outil de synthèse de population rapide en 1D pour la formation dynamique de binaires compactes. L'étude introduit une nouvelle capacité à générer de manière auto-cohérente des distributions de luminosités EM bolométriques (jets et chocs) correspondant aux événements GW prédits et détectables par LVK.

La méthodologie comprend :

1. Synthèse de population : Simulation de fusions de BBH au sein d'un modèle de disque d'AGN de Sirko-Goodman entourant un trou noir supermassif (SMBH) de $10^8 M_\odot$. Le code suit les fusions hiérarchiques, l'évolution du spin et la migration.
2. Modélisation de l'émission EM :
   • Luminosité de choc ($L_{shock}$) : Calculée sur la base de l'énergie dissipée lorsque le recul par recul (kick) de l'onde gravitationnelle déplace l'halos de Hill du rémanent dans le disque environnant (McKernan et al. 2019).
   • Luminosité de jet ($L_{jet}$) : Modélisée en supposant que le rémanent entre dans un état de disque magnétiquement arrêté (MAD), lançant un jet via le mécanisme de Blandford-Znajek. La puissance est paramétrée par le taux d'accrétion de Bondi-Hoyle-Lyttleton, le spin du trou noir ($\chi$) et l'efficacité du jet.
3. Critères d'observabilité : Un composant critique est le calcul du temps de rupture du jet ($t_{bo}$). Un jet est considéré comme potentiellement observable uniquement s'il brise le disque optiquement épais avant que le moteur ne s'éteigne ($t_{bo} < T_{jet}$). Si $t_{bo} > T_{jet}$, l'émission est « étouffée » et émerge sur une échelle de temps de diffusion de photons, la rendant indiscernable de la variabilité de l'AGN.
4. Exploration de l'espace des paramètres : L'étude fait varier plusieurs paramètres physiques pour évaluer leur impact sur les taux et les luminosités des contreforts EM, notamment :
   • La durée de vie du disque d'AGN ($\tau_{AGN}$).
   • L'étendue radiale du disque ($R_{outer}$).
   • La force de la turbulence du disque (migration stochastique).
   • La fonction de masse initiale (IMF) des trous noirs (variation des indices de la loi de puissance et des masses de pointe).
   • La largeur de la distribution de spin initiale.

Résultats clés

• Dominance des jets sur les chocs : Dans le modèle de disque de Sirko-Goodman fiduciel, les luminosités de choc sont généralement subdominantes par rapport à la luminosité de l'AGN hôte ($L_{AGN} \sim 10^{45}$ erg s$^{-1}$) et sont peu susceptibles d'être distinguables. À l'inverse, les jets alimentés par l'accrétion peuvent surpasser significativement l'AGN.
• Rôle des fusions hiérarchiques : Les régions denses du disque, spécifiquement les pièges de migration et les marécages (ex: à $R \sim 700 R_g$), favorisent efficacement les fusions hiérarchiques. Ces processus produisent des rémanents de haute masse et de haut spin (générations 2 et supérieures) capables d'alimenter des jets lumineux.
• Dépendance de la masse de chirp : La probabilité de produire un contrefort EM observable est fortement corrélée à la masse de chirp ($\mathcal{M}$). Pour le modèle fiduciel, la probabilité de produire une éruption observable augmente avec la masse, atteignant $\sim 100 \%$ pour des fusions avec $\mathcal{M} \in [44, 88] M_\odot$ (avant de prendre en compte l'obscuration géométrique).
• Sensibilité des paramètres :
  • Durée de vie du disque : Des disques plus durables ($\tau_{AGN} \ge 1$ Myr) permettent une fusion hiérarchique plus étendue, augmentant la population de rémanents à haut spin et haute masse, et donc de jets brillants.
  • IMF : Les IMF riches en masses élevées (ex: ségrégation de masse ou distributions de type pic autour de $34 M_\odot$) améliorent considérablement l'efficacité de la croissance hiérarchique, conduisant à des rémanents plus massifs et des jets plus brillants par rapport aux IMF standards en loi de puissance.
  • Turbulence : Une augmentation de la turbulence du disque lisse les pièges de migration, réduisant légèrement l'efficacité des fusions hiérarchiques et la variance de la luminosité des jets résultants.
• Contraintes géométriques : Même lorsqu'un jet brise le disque, des facteurs géométriques limitent l'observabilité. En supposant un nombre égal d'AGN de Type I et de Type II, seulement $\sim 1/4$ des événements de rupture sont géométriquement détectables (Type I face-à-face avec la rupture vers l'observateur).

Signification et revendications
L'article affirme présenter la première étude de synthèse de population des contreforts EM spécifiquement pour le canal AGN. Sa contribution primaire est le développement d'un cadre qui lie les observables GW (masse de chirp, spin) aux propriétés EM prédites (luminosité, temps de rupture) au sein d'un environnement de disque physiquement motivé.

Les auteurs affirment que :

1. Sélection des cibles : Les campagnes de suivi EM devraient prioriser les événements GW avec des masses de chirp $\mathcal{M} \gtrsim 44 M_\odot$, car ils présentent la probabilité statistique la plus élevée de produire des contreforts EM observables.
2. Échelles de temps : L'émission observable peut émerger aussi tard que $\sim 4$ mois après la fusion, ce qui suggère que les campagnes de surveillance doivent s'étendre sur des échelles de temps de plusieurs semaines à plusieurs mois.
3. Contraintes de modèle : La détection ou la non-détection de contreforts dans les futurs relevés (ex: LSST, Einstein Telescope) fournira des contraintes critiques sur les propriétés des disques d'AGN (densité, échelle de hauteur, durée de vie) et sur la nature de la population de BBH (IMF, alignement de spin).
4. Potentiel multi-messager : Des associations confirmées entre les événements GW et les éruptes EM dans les disques d'AGN offrent une voie pour contraindre indépendamment la constante de Hubble ($H_0$).

L'étude reste modeste quant aux garanties de détection spécifiques, reconnaissant que les incertitudes liées à la variabilité des AGN, à la physique du lancement des jets et aux angles de vue géométriques limitent actuellement la capacité à associer de manière définitive des événements GW spécifiques à des contreforts EM. Ce travail sert de guide pour les futures stratégies d'observation et les raffinements théoriques.