Comparing cosmic shear nulling methods for Stage-IV surveys

Resumen Técnico: Comparación de Métodos de Anulación de Cizalla Cósmica para Sondeos de la Fase IV

Enunciado del Problema
Las mediciones de cizalla cósmica procedentes de sondeos de la Fase IV (por ejemplo, Euclid, Observatorio Rubin) tienen como objetivo restringir la física de la Energía Oscura a nivel porcentual. Sin embargo, estas mediciones se ven contaminadas por procesos de retroalimentación bariónica (por ejemplo, retroalimentación de AGN, supernovas), que suprimen el espectro de potencia de la materia a pequeñas escalas ($k \gtrsim 10 \, \mathrm{Mpc}^{-1}$). Las simulaciones hidrodinámicas actuales predicen amplitudes de supresión que varían significativamente (del 10% a >30%), y las prescripciones de física subgrid permanecen inciertas. El análisis estándar mitiga esto aplicando cortes de escala angular (limitando $\ell_{\max}$), pero este enfoque es subóptimo porque una escala angular dada integra contribuciones de una amplia gama de escalas físicas 3D debido al amplio kernel de lentes radial. En consecuencia, se requieren cortes de escala agresivos para garantizar restricciones cosmológicas sin sesgo, lo que conduce a una pérdida significativa de precisión estadística. Este artículo investiga si las transformaciones de "anulación" aplicadas al vector de datos pueden aislar y eliminar más eficazmente las contribuciones de alto-$k$ (pequeñas escalas) en comparación con los cortes de escala estándar.

Metodología
Los autores realizan un análisis de pronóstico de Fisher comparando tres estrategias de anulación distintas frente a un análisis estándar de cizalla cósmica sin anulación. El análisis utiliza un vector de datos simulado basado en una configuración de sondeo similar a Euclid con seis bins de fuentes tomográficas, una densidad numérica de galaxias de 30 galaxias/min$^2$ y un área de cielo de 15.000 deg$^2$. La retroalimentación bariónica se modela utilizando HMCode2020 con un parámetro de calentamiento de AGN fiducial ($\Theta_{\mathrm{AGN}} = 7.8$). Los autores asumen una cosmología $w_0w_a$CDM y fijan todos los sistemáticos (por ejemplo, corrimientos al rojo fotométricos, sesgo de cizalla) a cero para aislar el impacto de la retroalimentación bariónica.

Los tres métodos comparados son:

1. Anulación LU (LUnul): Una nueva implementación de un método propuesto originalmente por Huterer & White (2005). Trata la integral de Limber como una multiplicación matricial discreta. Al realizar una factorización LU de la matriz del kernel, el método construye combinaciones lineales de los espectros de potencia tomográficos. Esto cancela efectivamente las contribuciones comunes de alto-$k$ (que son similares entre bins a bajo corrimiento al rojo) mientras retiene las diferencias a bajo $k$. El método permite un corte "suave" donde se descartan filas específicas del vector de datos transformado basándose en un parámetro de tolerancia ($tol=0.01$) para las contribuciones de alto-$k$.
2. Transformación BNT: La transformación de Bernardeau-Nishimichi-Taruya opera sobre el campo de lentes (mapas de cizalla). Utiliza una transformación lineal e invertible de los índices de los bins tomográficos para localizar los kernels de lentes en el espacio de corrimiento al rojo. Esto crea un mapeo más limpio entre el multipolo angular $\ell$ y el número de onda 3D $k$, suprimiendo las correlaciones cruzadas entre bins no adyacentes.
3. Correlación Cruzada: Adaptada de la literatura sobre lentes del CMB, este método utiliza una muestra de agrupamiento de galaxias de primer plano (bins de lentes de bajo corrimiento al rojo) para descorrelacionar la señal de cizalla. Al restar el campo de densidad proyectado (ponderado por el espectro de correlación cruzada) del campo de cizalla, se suprime la contribución de estructuras de bajo corrimiento al rojo y pequeña escala. En este estudio, los dos bins de corrimiento al rojo más bajos se reutilizan como lentes para "limpiar" los tres bins de fuentes de corrimiento al rojo más altos.

Contribuciones Clave

• Implementación de LUnul: Los autores proporcionan una implementación pública del método de factorización LU para cizalla cósmica, una técnica no explorada previamente en la literatura desde su propuesta inicial.
• Aplicación de Correlación Cruzada a la Cizalla Cósmica: El artículo adapta la técnica de anulación por correlación cruzada, utilizada previamente para lentes del CMB, al contexto de la cizalla cósmica de la Fase IV combinada con agrupamiento de galaxias fotométricas.
• Comparación Sistemática: El estudio ofrece una comparación directa de estos tres enfoques distintos (basado en campos, basado en estadísticas y basado en trazadores) dentro de un marco unificado de pronóstico de Fisher, evaluando su eficacia para reducir el sesgo mientras se retiene información.

Resultados
El análisis evalúa el compromiso entre el sesgo de parámetros (precisión) y las barras de error (precisión) a través de una gama de números de onda máximos ($k_{\max}$):

• Reducción del Sesgo: Los tres métodos reducen eficazmente el sesgo en los parámetros cosmológicos ($S_8, \Omega_m, w_0, w_a$) causado por la retroalimentación bariónica no modelada en comparación con el caso estándar.
• Rendimiento de LUnul: Este método es el más agresivo en la supresión de modos de alto-$k$. Logra el sesgo más bajo (reduciendo el sesgo de $S_8$ a $<0.5\sigma$ en $k_{\max} = 0.3 \, h/\mathrm{Mpc}$), superando a otros métodos en precisión. Sin embargo, esto conlleva un costo en precisión; el poder de restricción (Figura de Mérito, FoM) es significativamente menor que en el caso estándar en cortes de escala comparables porque descarta más información.
• Rendimiento de BNT: La transformación BNT preserva más información que LUnul y tiene una base teórica rigurosa. Sin embargo, en esta configuración específica con solo seis bins tomográficos, muestra solo una mejora marginal sobre el caso estándar. Los autores señalan que BNT es conocido por ser más efectivo con un mayor número de bins.
• Rendimiento de Correlación Cruzada: Este método es efectivo para reducir los sesgos absolutos en $S_8$, pero requiere información adicional del agrupamiento de galaxias. Aunque produce barras de error más ajustadas (debido a la inclusión de información del campo de densidad), el sesgo relativo (sesgo dividido por el error) permanece alto. Además, es menos restrictivo sobre los parámetros de Energía Oscura ($w_0, w_a$) en comparación con otros métodos. El método está limitado por la superposición entre muestras de lentes y fuentes en sondeos fotométricos.
• Sensibilidad de Parámetros: $S_8$ es consistentemente el parámetro más sesgado en todos los métodos debido a su sensibilidad a la amplitud del espectro de potencia.

Significado y Afirmaciones
El artículo concluye que, si bien las técnicas de anulación pueden mitigar significativamente el sesgo introducido por la retroalimentación bariónica, generalmente lo hacen a expensas de la precisión estadística. Los autores afirman que:

1. LUnul es el método más robusto para eliminar el sesgo y avanzar hacia escalas más pequeñas, aunque requiere un ajuste numérico cuidadoso y sacrifica precisión.
2. BNT es teóricamente sólido y preserva información, pero puede requerir más bins tomográficos (típicos de futuros lanzamientos de datos de la Fase IV) para realizar su pleno potencial sobre los cortes estándar.
3. Correlación cruzada ofrece una vía viable para reducir el sesgo utilizando datos existentes de 3$\times$2-pt, pero está limitada por la dificultad de separar perfectamente las poblaciones de lentes y fuentes en sondeos fotométricos.

Los autores enfatizan que estos métodos proporcionan una alternativa a la modelización compleja de la retroalimentación bariónica. Afirman que los sondeos de la Fase IV probablemente necesitarán emplear cortes de escala fuertes o depender de modelos bariónicos precisos independientemente de la elección de anulación, pero estas estrategias de anulación ofrecen una forma de construir vectores de datos que son más inmunes a las incertidumbres astrofísicas a pequeña escala. Los resultados refuerzan el compromiso establecido entre precisión (reducción del sesgo) y exactitud en los análisis de cizalla cósmica.