Diffuse gamma-ray emission in the vicinity of open cluster Berkeley 87

Resumen Técnico: Emisión Difusa de Rayos Gamma en las Cercanías del Cúmulo Abierto Berkeley 87

Planteamiento del Problema
Los cúmulos estelares masivos jóvenes (CEMJ) son sitios potenciales para la aceleración eficiente de partículas mediante aceleración difusiva por choques (DSA) impulsada por vientos estelares potentes de estrellas OB y Wolf-Rayet (WR). Si bien se han detectado emisiones de rayos gamma extendidas alrededor de varios CEMJ, la naturaleza específica de los mecanismos de emisión y el papel de los blancos de gas ambiental en el cúmulo Berkeley 87 aún deben caracterizarse completamente. Berkeley 87, un cúmulo masivo que alberga al menos 15 estrellas OB y la estrella WR 142, está incrustado en un medio interestelar denso. Observaciones anteriores, incluidos los catálogos de HEGRA y EGRET, arrojaron detecciones no concluyentes o nulas de emisión en el rango de TeV, mientras que se sugirieron correlaciones espaciales con fuentes EGRET no identificadas. Este estudio tiene como objetivo investigar la emisión difusa de rayos gamma alrededor de Berkeley 87 utilizando datos a largo plazo de Fermi-LAT, determinar sus propiedades espaciales y espectrales, y evaluar los escenarios hadrónicos versus leptónicos para su origen, abordando también la contribución del púlsar cercano PSR J2021+3651.

Metodología
Los autores analizaron 18 años de datos de Fermi-LAT dirigidos a la región de Berkeley 87 (Ascensión Recta = 305.396°, Declinación = 37.4206°).

• Selección de Datos: El análisis utilizó las funciones de respuesta instrumental P8_R3_v6, seleccionando fotones con energías >10 GeV para el análisis espacial con el fin de mitigar problemas de resolución angular a bajas energías y el fondo de albedo terrestre. La región de interés (ROI) fue un cuadrado de 15° × 15°.
• Puerta de Púlsar: Para aislar la emisión del cúmulo, los autores realizaron un análisis de puerta de púlsar en el cercano PSR J2021+3651 (asociado a la PWN Dragonfly). Utilizando efemérides del catálogo 3PC y Tempo2, definieron fases fuera del pulso para excluir la emisión magnetosférica directa del púlsar.
• Modelado Espacial: Se realizó un análisis de verosimilitud binned (gtlike). La significancia de la emisión extendida se probó utilizando la Estadística de Prueba para extensión ($TS_{ext}$) frente a modelos de fuente puntual. Se compararon modelos gaussianos radiales y modelos de disco radial.
• Análisis de Blancos de Gas: Para investigar el escenario hadrónico, los autores combinaron datos de CO (Dame et al. 2001) para hidrógeno molecular ($H_2$), la Encuesta HI Effelsberg-Bonn (EBHIS) para hidrógeno atómico neutro (HI) y datos de Planck a 353 GHz para hidrógeno ionizado (HII). Estos se utilizaron para construir mapas de densidad de columna de gas y calcular la masa total de gas y la densidad de protones.
• Modelado Espectral y Físico: La distribución espectral de energía (SED) se derivó utilizando un modelo de ley de potencias. Los autores emplearon el paquete NAIMA para ajustar tanto escenarios hadrónicos (interacción protón-protón) como leptónicos (Compton inverso y frenado) a la SED observada.

Resultados Clave

• Detección de Emisión Extendida: El estudio reporta una detección significativa de emisión extendida de rayos gamma centrada en Berkeley 87. La emisión se modela mejor con un perfil gaussiano radial con una extensión angular ($R_{68}$) de $0.36 \pm 0.07^\circ$ ($TS_{ext} = 27.6$).
• Propiedades Espectrales: El espectro de rayos gamma sigue una distribución de ley de potencias con un índice de fotones de $2.68 \pm 0.34$ en el rango de energía de 300 MeV a 500 GeV.
• Ambiente de Gas: La masa total de gas dentro de la región se calcula en $\sim 21.62 \times 10^2 M_\odot$, comprendiendo contribuciones de $H_2$, HI y HII. La densidad promedio de protones se estima en $\sim 368 \text{ cm}^{-3}$, proporcionando un blanco denso para interacciones de rayos cósmicos.
• Presupuesto Energético: Se estima que la tasa de inyección de energía cinética de los vientos estelares del cúmulo supera $1.0 \times 10^{36} \text{ erg s}^{-1}$. Asumiendo una eficiencia de aceleración de rayos cósmicos ($\eta_{CR}$) entre $10^{-2.0}$ y $10^{-2.5}$, el cúmulo puede suministrar energía suficiente para alimentar la emisión de rayos gamma observada.
• Mecanismo de Origen: Se favorece el escenario hadrónico. El índice de protones de mejor ajuste para el modelo hadrónico es $\alpha_p = 2.82 \pm 0.05$, con una energía total de protones ($W_p$) de $5.2 \times 10^{49} \text{ erg}$ por encima de 5 GeV. Este índice es consistente con la DSA en choques de terminación de vientos estelares y zonas de colisión viento-viento.
• Contribución del Púlsar: El análisis de PSR J2021+3651 en fases fuera del pulso arrojó una baja significancia de extensión ($TS_{ext} = 2.1$), y las pruebas de curvatura espectral sugirieron que la emisión fuera del pulso se origina en la magnetosfera del púlsar. Los autores concluyen que es poco probable que la emisión extendida alrededor de Berkeley 87 esté dominada por este púlsar o su PWN.

Significado y Afirmaciones
El artículo afirma que la emisión extendida de rayos gamma alrededor de Berkeley 87 es una detección robusta que respalda la hipótesis de que los CEMJ actúan como aceleradores de partículas eficientes. Los autores afirman que la emisión se explica mejor mediante procesos hadrónicos, donde los rayos cósmicos acelerados dentro del cúmulo interactúan con el gas ambiental denso (colisiones protón-protón que producen desintegración de piones). Esta conclusión se ve respaldada por la alta densidad de gas, la potente energía cinética del viento estelar del cúmulo y el índice espectral consistente con la aceleración por choques. El estudio refuerza la visión de que los vientos estelares colectivos en cúmulos masivos pueden impulsar una aceleración y transporte significativos de rayos cósmicos, contribuyendo a la población galáctica de rayos cósmicos. El trabajo también aclara la relación espacial entre el cúmulo y el púlsar cercano PSR J2021+3651, distinguiendo la emisión difusa del cúmulo de la actividad magnetosférica del púlsar.