Probing Cosmic Magnetism with Rotation Measure-Squared-Galaxy Cross-Correlations

Resumen Técnico: Sondeo del Magnetismo Cósmico mediante Correlaciones Cruzadas de la Medida de Rotación al Cuadrado y Galaxias

Planteamiento del Problema
Los campos magnéticos cósmicos permean el universo, desde los planetas hasta el medio intergaláctico (IGM), sin embargo, su origen, evolución y estructura espacial siguen siendo poco comprendidos. Una herramienta observacional primaria es la rotación de Faraday, que mide la Medida de Rotación (RM) de fuentes de fondo. No obstante, la RM proporciona solo una restricción integral a lo largo de la línea de visión, mezclando contribuciones de la fuente, la Vía Láctea y la estructura a gran escala (LSS) interviniente. Además, las correlaciones cruzadas directas entre la RM y la distribución de galaxias son problemáticas porque el campo de la RM puede ser positivo o negativo, lo que provoca cancelaciones que diluyen la señal. Los intentos previos para eludir esto mediante la correlación del valor absoluto de la RM ($|RM|$) con las galaxias introducen no diferenciabilidad en cero y, crucialmente, sufren de un sesgo de ruido, donde el valor esperado de $|RM + n|$ (donde $n$ es el ruido) difiere de $|RM|$.

Metodología
Los autores proponen un nuevo estadístico: la correlación cruzada entre el campo de la Medida de Rotación al cuadrado ($RM^2$) y la sobredensidad proyectada de galaxias ($g$), denotada como $\langle RM^2 \times g \rangle$. Este enfoque es análogo al estimador de "campos proyectados" utilizado para el efecto de Sunyaev-Zel'dovich cinético (kSZ) ($\langle kSZ^2 \times g \rangle$).

1. Marco Teórico: Los autores derivan que el espectro de potencia angular de $RM^2$ y $g$ se relaciona con un bispectro que involucra dos copias del campo magnético ponderado por la densidad de electrones ($\tilde{B}_{\perp,\parallel}$) y una copia de la sobredensidad de galaxias. En el espacio real, la función de correlación $w_{RM^2,g}(r_\perp)$ se muestra como una proyección a lo largo de la línea de visión de la correlación entre el campo magnético al cuadrado y las galaxias.
2. Aproximación Analítica: Se desarrolla un modelo analítico aproximado, que expresa la correlación como el producto del espectro de potencia cruzada entre electrones y galaxias y una integral sobre el espectro de potencia del campo magnético ponderado por la densidad de electrones. Este modelo sugiere que la forma de la señal está gobernada por el agrupamiento de electrones, mientras que su amplitud rastrea la fuerza del campo magnético proyectado.
3. Simulaciones: Los autores utilizan las simulaciones magnetohidrodinámicas (MHD) cosmológicas Illustris-TNG (específicamente la ejecución TNG300-3) para construir mapas de RM ficticios. Generan mapas de $RM^2$ integrando la densidad de electrones y el campo magnético a lo largo de las líneas de visión, y luego aplican un método de apilamiento (stacking) alrededor de halos de materia oscura para estimar la función de correlación cruzada $w_{RM^2,g}(r_\perp; z)$.
4. Análisis de Ruido: El artículo compara explícitamente el estimador $\langle RM^2 \times g \rangle$ frente al estimador $\langle |RM| \times g \rangle$. Demuestra matemáticamente y mediante simulaciones que $\langle RM^2 \times g \rangle$ es inmune al sesgo de ruido (asumiendo que el ruido no está correlacionado con las galaxias), mientras que $\langle |RM| \times g \rangle$ sufre de una supresión sistemática debido a la naturaleza no lineal de la operación de valor absoluto sobre datos con ruido.

Resultados Clave

• Detección de Señal en Simulaciones: Utilizando Illustris-TNG, los autores detectan un fuerte exceso de señal de $RM^2$ alrededor de los halos. La señal aumenta aproximadamente tres órdenes de magnitud desde $z \sim 3$ hasta $z \sim 0$.
• Evolución con el Redshift: La evolución es impulsada por la amplificación de los campos magnéticos en regiones sobredensas mediante dinamos turbulentos y el transporte de plasma magnetizado por salidas (outflows) de AGN y galácticas. La señal está dominada por regiones altamente magnetizadas dentro de los halos de materia oscura, en lugar del IGM difuso.
• Validación del Modelo Analítico: El modelo analítico aproximado (Eq. 3.3) reproduce con éxito la dependencia de la escala y la evolución del redshift de las señales simuladas, confirmando que la amplitud está determinada principalmente por la fuerza del campo magnético ponderado por la densidad de electrones.
• Dependencia de la Escala de Suavizado: La señal es altamente sensible a la escala de suavizado de las mediciones de RM. Debido a que el espectro de potencia del campo magnético contiene fuertes fluctuaciones a pequeña escala, un suavizado más grueso suprime significativamente la amplitud de la señal.
• Dependencia de la Resolución: Las simulaciones de mayor resolución (TNG300-1) producen amplitudes de señal más bajas en comparación con las ejecuciones de menor resolución (TNG300-3). Esto se atribuye a una mejor resolución de las inversiones del campo magnético, que causan cancelaciones en la proyección a lo largo de la línea de visión, una característica que no es capturada por las ejecuciones de menor resolución.
• Pronósticos: Los pronósticos para los actuales y futuros sondeos (p. ej., SKA, POSSUM, NVSS) sugieren que $\langle RM^2 \times g \rangle$ es detectable con alta significancia.
  • Un sondeo tipo SKA con $\sim 10^7$ fuentes de RM y $\sim 10^6$ galaxias por bin de redshift podría detectar la señal en un amplio rango de escalas ($0.5 - 50$ Mpc/h) y redshifts.
  • Incluso los conjuntos de datos actuales (p. ej., NVSS combinado con WISE) pueden permitir detecciones iniciales, particularmente a bajos redshifts, aunque con relaciones señal-ruido más bajas.
  • El estadístico $\langle RM^2 \times g \rangle$ produce una relación señal-ruido ligeramente superior que $\langle |RM| \times g \rangle$ y evita el sesgo de ruido de este último.

Significancia y Reivindicaciones
El artículo afirma proporcionar un método robusto y tomográfico para sondear la evolución de los campos magnéticos cósmicos a gran escala. Al utilizar $\langle RM^2 \times g \rangle$, el método:

1. Aísla la Evolución del Redshift: Extrae limpiamente las contribuciones de los campos magnéticos cerca de los redshifts de las galaxias trazadoras, filtrando las contribuciones intrínsecas de las fuentes y de la Vía Láctea, que actúan únicamente como varianza de ruido.
2. Evita Sesgos Sistemáticos: A diferencia de $\langle |RM| \times g \rangle$, el estimador al cuadrado es insesgado por el ruido de medición, lo que lo hace más fiable para el análisis cuantitativo.
3. Conecta con la Física Fundamental: El estadístico se vincula directamente con los bispectros subyacentes de los campos magnéticos y las densidades de electrones, ofreciendo un camino más claro para interpretar las propiedades físicas de la red cósmica magnetizada.
4. Habilita Ciencia Futura: Los autores argumentan que con los próximos datos del Square Kilometer Array (SKA) y grandes sondeos de galaxias (p. ej., DESI), esta técnica permitirá realizar un censo de la densidad de energía magnética cósmica a través del tiempo cósmico, abordando preguntas abiertas sobre la siembra y amplificación de campos magnéticos.

El artículo mantiene la modestia respecto a la interpretación del origen de la señal, señalando que, si bien la señal rastrea la evolución de la densidad de energía magnética, actualmente está dominada por el plasma magnetizado dentro de los halos. Se requiere trabajo futuro para refinar los modelos de galaxias, igualar con precisión las escalas de suavizado observacionales y tener en cuenta la varianza de la muestra y los términos de covarianza no gaussianos.