Drifters on the edge of town: $λ$ Boötis stars in clusters

Résumé technique : Des dériveurs en bordure de ville : les étoiles $\lambda$ Boötis dans les amas

Énoncé du problème
Les étoiles $\lambda$ Boötis ($\lambda$ Boo) sont des étoiles de type A chimiquement particulières, caractérisées par des abondances d'éléments légers (C, N, O, S) similaires à celles du Soleil, mais par des déficiences significatives en éléments du groupe du fer. L'explication théorique dominante postule que ces étoiles accrètent du gaz pur et pauvre en métaux provenant du milieu interstellaire (MIS). Cependant, ce scénario fait face à un défi majeur : les découvertes récentes d'étoiles $\lambda$ Boo au sein ou aux abords d'amas stellaires contredisent l'attente selon laquelle le rayonnement ionisant intense provenant d'étoiles intermédiaires et massives (> 5 $M_\odot$) dans de tels environnements empêcherait l'accrétion de gaz pur du MIS ou détruirait les disques circumstellaires nécessaires au processus. De plus, il existe une tension entre les taux d'accrétion requis pour produire les abondances $\lambda$ Boo ($10^{-14}$ à $10^{-9}$ $M_\odot$ an⁻¹) et les taux beaucoup plus faibles déduits des observations de naines blanches (< $10^{-17}$ $M_\odot$ an⁻¹). Le problème spécifique abordé est de savoir si les histoires dynamiques des étoiles de type A dans les amas en évolution leur permettent d'échapper au cœur irradié de l'amas, d'accréter du gaz pur, puis de revenir, réconciliant ainsi la présence d'étoiles $\lambda$ Boo dans les environnements d'amas avec l'hypothèse d'accrétion du MIS.

Méthodologie
Les auteurs utilisent des simulations N-corps purement gravitationnelles pour modéliser la formation et l'évolution dynamique de régions de formation d'étoiles sur 100 Myr.

• Configuration de la simulation : Quatre séries de simulations ont été menées avec des nombres totaux d'étoiles ($N_*$) de 150, 500, 1000 et 3000. Les conditions initiales comprenaient une sous-structure spatiale et cinématique générée par la méthode fractale de boîte (dimension fractale $D=1,6$) et un rapport virial de $\alpha_{vir} = 0,3$, conduisant à un effondrement initial et à la formation d'amas concentrés au centre. Les masses stellaires ont été échantillonnées à partir d'une fonction de masse initiale (IMF) de Maschberger (2013) avec une limite supérieure de 50 $M_\odot$. Les binaires primordiales ont été exclues.
• Analyse post-traitement : Les simulations ont été analysées a posteriori pour déterminer :
  1. Potentiel d'accrétion : Les étoiles de type A ont été suivies pour identifier quand elles se déplaçaient au-delà du rayon de Jacobi (tidal) ($r_J$), défini comme la limite où elles pourraient rencontrer du gaz pur du MIS. Lorsqu'elles se trouvaient au-delà de $r_J$, les taux d'accrétion de Bondi-Hoyle-Lyttleton (BHL) étaient calculés en fonction de la vitesse de l'étoile par rapport à un nuage de MIS stationnaire ($\rho = 10$ cm⁻³).
  2. Environnement de rayonnement : Les flux d'UV lointains (FUV) et d'UV extrêmes (EUV) incidents sur chaque étoile de type A ont été calculés en fonction de la distance aux étoiles massives (> 5 $M_\odot$ et > 20 $M_\odot$, respectivement).
  3. Survie du disque : En utilisant la grille FRIED de modèles de photoévaporation et des prescriptions analytiques, les auteurs ont modélisé la perte de masse des disques circumstellaires des étoiles de type A. Un disque était considéré comme « détruit » si la masse totale perdue dépassait la masse initiale du disque (supposée être de 10 % de la masse stellaire). Une étoile était considérée comme candidate à la formation de $\lambda$ Boo uniquement si elle accrétait du matériel du MIS tout en conservant son disque.
  4. Contraintes temporelles : Les calculs ont été répétés avec une contrainte empêchant l'accrétion après 40 Myr pour tenir compte de l'enrichissement potentiel du MIS par les supernovae de type Ia.

Résultats clés

• Histoires dynamiques : Dans toutes les réalisations de simulation, la majorité (90–98 %) des étoiles de type A ont voyagé au-delà du rayon de Jacobi à un moment donné au cours de l'évolution de 100 Myr. Cela s'est produit à la fois par éjection du cœur de l'amas et par « dérive », où les étoiles se sont déplacées à l'extérieur du rayon de marée pour ensuite revenir dans les régions internes.
• Taux d'accrétion : Les étoiles de type A se déplaçant au-delà du rayon de Jacobi ont atteint des taux d'accrétion BHL allant de $10^{-12}$ à $10^{-5}$ $M_\odot$ an⁻¹. Les étoiles éjectées à des vitesses élevées avaient des taux plus faibles ($10^{-12}–10^{-10}$ $M_\odot$ an⁻¹), tandis que celles sur des orbites centré sur l'amas qui dérivaient vers l'extérieur puis revenaient avaient des taux plus élevés ($10^{-9}–10^{-5}$ $M_\odot$ an⁻¹), chevauchant la plage requise pour expliquer les abondances $\lambda$ Boo.
• Survie du disque vs Rayonnement : La survie des disques circumstellaires dépendait fortement du nombre total d'étoiles (et donc du nombre de sources ionisantes) dans l'amas.
  • Dans les régions de faible masse ($N_* = 150$), environ 55 % des étoiles de type A ayant exploré au-delà du rayon de Jacobi ont conservé leurs disques pendant plus de 10 Myr.
  • Dans les régions de haute masse ($N_* = 3000$), 0 % des étoiles de type A ont conservé leurs disques en raison des champs FUV/EUV intenses, excluant de fait la formation d'étoiles $\lambda$ Boo dans ces environnements.
  • Dans les régions intermédiaires ($N_* = 500, 1000$), une fraction des étoiles de type A (8,8 % et 4,8 % respectivement) a conservé ses disques lors de l'accrétion.
• Systèmes binaires : L'étude a constaté que des systèmes binaires larges (> 100 au) contenant des étoiles de type A peuvent se former par capture lors de la dissolution des régions de formation d'étoiles. Crucialement, une étoile peut accréter du matériel pur du MIS, développer des abondances $\lambda$ Boo, et ensuite capturer un compagnon (ou vice versa). Ce mécanisme permet des systèmes binaires où les composantes présentent des abondances chimiques significativement différentes. La distribution de séparation simulée de ces binaires intersecte les binaires larges observées contenant des étoiles $\lambda$ Boo (par exemple, HD 87304 + CD-33 6615B).

Signification et affirmations
L'article affirme démontrer que la présence d'étoiles $\lambda$ Boötis à proximité ou dans des amas stellaires n'invalide pas intrinsèquement le scénario d'accrétion du MIS. Les auteurs soutiennent que :

1. L'évasion dynamique est suffisante : Les étoiles de type A n'ont pas besoin d'être éjectées de façon permanente pour accréter du gaz pur ; elles peuvent dériver au-delà du rayon de marée, accréter, puis revenir au centre de l'amas.
2. Seuils de rayonnement : La viabilité du scénario d'accrétion est fortement corrélée à la masse de l'amas. Les amas de faible masse (similaires à Theia 139 et HSC 1640, où des étoiles $\lambda$ Boo ont été récemment trouvées) possèdent des champs de rayonnement suffisamment faibles pour permettre la survie des disques, tandis que les amas massifs ne le permettent pas.
3. Diagnostics binaires : L'observation de systèmes binaires avec des abondances disparates entre les composantes n'invalide pas l'hypothèse d'accrétion. De tels systèmes peuvent se former dynamiquement après l'événement d'accrétion.

Les auteurs restent modestes quant à la validité absolue des taux d'accrétion, notant que la tension entre les taux requis et les limites des naines blanches reste non résolue. Ils concluent que, bien que leurs modèles dynamiques soutiennent la possibilité de la formation d'étoiles $\lambda$ Boo dans les amas, des recherches supplémentaires sont nécessaires pour valider les taux d'accrétion spécifiques et l'interaction physique entre les disques en accrétion et les champs de rayonnement.