Modelling dust coagulation, dynamical drag and turbulent mixing during star and disc formation

Sintesi Tecnica: Modellazione della Coagulazione della Polvere, del Drag Dinamico e del Mixing Turbolento durante la Formazione di Stelle e Dischi

Definizione del Problema
La polvere gioca un ruolo critico sia nella formazione delle stelle che in quella dei pianeti, influenzando il raffreddamento del gas, i livelli di ionizzazione e l'accoppiamento dei campi magnetici con il gas. Sebbene la coagulazione della polvere in grani più grandi (pebbles) sia essenziale per la formazione di planetesimali e la crescita dei nuclei, le precedenti simulazioni idrodinamiche tridimensionali della formazione stellare hanno ampiamente trascurato l'evoluzione dinamica della polvere rispetto al gas. I metodi precedenti, come quelli di Bate (2022), trattavano le popolazioni di polvere come se evolvessero indipendentemente in ogni elemento fluido (particella SPH), assumendo che la polvere fosse perfettamente accoppiata al gas. Tale approssimazione era valida per le fasi iniziali e a bassa densità del collasso, ma non riusciva a tenere conto della migrazione, del settling (assestamento) e del mixing della polvere su scale temporali più lunghe, particolarmente all'interno dei dischi protoplanetari dove le densità sono più elevate e la crescita dei grani è rapida. Inoltre, i modelli esistenti spesso mancavano di un trattamento auto-consistente di come la turbolenza guidi sia le velocità relative necessarie per la crescita dei grani, sia la diffusione turbolenta che mescola la polvere verticalmente.

Metodologia
Gli autori presentano una nuova implementazione all'interno del codice di smoothed particle hydrodynamics (SPH) SPHNG che combina tre distinti processi fisici:

1. Coagulazione della Polvere: Utilizzando il metodo discretizzato dell'equazione di Smoluchowski a campo medio (Hirashita & Yan 2009; Bate 2022), il codice traccia l'evoluzione di una popolazione di grani di polvere attraverso $N$ classi dimensionali (size bins) distribuite logaritmicamente per ogni particella di gas SPH. Il metodo considera collisioni binarie con una probabilità di adesione pari a unità al di sotto di una specifica soglia di velocità di coagulazione.
2. Drag Dinamico (MULTIGRAIN): Gli autori integrano l'algoritmo MULTIGRAIN (Hutchison et al. 2018), che risolve le equazioni del fluido continuo per una miscela di gas e molteplici specie di polvere nell'approssimazione della velocità terminale. Ciò consente il calcolo delle velocità differenziali macroscopiche tra polvere di diverse dimensioni e il gas a causa dei gradienti di pressione e della gravità.
3. Integrazione Implicita e Diffusione Turbolenta: Una chiave innovazione è l'uso di un metodo di integrazione implicita recentemente sviluppato (Elsender & Bate 2024) per risolvere le equazioni del drag della polvere. Questo evita i vincoli restrittivi del passo temporale (time-step) tipici dei metodi espliciti per i grani con tempi di arresto lunghi e previene frazioni di polvere negative non fisiche. Inoltre, gli autori introducono un nuovo metodo implicito per modellare la diffusione turbolenta. Trattano il rimescolamento turbolento della polvere come un processo di diffusione ($D_i = \alpha_{SS} c_s H / Sc$), dove il numero di Schmidt ($Sc$) tiene conto dell'efficienza di accoppiamento di diverse dimensioni dei grani. Questo termine di diffusione viene risolto utilizzando un approccio di operator-splitting di tipo Gauss-Seidel, permettendo alla polvere di mescolarsi tra gli elementi fluidi.

Il codice traccia l'evoluzione della polvere dal collasso di un nucleo di nube molecolare, attraverso la fase di primo nucleo idrostatico, fino alla formazione e all'evoluzione di un disco protoplanetario, utilizzando particelle "sink" per seguire il sistema oltre la formazione del nucleo stellare.

Contributi Chiave

• Prima Simulazione 3D Accoppiata: Questo lavoro rappresenta la prima simulazione idrodinamica tridimensionale a modellare simultaneamente la coagulazione della polvere, il drag dinamico del gas e il mixing turbolento durante le prime fasi di formazione di stelle e dei loro dischi.
• Schema di Diffusione Implicito: Lo sviluppo di un robusto solver implicito per la diffusione turbolenta della polvere consente simulazioni stabili del mixing della polvere senza i severi limiti di passo temporale associati agli schemi di diffusione espliciti.
• Trattamento Unificato: Il codice unifica il trattamento delle velocità relative (che guidano la crescita) e della diffusione turbolenta (che guida il mixing) all'interno di un unico framework, utilizzando lo stesso parametro di turbolenza ($\alpha_{SS}$) per entrambi i processi nei casi di test.

Risultati
Gli autori presentano i risultati sia dai calcoli di test (settling della polvere in una colonna verticale) che dalle simulazioni complete di formazione stellare:

• Calcoli di Test: Il codice recupera con successo le soluzioni analitiche per l'equilibrio tra il settling della polvere e la diffusione turbolenta. I test confermano che l'inclusione della diffusione turbolenta impedisce la formazione di strati di polvere irrealisticamente sottili e mantiene un'altezza di scala (scale height) della polvere finita.
• Simulazioni di Formazione Stellare:
  • Primo Nucleo Idrostatico: In linea con il lavoro precedente (Bate 2022), si verifica una rapida crescita della polvere all'interno del primo nucleo idrostatico, producendo popolazioni localmente monodisperse con dimensioni superiori a 100 $\mu$m all'interno del nucleo. L'inclusione di drag e diffusione ha avuto un impatto trascurabile su questa fase precoce a causa dei brevi tempi evolutivi.
  • Evoluzione del Disco Protoplanetario: Una volta formato il disco, l'inclusione della diffusione turbolenta altera significativamente l'evoluzione della polvere.
    • Crescita Potenziata: La diffusione turbolenta potenzia sostanzialmente il tasso di crescita dei grani di polvere. Fornendo continuamente grani piccoli e di medie dimensioni alle regioni ad alta densità (dove i grani grandi stanno crescendo), la diffusione impedisce l'esaurimento locale del "cibo" per la coagulazione.
    • Dipendenza dalla Turbolenza: Il livello di turbolenza ($\alpha_{SS}$) è critico. Una turbolenza elevata ($\alpha_{SS} = 10^{-3}$) porta a una rapida crescita verso popolazioni monodisperse di grani di $\sim30-100$ $\mu$m in tutto il disco. Una turbolenza inferiore ($\alpha_{SS} = 10^{-4}$) risulta in una crescita più lenta e in una popolazione di polvere più stratificata verticalmente, con i grani più grandi confinati più vicino al piano medio (mid-plane).
    • Settling della Polvere: Nelle simulazioni con turbolenza debole, il settling della polvere porta a un aumento significativo del rapporto polvere-gas nel piano medio (di oltre il 50%), mentre i grani piccoli rimangono impoveriti vicino al piano medio e confinati negli strati superficiali.
  • Drag vs. Diffusione: I confronti tra simulazioni con e senza specifici processi rivelano che il drag dinamico polvere-gas ha un effetto limitato sulle distribuzioni di polvere in queste scale temporali, poiché i grani rimangono ben accoppiati. Al contrario, l'inclusione della diffusione turbolenta ha un effetto enorme, livellando le distribuzioni di polvere e accelerando la coagulazione.

Significatività e Dichiarazioni
L'articolo sostiene che questa nuova metodologia consente, per la prima volta, la modellazione auto-consistente della crescita, del mixing e della migrazione della polvere durante la formazione di stelle e dei loro dischi. Gli autori sottolineano che, sebbene l'inclusione di drag e diffusione non alteri i risultati fondamentali riguardanti la crescita della polvere nel primo nucleo idrostatico, essa è cruciale per l'evoluzione successiva del disco protoplanetario. Nello specifico, affermano che la diffusione turbolenta della polvere fornisce un potenziamento sostanziale al tasso di crescita dei grani di polvere, agendo come una fonte di grani piccoli nelle regioni in cui i grani grandi si stanno formando.

Gli autori rimangono modesti riguardo all'ambito delle loro scoperte, osservando che:

• I risultati dipendono dal modello di turbolenza adottato (ad esempio, $\alpha$-disc vs. turbolenza MHD).
• Le simulazioni considerano attualmente la coagulazione ma non la frammentazione; diverse soglie di coagulazione potrebbero portare alla frammentazione e al ripopolamento di grani piccoli.
• I calcoli di trasferimento radiativo utilizzano ancora opacità interstellari standard invece di utilizzare le distribuzioni locali evolutive della dimensione dei grani, un passaggio identificato come una necessaria estensione futura.

In definitiva, il lavoro fornisce uno strumento per confrontare le simulazioni idrodinamiche 3D con la ricchezza di nuovi dati osservativi provenienti da strutture come ALMA, VLT-SPHERE e JWST, in particolare riguardo al settling della polvere e al rimescolamento turbolento nei dischi visti di taglio (edge-on).