Abbondanze cosmiche ed evoluzione chimica galattica

Misurare la composizione dell’universo

Per sapere quanto idrogeno, elio, carbonio o ferro c’è nell’universo non possiamo raccogliere campioni da una stella e portarli in laboratorio. Quello che possiamo fare è leggere le impronte digitali che gli atomi lasciano nella luce: ogni elemento, riscaldato nell’atmosfera di una stella, assorbe la luce a frequenze specifiche, formando uno spettro di assorbimento con righe a posizioni esattamente prevedibili dalla meccanica quantistica dei livelli atomici. Misurando profondità, forma e larghezza Doppler di queste righe, e confrontandole con calcoli accurati di trasferimento radiativo in un modello di atmosfera stellare, possiamo dedurre quanto di quell’elemento è presente nella fotosfera, e — sotto l’assunzione standard di mescolamento — nell’intera stella sotto la zona convettiva esterna.

La fotografia più dettagliata che abbiamo della composizione “media” della materia dell’universo è quella del Sole — la stella più vicina e meglio misurata — e quella delle meteoriti CI condritiche, i campioni più primitivi del Sistema Solare. Sole e CI condriti sono in eccellente accordo per gli elementi non volatili (entro $\pm 10\%$ per la maggior parte degli elementi), una convergenza che riflette la formazione del Sistema Solare da gas e polveri di composizione omogenea e che costituisce il punto di riferimento di tutta l’astrofisica. Le quattro CI condriti note — Ivuna, Alais, Orgueil, Tonk, cadute fra il 1806 e il 1938 — sono state analizzate con tecniche di precisione progressivamente migliori (ICP-MS multicollector, attivazione neutronica, secondary ion mass spectrometry, accelerator mass spectrometry) e costituiscono il “campione di riferimento” delle abbondanze del Sistema Solare per gli elementi pesanti. Le differenze sistematiche fra Sole e CI sono ben caratterizzate: il Sole è leggermente più povero di Li, Be e B perché questi elementi vengono distrutti dalla combustione protonica nello strato sotto la zona convettiva durante la sequenza principale, mentre le CI sono più povere di volatili (H, C, N, O, Ne, Ar) perché questi sono sfuggiti durante la formazione e l’evoluzione delle parent bodies asteroidali.

Scale di abbondanza

Le abbondanze cosmiche sono espresse in scale diverse a seconda del contesto. La scala logaritmica spettroscopica è il riferimento per il confronto fra stelle: $\log \epsilon(\mathrm{X}) \equiv \log_{10}(N_{\mathrm{X}}/N_{\mathrm{H}}) + 12$, normalizzata a $\log \epsilon(\mathrm{H}) = 12$ per definizione. La frazione in massa $X_i = M_i/M_{\mathrm{tot}}$, con $\sum_i X_i = 1$, è naturale per i modelli stellari: convenzionalmente $X = X_\mathrm{H}$, $Y = X_\mathrm{He}$, $Z = \sum_{i > 2} X_i$ è la cosiddetta “metallicità”, e per il Sole moderno $X_\odot \approx 0{,}74$, $Y_\odot \approx 0{,}24$, $Z_\odot \approx 0{,}014$. La frazione in numero $Y_i = X_i/A_i$ (a volte indicata anche $X_i/A_i$ a seconda della convenzione) è quella naturale per scrivere le reazioni nucleari. Per le stelle individuali, la metallicità relativa al Sole è $[\mathrm{Fe/H}] = \log(N_\mathrm{Fe}/N_\mathrm{H})_\star - \log(N_\mathrm{Fe}/N_\mathrm{H})_\odot$, e analoghi rapporti $[\mathrm{X/Y}]$ misurano le abbondanze relative confrontandole con la corrispondente proporzione solare.

Le compilazioni standard delle abbondanze solari di riferimento sono tre, prodotte da gruppi diversi con metodologie parzialmente sovrapposte. La compilazione AGSS09 (Asplund, Grevesse, Sauval, Scott 2009) è stata la prima a integrare in modo sistematico le simulazioni 3D radiative-idrodinamiche dell’atmosfera solare e l’analisi NLTE (non-Local-Thermodynamic-Equilibrium) delle righe principali, e ha prodotto $Z_\odot = 0{,}0134$. L’aggiornamento AAGS21 [Asplund et al. 2021] raffina ulteriormente i modelli atomici per Fe-peak e per gli elementi pesanti, portando $Z_\odot = 0{,}0139$. La compilazione Lodders (2020) [Lodders 2020] si focalizza sulla rivisitazione meteoritica per gli elementi non volatili e produce $Z_\odot = 0{,}0156$, leggermente più alta. La compilazione storica GS98 (Grevesse-Sauval) e i suoi aggiornamenti [Grevesse et al. 2010] , basata su atmosfere 1D LTE pre-3D-NLTE, dà $Z_\odot = 0{,}0170$ — circa $0{,}13$-$0{,}15$ dex più alta di AAGS21 per CNO. Le differenze fra queste compilazioni — apparentemente piccole numericamente — hanno conseguenze profonde sui modelli solari interni, e sono al cuore del cosiddetto Solar Modeling Problem discusso nella sezione successiva.

La compilazione solare

Le abbondanze solari fotosferiche misurate nello strato esterno della stella, attraverso lo spettro che ci arriva sulla Terra, costituiscono il riferimento moderno dell’astrofisica. Ogni nuovo metodo di analisi le ha rivedute progressivamente al ribasso di una frazione, e dopo l’introduzione delle simulazioni 3D della convezione solare e dell’analisi NLTE delle righe atomiche e molecolari (Asplund, Grevesse, Sauval e collaboratori, anni 2000), le abbondanze più importanti — C, N, O in particolare — sono state ridotte di una frazione consistente rispetto ai valori usati per decenni nelle simulazioni stellari. Questa revisione è importante perché le abbondanze fotosferiche entrano direttamente come input dei modelli solari interni e dei modelli stellari in generale: cambiare $Z_\odot$ cambia l’opacità della zona radiativa, la struttura termica del core, e l’evoluzione predetta.

I valori delle abbondanze solari fotosferiche nelle tre compilazioni di riferimento per gli elementi principali sono riassunti nella tabella seguente:

La metallicità complessiva $Z$ è $0{,}0134$ in AGSS09, $0{,}0139$ in AAGS21, $0{,}0170$ in GS98. La riduzione delle abbondanze CNO è la principale conseguenza dell’analisi 3D NLTE: nei modelli 1D LTE precedenti, l’inomogeneità tridimensionale della convezione e la deviazione dalla distribuzione termica delle popolazioni atomiche nelle righe di assorbimento erano mascherate, portando sistematicamente a sovrastima delle abbondanze fotosferiche.

Le tecniche di misura variano significativamente da elemento a elemento. L’elio non è osservabile nella fotosfera solare (le sue righe sono troppo deboli a $T_{\mathrm{eff}} = 5777$ K) ed è misurato indirettamente dalla elioseismologia — l’inversione delle frequenze dei modi p di oscillazione solare osservati da reti come BiSON e GONG — oppure dalla spettroscopia delle prominenze coronali (a $T \gtrsim 10^{4}$ K). Carbonio, azoto e ossigeno sono misurati attraverso una combinazione di righe atomiche (CI, NI, OI permesse e proibite) e molecolari (CH, CN, CO, NH, OH); l’analisi richiede simulazioni 3D NLTE perché le righe sono saturate, sensibili al gradiente termico della fotosfera, e sensibili agli effetti di mescolamento turbolento. Il ferro fornisce la base della scala di metallicità grazie a centinaia di righe FeI e FeII disponibili nello spettro visibile. I gas nobili Ne e Ar non sono osservabili in fotosfera e sono misurati indirettamente dai venti solari o dalle regioni coronali, con coefficienti FIP (First Ionization Potential) da applicare per correggere il bias di abbondanza del vento. Per gli elementi a bassa abbondanza — Mn, V, Sc, e gli s/r-process come Eu, Ba — si usano poche righe singole o data points con incertezze maggiori ($\pm 0{,}05$-$0{,}15$ dex).

Una controversia tecnica resta aperta sull’abbondanza solare di ossigeno. AGSS09 e AAGS21 convergono su $\log\epsilon(\mathrm{O}) = 8{,}69$, in tensione con i requisiti del modello solare standard (vedi sezione seguente). Analisi alternative (Caffau et al. 2011 con CO5BOLD, Bergemann et al. 2021 con metodi NLTE indipendenti) suggeriscono valori $\sim 0{,}05$ dex più alti che ridurrebbero la tensione. Il consenso della comunità non è ancora consolidato, e questa ambiguità di $\sim 0{,}05$ dex sull’ossigeno solare è — per ironia del destino — uno dei numeri più importanti dell’astrofisica stellare quantitativa.

Il Solar Modeling Problem

C’è una tensione di vent’anni che ha animato il dibattito della comunità solare e stellare: il modello solare standard (SSM), basato sulle abbondanze AGSS09/AAGS21, prevede una struttura interna che non è in pieno accordo con i vincoli derivati dall’elioseismologia, ossia dalle frequenze di milioni di modi acustici di oscillazione osservati con precisione $< 10^{-4}$ sulla superficie del Sole. I modelli SSM che usano le vecchie abbondanze GS98 (con CNO più alto) si accordano significativamente meglio con tutti i diagnostici elioseismologici. La domanda aperta è: l’analisi 3D NLTE è sbagliata, o il modello solare è incompleto?

Il Solar Modeling Problem si manifesta in modo concordante su tre diagnostici indipendenti. La profondità della zona convettiva è $r_{\mathrm{cz}}/R_\odot = 0{,}7133 \pm 0{,}0005$ dall’elioseismologia; SSM AGSS09 predicono $0{,}722$, in tensione di $\sim 1{,}5\sigma$ rispetto al valore osservato, mentre SSM GS98 producono $0{,}7133$ in accordo perfetto. L’abbondanza di elio nella superficie convettiva $Y_s$ è $0{,}2485 \pm 0{,}0035$ dall’elioseismologia (inversione del small frequency separation); SSM AGSS09 predicono $0{,}2426$, mentre SSM GS98 predicono $0{,}2485$. Il profilo della velocità del suono $c(r)$ — ricavato da inversione globale di $\sim 5000$ frequenze p-mode della spettroscopia Sun-as-a-star — mostra differenze fino a $\sim 1\%$ fra il modello AGSS09 e i dati, soprattutto nella regione subito sotto la zona convettiva, mentre GS98 è coerente entro la precisione. SSM AAGS21 (con metallicità leggermente più alta di AGSS09) riduce la tensione ma non la risolve. SSM aggiornati (Vinyoles et al. 2017, Magg et al. 2022) sintetizzano questi vincoli sismici.

Le soluzioni proposte si articolano su tre direzioni complementari. La revisione delle opacità è il candidato principale: misure dirette dell’opacità del Fe a temperature solari da parte di Bailey et al. (2015) con la Z-machine di Sandia hanno mostrato valori $\sim 25\%$ più alti dei calcoli teorici OP e OPAL usati standardmente nei codici SSM. Una revisione completa delle opacità nella zona radiativa potrebbe assorbire larga parte della tensione, riducendo la sensibilità del SSM alle abbondanze esatte di CNO. Misurazioni indipendenti di opacità a Z-machine e in altri impianti high-energy-density sono in corso. La revisione dei coefficienti di diffusione atomica e di settling nei modelli stellari modifica la depleszione predetta di elio in superficie e potrebbe contribuire al gap di $Y_s$. L’aggiunta di overshooting convettivo alla base della zona convettiva sposta marginalmente $r_{\mathrm{cz}}$ predetta. E infine, ulteriori revisioni delle abbondanze fotosferiche di C, N, O nella direzione di valori più alti — come suggerito dalle analisi indipendenti di Caffau e Bergemann — potrebbero ridurre la tensione senza richiedere nuove fisiche.

Il diagnostico nucleare indipendente più forte sul Solar Modeling Problem viene dal flusso di neutrini CNO solari misurato da Borexino al Gran Sasso. La rivelazione storica del 2020 [Collaboration 2020] ha fornito la prima misura diretta dei neutrini prodotti dal ciclo CNO nel core solare, con flusso integrato $(6{,}6^{+2{,}0}_{-0{,}9}) \times 10^{8}$ cm$^{-2}$s$^{-1}$ al 68% C.L. — compatibile con SSM GS98 e in marginale tensione (circa $25\%$) con SSM AGSS09/AAGS21. La misura aggiornata di Borexino del 2022 [Collaboration 2022] con statistica accresciuta e analisi raffinata sposta il flusso a $6{,}7^{+1{,}2}_{-0{,}8} \times 10^{8}$ cm$^{-2}$s$^{-1}$ e conferma la preferenza per metallicità solare alta, pur senza essere decisiva. Il flusso CNO è direttamente proporzionale all’abbondanza di C+N nel core solare ed è dunque un sondaggio quasi diretto di $Z$ nella regione di produzione di energia. Future misure di precisione da Hyper-Kamiokande e da JUNO ridurranno ulteriormente le incertezze e contribuiranno a chiarire definitivamente la questione.

La risoluzione finale del Solar Modeling Problem richiede plausibilmente una combinazione di abbondanze parzialmente revisite verso valori più alti di CNO, nuove opacità che incorporino misure di laboratorio aggiornate, e raffinamenti minori della modellistica di mixing e diffusione. I prossimi cinque-dieci anni — con i risultati di Hyper-Kamiokande e JUNO sul fronte neutrinico, le campagne high-EDX di laboratorio sul fronte opacità, e le analisi di nuova generazione 3D NLTE sul fronte spettroscopia — promettono di chiudere il capitolo aperto da AGSS09 nel 2009.

Abbondanze in stelle e ISM

Una volta nota la composizione del Sole, possiamo confrontarla con quella di altre stelle della Galassia. Le stelle più vecchie e più povere di metalli — vissute miliardi di anni fa, in epoche in cui la Via Lattea era più giovane e meno arricchita dalle generazioni stellari precedenti — ci dicono come si è evoluta nel tempo la composizione chimica galattica. Stelle in regioni della Galassia con storie di formazione stellare diverse (alone, disco spesso, disco sottile, bulge, ammassi globulari) ci dicono dove e quando la Galassia ha accumulato i suoi metalli. Le abbondanze del gas interstellare attuale completano il quadro: confrontandole con quelle stellari, possiamo stimare quanta materia “leggera” e “pesante” è ancora libera nello spazio rispetto a quanta è stata incorporata in stelle ormai inerti o in resti compatti.

Le principali componenti chimiche della Via Lattea sono caratterizzate dai loro range di $[\mathrm{Fe/H}]$ e di $[\alpha/\mathrm{Fe}]$, e da età stellare media, secondo il quadro schematico seguente:

L’evoluzione del rapporto $[\alpha/\mathrm{Fe}]$ con $[\mathrm{Fe/H}]$ — il classico “ginocchio” della cronografia chimica galattica — è discussa in dettaglio nella seconda parte di questo capitolo, dedicata all’evoluzione chimica della Via Lattea. Qui ci limitiamo a notare che la mappa chimica della Galassia è oggi prodotta da survey spettroscopici di larga scala con risoluzione e numero di stelle senza precedenti. APOGEE (Sloan Digital Sky Survey) ha caratterizzato $\sim 10^{5}$ stelle in banda H con risoluzione $R \sim 22\,500$, coprendo il disco e il bulge. GALAH all’AAT ha caratterizzato $\sim 10^{5}$ stelle nell’ottico a $R \sim 28\,000$, focalizzato sul disco. La Gaia-ESO Survey ha caratterizzato $\sim 10^{4}$ stelle in ammassi a risoluzione alta. Le survey di prossima generazione — WEAVE al WHT, 4MOST al VISTA, PFS al Subaru — produrranno fra $10^{6}$ e $10^{7}$ stelle a $R = 5000$-$20000$ nel prossimo quinquennio. Per il gas interstellare, i sondaggi principali sono la spettroscopia di assorbimento delle DLA (sistemi Lyman-$\alpha$ damped a $z = 2$-$4$ contro la luce di QSO di sfondo), la spettroscopia di emissione delle regioni HII (nella Via Lattea e nelle galassie vicine), le diffuse interstellar bands per le molecole, e l’emissione X del gas caldo della corona galattica e dell’alone.

Alcuni dati specifici hanno un ruolo diagnostico particolare. Il plateau di Spite del litio in stelle metal-poor di alone — l’osservazione che le stelle EMP con $T_{\mathrm{eff}} > 5800$ K hanno tutte $A(\mathrm{Li}) \approx 2{,}2$, indipendentemente dalla metallicità — è in tensione con le predizioni di nucleosintesi del Big Bang combinate con i parametri cosmologici Planck (che predirebbero $A(\mathrm{Li}) \approx 2{,}7$): è il cosiddetto problema del litio cosmologico, discusso nei capitoli 2 e 8. Le stelle r-only, come la classica CS 22892-052 di Sneden, mostrano un pattern r-process universale per $Z = 56$-$80$ scalabile entro $\pm 0{,}1$ dex al pattern solare (vedi capitolo 6), e sono oggi catalogate in oltre 200 esemplari dalla R-Process Alliance. Le stelle Carbon-Enhanced Metal-Poor (CEMP) rappresentano circa il $25\%$ delle stelle con $[\mathrm{Fe/H}] < -2$ e si suddividono in sottotipi CEMP-no (probabilmente arricchite da supernovae Pop III con fallback parziale), CEMP-s (arricchite per mass transfer da una compagna AGB), CEMP-r (singolo evento r-process ancestrale) e CEMP-rs (combinazione complicata). I solar twins — stelle simili al Sole entro $\pm 100$ K in $T_{\mathrm{eff}}$, $\pm 0{,}1$ dex in $\log g$ e $\pm 0{,}1$ dex in $[\mathrm{Fe/H}]$ — permettono misure differenziali ad altissima precisione ($< 0{,}01$ dex) e hanno rivelato la cosiddetta Sun’s chemical anomaly (Meléndez et al. 2009): il Sole risulta marginalmente povero di elementi refrattari (Fe, Si, Mg, Al, Ti) rispetto ai solar twins di pari $[\mathrm{Fe/H}]$, di $\sim 0{,}04$ dex. L’interpretazione preferita è che i pianeti terrestri abbiano sequestrato materia refrattaria dalla nebulosa prima della completa formazione del Sole, lasciando il gas residuo che ha alimentato la fotosfera leggermente impoverito — una firma piccola ma misurabile della formazione planetaria.

Stelle metal-poor estreme e Pop III

Le stelle più antiche e più povere di metalli che riusciamo a osservare oggi ci raccontano una storia unica: probabilmente sono nate dal materiale arricchito da una singola supernova primordiale (o da un piccolo numero di esse), e portano la firma chimica diretta di quegli eventi nei loro spettri. Studiandole, ricostruiamo proprietà delle prime stelle dell’universo — le Population III — che non possiamo vedere direttamente perché vissute e morte miliardi di anni fa, quando l’universo aveva centinaia di milioni di anni e l’arricchimento chimico era appena iniziato.

Una manciata di stelle Extremely Metal-Poor (EMP, $[\mathrm{Fe/H}] < -3$) e Ultra Metal-Poor (UMP, $[\mathrm{Fe/H}] < -4$) sono state identificate negli ultimi vent’anni, principalmente dalla survey Hamburg/ESO Survey, dalla SDSS/SEGUE, e dalle survey australiane SkyMapper. Le più estreme sono riassunte nella tabella seguente:

I pattern di abbondanze di queste stelle suggeriscono progenitori SN core-collapse di massa $\sim 25$-$60\,M_\odot$ con mixing-fallback model (Umeda-Nomoto, Heger-Woosley): durante l’esplosione, i layer profondi di Si-burning vengono parzialmente mescolati prima di cadere nel buco nero residuo, ed espellono solo gli $\alpha$-elements più esterni con pochissimo Fe — coerente con $[\mathrm{Fe/H}] \ll -3$ ma $[\mathrm{X/Fe}]$ moderato o elevato per C, N, O, Mg. La firma chimica include tipicamente $[\mathrm{C/Fe}] \gtrsim +1$ in CEMP-no (dominante in UMP), assenza di arricchimenti r-process o s-process significativi, e un eccesso variabile di N e Mg compatibile con specifici modelli Pop III. I best-fit yield della stella SMSS 0313 — l’oggetto con il limite di metallicità più stringente noto — suggeriscono un progenitore Pop III di $\sim 50\,M_\odot$ in regime mixing-fallback con energia di esplosione intermedia.

L’evidenza di stelle Pop III pure (a metallicità rigorosamente zero) ancora vive resta incerta. Le Pop III di alta massa ($M > 10\,M_\odot$), nate dalla frammentazione di nubi primordiali HI-H$_2$ senza polvere e senza metalli, sono già morte da $> 13$ Gyr e contribuiscono solo come progenitori delle stelle EMP osservate. Le Pop III di bassa massa ($M < 0{,}8\,M_\odot$), se si fossero formate, sarebbero ancora vive oggi: la loro ricerca attraverso survey progressivamente più profonde non ha finora prodotto candidati. La SMSS 0313, con $[\mathrm{Fe/H}] < -7{,}3$, ha pur sempre carbonio e azoto rivelati spettroscopicamente, e quindi non è Pop III pura. L’interpretazione preferita è che le Pop III di bassa massa non si siano mai formate in numero significativo, perché la frammentazione delle nubi primordiali a basse masse richiede un meccanismo di raffreddamento efficiente (polvere, metalli, molecole pesanti) che era assente prima del primo arricchimento — il cosiddetto fragmentation threshold. JWST sta cercando galassie a $z = 10$-$15$ con firme spettrali compatibili con popolazioni Pop III dominanti, e alcune candidate sono state riportate con interpretazione preliminare; il quadro definitivo richiederà ulteriori spettri ad alta risoluzione nei prossimi anni.

Stato della questione: abbondanze osservate

La caratterizzazione delle abbondanze cosmiche è oggi un programma maturo, con compilazioni di riferimento stabili al livello dell’incertezza statistica e con sistematiche residue dell’ordine di $0{,}05$-$0{,}1$ dex su molti elementi chiave. Restano aperte tre questioni quantitative principali. La prima è la convergenza definitiva del $Z_\odot$: la riconciliazione fra AAGS21, Lodders 2020, e i requisiti SSM si gioca su una decina di centesimi di dex nelle abbondanze CNO, ed è il cuore del Solar Modeling Problem. La seconda è il completamento del catalogo delle stelle EMP/UMP: identificare le poche centinaia di oggetti più estremi nell’alone galattico, nelle dwarf satellite e nei flussi stellari, è essenziale per ricostruire la funzione di massa iniziale e gli yield della prima generazione stellare. La terza è la mappatura chimica completa della Galassia: le survey 4MOST e WEAVE producano nei prossimi cinque anni un catalogo di $10^{7}$ stelle con $> 20$ elementi misurati per ciascuna, che permetterà la ricostruzione data-driven della storia chimica e dinamica della Via Lattea.

Le prospettive a 5-10 anni sono concrete e convergenti. Sul fronte solare, Hyper-Kamiokande (presa dati dal 2028) e JUNO (dal 2025) misureranno con precisione il flusso di neutrini $^{8}\mathrm{B}$ e CNO solari, fornendo vincoli diretti sul $Z$ del core solare e contribuendo alla risoluzione del SMP. Sul fronte spettroscopico, ELT, GMT e TMT permetteranno la spettroscopia ad altissima risoluzione di stelle EMP a $V \sim 18$-$20$ — un ordine di grandezza più deboli di quelle accessibili oggi — e l’identificazione di nuove stelle UMP candidate Pop III a redshift $z \sim 0$. Sul fronte spaziale, Gaia DR4 (atteso 2026-2027) completerà la cinematica 6D di $10^{9}$ stelle in combinazione con i metalli da APOGEE/GALAH/4MOST/WEAVE, fornendo la base per la chemodynamics auto-consistente della Galassia. JWST e i futuri telescopi di nuova generazione caratterizzeranno spettralmente galassie a $z = 8$-$15$ e potenzialmente identificheranno popolazioni Pop III ancora vive nei sistemi più primitivi accessibili. La combinazione di questi fronti promette di portare la cosmochimica osservativa da disciplina di rifinitura a disciplina di mappatura completa nel prossimo decennio.

Le abbondanze cosmiche misurate nel Sole, nelle meteoriti, nelle stelle e nel gas interstellare sono il banco di prova finale di tutta la nucleosintesi: i meccanismi descritti nei capitoli precedenti — combustioni quiescenti, processi s e r, p, nucleosintesi esplosiva di supernovae, kilonove di NSM — devono insieme riprodurre quantitativamente la curva delle abbondanze solari e la sua evoluzione nel tempo galattico. La seconda parte di questo capitolo è dedicata a come questi siti vengono integrati in modelli di evoluzione chimica galattica che fanno predizioni dirette sui rapporti $[\mathrm{X/H}]$ vs $[\mathrm{Fe/H}]$ osservati, e su come la combinazione di yield stellari, IMF, SFR e delay time distribution delle SN Ia plasma la firma chimica della Via Lattea che vediamo oggi.

La galassia come storia di generazioni

Una galassia non ha sempre avuto la composizione che ha oggi. Quando la Via Lattea era giovane, il suo gas conteneva quasi soltanto idrogeno ed elio, con tracce minime di litio dalla nucleosintesi primordiale (capitolo 2): tutto il resto della tavola periodica è stato cucinato e disperso dalle stelle nei tredici miliardi di anni successivi. Ogni generazione di stelle che nasce, vive e muore arricchisce un po’ di più il gas interstellare di “metalli” — nel gergo astrofisico, qualsiasi elemento più pesante dell’elio. Le stelle della generazione successiva nascono dal gas più ricco e quindi sono più metalliche alla nascita. Studiando la composizione di stelle di età diverse e di regioni diverse della Galassia possiamo leggere questa storia direttamente nelle proprietà spettroscopiche delle loro fotosfere: vediamo l’evoluzione chimica della Via Lattea incisa nelle abbondanze chimiche delle sue stelle, dal Sistema Solare attuale fino alle stelle EMP più antiche dell’alone con $[\mathrm{Fe/H}] < -4$.

Le equazioni della GCE

L’evoluzione chimica galattica (GCE, Galactic Chemical Evolution) è descritta in forma quantitativa da un set di equazioni di conservazione su una regione galattica considerata come “scatola” o, nei modelli più sofisticati, come collezione di “zone” connesse fra loro. Nel modello a singola zona, le equazioni base sono:

$$\frac{d \sigma_{\mathrm{gas}}}{dt} = -\psi(t) + E(t) + I(t) - W(t)$$ $$\frac{d \sigma_\star}{dt} = \psi(t) - E(t)$$

dove $\sigma_{\mathrm{gas}}$ è la densità superficiale di gas, $\sigma_\star$ di stelle, $\psi(t)$ il tasso di formazione stellare (SFR), $E(t)$ il rate di restituzione di massa dalle stelle morte (ejecta SN, venti AGB), $I(t)$ il infall di gas (accrezione cosmologica e satellite), $W(t)$ il outflow (galactic winds da SN feedback). Per ciascun elemento $i$ tracciato, l’evoluzione obbedisce a un’equazione analoga ma pesata per la frazione in massa:

$$\frac{d (\sigma_{\mathrm{gas}} X_i)}{dt} = -\psi(t) X_i(t) + E_i(t) + I_i(t) X_{i,\mathrm{infall}} - W_i(t) X_i(t),$$

dove il termine di ejecta chimica è

$$E_i(t) = \int_{M_\ast}^{M_{\mathrm{up}}} (M - M_R)\,Y_i(M)\,\psi(t - \tau(M))\,\xi(M)\,dM,$$

con $Y_i(M)$ yield della stella di massa $M$ per l’elemento $i$, $\xi(M)$ funzione di massa iniziale (IMF), $\tau(M)$ lifetime stellare, e $M_R$ massa del residuo compatto. L’integrale è il cuore quantitativo del modello: trasforma la storia di formazione stellare e la fisica stellare individuale in una predizione misurabile sull’evoluzione chimica del gas. Gli ingredienti principali sono l’IMF (Salpeter $\xi \propto M^{-2{,}35}$, Kroupa, Chabrier — quest’ultima usata standardmente nei modelli moderni); la legge SFR-storia (per esempio Schmidt-Kennicutt $\psi \propto \sigma_{\mathrm{gas}}^{1{,}4}$ o assunzioni più dirette su SFH); i yield stellari (Sukhbold-Woosley per SN core-collapse, Nomoto-Kobayashi per SN Ia, Karakas-Lugaro [Karakas & Lattanzio 2014] e Cristallo et al. [Cristallo et al. 2015] per AGB); i modelli di infall e outflow; e la delay time distribution (DTD) delle SN Ia. Per una trattazione formale e completa si rinvia alle monografie di Pagel [Pagel 2009] e Matteucci [Matteucci 2021] .

I codici GCE pubblici di riferimento includono NuPyCEE (Côté et al., Python single-zone con NSM, novae, super-AGB), OMEGA (Côté et al., multi-zona con database SFH), gce4py (Andrews et al., Python con yield personalizzabili) e Chempy (Rybizki et al., framework Bayesiano per parametri GCE). Tipicamente integrano $\sim 50$-$1000$ specie chimiche tracciate su mesh di $10$-$100$ zone, con runtime di poche ore per modello completo della Via Lattea. La rassegna metodologica recente della disciplina è Matteucci (2021) [Matteucci 2021] , mentre il testo di riferimento didattico resta Pagel (2009) [Pagel 2009].

Dalle scatole chiuse alle galassie reali

I primi modelli analitici usavano una scatola chiusa: niente infall, niente outflow, mescolamento istantaneo, e talvolta l’approssimazione di riciclo istantaneo. Restano utili perché mostrano come apparirebbe l’evoluzione chimica se l’arricchimento fosse puramente interno — e falliscono in modi istruttivi. Una scatola chiusa predice troppe stelle longeve a bassa metallicità rispetto a quante se ne osservano nei dintorni solari: è il classico G-dwarf problem. La soluzione non è una piccola correzione agli yield; richiede l’infall di gas, che diluisce il mezzo interstellare e sostiene la formazione stellare nel tempo.

Gli outflow sono altrettanto importanti nelle galassie di piccola massa. Il feedback delle supernovae può espellere gas arricchito dai potenziali gravitazionali poco profondi: è per questo che le galassie nane hanno metallicità più basse delle grandi spirali a parità di frazione di gas, e che i loro ginocchi $[\alpha/\mathrm{Fe}]$ cadono a $[\mathrm{Fe/H}]$ più basso. Un modello che ignora i venti può adattarsi male al disco della Via Lattea e fallire del tutto per le sferoidali nane.

Il ginocchio $[\alpha/\mathrm{Fe}]$

Una delle firme più caratteristiche dell’evoluzione chimica galattica è la curva del rapporto $\alpha$-elements/Fe in funzione della metallicità. Per stelle molto antiche e povere di metalli ($[\mathrm{Fe/H}] \ll -1$), il rapporto $[\alpha/\mathrm{Fe}]$ è elevato — circa $+0{,}3$-$+0{,}4$ — perché soltanto le supernovae core-collapse hanno avuto il tempo di esplodere e arricchire il gas con la loro firma chimica $\alpha$-ricca (capitolo 5). Per stelle progressivamente più giovani e più metalliche, $[\alpha/\mathrm{Fe}]$ scende monotonicamente verso zero, perché le SN Ia — che producono molto Fe ma pochi $\alpha$-elements — entrano in scena con un delay time di alcune centinaia di milioni di anni e contribuiscono progressivamente di più al budget di Fe del gas interstellare. Il punto in cui la curva piega — il ginocchio o knee — segna l’epoca in cui le SN Ia hanno iniziato a contribuire al pool di Fe, e la sua posizione esatta è un cronometro chimico diretto della storia di formazione stellare.

Quantitativamente, il rapporto $[\alpha/\mathrm{Fe}]$ funziona così. A $[\mathrm{Fe/H}] \ll [\mathrm{Fe/H}]_{\mathrm{knee}}$ solo le SN core-collapse hanno contribuito e $[\alpha/\mathrm{Fe}] \approx +0{,}3$-$+0{,}4$ forma un plateau. A $[\mathrm{Fe/H}] \gg [\mathrm{Fe/H}]_{\mathrm{knee}}$ le SN Ia dominano la produzione di Fe e $[\alpha/\mathrm{Fe}] \to 0$, raggiungendo il valore solare. La posizione del knee $[\mathrm{Fe/H}]_{\mathrm{knee}}$ dipende dalla SFR storica dell’ambiente: una SFR alta significa che molte SN core-collapse hanno arricchito velocemente il gas in Fe prima che le SN Ia entrassero in scena, e dunque il knee cade a $[\mathrm{Fe/H}]$ più alto; una SFR bassa, al contrario, lascia tempo alle SN Ia di iniziare a contribuire mentre la metallicità è ancora bassa, e il knee cade a $[\mathrm{Fe/H}]$ basso. La posizione osservata del knee in diversi ambienti galattici è dunque un indicatore quantitativo della SFR media:

Il delay time distribution (DTD) delle SN Ia è il parametro fisico chiave che determina la posizione del knee per una data SFR. Le osservazioni in cluster di galassie (Maoz, Mannucci 2012) e l’analisi statistica di SN Ia in galassie host con SFR misurata supportano una DTD a power-law $\Psi(\tau) \propto \tau^{-1}$ per $\tau > \tau_{\min} \approx 40$ Myr, motivata fisicamente dall’inspiral gravitazionale di binarie WD nello scenario Double Degenerate (capitolo 5). I modelli GCE che riproducono quantitativamente il knee della Via Lattea richiedono questa DTD. Gli yield delle SN Ia restano incerti di $\pm 30\%$ per Fe e $\pm 50\%$ per Mn, Cr, e i diversi scenari di esplosione (Chandrasekhar W7 deflagrazione, Chandra DDT, sub-Chandrasekhar double detonation) producono yield diversi. Il rapporto Mn/Fe è il diagnostico più sensibile a separare i canali: alto per Chandra DDT (cattura $e^{-}$ efficiente a $\rho > 10^{9}$ g/cm$^{3}$), basso per sub-Chandra a densità centrale minore — e il pattern Mn/Fe osservato in funzione di $[\mathrm{Fe/H}]$ nelle stelle di disco e di alone vincola progressivamente la partizione fra canali.

Yield nucleari come input

Il modello GCE è accurato soltanto quanto sono accurati i yield di ciascun tipo di stella morente: SN core-collapse, SN Ia, AGB, NSM, novae. Questi yield sono calcolati dai modelli di nucleosintesi stellare discussi nei capitoli precedenti, e le incertezze sugli yield sono — di gran lunga — le incertezze dominanti della GCE moderna. Le tabelle yield in uso nei modelli GCE attuali coprono i contributi principali:

• SN core-collapse: Sukhbold-Woosley (2018) per $13$-$120\,M_\odot$ via codice KEPLER; Limongi-Chieffi (2018) FRANEC con rotazione differenziale; Pignatari-Herwig (2016) NuGrid.
• SN Ia: Iwamoto W7 (1999) per Chandrasekhar deflagrazione; Travaglio et al. (2004); Seitenzahl 3D-DDT (2013); Shen-Townsley per sub-Chandra.
• AGB: Karakas-Lugaro (2014) Monash [Karakas & Lattanzio 2014] ; Cristallo et al. (2015) FRUITY [Cristallo et al. 2015] ; Ventura-D’Antona (2018) ATON.
• NSM: Wanajo et al. (2014) per kilonova dettagliata; varianti su distribuzione $Y_e$ dell’ejecta.
• Novae: José-Hernanz (2011); Denissenkov et al. (2014).

Gli yield integrati su IMF (Kroupa o Salpeter) per generazione stellare unitaria — definiti come $\eta_X = \int Y_X(M)\,\xi(M)\,dM / \int M\,\xi(M)\,dM$ — sono tabulati in Kobayashi, Karakas e Lugaro (2020) [Kobayashi et al. 2020] e sono input diretto per la GCE moderna. La partizione fra contributori dominanti per gli elementi principali è riassunta nella tabella seguente:

La sensitività dei modelli GCE agli yield individuali è significativa e variabile da elemento a elemento. Le yield di C, N, O sono sensibili al trattamento di mass loss in stelle massicce e di Hot Bottom Burning (HBB) nelle AGB di massa intermedia, con variazione $\pm 30\%$ fra modelli. La yield di Mn dipende dal canale SN Ia dominante e produce un rapporto Mn/Fe osservabile distinguibile fra stelle Pop II (basso, sub-Chandra dominante) e Pop I (alto, Chandra DDT contribuisce). La yield di Eu dipende dal mix NSM/collapsar e dal DTD delle NSM: lo scatter di $[\mathrm{Eu/Fe}]$ vs $[\mathrm{Fe/H}]$ a basse metallicità è diagnostico diretto della rarità degli eventi r-process individuali. La yield di Pb (componente s strong) dipende dalla mass loss AGB a bassa metallicità e dal $^{13}\mathrm{C}$ pocket (capitolo 4). La procedura iterativa standard di calibrazione è: yield come input → modello GCE → confronto con $[\mathrm{X/Fe}]$ vs $[\mathrm{Fe/H}]$ osservato in survey large → revisione di yield ambigui → iterazione. Il vincolo è progressivamente raffinato dall’aggiunta di dati indipendenti come quelli dei grani presolari (capitolo 4) e della spettroscopia ad alta risoluzione di stelle vicine.

Modelli multi-zona della Via Lattea

I primi modelli di evoluzione chimica trattavano la Via Lattea come una scatola omogenea — un’idealizzazione utile per i primi calcoli pioneristici (Tinsley negli anni ‘70-‘80), ma manifestamente inadeguata alla galassia reale. Oggi sappiamo che la nostra galassia ha una struttura ricca: un disco sottile e un disco spesso con cinematica e chimica diverse, un alone interno e uno esterno, un bulge con la propria storia evolutiva, ammassi globulari con composizioni a multipla popolazione, e una serie di galassie satellite in vari stadi di accrescimento. Ognuno di questi ambienti ha una storia di formazione stellare diversa e una composizione chimica diversa di conseguenza. I modelli moderni dividono la Galassia in zone multiple, ciascuna con il proprio bilancio gas-stelle e con scambi (infall di gas, outflow da SN, mixing radiale per migrazione stellare) fra zone vicine.

Il modello multi-zona standard della Via Lattea (Matteucci-Romano 2002 e successori) suddivide la Galassia in halo + bulge + disco con SFR e DTI distinti, e il disco in $\sim 20$ anelli concentrici fra $R = 2$ kpc e $R = 22$ kpc, ciascuno con SFR e infall tempo-dipendenti. Il mixing radiale per migrazione stellare è implementato in modo ad-hoc, calibrato sui gradienti chimici osservati. I risultati notevoli del modello includono il gradiente radiale di metallicità osservato $d[\mathrm{Fe/H}]/dR \approx -0{,}05$ dex/kpc nel disco esterno — spiegato naturalmente da SFR e infall dipendenti dal raggio, con il disco interno più rapido nella formazione e dunque chimicamente più evoluto. L’inside-out formation vede il bulge e il disco interno formarsi per primi su tempo scala di $\sim 1$ Gyr, il disco esterno per ultimo su tempo scala di $\sim 7$ Gyr, in accordo con le età stellari mappate da Gaia DR3 e da APOGEE. Il modello two-infall prevede una prima fase rapida di accrescimento (formazione halo + thick disk in $\tau \sim 1$ Gyr) seguita da una seconda fase lenta (formazione thin disk in $\tau \sim 7$ Gyr), e spiega quantitativamente il “salto” osservato di $[\alpha/\mathrm{Fe}]$ fra thick e thin disk e la doppia sequenza in $[\alpha/\mathrm{Fe}]$ vs $[\mathrm{Fe/H}]$ vista in APOGEE e GALAH.

Sul fronte cosmologico, le simulazioni idrodinamiche zoom-in di galassie tipo Via Lattea — EAGLE, FIRE, IllustrisTNG, AURIGA, NIHAO — integrano N-body, idrodinamica SPH/AMR/MFM, e chimica sub-grid che traccia $\sim 10$-$20$ specie chimiche con yield integrati IMF, su volumi cosmologici dal collasso primordiale fino a oggi. Riproducono qualitativamente l’evoluzione chimica galattica osservata: l’origine del thick disk come prodotto di early star formation interna riscaldata da merger gas-rich o instabilità del disco; l’inomogeneità chimica attesa a basse $[\mathrm{Fe/H}]$ con scatter ampio coerente con le osservazioni di stelle EMP; la radial migration in cui stelle si muovono in raggio nel tempo per scattering con braccia spirali e barra, mescolando popolazioni chimiche di raggi diversi e producendo dispersione $[\mathrm{Fe/H}]$ a età fissata. Le simulazioni di alta risoluzione (FIRE-3, AURIGA Plus) raggiungono risoluzione di $\sim 100$ pc nel disco galattico e tracciano dettagliatamente la struttura chimica del disco — un livello di realismo che permette di confrontare quantitativamente con le mappe Gaia + APOGEE/GALAH della Via Lattea attuale.

Galassie satellite e flussi di accrezione

La Via Lattea non si è formata in isolamento: nei suoi tredici miliardi di anni di vita ha inghiottito decine di galassie più piccole, accrescendone le stelle e (in epoche più antiche) il gas. Studiando le firme chimiche e cinematiche delle stelle nell’alone — combinando spettroscopia ad alta risoluzione con i moti propri e le parallassi misurate da Gaia — gli astronomi possono oggi distinguere le stelle “indigene” nate nella Via Lattea originaria da quelle “immigrate” accresciute da galassie satellite ormai dissolte o ancora in fase di mixing dinamico. Questa decomposizione, impensabile prima di Gaia DR2 (2018), è una delle frontiere più attive dell’archeologia galattica.

Le galassie satellite identificate e parzialmente accresciute includono diversi sistemi notevoli. Sagittarius dwarf spheroidal (Sgr dSph) è in fase di accrescimento attivo, con stream stellari ben visibili che attraversano l’alone in più gradi del cielo. Gaia-Enceladus (o Sausage) è una galassia accresciuta circa 10 Gyr fa, identificata da Belokurov e Helmi nel 2018 attraverso la combinazione di pattern $[\alpha/\mathrm{Fe}]$-$[\mathrm{Fe/H}]$ e cinematica retrograda nell’alone interno; le sue stelle contribuiscono significativamente alla popolazione del thick disk e dell’alone interno, e il sistema progenitore aveva probabilmente massa $\sim 10^{10}\,M_\odot$ all’epoca di accrescimento — comparabile alla LMC attuale. Sequoia, Thamnos, Helmi streams e altri stream più piccoli sono stati successivamente identificati con tecniche analoghe. Le splashed stars del thick disk — stelle inizialmente formate nel disco e poi diffuse verticalmente in orbite di alone da un evento di accrescimento massiccio — completano il quadro della popolazione dell’alone interno.

Le implicazioni per la GCE sono profonde. La composizione chimica del thick disk e dell’alone interno non è quella di una popolazione formata in situ in un’evoluzione singola, ma una mescolanza di popolazioni eterogenee con storie di formazione stellare diverse. Gli yield integrati che osserviamo riflettono dunque la combinazione di SFH della Via Lattea originaria e SFH delle galassie satellite accrese. Le stelle Gaia-Enceladus, per esempio, mostrano $[\mathrm{Fe/H}] \approx -1{,}3$ medio (più povero del thick disk in situ), $[\alpha/\mathrm{Fe}] \approx +0{,}3$ (coerente con formazione rapida prima del knee SN Ia), e un pattern $[\mathrm{Eu/Fe}]$ specifico distinguibile da quello del thick disk in situ — un quadro chimicamente coerente con un sistema dwarf di massa intermedia comparabile alle dSph sopravvissute (Sculptor, Fornax). L’approccio del chemical tagging — identificazione di gruppi di stelle co-formate attraverso la firma chimica multi-elementale ($> 10$ elementi) — è la frontiera promettente per disaggregare la storia di accrescimento della Galassia, e le survey APOGEE e GALAH dedicano una frazione significativa del loro tempo osservativo proprio a questo programma.

I problemi aperti

Il quadro generale dell’evoluzione chimica galattica è oggi consolidato nelle sue linee maestre, ma resta una serie di problemi quantitativi aperti che animano la ricerca corrente. Tre famiglie di questioni sono particolarmente vive.

La prima è il problema dell’azoto primario. L’azoto è prodotto principalmente come primario (cioè da carbonio e ossigeno sintetizzati nella stessa stella) attraverso il Hot Bottom Burning nelle AGB di massa intermedia ($M = 4$-$7\,M_\odot$), che ha tempo scala di $\sim 100$-$300$ Myr. A bassissima metallicità ($[\mathrm{Fe/H}] < -3$), queste AGB non hanno ancora avuto tempo di evolvere e contribuire al gas interstellare quando si formano le prime stelle metal-poor, e quindi $[\mathrm{N/Fe}]$ dovrebbe essere molto basso o addirittura $-\infty$. Le osservazioni invece mostrano $[\mathrm{N/Fe}] \approx 0$ in molte stelle EMP, con scatter significativo. Le soluzioni proposte includono la rotazione veloce in stelle massicce Pop III (Hirschi 2007 e successori), che attiva mixing rotazionale di N attraverso il CNO-cycle nelle shell di He durante la sequenza principale, e la frammentazione di nubi primordiali con segnatura specifica (Sharda 2021). Il problema non è ancora risolto in modo definitivo.

La seconda è il problema dei plateau di $[\alpha/\mathrm{Fe}]$. I rapporti $[\mathrm{O/Fe}]$, $[\mathrm{Mg/Fe}]$, $[\mathrm{Si/Fe}]$ osservati nelle stelle metal-poor formano plateau notevolmente piatti intorno a $+0{,}4$ per $-3 < [\mathrm{Fe/H}] < -1$, con scatter osservato $< 0{,}1$ dex. Questa piattezza richiede o yield SN core-collapse molto stabili in funzione della massa progenitrice (in tensione con quanto predetto dai modelli stellari, che mostrano variazione $\sim 0{,}3$ dex), o un mescolamento efficace del gas su scale large che cancelli le fluttuazioni stocastiche evento per evento. A $[\mathrm{Fe/H}] < -3$, lo scatter cresce significativamente ($\pm 0{,}3$ dex), suggerendo che le stelle EMP più estreme abbiano effettivamente visto pochi eventi di nucleosintesi — talvolta una singola SN II progenitrice — la cui firma chimica individuale è ancora leggibile nello spettro.

La terza è il problema delle stelle Pop III mancanti. I modelli teorici di formazione delle prime stelle (Bromm, Yoshida, Hirano, Stacy) predicono che le Pop III avessero IMF top-heavy con masse caratteristiche $M_\ast \sim 10$-$300\,M_\odot$ in conseguenza dell’assenza di metalli e polvere che innesca la frammentazione delle nubi primordiali. Le Pop III di altissima massa ($140 < M < 260\,M_\odot$) terminano la loro vita come pair-instability supernovae (PISN) con yield molto distintivi: nessun elemento del picco Fe più pesante del Ni, $^{56}\mathrm{Ni}$ prodotto fino a $\sim 10\,M_\odot$ per evento, pattern $[\mathrm{Mg/Ca}]$ anomalo e $[\mathrm{Cr/Mn}]$ alto, assenza completa di processo r. Nessuna stella metal-poor finora osservata mostra in modo convincente questa firma chimica: la candidata storica era la stella UMP SDSS J102915+172927 (“Caffau star”), ma analisi successive hanno mostrato che il suo pattern è più coerente con un progenitore SN core-collapse standard di $\sim 30\,M_\odot$ con fallback. Possibili spiegazioni includono PISN molto rare (forse $< 1\%$ dei progenitori EMP), o Pop III dominate dalla massa intermedia ($10$-$40\,M_\odot$) che ha fatto la maggior parte del lavoro di arricchimento primordiale, o ancora la completa assenza di Pop III di bassa massa dal sample di stelle vive osservate (la massa minima per sopravvivere a 13 Gyr è $M < 0{,}8\,M_\odot$, e le simulazioni di frammentazione primordiale non producono routinariamente stelle a queste masse).

Stato della questione: evoluzione chimica

Il quadro della GCE moderna è solido al primo ordine: i meccanismi di nucleosintesi sono identificati e i loro yield calcolati, l’evoluzione del gas è ben descritta dalle equazioni di conservazione, le firme chimiche osservate nelle stelle riflettono coerentemente la storia di formazione stellare degli ambienti galattici noti. La sintesi metodologica di riferimento è oggi Nomoto, Kobayashi e Tominaga (2013) [Nomoto et al. 2013] per la parte yield, e Matteucci (2021) [Matteucci 2021] per la parte modellistica complessiva. Restano aperte questioni quantitative significative nei tre problemi discussi sopra, e una serie di raffinamenti su rotazione stellare, binarità, IMF dipendente da metallicità, e mixing 3D in atmosfere stellari che modificano marginalmente le predizioni.

Le prospettive a 5-10 anni sono concrete e multidirezionali. Sul fronte osservativo, le survey 4MOST, WEAVE, PFS completeranno il mapping spettroscopico della Via Lattea con $10^{6}$-$10^{7}$ stelle e $> 20$ elementi per stella, raffinando i vincoli GCE su scale dell’intera Galassia e identificando migliaia di nuove candidate EMP/UMP per archeologia stellare. Gaia DR4 (atteso 2026-2027) fornirà la cinematica 6D di precisione per tutto il dataset, abilitando la chemodynamics auto-consistente della popolazione galattica. JWST sta cercando galassie a $z = 10$-$15$ con firme spettrali compatibili con popolazioni Pop III dominanti — candidate come GN-z11 e i sistemi rivelati nei survey JADES e CEERS sono in fase di analisi spettroscopica dettagliata. I telescopi 30-metro di nuova generazione ELT, GMT e TMT permetteranno spettroscopia ad altissima risoluzione di stelle Pop II in galassie satellite della Via Lattea e in ammassi globulari, accedendo a metallicità un ordine di grandezza più basse di quelle accessibili oggi. Sul fronte teorico, i modelli idrodinamici di nuova generazione (FIRE-3, AURIGA Plus, IllustrisTNG-50) raggiungono risoluzione tale da risolvere chimicamente la struttura del disco e la storia di accrescimento delle galassie satellite, e i modelli stellari di nuova generazione con rotazione e binarità self-consistent producono yield più affidabili. La combinazione di questi fronti promette di portare la GCE da disciplina di sintesi qualitativa a disciplina di archeologia quantitativa nel prossimo decennio, capace di ricostruire eventi specifici e popolazioni progenitrici dalla firma chimica delle stelle osservate.

L’evoluzione chimica galattica costituisce il quadro statistico in cui i singoli siti di nucleosintesi — combustioni quiescenti, processi s e r, nucleosintesi esplosiva, kilonove — convergono nella firma chimica delle stelle e del gas che osserviamo nella Via Lattea attuale. Resta da raccontare cosa resta da capire nella disciplina nel suo complesso: le frontiere aperte della nucleosintesi stellare, dalla fisica nucleare oltre il drip-line alle applicazioni cosmologiche, sono il tema del prossimo capitolo che conclude la prima parte sistematica del libro.