A Consistent Implementation of Cluster Strong Lensing in Cosmological Simulation Light Cones

Sintesi Tecnica: Un'implementazione Coerente del Lensing Forte di Cluster in Light Cone di Simulazioni Cosmologiche

Problematica
Sebbene il lensing gravitazionale forte dei cluster di galassie sia un pilastro della cosmologia di precisione, le robuste previsioni teoriche sono storicamente rimaste indietro rispetto alle osservazioni di alta qualità. Questo divario deriva da due problemi primari: il disallineamento geometrico tra gli standard cubic box delle simulazioni e la geometria allungata richiesta per i light cone, e la tensione fondamentale tra il grande volume necessario per ospitare cluster massicci e l'alta risoluzione richiesta per risolvere le galassie sorgente di fondo e la sottostruttura.

I metodi di previsione attuali si affidano spesso ad approcci ibridi, estraendo singoli lenti dalle simulazioni e combinandoli con popolazioni di sorgenti analitiche o osservate. Questi metodi spesso non riescono a preservare le correlazioni attraverso il redshift o omettono la struttura correlata e non correlata lungo la linea di vista (LoS). Di conseguenza, vi è una mancanza di immagini di strong lensing simulate generate coerentemente dagli output delle simulazioni, dove la lente primaria, le sorgenti e tutta la materia interveniente siano estratti dallo stesso volume cosmologico in modo fedele al paradigma standard della materia oscura fredda (CDM).

Metodologia
Gli autori presentano una procedura completamente basata sulla simulazione per generare immagini di strong lensing direttamente dai dati delle particelle, utilizzando le simulazioni idrodinamiche cosmologiche IllustrisTNG. La metodologia consiste in tre componenti principali:

1. Remapping del Volume: Per affrontare i vincoli geometrici, gli autori impiegano una tecnica di remapping del box (utilizzando il codice BoxRemap) per trasformare il volume di simulazione cubico in un prisma rettangolare a base quadrata. Questa geometria presenta due assi corti (che rappresentano il piano del cielo) e un asse lungo (che rappresenta la linea di vista). Questa trasformazione preserva il volume totale e la struttura locale senza introdurre strutture ripetute all'interno di un singolo light cone, a patto che il campo visivo sia scelto in modo appropriato rispetto alla dimensione del box. Il remapping è applicato coerentemente a ogni snapshot della simulazione.
2. Estrazione dei Piani e Smoothing: Dai box rimappati, gli autori estraggono informazioni di massa e luce per creare una serie di piani lente e sorgente corrispondenti agli snapshot della simulazione.
   • Piani di Massa: Le particelle sono divise in gusci radiali basati sui redshift degli snapshot. Le mappe di massa sono generate effettuando lo smoothing delle particelle di materia oscura, stelle e gas utilizzando un kernel di vicinato più prossimo (32 vicini, kernel Wendland-C2). Per tenere conto della massa diffusa nelle regioni sparse ed evitare il rumore di shot, la massa totale del guscio viene aggiunta come un foglio uniforme, e la densità di materia media dell'universo viene sottratta per creare una mappa di contrasto di massa.
   • Piani Sorgente: Le particelle stellari sono trasformate in mappe di flusso per specifici filtri osservativi (ad esempio, JWST f200w). Il flusso è stimato in base all'età stellare, alla massa iniziale e alla metallicità, redshiftato al piano sorgente e convertito in magnitudine apparente. I subhalo con un numero insufficiente di particelle vengono sostituiti con profili di Sersic analitici per garantire la completezza.
3. Ray Tracing Multi-Piano: I piani estratti vengono utilizzati in un algoritmo di ray tracing multi-piano (utilizzando il codice GLAMER). La luce proveniente dai piani sorgente al redshift $z_s$ viene tracciata all'indietro attraverso tutti i piani lente ai redshift $z_l < z_s$. Le immagini risultanti vengono sommate per produrre una previsione teorica finale. La pipeline non richiede l'iniezione manuale di sorgenti o profili analitici; la cosmologia e il modello di formazione delle galassie della simulazione fungono da prior.

Contributi Chiave e Risultati
Utilizzando la simulazione TNG300-1, gli autori hanno generato light cone centrati su massicci gruppi di galassie ($M_{200m} \ge 5 \times 10^{14} M_\odot$) a un redshift della lente $z_G \approx 0.38$, estendendosi fino a un redshift sorgente di $z_{max} \approx 4.18$.

• Generazione di Immagini: La pipeline produce con successo caratteristiche di strong lensing, inclusi configurazioni di immagini multiple, archi e curve critiche, direttamente dai dati della simulazione senza iniezione ad hoc di sorgenti. Le immagini risultanti mostrano caratteristiche tipiche come archi tangenziali e morfologie complesse delle curve critiche guidate da merger di cluster e sottostrutture.
• Impatto della Struttura lungo la Linea di Vista (LoS): Lo studio quantifica l'impatto della struttura LoS non correlata confrontando il ray tracing completo del light cone con il ray tracing attraverso il solo piano della lente primaria.
  • Posizioni delle Immagini: La struttura LoS non correlata può spostare le posizioni relative delle immagini lenti di diversi arcsecondi (fino a $\sim 10''$ in casi estremi), superando significativamente i tipici residui di modellazione della massa delle lenti forti ($\lesssim 1''$).
  • Morfologia della Curva Critica: L'inclusione del materiale LoS introduce uno scatter di circa il $6\%$ nell'area della curva critica primaria di un cluster. In alcuni casi, l'area della curva critica primaria è stata aumentata del 30% o diminuita del 50%.
  • Area Critica Totale: Entro $100''$ dal minimo del potenziale del cluster, l'area critica totale a $z_s = 4$ cambia del $16^{+20}_{-14}\%$ quando viene inclusa la struttura LoS. Questo incremento deriva sia dalla creazione di nuove curve critiche da parte degli oggetti LoS, sia dalla modifica delle curve esistenti.

Significato e Limitazioni
L'articolo sostiene di fornire la prima dimostrazione di immagini di cluster strong lensing generate in modo auto-coerente da un light cone di una simulazione idrodinamica cosmologica. Questo approccio offre un metodo più automatizzato e "hands-off" per produrre grandi cataloghi di osservazioni mock, che possono essere utilizzati per studiare i sistemi nelle modellazioni delle lenti, i bias nell'inferenza della massa e la tensione della sezione d'urto del galaxy-galaxy strong lensing (GGSL).

Gli autori evidenziano diverse limitazioni:

• Potenza su Grande Scala: Il metodo non introduce modi di potenza su grande scala assenti nel box di simulazione iniziale, il che potrebbe influenzare il campionamento dello spettro di potenza su larga scala.
• Risoluzione: Il metodo è limitato dalla risoluzione di massa della simulazione, che fatica a risolvere simultaneamente cluster massicci e sottostrutture/sorgenti a bassa massa.
• Massa Diffusa: La massa diffusa è approssimata come un foglio uniforme piuttosto che preservarne la specifica struttura, sebbene ciò sia necessario per evitare il rumore di shot.
• Costo Computazionale: Generare un singolo light cone richiede l'elaborazione di ogni particella nel box di simulazione per ogni snapshot, il che è computazionalmente costoso e genera requisiti di archiviazione dati significativi.
• Shear Esterno: L'attuale implementazione non include lo shear esterno dalla struttura su larga scala al di fuori del light cone, sebbene l'effetto atteso sia stimato come piccolo ($\gamma_{ext} \sim 0.02$).

Gli autori concludono che incorporare l'informazione completa del light cone è cruciale per le previsioni teoriche e che è necessario tenere conto con cura della struttura LoS quando si interpretano le osservazioni, poiché la struttura non correlata può alterare significativamente le configurazioni delle immagini e le proprietà delle curve critiche.