Testing Dark Matter with Generative Models for Extragalactic Stellar Streams

1. Il Problema Scientifico

La natura della materia oscura (DM) rimane una delle questioni aperte più importanti della cosmologia. Mentre il modello $\Lambda$CDM predice profili di densità degli aloni di materia oscura con un "cuspo" centrale ripido (profilo NFW), modelli alternativi come la Materia Oscura Auto-Interagente (SIDM), la Materia Oscura Calda (WDM) e la Materia Oscura "Fuzzy" prevedono la formazione di "core" centrali a densità costante o pendenze interne piatte.

Attualmente, misurare questi profili di densità, specialmente nelle regioni interne e esterne degli aloni galattici, è estremamente difficile. Le tecniche tradizionali basate sulla cinetica stellare o sul gas HI sono limitate alle galassie vicine e spesso richiedono assunzioni sull'equilibrio dinamico. Inoltre, le osservazioni dei profili esterni degli aloni (oltre il raggio di scala) sono scarsamente esplorate, sebbene siano cruciali per comprendere la storia di accrescimento e le fusioni delle galassie.

Il lavoro si propone di sfruttare le correnti stellari extragalattiche (tidal streams) come traccianti gravitazionali indipendenti. La morfologia di queste correnti (lunghezza, curvatura, punti di svolta) è sensibile al potenziale gravitazionale dell'alone ospite. Tuttavia, l'inferenza dei parametri del potenziale dalle sole immagini 2D (senza dati cinematici 3D) è un problema altamente degenerato e computazionalmente costoso.

2. Metodologia: Il Modello Generativo "X-Stream"

Gli autori sviluppano un nuovo approccio generativo chiamato X-Stream, progettato per vincolare il profilo di densità radiale completo degli aloni di materia oscura partendo esclusivamente dalle immagini delle correnti stellari.

• Simulazione Accelerata: Il metodo si basa sul codice GPU-accelerato streamsculptor, implementato nel framework JAX. Questo permette di generare migliaia di realizzazioni di correnti stellari in parallelo. Vengono simulati flussi di particelle (metodo "particle-spray") in potenziali gravitazionali variabili.
• Parametri del Potenziale: Gli aloni sono modellati con un profilo di densità generalizzato NFW:
  $$\rho(r) = \frac{\rho_0}{(r/r_s)^\gamma (1 + r/r_s)^{\beta-\gamma}} \exp\left(-\left(\frac{r}{r_{cut}}\right)^2\right)$$
  I parametri liberi includono la massa dell'alone ($M_{halo}$), il raggio di scala ($r_s$), la pendenza interna ($\gamma$) e la pendenza esterna ($\beta$).
• Osservazioni Mock e Punti di Controllo: Per simulare dati osservativi reali (es. da Euclid, Rubin, Roman), le correnti simulate vengono convertite in istogrammi di densità. Vengono estratti "punti di controllo" (control points) dalle regioni più dense per rappresentare la morfologia osservata, ignorando la luminosità superficiale per focalizzarsi sulla geometria.
• Funzione di Costo e Campionamento:
  • Viene utilizzata la Divergenza di Kullback-Leibler (KL) come funzione di costo per misurare la discrepanza tra la distribuzione dei punti di controllo osservati e quella del modello generato.
  • Per affrontare la superficie di verosimiglianza estremamente stretta e piccata, viene impiegato un approccio di likelihood tempering (modificando la funzione di costo) combinato con il campionamento Nested Sampling (tramite il pacchetto nautilus).
  • L'uso di vmap in JAX permette di simulare 500 set di parametri simultaneamente su una singola GPU, riducendo i tempi di esecuzione da oltre 10 ore a poche ore.
• Inferenza Statistica: I campioni ottenuti vengono convertiti in regioni di confidenza (68% e 95%) utilizzando una statistica basata sul rapporto di densità gaussiana, validata per gestire le degenerazioni fisiche.

3. Contributi Chiave

1. Metodo Generativo Completo: X-Stream è il primo metodo in grado di inferire l'intero profilo di densità radiale (sia interno che esterno) di un alone di materia oscura sferico utilizzando solo la morfologia proiettata di correnti stellari extragalattiche.
2. Risoluzione delle Degenerazioni: Il metodo identifica e gestisce nuove degenerazioni fisiche, in particolare tra la massa dell'alone, le pendenze di densità ($\gamma, \beta$) e l'età dinamica della corrente.
3. Combinazione di Correnti Multiple: Dimostra che l'uso combinato di correnti situate a diverse distanze dal centro galattico (una interna, una esterna) permette di "cucire" insieme il profilo di densità completo, vincolando sia il cuspo interno che la pendenza esterna.
4. Test per Modelli Alternativi: Il metodo è stato validato per distinguere tra profili NFW (cuspidali) e profili con "core" (tipici della SIDM), anche in assenza di dati cinematici.

4. Risultati Principali

• Recupero dei Parametri: Applicando il metodo a due correnti simulate (una interna, $r \lesssim 3r_s$, e una esterna, $r \gtrsim 3r_s$) in un alone NFW, il sampler recupera con successo tutti i parametri veri (massa, velocità, età, pendenze) all'interno della regione di credibilità del 68%.
• Vincoli Radiali:
  • La corrente interna è altamente sensibile alla pendenza interna ($\gamma$), permettendo di vincolare fortemente il cuspo centrale.
  • La corrente esterna è più sensibile alla pendenza esterna ($\beta$) e alla massa totale, fornendo vincoli sulla regione esterna dell'alone.
  • L'incertezza frazionaria sulla densità è circa 5 volte più stretta per la corrente interna nelle regioni centrali rispetto a quella esterna.
• Rilevamento di "Core": In un test dove la corrente è evoluta in un potenziale con un "core" interno ($\gamma = 0.2$), il metodo recupera correttamente la pendenza e scarta il profilo NFW puro a un livello di confidenza di $2\sigma$.
• Degenerazioni Identificate:
  • Massa vs. Età: Una massa del progenitore più alta può essere compensata da un'età dinamica più bassa per produrre la stessa lunghezza di corrente.
  • Pendenza vs. Massa: Pendenze interne più ripide ($\gamma$ alto) richiedono masse di alone più basse per mantenere la stessa morfologia, mentre pendenze esterne più ripide ($\beta$ alto) richiedono masse più elevate.
  • Posizione lungo la linea di vista: Esiste una degenerazione bimodale tra la posizione del progenitore lungo la linea di vista e la direzione della velocità radiale (progenitore davanti o dietro la galassia ospite).

5. Significato e Prospettive Future

Questo lavoro rappresenta un passo fondamentale verso l'utilizzo delle future survey di basso contrasto superficiale (Euclid, Rubin Observatory, Roman, ARRAKIHS) per testare la fisica della materia oscura su larga scala.

• Scalabilità: Sebbene il metodo attuale richieda alcune ore di calcolo per galassia (limitato dall'integrazione orbitale completa), è parallelizzabile su migliaia di galassie.
• Test di Fisica Fondamentale: X-Stream offre un test indipendente per modelli di materia oscura come la SIDM, permettendo di cercare "core" in galassie extragalattiche senza la necessità di dati cinematici, che sono spesso indisponibili per oggetti distanti.
• Mappatura degli Aloni: Permette di mappare la struttura degli aloni di materia oscura dal raggio di scala fino al raggio viriale, una regione inaccessibile alle tecniche cinetiche tradizionali.

In sintesi, X-Stream trasforma le immagini di correnti stellari in potenti sonde gravitazionali, aprendo la strada a una mappatura dettagliata della distribuzione della materia oscura in migliaia di galassie a diversi redshift, fornendo vincoli stringenti sulla natura delle particelle di materia oscura.