A Boundary-Consistent Two-Zone Electron Kernel for Distant Pulsar Contributions to Positron Flux and Anisotropy

Resumo Técnico: Um Núcleo de Elétrons de Duas Zonas com Consistência de Fronteira para Contribuições de Pulsares Distantes no Fluxo de Pósitrons e Anisotropia

Enunciado do Problema
Os espectros locais de elétrons e pósitrons de raios cósmicos servem como sondas de precisão para aceleradores de alta energia próximos, particularmente pulsares. Enquanto pulsares maduros como Geminga e Monogem são candidatos naturais para o fluxo de léptons na escala de TeV, uma questão crítica permanece: os fluxos de léptons locais devem ser dominados pelos pulsares mais próximos, ou pulsares mais distantes podem contribuir coletivamente? Esta investigação é complicada por dois fatores: (1) o rápido resfriamento radiativo (síncrotron e Compton inverso) limita o horizonte de propagação de partículas de alta energia, e (2) os limites de anisotropia da Fermi-LAT e AMS-02 restringem modelos onde um único alvo próximo domina de forma muito ampla. Estudos anteriores sugeriram que fontes distantes podem reproduzir elétrons de baixa energia, mas carecia-se de um arcabouço semi-analítico estável para quantificar o fator geométrico de contagem de fontes na presença de difusão lenta próxima à fonte. Além disso, os primeiros centenas de parsecs da Galáxia (a Bolha Local) apresentam uma descontinuidade física onde uma densidade contínua e suave de pulsares é suspeita devido à baixa densidade de gás molecular e aos chutes natais (natal kicks) dos pulsares.

Metodologia
Os autores desenvolvem uma solução de série semi-analítica para a propagação de elétrons e pósitrons em um modelo de difusão de duas zonas esféricas.

• Núcleo de Duas Zonas: O modelo define uma região de difusão lenta de raio $R_1$ ao redor da fonte, onde o coeficiente de difusão é $D_1(E)$, e uma região interestelar padrão ($R_1 < r < R_2$) com coeficiente $D_2(E)$. Os coeficientes compartilham a mesma dependência de energia ($D \propto f_D(E)$), permitindo soluções de autofunções analíticas.
• Condições de Contorno: A solução impõe continuidade de densidade e fluxo difusivo em $R_1$ e uma fronteira absorvente em $R_2$ (definida como 20 kpc). A formulação evita os cancelamentos oscilatórios inerentes às avaliações integrais diretas de modelos de duas zonas, preservando fronteiras nítidas de resfriamento radiativo.
• Resfriamento: As perdas radiativas (síncrotron e Compton inverso de Klein–Nishina sobre campos de CMB, infravermelho e óptico) são tratadas via características de energia, vinculando a energia de injeção à idade da fonte.
• Modelagem de População: O núcleo é aplicado a uma população contínua de pulsares em disco. Os autores impõem uma cavidade vazia para $R < 0,1$ kpc para contabilizar a Bolha Local e os movimentos próprios dos pulsares.
• Análise Estatística: Uma análise de Monte Carlo via Cadeia de Markov (MCMC) com "solucionador no loop" é realizada para ajustar o modelo aos dados de fluxo de pósitrons da AMS-02, medições de dipolo de anisotropia de todos os elétrons da CALET e limites superiores de anisotropia da Fermi-LAT. O ajuste inclui histórias de injeção de spin-down contínuas, um espectro de injeção de pares de lei de potência interrompida e perfis sobre amplitudes lineares não negativas para uma base de pósitrons do tipo secundário e um fundo de lei de potência suave.

Contribuições Principais

1. Solução de Série Estável: O artigo fornece uma solução de série semi-analítica numericamente estável para o problema de difusão de duas zonas com resfriamento radiativo. Este método é computacionalmente eficiente e evita as instabilidades numéricas de métodos de volume finito de força bruta ou integrais oscilatórios, particularmente no regime de distância média.
2. Quantificação Geométrica: O estudo quantifica como o "fator de casca" (contagem de fontes vs. distância) interage com os núcleos de propagação. Demonstra que, enquanto uma distribuição 3D uniforme é dominada por cascas de escala de quiloparsecs a 10–100 GeV, uma geometria de disco ($h=0,2$ kpc) desloca a contribuição para dentro, porém fontes além de 1 kpc ainda fornecem 37–47% do fluxo de 10–100 GeV.
3. Papel do Halo Lento: Os autores esclarecem que a inclusão de um halo de difusão lenta não aumenta o fluxo de fontes distantes em relação a um modelo de uma zona; em vez disso, introduz uma penalidade de resfriamento dependente da energia que reduz o fluxo absoluto de todas as distâncias, com uma supressão mais forte em energias mais altas.
4. Restrições de Anisotropia: O trabalho integra restrições de anisotropia diretamente no ajuste da população, usando amplitudes de dipolo da CALET e limites superiores da Fermi-LAT para testar a viabilidade de parâmetros de halo comum.

Resultos

• Decomposição de Fluxo: Na população de pulsares contínua ajustada (excluindo $R < 0,1$ kpc), a contribuição de fontes além de 0,3 kpc domina o anel local de 0,1–0,3 kpc em todas as energias de 10 GeV até a escala de TeV. A 10 GeV, a contribuição de $R > 0,3$ kpc é aproximadamente uma ordem de magnitude maior que a do anel local.
• Parâmetros do Halo: A análise MCMC produz uma distribuição posterior para o raio do halo lento $R_1$ centrada em $0,0709^{+0,0510}_{-0,0407}$ kpc ($\approx 71$ pc) e um coeficiente de difusão $D_1(100 \text{ TeV}) \approx 2,08 \times 10^{28}$ cm$^2$ s$^{-1}$. Este resultado é consistente com estudos de morfologia de halos de TeV de Geminga e Monogem.
• Compatibilidade de Anisotropia: O modelo ajustado, com os parâmetros de halo derivados, permanece compatível com os limites atuais de anisotropia (CALET e Fermi-LAT). As previsões de amplitude de dipolo situam-se dentro dos limites superiores de 95%.
• Dependência de Energia: Embora fontes distantes dominem o fluxo sub-TeV, o resfriamento radiativo faz com que a contribuição se localize próximo à escala de TeV. A 1 TeV, aproximadamente 90% do fluxo no modelo de disco provém de dentro de 1 kpc.

Significância e Alegações
O artigo conclui que uma contribuição coletiva de pulsares distantes (além da vizinhança imediata da Bolha Local) é um contribuinte plausível e dominante para o fluxo de pósitrons e elétrons abaixo de algumas centenas de GeV. O componente de pulsar ajustado não é dominado pelas fontes mais próximas, mas por uma população que se estende além de 0,3 kpc.

No entanto, os autores mantêm uma postura modesta quanto à unicidade de suas descobertas:

• Os dados de fluxo e de anisotropia sozinhos não determinam unicamente o tamanho do halo ($R_1$) ou o coeficiente de difusão ($D_1$); existem degenerações com formas de injeção, modelagem de fundo, histórias de liberação e incompletude de catálogos.
• O tamanho de halo derivado "escala Geminga" é uma restrição condicional baseada no modelo de população e conjunto de dados específicos utilizados, não uma medição definitiva de halos de sub-difusão físicos.
• Para passar de um teste de consistência para uma determinação da região de difusão lenta, são necessários observáveis adicionais, especificamente perfis de brilho superficial de Compton inverso em TeV de halos individuais, dados detalhados de população de pulsares (movimentos próprios, locais de nascimento) e likelihoods de anisotropia vetorial completos.

O trabalho serve como uma demonstração controlada pela propagação para testar o espaço de parâmetros de halo comum, fornecendo um núcleo estável para futuras inferências de população, ao mesmo tempo em que destaca a necessidade de entradas multimensageiras para uma solução completa.