Tomographic constraints on the high-energy cosmic neutrino emission rate

Resumen Técnico: Restricciones Tomográficas sobre la Tasa de Emisión de Neutrinos Cósmicos de Alta Energía

Planteamiento del Problema
A pesar de la detección de un flujo difuso de neutrinos astrofísicos de alta energía por el Observatorio IceCube, los orígenes específicos de la mayoría de estos eventos permanecen sin identificar. Si bien se han asociado fuentes específicas con eventos de neutrinos individuales, estas representan solo una fracción menor del flujo total. La población dominante parece ser isotrópica, lo que deja indeterminadas las características de las poblaciones de fuentes. La identificación de estas fuentes es un desafío primordial en la astrofísica contemporánea. Un enfoque prometedor consiste en buscar anisotropías estadísticas en las direcciones de llegada de los neutrinos, bajo la hipótesis de que estos trazan la estructura a gran escala (LSS, por sus siglas en inglés) del Universo. Sin embargo, los intentos previos de correlacionar los mapas de neutrinos con sondeos de galaxias han arrojado resultados mixtos o no significativos, careciendo a menudo de la capacidad para restringir los parámetros físicos absolutos de la tasa de emisión.

Metodología
Este trabajo emplea un análisis de correlación cruzada entre el conjunto de datos de 10 años de trazas de muones de IceCube y mapas tomográficos de sobredensidad de galaxias para restringir la emisividad de neutrinos de alta energía ponderada por el sesgo (bias). La metodología procede de la siguiente manera:

1. Construcción de Datos:

   • Datos de Neutrinos: Se utilizaron 10 años de datos de trazas de muones de IceCube (403,529 eventos) en el rango de energía de $10^3$ a $10^6$ GeV. Se construyeron mapas tanto para la densidad numérica ($I_N$) como para la densidad de energía ($I_E$). Para mitigar la contaminación por neutrinos atmosféricos, se excluyeron los eventos con declinaciones $\delta < -5^\circ$, y se restaron los componentes monopolo y dipolo de los mapas.
   • Datos de Galaxias: Se utilizaron cuatro catálogos distintos de galaxias y cuásares para cubrir un amplio rango de desplazamiento al rojo ($z \lesssim 3$): 2MPZ ($z \sim 0.06$), WI$\times$SC ($z \sim 0.23$), DECaLS ($z \sim 0.50$) y Quaia ($z \sim 1.72$). Estos fueron procesados en mapas HEALPix ($N_{side}=256$) aplicando máscaras para eliminar la contaminación del plano galáctico.

2. Marco Teórico:

   • El análisis define la intensidad de los neutrinos en términos de la emisividad de neutrinos comóvil ($\dot{n}_\nu$) y la emisividad de energía de neutrinos comóvil ($\dot{\rho}_\nu$), ponderadas por el sesgo lineal de las fuentes ($b_\nu$).
   • Los espectros de potencia angular de correlación cruzada ($C_\ell^{NG}$) entre los mapas de neutrinos y las sobredensidades de galaxias se calculan utilizando la aproximación de Limber. Los espectros teóricos dependen del producto del sesgo de la galaxia ($b_g$), el sesgo del neutrino ($b_\nu$), la emisividad media y el espectro de potencia de la materia.
   • Los sesgos de las galaxias ($b_g$) se determinaron independientemente mediante un análisis de autocorrelación y se trataron como parámetros fijos en la verosimilitud de la correlación cruzada.

3. Análisis Estadístico:

   • Se empleó un enfoque de inferencia bayesiana utilizando muestreo anidado (UltraNest) para restringir los parámetros del modelo.
   • Se probaron tres modelos para la evolución con el desplazamiento al rojo de la tasa de emisión ponderada por el sesgo ($b_\nu \dot{n}_\nu$ y $b_\nu \dot{\rho}_\nu$):
     • Ley de Potencia: $\propto (1+z)^a$.
     • Modelo de Pico: Siguiendo la historia de la formación estelar (Madau & Dickinson 2014).
     • Tomográfico (No paramétrico): Restringiendo la tasa de emisión de forma independiente en los cuatro intervalos de desplazamiento al rojo medios de los catálogos de galaxias.
   • Se calcularon las relaciones señal-ruido (SNR) para evaluar la significancia de la detección.

Resultados Clave

• Significancia de la Correlación: El análisis no encontró una correlación significativa entre los mapas de neutrinos de IceCube y los trazadores de la estructura a gran escala. La desviación más fuerte de un modelo de señal cero fue una detección de $1.9\sigma$ (específicamente para la muestra WI$\times$SC en el mapa ponderado por energía). Todas las demás correlaciones cruzadas fueron consistentes con cero dentro de las incertidumbres estadísticas.
• Restricciones de Parámetros:
  • Modelos Paramétricos: Para el modelo de ley de potencia, la amplitud del parámetro $b_\nu \dot{n}_\nu$ se restringió a $5.8 \pm 3.6 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ (compatible con cero a $1.6\sigma$). Para el modelo de pico, el valor fue $4.0 \pm 2.1 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ (compatible con cero a $1.9\sigma$). El índice de la ley de potencia $a$ permaneció sin restringir dentro del rango de la distribución previa (prior).
  • Restricciones Tomográficas: En un enfoque independiente del modelo, se establecieron límites superiores en la tasa de emisión en cuatro intervalos de desplazamiento al rojo. Aunque los valores centrales fueron positivos, permanecieron consistentes con cero dentro de $2\sigma$ en todos los intervalos (por ejemplo, $4.7 \pm 3.5 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ en $z \approx 0.06$ y $310 \pm 250 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ en $z \approx 1.7$).
• Pruebas de Robustez: Los resultados fueron verificados frente a variaciones en los intervalos de energía, cortes de declinación ($\delta \ge 0^\circ$ y $\delta \ge 5^\circ$) y la sustracción de una plantilla de neutrinos atmosféricos. Estas pruebas confirmaron que la falta de detección significativa es robusta y no un artefacto de la contaminación atmosférica o de las elecciones del análisis.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona las primeras restricciones tomográficas sobre la emisividad del fondo difuso astrofísico de neutrinos de alta energía. A diferencia de estudios previos que se centraron en anisotropías relativas o mapas normalizados, este análisis preserva la escala de flujo absoluto, lo que permite la derivación directa de límites superiores sobre la emisividad física de neutrinos comóvil ($\dot{n}_\nu$) y la emisividad de energía ($\dot{\rho}_\nu$), aunque sea degenerada con el sesgo de la fuente desconocido $b_\nu$.

El artículo posiciona este análisis como una nueva sonda para probar modelos de fuentes de neutrinos extragalácticos. Si bien los datos actuales no producen una detección definitiva, la metodología establece un marco para establecer límites cuantitativos superiores en las tasas de emisión como función del desplazamiento al rojo. Los autores señalan que se requerirán conjuntos de datos futuros de los observatorios de neutrinos de próxima generación (IceCube-Gen2, KM3NeT) y sondeos de galaxias profundos (LSST, Euclid) para aumentar el poder de restricción y potencialmente lograr una detección de esta correlación.