The dispersal of compact protoplanetary discs

Résumé technique : La dispersion des disques protoplanétaires compacts

Énoncé du problème
Les disques protoplanétaires sont les sites de la formation des planètes, mais leurs mécanismes de dispersion restent un sujet d'enquête active. Les preuves observationnelles, notamment issues des études de l'évolution couleur-couleur des populations stellaires jeunes, indiquent que la dispersion des disques progresse préférentiellement de l'intérieur vers l'extérieur, créant une cavité centrale avant que le disque externe ne soit érodé. Cependant, les modèles théoriques standards supposent souvent des disques de grands rayons externes (centaines d'unités astronomiques), bien au-delà du rayon gravitationnel ($R_g$) où la photoévaporation est la plus efficace. Des enquêtes récentes (par exemple, ALMA) révèlent qu'une fraction substantielle de disques est « compacte », avec des rayons externes comparables ou inférieurs à $R_g$. Les modèles existants, qui supposent généralement des profils de perte de masse non contraints s'étendant jusqu'à de grands rayons, ne parviennent pas à reproduire avec précision les trajectoires évolutives et les échelles de temps de dispersion de ces systèmes compacts. Plus précisément, les prescriptions standards prédisent souvent une érosion de l'extérieur vers l'intérieur pour les disques compacts, ce qui contredit les diagnostics observationnels.

Méthodologie
Les auteurs ont employé une approche multi-étapes combinant des simulations hydrodynamiques de rayonnement haute résolution avec des modèles de synthèse de population visqueuse 1D :

1. Simulations hydrodynamiques de rayonnement : En utilisant le code PLUTO couplé au cadre de transfert radiatif MOCASSIN, les auteurs ont simulé la photoévaporation induite par les rayons X pour une étoile de $1 M_\odot$. Ils ont exploré quatre rayons de coupure initiaux de disque distincts ($R_{cut}$) : 10, 50, 100 et 200 ua. Les simulations ont suivi la structure thermique, la densité et les champs de vitesse jusqu'à ce qu'un état quasi-stationnaire soit atteint. Le rayon gravitationnel a été adopté comme $R_g = 5$ ua.
2. Déduction de relations d'échelle : À partir des simulations hydrodynamiques, les auteurs ont analysé les distributions cumulatives du taux de perte de masse. Ils ont déterminé comment le taux total de perte de masse par vent dépend du rayon de coupure externe du disque.
3. Synthèse de population : Ces résultats ont été intégrés dans un code d'évolution visqueuse 1D (SPOCK). Les auteurs ont réalisé une synthèse de population de 4 096 disques synthétiques en utilisant une stratégie d'échantillonnage Quasi-Monte Carlo (séquences de Sobol) pour explorer un espace de paramètres à 9 dimensions (masse stellaire, luminosité en rayons X, masse initiale du disque, rayon de coupure, viscosité, champ FUV externe, etc.).
4. Prescriptions testées : Les modèles ont comparé deux scénarios :
   • Prescription standard : Taux de perte de masse par photoévaporation interne appliqués sans égard au rayon externe fini du disque (supposant effectivement que le disque s'étend à l'infini ou au-delà du pic de perte de masse).
   • Prescription dépendante de la coupure : Le taux de perte de masse est tronqué et normalisé en fonction du rayon de coupure spécifique du disque, dérivé des relations d'échelle hydrodynamiques.
   • Photoévaporation externe : Certains modèles incluaient une photoévaporation externe FUV modérée (en utilisant la grille FRIED) pour simuler des étoiles OB voisines.

Contributions et résultats clés

• Échelle hydrodynamique : Les simulations ont révélé que les profils de perte de masse de surface sont presque indépendants de la taille du disque ; la morphologie du vent et les taux de perte de masse locaux dans les régions internes restent cohérents indépendamment de la coupure externe. Cependant, le taux total de perte de masse intégré est considérablement réduit pour les disques compacts car le vent est tronqué au rayon de coupure. Le taux de perte de masse cumulatif évolue avec le rayon, croissant comme $\dot{M}_w \propto R^3$ près du rayon gravitationnel et s'aplatissant vers $\dot{M}_w \propto R^{0.5}$ à des rayons plus grands.
• Trajectoires évolutives :
  • Prescription standard : Lorsqu'elle est appliquée à des disques compacts ($R_{cut} < 20-30$ ua), le modèle standard prédit que le taux de perte de masse par vent dépasse le taux d'étalement visqueux en dehors du rayon de coupure. Cela conduit à une dispersion de l'extérieur vers l'intérieur, où les bords externes sont érodés en premier, contredisant les preuves observationnelles.
  • Prescription dépendante de la coupure : Lorsque le taux de perte de masse est mis à l'échelle du rayon de coupure, le taux total de vent est réduit suffisamment pour permettre à l'étalement visqueux de dominer initialement. Cela résulte en une dispersion de l'intérieur vers l'extérieur, où une ouverture se forme près du rayon gravitationnel et s'étend vers l'extérieur, cohérent avec les diagnostics observationnels.
• Résultats de la synthèse de population :
  • La prescription dépendante de la coupure, combinée à un taux de photoévaporation interne réduit (mis à l'échelle par un facteur de 10 pour s'aligner avec les modèles chimiques récents), reproduit les fractions de disques observées et les taux d'accrétion en fonction de l'âge de l'amas sans ajustement de paramètres.
  • La durée de vie médiane du disque dans le modèle amélioré est de $2,70 \pm 0,13$ Myr, ce qui correspond bien aux estimations observationnelles ($2,24 \pm 0,35$ Myr) et est significativement plus longue que les prédictions des modèles sans réduction de taux ($1,13 \pm 0,05$ Myr).
  • L'inclusion d'une photoévaporation externe modérée n'empêche pas l'étalement du disque lorsque la prescription du rayon de coupure est utilisée, mais entraîne une augmentation plus douce des rayons des disques au fil du temps et une fraction plus élevée de disques compacts persistant jusqu'à des âges plus avancés.

Importance
L'article soutient que la dispersion des disques protoplanétaires compacts ne peut pas être correctement capturée par des modèles standards de photoévaporation non mis à l'échelle. Les auteurs proposent une relation d'échelle simple et motivée physiquement qui tient compte du rayon externe fini du disque. En intégrant cette dépendance, les modèles parviennent à réconcilier avec succès les prédictions théoriques avec le nettoyage observé de l'intérieur vers l'extérieur des populations stellaires jeunes et la prévalence des disques compacts. L'étude conclut que les prescriptions de dispersion doivent inclure explicitement la dépendance au rayon de coupure du disque pour prédire correctement l'évolution et la durée de vie des disques compacts, fournissant une correction nécessaire aux modèles actuels de synthèse de population.