Can photoevaporation open gaps in protoplanetary discs?

Problematica
La dispersione dei dischi protoplanetari è un fattore critico nel determinare la tempistica e le condizioni della formazione planetaria. Sebbene la fotoevaporazione — la rimozione di gas tramite venti guidati termicamente e lanciati da radiazione stellare ad alta energia — sia ampiamente riconosciuta come un meccanio chiave di dispersione, la base teorica su come essa apra e mantenga i gap rimane incerta. I modelli standard si basano tipicamente su calcoli di evoluzione viscosa monodimensionali (1D), in cui il tasso di perdita di massa fotoevaporativa, $\dot{\Sigma}(r)$, è trattato come un termine di pozzo fisso derivato da simulazioni bidimensionali (2D) statiche. Questo approccio assume implicitamente che il profilo di perdita di massa rimanga costante anche mentre la struttura del disco evolve. Gli autori sostengono che ciò trascuri un cruciale ciclo di feedback: man mano che un gap inizia a formarsi, la struttura locale del disco cambia, alterando potenzialmente le condizioni di lancio del vento, riducendo i tassi di perdita di massa locali e permettendo l'afflusso viscoso di riempire la regione impoverita. Se la fotoevaporazione è auto-limitante in questo modo, potrebbe essere meno efficiente nel liberare il gas ai raggi in cui si formano i pianeti giganti rispetto a quanto precedentemente assunto, con implicazioni significative per i modelli di accrezione del nucleo e migrazione.

Metodologia
Per investigare se la fotoevaporazione da sola possa aprire e mantenere i gap quando il flusso risponde dinamicamente alla struttura del disco, gli autori hanno impiegato un approccio computazionale in due fasi:

1. Fase I (Evoluzione 1D): Hanno utilizzato il codice di evoluzione viscosa spock per simulare un disco soggetto a fotoevaporazione da raggi X fino all'inizio della formazione del gap. La simulazione ha evoluto un profilo di densità superficiale del gas ($\Sigma$) su 2,1 Myr, incorporando un profilo di perdita di massa da Picogna et al. (2019). Questa fase ha stabilito le condizioni iniziali per le simulazioni 2D, specificamente uno stato in cui si era sviluppato un parziale impoverimento nella densità superficiale.
2. Fase II (Idrodinamica Radiativa 2D): Il profilo di densità superficiale dalla simulazione 1D a 2,1 Myr è stato mappato in un sistema di coordinate sferiche assialsimmetriche 2D ($r-\theta$) utilizzando il codice pluto. Queste simulazioni hanno modellato l'uscita fotoevaporativa e la sua interazione con il disco.
   • Modello p-g: Ha utilizzato la densità superficiale evoluta dalla Fase I, che includeva un preesistente impoverimento di densità.
   • Modello p-ng: Ha utilizzato un profilo di densità superficiale iniziale modificato in cui l'impoverimento era soppresso (sostituito da una legge di potenza) per valutare l'impatto della morfologia iniziale sull'apertura del gap.
   • Le simulazioni hanno impiegato un'approssimazione di riscaldamento lungo la linea di vista e hanno assunto l'equilibrio radiativo, disaccoppiando il trasferimento radiativo dall'idrodinamica per ridurre il costo computazionale.

Risultati Chiave

• Formazione di Gap Auto-Limitante: Nelle simulazioni 2D, una volta formata una depressione di densità, il tasso di perdita di massa locale è diminuito bruscamente all'interno della regione del gap. Ciò è accaduto perché la regione impoverita è rimasta nell'ombra del disco interno, riducendo l'esposizione all'irraggiamento stellare. Di conseguenza, il gap non si è approfondito in una cavità completamente evacuata.
• Riempimento Parziale: L'afflusso viscoso e il trasporto di massa radiale lungo la superficie del disco (e attraverso il vento) hanno agito per riempire parzialmente la regione impoverita. Il sistema si è evoluto verso una zona persistente e parzialmente impoverita piuttosto che verso un foro interno pulito.
• Insensibilità alle Condizioni Iniziali: Il confronto tra i modelli con (p-g) e senza (p-ng) un impoverimento iniziale ha mostrato che l'esito evolutivo a lungo termine era ampiamente indipendente dalla morfologia iniziale del disco. In entrambi i casi, il sistema si è autoregolato verso uno stato di parziale impoverimento.
• Limitazioni delle Prescrizioni 1D: Gli autori hanno testato una prescrizione di primo ordine per approssimare il comportamento 2D nei modelli 1D scalando il profilo di perdita di massa primordiale in base all'impoverimento della densità superficiale locale. Sebbene ciò abbia migliorato l'accordo rispetto alle prescrizioni statiche, i modelli 1D sono evoluti ancora troppo rapidamente e hanno prodotto gap più profondi rispetto alle simulazioni 2D. Questa discrepanza evidenzia l'incapacità dei modelli 1D di catturare il trasporto di massa radiale lungo la superficie del disco e la complessa interazione dell'ombreggiamento.
• Firme Osservative: Sebbene il gap di gas rimanga parzialmente pieno, il massimo di pressione al bordo esterno della depressione può ancora agire come una trappola efficiente per i grani di polvere. Ciò potrebbe produrre firme osservative simili ai dischi di transizione (cavità interne nell'emissione del continuo) anche senza una completa rimozione del gas.

Significato e Rivendicazioni
L'articolo sfida il paradigma standard secondo cui la fotoevaporazione scava efficientemente cavità interne pulite e produce direttamente dischi di transizione. Invece, gli autori propongono che la fotoevaporazione possa portare a uno stato quasi stazionario caratterizzato da una zona persistente e parzialmente impoverita. Ciò suggerisce che i dischi che esibiscono un moderato impoverimento senza una completa rimozione interna possano rappresentare sistemi influenzati dalla fotoevaporazione ma non completamente dispersi.

Gli autori forniscono una modesta prescrizione di primo ordine per approssimare questo comportamento autoregolato nei modelli di evoluzione del disco 1D, che è adatta per l'uso negli studi di formazione planetaria e sintesi di popolazione. Tuttavia, dichiarano esplicitamente che questa prescrizione rimane approssimativa e non può catturare completamente la natura multidimensionale del flusso, in particolare la ridistribuzione della massa radiale. L'articolo conclude che sono necessarie ulteriori simulazioni multidimensionali ed esplorazioni dello spazio dei parametri per derivare ricette robuste per i modelli globali di disco e di popolazione planetaria. Il lavoro sottolinea la necessità di andare oltre i trattamenti statici della perdita di massa per comprendere l'interazione tra i venti interni e i percorsi di formazione planetaria.