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

Technische Samenvatting: Modellering van stofcoagulatie, dynamische weerstand en turbulente menging tijdens ster- en schijfvorming

Probleemstelling
Stof speelt een cruciale rol bij zowel de vorming van sterren als planeten, waarbij het de koeling van gas, ionisatiegraden en de koppeling van magnetische velden aan het gas beïnvloedt. Hoewel de coagulatie van stof tot grotere korrels (kiezelstenen) essentieel is voor de vorming van planetesimalen en de groei van kernen, hebben eerdere driedimensionale hydrodynamische simulaties van stervorming de dynamische evolutie van stof ten opzichte van het gas grotendeels genegeerd. Eerdere methoden, zoals die van Bate (2022), behandelden stofpopulaties als onafhankelijk evoluerend binnen elk fluïdumelement (SPH-deeltje), uitgaande van een perfecte koppeling tussen stof en gas. Deze benadering hield stand tijdens de vroege, laag-densiteitsfasen van de instorting, maar schoot tekort bij het verklaren van stofmigratie, sedimentatie en menging over langere tijdschalen, met name binnen protoplanetaire schijven waar de dichtheden hoger zijn en korrelgroei snel verloopt. Bovendien ontbraken in bestaande modellen vaak een zelfconsistente behandeling van hoe turbulentie zowel de relatieve snelheden nodig voor korrelgroei als de turbulente diffusie die stof verticaal mengt, aanstuurt.

Methodologie
De auteurs presenteren een nieuwe implementatie binnen de smoothed particle hydrodynamics (SPH) code SPHNG, die drie verschillende fysische processen combineert:

1. Stofcoagulatie: Door gebruik te maken van de gediscretiseerde mean-field Smoluchowski-vergelijking methode (Hirashita & Yan 2009; Bate 2022), volgt de code de evolutie van een populatie stofkorrels over $N$ logaritmisch gespatieerde grootteklassen voor elk SPH-gasdeeltje. De methode houdt rekening met binaire botsingen met een plakkingskans van één onder een specifieke coagulatie-drempelsnelheid.
2. Dynamische Weerstand (MULTIGRAIN): De auteurs integreren het MULTIGRAIN-algoritme (Hutchison et al. 2018), dat de continuüm-fluïdumvergelijkingen oplost voor een mengsel van gas en meerdere stofsoorten in de eindsnelheidsbenadering (terminal velocity approximation). Dit maakt de berekening mogelijk van bulk differentiële snelheden tussen verschillende stofgroottes en het gas als gevolg van drukgradiënten en zwaartekracht.
3. Impliciete Integratie en Turbulente Diffusie: Een belangrijke innovatie is het gebruik van een recent ontwikkelde impliciete integratiemethode (Elsender & Bate 2024) om de stofweerstandvergelijkingen op te lossen. Dit voorkomt de restrictieve tijdstap-beperkingen van expliciete methoden voor korrels met lange stop-tijden en voorkomt onfysische negatieve stoffracties. Daarnaast introduceren de auteurs een nieuwe impliciete methode om turbulente diffusie te modelleren. Zij behandelen de turbulente menging van stof als een diffusieproces ($D_i = \alpha_{SS} c_s H / Sc$), waarbij het Schmidt-getal ($Sc$) de koppelingsefficiëntie van verschillende korrelgroottes bepaalt. Deze diffusieterm wordt opgelost met een Gauss-Seidel operator-splitting benadering, waardoor stof tussen fluïdumelementen kan mengen.

De code volgt de evolutie van stof vanaf de instorting van een moleculaire wolk-kern, door de eerste hydrostatische kernfase, tot de vorming en evolutie van een protoplanetaire schijf, waarbij sink-deeltjes worden gebruikt om het systeem te volgen voorbij de vorming van de stellaire kern.

Belangrijkste Bijdragen

• Eerste 3D Gekoppelde Simulatie: Dit werk vertegenwoordigt de eerste driedimensionale hydrodynamische simulatie die gelijktijdig stofcoagulatie, dynamische gasweerstand en turbulente menging tijdens de vroege stadia van ster- en schijfvorming modelleert.
• Impliciet Diffusieschema: De ontwikkeling van een robuuste impliciete solver voor turbulente stofdiffusie maakt stabiele simulaties van stofmenging mogelijk zonder de zware tijdstap-beperkingen die geassocieerd worden met expliciete diffusieschema's.
• Verenigde Behandeling: De code verenigt de behandeling van relatieve snelheden (die groei aandrijft) en turbulente diffusie (die menging aandrijft) binnen één enkel kader, waarbij dezelfde turbulentieparameter ($\alpha_{SS}$) wordt gebruikt voor beide processen in de testgevallen.

Resultaten
De auteurs presenteren resultaten van zowel testberekeningen (stofsedimentatie in een verticale kolom) als volledige stervormingssimulaties:

• Testberekeningen: De code herstelt succesvol analytische oplossingen voor het evenwicht tussen stofsedimentatie en turbulente diffusie. Tests bevestigen dat de inclusie van turbulente diffusie de vorming van onrealistisch dunne stoflagen voorkomt en een eindige stof-schaalhoogte handhaaft.
• Stervormingssimulaties:
  • Eerste Hydrostatische Kern: Consistent met eerder werk (Bate 2022), vindt er snelle stofgroei plaats binnen de eerste hydrostatische kern, wat lokaal mono-dispers verspreide populaties produceert met groottes die de 100 $\mu$m overschrijden diep in de kern. De inclusie van weerstand en diffusie had een verwaarloosbare impact op deze vroege fase vanwege de korte evolutionaire tijdschalen.
  • Protoplanetaire Schijf Evolutie: Zodra de schijf gevormd wordt, verandert de inclusie van turbulente diffusie de stofevolutie aanzienlijk.
    • Versterkte Groei: Turbulente diffusie versterkt de snelheid van stofkorrelgroei aanzienlijk. Door continu kleine en middelgrote korrels aan te voeren naar gebieden met een hoge dichtheid (waar grote korrels groeien), voorkomt diffusie de lokale uitputting van "voedsel" voor coagulatie.
    • Afhankelijkheid van Turbulentie: Het niveau van turbulentie ($\alpha_{SS}$) is cruciaal. Hoge turbulentie ($\alpha_{SS} = 10^{-3}$) leidt tot snelle groei naar mono-dispers verspreide populaties van $\sim30-100$ $\mu$m korrels door de hele schijf. Lagere turbulentie ($\alpha_{SS} = 10^{-4}$) resulteert in tragere groei en een meer verticaal gestratificeerde stofpopulatie, waarbij grotere korrels dichter bij het mid-vlak geconcentreerd blijven.
    • Stofsedimentatie: In simulaties met zwakke turbulentie leidt stofsedimentatie tot een significante toename van de stof-gasverhouding in het mid-vlak (met >50%), terwijl kleine korrels uitgeput blijven nabij het mid-vlak en beperkt blijven tot de oppervlaktelagen.
  • Weerstand versus Diffusie: Vergelijkingen van simulaties met en zonder specifieke processen laten zien dat de dynamische stof-gasweerstand weinig effect heeft op de stofdistributies op deze tijdschalen, aangezien korrels goed gekoppeld blijven. In tegen tegenstelling hiertoe heeft de inclusie van turbulente diffusie een enorm effect, waarbij het stofdistributies gladstrijkt en coagulatie versnelt.

Betekenis en Claims
Het artikel beweert dat deze nieuwe methodologie voor het eerst een zelfconsistente modellering van stofgroei, menging en migratie tijdens de vorming van sterren en hun schijven mogelijk maakt. De auteurs benadrukken dat hoewel de inclusie van weerstand en diffusie de fundamentele bevindingen over stofgroei in de eerste hydrostatische kern niet verandert, het cruciaal is voor de daaropvolgende evolutie van de protoplanetaire schijf. Specifiek claimen zij dat turbulente stofdiffusie zorgt voor een substantiële versterking van de snelheid van stofkorrelgroei door te fungeren als een bron van kleine korrels voor de regio's waar grote korrels worden gevormd.

De auteurs blijven bescheiden over de reikwijdte van hun bevindingen, waarbij zij opmerken dat:

• De resultaten afhangen van het gehanteerde turbulentiemodel (bijv. $\alpha$-schijf versus MHD-turbulentie).
• De simulaties momenteel coagulatie overwegen maar geen fragmentatie; verschillende coagulatiedrempels zouden tot fragmentatie en herpopulatie van kleine korrels kunnen leiden.
• De radiatieve transferberekeningen nog steeds standaard interstellaire opaciteiten aannemen in plaats van de evoluerende lokale korrelgrootteverdelingen, een stap die geïdentificeerd is als een noodzakelijke toekomstige uitbreiding.

Uiteindelijk biedt dit werk een instrument om 3D hydrodynamische simulaties te vergelijken met de rijkdom aan nieuwe observationele data van faciliteiten zoals ALMA, VLT-SPHERE en JWST, met name met betrekking tot stofsedimentatie en turbulente menging in edge-on schijven.