Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment

Resumo Técnico: Implementação do Algoritmo CRESP no Código MHD PIERNIK – Propagação de Núcleos Primários e Secundários

1. Problema e Contexto

Os Raios Cósmicos (RCs) são uma componente essencial do Meio Interestelar (MIS) e desempenham um papel crucial na dinâmica galáctica, incluindo a produção de raios gama e radiação síncrotron. Compreender a sua propagação é um desafio complexo que envolve teoria, experimentos (como o AMS-02) e simulações numéricas.
Existem duas abordagens principais para modelar a propagação de RCs:

1. Modelos Fenomenológicos (ex: GALPROP): Baseiam-se na equação de Fokker-Planck em um ambiente interestelar estacionário.
2. Modelos Autoconsistentes: Acoplam a propagação de RCs às equações de Magnetohidrodinâmica (MHD) dependentes do tempo, considerando a interação dinâmica entre RCs, gás térmico e campos magnéticos.

O problema central abordado neste trabalho é a extensão de um código MHD autoconsistente (PIERNIK) para incluir a propagação de núcleos de raios cósmicos espectro-resolvidos (não apenas prótons ou elétrons), permitindo a simulação simultânea de núcleos primários e secundários, incluindo processos de espalação e decaimento radioativo, em um ambiente dinâmico.

2. Metodologia

Os autores implementaram uma extensão do algoritmo CRESP (Cosmic Ray Energy SPectrum) no código PIERNIK, um código MHD multifluido baseado em grades.

• Equações Fundamentais: O sistema utiliza equações de Fokker-Planck para descrever a distribuição de função dos RCs, acopladas dinamicamente às equações de MHD do gás interestelar.
• Abordagem de Dois Momentos: Para resolver as equações de forma eficiente, o espaço de momento é dividido em "bins" (intervalos finitos). O algoritmo calcula dois momentos da função de distribuição: densidade numérica ($n_a$) e densidade de energia cinética ($e_a$).
• Aproximação de Lei de Potência por Partes: Utiliza-se uma distribuição isotrópica de lei de potência por partes dentro de cada bin para integrar as equações de transporte, permitindo lidar com regimes não-relativísticos, trans-relativísticos e ultra-relativísticos.
• Espécies Simuladas:
  • Primários: Carbono (C), Nitrogênio (N) e Oxigênio (O), acelerados em remanescentes de supernovas (SN).
  • Secundários: Lítio (Li), Berílio (Be) e Boro (B), produzidos via espalação de núcleos primários contra prótons térmicos do MIS.
  • Prótons: Tratados em aproximação "cinza" (espectro não resolvido) para atuar como motor dinâmico do sistema.
• Configuração de Simulação:
  • Caixa 3D estratificada por gravidade (0,5 x 1 x 8 kpc), representando o ambiente local do Sistema Solar.
  • Injeção aleatória de RCs a partir de supernovas.
  • Inclusão de processos físicos: advecção, difusão anisotrópica (dependente da rigidez), perdas adiabáticas e perdas de Coulomb.
• Estudo de Parâmetros: Foram realizadas simulações variando três parâmetros chave:
  1. Taxa de Supernovas (SN).
  2. Coeficiente de difusão paralelo ($D_{\parallel}$).
  3. Índice de potência da dependência da rigidez ($\delta$).

3. Contribuições Principais

• Novo Algoritmo Multi-Espécie: Validação da capacidade do código PIERNIK de modelar simultaneamente a produção e propagação de múltiplas espécies de núcleos hadrônicos (primários e secundários) com resolução espectral.
• Acoplamento Dinâmico Completo: Demonstração de como a propagação de RCs afeta e é afetada pela evolução do gás e do campo magnético em tempo real, superando a limitação de modelos estacionários.
• Descoberta de um Novo Mecanismo de Transporte: Identificação de que a razão de fluxo secundário/primário (B/C) aumenta com o aumento do coeficiente de difusão, um comportamento contrário às previsões de modelos de "caixa vazada" (leaky-box) tradicionais.

4. Resultados Chave

• Dinâmica do Gás e Campo Magnético: A injeção de RCs gera instabilidades de Parker e ventos galácticos verticais. A taxa de SN controla a força desses ventos; taxas mais altas geram ventos mais fortes, reduzindo o tempo de residência das partículas primárias no disco.
• Razão B/C (Boro/Carbono):
  • Efeito da Taxa de SN: Uma taxa de SN mais alta leva a uma menor razão B/C. Isso ocorre porque ventos mais fortes expulsam as partículas primárias do disco mais rapidamente, reduzindo a probabilidade de espalação e produção de secundários.
  • Efeito do Coeficiente de Difusão ($D_{\parallel}$): Contrariando a intuição clássica, uma maior difusão resulta em uma maior razão B/C.
    • Explicação: Maior difusão distribui os RCs de forma mais homogênea ao longo das linhas de campo horizontais, reduzindo a flutuação local de pressão que gera campos magnéticos verticais (efeito de flutuação/empuxo dos RCs). Com um campo vertical mais fraco, a difusão vertical efetiva diminui, aprisionando as partículas primárias no disco por mais tempo e aumentando a produção de secundários.
• Razão Isotópica $^{10}\text{Be}/^{9}\text{Be}$: A razão de isótopos instáveis/estáveis é sensível ao tempo de residência. Os resultados mostram que a razão varia com a altitude e os parâmetros de transporte, sendo um bom indicador do tempo de escape.
• Espectros de Energia: Os espectros evoluídos mostram inclinações (slopes) que diferem do espectro de injeção original devido à difusão dependente da rigidez e perdas de energia. As inclinações em energias não-relativísticas são mais íngremes no disco devido ao confinamento e resfriamento.

5. Significado e Conclusões

Este trabalho representa um avanço significativo na astrofísica computacional de raios cósmicos ao integrar a física de espalação nuclear em simulações MHD autoconsistentes.

• Validação do Método: O modelo consegue reproduzir qualitativamente as tendências observadas nos dados do AMS-02 e do GALPROP, validando a abordagem numérica.
• Implicações Físicas: O estudo revela que a relação entre parâmetros de transporte (como difusão) e observáveis (como B/C) não é direta em ambientes dinâmicos. A geometria do campo magnético, modificada pelos próprios RCs (via empuxo), desempenha um papel crítico na determinação do tempo de residência e, consequentemente, nas abundâncias de elementos secundários.
• Limitações e Futuro: O modelo atual utiliza uma caixa de baixa resolução e não inclui modulação solar, resfriamento térmico detalhado ou o modo de "streaming" (corrente) de RCs. No entanto, estabelece a base para simulações galácticas mais realistas que podem refinar a determinação dos parâmetros de transporte de raios cósmicos na Via Láctea.

Em suma, o artigo demonstra que a interação dinâmica entre RCs, gás e campos magnéticos é fundamental para interpretar corretamente os dados de raios cósmicos, e que modelos estáticos podem falhar em capturar comportamentos contra-intuitivos, como o aumento da razão B/C com o aumento da difusão.