Can photoevaporation open gaps in protoplanetary discs?

Probleemstelling
De dispersie van protoplanetaire schijven is een cruciale factor bij het bepalen van de timing en de omstandigheden van planeetvorming. Hoewel fotoverdamping — het verwijderen van gas via thermisch gedreven winden die worden gelanceerd door hoogenergetische stellaire straling — algemeen wordt erkend als een belangrijk dispersiemechanisme, blijft de theoretische basis voor hoe dit gaten opent en in stand houdt onzeker. Standaardmodellen vertrouwen doorgaans op eendimensionale (1D) viskeuze evolutieberekeningen waarbij de fotoverdampingsmassaverliesvoet, $\dot{\Sigma}(r)$, wordt behandeld als een vaste sinkterm afgeleid van statische tweedimensionale (2D) simulaties. Deze benadering gaat er impliciet van uit dat het massaverliesprofiel constant blijft, zelfs wanneer de schijfstructuur evolueert. De auteurs stellen dat dit een cruciale feedbackloop negeert: naarmate een gat begint te vormen, verandert de lokale schijfstructuur, wat potentieel de condities voor het lanceren van de wind verandert, lokale massaverliessnelheden vermindert en viskeuze instroom mogelijk maakt om de uitgeputte regio aan te vullen. Als fotoverdamping op deze wijze zelflimiterend is, kan het minder efficiënt zijn in het leegmaken van gas op de radii waar reuzenplaneten ontstaan dan voorheen aangenomen, met significante implicaties voor kernaccretie- en migratiemodellen.

Methodologie
Om te onderzoeken of fotoverdamping alleen gaten kan openen en in stand houden wanneer de stroming dynamisch reageert op de schijfstructuur, gebruikten de auteurs een computationele aanpak in twee fasen:

1. Fase I (1D Evolutie): Ze gebruikten de viskeuze evolutiecode spock om een schijf te simuleren die onderhevig is aan röntgenfotoverdamping tot het begin van gatvorming. De simulatie evolueerde een gasoppervlakte-dichtheidsprofiel ($\Sigma$) over 2,1 Myr, waarbij een massaverliesprofiel van Picogna et al. (2019) werd opgenomen. Deze fase stelde de initiële condities vast voor de 2D-simulaties, specifiek een staat waarin een gedeeltelijke uitputting in de oppervlakte-dichtheid was ontstaan.
2. Fase II (2D Radiatieve Hydrodynamica): Het oppervlakte-dichtheidsprofiel uit de 1D-simulatie bij 2,1 Myr werd in een axiaal symmetrisch 2D ($r-\theta$) sferisch coördinatenstelsel gemapt met behulp van de pluto-code. Deze simulaties modelleerden de fotoverdampingsuitstroom en de interactie ervan met de schijf.
   • Model p-g: Maakte gebruik van de geëvolueerde oppervlakte-dichtheid uit Fase I, die een reeds bestaande dichtheidsuitputting bevatte.
   • Model p-ng: Maakte gebruik van een aangepast initieel oppervlakte-dichtheidsprofiel waarbij de uitputting werd onderdrukt (vervangen door een machtswet) om de impact van de initiële morfologie op gatopening te beoordelen.
   • De simulaties maakten gebruik van een line-of-sight verwaringsbenadering en namen radiatief evenwicht aan, waarbij radiatieve transfer ontkoppeld werd van hydrodynamica om de computationele kosten te verlagen.

Belangrijkste Resultaten

• Zelflimiterende Gatvorming: In de 2D-simulaties nam, zodra een dichtheidsdepressie ontstond, de lokale massaverliesvoet scherp af binnen de gatregio. Dit gebeurde omdat de uitgeputte regio in de schaduw van de binnenste schijf bleef, wat de blootstelling aan stellaire bestraling verminderde. Als gevolg hiervan verdiepte het gat zich niet tot een volledig uitgeholde holte.
• Gedeeltelijke Aanvulling: Viskeuze instroom en radiale massatransport langs het schijfoppervlak (en door de wind) werkten deels als aanvulling op de uitgeputte regio. Het systeem evolueerde naar een persistente, gedeeltelijk uitgeputte zone in plaats van een schone binnenste holte.
• Ongevoeligheid voor Initiële Condities: Het vergelijken van modellen met (p-g) en zonder (p-ng) een initiële uitputting toonde aan dat de langetermijn evolutie grotendeels onafhankelijk was van de initiële schijfmorfologie. In beide gevallen reguleerde het systeem zichzelf naar een staat van gedeeltelijke uitputting.
• Beperkingen van 1D Voorschriften: De auteurs testten een eerstegraads voorschrift om het 2D-gedrag te benaderen in 1D-modellen door het primordiale massaverliesprofiel te schalen op basis van de lokale oppervlakte-dichtheidsuitputting. Hoewel dit de overeenkomst verbeterde ten opzichte van statische voorschriften, evolueerden de 1D-modellen nog steeds te snel en produceerden ze diepere gaten dan de 2D-simulaties. Deze discrepantie benadrukt het onvermogen van 1D-modellen om radiaal massatransport langs het schijfoppervlak en de complexe interactie van schaduwvorming te vatten.
• Observationele Signaturen: Hoewel het gasgat gedeeltelijk gevuld blijft, kan de drukmaximum aan de buitenrand van de depressie nog steeds fungeren als een efficiënte val voor stofkorrels. Dit zou observationele signaturen kunnen produceren die lijken op transitieschijven (binnenholtes in continuümemissie), zelfs zonder volledige gasverwijdering.

Betekenis en Claims
Het artikel daagt het standaardparadigma uit dat fotoverdamping efficiënt schone binnenste holtes uitslijpt en direct transitieschijven produceert. In plaats daarvan stellen de auteurs voor dat fotoverdamping kan leiden tot een quasi-stationaire toestand die wordt gekenmerkt door een persistente, gedeeltelijk uitgeputte zone. Dit suggereert dat schijven die een matige uitputting vertonen zonder volledige binnenste verduistering, systemen kunnen vertegenwoordigen die beïnvloed worden door fotoverdamping maar niet volledig zijn gedispergeerd.

De auteurs bieden een bescheiden, eerstegraads voorschrift om dit zelfregulerende gedrag te benaderen in 1D-schijfevolutiemodellen, wat geschikt is voor gebruik in studies naar planeetvorming en populatiesynthese. Ze stellen echter expliciet dat dit voorschrift slechts benaderend is en de multidimensionale aard van de stroming, met name de radiale massaverdeling, niet volledig kan vatten. Het artikel concludeert dat verdere multidimensionale simulaties en exploratie van de parameterruimte noodzakelijk zijn om robuuste recepten af te leiden voor globale schijf- en planeetpopulatiemodellen. Het werk onderstreept de noodzaak om verder te gaan dan statische massaverwerking om de interactie tussen interne winden en paden van planeetvorming te begrijpen.