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

Résumé Technique : Une Implémentation Cohérente du Lentillage Fort de Clusters dans les Cônes de Lumière des Simulations Cosmologiques

Énoncé du Problème
Bien que le lentillage gravitationnel fort par les amas de galaxies soit une pierre angulaire de la cosmologie de précision, les prédictions théoriques robustes ont historiquement accusé un retard par rapport aux observations de haute qualité. Cet écart provient de deux problèmes principaux : l'incompatibilité géométrique entre les boîtes de simulation cubiques standards et la géométrie allongée requise pour les cônes de lumière, et la tension fondamentale entre le grand volume nécessaire pour héberger des amas massifs et la haute résolution requise pour résoudre les galaxies sources en arrière-plan et la substructure.

Les méthodes de prévision actuelles reposent souvent sur des approches hybrides, extrayant des lentilles individuelles de simulations et les combinant avec des populations de sources analytiques ou observées. Ces méthodes échouent fréquemment à préserver les corrélations à travers le décalage vers le rouge (redshift) ou omettent les structures corrélées et non corrélées le long de la ligne de visée (LoS). Par conséquent, il existe un manque d'images de lentillage fort simulées, générées de manière cohérente à partir des données de simulation, où la lentille primaire, les sources et toute la matière intermédiaire sont tirées du même volume cosmologique d'une manière fidèle au paradigme standard de la matière noire froide.

Méthodologie
Les auteurs présentent une procédure entièrement basée sur la simulation pour générer des images de lentillage fort directement à partir de données de particules, en utilisant les simulations hydrodynamiques cosmologiques IllustrisTNG. La méthodologie se compose de trois composantes principales :

1. Remappage de Volume : Pour répondre aux contraintes géométriques, les auteurs emploient une technique de remappage de boîte (utilisant le code BoxRemap) pour transformer le volume de simulation cubique en un prisme rectangulaire à base carrée. Cette géométrie présente deux axes courts (représentant le plan du ciel) et un axe long (représentant la ligne de visée). Cette transformation préserve le volume total et la structure locale sans introduire de structures répétées au sein d'un seul cône de lumière, à condition que le champ de vision soit choisi de manière appropriée par rapport à la taille de la boîte. Le remappage est appliqué de manière cohérente à chaque instantané (snapshot) de la simulation.
2. Extraction de Plans et Lissage : À partir des boîtes remappées, les auteurs extraient les informations de masse et de lumière pour créer une série de plans de lentilles et de sources correspondant aux instantanés de la simulation.
   • Plans de Masse : Les particules sont divisées en couches radiales basées sur les redshifts des instantanés. Des cartes de masse sont générées en lissant les particules de matière noire, d'étoiles et de gaz à l'aide d'un noyau de plus proches voisins (32 voisins, noyau Wendland-C2). Pour tenir compte de la masse diffuse dans les régions éparses et éviter le bruit de tir (shot noise), la masse totale de la couche est ajoutée sous forme de feuille uniforme, et la densité de matière moyenne de l'univers est soustraite pour créer une carte de contraste de masse.
   • Plans de Sources : Les particules stellaires sont lissées en cartes de flux pour des filtres observationnels spécifiques (par exemple, JWST f200w). Le flux est estimé en fonction de l'âge stellaire, de la masse initiale et de la métallicité, décalé vers le rouge vers le plan de la source, et converti en magnitude apparente. Les sous-halos ayant un nombre de particules insuffisant sont remplacés par des profils de Sersic analytiques pour assurer la complétude.
3. Ray Tracing Multi-Plans : Ces plans extraits sont utilisés dans un algorithme de lancer de rayons multi-plans (utilisant le code GLAMER). La lumière provenant des plans de sources à un redshift $z_s$ est tracée à rebours à travers tous les plans de lentilles à des redshifts $z_l < z_s$. Les images résultantes sont cumulées pour produire une prédiction théorique finale. Le pipeline ne nécessite aucune injection manuelle de sources ou de profils analytiques ; la cosmologie et le modèle de formation des galaxies de la simulation servent de priorité (prior).

Contributions Clés et Résultats
En utilisant la simulation TNG300-1, les auteurs ont généré des cônes de lumière centrés sur des groupes de galaxies massifs ($M_{200m} \ge 5 \times 10^{14} M_\odot$) à un redshift de lentille $z_G \approx 0.38$, s'étendant jusqu'à un redshift de source $z_{max} \approx 4.18$.

• Génération d'Images : Le pipeline produit avec succès des caractéristiques de lentillage fort, incluant des configurations d'images multiples, des arcs et des courbes critiques, directement à partir des données de simulation sans injection de source ad hoc. Les images résultantes montrent des caractéristiques typiques telles que des arcs tangentiels et des morphologies de courbes critiques complexes pilotées par les fusions d'amas et la substructure.
• Impact de la Structure de la Ligne de Visée (LoS) : L'étude quantifie l'impact de la structure non corrélée de la LoS en comparant le lancer de rayons complet du cône de lumière contre le lancer de rayons à travers le plan de la lentille primaire uniquement.
  • Positions des Images : La structure non corrélée de la LoS peut déplacer les positions relatives des images lentillées de plusieurs secondes d'arc (jusqu'à $\sim 10''$ dans les cas extrêmes), dépassant significativement les résidus typiques de modélisation de masse des lentilles fortes ($\lesssim 1''$).
  • Morphologie de la Courbe Critique : L'inclusion de la matière de la LoS introduit une dispersion de $\sim 6\%$ dans l'aire de la courbe critique primaire d'un amas. Dans certains cas, l'aire de la courbe critique primaire a été augmentée de 30 % ou diminuée de 50 %.
  • Aire Critique Totale : À l'intérieur de $100''$ du minimum du potentiel de l'amas, l'aire critique totale à $z_s = 4$ change de $16^{+20}_{-14}\%$ lorsque la structure de la LoS est incluse. Cet accroissement provient à la fois de la création de nouvelles courbes critiques par les objets de la LoS et de la modification des courbes existantes.

Signification et Limitations
L'article affirme fournir la première démonstration d'images de lentillage fort de clusters générées de manière auto-cohérente à partir d'un cône de lumière de simulation hydrodynamique cosmologique. Cette approche offre une méthode plus automatisée et "sans intervention" pour produire de vastes catalogues d'observations simulées, qui peuvent être utilisés pour étudier les biais systématiques dans la modélisation de lentilles, les biais d'inférence de masse et la tension de la section efficace de lentillage fort galaxie-galaxie (GGSL).

Les auteurs notent plusieurs limitations :

• Puissance à Grande Échelle : La méthode n'introduit pas de modes de puissance à grande échelle absents de la boîte de simulation initiale, ce qui pourrait biaiser l'échantillonnage du spectre de puissance à grande échelle.
• Résolution : La méthode est limitée par la résolution de masse de la simulation, qui peine à résoudre simultanément les amas massifs et la substructure/sources de faible masse.
• Masse Diffuse : La masse diffuse est approximée comme une feuille uniforme plutôt que de préserver sa structure spécifique, bien que cela soit nécessaire pour éviter le bruit de tir.
• Coût Computationnel : Générer un seul cône de lumière nécessite de traiter chaque particule de la boîte de simulation à chaque instantané, ce qui est coûteux en calcul et génère des exigences de stockage de données importantes.
• Cisaillement Externe : L'implémentation actuelle n'inclut pas le cisaillement externe provenant de la structure à grande échelle en dehors du cône de lumière, bien que l'effet attendu soit estimé comme étant faible ($\gamma_{ext} \sim 0.02$).

Les auteurs concluent qu'incorporer l'information complète du cône de lumière est crucial pour les prédictions théoriques et que la prise en compte de la structure de la LoS est nécessaire lors de l'interprétation des observations, car la structure non corrélée peut modifier considérablement les configurations d'images et les propriétés des courbes critiques.