How X-rays heat the IGM in different 21-cm simulation codes: a comparison between Licorice and Beorn

Immaginate l'universo come un gigantesco oceano invisibile di gas, composto principalmente da idrogeno, che si estende in ogni direzione. Per i primi centinaia di milioni di anni dopo il Big Bang, questo oceano era buio, freddo e silenzioso. Poi, le prime stelle e galassie si sono accese come fari nella nebbia. La loro intensa radiazione ha iniziato a riscaldare il gas e a strappare gli elettroni dagli atomi, un processo chiamato "reionizzazione". Gli astronomi sono disperati nel voler vedere questo periodo svelarsi perché contiene i segreti di come il nostro universo sia cresciuto. Per farlo, stanno costruendo un enorme radiotelescopio chiamato Square Kilometre Array (SKA), che agisce come una gigantesca orecchia in ascolto di un debole sussurro: il "segnale a 21 cm". Questo segnale è un'onda radio specifica emessa dall'idrogeno neutro, e la sua intensità cambia a seconda di quanto è caldo il gas. Se il gas è freddo, il segnale appare diverso rispetto a quando è caldo. Misurando questo segnale, gli scienziati sperano di mappare la storia del riscaldamento dell'universo.

Tuttavia, c'è un problema. L'universo è troppo vasto e complesso per essere osservato direttamente in tempo reale; dobbiamo costruire modelli informatici per simulare ciò che è accaduto miliardi di anni fa. Queste simulazioni sono come macchine del tempo digitali. Ma proprio come ogni macchina del tempo, esse sono affidabili solo quanto le regole che seguono. Alcune simulazioni sono incredibilmente dettagliate, tracciando ogni singolo particella di gas e ogni raggio di luce, ma richiedono settimane di lavoro ai supercomputer. Altre sono più veloci e utilizzano scorciatoie, rendendole più facili da usare ma potenzialmente meno accurate. La grande domanda è: queste scorciatoie funzionano abbastanza bene da rendere affidabili i risultati? Se due diversi codici di simulazione ci raccontano storie differenti su quanto sia diventato caldo il gas, quale dovremmo credere? È l'enigma che un nuovo studio di Romain Merot e del suo team si propone di risolvere.

Il team ha deciso di mettere due "macchine del tempo" molto diverse testa a testa per vedere come gestiscono il riscaldamento del gas cosmico tramite i raggi X. Da una parte, avevano Licorice, una simulazione 3D ad alta definizione che agisce come un videogioco a tutto movimento, calcolando come il gas si muove e come la luce viaggia attraverso di esso in tre dimensioni. Dall'altra parte, avevano Beorn, un codice più veloce e semplice che utilizza un approccio "unidimensionale". Pensate a Licorice come a uno chef maestro che cucina uno stufato complesso, mescolando ogni singolo ingrediente, mentre Beorn è come un microonde che scalda l'intero pasto basandosi su una ricetta media. I ricercatori volevano vedere se il microonde potesse produrre un pasto che avesse lo stesso sapore dello stufato dello chef, specificamente riguardo a come i raggi X hanno riscaldato l'idrogeno del gas.

Per rendere la sfida equa, gli scienziati non si sono limitati a far girare i due codici con le proprie impostazioni casuali. Invezione, hanno preso gli "ingredienti" (le stelle e la loro emissione di luce) dalla dettagliata simulazione Licorice e li hanno inseriti in Beorn. Ciò ha garantito che entrambi i codici cercassero di cucinare esattamente lo stesso universo, pur usando metodi diversi. Hanno poi confrontato i risultati: la temperatura del gas, la luminosità del segnale a 21 cm e i modelli statistici di come quel segnale fluttua attraverso il cielo.

I risultati sono stati un misto di buone notizie e un avvertimento. Quando hanno guardato il quadro generale — la temperatura media del gas e la forza complessiva del segnale — i due codici concordavano molto bene. Era come se il microonde e lo chef fossero riusciti entrambi a portare lo stufato alla giusta temperatura media. Tuttavia, quando i ricercatori hanno guardato più da vicino i dettagli, le differenze hanno iniziato a emergere. La distribuzione del calore non era la stessa. Nella simulazione 3D dettagliata, il calore era a chiazze, con alcuni punti molto più caldi e altri più freddi, creando una texture complaesa. Nella simulazione 1D più veloce, il calore era più fluido e uniforme.

Queste piccole differenze nella "texture" del calore si sono rivelate sorprendentemente importanti. Quando il team ha calcolato lo "spettro di potenza" — un modo elaborato per misurare quanto il segnale fluttua a diverse scale — ha trovato una discrepanza di circa il 30% tra i due codici. Per dare un termine di paragone, se cercaste di indovinare la ricetta di una torta basandovi su un singolo boccone, una differenza del 30% nel sapore significherebbe che indovinereste completamente gli ingredienti.

I ricercatori hanno poi scavato più a fondo per scoprire perché i codici discordassero. Hanno scoperto che le scorciatoie utilizzate dal codice più veloce erano la causa del problema. Nello specifico, il codice 1D assumeva che la densità del gas fosse la stessa ovunque (come assumere che l'oceano sia perfettamente piatto) e utilizzava una formula semplificata per quanto calore i raggi X effettivamente depositassero. In realtà, il gas è grumoso, e il riscaldamento dipende esattamente da dove si trova il gas. Quando il team ha costretto il codice 3D dettagliato a utilizzare queste stesse scorciatoie semplificate, i risultati si sono avvicinati al codice 1D, confermando che queste approssimazioni erano la fonte dell'errore.

Infine, il team si è posto la domanda suprema: "Se usiamo queste diverse simulazioni per analizzare i dati del futuro telescopio SKA, otterremo le risposte sbagliate?". Hanno simulato uno scenario in cui il telescopio avrebbe osservato per 100 ore. Hanno scoperto che le differenze tra i codici avrebbero portato a un bias significativo nelle proprietà stimate delle prime stelle. In termini statistici, l'errore era abbastanza grande da essere maggiore di 1 sigma (una misura di confidenza), e per alcuni parametri, era addirittura superiore. Ciò significa che se gli astronomi avessero usato il codice più veloce e semplice per interpretare i dati reali, avrebbero potuto concludere con sicurezza che le prime stelle fossero diverse da quelle reali.

Lo studio conclude che, sebbene i codici più veloci siano utili e riescano a comprendere la storia "media", non sono ancora abbastanza precisi per l'era dello SKA. La differenza del 30% nello spettro di potenza e il conseguente bias nei dati suggeriscono che non possiamo semplicemente fidarci delle scorciatoie quando cerchiamo di misurare l'universo con estrema precisione. Gli autori suggeriscono che, sebbene i codici siano abbastanza simili per gli attuali telescopi meno sensibili, è necessario molto lavoro per perfezionare i codici più veloci o sviluppare nuovi modi per correggere i loro errori prima che lo SKA inizi ad ascoltare. È un promemoria del fatto che, nella ricerca di una comprensione dell'alba dell'universo, il diavolo sta davvero nei dettagli.