Analysis and implications of the spatio-spectral morphology of the Fermi Bubbles

Resumo Técnico: Análise e Implicações da Morfologia Espácio-Espectral das Bolhas de Fermi

Enunciado do Problema
A origem e a fonte de potência das Bolhas de Fermi (FBs)—estruturas gigantes de raios $\gamma$ que se estendem $\sim 50^\circ$ acima e abaixo do plano Galáctico—permanecem incertas. Embora sua morfologia sugira uma origem no Centro Galáctico (GC), o mecanismo microfísico que impulsiona a emissão de raios $\gamma$ é debatido entre dois canais primários: o canal hadrônico (colisões inelásticas de prótons de raios cósmicos com o gás ambiente) e o canal leptônico (espalhamento Compton inverso de fótons do campo de radiação interestelar por elétrons de raios cósmicos). Análises anteriores basearam-se em representações espaciais grosseiras (ex: fatias com média de latitude), o que pode obscurecer variações espectrais dependentes da posição. Este artigo visa resolver essas incertezas realizando uma análise espectral de alta resolução, pixel a pixel, para restringir as populações de raios cósmicos (CR) parentais necessárias para reproduzir os espectros de raios $\gamma$ observados.

Metodologia
Os autores utilizam uma reconstrução livre de modelos (template-free) de dez anos de dados de todo o céu do Fermi/LAT fornecidos por Platz et al. (2023). Especificamente, eles empregam o componente "não-poeira" ($I_{nd}$) do modelo M2, que isola a emissão não correlacionada com mapas de poeira térmica, revelando as FBs claramente contra o primeiro plano Galáctico. O céu é pixelizado usando o esquema HEALPix ($N_{side}=128$), resultando em aproximadamente 20.000 pixels cobrindo as FBs.

Para modelar a emissão, os autores:

1. Constroem uma Geometria 3D: Eles aproximam as FBs como figuras de rotação geradas pela rotação de um contorno traçado à mão em torno do eixo $z$ Galáctico. Isso permite o cálculo das profundidades na linha de visada e a integração de quantidades físicas ao longo de cada raio.
2. Modelam o Campo de Radiação Interestelar (ISRF): Utilizam o modelo de transferência radiativa axissimétrico de Popescu et al. (2017) para determinar o ISRF dependente da posição (componentes estelar, de poeira e CMB) ao longo de cada linha de visada. Para regiões além da grade do modelo ($z > 10$ kpc), eles extrapolam usando escalonamento de lei de potência ou $z^{-1}$.
3. Modelagem Direta (Forward Modeling): Usando o pacote Python Naima, eles ajustam os espectros de raios $\gamma$ observados em cada pixel contra três modelos de população parental: Lei de Potência (PL), Lei de Potência com Corte Exponencial (EPL) e Lei de Potência Quebrada (BPL).
4. Hadrônico vs. Leptônico: Para modelos hadrônicos, assumem uma densidade numérica de hidrogênio uniforme de $n_H = 10^{-3} \text{ cm}^{-3}$ e normalizam os resultados de acordo. Para modelos leptônicos, calculam o espalhamento Compton Inverso (IC) sobre o ISRF médio derivado para cada pixel.
5. Bondade do Ajuste: Eles avaliam os modelos usando $\chi^2$ reduzido, observando que, devido às incertezas correlacionadas na decomposição dos dados, esta métrica serve para comparar o desempenho relativo dos modelos, em vez de significância estatística absoluta.

Resultos Principais

• Forma Espectral: Modelos simples de lei de potência (PL) tanto para prótons quanto para elétrons proporcionam ajustes pobres aos dados e são descartados. Ambos os modelos EPL e BPL fornecem ajustes igualmente bons para os cenários hadrônicos e leptônicos.
• Índices Espectrais e Endurecimento: Os melhores índices espectrais em baixas energias mostram um endurecimento distinto em direção à extremidade sul das bolhas tanto para prótons de CR (CRp) quanto para elétrons de CR (CRe). Essa tendência é consistente com achados anteriores, mas estende-se a latitudes mais altas do que as analisadas anteriormente.
• Energias de Quebra/Corte:
  • Para modelos EPL, as energias de corte são da ordem de alguns TeV e mostram pouca variação com a latitude.
  • Para modelos BPL, as energias de quebra são de $\sim 100$ GeV com variação latitudinal mínima.
  • A mudança no índice espectral através da quebra em modelos BPL é $\lesssim 0,5$. Isso é significativamente menor do que o $\Delta \alpha \approx 1$ esperado para uma quebra de resfriamento devido a perdas de sincrotron ou IC, sugerindo que as quebras espectrais não são características de resfriamento.
• Contribuições do Campo de Radiação: Em modelos leptônicos, o componente de poeira do ISRF domina a produção de raios $\gamma$ em todas as energias do Fermi/LAT, seguido pelo componente estelar (que é suprimido por efeitos de Klein-Nishina em altas energias) e pela CMB. Consequentemente, as perdas de IC podem ser aproximadas como ocorrendo no regime de Thomson em todas as FBs.
• Densidade de Energia e Pressão:
  • Hadrônico: A energia total de CR é estimada em $\sim 1,7 \times 10^{57} / n_{H,-3}$ erg. A pressão de CR é provavelmente dominante sobre a pressão do gás térmico, mas apenas por um fator de poucas unidades, dependendo fortemente da densidade de gás alvo assumida.
  • Leptônico: A densidade de energia de CRe aumenta com a distância do plano Galáctico e atinge o pico em direção às bordas das bolhas, qualitativamente consistente com o brilho de borda (limb brightening).
• Escalas de Tempo de Resfriamento:
  • Hadrônico: A escala de tempo de resfriamento é extremamente longa ($\sim 100$ Gyr para $n_H = 10^{-3} \text{ cm}^{-3}$), implicando que as bolhas podem ser muito antigas ou estar longe do estado estacionário.
  • Leptônico: Tempos de resfriamento para elétrons de alta energia são extremamente curtos ($\lesssim 1$ Myr nos capuzes das bolhas, caindo para $\sim 0,1$ Myr perto do GC).

Significância e Implicações
O artigo conclui que nem os modelos hadrônicos nem os leptônicos podem ser definitivamente descartados baseando-se apenas na forma espectral, pois ambos os modelos EPL e BPL se ajustam bem aos dados para ambos os tipos de partículas. No entanto, os resultados impõem restrições estritas aos cenários físicos:

1. Transporte Leptônico vs. Aceleração In-Situ: A combinação de tempos de resfriamento curtos e o aumento observado na densidade de energia de CRe em direção às bordas das bolhas apresenta um desafio severo aos modelos onde os elétrons são injetados no Centro Galáctico e transportados para fora. Transportar elétrons até os capuzes das bolhas ($\sim 8$ kpc) dentro de seu tempo de resfriamento de $\sim 1$ Myr exigiria velocidades de expansão ($\sim 8000$ km/s) inconsistentes com a fenomenologia de raios X. Os autores sugerem que a reaceleração in-situ de léptons na periferia da bolha é uma explicação mais natural para a morfologia e distribuição de densidade de energia observadas.
2. Restrições de Idade Hadrônica: A longa escala de tempo de resfriamento hadrônico implica que, se as bolhas forem jovens (ex: $< 100$ Myr), o orçamento de energia necessário exigiria uma origem de AGN em vez de formação estelar. Por outro outro lado, se as bolhas forem velhas (multi-Gyr), uma origem de formação estelar permanece plausível, desde que a população de CR tenha tido tempo suficiente para se acumular.
3. Endurecimento Espectral: O endurecimento espectral observado em direção à extremidade sul e a falta de uma quebra de resfriamento ($\Delta \alpha < 0,5$) sugerem que as características espectrais são impulsionadas pelo mecanismo de injeção ou aceleração, em vez de perdas de resfriamento radiativo.

Em resumo, este trabalho fornece as primeiras restrições espectrais pixel a pixel sobre as Bolhas de Fermi, demonstrando que, embora leis de potência simples sejam insuficientes, modelos com cortes exponenciais ou quebras se ajustam bem aos dados. Os achados favorecem fortemente cenários envolvendo ou o transporte relativístico do núcleo ou, mais provavelmente, a aceleração distribuída in-situ de raios cósmicos nas bordas das bolhas, particularmente para o canal leptônico.