The dispersal of compact protoplanetary discs

Technische Samenvatting: De Verspreiding van Compacte Protoplanetaire Schijven

Probleemstelling
Protoplanetaire schijven zijn de locaties van planeetvorming, maar hun verspreidingsmechanismen blijven onderwerp van actief onderzoek. Observationeel bewijs, met name uit kleurencolor-evolutiestudies van jonge sterpopulaties, geeft aan dat schijfverspreiding bij voorkeur van binnen naar buiten verloopt, waarbij een centrale holte ontstaat voordat de buitenste schijf wordt geërodeerd. Echter, standaard theoretische modellen gaan vaak uit van schijven met grote buitenstralen (honderden astronomische eenheden), ver buiten de gravitationele straal ($R_g$) waar fotoverdamping het meest efficiënt is. Recente surveys (bijv. ALMA) tonen aan dat een aanzienlijk deel van de schijven "compact" is, met buitenstralen die vergelijkbaar zijn met of kleiner dan $R_g$. Bestaande modellen, die doorgaans uitgaan van onbeperkte massa-verliesprofielen die zich uitstrekken tot grote stralen, slagen er niet in om de evolutionaire banen en verspreidingstijdschalen van deze compacte systemen nauwkeurig weer te geven. Specifiek voorspellen standaard voorschriften vaak erosie van buiten naar binnen voor compacte schijven, wat in strijd is met observationele diagnostiek.

Methodologie
De auteurs hanteerden een meerstapsaanpak waarbij hoogresolutie stralings-hydrodynamische simulaties werden gecombineerd met 1D viskeuze populatiesynthesemodellen:

1. Stralings-hydrodynamische Simulaties: Met behulp van de PLUTO-code gekoppeld aan het MOCASSIN-stralingstransportframework simuleerden de auteurs door röntgenstraling aangedreven fotoverdamping voor een ster van $1 M_\odot$. Ze onderzochten vier verschillende initiële schijfafsnijstralen ($R_{cut}$): 10, 50, 100 en 200 au. De simulaties volgden de thermische structuur, dichtheid en snelheidsvelden totdat een quasi-stationaire toestand was bereikt. De gravitationele straal werd aangenomen als $R_g = 5$ au.
2. Afleiding van Schalingsrelaties: Uit de hydrodynamische simulaties analyseerden de auteurs de cumulatieve verdelingen van het massa-verlies. Ze bepaalden hoe het totale wind-massa-verlies afhangt van de buitenste afsnijstraal van de schijf.
3. Populatiesynthese: Deze bevindingen werden opgenomen in een 1D viskeus evolutiecode (SPOCK). De auteurs voerden een populatiesynthese uit van 4.096 synthetische schijven met behulp van een Quasi-Monte Carlo-steekproefstrategie (Sobol-sequenties) om een 9-dimensionale parameterruimte te verkennen (stermassa, röntgenluminositeit, initiële schijfmassa, afsnijstraal, viscositeit, extern FUV-veld, enz.).
4. Geteste Voorschriften: De modellen vergeleken twee scenario's:
   • Standaard Voorschrift: Interne fotoverdampings-massa-verliespercentages toegepast zonder rekening te houden met de eindige buitenstraal van de schijf (effectief aannemende dat de schijf zich uitstrekt tot oneindig of voorbij het piek-massa-verlies).
   • Afsnijstraal-afhankelijk Voorschrift: Het massa-verliespercentage wordt afgekapt en genormaliseerd op basis van de specifieke afsnijstraal van de schijf, afgeleid uit de hydrodynamische schalingsrelaties.
   • Externe Fotoverdamping: Sommige modellen omvatten milde externe FUV-fotoverdamping (met behulp van het FRIED-rooster) om nabijgelegen OB-sterren te simuleren.

Belangrijkste Bijdragen en Resultaten

• Hydrodynamische Schalings: De simulaties onthulden dat de oppervlakte-massa-verliesprofielen bijna onafhankelijk zijn van de schijfgrootte; de windmorfologie en lokale massa-verliespercentages in de binnenste gebieden blijven consistent ongeacht de buitenste afsnijstraal. Het geïntegreerde totale massa-verliespercentage is echter aanzienlijk verminderd voor compacte schijven omdat de wind wordt afgekapt bij de afsnijstraal. Het cumulatieve massa-verliespercentage schaalt met de straal en groeit als $\dot{M}_w \propto R^3$ nabij de gravitationele straal en vloeit af naar $\dot{M}_w \propto R^{0.5}$ bij grotere stralen.
• Evolutionaire Banen:
  • Standaard Voorschrift: Wanneer toegepast op compacte schijven ($R_{cut} < 20-30$ au), voorspelt het standaardmodel dat het wind-massa-verliespercentage de viskeuze spreidingsnelheid buiten de afsnijstraal overtreft. Dit leidt tot verspreiding van buiten naar binnen, waarbij de buitenste randen eerst worden geërodeerd, in strijd met observationele bewijslast.
  • Afsnijstraal-afhankelijk Voorschrift: Wanneer het massa-verliespercentage wordt geschaald naar de afsnijstraal, wordt het totale windpercentage voldoende verminderd om viskeuze spreiding aanvankelijk te laten domineren. Dit resulteert in verspreiding van binnen naar buiten, waarbij een spleet opent nabij de gravitationele straal en zich naar buiten uitbreidt, in overeenstemming met observationele diagnostiek.
• Resultaten van Populatiesynthese:
  • Het afsnijstraal-afhankelijke voorschrift, gecombineerd met een verminderd intern fotoverdampingspercentage (geschaald met een factor 10 om overeen te komen met recente chemische modellen), reproduceert waargenomen schijffracties en accretiepercentages als functie van de clusterleeftijd zonder parameteraanpassing.
  • De mediale schijflevensduur in het verbeterde model bedraagt $2,70 \pm 0,13$ Myr, wat goed overeenkomt met observationele schattingen ($2,24 \pm 0,35$ Myr) en aanzienlijk langer is dan voorspellingen van modellen zonder de reductie van het percentage ($1,13 \pm 0,05$ Myr).
  • De opname van milde externe fotoverdamping verhindert de schijfspreiding niet wanneer het afsnijstraal-voorschrift wordt gebruikt, maar resulteert in een zachtere toename van schijfstralen in de loop van de tijd en een hoger percentage compacte schijven dat tot latere leeftijden blijft bestaan.

Betekenis
Het artikel betoogt dat de verspreiding van compacte protoplanetaire schijven niet nauwkeurig kan worden vastgelegd door standaard, niet-geschaalde fotoverdampingsmodellen. De auteurs stellen een eenvoudige, fysiek gemotiveerde schalingsrelatie voor die rekening houdt met de eindige buitenstraal van de schijf. Door deze afhankelijkheid op te nemen, slagen de modellen erin theoretische voorspellingen te verzoenen met de waargenomen van binnen naar buiten gerichte ruiming van jonge sterpopulaties en de prevalentie van compacte schijven. De studie concludeert dat verspreidingsvoorschriften expliciet de afhankelijkheid van de schijfafsnijstraal moeten opnemen om de evolutie en levensduur van compacte schijven correct te voorspellen, waardoor een noodzakelijke correctie wordt aangebracht aan huidige populatiesynthesemodellen.