Are Single-Zone Emission models Sufficient to Explain GRB 220426A and GRB 230812B?

Resumen Técnico: ¿Son suficientes los modelos de emisión de zona única para explicar GRB 220426A y GRB 230812B?

Planteamiento del Problema
El mecanismo de emisión prompt de los estallidos de rayos gamma (GRB) sigue siendo un tema controvertido en la astrofísica de altas energías. Aunque el modelado empírico suele utilizar la función de Band, el origen físico de la radiación se debate entre escenarios térmicos (bola de fuego) y no térmicos (sincrotrón). Estudios anteriores han dependido en gran medida del análisis espectral integrado en el tiempo y del índice espectral de baja energía ($\alpha$) para distinguir estos orígenes, invocando a menudo los límites de la "línea de muerte" (LOD) para la radiación de sincrotrón. Sin embargo, los análisis integrados en el tiempo pueden enmascarar una evolución espectral crucial, y las bajas estadísticas de fotones a altas energías a menudo permiten que múltiples modelos reproduzcan el mismo espectro, lo que genera ambigüedades en la interpretación. Además, trabajos teóricos recientes sugieren que el ancho espectral a media altura ($W$) en la representación $E^2N(E)$ es un discriminador crítico, pero su evolución temporal no ha sido estudiada sistemáticamente para estallidos brillantes específicos para probar los modelos de emisión de zona única.

Metodología
Los autores realizaron un análisis espectral resuelto en el tiempo de dos GRB excepcionalmente brillantes detectados por Fermi/GBM: GRB 220426A y GRB 230812B.

• Selección de Datos: Se analizaron datos de eventos con marca de tiempo (TTE) de los detectores NaI y BGO más brillantes. Se generaron intervalos de tiempo utilizando el algoritmo Bayesian Block, asegurando una significancia de detección superior a $8\sigma$ y un mínimo de 800 conteos netos para evitar degeneraciones de parámetros.
• Reestructuración del Modelo: Para sondear directamente el ancho espectral, los autores reestructuraron la función de Band estándar y la función de Ley de Potencia con Corte (CPL). En estas formas reestructuradas, el ancho espectral a media altura ($W$) se trató como un parámetro libre junto con el índice de baja energía ($\alpha$) y la energía de pico ($E_p$) o la energía de corte ($E_c$).
• Procedimiento de Ajuste: El ajuste espectral se realizó utilizando el marco de Máxima Verosimilitud Multi-Misión (3ML) con una verosimilitud Poisson-Gaussiana y el algoritmo DYNESTY para inferencia bayesiana. Se seleccionó el modelo con el Criterio de Información Bayesiano (BIC) más bajo y posteriores bien restringidas como el mejor ajuste para cada intervalo de tiempo.
• Análisis Comparativo: La evolución temporal de los parámetros de mejor ajuste ($W$, $E_p$, $\alpha$) se comparó con las predicciones teóricas para:
  1. Modelos de Bola de Fuego Térmica: Específicamente, bolas de fuego en expansión en regímenes dominados por materia o por radiación, donde se espera que $W$ disminuya a medida que la bola de fuego se expande y el espectro se acerca a una función de Planck instantánea.
  2. Modelos de Sincrotrón: Específicamente, emisión de sincrotrón de una distribución de electrones Maxwelliana relativista, donde se espera que el enfriamiento por sincrotrón cause que el ancho espectral $W$ disminuya con el tiempo a medida que se agotan los electrones de alta energía.

Contribuciones Clave

• Aplicación de una Métrica Novel: Este trabajo introduce la evolución temporal del ancho espectral ($W$) como una herramienta de diagnóstico principal para la emisión prompt de GRB, superando los ajustes espectales estáticos e integrados en el tiempo.
• Funciones Espectrales Reestructuradas: El artículo proporciona un marco matemático para ajustar funciones de Band y CPL donde $W$ es un parámetro libre directo, permitiendo una estimación estadística más robusta de la amplitud espectral sin depender de datos desenrollados en post-procesamiento.
• Restricciones Resueltas en el Tiempo: Al analizar 18 intervalos de tiempo para GRB 220426A y 16 para GRB 230812B, el estudio proporciona restricciones de alta resolución sobre cómo evolucionan los parámetros espectrales durante la duración del estallido, revelando comportamientos que contradicen los modelos estándar de zona única.

Resultados

• Evolución del Ancho Espectral ($W$): Contrario a las predicciones de los modelos de enfriamiento tanto de bola de fuego térmica como de sincrotrón (que predicen una disminución de $W$), el $W$ observado para ambos GRB 220426A y GRB 230812B aumenta con el tiempo.
  • Para GRB 230812B, $W$ aumenta de $\sim0.9$ a $\sim1.6$.
  • Para GRB 220426A, $W$ permanece relativamente constante alrededor de $\sim0.7$ pero muestra un aumento moderado durante la fase de decaimiento del segundo pico.
  • En ambos casos, los espectros permanecen más amplios que una función de Planck ($W \approx 0.54$) durante toda la duración.
• Energía de Pico ($E_p$) e Índice ($\alpha$):
  • $E_p$ disminuye monótonamente, lo que es consistente con algunos escenarios de enfriamiento por sincrotrón, pero también potencialmente consistente con la expansión de la bola de fuego.
  • El índice de baja energía $\alpha$ presenta una contradicción: para GRB 220426A, $\alpha$ permanece más duro que la LOD de sincrotrón durante todo el estallido, desfavoreciendo un origen puramente de sincrotrón. Para GRB 230812B, $\alpha$ se ablanda y cruza la LOD, sugiriendo una transición, sin embargo, la tasa de decaimiento de $E_p$ no concuerda con interpretaciones térmicas puras.
• Inconsistencia del Modelo: Ningún modelo de zona única (térmico puro o no térmico puro) puede explicar simultáneamente el aumento de $W$, la evolución específica de $E_p$ y el comportamiento de $\alpha$. El aumento de $W$ es particularmente problemático para los escenarios estándar de zona única.

Significado y Afirmaciones
El artículo concluye que los modelos de emisión de zona única son insuficientes para explicar la fase prompt de GRB 220426A y GRB 230812B. El aumento observado en el ancho espectral ($W$) a lo largo del tiempo sugiere fuertemente que la emisión proviene de múltiples zonas de emisión cuyas contribuciones relativas cambian dinámicamente durante el estallido.

Los autores proponen que esta complejidad podría surgir de:

1. Choques Internos: Formados por un motor central activado esporádicamente, donde el ensanchamiento del espectro sugiere un aumento en el número de choques, mientras que la disminución de $E_p$ indica un declive en la actividad del motor central.
2. Reconexión Turbulenta (ICMART): Un entorno de plasma turbulento que induce múltiples eventos de reconexión magnética, creando zonas de emisión distintas. La radiación acumulada de estas zonas amplía el espectro observado.

El estudio enfatiza que el análisis resuelto en el tiempo es esencial para descubrir estos componentes ocultos, ya que las observaciones integradas en el tiempo a menudo enmascaran la verdadera naturaleza del mecanismo de emisión. Los autores sugieren que futuros estudios de polarización (por ejemplo, a través de COSI o POLAR2) podrían elucidar aún más la reorientación del campo magnético asociada con estos escenarios de múltiples zonas.