Diffuse gamma-ray emissions around the stellar cluster Berkeley 59

Resumen Técnico: Emisiones difusas de rayos gamma alrededor del cúmulo estelar Berkeley 59

Problema y Motivación
El origen de los rayos cósmicos (RC) galácticos sigue siendo una pregunta abierta fundamental en astrofísica. Aunque los remanentes de supernova (RSN) han sido considerados durante mucho tiempo los sitios principales de aceleración, crecientes evidencias sugieren que los cúmulos estelares, particularmente aquellos que albergan estrellas masivas OB, actúan como fuentes significativas de energía para los RC de muy alta energía (VHE) y ultra alta energía (UHE). Estos cúmulos generan vientos estelares potentes y explosiones de supernova que crean superburbujas, acelerando potencialmente partículas a altas energías. Sin embargo, distinguir entre orígenes leptónicos y hadrónicos para las emisiones de alta energía en estos entornos es desafiante, y la distribución espacial de los RC dentro de los entornos de cúmulos a menudo se simplifica en exceso en los modelos de difusión. Este estudio investiga el cúmulo estelar joven Berkeley 59 (Ber 59) para determinar si actúa como una fuente de emisión difusa de rayos gamma y para restringir los mecanismos físicos de aceleración y propagación de RC dentro de dicho cúmulo.

Metodología
Los autores realizaron un análisis exhaustivo utilizando 16 años de datos de Fermi-LAT (Pass 8, recopilados desde agosto de 2008 hasta octubre de 2024) dentro de una región de interés de $15^\circ \times 15^\circ$ centrada en la fuente 4FGL J0002.1+6721, la cual está desplazada solo 4.8 minutos de arco de Ber 59.

1. Análisis de rayos gamma: El análisis utilizó los paquetes Fermipy y Fermi Science Tools. Se realizó un análisis de verosimilitud binned sobre el rango de energía de 300 MeV a 300 GeV. El equipo probó varios modelos espaciales (fuente puntual, gaussiana radial y disco radial) contra los modelos de fondo (emisión difusa galáctica gll_iem_v07 y emisión extragaláctica isotrópica P8_R3_V3). Emplearon la Estadística de Prueba (TS) y el Criterio de Información de Akaike (AIC) para evaluar la significancia de las extensiones espaciales y la preferencia por modelos espectrales específicos (Ley de Potencia, LogParábola y PLSuperExpCutoff).
2. Caracterización del contenido de gas: Para investigar el origen hadrónico, los autores mapearon la distribución del hidrógeno molecular ($H_2$), el hidrógeno atómico neutro ($HI$) y el hidrógeno ionizado ($HII$) alrededor de Ber 59.
   • La densidad de columna de $H_2$ se derivó de datos de CO (Dame et al. 2001) utilizando un factor de conversión $X_{CO} = 2.0 \times 10^{20} \text{ cm}^{-2} \text{ K}^{-1} \text{ km}^{-1} \text{ s}$.
   • Los datos de $HI$ se obtuvieron del sondeo de todo el cielo HI4PI.
   • La densidad de columna de $HII$ se derivó del mapa de radiación libre-libre de Planck (353 GHz).
3. Modelado de la inyección de RC: Los autores modelaron la distribución espectral de energía (SED) utilizando el paquete NAIMA. Probaron tanto escenarios hadrónicos (interacciones protón-protón que conducen al decaimiento de piones neutros) como escenarios leptónicos (dispersión Compton inversa y frenado). Además, derivaron la relación entre la eficiencia de aceleración de rayos cósmicos ($\eta_{CR}$) y el coeficiente de difusión ($D$) basada en la tasa de inyección de energía cinética de los vientos estelares.

Resultados Clave

• Extensión Espacial: El análisis reveló una emisión extendida significativa en GeV alrededor de Berkeley 59. La emisión se modela mejor con un disco radial de un radio de $1.02^\circ \pm 0.17^\circ$, arrojando una significancia de extensión de $10.6\sigma$ (TS de extensión = 97 para el modelo de disco). Esta extensión es robusta incluso al tener en cuenta fuentes cercanas como 4FGL J0005.6+6746c.
• Propiedades Espectrales: El espectro de rayos gamma se describe bien mediante una simple Ley de Potencia (PL) con un índice fotónico de $\Gamma = 2.88 \pm 0.06$ y una luminosidad de $L_\gamma = 1.15 \times 10^{33} \text{ erg/s}$. La falta de curvatura espectral significativa ($TS_{LP} = 0$, $TS_{PLSC} = 1.6$) favorece un modelo PL sobre el modelo LogParábola utilizado en el catálogo 4FGL para la fuente cercana.
• Masa y Densidad del Gas: La masa total de gas dentro de la región extendida (derivada de $H_2$ y $HII$) es aproximadamente $2.89 \times 10^2 M_\odot$, con una densidad numérica promedio de protones de $99.41 \text{ cm}^{-3}$. La masa estelar total del cúmulo se estima en $\sim 5.4 \times 10^4 M_\odot$.
• Eficiencia de Aceleración y Difusión: Asumiendo un origen hadrónico, el índice de protones de mejor ajuste es $\alpha_p = 2.91 \pm 0.05$, con una energía total de protones por encima de 1 GeV de $W_p = 1.51 \times 10^{49} \text{ erg}$. Los autores derivaron un coeficiente de difusión $D$ en el rango de $8.3 \times 10^{26}$ a $3.3 \times 10^{27} \text{ cm}^2 \text{ s}^{-1}$, correspondiente a una eficiencia de aceleración $\eta_{CR}$ de 0.1 a 0.4. Este coeficiente de difusión es menor que el promedio galáctico, lo que sugiere un confinamiento mejorado de RC dentro del entorno del cúmulo.
• Contribución Leptónica: Aunque se probó un escenario leptónico (dispersión Compton inversa y frenado) y arrojó un índice de electrones de mejor ajuste de $\alpha_e = 2.19 \pm 0.04$, los autores señalan que este origen no puede descartarse definitivamente, aunque el escenario hadrónico es consistente con la densidad de gas observada y el índice espectral.

Significancia y Afirmaciones
El artículo afirma que Berkeley 59 es una fuente de emisión extendida de rayos gamma originada por rayos cósmicos acelerados por vientos del cúmulo que interactúan con el gas circundante. El estudio proporciona una restricción detallada sobre la eficiencia de aceleración de RC y el coeficiente de difusión en un entorno de cúmulo estelar joven. El coeficiente de difusión derivado, al ser menor que el promedio galáctico, respalda la hipótesis de que los campos magnéticos turbulentos dentro del cúmulo (impulsados por vientos estelares y supernovas) mejoran el confinamiento de RC, facilitando una aceleración eficiente de partículas.

Los autores concluyen que la energía cinética total inyectada por los vientos estelares durante la vida de $\sim 2$ Myr del cúmulo ($\sim 6.07 \times 10^{49} \text{ erg}$) es suficiente para explicar las emisiones de rayos gamma observadas. Este trabajo contribuye al creciente cuerpo de evidencia de que los cúmulos estelares son contribuyentes significativos al presupuesto de rayos cósmicos galácticos, ofreciendo un laboratorio único para estudiar la interacción entre la formación estelar, la estructura del medio interestelar y la física de partículas de alta energía. Los autores reconocen que las incertidumbres en la densidad efectiva de gas y la eficiencia de aceleración requieren un refinamiento adicional mediante un mapeo de gas más detallado en observaciones futuras.