Infall onto the Protoplanetary Disk during the Gravitational Collapse of a Molecular Cloud

Titel: Inval op de Protoplanetaire Schijf tijdens de Gravitatie-instorting van een Moleculaire Wolk

Auteurs: E.R. Redkin en E.I. Vorobyov
Datum: 24 februari 2026 (gepubliceerd op arXiv)

---

1. Probleemstelling

De vorming van protoplanetaire schijven en jonge sterrenobjecten (YSO's) wordt vaak gemodelleerd via de gravitatie-instorting van koude moleculaire wolken. Een cruciale stap in deze numerieke studies is het modelleren van de instorting van materie vanuit de wolk naar de schijf (de zogenaamde "infall").

De meest gebruikte aanpak hiervoor is het model van Nakamoto en Nakagawa (1994) (hierna NN-model). Hoewel dit model populair is vanwege zijn eenvoudige numerieke implementatie in één-dimensionale viskeuze evolutiemodellen, vertoont het ernstige tekortkomingen:

1. Idealistische dichtheidsverdeling: Het gaat uit van een "singular isothermische sfeer" (SIS) met een wiskundige singulariteit in het centrum, wat fysisch onrealistisch is voor de centrale regio's van een wolk.
2. Gebrek aan wolksgrens: Het model beschouwt de wolk als oneindig groot, wat het onmogelijk maakt om parameters zoals de verhouding tussen rotatie- en gravitatie-energie ($\beta$) realistisch te berekenen.
3. Constante instortingsrate: Het model veronderstelt een constante massa-instortingsrate ($\dot{M}_c$) totdat een vooraf bepaald totaal massabudget is bereikt, waarna deze abrupt stopt. Dit negeert de fysica van een eindige wolk, waarbij een verdunningsgolf (rarefaction wave) de instortingsrate in de loop van de tijd doet afnemen.

Deze beperkingen leiden tot discrepanties tussen simulaties en observationele data, met name wat betreft de duur van de "ingebouwde fase" (embedded phase) van YSO's (Class 0 en Class I).

2. Methodologie

De auteurs analyseren het NN-model en stellen een reeks technische modificaties voor om de fysieke realisme te vergroten zonder de numerieke eenvoud te verliezen. De kern van de methode omvat:

• Definitie van een wolksgrens: In plaats van een oneindige wolk, wordt een buitenste straal ($r_{out}$) geïntroduceerd waarbinnen de totale massa gelijk is aan de wolkmassa ($M_c$). Dit maakt de berekening van de verhouding $\beta = E_{rot} / |E_{grav}|$ mogelijk.
• Verbeterde dichtheidsverdeling: De SIS-verdeling ($\rho \propto r^{-2}$) wordt vervangen door een profiel dat een centraal plateau heeft, beschreven door een Bonnor-Ebert-sfeer-achtige verdeling:
  $$\rho(r) = \frac{\rho_0}{1 + (r/r_0)^2}$$
  Dit profiel is fysisch realistischer voor de pre-stellare condensatie en convergeert asymptotisch naar de SIS op grote afstanden.
• Tijdsafhankelijke instortingsrate ($\dot{M}_c$): De constante instortingsrate wordt vervangen door een functie die rekening houdt met de eindige grootte van de wolk. De rate blijft aanvankelijk constant totdat een verdunningsgolf vanuit de buitenste rand de instortingsfront bereikt (tijd $t_d$), waarna de rate exponentieel afneemt. De auteurs stellen een nieuwe formule (Modificatie M2) voor die garandeert dat de totale geaccreteerde massa exact gelijk is aan de initiële wolkmassa ($M_c$), ongeacht de parameters.
• Aanpassing van de schijf-evolutievergelijking: Er wordt erkend dat materie met een sub-Kepleriaans snelheid op de schijf valt. Dit vereist een extra term in de continuïteitsvergelijking voor de oppervlaktedichtheid ($\Sigma$) om het verschil in hoekmoment tussen de invallende materie en de schijfmaterie correct te modelleren.
• Disc-gemedieerde accretie: Voor het geval de instorting alleen op de schijf gebeurt (en niet direct op de ster), wordt een correctiefactor ($\xi$) ingevoerd in de instortingsformule om de binnenste rand van de schijf ($R_{in}$) te respecteren.

3. Belangrijkste Bijdragen

De paper introduceert vier hoofdverbeteringen aan het bestaande NN-model:

1. Berekening van de buitenste grens: Hierdoor kunnen sleutelparameters van moleculaire wolken (zoals $\beta$) direct worden gekoppeld aan de eigenschappen van de resulterende protoplanetaire schijf.
2. Realistische dichtheidsprofielen: De vervanging van de SIS door een profiel met een centraal plateau (Bonnor-Ebert) verbetert de fysieke nauwkeurigheid in de centrale regio's.
3. Analytische uitdrukking voor afnemende instortingsrate: Een nieuwe formule voor $\dot{M}_c(t)$ die de fysica van een eindige wolk weergeeft en zorgt voor een massabalans die consistent is met de initiële condities.
4. Veralgemening van de instortingsvergelijking: Een meer algemene vorm van de vergelijking voor $\dot{\Sigma}_{infall}$ die de binnenste rand van de schijf expliciet meeneemt, essentieel voor het modelleren van disk-gemedieerde accretie na de vorming van de protoster.

4. Resultaten

De auteurs voerden simulaties uit met de gemodificeerde modellen (NN-1 tot NN-7) en vergeleken deze met observationele data:

• Invloed van parameters: De duur van de instorting en de evolutie van de YSO hangen sterk af van de temperatuur van de wolk ($T_c$) en de amplitude van de dichtheidsverstoring ($A$).
• Duur van de ingebouwde fase: De gemodificeerde modellen voorspellen over het algemeen een snellere overgang van Class 0 naar Class I en van Class I naar Class II dan waargenomen in de praktijk.
  • Bij hoge waarden van $A$ (sterke dichtheidsverstoring) is de instortingsrate zeer hoog, wat leidt tot een zeer korte duur van de Class 0-fase (soms < 100 kyr).
  • Alleen bij lage temperaturen ($\approx 7$ K) of zeer lage verstoringen ($A \approx 2.2$) komen de simulatieresultaten dichter bij de observationele schattingen (Class 0: 0.15–0.24 Myr; Class I: 0.31–0.48 Myr).
• Oorzaak van discrepantie: De auteurs stellen dat de snellere evolutie in de modellen mogelijk te wijten is aan:
  1. Observatiebias (bijv. edge-on systemen worden verkeerd geclassificeerd als Class 0).
  2. Beperkingen van het model zelf, dat een geïsoleerde wolk veronderstelt zonder massatoevoer vanuit omringende filamenten (wat in de realiteit de instortingsrate kan verlengen).

5. Betekenis en Conclusie

Deze studie is van groot belang voor de numerieke astrofysica omdat het een veelgebruikt, maar onvolmaakt model (Nakamoto & Nakagawa) verfijnt. De voorgestelde modificaties verhogen de fysieke betrouwbaarheid aanzienlijk zonder de rekencomplexiteit te verhogen, waardoor het model geschikt blijft voor gebruik in één-dimensionale viskeuze evolutiesimulaties.

De resultaten benadrukken dat de initiële condities van de moleculaire wolk (temperatuur en dichtheidsverstoring) een doorslaggevende rol spelen in de tijdschaal van de ster- en schijfvorming. De discrepantie met observationele data suggereert dat toekomstige modellen rekening moeten houden met externe massatoevoer (via filamenten) en complexere classificatiemethoden voor YSO's.

Deze paper vormt het eerste deel van een serie studies; de volgende publicatie zal de schijfmodellen en de resultaten van de numerieke simulaties van de evolutie zelf presenteren.