EMU: Cross-correlating EMU Pilot Survey 1 with Dark Energy Survey to validate the radio galaxy bias and redshift distribution

Resumen Técnico: Correlación Cruzada de EMU con DES para Validar el Sesgo de Galaxias de Radio y la Distribución de Redshift

Planteamiento del Problema
Los sondeos de continuo de radio ofrecen un método poderoso para mapear la estructura a gran escala (LSS) del Universo debido a su capacidad para cubrir grandes áreas del cielo y su insensibilidad a la extinción por polvo. Sin embargo, una limitación crítica para los análisis cosmológicos utilizando estos sondeos es la falta de información directa de redshift (desplazamiento al rojo) para las fuentes de radio. El mecanismo de emisión dominante, la radiación sincrotrón, produce un espectro de ley de potencia sin rasgos distintivos, lo que hace que la estimación del redshift sea extremadamente desafiante. Sin una distribución de redshift ($dN/dz$) precisa y un modelo robusto para el sesgo de las galaxias, los errores sistemáticos pueden propagarse hacia las restricciones de los parámetros cosmológicos, como $\sigma_8$. Si bien estudios previos han dependido del emparejamiento con sondeos ópticos (lo que proporciona redshifts fotométricos solo para una fracción de las fuentes) o simulaciones (por ejemplo, TRECS, SKADS), existe la necesidad de un método que utilice toda el área del sondeo para inferir estadísticamente la distribución de redshift sin depender de la identificación de fuentes individuales.

Metodología
Este estudio emplea una técnica de correlación cruzada entre el mapa de continuo de radio bidimensional del Evolutionary Map of the Universe Pilot Survey 1 (EMU-PS1) y la muestra de galaxias fotométricas ópticas del Dark Energy Survey (DES). El análisis procede de la siguiente manera:

1. Preparación de Datos:

   • EMU-PS1: Utiliza un catálogo generado por el buscador de fuentes PyBDSF, que contiene $\sim$224,100 fuentes. El análisis aplica cortes de densidad de flujo en 180 $\mu$Jy y 400 $\mu$Jy para crear muestras conservadoras. Se construyen mapas de peso para tener en cuenta el ruido observacional y las variaciones en el recuento de fuentes.
   • DES: Utiliza la muestra de galaxias limitada por magnitud (MagLim) dividida en seis bins tomográficos de redshift ($0.20 \le z < 1.05$). Los redshifts fotométricos se derivan mediante el algoritmo Directional Neighbourhood Fitting (DNF).
   • Planck: Incluye el mapa de convergencia de lente de la CMB de Planck ($\kappa$) para proporcionar una validación independiente del sesgo de radio y la distribución de redshift.

2. Estimación del Espectro de Potencia:

   • Los espectros de potencia angular ($C_\ell$) se computan utilizando un estimador de pseudo-$C_\ell$ (NaMaster) para manejar la cobertura incompleta del cielo.
   • El análisis mide el espectro de potencia auto-espectro de EMU-PS1 ($C^{EMU}_\ell$), los auto-espectros de los bins de DES ($C^{DES}_\ell$), los espectros de potencia cruzados entre EMU-PS1 y los bins de DES ($C^{EMU \times DES}_\ell$), y el espectro de potencia cruzado entre EMU-PS1 y el lensing de Planck ($C^{EMU \times \kappa}_\ell$).
   • Se estima una matriz de covarianza utilizando 1,000 realizaciones ficticias (mocks) generadas mediante GLASS, incorporando las huellas del sondeo y el ruido de Poisson.

3. Modelado e Inferencia:

   • Distribución de Redshift ($n(z)$): Se prueban dos modelos. Primero, la distribución de redshift del Tiered Radio Extragalactic Continuum Simulation (TRECS), que incluye contribuciones de Núcleos Activos de Galaxias (AGNs) y Galaxias con Formación Estelar (SFGs). Segundo, un modelo paramétrico (adaptado de Nakoneczny et al. 2024) con tres parámetros libres ($z_0$, $r$, $\gamma$) diseñado para motivar físicamente el pico y la cola de la distribución.
   • Sesgo de Galaxia: Se adoptan modelos de sesgo lineal e independiente de la escala. Para el $n(z)$ paramétrico, se utiliza un modelo de amplitud constante para las SFGs y una evolución polinómica para los AGNs (o un sesgo efectivo combinado).
   • Estimación de Parámetros: Se emplea un marco de máxima verosimilitud (usando PyMC) para restringir los parámetros de sesgo y los parámetros de la distribución de redshift mediante el ajuste del vector de datos conjunto de los espectros de potencia cruzados.

4. Validación:

   • El flujo de trabajo de análisis se valida utilizando catálogos ficticios donde el $n(z)$ y el sesgo reales son conocidos. Esto confirma la capacidad del método para recuperar la ubicación del pico de la distribución de redshift.
   • Se evalúa el impacto de la incompleta residual en el catálogo de radio, mostrando que, si bien los auto-espectros de potencia son sensibles a la incompleta, los espectros de potencia cruzados utilizados para la inferencia siguen siendo robustos.

Resultados Clave

• Significancia de Detección: La correlación cruzada entre EMU-PS1 y los seis bins tomográficos de DES se detecta con significancias que oscilan entre $3.5\sigma$ y $5.5\sigma$ en el rango de multipolos $50 < \ell < 250$.
• Validación del Sesgo de DES: El sesgo de galaxia estimado para los bins de DES utilizando únicamente mediciones de agrupamiento es consistente con los resultados previos de DES (por ejemplo, Abbott et al. 2022) dentro de $1\sigma$, validando el flujo de trabajo de análisis.
• Recuperación de la Distribución de Redshift:
  • Modelo TRECS: El $n(z)$ simulado por TRECS, cuando se combina con un modelo de mezcla para el sesgo de galaxia, se ajusta a los espectros de potencia cruzados, al auto-espectro de potencia de EMU-PS1 y a la correlación cruzada con el lensing de la CMB. Los valores de $\chi^2$ totales son 7.78 (180 $\mu$Jy) y 3.15 (400 $\mu$Jy) para 20 grados de libertad.
  • Modelo Paramétrico: El modelo paramétrico (Ec. 11) también proporciona un buen ajuste a los espectros de potencia cruzados sin requerir conocimiento previo de la mezcla AGN/SFG. La distribución de mejor ajuste tiene un pico en $z \sim 0.85$. Este pico es superior al predicho por TRECS y los campos profundos de LOFAR, pero es consistente con los sondeos COSMOS 3GHz y MIGHTEE.
• Evolución del Sesgo: El sesgo de galaxia recuperado para la evolución de EMU-PS1 es consistente tanto con las predicciones de TRECS como con estimaciones independientes de la literatura, incluso cuando se utiliza el modelo paramétrico que no separa explícitamente las poblaciones de fuentes.
• Sistemáticas: El estudio señala una discrepancia donde el modelo paramétrico no logra ajustar el auto-espectro de potencia de EMU-PS1 en el corte de 180 $\mu$Jy, lo que sugiere sistemáticas residuales en los datos de radio (potencialmente incompleta) que no afectan los resultados de la correlación cruzada.

Significancia y Reivindicaciones
Los autores afirman que este trabajo demuestra la viabilidad de estimar las distribuciones de redshift y el sesgo de galaxias para sondeos de continuo de radio de gran área directamente a partir de correlaciones cruzadas con sondeos ópticos, sin depender de los redshifts de fuentes individuales o únicamente de simulaciones.

• Validación del Método: La técnica de correlación cruzada recupera con éxito la forma de la distribución de redshift y valida el modelo de sesgo, proporcionando un enfoque complementario al emparejamiento óptico directo.
• Utilidad Cosmológica: Los resultados validan el uso del modelo $n(z)$ paramétrico y la simulación TRECS para futuros análisis cosmológicos. Esto es significativo para los próximos sondeos de gran área como el sondeo EMU completo y el Observatorio del Square Kilometre Array (SKAO), donde la medición directa de redshift para todas las fuentes será imposible.
• Robustez: Se demuestra que la metodología es robusta frente a las sistemáticas internas y la incompleta residual en el catálogo de radio cuando se utilizan espectros de potencia cruzados, aunque los autores advierten que el pequeño área del piloto ($\sim$270 deg$^2$) limita la precisión de las restricciones en comparación con lo que será alcanzable con el sondeo principal.

El artículo concluye que la correlación cruzada con sondeos ópticos es una herramienta viable y necesaria para la era de la cosmología de precisión de la radioastronomía, permitiendo la mitigación de las incertidumbres en la distribución de redshift que actualmente limitan el retorno científico de los datos de continuo de radio.