Multi-Messenger Modeling of Low-Luminosity Gamma-Ray Bursts

Resumen Técnico: Modelado Multimensajero de Brotes de Rayos Gamma de Baja Luminosidad

Planteamiento del Problema
El origen astrofísico del flujo difuso de neutrinos astrofísicos de alta energía (rango TeV–PeV) permanece en gran medida sin identificar. Aunque los Núcleos Galácticos Activos (AGN), como NGC 1068 y TXS 0506+056, han sido identificados como fuentes, estos no pueden explicar la totalidad del flujo difuso observado, particularmente en energías $\gtrsim 0.1$ PeV. Los brotes de rayos gamma (GRBs) han sido propuestos durante mucho tiempo como candidatos a fuentes; sin embargo, las estrictas restricciones de IceCube sobre los típicos GRBs de alta luminosidad (HL) limitan su contribución al flujo difuso a $\lesssim 1\%$. Los GRBs de baja luminosidad (LL GRBs), que son de 10 a 1000 veces más abundantes que los HL GRBs y presentan eyecciones trans-relativistas, siguen siendo candidatos prometedores, pero se encuentran insuficientemente caracterizados. Sus orígenes físicos son diversos, surgiendo potencialmente de escenarios de ruptura de choque (SBO), chorros obstruidos (choked jets) o chorros emergentes ligeramente relativistas; sin embargo, sus contribuciones específicas al flujo de neutrinos y a los factores de carga de rayos cósmicos (CR) no han sido cuantificadas exhaustivamente utilizando datos multimensajeros.

Metodología
Los autores presentan un enfoque de modelado leptohadrónico unificado y novedoso aplicado a una muestra de siete LL GRBs (GRB 050826, GRB 060218, GRB 100316D, GRB 120422A, GRB 171205A, GRB 190829A y GRB 201015A) observados por Swift-BAT.

1. Modelo Físico: El estudio utiliza el Astrophysical Multimessenger Emission Simulator (AMES) para modelar la emisión prompt. El marco acomoda tanto escenarios de choque interno (IS) como de ruptura de choque (SBO) al tratar el radio de disipación ($R$) y el factor de Lorentz de movimiento global ($\Gamma$) como parámetros libres que abarcan rangos amplios ($R \sim 10^{11.5}–10^{16}$ cm, $\Gamma \sim 2–20$).
2. Física de Partículas: El modelo asume la aceleración de choque difusiva (DSA) para protones y electrones. Resuelve las ecuaciones de transporte para las densidades numéricas de partículas y la redistribución de energía, contabilizando el enfriamiento por sincrotrón, la dispersión Compton inversa (IC), la producción de pares por Bethe-Heitler (BH) e interacciones fotomesónicas ($p\gamma$). Se realiza un seguimiento de los piones y muones secundarios para derivar los espectros de neutrinos y fotones de alta energía.
3. Análisis de Datos: Los autores realizan un ajuste de Monte Carlo de Cadenas de Markov (MCMC) utilizando el paquete emcee. El modelo está restringido por:
   • Espectros de rayos X de Swift-BAT durante la duración $T_{90}$.
   • Límites superiores de rayos gamma de Fermi-LAT (donde estén disponibles) durante la fase prompt.
   • Redshifts observados y luminosidades isotrópicas equivalentes ($L_{\gamma, \text{iso}}$).
   • Los límites superiores de neutrinos de IceCube se mencionan pero no se utilizan como restricciones, ya que exceden las predicciones del modelo.
4. Análisis Estadístico: Las distribuciones posteriores para los parámetros físicos (fracción de energía magnética $\xi_B$, fracción de energía electrónica $\xi_e$, factor de carga de CR $\xi_p$, radio de emisión $R$, factor de Lorentz del chorro $\Gamma_j$, índice espectral electrónico $s_e$ y partición de energía electrónica $\epsilon_e$) son derivadas. Se aplica el Análisis de Componentes Principales (PCA) a las muestras de MCMC para visualizar la diversidad y el agrupamiento de los siete eventos en el espacio de parámetros.

Contribuciones Clave y Resultados

• Restricciones de Parámetros: El estudio logra restringir el factor de carga de CR ($\xi_p$) a través de la muestra, encontrando medianas posteriores en el rango $\xi_p \sim 0.2–1.6$. Los intervalos de confianza de 1$\sigma$ abarcan $\log(\xi_p) \sim -2.03$ a $2.00$.
• Diversidad de Mecanismos de Emisión: El análisis de PCA revela que los LL GRBs no son una población homogénea.
  • Candidatos a SBO: GRB 060218 y GRB 100316D (los eventos de menor luminosidad, $L_{\gamma, \text{iso}} \sim 10^{46-47}$ erg s$^{-1}$) exhiben los factores de carga de CR más altos ($\xi_p \sim 1.3–1.6$) y son consistentes con escenarios de SBO. Muestran una partición de energía equilibrada entre campos magnéticos, electrones y protones.
  • Eventos Dominados por IC: GRB 050826 y GRB 171205A muestran espectros de rayos X duros dominados por la dispersión IC, lo que implica campos magnéticos más débiles y menores requerimientos de carga de CR.
  • Eventos Dominados por Sincrotrón: GRB 190829A y GRB 201015A se caracterizan por componentes de sincrotrón fuertes y altas fracciones de energía magnética ($\xi_B$), lo que suprime las contribuciones de cascada y limita la eficiencia de producción de neutrinos.
• Predicciones de Neutrinos: El modelo predice la fluencia de neutrinos esperada para cada evento. Si bien las instalaciones actuales (IceCube, KM3NeT) son improbables de detectar LL GRBs individuales (eventos esperados $\lesssim 0.01$ por brote), los observatorios de próxima generación (IceCube-Gen2, KM3NeT, P-ONE, HUNT, TRIDENT) podrían alcanzar $\sim 0.1$ eventos por cada LL GRB cercano.
• Potencial de Apilamiento (Stacking): Los autores sugieren que los análisis de apilamiento, particularmente aprovechando muestras identificadas por misiones de campo amplio como el Einstein Probe (EP), podrían producir señales estadísticamente significativas.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera caracterización cuantitativa y exhaustiva de las propiedades de los LL GRBs utilizando un modelo leptohadrónico unificado restringido por datos multilongitud de onda. Las reivindicaciones clave son:

1. Los LL GRBs como Fuentes de Neutrinos Viables: A pesar de las restricciones de IceCube sobre los HL GRBs, los LL GRBs siguen siendo fuentes prometedoras de neutrinos astrofísicos de alta energía, particularmente aquellos con altos factores de carga de CR y duraciones largas (por ejemplo, GRB 060218, GRB 100316D).
2. Diversidad Física: El estudio demuestra que los LL GRBs abarcan un continuo de mecanismos físicos (desde chorros obstrucidos/SBO hasta chorros emergentes exitosos), lo que requiere un modelado individual en lugar de un promedio de población único.
3. Marco Metodológico: El uso combinado de ajuste MCMC y PCA ofrece un método robusto para distinguir subpoblaciones e identificar los candidatos más prometedores para la producción de neutrinos.
4. Perspectiva Futura: Los autores enfatizan que las instalaciones de rayos X y MeV de próxima generación son críticas para identificar más LL GRBs y refinar las restricciones de los parámetros. Abogan por búsquedas en tiempo real de eventos coincidentes de LL GRB y neutrinos, señalando que incluso las no detecciones proporcionan restricciones valiosas sobre los parámetros físicos. El enfoque de modelado se presenta como ampliamente aplicable a otros transitorios astrofísicos multimensajeros.