ECHO21: a tool for modelling global 21-cm signal from dark ages to reionization

Sintesi Tecnica di "ECHO21: Uno strumento per modellare il segnale globale a 21 cm dai Dark Ages alla riionizzazione"

Problema e Motivazione
I "Dark Ages" (Età Oscure) e il "Cosmic Dawn" (Alba Cosmica) sono periodi del cronogramma cosmico scarsamente compresi a causa della mancanza di sonde osservative dirette. La transizione iperfine a 21 cm dell'idrogeno neutro offre una potenziale sonda per svelare la fisica delle prime stelle, del mezzo intergalattico (IGM) e dei parametri cosmologici. Sebbene i codici numerici e semi-numerici come 21cmFAST e BEORN forniscano previsioni accurate e fisicamente ricche, essi sono computazionalmente lenti, ostacolando un'efficiente inferenza astrofisica o cosmologica da grandi dataset. Al contrario, i codici analitici esistenti come ARES e Zeus21 offrono velocità ma mancano di certi componenti fisici, specificamente gli effetti di riscaldamento dei fotoni Lyman-$\alpha$ (Ly$\alpha$). Gli autori identificano che l'omissione del riscaldamento da Ly$\alpha$ può portare a forti bias nella modellazione del segnale. Per rispondere alla necessità di uno strumento veloce, sufficientemente realistico e flessibile, capace di gestire ampi spazi di parametri per l'inferenza, gli autori introducono ECHO21 (Exploring the Cosmos with Hydrogen Observation).

Metodologia
ECHO21 è un pacchetto Python che utilizza un framework analitico per modellare il segnale globale a 21 cm da un redshift $z \approx 1500$ fino al presente. Il codice opera su un modello a due zone dell'IGM, distinguendo tra l'IGM bulk e le regioni H II completamente ionizzate.

• Framework Fisico: Il codice risolve equazioni differenziali accoppiate per la temperatura cinetica del gas ($T_k$), la frazione elettronica dell'IGM bulk ($x_e$) e il fattore di riempimento volumetrico delle regioni H II ($Q$).

  • Storia Termica: L'evoluzione di $T_k$ tiene conto del raffreddamento adiabatico, del riscaldamento Compton, del riscaldamento X e, crucialmente, del riscaldamento da Ly$\alpha$. Gli autori notano che, mentre i precedenti codici analitici ignoravano il riscaldamento da Ly$\alpha$, ECHO21 lo include, calcolando il tasso di riscaldamento basandosi sull'intensità specifica di Ly$\alpha$ ($J_{Ly}$) derivata dalla densità del tasso di formazione stellare (SFRD).
  • Storia di Ionizzazione: Il codice traccia la ricombinazione dell'IGM bulk e la ionizzazione guidata da raggi X e fotoni UV. Modella il fattore di riempimento volumetrico delle regioni H II confrontando il tasso di emissione di fotoni ionizzanti con le ricombinazioni radiative.
  • Densità del Tasso di Formazione Stellare (SFRD): Una caratteristica innovativa di ECHO21 è l'inclusione di tre distinti modelli di SFRD:
    1. Fisicamente motivato: Basato sul collasso dei barioni nei aloni di materia oscura (utilizzando le funzioni di massa degli aloni di Press-Schechter, Sheth-Tormen o Tinker con un'efficienza di formazione stellare costante).
    2. Semi-empirico: Una variante del modello fisico che incorpora una scala temporale caratteristica della formazione stellare e un taglio netto a una temperatura viriale minima.
    3. Empiricamente motivato: Un modello vincolato dai dati della funzione di luminosità UV (forma di Madau-Dickinson a basso $z$, decadimento esponenziale ad alto $z$).

• Parametri: Il modello è controllato da sei parametri astrofisici liberi: normalizzazione di Ly$\alpha$ ($f_{Ly}$), indice spettrale della serie di Lyman ($s$), normalizzazione X ($f_X$), indice spettrale X ($w$), frazione di fuga dei fotoni ionizzanti ($f_{esc}$) e temperatura viriale minima ($min(T_{vir})$). Inoltre, il codice permette la variazione dei parametri cosmologici standard ($H_0, \Omega_m, \Omega_b, \sigma_8, n_s, T_{\gamma0}, Y_p$).

• Performance: L'approccio analitico consente a ECHO21 di generare un singolo caso di segnale globale a 21 cm in circa $O(1)$ secondo su una singola CPU. Il codice è progettato per la parallelizzazione MPI, consentendo la generazione di migliaia di modelli in pochi minuti.

Contributi Chiave

1. Inclusione del Riscaldamento da Ly$\alpha$: ECHO21 è il primo codice analitico a includere esplicitamente il riscaldamento da Ly$\alpha$ nell'evoluzione termica dell'IGM. Gli autori dimostrano che trascurare questo termine può introdurre bias significativi, particolarmente in scenari con bassa emissività X e forte emissione di Ly$\alpha$.
2. Diversità dei Modelli SFRD: Il codice offre unicamente tre diversi approcci per modellare la SFRD (fisico, semi-empirico ed empirico), permettendo agli utenti di valutare come la scelta di una prescrizione di formazione stellare impatti il segnale a 21 cm risultante.
3. Flessibilità Cosmologica: A differenza di molti codici specializzati, ECHO21 tratta i parametri cosmologici con la stessa facilità dei parametri astrofisici, facilitando lo studio delle incertezze cosmologiche (es. la tensione di Hubble) sul segnale a 21 cm.
4. Velocità e Scalabilità: Il codice è ottimizzato per un'inferenza rapida, capace di generare 15.625 modelli in uno spazio di parametri astrofisici a sei dimensioni in circa 17 minuti utilizzando 50 CPU.

Risultati

• Modello Fiduciale: Utilizzando una SFRD fisicamente motivata e una funzione di massa degli aloni di Press-Schechter, il modello fiduciale predice una storia di riionizzazione completata entro $z \approx 5.3$ con un redshift di riionizzazione del 50% di $z_{re} = 7.42$. L'opacità di scattering degli elettroni è calcolata come $\tau_e \approx 0.0592$, coerente con i vincoli di Planck 2018.
• Caratteristiche del Segnale: Il segnale globale a 21 cm fiduciale esibisce un primo picco di assorbimento di $-41.1$ mK a $z=87.5$ (Dark Ages) e un secondo picco di assorbimento più profondo di $-142.4$ mK a $z=19.2$ (Cosmic Dawn), seguito da un picco di emissione di $22.2$ mK a $z=11.2$.
• Impatto del Riscaldamento da Ly$\alpha$: I confronti con i modelli che escludono il riscaldamento da Ly$\alpha$ mostrano che l'omissione porta a un segnale di assorbimento a 21 cm più forte (più profondo). Gli autori sottolineano che il riscaldamento da Ly$\alpha$ è un componente cruciale per una modellazione accurata, poiché può contrastare gli effetti di raffreddamento e alterare la forza del segnale in modo non monotono a seconda dei parametri X.
• Impatto dei Modelli SFRD: La scelta del modello SFRD altera significativamente il segnale. I modelli empiricamente motivati, che predicono un inizio della formazione stellare più tardivo rispetto ai modelli fisicamente motivati, risultano in segnali ritardati di circa $\sim 8.8$ in redshift e più superficiali di $\sim 39.4$ mK. Allo stesso modo, la scelta della Funzione di Massa degli Aloni (HMF) influenza la tempistica e la profondità delle caratteristiche del segnale.
• Sensibilità ai Parametri: Il segnale è più sensibile a $f_{Ly}$, $f_X$, $w$ e $min(T_{vir})$. Le variazioni dei parametri cosmologici entro le attuali incertezze osservative risultano in variazioni relativamente piccole dell'ampiezza del segnale (ad esempio, un intervallo di $\sim 18$ mK nella profondità del picco di assorbimento).

Significatività e Rivendicazioni
Gli autori posizionano ECHO21 come uno strumento open-source flessibile ed estensibile, progettato per effettuare "inferenze astrofisiche rapide ma sufficientemente realistiche" dai futuri dataset a 21 cm (es. da REACH, SARAS o futuri esperimenti spaziali). Essi affermano che il codice è particolarmente adatto per:

• Eseguire robuste inferenze astrofisiche e cosmologiche tramite l'esplorazione di parametri su larga scala.
• Testare le pipeline di analisi dei dati (specificamente per la collaborazione REACH) con segnali fisicamente motivati piuttosto che con semplici fit empirici.
• Costruire sinergie tra le osservazioni a 21 cm e altre sonde cosmologiche, come le previsioni BBN, le anisotropie del CMB e la tensione di Hubble.

Il paper riconosce con modestia che, sebbene ECHO21 sia meno accurato dei modelli semi-numerici a causa delle sue approssimazioni analitiche (es. trascurare i molteplici scattering di Ly$\alpha$), la sua velocità e modularità lo rendono uno strumento ideale per le inferenze iniziali e per incorporare nuova fisica nelle versioni future. Il codice è pubblicamente disponibile tramite PyPI e GitHub.