The Shape of the Vertical Action Distribution Locates the Scatterers that Heat the Galactic Disc

Résumé technique : La forme de la distribution de l'action verticale localise les diffuseurs qui chauffent le disque galactique

Énoncé du problème
La dispersion des vitesses verticales des étoiles dans la Voie lactée augmente avec l'âge, un phénomène connu sous le nom de chauffage du disque. Bien que l'existence de ce chauffage soit bien établie, l'identité et la distribution spatiale des « diffuseurs » responsables restent non résolues. Deux classes principales de modèles existent :

1. Diffuseurs remplissant le volume : Des perturbations telles que la structure spirale, la barre galactique ou les sous-halos de matière noire qui agissent sur une étoile tout au long de son orbite.
2. Diffuseurs confinés au plan moyen : Des perturbations, telles que les nuages moléculaires géants (GMC), qui sont confinées dans une couche mince autour du plan galactique moyen et n'interagissent avec les étoiles que lors de leurs passages verticaux.

Les tentatives précédentes pour distinguer ces mécanismes en utilisant la relation entre l'âge et la dispersion des vitesses ou la forme de l'ellipsoïde de vitesse se sont révélées peu concluantes en raison des dégénérescences entre l'abondance, la masse des diffuseurs et l'histoire du chauffage. L'article soutient que l'action verticale ($J_z$), un invariant adiabatique, offre un diagnostic supérieur car elle étiquette l'orbite plutôt que la vitesse instantanée et est insensible à la croissance du potentiel du disque au fil du temps.

Méthodologie
Les auteurs dérivent la distribution théorique de l'action verticale pour une population coevalle d'étoiles chauffées par des diffuseurs possédant une distribution de hauteur verticale arbitraire.

• Dérivation théorique :

  • Le processus de chauffage est modélisé comme une diffusion dans l'espace de l'action verticale en utilisant l'équation de Fokker–Planck.
  • Les auteurs dérivent la diffusivité $D(J_z)$ en fonction du profil de hauteur $q(z)$ du diffuseur. Ils montrent que la diffusivité suit une loi de puissance $D \propto J_z^b$, où l'exposant $b$ est déterminé uniquement par la géométrie des diffuseurs.
  • Cas de remplissage de volume : Si les diffuseurs remplissent le volume, la diffusivité croît linéairement avec l'action ($b=1$), entraînant une distribution exponentielle $p(J_z) \propto \exp(-J_z/J_0)$.
  • Cas de la couche du plan moyen : Si les diffuseurs sont confinés dans une couche mince, le temps d'interaction diminue à mesure que la vitesse de passage de l'étoile augmente. Cela réduit la croissance de la diffusivité à $D \propto J_z^{1/2}$ ($b=1/2$), produisant une distribution d'« exponentielle étirée » $p(J_z) \propto \exp[-(J_z/J_0)^{3/2}]$.
  • La solution générale pour une population coevallle est $p(J_z) \propto \exp[-(J_z/J_0)^{2-b}]$, où la forme (l'exposant de la queue de distribution) dépend uniquement de la géométrie ($b$), tandis que l'échelle ($J_0$) encode l'amplitude et l'histoire du chauffage.

• Données et ajustement :

  • Le modèle est ajusté à un échantillon de 7 589 étoiles de la branche de la fourche rouge (red clump) à faible teneur en $\alpha$ provenant du relevé APOGEE DR14, croisé avec l'astrométrie Gaia DR2.
  • L'échantillon est restreint aux étoiles ayant des rayons effectifs compris entre $5 \le R_{\text{eff}} \le 10$ kpc et des âges compris entre $2 \le \tau \le 8$ Gyr.
  • Une composante critique de la méthodologie est le traitement rigoureux de la fonction de sélection, qui manque préférentiellement les étoiles à haute action proches de la limite de magnitude du relevé. Les auteurs appliquent une fonction de sélection par étoile dérivée de Ting & Rix (2019) directement à la vraisemblance du modèle.
  • La procédure d'ajustement marginalise les erreurs de mesure sur l'âge, l'action et le rayon en utilisant la quadrature de Gauss. Crucialement, le paramètre d'échelle $J_0$ est laissé libre pour chaque âge (paramétré par quatre nœuds et une pente radiale). Cela permet de mesurer l'exposant de forme $b$ sans supposer une histoire de chauffage ou une dispersion de naissance spécifique, permettant ainsi de découpler efficacement la géométrie du taux de chauffage.

Résultats clés

1. Exposant mesuré : L'ajustement produit une pente de diffusivité de $b = 0,51^{+0,06}_{-0,07}$.

   • Cette valeur est cohérente avec la prédiction pour une couche confinée au plan moyen ($b=0,5$) mais est fortement incompatible avec la prédiction de remplissage de volume ($b=1$).
   • Des tests synthétiques confirment que la méthode peut récupérer l'exposant correct sans biais et sépare les deux géométries avec une signification de $5,6\sigma$.
   • Le résultat est robuste à travers différents intervalles radiaux et coupes d'âge, et reste stable face aux variations des erreurs d'action supposées.

2. Masse effective des diffuseurs :

   • En utilisant l'exposant $b$ mesuré et les échelles libres ($J_0$), les auteurs calculent l'amplitude du chauffage.
   • Pour le groupe d'âge 2–4 Gyr, en comparant l'amplitude du chauffage à la densité de surface moléculaire actuelle, on obtient une masse effective des diffuseurs de $2,7 \times 10^6 M_\odot$.
   • Cette valeur se situe dans la plage ($1,4\text{--}3,4 \times 10^6 M_\odot$) dérivée de catalogues de nuages indépendants et de fonctions de masse (par exemple, Miville-Deschênes et al. 2017).
   • Pour les étoiles plus âgées, la masse effective apparente augmente si on la compare à la densité de gaz actuelle. Cependant, lorsqu'elle est corrigée par l'histoire décroissante de la formation stellaire (et donc la diminution de la densité de gaz au cours des derniers gigannées), la masse effective pour tous les groupes d'âge converge vers la plage $1,9\text{--}2,8 \times 10^6 M_\odot$.

3. Histoire du chauffage :

   • Les échelles libres croissent avec l'âge selon $J_0 \propto \tau^{1,12}$. C'est plus rapide que le $\tau^{0,67}$ attendu pour un taux de chauffage constant avec une naissance « froide », ce qui indique que le disque a été plus fortement chauffé dans le passé.

Signification et affirmations
L'article affirme résoudre le débat de longue date concernant la localisation des diffuseurs du disque galactique en utilisant la forme de la distribution de l'action verticale plutôt que simplement sa moyenne ou sa dispersion de vitesse.

• Localisation des diffuseurs : L'exposant mesuré $b \approx 0,5$ fournit une preuve directe que le mécanisme principal de chauffage vertical est localisé près du plan galactique. Cela exclut les mécanismes de remplissage de volume (tels que les sous-halos de matière noire, les ondes de flexion ou les perturbations directes des structures spirales) comme source dominante de chauffage diffusif. Les auteurs estiment que toute composante de remplissage de volume contribue à moins de 7 % au chauffage.
• Identification des diffuseurs : La combinaison de la géométrie du plan moyen et de la masse effective dérivée ($\sim 10^6 M_\odot$) identifie les nuages moléculaires géants (GMC) comme les diffuseurs. La masse est cohérente avec le haut de la fonction de masse des GMC, et la géométrie correspond à la fine couche de gaz moléculaire.
• Résolution des divergences antérieures : L'article explique pourquoi les études précédentes (par exemple, Binney & Lacey 1988) ont déduit $b=1$ (remplissage de volume). Ces études analysaient des populations d'âges mixtes ou des dispersions de vitesses sans résoudre la forme à âge fixe. Les auteurs démontrent que le mélange de populations avec différentes échelles de chauffage aplatit artificiellement la queue de la distribution, simulant une géométrie de remplissage de volume ($b \approx 1$), même si la véritable géométrie sous-jacente est une couche mince ($b=0,5$).
• Mécanisme vs Source d'énergie : Les auteurs précisent que, bien que la diffusion se produise au plan moyen (via les nuages), l' énergie peut provenir de mouvements dans le plan (par exemple, pilotés par la structure spirale) qui sont redirigés verticalement lors des passages dans les nuages. L'exposant mesure le lieu de la redirection, et non la source ultime de l'énergie.

En conclusion, l'article établit que le disque de la Voie lactée est chauffé près du plan par une population évolutive d'objets ayant la masse des nuages moléculaires géants, et que la forme de la distribution de l'action verticale est un diagnostic robuste et indépendant du modèle pour la distribution spatiale des agents de diffusion.