FORGE'd in the Early Universe: The Effect of Protostellar Outflows on Pop III Accretion

Sintesi Tecnica: Simulazione FORGE'd in FIRE di una (proto)stella di Pop. III

Problema
La formazione e l'evoluzione della prima generazione di stelle (Popolazione III) rimangono scarsamente vincolate a causa della mancanza di osservazioni dirette e dell'immensa gamma dinamica richiesta per connettere la formazione delle strutture cosmologiche alla dettagliata fisica protostellare. I precedenti modelli teorici si basavano spesso su condizioni iniziali idealizzate, omettevano la fisica critica (come campi magnetici, radiazione e getti), o non riuscivano a produrre funzioni di massa iniziale (IMF) e molteplicità stellari realistiche osservate negli ambienti ricchi di metalli. Esiste un dibattito specifico sul fatto che le stelle di Pop. III possano lanciare getti protostellari, un processo che dipende dall'amplificazione di campi magnetici traccia e dalla presenza di un disco circumstellare rotante e magnetizzato.

Metodologia
Questo studio presenta una simulazione di radiazione magneto-idrodinamica (RMHD) cosmologica "zoom-in" utilizzando il framework FORGE'd in FIRE. La simulazione traccia la formazione e l'evoluzione di una singola (proto)stella di Pop. III priva di metalli a redshift $z \sim 14$ all'interno di un contesto completamente cosmologico.

• Framework di Simulazione: Viene utilizzato il codice GIZMO, che impiega il metodo meshless finite-mass (MFM) per risolvere le equazioni accoppiate di gravità, magneto-idrodinamica (MHD) e trasporto della radiazione.
• Fisica Inclusa: La simulazione incorpora il modello fisico FIRE-3, che include:
  • Radiazione: Trasporto multifrequenza (ionizzante, Lyman-Werner, ecc.) accoppiato in modo auto-consistente con gas e polvere.
  • Termochimica: Estesa a temperature di $1\text{--}10^{10}$ K e densità di $10^{-8}\text{--}10^{24}$ cm$^{-3}$, tracciando le abbondanze di H, He, HD e polvere.
  • Campi Magnetici: Inizializzati con traccia campi semente cosmologici ($\sim 10^{-15}$ G) che si amplificano in modo auto-consistente. Gli effetti di MHD non ideale (resistività ohmica, effetto Hall, diffusione ambipolare) sono trascurati, giustificato dal confronto delle scale temporali nel post-processing che mostra come siano trascurabili rispetto alle scale temporali dinamiche nel gas privo di metalli.
  • Evoluzione Protostellare: Particelle "sink" rappresentano le (proto)stelle, evolvendo attraverso l'accrezione, i percorsi di sequenza principale e i meccanismi di feedback.
• Modellazione di Jet e Venti: I getti protostellari sono modellati seguendo Cunningham et al., deviando una frazione ($f_W = 0.3$) della massa accresciuta in deflussi bipolari collimati lanciati a una frazione ($f_K = 0.3$) della velocità Kepleriana. I venti stellari sono trascurati per le stelle prive di metalli poiché scalano con la metallicità.
• Risoluzione: La simulazione utilizza una "iper-raffinazione", risolvendo le celle di gas fino a $\sim 10^{-3} M_\odot$ e scale spaziali di $\sim 10^{-4}$ pc ($< 100$ au) vicino alla (proto)stella, mantenendo al contempo il contesto cosmologico da $z \sim 100$.

Risultati Chiave

1. Crescita e Massa Stellare:

   • La (proto)stella cresce da un tasso di accrezione iniziale di $\sim 5 \times 10^{-3} M_\odot \text{yr}^{-1}$ a una massa finale di $\sim 27 M_\odot$ in 31.000 anni.
   • Il tasso di crescita si stabilizza mentre la stella entra nella fase di "combustione del guscio di He" (Shell D burning). La massa finale è regolata dal feedback (getti e radiazione) piuttosto che dall'esaurimento del serbatoio circostante.
   • Al termine della simulazione, la stella ha un raggio di $\sim 10 R_\odot$, una temperatura effettiva di $\sim 3.5 \times 10^4$ K e una luminosità di $\sim 1.3 \times 10^5 L_\odot$. La scala temporale di Kelvin-Helmholtz ($t_{KH}$) diventa più breve della scala temporale di accrezione ($t_{acc}$), indicando che la stella raggiungerà la sequenza principale prima di accrescere una massa significativa.

2. Disco Circumstellare:

   • Si forma un disco supportato dalla rotazione, che si estende fino a $\sim 0.01$ pc ($\sim 2000$ au).
   • Il disco rimane gravitazionalmente stabile durante tutta la simulazione, con un parametro di Toomre $Q$ costantemente $>> 1$.
   • La stabilità è mantenuta da una combinazione di supporto termico (nelle regioni interne, $T \sim 10^4$ K) e supporto turbolento (nelle regioni esterne, $R > 1000$ au).
   • Il profilo della densità superficiale transita da $\Sigma \propto R^{-1}$ nel disco esterno a un profilo a disco-$\alpha$ nelle regioni interne, con un brusco aumento della densità a $\sim 20$ au.

3. Jetti Protostellari e Outflow:

   • Lancio dei Jetti: I getti si formano perché il rapporto massa-flusso magnetico ($\mu$) cade nell'intervallo ($\mu \sim 10^{-1} \text{--} 10^1$) che permette il lancio di getti, come definito da Machida e Doi. Ciò avviene nonostante l'uso di campi magnetici semente cosmologici traccia.
   • Dinamica: I getti si propagano fino a $\sim 1$ pc, collidendo con il gas in caduta per formare shock di terminazione.
   • Regolazione: I getti guidano un feedback di quantità di moto sostanziale, spazzando via il gas e riducendo il tasso netto di accrezione verso l'interno di circa un ordine di grandezza alle grandi scale radiali. Questo feedback guidato dalla quantità di moto regola la massa stellare finale.
   • Il tasso di outflow raggiunge il picco a $\sim 10^{-3} M_\odot \text{yr}^{-1}$ entro i primi 1.000 anni prima di declinare, mentre il tasso di afflusso si stabilizza intorno a $3\text{--}5 \times 10^{-4} M_\odot \text{yr}^{-1}$.

4. Struttura Chimica e Termica:

   • Nel disco interno, le alte temperature ($T \sim 10^4$ K) sopprimono l'idrogeno molecolare ($H_2$) a favore dell'idrogeno ionizzato, spostando i meccanismi di raffreddamento da $H_2$ al raffreddamento di Lyman.
   • Il disco rimane otticamente sottile alla sua radiazione di raffreddamento dominante, permettendo ai fotoni IR di sfuggire efficientemente.

Significato e Rivendicazioni
Il documento rivendica di presentare la prima applicazione del framework FORGE'd in FIRE a una (proto)stella di Pop. III risolta individualmente all'interno di un setup cosmologico "zoom-in" completo. La significatività primaria risiede nel dimostrare che:

• I getti sono vitali nella formazione di Pop. III: Anche con campi magnetici semente cosmologici traccia, le condizioni per il lancio di forti getti sono soddisfatte, sfidando le precedenti assunzioni secondo cui i getti potrebbero essere assenti in ambienti privi di metalli.
• Il feedback regola la massa: I getti protostellari giocano un ruolo critico nel plasmare la massa finale delle stelle di Pop. III, regolando l'accrezione attraverso deflussi guidati dalla quantità di moto, piuttosto che attraverso il semplice consumo dell'intero serbatoio da parte della stella.
• Stabilità del Disco: La formazione di un disco circumstellare stabile e non frammentante è possibile in questo regime, mantenuta da supporto termico e turbolento.

Questi risultati fornisiscono un quadro fisicamente più completo delle prime regioni di formazione stellare, offrendo spunti sulla Funzione di Massa Iniziale (IMF) di Pop. III e sui processi di feedback. Gli autori notano che il lavoro futuro espanderà questo approccio a un gruppo di stelle di Pop. III per caratterizzare pienamente la prevista IMF e la molteplicità stellare.