Formation of Isotopically Heterogeneous Molecular Cloud Cores in Filamentary Molecular Clouds

Résumé technique : Formation de cœurs de nuages moléculaires isotopiquement hétérogènes dans les nuages moléculaires filamentaires

Énoncé du problème
L'analyse des météorites révèle que les solides du système solaire primitif, tels que les inclusions riches en calcium et aluminium (CAI), présentent une inhomogénéité isotopique (par exemple, des variations des rapports $^{54}\text{Cr}/^{52}\text{Cr}$ et $^{50}\text{Ti}/^{47}\text{Ti}$). Bien que le traitement thermique dans les disques protoplanétaires et l'héritage d'inhomogénéités des cœurs de nuages moléculaires parents soient des explications proposées, l'étendue de la survie de l'hétérogénéité isotopique lors de l'effondrement des nuages moléculaires filamentaires en cœurs denses reste incertaine. Plus précisément, il est nécessaire de quantifier quelle proportion de l'hétérogénéité isotopique initiale présente dans un filament est préservée ou lissée lors de la fragmentation gravitationnelle et de la formation des cœurs, et comment cette hétérogénéité se traduit dans le disque circumstellaire résultant.

Méthodologie
Les auteurs ont utilisé des simulations de la méthode de la SPH (Smoothed Particle Hydrodynamics), employant spécifiquement la méthode GSPM (Godunov Smoothed Particle Magnetohydrodynamics) pour modéliser la fragmentation gravitationnelle de nuages moléculaires filamentaires.

• Configuration de la simulation : Les simulations ont modélisé des filaments avec une masse linéique critique ($M_{\text{line,crit}} \approx 18 M_\odot \text{pc}^{-1}$) et une largeur d'environ $0,1$ pc, conformément aux observations. Les conditions initiales comprenaient un profil de densité en équilibre hydrostatique, une turbulence de Kolmogorov ($\sigma = 2c_s$), et un champ magnétique initial de $10 \mu\text{G}$ perpendiculaire à l'axe du filament. Des cas purement hydrodynamiques ont également été exécutés à des fins de comparaison.
• Traçage isotopique : Les espèces isotopiques ont été traitées comme des traceurs passifs. Les auteurs ont construit six modèles distincts (Tableau 1) pour imposer des gradients de rapport isotopique initiaux ($C_0$) selon différents axes : l'axe longitudinal ($z$), les axes mineurs ($x, y$), et des combinaisons de ceux-ci. Le rapport isotopique a été normalisé entre 0 et 1 sur le domaine de calcul.
• Analyse : L'évolution du rapport isotopique a été suivie en suivant les particules SPH. Les auteurs ont défini des couches sphériques autour du pic de densité des cœurs en formation pour calculer le rapport isotopique moyen ($\bar{C}$) et la différence par rapport au centre du cœur ($|\Delta C|$). Ils ont également employé une approche semi-analytique pour estimer le déplacement initial du centre de masse des couches ($\delta z_{\text{ini}}$) basé sur les champs de vitesse turbulents et les spectres de puissance de densité.

Résultats clés

1. Survie de l'hétérogénéité : Les simulations démontrent que les inhomogénéités isotopiques ne sont pas complètement lissées lors de la fragmentation du filament. Les cœurs conservent environ 0,1 % à 1 % de l'inhomogénéité isotopique initiale présente sur la longueur initiale du filament de 1 pc.
2. Dépendance géométrique : L'ampleur de l'hétérogénéité conservée dépend de l'orientation du gradient isotopique initial par rapport à la géométrie du filament.
   • Les gradients le long de l'axe longitudinal (direction $z$, Modèle $M_z$) entraînent les plus grandes différences isotopiques entre les couches ($|\Delta C| \sim 0,01$).
   • Les gradients le long des axes mineurs ($x$ ou $y$) entraînent des différences nettement plus faibles. Cela est attribué à la géométrie du filament, qui restreint le déplacement des éléments fluides le long des axes mineurs, limitant ainsi le mélange de matériaux ayant des rapports isotopiques différents.
3. Mécanisme d'inhomogénéité : Le principal moteur de la différence isotopique résiduelle entre les couches est la différence initiale de centre de masse des couches, causée par le champ de vitesse turbulent. Les estimations semi-analytiques, qui tiennent compte du déplacement initial des éléments fluides, concordent avec les résultats de la simulation à un facteur 2, confirmant que la configuration turbulente initiale dicte le profil isotopique final.
4. Influence du champ magnétique : La présence d'un champ magnétique ($10 \mu\text{G}$) n'altère pas de manière significative le profil interne du rapport isotopique par rapport au cas purement hydrodynamique.
5. Implications pour le disque circumstellaire : En estimant le rayon centrifuge du matériel accrétant, les auteurs ont trouvé que l'hétérogénéité isotopique persiste dans le disque circumstellaire. Bien que le disque hérite de l'hétérogénéité du cœur, la distribution anisotrope du moment angulaire signifie que le gaz accrète à partir de diverses couches sphériques. Les résultats suggèrent qu'environ 1 % de l'hétérogénéité initiale du filament subsiste dans le disque.

Signification et affirmations
L'article affirme fournir un lien quantitatif entre les anomalies isotopiques observées dans les météorites et les conditions astrophysiques du nuage moléculaire parent. En combinant leurs résultats avec le facteur de réduction rapporté par Jacquet et al. (2019) (qui suggère une réduction d'un facteur $\sim 3$ du cœur au disque), les auteurs infèrent que l'hétérogénéité résiduelle dans le disque circumstellaire final est d'environ 0,03 % à 0,3 % de la variation initiale du filament.

Les auteurs affirment que ce niveau d'hétérogénéité résiduelle est suffisant pour expliquer les anomalies isotopiques observées dans les matériaux du système solaire (par exemple, les variations de $^{54}\text{Cr}/^{52}\text{Cr}$ et $^{50}\text{Ti}/^{47}\text{Ti}$ de 1 à 10 parties pour 10 000) sans nécessiter de processus thermiques extrêmes ou de scénarios d'injection nucléosynthétique exotiques post-effondrement. L'étude conclut que l'inhomogénéité isotopique des blocs de construction du système solaire pourrait naturellement survivre à l'ensemble du processus, de la fragmentation du filament à la formation du disque, à condition que le nuage parent possède un niveau approprié d'inhomogénéité initiale. Les auteurs notent que des simulations à haute résolution sont nécessaires pour les travaux futurs afin de résoudre pleinement les propriétés réalistes des rapports isotopiques dans le disque circumstellaire, particulièrement concernant l'infall et les flux sortants anisotropes.