Characterizing the Scale Height and Filamentary Structure of Radiatively Cooled MADs

Résumé technique : Caractérisation de l'échelle de hauteur et de la structure filamentaire des MAD refroidis par rayonnement

Énoncé du problème
Le refroidissement radiatif est un composant fondamental de la dynamique réaliste des disques d'accrétion autour des trous noirs, mais son impact spécifique sur la structure des disques arrêtés magnétiquement (MAD) reste sous-exploré. Bien que les simulations antérieures de magnétohydrodynamique en relativité générale (GR-MHD) aient étudié les MAD de manière approfondie, la plupart se sont concentrées sur des écoulements idéaux, non radiatifs, adaptés aux sources de luminosité extrêmement faible. Les études existantes intégrant le rayonnement ont souvent été limitées par la complexité computationnelle et la dimensionalité de l'espace des paramètres. De plus, un problème critique surgit dans la caractérisation de l'épaisseur verticale des MAD : la définition traditionnelle de l'échelle de hauteur du disque suppose une distribution de densité symétrique par rapport au plan équatorial. Cependant, dans les états MAD, l'accumulation de flux magnétique intense entraîne des instabilités d'échange, créant des courants d'accrétion discrets et filamentaires qui s'écartent souvent significativement du plan équatorial. Par conséquent, les métriques standard peuvent produire des résultats trompeurs concernant l'épaisseur géométrique réelle et l'état dynamique du disque, en particulier lorsque le refroidissement radiatif modifie le soutien de la pression thermique.

Méthodologie
Les auteurs ont réalisé une série de simulations GR-MHD tridimensionnelles à l'aide du code accéléré par GPU cuHARM autour d'un trou noir non tournant (Schwarzschild) de masse $M = 4 \times 10^6 M_\odot$. Les simulations ont modélisé des états MAD avec des taux d'accrétion de masse ($\dot{M}$) variant de $10^{-7}$ à $10^{-4} \dot{M}_{\text{Edd}}$.

Les composants méthodologiques clés incluent :

• Refroidissement radiatif : Les simulations ont systématiquement inclus les pertes radiatives via l'émission synchrotron et le rayonnement de freinage. Le refroidissement Compton a été exclu car il a été jugé subdominant à ces taux d'accrétion et nécessite un transport complet du rayonnement.
• Thermodynamique : Une prescription à température unique a été adoptée, supposant un rapport fixe entre les températures des protons et des électrons ($T_p/T_e = 3$), impliquant un couplage fort entre les espèces. Le taux de refroidissement $\Lambda$ a été calculé en fonction de la densité locale du gaz, de la température des électrons et de la densité d'énergie du champ magnétique.
• Conditions initiales : Les simulations ont commencé avec un tore axisymétrique non magnétisé en équilibre hydrostatique, perturbé pour déclencher l'instabilité magnéto-rotationnelle (MRI). Une configuration spécifique de champ magnétique (potentiel vecteur en boucle simple) a été utilisée pour assurer le développement d'un état MAD.
• Analyse : Des diagnostics moyennés dans le temps ont été calculés sur l'intervalle $1.5 \times 10^4 t_g < t < 3.0 \times 10^4 t_g$. Les auteurs ont comparé l'échelle de hauteur angulaire pondérée par la densité standard à une nouvelle métrique proposée basée sur la position polaire des maxima de densité locaux afin de mieux capturer les structures filamentaires.

Contributions et résultats clés

1. Domination du refroidissement synchrotron :
   L'étude confirme que l'émission synchrotron est le mécanisme de refroidissement radiatif dominant sur toute la gamme des taux d'accrétion simulés, dépassant les pertes par rayonnement de freinage de plusieurs ordres de grandeur. Le refroidissement synchrotron atteint son maximum dans les régions internes du disque près du trou noir et le long du plan équatorial, entraîné par des champs magnétiques intenses, des densités et des températures élevées.

2. Efficacité radiative et structure du disque :
   L'efficacité radiative ($\eta_{\text{rad}}$) augmente de manière non linéaire avec le taux d'accrétion de masse. À de faibles taux ($\dot{M} \lesssim 10^{-6} \dot{M}_{\text{Edd}}$), l'écoulement se comporte comme un écoulement d'accrétion radiativement inefficace (RIAF). Cependant, à des taux plus élevés (approchant $10^{-4} \dot{M}_{\text{Edd}}$), le refroidissement devient suffisamment efficace pour réduire significativement la pression thermique. Cela conduit à une transition où le disque devient radiativement efficace (atteignant une efficacité d'environ 30 % à $r=5r_g$ pour le taux le plus élevé) et où le gaz est comprimé en filaments plus minces et plus denses.

3. Échec des définitions traditionnelles de l'échelle de hauteur :
   Les auteurs démontrent que la définition conventionnelle de l'échelle de hauteur pondérée par la densité (Équation 14) est trompeuse pour les MAD. Parce que les écoulements MAD présentent des structures filamentaires déplacées du plan équatorial en raison des instabilités d'échange magnétique, la moyenne sur le plan équatorial gonfle artificiellement la hauteur mesurée du disque. La métrique traditionnelle échoue à distinguer entre l'épaisseur physique du gaz accrétant et le déplacement vertical des filaments.

4. Proposition d'une nouvelle métrique d'échelle de hauteur :
   Pour remédier à cela, l'article introduit une définition alternative de la hauteur du disque basée sur l'emplacement des maxima de densité locaux au sein de tranches azimutales (Équation 15). Cette méthode identifie des « îlots » cohérents de haute densité, exclut les régions coronales de faible densité et calcule l'épaisseur de ces filaments par rapport à leurs centres locaux.

   • Résultat : En utilisant cette nouvelle métrique, la hauteur du disque est montrée comme diminuant de manière monotone vers le trou noir, s'approchant de zéro à l'horizon. En revanche, la métrique traditionnelle montrait une augmentation artificielle de l'épaisseur près de l'horizon.
   • Impact du refroidissement : La nouvelle métrique révèle que pour des taux d'accrétion $\dot{M} \gtrsim 10^{-4} \dot{M}_{\text{Edd}}$, le refroidissement radiatif provoque une réduction substantielle de l'échelle de hauteur du disque à l'intérieur des $\sim 30-40 r_g$ internes, tandis que les taux d'accrétion plus faibles ne montrent aucun tel amincissement.

5. Taux d'accrétion critique :
   Les auteurs dérivent une condition analytique pour le taux d'accrétion de masse critique au-dessus duquel le refroidissement radiatif devient dynamiquement significatif (plus rapide que le chauffage). En équilibrant le refroidissement synchrotron avec le chauffage magnétique (médié par la turbulence et la reconnexion), ils estiment que cette transition se produit autour de $\dot{M} \sim 2.17 \times 10^{-4} \dot{M}_{\text{Edd}}$ à $r=30r_g$. Cette estimation analytique correspond à l'observation numérique des changements structurels dans la simulation S-4 ($10^{-4} \dot{M}_{\text{Edd}}$).

Signification et affirmations
L'article affirme que le refroidissement radiatif modifie fondamentalement la structure interne des MAD en réduisant le soutien de la pression thermique, comprimant ainsi les filaments d'accrétion et réduisant l'échelle de hauteur du disque. La signification principale de ce travail réside dans :

• Correction du biais diagnostique : Il établit que les mesures conventionnelles de l'épaisseur du disque sont inadéquates pour les MAD, en particulier en présence d'un refroidissement intense, et propose une définition plus robuste, consciente des filaments.
• Identification d'un régime de transition : Il identifie un seuil de taux d'accrétion critique où le refroidissement passe d'une perturbation mineure à une force dominante remodelant la géométrie et la dynamique du disque.
• Implications observationnelles : Les auteurs suggèrent que ces changements structurels induits par le refroidissement (filaments plus minces et plus denses) devraient se marquer sur les propriétés observables, telles que les distributions d'énergie spectrales et les fluctuations des courbes de lumière, offrant potentiellement une nouvelle méthode pour déduire les taux d'accrétion de masse dans des systèmes comme Sgr A* ou les binaires X.

L'étude reste modeste dans son ampleur, limitant l'analyse aux trous noirs non tournants et supposant un rapport fixe entre les températures des ions et des électrons, notant que des travaux futurs sont nécessaires pour explorer l'impact de la rotation des trous noirs et de modèles de couplage électron-proton plus complexes à des taux d'accrétion plus élevés.