Impact and interplay of $Λ$CDM analysis choices for LSST cosmic shear

Resumen Técnico: Impacto e Interacción de las Elecciones de Análisis ΛCDM para el Cizallamiento Cósmico de LSST

1. El Problema
El Observatorio Vera C. Rubin (LSST) y la misión Euclid están diseñados para realizar mediciones de precisión del modelo cosmológico estándar ($\Lambda$CDM) y de la energía oscura mediante el cizallamiento cósmico (lente gravitacional débil). Sin embargo, existe una tensión crítica entre las proyecciones teóricas optimistas y la realidad de las incertidumbres sistemáticas.

El documento identifica dos fuentes principales de degradación en la precisión de los parámetros cosmológicos, específicamente en el parámetro compuesto $S_8 = \sigma_8 \sqrt{\Omega_m/0.3}$:

• Retroalimentación Bariónica (Baryon Feedback): La incertidumbre sobre cómo los procesos astrofísicos (principalmente la retroalimentación de núcleos galácticos activos y formación estelar) redistribuyen la materia en escalas no lineales. Los modelos actuales sugieren que este efecto puede suprimir la agrupación de materia en escalas mucho más grandes de lo que se creía anteriormente (hasta $k \sim 0.15 h \text{Mpc}^{-1}$), afectando significativamente el espectro de potencia de cizallamiento.
• Calibración de Corrimiento al Rojo Fotométrico ($z_{phot}$): Para el primer año de observaciones (Y1), existen muestras espectroscópicas representativas para calibrar los $z_{phot}$. Sin embargo, para el décimo año (Y10), la profundidad de LSST excederá la capacidad de las muestras espectroscópicas actuales. Esto obligará a utilizar técnicas de correlación cruzada, que actualmente tienen una precisión un orden de magnitud menor que la requerida por los documentos de requisitos científicos (DESC-SRD).

El riesgo es que, si no se modelan correctamente estas sistemáticas, se introduzcan sesgos (bias) en los parámetros cosmológicos recuperados, invalidando la capacidad de LSST para resolver tensiones actuales en cosmología o detectar nueva física.

2. Metodología
Los autores desarrollan un marco de análisis robusto utilizando simulaciones de datos ficticios (mocks) para LSST-Y1 y Y10, basándose en el modelo cosmológico $\Lambda$CDM plano.

• Simulaciones y Datos Ficticios: Utilizan las simulaciones hidrodinámicas FLAMINGO para generar los datos de "verdad" (mocks), que incluyen efectos realistas de retroalimentación bariónica. Comparan estos datos contra modelos teóricos de análisis.
• Modelado de Sistemáticas:
  • Retroalimentación Bariónica: Evalúan dos enfoques principales: HMCODE2020 (un modelo de halo con un parámetro libre $\Theta_{AGN}$) y SP(K) (un emulador flexible basado en la fracción bariónica $f_b(M_{halo})$ entrenado con las simulaciones ANTILLES).
  • Corrimiento al Rojo: Modelan la incertidumbre en el desplazamiento medio del corrimiento al rojo ($\delta_z$). Para Y1 asumen una precisión de $\delta_z \sim 0.002(1+z)$ (basada en la encuesta C3R2), mientras que para Y10 asumen una degradación a $\delta_z \sim 0.015$ debido a la falta de espectroscopía representativa.
  • Alineación Intrínseca (IA): Comparan el modelo NLA-z (dos parámetros) con el modelo TATT (cinco parámetros).
• Estrategia de Análisis:
  • Implementan cortes de escala ($K_{max}$) para eliminar datos en escalas altamente no lineales donde la física es incierta.
  • Marginalizan sobre los parámetros de sistemáticas utilizando priores informativos (basados en simulaciones FLAMINGO) y no informativos.
  • Utilizan una métrica de éxito estricta: recuperar el valor de verdad de $S_8$ dentro de un sesgo de $< 0.5\sigma$ respecto al valor real, corrigiendo los efectos de proyección (projection bias) mediante análisis de "verdad" y "prueba" personalizados.
• Herramientas: Se utiliza el código COSMOSIS y el muestreador NAUTILUS para explorar el espacio de parámetros, evitando la aproximación de la matriz de Fisher tradicional para capturar no linealidades.

3. Contribuciones Clave

• Reevaluación de la Precisión de Y10: Demuestran que, bajo las condiciones actuales de calibración de $z_{phot}$ para profundidades extremas, las restricciones de cizallamiento cósmico de LSST en el año 10 (Y10) apenas superan las del año 1 (Y1), a pesar de la mayor densidad de fuentes.
• Validación de Modelos de Retroalimentación: Comparan la eficacia de HMCODE2020 frente a emuladores más complejos como SP(K), mostrando que la flexibilidad del modelo no compensa automáticamente la falta de priores observacionales estrictos.
• Análisis de Sesgo de Proyección: Introducen una metodología rigurosa para cuantificar y corregir el "sesgo de proyección" que surge de la interacción entre priores asimétricos y degeneraciones no lineales, asegurando que las métricas de éxito no sean engañosas.
• Contraste con Estudios Previos: Explican por qué sus resultados difieren de proyecciones anteriores (como Boruah et al. 2024), atribuyendo la discrepancia a la diferencia entre analizar el modelo $\Lambda$CDM (donde $S_8$ es crítico) frente a modelos de energía oscura ($w_0w_a$CDM) en un análisis combinado de 3x2 puntos.

4. Resultados Principales

• Impacto de la Retroalimentación Bariónica: Cuando se tiene en cuenta la incertidumbre actual sobre la supresión de la retroalimentación bariónica, el error en $S_8$ se casi duplica en comparación con un análisis que ignora este fenómeno o asume un modelo perfecto.
• Compromiso entre Precisión y Exactitud:
  • Un análisis de "solo cortes de escala" (sin modelar la retroalimentación) requiere cortes muy agresivos ($K_{max} \approx 0.1$) para evitar sesgos, lo que triplica el error en $S_8$.
  • Un análisis que marginaliza sobre la retroalimentación con un prior informativo basado en FLAMINGO ($K_{max} \approx 0.5$) logra un equilibrio: reduce el sesgo a niveles aceptables ($<0.5\sigma$) con una pérdida de precisión moderada (el error es 1.5 veces mayor que en el caso sesgado).
• El Dilema de Y10: La degradación en la calibración de $z_{phot}$ para Y10 ($\delta_z = 0.015$) anula las ganancias estadísticas de la mayor profundidad. El análisis de Y10 con las sistemáticas actuales ofrece restricciones en $S_8$ casi idénticas a las de Y1.
• Modelos de Alineación Intrínseca: El cambio del modelo NLA-z al modelo TATT no afecta significativamente las restricciones en este escenario de simulación, aunque los autores advierten que esto depende de la verdad subyacente de las alineaciones.
• Potencial Futuro: Si se logran priores estrictos sobre la retroalimentación bariónica (mediante datos externos como DESxROSAT) y se mejora la calibración de $z_{phot}$ para Y10 a niveles de precisión de Y1, las restricciones de LSST en $S_8$ podrían ser más de cinco veces más precisas que las encuestas actuales (como DES y KiDS).

5. Significado y Conclusión
El estudio concluye que la comunidad cosmológica no puede depender únicamente de cortes de escala o de modelos teóricos imperfectos para mitigar la retroalimentación bariónica. Para que LSST y Euclid cumplan sus objetivos científicos, es imperativo:

1. Desarrollar metodologías observacionales y teóricas para restringir los modelos de retroalimentación bariónica con mayor precisión.
2. Avanzar en las técnicas de calibración de corrimiento al rojo fotométrico para profundidades donde la espectroscopía no es representativa.

El trabajo enfatiza que, sin abordar estas sistemáticas, el potencial estadístico masivo de LSST no se traducirá en una mejora real en la precisión cosmológica, y las tensiones actuales en los parámetros cosmológicos podrían seguir siendo ambiguas debido a errores sistemáticos no modelados. La prioridad debe ser entender la física de la retroalimentación bariónica para desbloquear la verdadera potencia de las encuestas de próxima generación.