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

Resumo Técnico: Emissão Difusa de Raios Gama nas Proximidades do Aglomerado Aberto Berkeley 87

Enunciado do Problema
Aglomerados estelares massivos jovens (YMSCs) são locais potenciais para aceleração eficiente de partículas via aceleração difusiva por choque (DSA) impulsionada por ventos estelares poderosos de estrelas OB e Wolf-Rayet (WR). Embora emissões de raios gama estendidas tenham sido detectadas ao redor de vários YMSCs, a natureza específica dos mecanismos de emissão e o papel dos alvos de gás ambiente no aglomerado Berkeley 87 permanecem a serem totalmente caracterizados. Berkeley 87, um aglomerado massivo que abriga pelo menos 15 estrelas OB e a estrela WR 142, está embutido em um meio interestelar denso. Observações anteriores, incluindo os catálogos HEGRA e EGRET, produziram detecções inconclusivas ou não detecções de emissão em TeV, enquanto correlações espaciais com fontes EGRET não identificadas foram sugeridas. Este estudo visa investigar a emissão difusa de raios gama ao redor de Berkeley 87 usando dados de longo prazo do Fermi-LAT, determinar suas propriedades espaciais e espectrais e avaliar os cenários hadrônicos versus leptônicos para sua origem, abordando também a contribuição da pulsar próxima PSR J2021+3651.

Metodologia
Os autores analisaram 18 anos de dados do Fermi-LAT direcionados à região de Berkeley 87 (Ascensão Reta = 305,396°, Declinação = 37,4206°).

• Seleção de Dados: A análise utilizou as funções de resposta instrumental P8_R3_v6, selecionando fótons com energias >10 GeV para análise espacial, a fim de mitigar problemas de resolução angular de baixa energia e o fundo de albedo terrestre. A região de interesse (ROI) foi um quadrado de 15° × 15°.
• Portaria da Pulsar: Para isolar a emissão do aglomerado, os autores realizaram uma análise de portaria de pulsar na pulsar próxima PSR J2021+3651 (associada ao PWN Dragonfly). Usando efemérides do catálogo 3PC e o Tempo2, definiram fases fora do pulso para excluir a emissão direta da magnetosfera da pulsar.
• Modelagem Espacial: Foi realizada uma análise de verossimilhança binned (gtlike). A significância da emissão estendida foi testada usando a Estatística de Teste para extensão ($TS_{ext}$) contra modelos de fonte pontual. Modelos gaussianos radiais e modelos de disco radial foram comparados.
• Análise de Alvos de Gás: Para investigar o cenário hadrônico, os autores combinaram dados de CO (Dame et al. 2001) para hidrogênio molecular ($H_2$), o Levantamento HI Effelsberg-Bonn (EBHIS) para hidrogênio atômico neutro (HI) e dados do Planck a 353 GHz para hidrogênio ionizado (HII). Estes foram usados para construir mapas de densidade de coluna de gás e calcular a massa total de gás e a densidade de prótons.
• Modelagem Espectral e Física: A distribuição de energia espectral (SED) foi derivada usando um modelo de lei de potência. Os autores empregaram o pacote NAIMA para ajustar tanto cenários hadrônicos (interação próton-próton) quanto leptônicos (Compton inverso e bremsstrahlung) à SED observada.

Principais Resultados

• Detecção de Emissão Estendida: O estudo relata uma detecção significativa de emissão estendida de raios gama centrada em Berkeley 87. A emissão é melhor modelada por um perfil gaussiano radial com uma extensão angular ($R_{68}$) de $0,36 \pm 0,07^\circ$ ($TS_{ext} = 27,6$).
• Propriedades Espectrais: O espectro de raios gama segue uma distribuição de lei de potência com um índice de fótons de $2,68 \pm 0,34$ na faixa de energia de 300 MeV a 500 GeV.
• Ambiente de Gás: A massa total de gás dentro da região é calculada em $\sim 21,62 \times 10^2 M_\odot$, compreendendo contribuições de $H_2$, HI e HII. A densidade média de prótons é estimada em $\sim 368 \text{ cm}^{-3}$, fornecendo um alvo denso para interações de raios cósmicos.
• Orçamento de Energia: A taxa de injeção de energia cinética dos ventos estelares do aglomerado é estimada em exceder $1,0 \times 10^{36} \text{ erg s}^{-1}$. Assumindo uma eficiência de aceleração de raios cósmicos (CR) ($\eta_{CR}$) entre $10^{-2,0}$ e $10^{-2,5}$, o aglomerado pode fornecer energia suficiente para alimentar a emissão de raios gama observada.
• Mecanismo de Origem: O cenário hadrônico é favorecido. O índice de prótons de melhor ajuste para o modelo hadrônico é $\alpha_p = 2,82 \pm 0,05$, com uma energia total de prótons ($W_p$) de $5,2 \times 10^{49} \text{ erg}$ acima de 5 GeV. Este índice é consistente com DSA em choques de terminação de ventos estelares e zonas de colisão vento-vento.
• Contribuição da Pulsar: A análise da PSR J2021+3651 em fases fora do pulso resultou em uma baixa significância de extensão ($TS_{ext} = 2,1$), e testes de curvatura espectral sugeriram que a emissão fora do pulso origina-se da magnetosfera da pulsar. Os autores concluem que a emissão estendida ao redor de Berkeley 87 dificilmente é dominada por esta pulsar ou seu PWN.

Significado e Alegações
O artigo alega que a emissão estendida de raios gama ao redor de Berkeley 87 é uma detecção robusta que apoia a hipótese de YMSCs atuando como aceleradores de partículas eficientes. Os autores afirmam que a emissão é melhor explicada por processos hadrônicos, onde raios cósmicos acelerados dentro do aglomerado interagem com o gás ambiente denso (colisões próton-próton produzindo decaimento de píons). Esta conclusão é apoiada pela alta densidade de gás, pela forte potência cinética do vento estelar do aglomerado e pelo índice espectral consistente com aceleração por choque. O estudo reforça a visão de que ventos estelares coletivos em aglomerados massivos podem impulsionar aceleração e transporte significativos de CR, contribuindo para a população de raios cósmicos galácticos. O trabalho também esclarece a relação espacial entre o aglomerado e a pulsar próxima PSR J2021+3651, distinguindo a emissão difusa do aglomerado da atividade magnetosférica da pulsar.