The Effects of Cosmic Ray Protons on Galactic Nonthermal Filaments

Resumo Técnico: Os Efeitos de Prótons de Raios Cósmicos em Filamentos Não Térmicos Galácticos

Problema
O Centro Galáctico (CG) abriga uma população de filamentos não térmicos (NTFs) que se estendem por dezenas de parsecs de comprimento, iluminados por radiação síncrotron de raios cósmicos (CRs). A origem física dessas estruturas permanece sem resolução. Dois mecanismos primários de injeção são propostos: (1) os NTFs são alimentados por jatos de nebulosas de vento de pulsar, resultando em uma população dominada por léptons; ou (2) os NTFs são alimentados por partículas aceleradas de choques interestelares, resultando em uma população dominada por prótons (com uma razão próton-lépton de aproximadamente 100:1). Embora a população de CRs galácticos exiba uma distribuição de energia em lei de potência consistente com as observações dos NTFs, a conexão causal específica entre a fonte de injeção e a morfologia do filamento permanece incerta. Além disso, a topologia do campo magnético poloidal dentro da Zona Molecular Central (CMZ) é complexa, com muitos filamentos parecendo não lineares, sugerindo o papel potencial da turbulência interestelar. Este estudo visa distinguir entre os cenários dominados por léptons e dominados por prótons simulando a propagação de CRs, o aquecimento e o fluxo de plasma sob variadas condições do Centro Galáctico.

Metodologia
Os autores empregaram o código de magnetohidrodinâmica (MHD) Athena++, modificado para incluir termos de perda radiativa e colisional, para simular a propagação de CRs em um tubo de fluxo 1D representando um filamento. A simulação modela um meio de fundo uniforme perturbado pela injeção de energia de CR no ponto médio do filamento.

• Processos Físicos: O modelo incorpora difusão (espalhamento em ondas de Alfvén) e streaming (convecção na velocidade de Alfvén). Os mecanismos de perda para léptons incluem radiação síncrotron, espalhamento de Compton inverso, perdas de Coulomb e bremsstrahlung. A produção secundária é modelada onde prótons primários interagem com prótons ambientes para gerar elétrons e pósitrons secundários.
• Parâmetros: O estudo varia parâmetros fundamentais incluindo a força do campo magnético ($B$), a densidade do plasma ($n_p$), o coeficiente de difusão de CR ($\kappa$) e a energia total injetada ($E_{\pm}$).
• Conjuntos de Simulação: Cada conjunto de parâmetros foi simulado duas vezes: uma vez apenas com léptons e outra com uma mistura de prótons e léptons (escalonando a taxa de injeção total por 100 para contabilizar a razão próton-lépton). Os parâmetros fiduciais adotados foram $B = 200\,\mu\text{G}$, $n_p = 100\,\text{cm}^{-3}$ e $\kappa = 1 \times 10^{26}\,\text{cm}^2\,\text{s}^{-1}$.
• Análise: Os autores analisaram os perfis espaciais de densidade de energia, fluxos de gás induzidos, taxas de aquecimento e o brilho síncrotron e índices espectrais resultantes ao longo do tempo (até 4 de 40 kyr).

Resultos Principais

1. Indistinguibilidade dos Mecanismos de Injeção: As simulações revelaram poucas diferenças observáveis entre os casos dominados por prótons e dominados por léptons. Embora a injeção de prótons impulsione fluxos de gás mais fortes e proporcione mais aquecimento devido à maior densidade de energia total, esses efeitos permanecem abaixo dos limiares detectáveis. Mesmo nas densidades mais baixas consideradas, as velocidades induzidas ($\sim 2,5\,\text{km s}^{-1}$) e as taxas de aquecimento são negligenciáveis em comparação com o resfriamento radiativo e as larguras típicas de linhas térmicas.
2. Restrições de Propagação: O brilho de superfície observado, constante e uniforme dos NTFs sugere um coeficiente de difusão significativamente menor do que o valor galáctico típico ($\kappa \sim 10^{28}\,\text{cm}^2\,\text{s}^{-1}$). Os autores descobrem que um coeficiente de difusão mais baixo ($\kappa \sim 10^{26}\,\text{cm}^2\,\text{s}^{-1}$) é necessário para reproduzir os perfis de luminosidade observados.
3. Comprimento e Idade do Filamento: As simulações indicam que, sem um mecanismo para truncar a propagação, os filamentos continuariam a se expandir indefinidamente. O comprimento de um NTF é, portanto, dependente de sua idade. Perdas radiativas e colisionais, mesmo quando combinadas com a produção secundária, não criam um comprimento de estado estacionário porque os léptons secundários nascem mais longe do local de injeção e difundem-se de forma semelhante às partículas primárias.
4. Variabilidade do Índice Espectral: O estudo explorou mecanismos para a variação observada nos índices espectrais de síncrotron (variando de -2 a +0,5).
   • Difusão dependente de energia: Se a difusividade escala com a energia ($\kappa \propto \gamma^q$), o índice espectral varia ao longo do filamento, geralmente endurecendo ao se afastar do local de injeção.
   • Envelhecimento de síncrotron: Partículas de alta energia perdem energia mais rapidamente, levando ao steepening (suavização) espectral ao longo do tempo.
   • Produção secundária: Suaviza o espectro em tempos tardios.
   • Apesar desses mecanismos, os autores observam que nenhuma combinação única desses efeitos explica totalmente as variações de índice espectral irregulares e assimétricas observadas em NTFs reais.

Significância e Conclusões
A principal contribuição deste trabalho é a demonstração de que a presença de uma grande população de prótons de raios cósmicos não produz assinaturas distintas e observáveis na morfologia, fluxo ou aquecimento dos filamentos que permitiriam aos astrônomos distinguir entre origens impulsionadas por pulsares (léptons) ou por choques (prótons) baseando-se em simulações MHD 1D atuais.

Consequentemente, os autores argumentam que a incapacidade de distinguir entre esses dois modelos padrão de injeção, combinada com as imperfeições estruturais observadas nos NTFs (dobras, divisões e alinhamentos não lineares), motiva a consideração de um terceiro mecanismo de formação. O artigo propõe que os NTFs podem surgir de estruturas intermitentes na turbulência do Centro Galáctico. Neste cenário, a turbulência MHD não homogênea poderia amplificar campos magnéticos e esticar redemoinhos em filamentos finos, criando as estruturas observadas sem exigir uma fonte de injeção de CR específica e localizada. Os autores concluem que trabalhos futuros devem focar em simular esses efeitos turbulentos para compreender melhor a origem dos NTFs.