Exploring selection biases in FRB dispersion-galaxy cross-correlations with magnetohydrodynamical simulations

Título: Explorando sesgos de selección en las correlaciones cruzadas de dispersión de FRB y galaxias con simulaciones magnetohidrodinámicas

1. El Problema

Las ráfagas rápidas de radio (FRB) son transitorios de radio extragalácticos cuya medida de dispersión (DM) es una herramienta poderosa para trazar la distribución de plasma extragaláctico y, por ende, la materia bariónica del universo. Una técnica prometedora para medir la distribución de electrones libres es la correlación cruzada angular entre las DM de los FRB y las posiciones de las galaxias, lo que permite extraer el espectro de potencia de electrones-galaxias ($P_{eg}$).

Sin embargo, la precisión de esta medición depende críticamente de la comprensión de los efectos de selección observacional. La probabilidad de detectar un FRB y obtener su redshift (y por tanto su galaxia anfitriona) no es uniforme; depende de propiedades intrínsecas como la DM, el brillo de la galaxia anfitriona y efectos de propagación como la dispersión (scattering). Si estos sesgos no se cuantifican y corrigen, pueden introducir errores sistemáticos significativos en la estimación de la distribución de bariones, subestimando o sobreestimando la señal de correlación.

2. Metodología

Los autores utilizaron simulaciones de magnetohidrodinámica (MHD) a gran escala del proyecto IllustrisTNG (específicamente la caja de simulación TNG300-1 de 300 Mpc) para realizar un seguimiento de rayos (ray tracing) y estudiar el impacto de diversos sesgos de selección.

• Simulación y Muestreo:

  • Se simularon 3000 FRB en un rango de redshift $0.3 \le z \le 0.4$, distribuidos aleatoriamente en galaxias ponderadas por su tasa de formación estelar (SFR).
  • Se utilizaron galaxias de fondo en el rango $0.2 \le z \le 0.3$ para la correlación cruzada.
  • La densidad de electrones se mapeó en una cuadrícula de 500 ckpc/h para permitir el cálculo eficiente de la DM a lo largo de líneas de visión arbitrarias.
  • Se definieron 48 parches de cielo independientes ($0.1 \times 0.1$ rad) dentro de "regiones buenas" para evitar artefactos geométricos causados por las condiciones de frontera periódicas de la simulación.

• Estimador de Espectro de Potencia:

  • Se empleó un Estimador Cuadrático Óptimo (OQE) para calcular el espectro de potencia cruzada angular ($C^{Dg}_\ell$) entre el campo de DM y el campo de densidad de galaxias. Este método es robusto frente a ventanas de observación incompletas (parches de cielo no cubiertos uniformemente).
  • Se compararon los resultados simulados con el espectro de potencia "verdadero" obtenido directamente de la densidad de electrones y el catálogo de halos en la simulación.

• Pruebas de Sesgos de Selección:
  Se introdujeron artificialmente varios efectos de selección para medir su impacto en la señal de correlación:

  1. Propiedades de la galaxia anfitriona: Selección basada en la magnitud aparente (galaxias tenues) y la contribución de la DM de la propia galaxia anfitriona (inyectada mediante una distribución log-normal).
  2. Efectos de propagación (DM): Funciones de selección dependientes de la DM, modeladas basándose en las limitaciones del sondeo CHIME/FRB (función log-normal y cortes abruptos).
  3. Dispersión (Scattering): Modelado de la probabilidad de no detección debido a la dispersión temporal causada por halos intermedios.

3. Contribuciones Clave

• Marco de Simulación: Desarrollo de un pipeline de ray tracing sobre datos de IllustrisTNG que permite inyectar efectos de selección observacional realistas y medir su impacto directo en estadísticas de dos puntos (correlaciones cruzadas).
• Validación de Robustez: Demostración cuantitativa de que las correlaciones cruzadas FRB-galaxia son inherentemente robustas frente a la mayoría de los sesgos locales (propiedades de la anfitriona, DM local), a diferencia de las estadísticas de "un punto" (distribución de DM).
• Identificación de Riesgos Críticos: Identificación de que los cortes severos en la DM (específicamente aquellos que eliminan los FRB más dispersos) son la principal fuente de sesgo sistemático, capaz de alterar drásticamente la amplitud de la señal cosmológica.

4. Resultados Principales

• Robustez frente a propiedades de la anfitriona:
  • La selección de galaxias anfitrionas basada en su brillo (magnitud aparente) no sesga significativamente la señal de correlación.
  • La contribución de la DM de la propia galaxia anfitriona (un término de ruido) aumenta la varianza de la medición pero no introduce un sesgo en la amplitud del espectro de potencia cruzada, ya que la DM local no está correlacionada con las galaxias de primer plano.
• Robustez frente a efectos de CHIME/FRB suaves:
  • Una función de selección dependiente de la DM suave (similar a la sensibilidad reducida de CHIME a altas DMs, escalada al rango de redshift estudiado) afecta mínimamente la señal de correlación.
• Sesgo Crítico por Cortes de DM:
  • Se encontró que un efecto de selección que elimina abruptamente el 10% de los FRB con mayor DM (un corte de tipo "step-function") reduce la amplitud de la señal de correlación cruzada en más del 50% en escalas angulares de $\sim 0.1^\circ$ (correspondientes a escalas físicas de $\sim$ Mpc).
  • Esto se debe a que los FRB con mayor DM son los que trazan las mayores sobredensidades de electrones (cúmulos de galaxias); eliminarlos elimina la mayor parte de la señal cosmológica de interés.
• Efectos de Dispersión (Scattering):
  • Para el rango de redshift local estudiado ($z < 0.4$), la dispersión causada por halos intermedios no afecta significativamente la señal de correlación. Sin embargo, se advierte que este efecto podría volverse relevante a redshifts más altos.

5. Significado e Implicaciones

Este trabajo es fundamental para la próxima generación de estudios cosmológicos con FRB (con telescopios como CHORD, DSA-2000 y BURSTT).

• Validación del Método: Confirma que la correlación cruzada FRB-galaxia es una técnica viable y robusta para mapear la materia bariónica, superando muchas de las limitaciones de las estadísticas de un punto.
• Guía para Diseño de Observatorios: Destaca la necesidad crítica de diseñar instrumentos y pipelines de detección que minimicen los cortes abruptos en la DM. Para evitar sesgos del 50% en la medición de la estructura a gran escala, los futuros sondeos deben ser capaces de detectar y localizar FRB con altas medidas de dispersión, o bien deben modelar y corregir rigurosamente estos sesgos utilizando simulaciones como las presentadas aquí.
• Planificación de Análisis: Proporciona un marco metodológico para que los futuros análisis de datos reales incluyan simulaciones de selección específicas para sus instrumentos, asegurando que las mediciones de la distribución de bariones no estén contaminadas por efectos instrumentales.