Probing Evolution of Long Gamma-Ray Burst Properties through Their Cosmic Formation History

Resumen Técnico: Sondeo de la Evolución de las Propiedades de los GRB Largos a Través de su Historia de Formación Cósmica

Enunciado del Problema
La astrofísica de los progenitores de los Estallidos de Rayos Gamma Largos (LGRB) y la posible evolución cosmológica de sus propiedades siguen siendo temas de significativo debate. Mientras que estudios previos sugieren que la densidad de tasa de LGRB (LGRB-RD) rastrea la densidad de tasa de formación estelar (SFRD) a altos corrimientos al rojo ($z$), atribuyendo esto a la evolución de la metalicidad, la relación a bajos corrimientos al rojo es controvertida. Algunos estudios indican un "repunte" o exceso en la LGRB-RD a bajo $z$ en comparación con la SFRD, potencialmente causado por una clase distinta de GRB de baja luminosidad (LLGRB), contaminación de la muestra o diferentes subclases de progenitores. Además, la literatura existente ofrece visiones conflictivas sobre si los LGRB son trazadores imparciales de la formación estelar o trazadores sesgados influenciados por la metalicidad y la función de masa inicial (IMF). Este estudio tiene como objetivo resolver estas discrepancias investigando la LGRB-RD y su mapeo a la SFRD utilizando una muestra homogénea de LGRB con características específicas de meseta en rayos X.

Metodología
Los autores utilizan una muestra de 242 LGRB detectados entre marzo de 2005 y diciembre de 2023, caracterizados por mesetas en rayos X en sus curvas de luz de remanente observadas por el Observatorio Neil Gehrels Swift. La selección de esta submuestra específica está motivada por la homogeneidad de los datos y la reducción de la dispersión en las correlaciones en comparación con los estudios de emisión de impulso.

El análisis procede a través de los siguientes pasos:

1. Límites de Luminosidad y Flujo: Se calculan las luminosidades observadas al final de la meseta ($L_{a,obs}$) y las luminosidades observables mínimas ($L_{min,obs}$). Se aplica un corte de límite de flujo ($\log F_{lim} = -11.9$), eliminando 17 estallidos para mitigar el sesgo de Malmquist.
2. Corrección No Paramétrica: Para abordar la evolución intrínseca de la luminosidad y los sesgos de selección, el estudio emplea el método no paramétrico de Efron & Petrosian (EP). Esta técnica prueba la independencia de la luminosidad ($L$) y el corrimiento al rojo ($1+z$) utilizando la estadística $\tau$ de Kendall.
3. Determinación del Factor de Evolución: El análisis determina un factor de evolución de luminosidad, $(1+z)^k$, donde $k = 2.1^{+0.5}_{-0.7}$. Este factor se utiliza para corregir la densidad de tasa acumulada ($\sigma$) y la función de luminosidad ($\Phi$), produciendo distribuciones corregidas ($\sigma_{corr}$ y $\Phi_{corr}$).
4. Cálculo de la Densidad de Tasa: La LGRB-RD diferencial ($\rho(z)$) se deriva de la distribución acumulada corregida. Los autores aplican factores de corrección ($F_1$ y $F_2$) para contabilizar los estallidos excluidos (debido a límites de flujo o conjuntos asociados nulos) y las limitaciones instrumentales (campo de visión de Swift-BAT, tiempo de actividad y años de operación) para estimar la verdadera LGRB-RD.
5. Comparación Teórica: La LGRB-RD corregida ($\rho_{corr}$) se compara contra una LGRB-RD teórica ($\rho_{theor}$) derivada de la SFRD (Madau & Dickinson 2014) utilizando un enfoque "similar a la ecuación de Drake". Este modelo teórico incorpora factores para la fracción de estrellas masivas, supernovas de colapso del núcleo, supernovas Tipo Ic, chorro (beaming) y metalicidad.

Contribuciones Clave

• Función de Luminosidad Actualizada: El estudio actualiza la función de luminosidad (LF) de LGRB para LGRB con mesetas en rayos X utilizando una muestra más reciente y mayor que trabajos anteriores (por ejemplo, Dainotti et al. 2013).
• Primera Análisis de RD para Muestra de Meseta en Rayos X: Este es el primer estudio que deriva la LGRB-RD específicamente para la muestra de meseta en rayos X y la compara directamente con modelos de SFRD.
• Restricciones Evolutivas: El artículo proporciona una restricción cuantitativa sobre el factor de evolución de luminosidad ($k \approx 2.1$) requerido para hacer independientes las distribuciones de luminosidad y corrimiento al rojo, corrigiendo el sesgo de Malmquist.
• Análisis Comparativo: El trabajo compara la $\rho_{corr}$ derivada tanto con modelos teóricos como con resultados observacionales previos (incluyendo una muestra de meseta óptica de Dainotti et al. 2024) para aislar efectos de selección y tendencias evolutivas.

Resultados

• Evolución de Luminosidad: El análisis confirma una evolución significativa de la luminosidad, siendo el factor de evolución de mejor ajuste $(1+z)^{2.1}$. Sin esta corrección, la naturaleza intrínseca de la LGRB-RD queda oscurecida.
• Comportamiento de la Densidad de Tasa: La LGRB-RD corregida ($\rho_{corr}$) exhibe una distribución relativamente plana a través de la mayoría de los corrimientos al rojo con un pico sutil alrededor de $(1+z) \approx 2.5$.
• Desajuste con la SFRD:
  • Exceso a Bajo-$z$: El estudio confirma un exceso de LGRB a bajos corrimientos al rojo ($z < 1$) en comparación con la SFRD, consistente con los hallazgos de Yu et al. (2015), Petrosian et al. (2015) y otros. Este exceso disminuye en $z \sim 0.1$, probablemente debido a estadísticas de números pequeños.
  • Discrepancia a Alto-$z$: A altos corrimientos al rojo, la $\rho_{corr}$ muestra un desajuste con la mayoría de los modelos de SFRD, aunque se alinea parcialmente con el modelo de Hopkins & Beacom (2006). El pico de la SFRD ocurre a un corrimiento al rojo ligeramente mayor que el pico de la LGRB-RD.
• Prueba de Hipótesis de Evolución: Los autores prueban tres escenarios para mapear $\rho_{theor}$ a $\rho_{corr}$:
  1. Sin Evolución: Descartado ya que las distribuciones son incompatibles.
  2. Evolución del Ángulo de Chorro: Un factor de $(1+z)^{1.6}$ aplanaría la pendiente de alto-$z$ de la curva teórica pero falla en explicar el comportamiento a bajo-$z$.
  3. Evolución Simple de Ley de Potencia: Un único factor de ley de potencia $(1+z)^\delta$ no puede coincidir simultáneamente con los extremos de bajo-$z$ y alto-$z$ de la densidad de tasa observada. Valores negativos de $\delta$ ajustan el extremo de bajo-$z$, mientras que valores positivos ajustan el extremo de alto-$z$, pero ninguno funciona para todo el rango.

Significado y Afirmaciones
El artículo concluye que la hipótesis de "sin evolución" es insuficiente para explicar la relación entre la formación de LGRB y la formación estelar. Además, los modelos evolutivos simples que involucran solo el ángulo de chorro o una evolución de ley de potencia simple de las propiedades son inadecuados para mapear la LGRB-RD a la SFRD a través de todo el rango de corrimientos al rojo.

Los autores afirman que coincidir con la LGRB-RD y la SFRD en su totalidad probablemente requiere la inclusión de múltiples propiedades evolutivas de los LGRB simultáneamente. Las discrepancias observadas sugieren que deben tenerse en cuenta los sesgos de selección (como la contaminación por LLGRB o GRB impulsados por fusiones mal clasificados como LGRB) y factores evolutivos complejos (como la dependencia de la metalicidad, IMFs evolutivas o la evolución del sistema progenitor). El estudio enfatiza que, aunque la $\rho_{corr}(z)$ derivada es una herramienta poderosa para sondear sistemas progenitores y poblaciones de agujeros negros de masa estelar, es necesaria una función de evolución de corrimiento al rojo más elaborada para comprender plenamente la historia de formación cósmica de los LGRB. Los autores sugieren modestamente que se requerirán datos futuros, potencialmente ayudados por aprendizaje automático para la inferencia de corrimientos al rojo, para abordar aún más estos problemas complejos.