Connecting the long-term variability behaviour of active galactic nuclei to their central engines

Résumé technique : Relier le comportement de variabilité à long terme des noyaux actifs de galaxies à leurs moteurs centraux

Énoncé du problème
Les Noyaux Actifs de Galaxies (NAG) présentent une variabilité multi-longueurs d'onde pilotée par des jets relativistes émis par des trous noirs supermassifs. Bien que l'analyse de la Densité Spectrale de Puissance (DSP) soit une méthode standard pour identifier les échelles de temps caractéristiques de la variabilité (spécifiquement le coude à basse fréquence dans le spectre de bruit rouge), l'interprétation physique de ces échelles de temps fait toujours l'objet de débats. Des études antérieures sur les binaires X et les NAG suggèrent que les coudes de la DSP évoluent avec la masse du trou noir ($M_{BH}$) et le taux d'accrétion ($\dot{M}$), mais la caractéristique spécifique de la courbe de lumière correspondant au coude de la DSP radio reste floue. Il est incertain si l'inverse de la fréquence de coude correspond à la durée des éruptions individuelles ou à la séparation temporelle entre elles. De plus, la relation entre ces échelles de temps de variabilité et la taille physique de la région d'émission (souvent associée à la région de transition du jet) et aux paramètres du moteur central nécessite une investigation supplémentaire utilisant un échantillon plus vaste que ceux précédemment disponibles.

Méthodologie
Les auteurs ont analysé des courbes de lumière radio à long terme de 54 NAG surveillés à l'Observatoire Radio Metsähovi de l'Université Aalto (MRO) à 37 GHz. L'échantillon comprenait 11 sources dont les échelles de temps de la DSP étaient déjà contraintes (d'après Kankkunen et al. 2025, « Article I ») et 43 sources supplémentaires (incluant 3 cas limites) sélectionnées pour leur forte densité de données et leurs éruptions marquées.

La méthodologie centrale a impliqué :

1. Décomposition de la courbe de lumière : Le code Magnetron a été employé pour décomposer les courbes de lumière en éruptions individuelles en utilisant une modélisation hiérarchique bayésienne et un échantillonnage imbriqué diffusif. Un processus d'Ornstein-Uhlenbeck (OU) a été utilisé pour modéliser la variabilité de fond sous-jacente, empêchant la mauvaise identification de variations turbulentes mineures comme des éruptions.
2. Extraction des échelles de temps : Pour chaque source, les auteurs ont calculé la durée moyenne et la séparation moyenne des pics des éruptions les plus brillantes (définies comme celles dépassant l'amplitude moyenne dans un tirage posterior). Les temps de montée ont également été extraits pour les éruptions les plus brillantes (amplitude > 50 % du maximum) afin de servir de substituts pour la taille de la région d'émission.
3. Estimation des paramètres :
   • Facteurs Doppler ($\delta_{var}$) : Estimés en utilisant les températures de brillance de variabilité dérivées des temps de montée et des amplitudes des éruptions, comparées à une distribution gaussienne de températures de brillance intrinsèques.
   • Facteurs de Lorentz ($\Gamma$) : Calculés en utilisant $\delta_{var}$ et les vitesses apparentes des jets ($\beta_{app}$) issues de l'enquête MOJAVE.
   • Paramètres du moteur central : Les masses des trous noirs ($M_{BH}$) et les luminosités des disques d'accrétion ($L_{acc}$) ont été estimées en utilisant les luminosités des raies larges (H$\beta$, MgII ou [OIII]) et les valeurs de largeur à mi-hauteur (FWHM). Les taux d'accrétion de masse normalisés ($\dot{m}$) ont été dérivés en supposant une efficacité radiative de $\eta = 0,4$.
4. Analyse statistique : Des tests de corrélation de rang de Kendall ont été effectués pour évaluer les associations entre les échelles de temps dérivées (durées d'éruptions, séparations et temps de montée) et les paramètres physiques ($M_{BH}$, $\dot{m}$, $L_{acc}$, $\Gamma$). Des analyses de corrélation partielle ont été menées pour tenir compte des dépendances potentielles au décalage vers le rouge ($z$).

Résultats clés

• Échelle de temps de la DSP vs Caractéristiques de la courbe de lumière : L'inverse de la fréquence de coude de la DSP correspond le mieux à la durée moyenne des éruptions les plus brillantes. Bien que les séparations moyennes des éruptions soient également corrélées à l'échelle de temps de la DSP, les durées d'éruptions montrent un accord un-à-un plus fort, en particulier pour les 11 sources contraintes.
• Lien avec le moteur central : En utilisant les durées et séparations d'éruptions comme substituts pour l'échelle de temps de la DSP, les auteurs ont trouvé une corrélation positive modérée avec le rapport de la masse du trou noir sur le taux d'accrétion de masse normalisé ($M_{BH}/\dot{m}$). Cette corrélation est significative dans le repère de l'observateur mais s'affaiblit ou devient non significative lorsqu'elle est corrigée pour le repère d'émission. La corrélation semble principalement pilotée par $M_{BH}$, bien que le rôle de $\dot{m}$ ne puisse être définitivement exclu en raison des limitations de la taille de l'échantillon.
• Taille de la région d'émission : Les temps de montée des éruptions les plus brillantes, utilisés comme substituts pour la largeur de la région de transition du jet, ont montré des corrélations faibles et préliminaires avec $M_{BH}$, $L_{acc}$ et $\Gamma$. Cependant, ces associations se sont révélées être pilotées par une dépendance commune au décalage vers le rouge. Lorsque le décalage vers le rouge a été pris en compte via une corrélation partielle, la significativité des relations entre les temps de montée et les paramètres du moteur central a disparu.
• Dépendance au décalage vers le rouge : Une corrélation positive significative a été trouvée entre le décalage vers le rouge et les durées d'éruptions dans le repère d'émission, tandis qu'une corrélation négative existait entre le décalage vers le rouge et les séparations d'éruptions dans le repère d'émission (impliquant un taux d'éjection plus élevé pour les sources à plus haut décalage vers le rouge). Ces tendances opposées peuvent expliquer l'absence de corrélation entre les durées/séparations d'éruptions observées et les facteurs Doppler.

Signification et affirmations
L'article affirme que les échelles de temps caractéristiques dérivées des coudes de la DSP radio sont physiquement associées au moteur central des NAG, évoluant spécifiquement avec $M_{BH}/\dot{m}$. Cela soutient l'hypothèse selon laquelle le moteur central module la variabilité du jet, similaire aux résultats trouvés dans le domaine des rayons X. L'étude suggère que la durée moyenne des éruptions les plus brillantes est le substitut de courbe de lumière le plus approprié pour l'échelle de temps du coude de la DSP.

Les auteurs notent modestement que, bien que des corrélations avec les paramètres du moteur central existent, elles sont souvent faibles et potentiellement faussées par des effets de décalage vers le rouge et des biais d'échantillon (par exemple, la sélection de sources avec des éruptions discernables). Les résultats concernant la taille de la région d'émission et le facteur de Lorentz sont décrits comme des « corrélations faibles et préliminaires » qui pourraient ne pas tenir une fois les dépendances au décalage vers le rouge éliminées. Par conséquent, l'article conclut que, bien qu'un lien avec le moteur central soit suggéré, des échantillons plus vastes et des estimations de paramètres plus précises sont nécessaires pour caractériser définitivement ces relations. Le travail ne propose pas de nouveaux dispositifs expérimentaux mais valide plutôt l'utilisation de la décomposition de courbes de lumière comme méthode pour étendre les analyses basées sur la DSP à un échantillon plus large de NAG.