Formation of Isotopically Heterogeneous Molecular Cloud Cores in Filamentary Molecular Clouds

Technische Samenvatting: Vorming van Isotopisch Heterogene Moleculaire Wolkkeren in Filamentaire Moleculaire Wolken

Probleemstelling
Meteorietanalyse onthult dat vaste stoffen uit het vroege zonnestelsel, zoals calcium-aluminium-rijke inclusies (CAI's), isotopische inhomogeniteit vertonen (bijv. variaties in $^{54}\text{Cr}/^{52}\text{Cr}$ en $^{50}\text{Ti}/^{47}\text{Ti}$-verhoudingen). Hoewel thermische verwerking in protoplanetaire schijven en overerving van inhomogene ouderlijke moleculaire wolkkeren als verklaringen worden voorgesteld, blijft de mate waarin isotopische heterogeniteit de ineenstorting van filamentaire moleculaire wolken naar dichte kernen overleeft, onduidelijk. Specifiek is het noodzakelijk om te kwantificeren hoeveel van de initiële isotopische variatie die aanwezig is in een filamentaire wolk behouden blijft of wordt uitgevlakt tijdens gravitationele fragmentatie en kernvorming, en hoe deze heterogeniteit zich vertaals naar de resulterende circumstellaire schijf.

Methodologie
De auteurs maakten gebruik van Smoothed Particle Hydrodynamics (SPH) simulaties, specifiek gebruikmakend van de Godunov smoothed particle magnetohydrodynamische (GSPM) methode, om de gravitationele fragmentatie van filamentaire moleculaire wolken te modelleren.

• Simulatie-opstelling: De simulaties modelleerden filamenten met een kritische lijnmassa ($M_{\text{line,crit}} \approx 18 M_\odot \text{pc}^{-1}$) en een breedte van $\sim 0,1$ pc, consistent met waarnemingen. De initiële condities omvatten een hydrostatisch evenwichtig dichtheidsprofiel, Kolmogorov-turbulentie ($\sigma = 2c_s$), en een initieel magnetisch veld van $10 \mu\text{G}$ loodrecht op de filamentas. Alleen-hydrodynamische gevallen werden ook uitgevoerd voor vergelijking.
• Isotopische Sporing: Isotopische soorten werden behandeld als passieve tracers. De auteurs construeerden zes afzonderlijke modellen (Tabel 1) om initiële isotopische ratio-gradiënten ($C_0$) langs verschillende assen op te leggen: de longitudinale as ($z$), de minore assen ($x, y$), en combinaties daarvan. De isotopische ratio werd genormaliseerd tussen 0 en 1 over het computationele domein.
• Analyse: De evolutie van de isotopische ratio werd gevolgd door SPH-deeltjes te volgen. De auteurs definieerden sferische schillen rond de dichtheidspiek van vormende kernen om de gemiddelde isotopische ratio ($\bar{C}$) en het verschil met de kerncentra ($|\Delta C|$) te berekenen. Ze gebruikten ook een semi-analytische benadering om de initiële zwaartepuntverplaatsing van schillen ($\delta z_{\text{ini}}$) te schatten op basis van turbulente snelheidvelden en dichtheidsspectra.

Belangrijkste Resultaten

1. Overleving van Heterogeniteit: De simulaties demonstreren dat isotopische inhomogeniteiten niet volledig worden uitgevlakt tijdens de fragmentatie van het filament. Kernen behouden ongeveer 0,1% tot 1% van de initiële isotopische inhomogeniteit die aanwezig was over de initiële filamentlengte van 1 pc.
2. Geometrische Afhankelijkheid: De omvang van de behouden heterogeniteit hangt af van de oriëntatie van de initiële isotopische gradiënt ten opzichte van de filamentgeometrie.
   • Gradiënten langs de longitudinale as ($z$-richting, Model $M_z$) resulteren in de grootste isotopische verschillen tussen schillen ($|\Delta C| \sim 0,01$).
   • Gradiënten langs de minore assen ($x$ of $y$) resulteren in aanzienlijk kleinere verschillen. Dit wordt toegeschreven aan de geometrie van het filament, die de verplaatsing van fluïdum-elementen langs de minore assen beperkt, waardoor het mengen van materiaal met verschillende isotopische ratio's wordt beperkt.
3. Mechanisme van Inhomogeniteit: De primaire drijfveer van het resterende isotopische verschil tussen schillen is het initiële verschil in het zwaartepunt van de schillen, veroorzaakt door het turbulente snelheidveld. De semi-analytische schattingen, die rekening houden met de initiële verplaatsing van fluïdum-elementen, komen binnen een factor 2 overeen met de simulatie-resultaten, wat bevestigt dat de initiële turbulente configuratie het uiteindelijke isotopische profiel bepaalt.
4. Invloed van het Magnetisch Veld: De aanwezigheid van een magnetisch veld ($10 \mu\text{G}$) verandert het interne profiel van de isotopische ratio niet significant vergeleken met het puur hydrodynamische geval.
5. Implicaties voor de Circumstellaire Schijf: Door de centrifugale straal van accreterend materiaal te schatten, vonden de auteurs dat isotopische heterogeniteit persisteert in de circumstellaire schijf. Hoewel de schijf de heterogeniteit van de kern overneemt, betekent de anisotrope distributie van impulsmoment dat gas vanuit verschillende sferische schillen accreteert. De resultaten suggereren dat ongeveer 1% van de initiële filamentaire heterogeniteit in de schijf overblijft.

Betekenis en Claims
Het artikel claimt een kwantitatieve link te bieden tussen de isotopische anomalieën waargenomen in meteorieten en de astrofysische condities van de ouderlijke moleculaire wolk. Door hun bevindingen te combineren met de reductiefactor gerapporteerd door Jacquet et al. (2019) (die een reductie van factor $\sim 3$ suggereert van kern naar schijf), leiden de auteurs af dat de resterende heterogeniteit in de uiteindelijke circumstellaire schijf ongeveer 0,03% tot 0,3% is van de initiële filamentaire variatie.

De auteurs stellen dat dit niveau van resterende heterogeniteit voldoende is om de waargenomen isotopische anomalieën te verklaren (bijv. $^{54}\text{Cr}/^{52}\text{Cr}$ en $^{50}\text{Ti}/^{47}\text{Ti}$ variaties van 1–10 delen per 10.000) zonder dat daar extreme thermische verwerking of exotische nucleosynthetische injectiescenario's na de ineenstorting voor nodig zijn. De studie concludeert dat de isotopische inhomogeniteit van de bouwstenen van het zonnestelsel op natuurlijke wijze het hele proces van filamentfragmentatie tot schijfvorming kan overleven, mits de ouderlijke wolk een gepaste mate van initiële inhomogeniteit bezat. De auteurs merken op dat simulaties met een hogere resolutie vereist zijn voor toekomstig werk om de realistische eigenschappen van isotopische ratio's in de circumstellaire schijf volledig op te lossen, met name met betrekking tot anisotrope instroom en uitstroom.