Polycyclic aromatic hydrocarbon (PAH) abundances in the disk around T Chamaeleontis (T Cha): PAH sizes, ionization fraction, and mass during JWST observations

Résumé technique : Abondances des hydrocarbures aromatiques polycycliques (HAP) dans le disque autour de T Chamaeleontis

Problématique et contexte
Les hydrocarbures aromatiques polycycliques (HAP) sont largement acceptés comme étant les porteurs des bandes infrarouges aromatiques (BIA) observées dans les environnements astrophysiques. Dans les disques protoplanétaires (DPP), les HAP sont dynamiquement couplés au gaz et aux surfaces de bord du disque, jouant un rôle critique dans le chauffage du gaz, la thermochimie et la photoévaporation. Cependant, leurs abondances dans les DPP restent mal contraintes. Bien que les BIA soient fréquemment détectées dans les disques d'Herbig Ae/Be, elles sont rares dans les disques de type T Tauri, ce qui est souvent attribué au flux UV plus faible des étoiles T Tauri ou à l'appauvrissement des HAP via la formation d'amas et l'adsorption sur des grains plus gros. T Chamaeleontis (T Cha), une étoile T Tauri présentant un disque de transition connu avec un gap de poussière, constitue un cas unique où des BIA proéminentes ont été détectées lors de récentes observations du télescope spatial James Webb (JWST). Le défi consiste à caractériser la population de HAP (tailles, fraction d'ionisation et masse) responsable de ces caractéristiques et à comprendre leur distribution au sein de la structure du disque.

Méthodologie
Les auteurs ont utilisé des données d'archives de l'instrument de spectroscopie à moyenne résolution (MRS) du MIRI du JWST pour T Cha, ainsi que la photométrie optique et les flux dans les bandes millimétriques provenant d'ALMA, SEST et ATCA. L'étude a employé un modèle de transfert radiatif de poussière 3D paramétrique à l'aide du code MCFOST. Les étapes méthodologiques clés comprenaient :

1. Décomposition spectrale : Les BIA observées à 6,2, 8,1 et 11,3 µm ont été ajustées à l'aide de profils de Drude afin d'identifier les composantes spectrales sous-jacentes. Cette analyse a révélé un profil de classe C avec certaines contributions de classe B, distinct des profils de classe A souvent utilisés dans les modèles standards.
2. Modification de l'opacité : Les opacités « astro-HAP » standard (Li & Draine 2001), qui reproduisent principalement les profils de classe A/B, étaient insuffisantes. Les auteurs ont modifié le formalisme de l'opacité des HAP en se basant sur leur analyse de profil de Drude pour reproduire les caractéristiques spécifiques de classe C observées dans T Cha.
3. Modélisation du disque : Un modèle de disque à plusieurs composantes a été construit, incluant un disque interne (avec un bord bombé et une région ombragée), un gap de poussière et un disque externe. Le disque externe a été subdivisé en un « mur » (bord interne) et une « région externe » pour tenir compte des variations de taille de grains.
4. Chauffage stochastique : Les HAP ont été traitées comme des grains de poussière subissant un chauffage stochastique. Le modèle faisait varier les distributions de taille des HAP (nombre d'atomes de carbone, $N_C$), les fractions d'ionisation ($\phi$) et le rapport de masse HAP/poussière ($M_{HAP}/M_{poussière}$) pour ajuster simultanément le continuum infrarouge moyen (MIR), les flux/rapports des BIA, ainsi que les flux optiques et millimétriques.
5. Optimisation : Une grille de modèles a été explorée en utilisant la minimisation du $\chi^2$ et l'inspection visuelle pour contraindre les paramètres, réduisant les dégénérescences en intégrant les contraintes issues des études antérieures par Spitzer, VLT et ALMA.

Contributions clés et résultats
L'étude a réussi à reproduire le spectre de T Cha, incluant le continuum MIR et les BIA proéminentes, en obtenant les résultats spécifiques suivants :

• Caractéristiques de la population de HAP : Le meilleur modèle prédit une population de petits HAP contenant environ 10 à 26 atomes de carbone ($N_C \lesssim 26$). La fraction d'ionisation est estimée à $\phi \approx 0,15$, indiquant que la population de HAP est majoritairement neutre (~85 %).
• Abondance de HAP : Le rapport de masse HAP/poussière dérivé est de $\sim 7 \times 10^{-3}$, ce qui correspond à environ 17 % de la valeur du milieu interstellaire (ISM). Il s'agit d'une valeur élevée comparée à d'autres disques modélisés.
• Morphologie du disque et distribution des grains :
  • Structure du disque externe : Le modèle nécessite un mur de disque externe (31–36 au) composé de grains de poussière plus petits (jusqu'à 4 µm) pour ajuster la pente raide du continuum entre 14 et 25 µm. La région externe (36–50 au) contient des grains plus gros (jusqu'à 5 mm) avec une distribution de taille plus faible ($p = -3,3$), suggérant une croissance efficace des grains.
  • Remplissage du gap : Pour expliquer le plateau de continuum plus faible autour de 10 µm dans le spectre du JWST (comparé au creux plus profond des données de Spitzer), le modèle propose la présence de grains de poussière submicroniques (0,01–1 µm) au sein du gap de poussière (1–31 au). Ces grains sont probablement couplés au gaz en accrétion.
  • Disque interne : La masse de poussière du disque interne a été réduite par rapport aux modèles précédents pour correspondre à la diminution du continuum de 5–9 µm, ce qui est cohérent avec la variabilité de type « see-saw » observée pour T Cha.
• Classification spectrale : Le profil des BIA dans T Cha est qualitativement caractérisé comme un profil de classe C avec certaines signatures de classe B, impliquant un rapport aliphatique-sur-aromatique plus élevé ou des contributions spécifiques de sous-groupes aliphatiques par rapport aux environnements hautement traités par les UV.

Signification et affirmations
Les auteurs affirment que leur modélisation fournit une estimation cohérente de la présence de HAP dans un disque de type T Tauri, comblant une lacune dans la littérature où de telles caractérisations détaillées sont rares.

• Implications physiques : Le rapport de masse HAP/poussière élevé et la présence de petits HAP majoritairement neutres à la surface du disque suggèrent que la photoévaporation par FUV est très efficace dans le système T Cha. Les auteurs proposent que cela pourrait engendrer des taux de perte de masse de l'ordre de $10^{-9}$ à $10^{-8} M_{\odot} \text{an}^{-1}$.
• Évolution du disque : La présence de petits grains dans le gap soutient les théories de réapprovisionnement du disque interne suite à des événements de destruction, tandis que la stratification de la taille des grains dans le disque externe indique une croissance active des grains.
• Limites : L'article reconnaît modestement que les modèles présentent des dégénérescences, particulièrement concernant l'indice d'épanouissement (flaring index) et la distribution exacte de la taille des HAP. Les auteurs notent que l'estimation des abondances de HAP repose sur l'hypothèse que les opacités modifiées capturent avec précision la physique spécifique de la classe C et que l'ajustement du continuum dans la plage 9–14 µm n'était pas parfait. Ils concluent que bien que leurs résultats soient robustes dans l'espace de paramètres actuel, une modélisation simultanée future de l'émission des raies et de la poussière pourrait réduire davantage ces dégénérescences.