A Multiwavelength Interpretation of HESS J1857+026 Emission Using the Fermi-LAT, VERITAS, and HAWC Observatories

Resumo Técnico: Uma Interpretação Multicomprimento de Onda da Emissão de HESS J1857+026

Definição do Problema
HESS J1857+026 é uma fonte de raios gama de muito alta energia (VHE) estendida e não identificada, coincidente com o pulsar energético PSR J1856+0245. Embora a fonte seja uma candidata a uma nebulosa de vento de pulsar (PWN), sua natureza multicomprimento de onda permanece complexa. Observações anteriores revelaram um descompasso espacial entre as extensões das emissões em GeV e TeV e uma falta de contrapartidas de raios X confirmadas, levantando questões sobre o mecanismo de emissão dominante. Especificamente, não está claro se a emissão é puramente leptônica (espalhamento Compton inverso de uma PWN), puramente hadrônica (decaimento de pioneiro de uma remanescente de supernova interagindo com gás molecular), ou uma combinação de ambas. Além disso, as propriedades de difusão de partículas no ambiente imediato de PWNe evoluídas permanecem mal restringidas quando comparadas ao meio interestelar (ISM) médio.

Metodologia
Este estudo apresenta uma análise abrangente multi-instrumental de HESS J1857+026 utilizando 15 anos de dados do Fermi-LAT (300 MeV–2 TeV), dados do VERITAS (0,3–10 TeV) e novos dados do HAWC (1–316 TeV). A análise procede através de várias etapas:

1. Análise de Dados:

   • Fermi-LAT: Uma análise de verossimilhança agrupada (binned likelihood analysis) foi realizada em dados da classe Pass 8 SOURCE usando templates de fontes difusas e o catálogo 4FGL-DR3. Testes de extensão espacial compararam templates de ponto, disco radial e Gaussiana radial.
   • VERITAS: Uma análise de verossimilhança máxima 3D agrupada foi conduzida em 30 horas de exposição usando o framework Gammapy, modelando fontes como leis de potência com formas espaciais Gaussianas simétricas.
   • HAWC: Uma análise de verossimilhança máxima de forward-folding foi aplicada ao conjunto de dados Pass 5 (2860 dias). O estudo testou múltiplos templates espaciais (Gaussiano vs. baseado em difusão) e formas espectrais (lei de potência, log-parábola, lei de potência com corte exponencial) usando critérios de informação Bayesiano e de Akaike (BIC/AIC) para selecionar os melhores modelos de ajuste.

2. Modelagem Radiativa:

   • Independente do Tempo: Distribuições de energia espectral (SEDs) de banda larga foram ajustadas usando o pacote NAIMA para testar cenários puramente hadrônicos (SNR) e lepto-hadrônicos (PWN/SNR).
   • Dependente do Tempo: Um modelo evolutivo (semelhante a Gelfand et al. 2009) foi empregado para simular a expansão e o resfriamento da PWN dentro de uma casca de SNR. Um procedimento de Monte Carlo via Cadeias de Markov (MCMC) restringiu parâmetros físicos como a idade do sistema, a força do campo magnético e os espectros de injeção de partículas.

3. Análise de Difusão:

   • Perfis de brilho superficial radiais foram extraídos de todos os três instrumentos. Esses perfis foram ajustados com um modelo analítico baseado em difusão (Equação 16) para restringir o coeficiente de difusão de partículas, $D(E)$, e a supressão da difusão em relação ao ISM.

Principais Resultados

• Características Espectrais e Espaciais:

  • Os dados do Fermi-LAT revelam uma fonte estendida que é melhor ajustada por uma Gaussiana radial com um raio de $0,38^\circ \pm 0,04^\circ$ e um índice espectral $\Gamma \approx 2,07$. Um quebra espectral é observado, com um índice mais suave ($\Gamma \approx 2,63$) abaixo de 5 GeV e um índice mais duro ($\Gamma \approx 1,99$) acima de 5 GeV.
  • O VERITAS detecta a fonte com uma significância de $15,2\sigma$, medindo um índice espectral de $\Gamma = 2,48 \pm 0,06$ e uma largura Gaussiana de $0,236^\circ \pm 0,017^\circ$.
  • Os dados do HAWC mostram que a fonte se estende até $\sim 37$ TeV. O melhor modelo espacial de ajuste favorece um template de difusão em vez de um simples Gaussiano para um campo magnético de $1 \mu$G, com uma energia de corte $E_c \approx 14$ TeV.
  • O tamanho da fonte diminui com o aumento da energia (de $0,38^\circ$ em energias de GeV para $0,26^\circ$ em energias de HAWC), indicando uma morfologia dependente da energia consistente com uma PWN evoluída.

• Origem Radiativa:

  • Lepto-hadrônica vs. Hadrônica: Embora um modelo puramente hadrônico não possa ser inteiramente descartado para o componente de baixa energia ($<10$ GeV), o modelo lepto-hadrônico fornece um ajuste superior aos dados de banda larga. A falta de uma detecção da casca de SNR e a baixa densidade de gás ambiente ($n_0 \lesssim 0,13$ cm$^{-3}$) derivada do modelo evolutivo desfavorecem uma origem hadrônica dominante para a emissão de TeV.
  • Evolução da PWN: O modelo evolutivo sugere que o sistema é uma PWN evoluída com uma idade de $\tau = [16, 21]$ kyr e uma baixa força de campo magnético de $B = [0,4, 1,6] \mu$G. Esses parâmetros explicam a fraca emissão de raios X (devido ao resfriamento por sincrotron) e a morfologia dependente da energia. O baixo campo magnético também sustenta a dominância do espalhamento Compton inverso (ICS) sobre campos de fótons locais para a emissão de TeV.

• Difusão de Partículas:

  • O coeficiente de difusão ao redor de HESS J1857+026 foi encontrado como suprimido em relação ao ISM galáctico médio. Em 50 TeV, o coeficiente de difusão é $D_{PWN} \sim 2,14 \times 10^{28}$ cm$^2$ s$^{-1}$ (para $B=1 \mu$G), o que é aproximadamente dois ordens de magnitude inferior ao valor do ISM ($D_{ISM} \sim 1,2 \times 10^{30}$ cm$^2$ s$^{-1}$).
  • A supressão ocorre principalmente na normalização ($D_0$) e não na dependência de rigidez (índice espectral $\delta$), que permanece consistente com modelos de turbulência Kolmogorov ($\delta=1/3$) ou Kraichnan ($\delta=1/2$). Isso apoia o cenário de turbulência magnética localmente aumentada perto da fonte.

Significância e Alegações
O artigo afirma que os dados combinados de múltiplos comprimentos de onda fornecem uma ferramenta poderosa para restringir as propriedades de difusão de partículas nas proximidades de PWNe evoluídas. A conclusão primária é que HESS J1857+026 é muito provavelmente uma PWN alimentada por PSR J1856+0245, onde a emissão de raios gama observada é dominada por elétrons relativísticos sofrendo espalhamento Compton inverso.

O estudo destaca que o componente espectral de baixa energia ($E < 10$ GeV) ainda pode abrigar uma contribuição hadrônica de uma SNR hospedeira ou próxima, ou alternativamente, uma população distinta de elétrons dentro da mesma PWN. No entanto, a emissão de alta energia é firmemente leptônica. O fator de supressão de difusão derivado é consistente com outras PWNe de TeV (ex: Geminga, Monogem), reforçando a hipótese de que o transporte de partículas perto dessas fontes é significativamente dificultado pela turbulência magnética local em comparação com o ISM geral. Os autores observam que, embora um excesso de neutrinos tenha sido relatado com baixa significância, o grande limite superior da região de emissão de neutrinos em comparação com a região de raios gama sugere uma origem complexa ou não relacionada para qualquer potencial sinal de neutrinos.