KiDS-1000: Improved constraints on cosmology, intrinsic alignments and baryonic feedback from clipped cosmic shear

Resumen Técnico: KiDS-1000: Restricciones mejoradas sobre cosmología, alineaciones intrínsecas y retroalimentación bariónica mediante cizalladura cósmica recortada (clipped cosmic shear)

Planteamiento del Problema
El Modelo Estándar de Cosmología ($\Lambda$CDM) enfrenta tensiones persistentes con respecto a los valores de parámetros clave, particularmente la constante de Hubble ($H_0$) y la amplitud de agrupamiento de materia ($S_8$). Si bien el lente gravitacional débil (cizalladura cósmica o cosmic shear) es una sonda primaria para estos parámetros, los análisis tradicionales que dependen de estadísticas de dos puntos (por ejemplo, funciones de correlación de cizalladura) son subóptimos porque extraen información únicamente hasta el segundo orden. La distribución de materia en tiempos tardíos que sondea el lensing es altamente no gaussiana, conteniendo información de orden superior que las estadísticas de dos puntos no logran aprovechar plenamente. Además, las escalas pequeñas donde reside esta información no gaussiana están fuertemente contaminadas por la retroalimentación bariónica (baryonic feedback) y las alineaciones intrínsecas (IA), lo que dificulta un modelado robusto. Los intentos previos de utilizar estadísticas de orden superior han tenido dificultades para modelar conjuntamente estos sistemas con la precisión suficiente. Adicionalmente, aunque el "recorte" (clipping —eliminar regiones de alta densidad del campo proyectado—) ha sido propuesto para decorrelacionar los modos de Fourier y mejorar la extracción de información de las estadísticas de dos puntos, aplicaciones previas a datos reales de lensing (por ejemplo, Giblin et al. 2018) carecieron de un modelado exhaustivo de sistemas como las alineaciones intrínsecas, la retroalimentación bariónica y las incertidumbres en los desplazamientos al rojo fotométricos.

Metodología
Este artículo presenta las primeras restricciones cosmológicas del sondeo KiDS-1000 utilizando funciones de correlación de cizalladura recortadas, empleando un enfoque de modelado directo (forward-modeling) basado en simulaciones de $N$-cuerpos y emuladores de Procesos Gaussianos (GP).

• Datos y Simulaciones: El análisis utiliza el cuarto lanzamiento de datos del Kilo-Degree Survey (KiDS-1000), que cubre 777 deg$^2$ con 21 millones de galaxias. Los autores se apoyan en las suites de simulaciones de materia oscura únicamente cosmoSLICS y SLICS. Estas simulaciones ofrecen volúmenes significativamente mayores y una mayor resolución que los conjuntos de datos utilizados anteriormente para el lensing recortado.
• Procedimiento de Recorte (Clipping): La metodología sigue el enfoque de Giblin et al. (2018). Se generan mapas de convergencia ($\kappa$) a partir de mapas de elipticidad mediante la técnica de mapeo de masa de Kaiser-Squires. Las regiones que exceden un umbral $\kappa_c = 0.010$ son "recortadas" (ajustadas al valor del umbral), creando un mapa recortado ($\kappa_{clip}$) y un mapa residual ($\Delta\kappa$). El mapa residual se invierte para obtener elisticidades residuales, que se restan de las elipticidades observadas para producir elisticidades recortadas.
• Modelado de Sistemas: Un avance crítico es el modelado exhaustivo de los sistemas. Los autores entrenan emuladores de GP sobre variaciones de las simulaciones para predecir la dependencia de las funciones de correlación de cizalladura tanto recortadas como no recortadas en función de:
  • Cosmología: Un espacio de parámetros $w$CDM de 4 dimensiones ($\Omega_m, S_8, h, w_0$).
  • Alineaciones Intrínsecas (IA): Modeladas mediante el modelo de alineación lineal no lineal (NLA) con amplitud $A_{IA}$.
  • Retroalimentación Bariónica (BF): Modelada utilizando simulaciones hidrodinámicas magneticum, escalada por un parámetro de molestia libre $b_{bary}$.
  • Desplazamientos al rojo fotométricos ($z$): Modelados mediante desplazachientos ($\Delta z$) en las distribuciones de desplazamiento al rojo.
  • Agrupamiento Fuente-Lente (SLC): Se investigó y se encontró que es insignificante para este estadístico.
• Estrategia de Emulación: Los emuladores de GP se entrenan para predecir los sesgos cosmológicos y de sistemas con respecto a un modelo fiduciario. Los emuladores se validan mediante validación cruzada de dejar uno fuera (leave-one-out), mostrando errores subdominantes comparados con el ruido estadístico. La incertidumbre del emulador se incluye en la matriz de covarianza para evitar sesgos.
• Inferencia de Parámetros: El análisis utiliza una verosimilitud gaussiana (con corrección de Hartlap para la inversa de la covarianza) y el muestreador nautilus para explorar el espacio de parámetros de 11 dimensiones (4 cosmológicos + 7 de molestia). El vector de datos consiste en auto-correlaciones y correlaciones cruzadas para cinco bins de desplazamiento al rojo tomográficos, midiendo tanto $\xi_+$ como $\xi_-$ recortados, resultando en un vector de longitud 540.

Contribuciones Clave

1. Primer Análisis de Recorte Controlado por Sistemas: Esta es la primera aplicación de recorte a datos reales de lensing débil que modela simultáneamente las alineaciones intrínsecas, la retroalimentación bariónica, las incertidumbres de desplazamiento al rojo fotométrico y el agrupamiento fuente-lente.
2. Modelado Directo con Simulaciones: A diferencia de análisis previos que dependían de prescripciones teóricas (por ejemplo, hmcode) para las estadísticas no recortadas, este trabajo modela directamente tanto las estadísticas recortadas como las no recortadas utilizando simulaciones de $N$-cuerpos, proporcionando una validación independiente de los resultados de cizalladura cósmica de KiDS-1000 (Asgari et al. 2021).
3. Recorte Tomográfico: El análisis introduce el agrupamiento de desplazamiento al rojo tomográfico dentro de un marco de lensing recortado por primera vez, expandiendo significativamente el vector de datos en comparación con estudios previos de prueba de concepto.
4. Manejo del Sesgo de Máscara: Al replicar directamente la huella del sondeo KiDS-1000 (incluyendo la máscara) en los mosaicos de simulación, los autores modelan directamente el sesgo de máscara, permitiendo el uso de ambos componentes $\xi_+$ y $\xi_-$ del estadístico recortado, que anteriormente estaban sesgados por la máscara.

Resultos
La combinación de sondas recortadas y no recortadas produce restricciones más estrictas sobre los parámetros cosmológicos y de molestia en comparación con el enfoque no recortado por sí solo:

• Restricciones Cosmológicas: El análisis arroja $\Omega_m = 0.263^{+0.035}_{-0.038}$, $S_8 = 0.724^{+0.027}_{-0.027}$, y $w_0 = -1.24^{+0.26}_{-0.28}$.
  • La incertidumbre sobre $S_8$ se reduce en un 16%.
  • La incertidumbre sobre $w_0$ se reduce en un 24%.
  • La Figura de Mérito (FoM) en el plano $\Omega_m-S_8$ mejora por un factor de 1.2.
• Restricciones de Sistemas:
  • La restricción de la amplitud de alineación intrínseca ($A_{IA}$) mejora en un 27%, alcanzando una precisión de ~39% (comparado con ~90% en el análisis no recortado).
  • El parámetro de retroalimentación bariónica ($b_{bary}$), que no está restringido en el análisis no recortado, se restringe a un límite superior de $b_{bary} < 0.97$ en el análisis combinado.
• Consistencia: Los resultados son consistentes con el modelo $\Lambda$CDM y el análisis principal de cizalladura cósmica de KiDS-1000 (Asgari et al. 2021), a pesar de utilizar un marco de simulación y emulador completamente independiente. El ajuste (goodness-of-fit, $p$-valores) es aceptable para todos los análisis.
• Tensiones: La restricción de $S_8$ presenta una tensión leve (1.8$\sigma$) con el resultado de A21 y una de 3.0$\sigma$ con los resultados de Planck CMB, aunque los autores señalan que las diferencias en la calibración de desplazamiento al rojo y las especificaciones del sondeo en los datos más nuevos (KiDS-Legacy) pueden desplazar estos valores.

Significado y Afirmaciones
El artículo afirma que el recorte transforma la cizalladura cósmica de una técnica de "prueba de concepto" en un estadístico de lensing competitivo y robusto. Al filtrar las regiones de alta densidad, el recorte permite que las estadísticas de dos puntos extraigan información cosmológica adicional de escalas más grandes y lineales donde la contaminación bariónica es más débil, mientras que la información complementaria de la estadística no recortada ayuda a romper las degeneraciones.

Los autores enfatizan que la naturaleza complementaria de las sondas recortadas y no recortadas permite el desacoplamiento de los procesos dependientes de la escala. Esto se evidencia por la mejora significativa en la restricción de la amplitud de alineación intrínseca y en establecer un límite superior para la fuerza de la retroalimentación bariónica, parámetros que son difíciles de restringir con estadísticas no recortadas por sí solas. El trabajo demuestra que el modelado directo basado en simulaciones puede controlar con éxito los sistemas en estadísticas de lensing de orden superior o modificadas, validando el enfoque para futuros sondeos de Etapa IV como Euclid y LSST. Los autores concluyen modestamente que, si bien persisten desafíos (por ejemplo, imprecisiones del emulador en escalas pequeñas), los obstáculos restantes para utilizar plenamente el recorte son probablemente de corta duración dados los avances actuales en la calidad de los datos y los recursos de simulación.