Resumen Técnico: Reconstrucción de las Escalas Más Grandes del Universo mediante la Inferencia a Nivel de Campo Aplicada al Catálogo de Cuásares Quaia

1. Planteamiento del Problema

El mapeo de la estructura a gran escala (LSS, por sus siglas en inglés) del Universo es esencial para extraer física fundamental; sin embargo, los análisis tradicionales que dependen de estadísticas de resumen (por ejemplo, espectros de potencia, bispectros) sufren de una compresión de datos subóptima, al no lograr explotar plenamente la información de alta dimensión. Si bien los cuásares son trazadores únicos capaces de sondear el Universo hasta desplazamientos al rojo $z \sim 7$ y cubrir escalas ultra grandes, su análisis se complica por las bajas densidades numéricas (alto ruido de disparo o shot noise), la contaminación de primer plano y las incertidumbres en el desplazamiento al rojo (particularmente para muestras fotométricas).

El recientemente publicado catálogo de cuásares Quaia, derivado de los datos de Gaia y unWISE, ofrece una cobertura de todo el cielo ($f_{sky} \approx 0.71$) y un amplio rango de desplazamiento al rojo ($0.084 < z \lesssim 4.0$). No obstante, aprovechar este conjunto de datos para una reconstrucción tridimensional completa de la distribución de materia requiere un método que pueda manejar simultáneamente los efectos de selección del sondeo, las distorsiones en el espacio de redshift y el sesgo (bias) de los cuásares, así como la inherente escasez de la población de trazadores. Los métodos existentes suelen tener dificultades para integrar estos sistemas en un modelo físico coherente sin perder información.

2. Metodología

Este trabajo aplica el algoritmo BORG (Bayesian Origin Reconstruction from Galaxies), un marco de inferencia a nivel de campo (FLI, por sus siglas en inglés), a los catálogos de cuásares Quaia. A diferencia de los enfoques tradicionales que comprimen los datos en estadísticas de resumen, BORG realiza un modelado bayesiano hacia adelante (forward-modeling) de todo el campo galáctico observado.

2.1 Muestras de Datos

El análisis utiliza dos muestras distintas del catálogo Quaia (Storey-Fisher et al. 2024):

• Quaia Clean: Una muestra limitada por magnitud ($G < 20.0$) con $\sim 755,850$ fuentes.
• Quaia Deep Cut: Una muestra extendida ($G < 20.5$) con $\sim 1,295,502$ fuentes. Esta muestra se somete a una selección angular más estricta, excluyendo regiones cercanas a las Nubes de Magallanes Grande y Pequeña (LMC/SMC) e imponiendo un umbral de completitud $> 0.5$ para mitigar la contaminación de primer plano.

Ambas muestras se dividen en ocho contenedores radiales para sondear diferentes épocas cósmicas manteniendo el poder estadístico. Los desplazamientos al rojo se transforman al marco de reposo del CMB para eliminar el dipolo inducido por el movimiento peculiar del Sol.

2.2 Modelo de Propagación (Forward Model)

El algoritmo BORG emplea un modelo de propagación basado en la física para evolucionar las condiciones iniciales hacia los observables actuales:

1. Condiciones Iniciales: Un campo de ruido blanco gaussiano es modulado por un espectro de potencia primordial (fijado a la cosmología de Planck 2018) para generar el potencial gravitatorio inicial.
2. Evolución Gravitatoria: El campo de densidad se evoluciona utilizando la Teoría de Perturbaciones Lagrangiana (LPT) (de primer orden) para modelar la formación de estructuras desde las condiciones iniciales hasta el presente. Este enfoque es computacionalmente eficiente para las grandes escalas y la resolución gruesa ($39.1 \, h^{-1}\text{Mpc}$) considerada.
3. Efectos Observacionales: El modelo incorpora distorsiones en el espacio de redshift, efectos de cono de luz y la función de ventana del sondeo (incluyendo funciones de selección radial y máscaras de completitud angular).
4. Modelado de Sesgo (Bias): Se utiliza un modelo de sesgo de potencia de ley de Poisson para relacionar la densidad de materia oscura subyacente con la distribución observada de los cuásares.

2.3 Motor de Inferencia

El algoritmo explora la distribución posterior de alta dimensión de las condiciones iniciales y los campos tardíos utilizando Monte Carlo Hamiltoniano (HMC) combinado con muestreadores de rebanada (slice samplers). Esto genera un conjunto de realizaciones físicamente plausibles de:

• Las condiciones iniciales primordiales.
• El campo de densidad de materia oscura actual.
• Los campos de velocidad a gran escala (radiales y absolutos).
• La distribución de cuásares predicha.

2.4 Estrategia de Validación

Para validar las reconstrucciones, los autores realizan:

• Verificaciones de Consistencia Interna: Comparación de los espectros de potencia y bispectros recuperados contra las predicciones teóricas (CLASS) y comprobaciones de no gaussianidades en las condiciones iniciales.
• Correlación Cruzada Externa: Un paso crítico de validación implica la correlación cruzada del campo de densidad de materia inferido por BORG con el mapa de convergencia de lente gravitacional del CMB de Planck ($\kappa_{Planck}$). Esto pone a prueba si las estructuras de gran escala reconstruidas se alinean con la distribución de materia integrada observada mediante el efecto de lente gravitacional del CMB.

3. Contribuciones Clave y Resultados

3.1 Volumen de Reconstrucción Sin Precedentes

El estudio produce la reconstrucción de campo más grande del Universo observable hasta la fecha en términos de volumen comóvil.

• Volumen: $(10 \, h^{-1}\text{Gpc})^3$.
• Resolución: Resolución espacial máxima de $39.1 \, h^{-1}\text{Mpc}$ (tamaño del vóxel).
• Alcance: La reconstrucción cubre el volumen completo del sondeo Quaia, proporcionando mapas posteriores de las condiciones iniciales, la densidad de materia oscura actual y los campos de velocidad.

3.2 Campos Inferidos

• Campo de Densidad: El campo de densidad de materia oscura reconstruido muestra una fuerte consistencia cualitativa con la distribución observada de cuásares. Captura estructuras coherentes a gran escala (filamentos y vacíos) y se extiende más allá de los límites del sondeo mediante el acoplamiento gravitatorio, con incertidumbres que aumentan naturalmente en regiones de baja densidad de trazadores.
• Condiciones Iniciales: Las fluctuaciones primordiales inferidas son consistes con estadísticas gaussianas (media $\approx 0$, asimetría y curtosis consistentes con cero), lo que confirma la validez del modelo de propagación y la ausencia de no gaussianidades artificiales introducidas por el proceso de inferencia.
• Campos de Velocidad: Los campos de velocidad reconstruidos revelan flujos masivos coherentes a gran escala, con entradas y salidas que trazan la distribución de densidad subyacente. La relación señal-ruido sigue la función de ventana del sondeo, resaltando la capacidad del algoritmo para dar cuenta de las limitaciones observacionales.

3.3 Validación Estadística

• Espectro de Potencia y Bispectro: El espectro de potencia de materia inferido concuerda con las predicciones teóricas de CLASS en una amplia gama de escalas. El bispectro reducido es consistente con cero, lo que indica la ausencia de contribuciones no lineales espurias.
• Consistencia de Subvolúmenes: Los espectros de potencia calculados en subregiones del volumen de reconstrucción muestran consistencia estadística, descartando desajustes espaciales significativos del modelo o contaminación localizada.

3.4 Correlación Cruzada con Lente de CMB

La validación externa más significativa es la detección de una señal de correlación cruzada entre el campo de densidad reconstruido por BORG y el mapa de lentes de Planck.

• Significancia: Para la muestra Quaia Deep Cut con $N=256$ de resolución, la significancia de la detección alcanza aproximadamente $\sim 4\sigma$ (específicamente $3.96\sigma$ mediante ajuste de amplitud y $3.0\sigma$ mediante alineación de fase).
• Amplitud: El parámetro de amplitud ajustado es $A = 0.70 \pm 0.17$, consistente con la expectativa teórica ($A=1$) dentro de las incertidumbres.
• Implicación: Esta detección confirma que el marco de BORG recupera con éxito la verdadera distribución de materia a gran escala trazada por los cuásares, validando la reconstrucción frente a una sonda cosmológica independiente.

4. Significado y Afirmaciones

El artículo afirma establecer un marco robusto para explotar los sondeos de cuásares en la cosmología de campo. Su principal importancia radica en:

1. Primera Aplicación a Cuásares: Esta es la primera aplicación de la inferencia a nivel de campo a un conjunto de datos de cuásares, demostrando que los trazadores dispersos y ruidosos pueden utilizarse para reconstruir la distribución de materia 3D y las condiciones iniciales.
2. Manejo de Sistemas de Error (Systematics): El trabajo demuestra que un enfoque bayesiano jerárquico de propagación completa puede dar cuenta de manera efectiva de los complejos sistemas de error observacionales (efectos de selección, errores de redshift, sesgo) sin perder información debido a la compresión de datos.
3. Validación de la Física a Gran Escala: La detección de la señal de lente de CMB de $\sim 4\sigma$ proporciona una confirmación independiente de que los campos reconstruidos reflejan con precisión la estructura de materia a gran escala.
4. Base para la Cosmología Futura: Los autores postulan que este marco permite estudios de precisión de la física a escalas ultra grandes, incluyendo restricciones sobre la no gaussianidad primordial ($f_{NL}$), pruebas de la relatividad general e investigaciones del efecto de Sachs-Wolfe Integrado (ISW). Los campos reconstruidos sirven como un "gemelo digital" del Universo local, proporcionando productos de datos (mapas posteriores) para simulaciones con restricciones y análisis cosmológicos adicionales.

Los autores señalan que la resolución de la reconstrucción está actualmente limitada por el ruido de disparo y las incertididades en el desplazamiento al rojo espectrofotométrico, lo que impide una inferencia fiable en escalas menores. Sin embargo, los resultados confirman que la inferencia a nivel de campo es una herramienta viable y poderosa para transformar las observaciones dispersas de cuásares en mapas tridimensionales estadísticamente robustos de la estructura a gran escala del Universo.