The FUor Mass Distribution Matches the Solar Neighborhood IMF: Evidence for a Universal Eruptive Phase

1. Problématique et Contexte

La formation des étoiles implique une phase d'accrétion de masse, mais la chronologie de ce processus reste mal comprise. Les taux d'accrétion observés pour les objets stellaires jeunes (YSO) de classes I et II sont généralement trop faibles ($10^{-8} - 10^{-7} M_{\odot} \text{an}^{-1}$) pour expliquer la masse finale des étoiles sur la séquence principale, compte tenu de la durée de vie limitée des disques protostellaires (< $10^7$ ans).

Les objets FU Orionis (FUors) représentent une solution potentielle à ce paradoxe. Ce sont des étoiles en éruption qui subissent des augmentations extrêmes de leur taux d'accrétion ($10^{-5} - 10^{-4} M_{\odot} \text{an}^{-1}$) durant des décennies. Cependant, durant ces éruptions, le disque d'accrétion est environ 100 fois plus lumineux que l'étoile centrale, rendant l'observation directe de la masse stellaire impossible. La question centrale est donc de déterminer si toutes les jeunes étoiles subissent une phase d'éruption de type FUor, ce qui résoudrait la tension entre les taux d'accrétion observés et les masses finales attendues. Pour cela, il est nécessaire de caractériser la distribution de masse de la population des FUors.

2. Méthodologie

L'étude repose sur une analyse spectrale à haute résolution dans le proche infrarouge (NIR) d'un échantillon de 28 FUors.

• Données Observations :

  • Keck/NIRSPEC : 23 FUors observés entre 2022 et 2023, couvrant les bandes Y, J, H et K (1,0–2,4 µm) avec une résolution $R \approx 15\,000 - 25\,000$.
  • IGRINS (Archives) : 15 spectres supplémentaires (bandes H et K) pour 8 sources, avec une résolution plus élevée ($R \approx 45\,000$).
  • L'échantillon total permet une couverture spectrale de 1,2 à 2,4 µm pour la majorité des objets.

• Analyse Spectrale et Rotation Keplerienne :

  • Les auteurs exploitent le fait que le disque interne des FUors est en rotation keplerienne ($v \propto r^{-1/2}$). En raison du profil de température du disque ($T \propto r^{-3/4}$), les longueurs d'onde plus rouges (K) proviennent de régions plus froides et plus éloignées (vitesses plus faibles), tandis que les longueurs d'onde plus bleues (Y/J) proviennent de régions plus chaudes et proches (vitesses plus élevées).
  • Technique de Corrélation Croisée (CCF) : Au lieu d'analyser des raies individuelles (souvent confondues par des effets de vent ou de température), les auteurs utilisent la fonction de corrélation croisée (CCF) sur l'ensemble des bandes spectrales. Ils comparent les spectres observés à un modèle de disque d'accrétion de référence (basé sur FU Ori) généré avec une grille de vitesses maximales ($v_{max}$) et d'inclinaisons ($i$).
  • Focus sur la bande H : La bande H (1,5–1,8 µm) est identifiée comme optimale car elle offre un bon compromis entre la sensibilité aux objets enfouis, la richesse des raies atomiques et moléculaires, et une moindre contamination par l'absorption des vents ou l'extinction par rapport aux bandes Y/J/K.

• Simulation de la Distribution de Masse :

  • Les auteurs simulent la distribution attendue de $v_{max}$ pour différents scénarios de fonction de masse initiale (IMF) :
    1. Une IMF « plate » (uniforme).
    2. L'IMF du voisinage solaire (Kirkpatrick et al. 2024, K24), testée avec et sans naines brunes (BD).
  • Les simulations intègrent une isochrone de 1 Myr (PARSEC) pour relier la masse stellaire au rayon, et supposent que le rayon d'interruption interne du disque ($R_{inner}$) est égal au rayon stellaire ($R_*$) durant l'éruption.

3. Résultats Clés

• Preuve de la Rotation Keplerienne Différentielle :
  L'analyse des CCFs démontre clairement que, pour la majorité des FUors, les profils de raies dans la bande K sont plus étroits que dans les bandes Y et J. Cela confirme que le disque interne est en rotation keplerienne et que les différentes longueurs d'onde sondent des rayons différents du disque. Des exceptions (absence de différentiel) sont expliquées par des modèles de disques avec un rayon externe ($R_{outer}$) très petit, ce qui limite la gamme de vitesses sondée.

• Mesure de $v_{max}$ :
  Les auteurs ont mesuré la vitesse maximale projetée ($v_{max} = \sqrt{GM_*/R_{inner}} \sin i$) pour 28 FUors. La distribution observée présente une bimodalité avec deux pics distincts : l'un autour de $\sim 30 \text{ km s}^{-1}$ et l'autre autour de $\sim 120 \text{ km s}^{-1}$.

• Comparaison avec les IMF Simulées :

  • L'IMF « plate » est rejetée car elle surestime le nombre d'objets à très haute vitesse.
  • L'IMF K24 incluant les naines brunes surestime les objets à très basse vitesse (ce qui suggère une incompréhension observationnelle des FUors de très faible masse).
  • Résultat Principal : L'IMF K24 restreinte à la gamme stellaire ($0,1 - 2 M_{\odot}$) correspond remarquablement bien à la distribution observée de $v_{max}$, reproduisant à la fois la bimodalité et la densité d'objets.

• Rayons Stellaires :
  Les simulations montrent que l'hypothèse d'étoiles « gonflées » (rayons augmentés de 20 % en raison de l'accrétion intense) n'améliore pas significativement l'ajustement par rapport aux modèles standards, suggérant que tous les FUors ne sont pas nécessairement fortement gonflés.

4. Contributions et Signification

• Validation de l'IMF Universelle pour les FUors :
  La découverte que la distribution de masse des FUors correspond à celle du voisinage solaire (IMF standard) est une preuve forte que toutes les jeunes étoiles, et pas seulement une sous-classe, subissent une phase d'éruption de type FUor au début de leur évolution pré-séquence principale. Cela résout le problème de l'accrétion de masse : les étoiles accumulent la majeure partie de leur masse durant ces éruptions brèves mais intenses.

• Nouvelle Méthodologie Spectrale :
  L'article établit une technique robuste utilisant la corrélation croisée sur la bande H pour mesurer les vitesses de rotation dans les disques d'accrétion, contournant les difficultés liées à l'observation directe de l'étoile centrale et aux effets de vent.

• Implications pour les Naines Brunes :
  L'absence apparente de FUors de type naine brune dans l'échantillon actuel pourrait être due à un biais d'observation (faible luminosité d'éruption) plutôt qu'à une absence physique. Des relevés futurs à haute sensibilité pourraient révéler cette population.

• Perspectives Futures :
  L'étude souligne la nécessité d'observations à très haute résolution (avec des instruments comme CRIRES ou HISPEC) pour lever les ambiguïtés liées à la résolution instrumentale et à la microturbulence, permettant ainsi de mieux contraindre la masse des FUors de faible masse et d'inclinaison faible.

En conclusion, ce travail fournit la première contrainte statistique sur la masse des FUors, validant l'hypothèse d'une phase d'éruption universelle dans la formation stellaire et reliant directement les propriétés des disques d'accrétion en éruption à la fonction de masse initiale de la population stellaire.