Diffuse gamma-ray emissions around the stellar cluster Berkeley 59

Resumo Técnico: Emissões difusas de raios gama ao redor do aglomerado estelar Berkeley 59

Problema e Motivação
A origem dos raios cósmicos (RC) galácticos permanece uma questão fundamental em aberto na astrofísica. Embora os remanescentes de supernova (SNRs) tenham sido há muito considerados os principais locais de aceleração, evidências crescentes sugerem que aglomerados estelares, particularmente aqueles que hospedam estrelas massivas do tipo OB, servem como fontes significativas de energia para raios cósmicos de muito alta energia (VHE) e ultra-alta energia (UHE). Esses aglomerados geram ventos estelares poderosos e explosões de supernova que criam superbolhas, potencialmente acelerando partículas a altas energias. No entanto, distinguir entre origens leptônicas e hadrônicas para emissões de alta energia nesses ambientes é desafiador, e a distribuição espacial dos raios cósmicos dentro de ambientes de aglomerados é frequentemente simplificada em modelos de difusão. Este estudo investiga o jovem aglomerado estelar Berkeley 59 (Ber 59) para determinar se ele atua como uma fonte de emissão difusa de raios gama e para restringir os mecanismos físicos de aceleração e propagação de raios cósmicos dentro de tal aglomerado.

Metodologia
Os autores conduziram uma análise abrangente utilizando 16 anos de dados do Fermi-LAT (Pass 8, coletados de agosto de 2008 a outubro de 2024) dentro de uma região de interesse de $15^\circ \times 15^\circ$ centrada na fonte 4FGL J0002.1+6721, que está deslocada por apenas 4,8 minutos de arco do Ber 59.

1. Análise de Raios Gama: A análise utilizou os pacotes Fermipy e Fermi Science Tools. Uma análise de verossimilhança em bins foi realizada na faixa de energia de 300 MeV a 300 GeV. A equipe testou vários modelos espaciais (fonte pontual, Gaussiana radial e disco radial) contra os modelos de fundo (emissão difusa galáctica gll_iem_v07 e emissão extragaláctica isotrópica P8_R3_V3). Eles empregaram o Estatístico de Teste (TS) e o Critério de Informação de Akaike (AIC) para avaliar a significância das extensões espaciais e a preferência por modelos espectrais específicos (Lei de Potência, LogParabola e PLSuperExpCutoff).
2. Caracterização do Conteúdo Gasoso: Para investigar a origem hadrônica, os autores mapearam a distribuição de hidrogênio molecular ($H_2$), hidrogênio atômico neutro ($HI$) e hidrogênio ionizado ($HII$) ao redor do Ber 59.
   • A densidade de coluna de $H_2$ foi derivada de dados de CO (Dame et al. 2001) usando um fator de conversão $X_{CO} = 2.0 \times 10^{20} \text{ cm}^{-2} \text{ K}^{-1} \text{ km}^{-1} \text{ s}$.
   • Dados de $HI$ foram obtidos do levantamento de todo o céu HI4PI.
   • A densidade de coluna de $HII$ foi derivada do mapa de free-free do Planck (353 GHz).
3. Modelagem da Injeção de RC: Os autores modelaram a distribuição espectral de energia (SED) usando o pacote NAIMA. Eles testaram tanto cenários hadrônicos (interações próton-próton levando ao decaimento de píons neutros) quanto cenários leptônicos (espalhamento Compton inverso e bremsstrahlung). Além disso, eles derivaram a relação entre a eficiência de aceleração de raios cósmicos ($\eta_{CR}$) e o coeficiente de difusão ($D$) com base na taxa de injeção de energia cinética dos ventos estelares.

Principais Resultados

• Extensão Espacial: A análise revelou uma emissão estendida significativa em GeV ao redor do Berkeley 59. A emissão é melhor modelada por um disco radial com um raio de $1.02^\circ \pm 0.17^\circ$, resultando em uma significância de extensão de $10.6\sigma$ (TS de extensão = 97 para o modelo de disco). Essa extensão é robusta mesmo ao levar em conta fontes próximas como 4FGL J0005.6+6746c.
• Propriedades Espectrais: O espectro de raios gama é bem descrito por uma simples Lei de Potência (PL) com um índice de fótons de $\Gamma = 2.88 \pm 0.06$ e uma luminosidade de $L_\gamma = 1.15 \times 10^{33} \text{ erg/s}$. A falta de curvatura espectral significativa ($TS_{LP} = 0$, $TS_{PLSC} = 1.6$) favorece um modelo PL sobre o modelo LogParabola usado no catálogo 4FGL para a fonte próxima.
• Massa e Densidade Gasosas: A massa total de gás dentro da região estendida (derivada de $H_2$ e $HII$) é de aproximadamente $2.89 \times 10^2 M_\odot$, com uma densidade numérica média de prótons de $99.41 \text{ cm}^{-3}$. A massa estelar total do aglomerado é estimada em $\sim 5.4 \times 10^4 M_\odot$.
• Eficiência de Aceleração e Difusão: Assumindo uma origem hadrônica, o melhor ajuste para o índice de prótons é $\alpha_p = 2.91 \pm 0.05$, com uma energia total de prótons acima de 1 GeV de $W_p = 1.51 \times 10^{49} \text{ erg}$. Os autores derivaram um coeficiente de difusão $D$ na faixa de $8.3 \times 10^{26}$ a $3.3 \times 10^{27} \text{ cm}^2 \text{ s}^{-1}$, correspondendo a uma eficiência de aceleração $\eta_{CR}$ de 0,1 a 0,4. Este coeficiente de difusão é menor que a média galáctica, sugerindo um confinamento aprimorado de RC dentro do ambiente do aglomerado.
• Contribuição Leptônica: Embora um cenário leptônico (espalhamento Compton inverso e bremsstrahlung) tenha sido testado e tenha produzido um melhor ajuste para o índice de elétrons de $\alpha_e = 2.19 \pm 0.04$, os autores observam que essa origem não pode ser definitivamente descartada, embora o cenário hadrônico seja consistente com a densidade de gás observada e o índice espectral.

Significância e Alegações
O artigo afirma que o Berkeley 59 é uma fonte de emissão estendida de raios gama originada de raios cósmicos acelerados por ventos do aglomerado interagindo com o gás circundante. O estudo fornece uma restrição detalhada sobre a eficiência de aceleração de RC e o coeficiente de difusão em um ambiente de aglomerado estelar jovem. O coeficiente de difusão derivado, sendo menor que a média galáctica, apoia a hipótese de que campos magnéticos turbulentos dentro do aglomerado (impulsionados por ventos estelares e supernovas) aprimoram o confinamento de RC, facilitando a aceleração eficiente de partículas.

Os autores concluem que a energia cinética total injetada pelos ventos estelares ao longo da vida útil de $\sim 2$ Myr do aglomerado ($\sim 6.07 \times 10^{49} \text{ erg}$) é suficiente para explicar as emissões de raios gama observadas. Este trabalho contribui para o crescente corpo de evidências de que aglomerados estelares são contribuintes significativos para o orçamento de raios cósmicos galácticos, oferecendo um laboratório único para estudar a inter-relação entre formação estelar, estrutura do meio interestelar e física de partículas de alta energia. Os autores reconhecem que incertezas na densidade efetiva de gás e na eficiência de aceleração exigem refinamento adicional por meio de mapeamento de gás mais detalhado em observações futuras.