Formation of Isotopically Heterogeneous Molecular Cloud Cores in Filamentary Molecular Clouds

Sintesi Tecnica: Formazione di nuclei di nubi molecolari isotopicamente eterogenei in nubi molecolari filamentose

Definizione del problema
L'analisi delle meteoriti rivela che i solidi del sistema solare primordiale, come le inclusioni calcio-alluminio ricche (CAI), esibiscono un'eterogeneità isotopica (ad esempio, variazioni nei rapporti $^{54}\text{Cr}/^{52}\text{Cr}$ e $^{50}\text{Ti}/^{47}\text{Ti}$). Sebbene siano state proposte spiegazioni basate sulla lavorazione termica nei dischi protoplanetari e l'ereditarietà da nuclei di nubi molecolari madri eterogenee, rimane incerto l'entità con cui l'eterogeneità isotopica sopravviva al collasso delle nubi molecolari filamentose in nuclei densi. Nello specifico, è necessario quantificare quanto dell'iniziale variazione isotopica presente in una nube filamentosa venga preservata o livellata durante la frammentazione gravitazionale e la formazione dei nuclei, e come questa eterogeneità si traduca nel conseguente disco circumstellare.

Metodologia
Gli autori hanno utilizzato simulazioni Smoothed Particle Hydrodynamics (SPH), impiegando specificamente il metodo magnetoidrodinamico Godunov smoothed particle (GSPM), per modellare la frammentazione gravitazionale di nubi molecolari filamentose.

• Configurazione della simulazione: Le simulazioni hanno modellato filamenti con una massa lineare critica ($M_{\text{line,crit}} \approx 18 M_\odot \text{pc}^{-1}$) e una larghezza di $\sim 0,1$ pc, coerenti con le osservazioni. Le condizioni iniziali includevano un profilo di densità in equilibrio idrostatico, turbolenza di Kolmogorov ($\sigma = 2c_s$) e un campo magnetico iniziale di $10 \mu\text{G}$ perpendicolare all'asse del filamento. Sono stati eseguiti anche casi puramente idrodinamici per confronto.
• Tracciamento isotopico: Le specie isotopiche sono state trattate come traccianti passivi. Gli autori hanno costruito sei modelli distinti (Tabella 1) per imporre gradienti di rapporto isotopico iniziale ($C_0$) lungo diversi assi: l'asse longitudinale ($z$), gli assi minori ($x, y$) e combinazioni di essi. Il rapporto isotopico è stato normalizzato tra 0 e 1 sull'intero dominio computazionale.
• Analisi: L'evoluzione del rapporto isotopico è stata monitorata seguendo le particelle SPH. Gli autori hanno definito gusci sferici attorno al picco di densità dei nuclei in formazione per calcolare il rapporto isotopico medio ($\bar{C}$) e la differenza dal centro del nucleo ($|\Delta C|$). Hanno inoltre impiegato un approccio semi-analitico per stimare lo spostamento iniziale del centro di massa dei gusci ($\delta z_{\text{ini}}$) basandosi sui campi di velocità turbolenta e sugli spettri di potenza della densità.

Risultati chiave

1. Sopravvivenza dell'eterogeneità: Le simulazioni dimostrano che le eterogeneità isotopiche non vengono completamente livellate durante la frammentazione del filamento. I nuclei conservano circa lo 0,1% - 1% dell'iniziale eterogeneità isotopica presente lungo la lunghezza iniziale del filamento di 1 pc.
2. Dipendenza geometrica: L'entità dell'eterogeneità conservata dipende dall'orientamento del gradiente isotopico iniziale rispetto alla geometria del filamento.
   • I gradienti lungo l'asse longitudinale (direzione $z$, Modello $M_z$) producono le maggiori differenze isotopiche tra i gusci ($|\Delta C| \sim 0,01$).
   • I gradienti lungo gli assi minori ($x$ o $y$) producono differenze significativamente minori. Ciò è attribuito alla geometria del filamento, che restringe lo spostamento degli elementi di fluido lungo gli assi minori, limitando così il mescolamento di materiale con diversi rapporti isotopici.
3. Meccanismo di eterogeneità: Il principale motore della residua differenza isotopica tra i gusci è la differenza iniziale nel centro di massa dei gusci, causata dal campo di velocità turbolenta. Le stime semi-analitiche, che tengono conto dello spostamento iniziale degli elementi di fluido, concordano con i risultati della simulazione entro un fattore di 2, confermando che la configurazione turbolenta iniziale determina il profilo isotopico finale.
4. Influenza del campo magnetico: La presenza di un campo magnetico ($10 \mu\text{G}$) non altera significativamente il profilo interno del rapporto isotopico rispetto al caso puramente idrodinamico.
5. Implicazioni per il disco circumstellare: Stimando il raggio centrifugo del materiale in accrezione, gli autori hanno scoperto che l'eterogeneità isotopica persiste nel disco circumstellare. Sebbene il disco erediti l'eterogeneità del nucleo, la distribuzione anisotropa del momento angolare implica che il gas accrezzi da vari gusci sferici. I risultati suggeriscono che circa l'1% dell'iniziale eterogeneità filamentosa rimanga nel disco.

Significatività e affermazioni
L'articolo sostiene di fornire un collegamento quantitativo tra le anomalie isotopiche osservate nelle meteoriti e le condizioni astrofisiche della nube molecolare madre. Combinando le proprie scoperte con il fattore di riduzione riportato da Jacquet et al. (2019) (il quale suggerisce una riduzione di un fattore $\sim 3$ dal nucleo al disco), gli autori inferiscono che l'eterogeneità residua nel disco circumstellare finale è approssimativamente dello 0,03% - 0,3% della variazione iniziale del filamento.

Gli autori affermano che questo livello di eterogeneità residua è sufficiente a spiegare le anomalie isotopiche osservate nei materiali del Sistema Solare (ad esempio, le variazioni di $^{54}\text{Cr}/^{52}\text{Cr}$ e $^{50}\text{Ti}/^{47}\text{Ti}$ di 1–10 parti su 10.000) senza richiedere processi termici estremi o scenari di iniezione nucleosintetica esotica post-collasso. Lo studio conclude che l'eterogeneità isotopica dei componenti del sistema solare potrebbe naturalmente sopravvivere all'intero processo dalla frammentazione del filamento alla formazione del disco, a condizione che la nube madre possedesse un'appropriata eterogeneità iniziale. Gli autori osservano che sono necessarie simulazioni ad alta risoluzione per il lavoro futuro per risolvere pienamente le proprietà realistiche dei rapporti isotopici nel disco circumstellare, in particolare riguardo all'infall e agli outflow anisotropi.