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

Resumen Técnico: Una Implementación Consistente de Lente Fuerte de Cúmulos en Conos de Luz de Simulaciones Cosmológicas

Planteamiento del Problema
Si bien el lente gravitacional fuerte de cúmulos de galaxias es un pilar de la cosmología de precisión, las predicciones teóricas robustas han quedado históricamente rezagadas respecto a las observaciones de alta calidad. Esta brecha surge de dos problemas principales: la discordancia geométrica entre las cajas cúbicas estándar de las simulaciones y la geometría alargada requerida para los conos de luz, y la tensión fundamental entre el gran volumen necesario para albergar cúmulos masivos y la alta resolución requerida para resolver galaxias fuentes de fondo y la subestructura.

Los métodos de pronóstico actuales suelen depender de enfoques híbridos, extrayendo lentes individuales de las simulaciones y combinándolos con poblaciones de fuentes analíticas u observadas. Estos métodos frecuentemente fallan al no preservar las correlaciones a través del desplazamiento al rojo (redshift) o al omitir la estructura correlacionada y no correlacionada de la línea de visión (LoS, por sus siglas en inglés). En consecuencia, existe una falta de imágenes de lente fuerte simuladas generadas de manera consistente a partir de los datos de la simulación, donde el lente primario, las fuentes y toda la materia interviniente se extraen del mismo volumen cosmológico de una manera fiel al paradigma estándar de la materia oscura fría.

Metodología
Los autores presentan un procedimiento basado totalmente en simulaciones para generar imágenes de lente fuerte directamente a partir de datos de partículas, utilizando las simulaciones hidrodinámicas cosmológicas IllustrisTNG. La metodología consta de tres componentes principales:

1. Remapeo de Volumen: Para abordar las restricciones geométricas, los autores emplean una técnica de remapeo de caja (utilizando el código BoxRemap) para transformar el volumen cúbico de la simulación en un prisma rectangular de base cuadrada. Esta geometría presenta dos ejes cortos (que representan el plano del cielo) y un eje largo (que representa la línea de visión). Esta transformación preserva el volumen total y la estructura local sin introducir estructuras repetidas dentro de un solo cono de luz, siempre que el campo de visión se elija adecuadamente en relación con el tamaño de la caja. El remapeo se aplica consistentemente a cada instantánea (snapshot) de la simulación.
2. Extracción de Planos y Suavizado: De las cajas remapeadas, los autores extraen información de masa y luz para crear una serie de planos de lente y de fuente correspondientes a las instantáneas de la simulación.
   • Planos de Masa: Las partículas se dividen en capas radiales basadas en los desplazamientos al rojo de las instantáneas. Los mapas de masa se generan suavizando las partículas de materia oscura, estrellas y gas mediante un núcleo de vecino más cercano (32 vecinos, núcleo Wendland-C2). Para contabilizar la masa difusa en regiones dispersas y evitar el ruido de disparo (shot noise), se añade la masa total de la capa como una lámina uniforme, y se resta la densidad de materia media del universo para crear un mapa de contraste de masa.
   • Planos de Fuente: Las partículas estelares se suavizan en mapas de flujo para filtros observacionales específicos (por ejemplo, JWST f200w). El flujo se estima basándose en la edad estelar, la masa inicial y la metalicidad, desplazado al rojo hacia el plano de la fuente, y convertido a magnitud aparente. Los subhalos con recuentos de partículas insuficientes se reemplazan con perfiles analíticos de Sersic para asegurar la completitud.
3. Trazado de Rayos Multi-plano: Los planos extraídos se utilizan en un algoritmo de trazado de rayos multi-plano (utilizando el código GLAMER). La luz de los planos de fuente en un desplazamiento al rojo $z_s$ se traza hacia atrás a través de todos los planos de lente en desplazamientos al rojo $z_l < z_s$. Las imágenes resultantes se apilan para producir una predicción teórica final. El proceso requiere la no inyección manual de fuentes o perfiles analíticos; el modelo de formación de galaxias y la cosmología de la simulación sirven como el prior.

Contribuciones Clave y Resultados
Utilizando la simulación TNG300-1, los autores generaron conos de luz centrados en grupos galácticos masivos ($M_{200m} \ge 5 \times 10^{14} M_\odot$) en un desplazamiento al rojo del lente de $z_G \approx 0.38$, extendiéndose hasta un desplazamiento al rojo de la fuente de $z_{max} \approx 4.18$.

• Generación de Imágenes: El proceso produce con éxito características de lente fuerte, incluyendo configuraciones de imágenes múltiples, arcos y curvas críticas, directamente desde los datos de la simulación sin inyección de fuentes ad-hoc. Las imágenes resultantes muestran características típicas como arcos tangenciales y morfologías de curvas críticas complejas impulsadas por fusiones de cúmulos y subestructura.
• Impacto de la Estructura de la Línea de Visión (LoS): El estudio cuantifica el impacto de la estructura de la LoS no correlacionada comparando el trazado de rayos de cono de luz completo contra el trazado de rayos a través únicamente del plano del lente primario.
  • Posiciones de las Imágenes: La estructura de la LoS no correlacionada puede desplazar las posiciones relativas de las imágenes lenticuladas por varios segundos de arco (hasta $\sim 10''$ en casos extremos), superando significativamente los residuos típicos de modelado de masa de lentes fuertes ($\lesssim 1''$).
  • Morfología de la Curva Crítica: La inclusión de material de la LoS introduce una dispersión de aproximadamente un $6\%$ en el área de la curva crítica primaria de un cúmulo. En algunos casos, el área de la curva crítica primaria se incrementó en un $30\%$ o disminuyó en un $50\%$.
  • Área Crítica Total: Dentro de $100''$ del potencial mínimo del cúmulo, el área crítica total en $z_s = 4$ cambia en $16^{+20}_{-14}\%$ cuando se incluye la estructura de la LoS. Este incremento surge tanto de la creación de nuevas curvas críticas por objetos de la LoS como de la modificación de las curvas existentes.

Significancia y Limitaciones
El artículo afirma proporcionar la primera demostración de imágenes de lente fuerte de cúmulos generadas de manera consistente a partir de un cono de luz de una simulación hidrodinámica cosmológica. Este enfoque ofrece un método más automatizado y "sin intervención manual" para producir grandes catálogos de observaciones simuladas, que pueden utilizarse para estudiar sistemáticas en el modelado de lentes, sesgos en la inferencia de masa y la tensión en la sección eficaz de lente fuerte galaxia-galaxia (GGSL).

Los autores señalan varias limitaciones:

• Potencia de Gran Escala: El método no introduce modos de potencia de gran escala ausentes en la caja de simulación inicial, lo que podría sesgar el muestreo del espectro de potencia de gran escala.
  ం
• Resolución: El método está limitado por la resolución de masa de la simulación, la cual tiene dificultades para resolver simultáneamente cúmulos masivos y subestructura/fuentes de baja masa.
• Masa Difusa: La masa difusa se aproxima como una lámina uniforme en lugar de preservar su estructura específica, aunque esto es necesario para evitar el ruido de disparo.
• Costo Computacional: Generar un solo cono de luz requiere procesar cada partícula en la caja de la simulación en cada instantánea, lo cual es computacionalmente costoso y genera requisitos significativos de almacenamiento de datos.
• Cizalladura Externa (External Shear): La implementación actual no incluye la cizalladura externa de la estructura de gran escala fuera del cono de luz, aunque se estima que su efecto es pequeño ($\gamma_{ext} \sim 0.02$).

Los autores concluyen que incorporar la información completa del cono de luz es crucial para las predicciones teóricas y que es necesario contabilizar cuidadosamente la estructura de la LoS al interpretar las observaciones, ya que la estructura no correlacionada puede alterar significativamente las configuraciones de las imágenes y las propiedades de las curvas críticas.