Detection of a dark matter subhalo in the strongly lensed system PJ011646

Sintesi Tecnica: Rilevamento di un Sottogalassia di Materia Oscura nel Sistema Fortemente Lenticolare PJ011646

Enunciato del Problema
La natura della materia oscura rimane un problema fondamentale irrisolto in cosmologia. Mentre il modello di Materia Oscura Fredda (CDM) prevede una formazione gerarchica di aloni auto-legati che si estendono fino a masse molto basse, i modelli di Materia Oscura Calda (WDM) prevedono una soppressione degli aloni a bassa massa al di sotto di una scala di taglio caratteristica. Distinguere tra questi modelli richiede misure statistiche della popolazione di sottogalassie, in particolare nel regime di massa in cui le previsioni divergono ($M_{200} \sim 10^7 - 10^9 M_\odot$). Una sfida primaria nel rilevare queste sottostrutture tramite lente gravitazionale forte è la "degenerazione" tra i segnali genuini di sottogalassie e la complessità angolare non modellata nella distribuzione di massa liscia della galassia lente. I modelli semplificati a legge di potenza ellittica (PL) spesso non riescono a catturare le deviazioni nelle galassie di tipo precoce (ETG), portando a rilevamenti spurii di sottogalassie o a mascherarne di reali. Inoltre, i rilevamenti precedenti sono stati limitati dalla risoluzione angolare dei telescopi ottici (ad es. HST) e dalla contaminazione della luce della galassia lente.

Metodologia
Gli autori presentano un'analisi di lente gravitazionale forte del sistema PJ011646 (lente $z=0.555$, sorgente $z=2.125$) utilizzando osservazioni ad alta risoluzione ($\sim 0.1$ arcsec) del continuo in polvere dall'Atacama Large Millimeter/submillimeter Array (ALMA) nella Banda 7. L'analisi utilizza la suite di software open-source PyAutoLens, PyAutoFit e PyAutoGalaxy.

La metodologia procede in tre fasi principali:

1. Modellizzazione Macroscopica: La distribuzione di massa della lente è modellata utilizzando un profilo a legge di potenza ellittica (EPL) con shear esterno. Per affrontare le limitazioni dei modelli ellittici semplici, gli autori incorporano multipoli di Fourier di ordine superiore (ordini $m=3$ e $m=4$, e inoltre $m=1$) per catturare le deviazioni angolari nella densità di massa. La convergenza di questi multipoli è parametrizzata per evitare degenerazioni associate alle fasi periodiche.
2. Ricostruzione della Sorgente: La luminosità superficiale della sorgente è ricostruita utilizzando un metodo di inversione bayesiana pixelizzata (Suyu et al. 2006) con pixelizzazione a mesh di Voronoi e interpolazione del vicino naturale. Ciò viene eseguito direttamente nel piano $uv$ per sfruttare i dati interferometrici completi.
3. Ricerca di Sottogalassie: Viene eseguita una ricerca basata su griglia nel piano immagine per rilevare perturbazioni. Un profilo sferico Navarro-Frenk-White (NFW) viene aggiunto come perturbatore in ogni cella della griglia, e la differenza di log-verosimiglianza bayesiana ($\Delta \ln Z$) viene calcolata rispetto al miglior modello macroscopico di adattamento. Viene adottata una soglia di $\Delta \ln Z > 10$ per identificare candidati robusti. Gli autori testano inoltre un profilo NFW generalizzato (gNFW) con pendenza interna libera ($\gamma$) per verificare deviazioni dal profilo CDM standard.

Risultati Chiave

• Miglioramento della Modellizzazione Macroscopica: L'inclusione di multipoli di terzo e quarto ordine ($m=3, 4$) migliora significativamente l'adattamento ai dati, producendo un aumento della verosimiglianza bayesiana di $\Delta \ln Z = 52.1$ rispetto a un modello PL puro. Le ampiezze di questi multipoli sono $\simeq 1.5\%$ della convergenza, coerenti con le osservazioni di ETG vicine. Un modello che include un multipolo di primo ordine ($m=1$) fornisce un miglioramento aggiuntivo marginale ($\Delta \ln Z = 6.5$), suggerendo una possibile compagna fisica ma non guidando il segnale primario.
• Rilevamento di Sottogalassie: Utilizzando il miglior modello macroscopico PL+MP di adattamento, la ricerca su griglia identifica due regioni nel piano immagine con $\Delta \ln Z > 10$. Tuttavia, un'analisi di follow-up rivela che queste due regioni mappano alla stessa posizione nel piano sorgente, indicando una degenerazione. Un modello a due sottogalassie non migliora la verosimiglianza rispetto a un modello a singola sottogalassia, e la distribuzione a posteriori per il secondo candidato produce concentrazioni non fisiche ($c_{200} > 1000$). Di conseguenza, gli autori concludono che solo una sottogalassia è fisicamente reale, situata a destra dell'arco settentrionale.
• Proprietà della Sottogalassia: La sottogalassia rilevata ha una massa $M_{200} = 2.78^{+0.43}_{-0.66} \times 10^{10} M_\odot$ e una concentrazione $c_{200} = 30^{+5}_{-7}$. Questa concentrazione è coerente con le aspettative per una sottogalassia NFW privata di massa all'interno di un alone genitore. La significatività del rilevamento è di $\sim 5.8\sigma$ rispetto al modello macroscopico PL+MP.
• Vincoli di Massa Robusti: A causa della degenerazione tra il raggio di scala NFW e la normalizzazione, la massa totale $M_{200}$ non è vincolata strettamente. Invece, gli autori definiscono un "raggio robusto" ($R_{sub} \approx 2$ kpc) dove l'incertezza frazionaria sulla massa racchiusa proiettata è minimizzata. All'interno di questo raggio, la massa racchiusa è $M_{sub} = 3.57^{+0.16}_{-0.14} \times 10^9 M_\odot$.
• Limiti di Sensibilità: Nelle celle della griglia in cui non viene rilevata alcuna sottogalassia ($\Delta \ln Z < 10$), gli autori derivano limiti superiori sulle masse di sottogalassie rilevabili. Nelle regioni più sensibili (lungo gli archi lenticolari luminosi), la massa minima rilevabile a una significatività di $3\sigma$ è $M_{200} \approx 8 \times 10^8 M_\odot$, assumendo la relazione massa-concentrazione di Ludlow et al. (2016).

Significato e Affermazioni
Il lavoro afferma che questa ricerca dimostra la capacità delle immagini in continuo di ALMA a risoluzione sub-arcosecondo di sondare la sottostruttura di materia oscura in un regime di massa ($M_{200} \sim 10^9 M_\odot$) in cui i modelli CDM e WDM prevedono abbondanze diverse. Modellando esplicitamente la complessità angolare della galassia lente tramite multipoli, gli autori mitigano le incertezze sistemiche che hanno storicamente afflitto le ricerche di sottogalassie. Il rilevamento di una sottogalassia con proprietà coerenti con le previsioni CDM, insieme alla derivazione di limiti superiori stringenti nelle regioni di non rilevamento, fornisce un test robusto della formazione di strutture su piccola scala. Gli autori notano che, sebbene la risoluzione angolare sia paragonabile a quella di HST, l'alto ingrandimento del sistema e l'assenza di contaminazione da luce di primo piano nei dati ALMA permettono di sondare masse inferiori rispetto a quanto precedentemente realizzabile con dati ottici. Lo studio si posiziona come un passo verso ricerche sistematiche e statisticamente significative di sottostrutture di materia oscura utilizzando il crescente campione di galassie a formazione stellare polverose fortemente lenticolari (DSFG) osservate con ALMA.