Tracing Large-scale Structure with the MeerKLASS On-the-Fly Survey: Angular Clustering of Radio Sources at 816 MHz

Resumen Técnico: Rastreo de la Estructura a Gran Escala con el Sondeo MeerKLASS On-the-Fly

Problema y Motivación
Los sondeos de continuo de radio ofrecen una ventana única, libre de la extinción por polvo, hacia la estructura a gran escala (LSS, por sus siglas en inglés) del Universo, rastreando el campo de densidad de materia a través de núcleos galácticos activos (AGN) radio-fuertes y galaxias con formación estelar (SFG). Sin embargo, las mediciones robustas de agrupamiento angular a partir de catálogos de radio no son triviales debido al agrupamiento espacialmente variable de la completitud impulsado por el ruido dependiente de la posición, la morfología compleja de las fuentes (donde una sola galaxia física se resuelve en múltiples componentes del catálogo) y la falta general de información directa de desplazamiento al rojo (redshift). Mediciones previas utilizando sondeos como FIRST y NVSS establecieron el agrupamiento de poblaciones brillantes dominadas por AGN, pero persiste la necesidad de mediciones de agrupamiento de alta fidelidad a frecuencias intermedias de $\sim$GHz utilizando interferómetros modernos con variaciones de profundidad bien caracterizadas y solapamiento con programas contemporáneos de espectroscopía óptica. Este artículo aborda estos desafíos presentando el primer análisis de agrupamiento angular del primer lanzamiento de datos (DR1) del sondeo on-the-fly (OTF) de la banda UHF del MeerKLASS (MeerKAT Large Area Synoptic Survey).

Metodología
El análisis utiliza el conjunto de datos MeerKLASS UHF DR1, que cubre $\sim$800 deg$^2$ dentro de la huella de DESI a una frecuencia de referencia de 816 MHz. El sondeo emplea una estrategia OTF de elevación constante con escaneos de ascenso y descenso entrelazados, lo que resulta en un patrón de profundidad no uniforme que debe ser modelado explícitamente.

1. Definición de la Muestra: Se construyó una muestra de agrupamiento conservadora a partir de los 95,483 fuentes únicas en el catálogo de fuentes (SRL) del DR1. La selección fiduciaria aplica:

   • Un umbral de densidad de flujo de $S_{816} > 2$ mJy.
   • Una máscara controlada por la profundidad que restringe el análisis a regiones con ruido RMS local $\sigma_{local} \le 250$ $\mu$Jy beam$^{-1}$.
   • Un criterio de "solo compactas" ($S\_Code = 'S'$ y $N\_Gaus = 1$) para mitigar el exceso de potencia a pequeña escala causado por la fragmentación de galaxias de radio de múltiples componentes.
   • Un corte conservador de borde en Ascensión Recta (RA) ($RA \le 184.07^\circ$) para asegurar la uniformidad entre los datos y los aleatorios.
     Esto resultó en un catálogo de agrupamiento final de 39,483 fuentes.

2. Máscara y Aleatorios: Se construyó una máscara de sondeo binaria a partir de la huella geométrica y el mapa de profundidad de RMS. Los catálogos de aleatorios se generaron uniformemente dentro de esta máscara para trazar la función de selección angular, validando la uniformidad de la relación de densidad superficial de datos a aleatorios.

3. Medición de Correlación: La función de correlación de dos puntos angular, $w(\theta)$, se midió utilizando el estimador de Landy–Szalay sobre 20 bins logarítmicamente espaciados ($0.02^\circ \le \theta \le 10^\circ$). Las incertidumbres y las covarianzas entre bins se estimaron utilizando una técnica de remuestreo de jackknife espacial con 30 parches. Se aplicó una corrección de restricción integral (IC) para tener en cuenta la huella finita del sondeo.

4. Modelado e Inferencia de Sesgo:

   • Ajuste Fenomenológico: La señal se caracterizó mediante un modelo de ley de potencia $w(\theta) \approx A(\theta/1^\circ)^{1-\gamma}$ con una pendiente fija $\gamma=1.8$.
   • Interpretación Física: El agrupamiento medido se interpretó ajustando plantillas de materia proyectada $\Lambda$CDM ($w_{DM}(\theta)$) computadas con CAMB y la aproximación de Limber. Dado que no se dispone de desplazamientos al rojo individuales, el análisis acotó la principal incertidumbre sistemática (la distribución de desplazamiento al rojo $N(z)$) utilizando dos priors derivados de simulaciones semiempíricas de T-RECS: un prior dominado por AGN y un prior TOTAL (AGN + SFG). Se infirió un sesgo efectivo a gran escala ($b_{eff}$) escalando la plantilla a los datos.

Resultados Clave

• Detección: Se detectó una señal de agrupamiento positiva en el rango de $0.02^\circ \lesssim \theta \lesssim 10^\circ$. La medición es estable ante variaciones en la máscara de profundidad y el umbral de flujo en escalas intermedias y grandes.
• Amplitud de la Ley de Potencia: El ajuste de la señal a escala intermedia ($0.112^\circ \le \theta \le 1.36^\circ$) con una pendiente fija $\gamma=1.8$ arroja una amplitud de $A(1^\circ) = (1.434 \pm 0.475) \times 10^{-3}$, correspondiente a $\log_{10} A = -2.843^{+0.124}_{-0.175}$.
• Sesgo Efectivo: El sesgo efectivo inferido depende significativamente del prior de la distribución de desplazamiento al rojo asumido:
  • Prior de AGN: $b_{eff} = 1.998 \pm 0.350$.
  • Prior TOTAL: $b_{eff} = 1.530 \pm 0.265$.
    La dispersión entre estos valores ($\Delta b_{eff} \approx 0.47$) se identifica como la incertidumbre de modelado dominante, derivada de la mezcla de poblaciones desconocida en la muestra limitada por flujo.
• Longitud de Correlación: La inversión de la amplitud de la ley de potencia mediante el método de Limber produce longitudes de correlación espacial de $r_0 = 6.18 \pm 1.13$ $h^{-1}$Mpc (prior de AGN) y $r_0 = 5.59 \pm 1.02$ $h^{-1}$Mpc (prior TOTAL).
• Sistemáticas: Relajar el criterio de compacidad resultó en un aumento pronunciado en $w(\theta)$ a escalas pequeñas ($\theta < 0.1^\circ$), confirmando la necesidad de la selección de solo compactas para mediciones robustas de agrupamiento a gran escala.

Significancia y Reivindicaciones
El artículo establece a MeerKLASS UHF DR1 como un nuevo conjunto de datos de amplio área y frecuencia intermedia para el agrupamiento de continuo de radio. Los autores afirman que sus resultados demuestran la viabilidad de la obtención de imágenes de continuo interferométrico OTF para mediciones de precisión de la estructura a gran escala. El estudio destaca que, si bien los resultados actuales están limitados por la incertidumbre en la distribución de desplazamiento al rojo ($N(z)$), el solapamiento con la huella de DESI proporciona un camino claro a seguir. Los autores declaran que futuros lanzamientos, incorporando espectroscopía óptica/IR superpuesta para una mejor calibración del desplazamiento al rojo, permitirán mediciones de agrupamiento divididas por población y mediciones de evolución del sesgo con una incertidumbre de modelado sustancialmente reducida. El presente trabajo sirve como un paso fundacional, proporcionando una medición conservadora y bien controlada que valida la estrategia de sondeo y los productos de datos para aplicaciones cosmológicas.