Modelling dust coagulation, dynamical drag and turbulent mixing during star and disc formation

Technische Zusammenfassung: Modellierung der Staubkoagulation, des dynamischen Drags und der turbulenten Vermischung während der Stern- und Scheibenbildung

Problemstellung
Staub spielt eine entscheidende Rolle sowohl bei der Stern- als auch bei der Planetenentstehung, da er die Kühlung des Gases, die Ionisationsgrade und die Kopplung magnetischer Felder an das Gas beeinflusst. Während die Koagulation von Staub zu größeren Körnern (Pebbles/Kieseln) essenziell für die Planetesimalbildung und das Kernwachstum ist, haben frühere dreidimensionale hydrodynamische Simulationen der Sternentstehung die dynamische Entwicklung des Staubs relativ zum Gas weitgehend vernachlässigt. Frühere Methoden, wie etwa die von Bate (2022), behandelten Staubpopulationen als unabhängig innerhalb jedes Fluidelements (SPH-Partikel) und nahmen an, dass Staub perfekt an das Gas gekoppelt sei. Diese Annahme galt für die frühen, niederdichten Phasen des Kollapses, versäumte es jedoch, die Migration, das Absetzen (Settling) und die Vermischung von Staub über längere Zeitskalen hinweg zu berücksichtigen – insbesondere innerhalb protoplanetarer Scheiben, in denen die Dichten höher und das Kornwachstum schneller ist. Darüber hinaus fehlte es bestehenden Modellen oft an einer konsistenten Behandlung der Frage, wie Turbulenz sowohl die relativen Geschwindigkeiten, die für das Kornwachstum notwendig sind, als auch die turbulente Diffusion, die den Staub vertikal vermischt, antreibt.

Methodik
Die Autoren präsentieren eine neue Implementierung innerhalb des Smoothed Particle Hydrodynamics (SPH)-Codes SPHNG, die drei unterschiedliche physikalische Prozesse kombiniert:

1. Staubkoagulation: Unter Verwendung der diskretisierten Mean-Field-Smoluchowski-Gleichung (Hirashita & Yan 2009; Bate 2022) verfolgt der Code die Entwicklung einer Staubpopulation über $N$ logarithmisch gestaffelte Größenklassen für jedes SPH-Gaspartikel. Die Methode berücksichtigt binäre Kollisionen mit einer Haftwahrscheinlichkeit von eins unterhalb einer spezifischen Koagulations-Schwellengeschwindigkeit.
2. Dynamischer Drag (MULTIGRAIN): Die Autoren integrieren den MULTIGRAIN-Algorithmus (Hutchison et al. 2018), der die Kontinuums-Flüssigkeitsgleichungen für ein Gemisch aus Gas und mehreren Staubarten unter der Annahme der Endgeschwindigkeit (Terminal Velocity Approximation) löst. Dies ermöglicht die Berechnung der mittleren Differenzgeschwindigkeiten zwischen verschiedenen Staubgrößen und dem Gas aufgrund von Druckgradienten und Gravitation.
3. Implizite Integration und turbulente Diffusion: Eine zentrale Innovation ist die Verwendung einer kürzlich entwickelten impliziten Integrationsmethode (Elsender & Bate 2024) zur Lösung der Staub-Drag-Gleichungen. Dies vermeidet die restriktiven Zeitschrittbeschränkungen expliziter Methoden für Körner mit langen Haltezeiten und verhindert unphysikalische negative Staubfraktionen. Zusätzlich führen die Autoren eine neue implizite Methode zur Modellierung der turbulenten Diffusion ein. Sie behandeln die turbulente Durchmischung des Staubs als Diffusionsprozess ($D_i = \alpha_{SS} c_s H / Sc$), wobei die Schmidt-Zahl ($Sc$) die Kopplungseffizienz verschiedener Korngrößen berücksichtigt. Dieser Diffusionsterm wird mittels eines Gauss-Seidel-Operator-Splitting-Ansatzes gelöst, was es dem Staub ermöglicht, sich zwischen Fluidelementen zu vermischen.

Der Code verfolgt die Entwicklung des Staubs vom Kollaps eines Molekülwolkenkerns über die Phase des ersten hydrostatischen Kerns bis hin zur Bildung und Entwicklung einer protoplanetaren Scheibe, wobei Sink-Partikel verwendet werden, um das System über die Sternkernbildung hinaus zu verfolgen.

Wesentliche Beiträge

• Erste gekoppelte 3D-Simulation: Diese Arbeit stellt die erste dreidimensionale hydrodynamische Simulation dar, die Staubkoagulation, dynamischen Gas-Drag und turbulente Vermischung während der frühen Stadien der Stern- und Scheibenbildung gleichzeitig modelliert.
• Implizites Diffusionsschema: Die Entwicklung eines robusten impliziten Lösers für die turbulente Staubdiffusion ermöglicht stabile Simulationen der Staubvermischung ohne die schweren Zeitschrittbeschränkungen, die mit expliziten Diffusionsschemata verbunden sind.
• Vereinheitlichte Behandlung: Der Code vereint die Behandlung von relativen Geschwindigkeiten (die das Wachstum antreiben) und turbulenter Diffusion (die die Vermischung antreibt) in einem einzigen Rahmenwerk, indem er in den Testfällen denselben Turbulenzparameter ($\alpha_{SS}$) für beide Prozesse verwendet.

Ergebnisse
Die Autoren präsentieren Ergebnisse sowohl aus Testrechnungen (Staubabsetzung in einer vertikalen Säule) als auch aus vollständigen Sternentstehungssimulationen:

• Testrechnungen: Der Code reproduziert erfolgreich analytische Lösungen für das Gleichgewicht zwischen Staubabsetzung und turbulenter Diffusion. Tests bestätigen, dass die Einbeziehung der turbulenten Diffusion die Bildung unrealistisch dünner Staubschichten verhindert und eine endliche Staub-Skalenhöhe aufrechterhält.
• Sternentstehungssimulationen:
  • Erster hydrostatischer Kern: Konsistent mit vorangegangenen Arbeiten (Bate 2022) findet innerhalb des ersten hydrostatischen Kerns ein schnelles Staubwachstum statt, das lokal monodisperse Populationen mit Größen von über 100 $\mu$m tief im Inneren des Kerns erzeugt. Die Einbeziehung von Drag und Diffusion hatte aufgrund der kurzen Evolutionszeitlinien einen vernachlässigbaren Einfluss auf diese frühe Phase.
  • Entwicklung der protoplanetaren Scheibe: Sobald sich die Scheibe bildet, verändert die Einbeziehung der turbulenten Diffusion die Staubentwicklung signifikant.
    • Verstärktes Wachstum: Die turbulente Diffusion verstärkt die Rate des Staubkornwachstums erheblich. Indem sie kontinuierlich kleine und mittelgroße Körner in Regionen hoher Dichte (in denen große Körner wachsen) liefert, verhindert die Diffusion die lokale Erschöpfung des „Nahrungsguts“ für die Koagulation.
    • Abhängigkeit von der Turbulenz: Das Niveau der Turbulenz ($\alpha_{SS}$) ist entscheidend. Hohe Turbulenz ($\alpha_{SS} = 10^{-3}$) führt zu einem schnellen Wachstum zu monodispersen Populationen von $\sim30-100$ $\mu$m großen Körnern im gesamten Disk. Niedrigere Turbulenz ($\alpha_{SS} = 10^{-4}$) führt zu langsamerem Wachstum und einer stärker vertikal geschichteten Staubpopulation, wobei größere Körner näher an der Mittelebene konzentriert sind.
    • Staubabsetzung (Settling): In Simulationen mit schwacher Turbulenz führt das Absetzen des Staubs zu einer signifikanten Erhöhung des Staub-Gas-Verhältnisses in der Mittelebene (um >50 %), während kleine Körner in der Nähe der Mittelebene verarmt sind und in den Oberflächenschichten verbleiben.
  • Drag vs. Diffusion: Vergleiche von Simulationen mit und ohne spezifische Prozesse zeigen, dass der dynamische Staub-Gas-Drag auf diesen Zeitskalen kaum Auswirkungen auf die Staubverteilungen hat, da die Körner gut gekoppelt bleiben. Im Gegensatz dazu hat die Einbeziehung der turbulenten Diffusion einen enormen Effekt, indem sie die Staubverteilungen glättet und die Koagulation beschleunigt.

Bedeutung und Ansprüche
Das Paper behauptet, dass diese neue Methodik erstmals eine selbstkonsistente Modellierung von Staubwachstum, Vermischung und Migration während der Entstehung von Sternen und deren Scheiben ermöglicht. Die Autoren betonen, dass die Einbeziehung von Drag und Diffusion zwar die grundlegenden Erkenntnisse bezüglich des Staubwachstums im ersten hydrostatischen Kern nicht verändert, aber entscheidend für die anschließende Entwicklung der protoplanetaren Scheibe ist. Insbesondere behaupten sie, dass die turbulente Staubdiffusion eine substanzielle Steigerung der Rate des Staubkornwachstums bewirkt, indem sie als Quelle für kleine Körner in Regionen fungiert, in denen große Körner entstehen.

Die Autoren bleiben hinsichtlich des Umfangs ihrer Ergebnisse bescheiden und merken an:

• Die Ergebnisse hängen vom gewählten Turbulenzmodell ab (z. B. $\alpha$-Disk vs. MHD-Turbulenz).
• Die Simulationen berücksichtigen derzeit die Koagulation, aber nicht die Fragmentierung; unterschiedliche Koagulationsschwellen könnten zu Fragmentierung und der Neubildung kleiner Körner führen.
• Die Strahlungstransferberechnungen verwenden weiterhin Standard-Interstellar-Opazitäten anstatt der sich entwickelnden lokalen Korngrößenverteilungen – ein Schritt, der als notwendige zukünftige Erweiterung identifiziert wurde.

Letztendlich bietet die Arbeit ein Werkzeug, um 3D-hydrodynamische Simulationen mit der Fülle neuer Beobachtungsdaten von Einrichtungen wie ALMA, VLT-SPHERE und JWST zu vergleichen, insbesondere im Hinblick auf Staubabsetzung und turbulente Durchmischung in Edge-on-Scheiben.