Formation of Isotopically Heterogeneous Molecular Cloud Cores in Filamentary Molecular Clouds

Technische Zusammenfassung: Bildung isotopisch heterogener Molekülwolkenkerne in filamentären Molekülwolken

Problemstellung
Die Analyse von Meteoriten zeigt, dass Feststoffe des frühen Sonnensystems, wie etwa calcium-aluminium-reiche Einschlüsse (CAIs), isotopische Inhomogenitäten aufweisen (z. B. Variationen in den $^{54}\text{Cr}/^{52}\text{Cr}$- und $^{50}\text{Ti}/^{47}\text{Ti}$-Verhältnissen). Während thermische Prozesse in protoplanetaren Scheiben und die Vererbung aus inhomogenen Muttermolekülwolkenkernen als Erklärungen vorgeschlagen werden, bleibt das Ausmaß, in dem isotopische Heterogenität während des Kollapses filamentärer Molekülwolken in dichte Kerne überlebt, unklar. Insbesondere ist es notwendig zu quantifizieren, wie viel der anfänglichen isotopischen Variation, die in einer filamentären Wolke vorhanden ist, während der gravitativen Fragmentierung und Kernbildung bewahrt oder geglättet wird, und wie sich diese Heterogenität auf die resultierende stellare Scheibe überträgt.

Methodik
Die Autoren verwendeten Smoothed Particle Hydrodynamics (SPH)-Simulationen, speziell die Godunov-geglättete magnetohydrodynamische (GSPM)-Methode, um die gravitative Fragmentierung filamentärer Molekülwolken zu modellieren.

• Simulationsaufbau: Die Simulationen modellierten Filamente mit einer kritischen Linienmasse ($M_{\text{line,crit}} \approx 18 M_\odot \text{pc}^{-1}$) und einer Breite von $\sim 0,1$ pc, was mit Beobachtungen übereinstimmt. Die Anfangsbedingungen beinhalteten ein hydrostatisches Gleichgewicht-Dichteprofil, Kolmogorov-Turbulenz ($\sigma = 2c_s$) und ein initiales Magnetfeld von $10 \mu\text{G}$, das senkrecht zur Filamentachse steht. Es wurden auch rein hydrodynamische Fälle für den Vergleich durchgeführt.
• Isotopenverfolgung: Isotopische Spezies wurden als passive Tracer behandelt. Die Autoren konstruierten sechs verschiedene Modelle (Tabelle 1), um initiale isotopische Verhältnisgradienten ($C_0$) entlang verschiedener Achsen aufzuerlegen: die longitudinale Achse ($z$), die Nebenachsen ($x, y$) sowie Kombinationen daraus. Das Isotopenverhältnis wurde über das Rechengebiet normalisiert zwischen 0 und 1.
• Analyse: Die Entwicklung des Isotopenverhältnisses wurde durch das Verfolgen von SPH-Partikeln nachverfolgt. Die Autoren definierten sphärische Schalen um das Dichtemaximum der sich bildenden Kerne, um das durchschnittliche Isotopenverhältnis ($\bar{C}$) und die Differenz zum Kernzentrum ($|\Delta C|$) zu berechnen. Zudem verwendeten sie einen semi-analytischen Ansatz, um die initiale Verschiebung des Massenschwerpunkts der Schalen ($\delta z_{\text{ini}}$) basierend auf turbulenten Geschwindigkeitsfeldern und Dichtespektren abzuschätzen.

Kernergebnisse

1. Überleben der Heterogenität: Die Simulationen zeigen, dass isotopische Inhomogenitäten während der Filamentfragmentierung nicht vollständig geglättet werden. Kerne bewahren etwa 0,1 % bis 1 % der anfänglichen isotopischen Inhomogenität, die über die initiale Filamentlänge von 1 pc vorhanden war.
2. Geometrische Abhängigkeit: Das Ausmaß der bewahrten Heterogenität hängt von der Orientierung des initialen Isotopengradienten relativ zur Filamentgeometrie ab.
   • Gradienten entlang der longitudinalen Achse ($z$-Richtung, Modell $M_z$) führen zu den größten Isotopenunterschieden zwischen den Schalen ($|\Delta C| \sim 0,01$).
   • Gradienten entlang der Nebenachsen ($x$ oder $y$) führen zu signifikant kleineren Unterschieden. Dies wird auf die Geometrie des Filaments zurückgeführt, welches die Verschiebung von Fluidlementen entlang der Nebenachsen einschränkt und somit die Vermischung von Material mit unterschiedlichen Isotopenverhältnissen begrenzt.
3. Mechanismus der Inhomogenität: Der primäre Treiber der verbleibenden Isotopenunterschiede zwischen den Schalen ist der initiale Unterschied im Massenschwerpunkt der Schalen, verursacht durch das turbulente Geschwindigkeitsfeld. Die semi-analytischen Schätzungen, welche die initiale Verschiebung der Fluid-Elemente berücksichtigen, stimmen innerhalb eines Faktors von 2 mit den Simulationsergebnissen überein und bestätigen, dass die initiale turbulente Konfiguration das finale Isotopenprofil bestimmt.
4. Einfluss des Magnetfeldes: Die Anwesenheit eines Magnetfeldes ($10 \mu\text{G}$) verändert das interne Profil des Isotopenverhältnisses im Vergleich zum rein hydrodynamischen Fall nicht signifikant.
5. Implikationen für die stellare Scheibe: Durch die Abschätzung des Zentrifugerradius des akkretierenden Materials fanden die Autoren heraus, dass die isotopische Heterogenität in der stellaren Scheibe fortbesteht. Während die Scheibe die Heterogenität des Kerns erbt, bedeutet die anisotrope Verteilung des Drehimpulses, dass Gas aus verschiedenen sphärischen Schalen akkretiert. Die Ergebnisse legen nahe, dass etwa 1 % der initialen filamentären Heterogenität in der Scheibe verbleibt.

Bedeutung und Behauptungen
Das Paper beansprucht, eine quantitative Verbindung zwischen den in Meteoriten beobachteten Isotopenanomalien und den astrophysikalischen Bedingungen der Muttermolekülwolke herzustellen. Durch die Kombination ihrer Ergebnisse mit dem von Jacquet et al. (2019) berichteten Reduktionsfaktor (welcher eine Reduktion um den Faktor $\sim 3$ vom Kern zur Scheibe nahelegt), schließen die Autoren daraus, dass die verbleibende Heterogenität in der finalen stellaren Scheibe etwa 0,03 % bis 0,3 % der initialen filamentären Variation beträgt.

Die Autoren behaupten, dass dieses Niveau der verbleibenden Heterogenität ausreicht, um die beobachteten Isotopenanomalien in Materialien des Sonnensystems (z. B. $^{54}\text{Cr}/^{52}\text{Cr}$- und $^{50}\text{Ti}/^{47}\text{Ti}$-Variationen von 1–10 Teilen auf 10.000) zu erklären, ohne dass extreme thermische Prozesse oder exotische nukleosynthetische Injektionsszenarien nach dem Kollaps erforderlich sind. Die Studie kommt zu dem Schluss, dass die isotopische Inhomogenität der Bausteine des Sonnensystems den gesamten Prozess von der Filamentfragmentierung bis zur Scheibenbildung auf natürliche Weise überleben könnte, sofern die Mutterwolke eine angemessene initiale Inhomogenität besaß. Die Autoren merken an, dass hochauflösende Simulationen für zukünftige Arbeiten erforderlich sind, um die realistischen Eigenschaften der Isotopenverhältnisse in der stellaren Scheibe vollständig aufzulösen, insbesondere im Hinblick auf anisotropen Einfall und Ausflüsse.