JWST lensed quasar dark matter survey IV: Stringent warm dark matter constraints from the joint reconstruction of extended lensed arcs and quasar flux ratios

Sintesi Tecnica: JWST Lensing Survey IV – Vincoli sulla Materia Oscura Calda tramite Ricostruzione Congiunta

Problematica
La natura della materia oscura (DM) rimane ignota, poiché il paradigma della Materia Oscura Fredda (CDM) prevede una gerarchia di aloni priva di scala fino a masse molto piccole. Teorie alternative, come la Materia Oscura Calda (WDM), prevedono una soppressione della struttura al di sotto di una specifica scala di free-streaming, caratterizzata dalla massa di half-mode ($m_{hm}$). Il lensing gravitazionale forte di quasar con immagini quadruplamente sdoppiate offre una sonda unica di questa sottostruttura, misurando le perturbazioni nelle posizioni delle immagini e nei rapporti di flusso. Tuttavia, le analisi precedenti, basate esclusivamente sulle posizioni di sorgenti puntiformi e sui rapporti di flusso, sono state limitate dalle degenerazioni tra il profilo di massa del deflettore principale (macromodello) e la popolazione di sottostrutture. Le incertezze nel macromodello (ad esempio, ellitticità, shear esterno e struttura angolare) possono simulare o mascherare le firme degli aloni di materia oscura, portando a vincoli più deboli sulle proprietà della DM. Inoltre, il microlensing stellare ha storicamente contaminato le misurazioni dei rapporti di flusso per sorgenti compatte, sebbene ciò sia stato mitigato dalle recenti osservazioni JWST delle regioni di polvere calda estese.

Metodologia
Questo articolo presenta un framework di forward-modeling che ricostruisce simultaneamente archi lenti estesi, posizioni delle immagini e rapporti di flusso per un campione di 28 quasar con immagini quadruplamente sdoppiate. L'analisi utilizza un framework di inferenza bayesiana per vincolare gli iper-parametri che descrivono la popolazione di sottostrutture di DM.

• Integrazione dei Dati: Lo studio combina i rapporti di flusso MIRI di JWST (dalla regione della polvere calda del quasar, evitando il microlensing stellare) con dati astrometrici e immagini di archi lenti estesi da JWST (NIRCam e MIRI) e HST.
• Modellazione Congiunta: A differenza dei lavori precedenti che trattavano separatamente gli archi lenti e i rapporti di flusso, questa analisi impiega pesi di campionamento per importanza ($w_{img}$) per propagare i vincoli degli archi estesi sui parametri del macromodello. Questo rompe la covarianza tra il macromodello e la popolazione di sottostrutture, stringendo i vincoli sui rapporti di flusso previsti dal modello.
• Modellazione Fisica:
  • Sottostruttura: Il framework genera popolazioni complete di subaloni, aloni lungo la linea di vista e ammassi globulari utilizzando il software pyHalo. Incorpora un modello di evoluzione mareale (basato su Du et al.) che predice i profili di densità dei subaloni in base alle proprietà di infall, tenendo conto della perdita di massa e dei cambiamenti di concentrazione.
  • Effetti WDM: Il modello modifica la funzione di massa degli aloni CDM e la relazione concentrazione-massa per tenere conto della soppressione da free-streaming in WDM, parametrizzata dalla massa di half-mode $m_{hm}$.
  • Macromodello: Il deflettore principale è modellato utilizzando un profilo a legge di potenza ellittica (EPL) aumentato da perturbazioni multipolari ellittiche (m=1, 3, 4) per tenere conto della struttura angolare (ad esempio, boxiness/diskiness) che altrimenti potrebbe mascherarsi da sottostruttura.
• Inferenza: La funzione di verosimiglianza viene valutata generando milioni di realizzazioni (oltre 174 milioni per l'intero campione) di modelli di lente. L'analisi marginalizza sui parametri di disturbo (macromodello, luce della sorgente, realizzazioni della sottostruttura) per inferire le distribuzioni a posteriori degli iper-parametri, specificamente la massa di half-mode $m_{hm}$ e l'ampiezza della funzione di massa dei subaloni $\Sigma_{sub}$.

Contributi Chiave

1. Rottura delle Degenerazioni: Il principale avanzamento metodologico è la modellazione simultanea di archi lenti estesi insieme a rapporti di flusso e posizioni delle immagini. Questo approccio riduce significativamente l'incertezza nel macromodello, rompendo così la covarianza tra l'abbondanza di subaloni ($\Sigma_{sub}$) e la scala di free-streaming ($m_{hm}$).
2. Migliore Modellazione dei Subaloni: Il documento implementa un modello di evoluzione mareale più rigoroso dal punto di vista fisico che predice le proprietà dei subaloni al momento dell'infall, piuttosto che alla redshift della lente. Ciò tiene conto della dipendenza dello stripping mareale dalla concentrazione di infall, il che è particolarmente rilevante per la WDM dove gli aloni sono meno concentrati.
3. Fonti di Perturbazione Complete: L'analisi include esplicitamente ammassi globulari e perturbazioni multipolari ellittiche come fonti di rumore nei rapporti di flusso, garantendo che i vincoli sulla WDM non siano distorti da strutture a piccola scala non modellate.
4. Ampiezza del Campione: Questo lavoro analizza il più grande campione di quasar con immagini quadruplamente sdoppiate ad oggi per l'inferenza della DM (28 sistemi), sfruttando tutte le capacità del survey di quasar lenti di JWST.

Risultati

• Vincoli sulla Materia Oscura Calda (WDM):
  • Assumendo un prior log-uniformo sull'ampiezza della funzione di massa dei subaloni, l'inclusione degli archi lenti migliora il limite sulla massa di half-mode di 0,6 dex rispetto all'uso dei soli rapporti di flusso e posizioni delle immagini.
  • Quando vengono applicati prior informativi su $\Sigma_{sub}$ derivati dal modello semi-analitico galacticus, l'analisi inferisce $m_{hm} < 10^{7.4} M_\odot$ (fattore di Bayes 10:1).
  • Applicando i prior basati su simulazioni N-body (Symphony suite), il vincolo si stringe a $m_{hm} < 10^{7.2} M_\odot$.
  • Questi limiti corrispondono a limiti inferiori sulla massa di una particella relitto termico WDM di 6.5 keV e 7.4 keV, rispettivamente.
  • I limiti di esclusione al 95% sono $m_{hm} < 10^{6.8} M_\odot$ e $m_{hm} < 10^{6.7} M_\odot$, corrispondenti a masse di particelle $> 9.6$ keV e $> 10.2$ keV.
• Abbondanza di Subaloni CDM:
  • Assumendo il limite CDM ($m_{hm} < 10^5 M_\odot$), lo studio misura la densità di massa superficiale proiettata dei subaloni pari a $1.7^{+2.6}_{-1.2} \times 10^7 M_\odot \text{kpc}^{-2}$ con una confidenza del 95%.
  • Ciò corrisponde a una frazione di massa proiettata del $3.2^{+5.1}_{-2.3}\%$.
  • Questa misurazione è statisticamente coerente con le previsioni di galacticus ma mostra una lieve tensione ($\sim 1.8\sigma$) con recenti simulazioni N-body, che predicono un'abbondanza di subaloni inferiore.

Significatività
L'articolo sostiene che questi risultati rappresentano alcuni dei vincoli più forti sulla WDM ad oggi e la misurazione più precisa dell'abbondanza di subaloni attorno alle lenti forti. La significatività risiede nella dimostrazione che la modellazione congiunta di archi lenti estesi e rapporti di flusso isola efficacemente il segnale della sottostruttura di materia oscura dalle incertezze del profilo di massa del deflettore principale. Gli autori osservano che i loro vincoli superano i precedenti limiti del lensing forte e sono competitivi o più forti di altre sonde a piccola scala come la foresta Lyman-$\alpha$ e i conteggi di galassie nane. Questo lavoro stabilisce un framework robusto per future inferenze sulla DM utilizzando i grandi campioni di lenti previsti dai prossimi survey (Euclid, Rubin, Roman), evidenziando il potenziale trasformativo della combinazione di imaging ad alta risoluzione dallo spazio con l'avanzata modellazione forward.