Balancing bias, baryons, and scale cuts in LSST 3x2pt analysis

Resumen Técnico: Equilibrio entre Sesgo, Bariones y Cortes de Escala en el Análisis $3 \times 2$pt de LSST

Planteamiento del Problema
Los sondeos cosmológicos de Etapa IV, como el Legacy Survey of Space and Time (LSST) del Observatorio Rubin, sondearán el régimen no lineal de la formación de estructuras donde los efectos sistemáticos de la formación de galaxias (bias) y la retroalimentación bariónica se vuelven dominantes. En un análisis $3 \times 2$pt (que combina el agrupamiento de galaxias, el cizallamiento cósmico y su correlación cruzada), estos efectos introducen parámetros de confusión (nuisance parameters) que a menudo están débilmente restringidos y pueden exhibir degeneraciones con los parámetros cosmológicos. Específicamente, la interacción entre el sesgo galáctico no lineal y la supresión bariónica en el espectro de potencia de la materia plantea un desafío: los modelos de sesgo puramente lineales son insuficientes en escalas pequeñas ($k \gtrsim 0.1 \, h/\text{Mpc}$), pero los modelos perturbativos más complejos introducen grados de libertad adicionales que pueden ser degenerados con la física bariónica. El problema central abordado es cómo equilibrar la complejidad del modelo frente a los cortes de escala para maximizar el poder de restricción sobre los parámetros cosmológicos ($\Omega_m, \sigma_8$) y sus extensiones (masa de neutrinos $M_\nu$) sin introducir sesgos significativos.

Metodología
Los autores presentan un análisis exhaustivo utilizando un nuevo canal de inferencia bayesiana de código abierto llamado MGL (Modified Gravity Lensing). El estudio genera vectores de datos simulados (mock data) que asemejan las especificaciones de LSST Año 1 (Y1) y Año 10 (Y10), incluyendo espectros de cizallamiento cósmico y agrupamiento de galaxias en 20 bins logarítmicos espaciados en $\ell \in [20, 5000]$.

Componentes metodológicos clave:

• Marco de Modelado: El análisis emplea el emulador BACCO para modelar tanto el espectro de potencia de la materia no lineal como la retroalimentación bariónica mediante un mecanismo de "barianificación". Para el sesgo de las galaxias, se comparan tres modelos:
  1. Sesgo Lineal: Válido solo en escalas grandes.
  2. Teoría de Campo Efectiva Híbrida (HEFT): Un enfoque perturbativo que combina la expansión de sesgo lagrangiano con la evolución de simulaciones N-body, válido hasta $k_{\text{max}} \approx 0.7 \, h/\text{Mpc}$.
  3. Sesgo Mínimo: Una variante simplificada de HEFT donde los términos de sesgo de orden superior ($b_{s^2}$ y $b_{\nabla^2}$) se fijan en cero.
• Tratamiento Bariónico: La retroalimentación bariónica se modela utilizando un factor de supresión $S_{mm}$ aplicado al espectro de potencia de la materia. Para el lensing galaxia-masa ($P_{gm}$), los autores adoptan la aproximación de Zennaro et al. (2024), donde $P_{gm}$ se corrige por $\sqrt{S_{mm}}$, mientras que $P_{gg}$ depende de la flexibilidad de la expansión del sesgo.
• Estrategia de Análisis: Los autores realizan análisis de verosimilitud sobre datos simulados generados con el modelo HEFT completo y parámetros fiduciales. Prueban diversas suposiciones de análisis, incluyendo modelos desajustados (por ejemplo, analizar datos generados con HEFT usando un modelo de sesgo lineal o mínimo) y la inclusión variable de parámetros bariónicos.
• Métricas: El rendimiento se cuantifica mediante la Figura de Mérito (FoM) para el poder de restricción y la Figura de Sesgo (FoB) para medir las desviaciones de los valores fiduciales.

Contribuciones Clave

1. Desarrollo del Pipeline MGL: El artículo presenta MGL, un pipeline independiente y de código abierto diseñado para sondeos de Etapa IV que utiliza emuladores para el cálculo eficiente de los espectros de potencia y soporta el análisis MCMC con el muestreador nautilus.
2. Comparación Sistemática de Modelos de Sesgo: El trabajo proporciona una comparación directa de los modelos de sesgo lineal, HEFT y mínimo en el contexto de LSST Y1 y Y10, evaluando específicamente su robustez frente a la retroalimentación bariónica y los cortes de escala.
3. Análisis de Degeneración: El estudio investiga explícitamente la degeneración entre los parámetros de sesgo de orden superior y la retroalimentación bariónica, así como la interacción entre el modelado del sesgo y la detección de neutrinos masivos.

Resultados

• Cortes de Escala y Modelos de Sesgo: Un modelo de sesgo lineal produce restricciones no sesgadas para $\Omega_m$ y $\sigma_8$ (FoB $< 0.5\sigma$) solo hasta $k_{\text{max}} \approx 0.1 \, h/\text{Mpc}$. Extenderse a escalas más pequeñas requiere un enfoque perturbativo.
• Eficacia del Sesgo Mínimo: El modelo de "sesgo mínimo" (fijando $b_{s^2}$ y $b_{\nabla^2}$ a cero) se comporta de manera comparable al modelo HEFT completo. Permanece sin sesgo para $\Omega_m$ y $\sigma_8$ incluso en $k_{\text{max}} = 0.7 \, h/\text{Mpc}$, siempre que los parámetros de sesgo fiduciales se encuentren dentro de los rangos estándar. Esto sugiere que la complejidad total de HEFT puede no ser necesaria para las restricciones cosmológicas estándar, ofreciendo una vía para reducir el espacio de parámetros.
• Interacción Sesgo-Barión: Los parámetros de sesgo de orden superior pueden imitar los efectos de supresión de la retroalimentación bariónica en el espectro de potencia galaxia-masa. Sin embargo, la retroalimentación bariónica no puede reproducir completamente el rango de comportamientos generados por los términos de sesgo de orden superior. Omitir las correcciones bariónicas en el componente de lensing galaxia-masa no sesga significativamente a $\Omega_m$ o $\sigma_8$, aunque puede desplazar el parámetro bariónico inferido ($\log_{10} M_c$).
• Detección de la Masa de Neutrinos: Una detección de la masa total de los neutrinos ($M_\nu$) es posible tanto para Y1 como para Y10 al extenderse a escalas no lineales ($k \geq 0.3 \, h/\text{Mpc}$). Sin embargo, el valor específico de $M_\nu$ inferido no es robusto en todos los escenarios plausibles. Cuando se analizan datos simulados con parámetros de sesgo de orden superior elevados utilizando el modelo de sesgo mínimo, ocurren sesgos significativos en la $M_\nu$ inferida (hasta $\sim 2\sigma$ para Y10). Esto indica una fuerte degeneración entre la masa de los neutrinos, los términos de sesgo de orden superior y los efectos de proyección.
• Y1 vs. Y10: Pasar de Y1 a Y10 mejora la FoM por un factor de $\sim 2.3$ en $k_{\text{max}} = 0.7 \, h/\text{Mpc}$. Mientras que $\Omega_m$ y $\sigma_8$ siguen siendo robustos, la fragilidad de las restricciones de $M_\nu$ bajo modelos de sesgo simplificados persiste.

Significancia y Reivindicaciones
El artículo afirma que para los parámetros estándar de $\Lambda$CDM ($\Omega_m, \sigma_8$), una parametrización de sesgo mínima es suficiente para extraer restricciones no sesgadas de los datos de LSST hasta el límite de validez de HEFT ($k_{\text{max}} = 0.7 \, h/\text{Mpc}$), incluso en presencia de retroalimentación bariónica. Esto permite la inclusión de escalas cuasi-no lineales sin la penalización computacional y estadística de una expansión de sesgo completa de 15 términos.

Sin embargo, los autores advierten modestamente que esta conclusión es sensible a los valores fiduciales específicos de los parámetros de sesgo. Aunque LSST debería ser capaz de una detección inequívoca de la masa de neutrinos no nula, la significancia y el valor preciso de esta detección dependen altamente de la precisión del modelado del sesgo. El estudio destaca que los enfoques de modelado simplificados (como el sesgo mínimo) pueden producir resultados robustos para la cosmología estándar, pero pueden introducir sesgos significativos para parámetros de extensión como la masa de los neutrinos si la física subyacente del sesgo galáctico se desvía del fiducial asumido. Por consiguiente, los autores abogan por realizar pruebas cuidadosas de los modelos simplificados contra casos científicos específicos antes de su adopción.