Time-dependent cosmic-ray escape from wind bubbles: hard spectra formation

Resumo Técnico: Escape de raios cósmicos dependente do tempo de bolhas de vento: formação de espectros duros

Problema
Bolhas impulsionadas por ventos, geradas por fontes como aglomerados estelares massivos jovens (YMSCs), estrelas Wolf–Rayet e objetos compactos com acreção super-Eddington, são candidatas promissoras para a aceleração de raios cósmicos (CRs) até energias de PeV (PeVatrons). Embora a aceleração por choque difusivo (DSA) no choque de terminação do vento interno seja amplamente aceita como o principal mecanismo de aceleração, o transporte dependente do tempo dessas partículas do choque de terminação para o choque frontal (a fronteira de escape para o meio interestelar) permanece pouco explorado. Estudos anteriores frequentemente basearam-se em descrições estacionárias ou estimativas de escalas de tempo de transporte, em vez de resolver a equação de transporte completa dependente do tempo. Este trabalho aborda a lacuna na compreensão de como a expansão dinâmica da bolha e a distância evolutiva no tempo entre o choque de terminação e o choque frontal impactam o espectro de energia dos raios cósmicos que escapam.

Metodologia
Os autores modelam a bolha de vento como um sistema esfericamente simétrico onde tanto o choque de terminação (TS) quanto o choque frontal (FS) se expandem com velocidades dependentes do tempo ($R_{ts} \propto t^{2/5}$ e $R_{fs} \propto t^{3/5}$ durante a fase de desaceleração). O transporte de CRs é governado por uma equação de transporte esférica unidimensional dependente do tempo (equação de Fokker-Planck), que inclui advecção pelo fluxo do vento chocado e difusão.

Para resolver esta equação, os autores a transformam em um conjunto de equações diferenciais estocásticas (SDEs) usando o léma de Itô. Estas SDEs são integradas utilizando uma versão modificada do framework de código aberto CRPropa3.2. Principais características metodológicas incluem:

• Injeção: Partículas são injetadas no TS com um espectro de lei de potência pré-acelerado ($dN/dE \propto E^{-2}$) e um corte de alta energia dependente do tempo determinado pelo limite de Hillas e pelas escalas de tempo de aceleração.
• Advecção: Assume-se um perfil de velocidade radial $u(r) \propto r^{-2}$ na região downstream, normalizado à velocidade pós-choque.
• Difusão: O estudo investiga três modelos de turbulência: Bohm ($\delta=1$), Kraichnan ($\delta=3/2$) e Kolmogorov ($\delta=5/3$). O coeficiente de difusão depende da força do campo magnético local e do espectro de potência da turbulência.
• Campo Magnético: O campo magnético é estimado com base em uma fração da pressão de radiação. O estudo testa várias dependências radiais do campo magnético downstream ($B \propto r^{-\alpha}$ com $\alpha = 0, -1, -2$) e níveis de magnetização ($\epsilon_B$).
• Parâmetros de Simulação: As simulações rodam por até 5 Myr, cobrindo a fase impulsionada pelo vento antes que as supernovas dominem o orçamento de energia.

Principais Resultados
O tratamento dependente do tempo revela desvios significativos da expectativa padrão dos espectros de escape:

1. Endurecimento Espectral: Em configurações físicas específicas, os espectros de CR que escapam podem ser mais duros que o espectro injetado de $E^{-2}$. Este endurecimento é mais pronunciado para modelos de turbulência com dependência de energia mais forte (ex: difusão de Bohm) e dependências de campo magnético mais fracas.
2. Supressão de Baixa Energia: Uma descoberta crítica é a supressão significativa de partículas de baixa energia no fluxo de escape. Como o choque frontal (fronteira de escape) se afasta da fonte ao longo do tempo, partículas de baixa energia com coeficientes de difusão pequenos não conseguem alcançar a fronteira de escape e permanecem presas na região downstream.
   • Para a difusão de Bohm, partículas de baixa energia ($E < 1$ TeV) são quase completamente depletadas do fluxo de escape.
   • Para a difusão de Kraichnan, observa-se um endurecimento moderado ($\Delta \gamma \approx 0,1-0,2$) e uma supressão parcial de baixa energia.
   • Para a difusão de Kolmogorov, a depleção é mínima e o espectro de escape assemelha-se ao espectro injetado, embora a energia máxima permaneça inferior ($\sim 20$ TeV).
3. Dependência de Parâmetros: O grau de depleção de baixa energia e o endurecimento espectral dependem do modelo de turbulência, da dependência radial do campo magnético (dependências mais fracas levam a espectros mais duros) e do nível de magnetização.
4. Implicações Observacionais: Quando propagados para a Terra (simulado para uma fonte a 3 kpc), a supressão de baixa energia persiste, embora o endurecimento seja parcialmente suavizado pela difusão galáctica ao longo de longas escalas de tempo.

Significância e Alegações
O artigo afirma que a natureza dependente do tempo da expansão da bolha de vento é um fator crucial na modelagem do espectro de CR observado, independentemente dos processos de perda de energia. Os autores enfatizam que:

• Formação de Espectros Duros: Bolhas de vento podem naturalmente produzir espectros de CR duros (mais duros que $E^{-2}$) sem requerer física de aceleração não convencional, desde que o transporte seja dependente do tempo.
• Corte de Baixa Energia: O "aprisionamento" de partículas de baixa energia dentro da bolha é uma previsão robusta deste modelo dependente do tempo. Isso pode levar a características observáveis em radiação multi-mensageira (raios gama e neutrinos) de dentro da bolha e afetar a gramática acumulada dos CRs.
• Distinção de SNRs: Os espectros duros resultantes e os cortes de baixa energia diferem daqueles tipicamente esperados de remanescentes de supernovas (SNRs), potencialmente ajudando a explicar o espectro de CR ao redor do "knee" ou observações específicas de raios gama de aglomerados estelares (ex: Westerlund 1).

Os autores concluem que, embora seu modelo foque na fase impulsionada pelo vento, essas partículas de baixa energia presas podem ser posteriormente re-aceleradas por choques de supernovas ou escapar quando o choque frontal cessa, podendo potencialmente contribuir para a população de CR de baixa energia com uma forma espectral modificada. Trabalhos futuros incorporarão interações hadrônicas para prever emissões de raios gama e neutrinos dependentes do tempo.