Constraining high-energy neutrinos from tidal disruption events with IceCube high-energy starting events

Resumen Técnico: Restricción de Neutrinos de Alta Energía de Eventos de Disrupción de Marea con IceCube HESE

Planteamiento del Problema
Los Eventos de Disrupción de Marea (TDE, por sus siglas en inglés), donde las estrellas son destrozadas por agujeros negros supermasivos, se proponen como candidatos para fuentes de neutrinos de alta energía. Los modelos teóricos sugieren que la producción de neutrinos puede ocurrir en diversas regiones, incluyendo chorros relativistas (tanto exitosos como sofocados), discos de acreción, coronas y flujos de salida. Si bien estudios previos identificaron posibles coincidencias entre TDEs específicos (por ejemplo, AT2019dsg) y eventos de neutrinos, mejoras recientes en la reconstrucción direccional de IceCube han demostrado que las direcciones reconstruidas de los neutrinos no se alinean con las localizaciones de los TDEs. En consecuencia, se requiere un enfoque estadístico para determinar si los TDEs, como población, contribuyen al flujo difuso de neutrinos de alta energía. El desafío principal radica en las bajas estadísticas tanto del catálogo de TDEs observado como del conjunto de datos de neutrinos de alta energía, lo que requiere un análisis de apilamiento (stacking) robusto para restringir los parámetros físicos que gobiernan la producción de neutrinos.

Metodología
Los autores realizan un análisis de máxima verosimilitud no agrupado (unbinned maximum likelihood analysis) para buscar correlaciones espaciales y temporales entre el conjunto de datos de 12.5 años de HESE de IceCube y un catálogo de 89 TDEs.

• Catálogos de Datos:

  • TDEs: Se seleccionó una muestra de 89 TDEs de un catálogo de 134 eventos [73]. Los criterios de selección requirieron la disponibilidad de fotometría de series temporales para determinar el tiempo de brillo máximo ($t_{peak}$) utilizando un algoritmo de Bloques Bayesianos. La muestra abarca varios desplazamientos al rojo e incluye eventos detectados en múltiples bandas ópticas (ZTF, ATLAS, ASASSN, etc.).
  • Neutrinos: El análisis utiliza 164 eventos HESE del conjunto de datos de 12.5 años de IceCube (junio 2010 – agosto 2022). El conjunto de datos incluye direcciones de llegada reconstruidas modeladas mediante una distribución Fisher-Bingham de 8 parámetros (FB8) para dar cuenta de las incertidumbres de localización anisotrópicas.

• Análisis de Verosimilitud:

  • El análisis emplea una función de verosimilitud no agrupada (Ec. 3.1) que apila las asociaciones entre cada TDE y todos los eventos de neutrinos en la vecindad temporal y espacial.
  • PDF de Señal ($S_i$): Modelada como el producto de componentes temporales y espaciales. El componente temporal asume una distribución Gaussiana centrada en un retraso de $\lambda_T = 100$ días respecto al pico óptico del TDE, con una dispersión de $\sigma_T = 10$ días, motivada por modelos de lanzamiento de chorros retardados (por ejemplo, estados de Disco Magnéticamente Arrestado). El componente espacial utiliza la densidad de probabilidad direccional FB8 proporcionada por IceCube.
  • PDF de Fondo ($B_i$): Asume una distribución uniforme en el tiempo (dentro de una ventana de $\sim$50 días) y una distribución uniforme sobre el cielo.
  • Estadístico de Prueba (TS): Definido como $TS = 2 \log(L(ns)/L(ns=0))$, donde $ns$ es el parámetro de intensidad de la señal. El $TS$ observado ($TS_{obs}$) se compara contra una distribución de $TS_{max}$ generada a partir de $10^5$ catálogos simulados con posiciones y tiempos de TDE generados aleatoriamente para establecer la significancia.

• Límites Superiores e Inyección de Señal:

  • Para restringir los parámetros físicos, los autores generan catálogos de TDE sintéticos y realizan una inyección de señal artificial. Varían dos parámetros clave: la fracción de TDEs que albergan chorros intrínsecos ($f_{jet}$) y la energía de rayos cósmicos isotrópicamente equivalente ($E_{CR}$).
  • Los catálogos simulados se construyen seleccionando una submuestra de TDEs con chorro basándose en $f_{jet}$, calculando los rendimientos de neutrinos esperados según $E_{CR}$, e inyectando pares de señal en el conjunto de datos con una correlación espacial perfecta y retrasos temporalmente correlacionados.

Resultos Clave

• Búsqueda de Correlación: El análisis de los datos observados arroja un estadístico de prueba máximo de $TS_{max} = 0$ en $ns = 0$. Los datos observados son consistentes con la hipótesis de solo fondo, mostrando ninguna correlación espacial o temporal significativa entre los 89 TDEs y los 164 neutrinos HESE.
• Restricciones sobre Parámetros Físicos: Utilizando el resultado nulo, los autores establecen límites superiores en el espacio de parámetros de $f_{jet}$ y $E_{CR}$.
  • Incluso con estadísticas limitadas, el análisis restringe la región donde $f_{jet} \gtrsim 0.6$ y $E_{CR} \gtrsim 3 \times 10^{53}$ erg con un nivel de confianza superior al 90% (C.L.).
  • Para el caso específico donde todos los TDEs son emisores de neutrinos ($f_{jet} = 1$), la energía de rayos cósmicos isotrópicamente equivalente se restringe a $E_{CR} < 2 \times 10^{53}$ erg al 90% C.L.
• Implicaciones Teóricas:
  • Los límites derivados restringen la eficiencia ($\epsilon_{CR}$) de la conversión de la energía de unión gravitacional en rayos cósmicos. Para una estrella estándar de $1 M_\odot$ y un agujero negro de $10^8 M_\odot$, los límites implican $\epsilon_{CR} \lesssim 0.5$ (para $f_{jet}=1$) o $\epsilon_{CR} \lesssim 0.75 f_b$ (para $f_{jet} < 0.6$).
  • Los resultados sugieren que, si bien los chorros con emisión de radio (que constituyen una fracción minúscula de los TDE) son exitosos, es poco probable que dominen el flujo difuso; los escenarios que involucran vientos "ocultos", chorros sofocados o emisión isotrópica de TDEs sin chorro siguen siendo viables pero están restringidos por el presupuesto energético.
  • El límite superior calculado para el flujo de neutrinos difuso de todo el cielo proveniente de los TDE es actualmente débil ($\sim 10^{-7}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$) comparado con el flujo observado de IceCube, pero los autores señalan que esto mejorará con datos futuros.

Significancia y Reivindicaciones
El artículo afirma ser la primera búsqueda de apilamiento que correlaciona el conjunto de datos HESE de IceCube con un catálogo específico de TDEs observados. La significancia principal radica en establecer que los datos actuales no apoyan una fuerte correlación, descartando así modelos donde una gran fracción de los TDEs ($f_{jet} \gtrsim 0.6$) produce rayos cósmicos de energía extremadamente alta ($E_{CR} \gtrsim 3 \times 10^{53}$ erg).

Los autores enfatizan que su metodología proporciona un marco para establecer restricciones estrictas sobre los parámetros físicos de los TDE a medida que las estadísticas mejoren. Expresan explícitamente que con la llegada del Observatorio Vera C. Rubin (LSST), que se espera detecte $\sim$1000 TDEs por año, y futuros observatorios de neutrinos (IceCube-Gen2, KM3NeT, GRAND, etc.), los límites actuales pueden restringirse significativamente. El artículo concluye modestamente que, aunque las estadísticas actuales son limitadas, el análisis demuestra el potencial para restringir el presupuesto energético y las fracciones de chorro de los TDE, ofreciendo un camino para comprender su papel en el cielo de neutrinos de alta energía.