Mirage Sources and Large TeV Halo-Pulsar Offsets: Exploring the Parameter Space

Resumo Técnico: Fontes Miragem e Deslocamentos de Halo de Pulsar em TeV

Enunciado do Problema
Observações recentes do Large High Altitude Air Shower Observatory (LHAASO) e do HAWC identificaram emissões estendidas de raios gama em TeV ("halos de pulsar") ao redor de pulsares de meia-idade. Um desafio significativo na interpretação dessas fontes é a observação frequente de grandes deslocamentos espaciais entre a posição do pulsar e o centroide da emissão de TeV. Modelos tradicionais que assumem a difusão isotrópica de pares de elétrons-pósitrons ($e^\pm$) no meio interestelar (ISM) não podem explicar naturalmente esses deslocamentos, pois a difusão isotrópica é apenas uma aproximação válida quando a escala física excede múltiplas escalas de coerência magnética. Consequentemente, para elétrons de ultra-alta energia ($>100$ TeV) com taxas de resfriamento rápidas, a escala de propagação pode ser menor que o comprimento de coerência do campo magnético turbulento do ISM. Como resultado, o transporte de partículas pode ser altamente assimétrico, seguindo as linhas de campo magnético locais em vez de difundir uniformemente. Este artigo investiga se tal propagação assimétrica pode gerar fontes "miragem" — picos de emissão aparentes deslocados do acelerador — e deslocamentos espaciais significativos.

Metodologia
Os autores utilizam uma simulação de partícula de teste acelerada por GPU, baseada em primeiros princípios, para rastrear as trajetórias de elétrons de 100 TeV injetados em um campo magnético turbulento.

• Geração do Campo Magnético: O campo magnético turbulento é modelado como uma superposição de quatro campos independentes de espectro de Kolmogorov ($n_{in} = -11/3$) gerados via um método de espaço de Fourier. Esta abordagem multiescala (usando espaçamentos de grade de $4 \times 10^{-9}$ a $5 \times 10^{-4}$ kpc) garante um intervalo inercial realista e campos divergência-livres ($\nabla \cdot B = 0$). O campo total inclui componentes regulares ($B_r$) e turbulentos ($B_t$).
• Propagação de Partículas: Os elétrons são propagados usando um algoritmo Boris pusher para resolver a equação de movimento relativística sob a força de Lorentz. A simulação contabiliza perdas de energia radiativa via radiação de sincrotron e espalhamento Compton Inverso (IC) sobre a Radiação Cósmica de Fundo (CMB) e o Campo de Radiação Interestelar (ISRF).
• Amostragem Estatística: Para amostrar eficientemente o espectro de energia inicial ($dN/dE \propto E^{-\alpha} \exp(-E/E_{cut})$), os autores utilizam uma abordagem de macropartícula ponderada. 8.192 macropartículas são injetadas, cada uma representando um intervalo logarítmico de energia com um peso proporcional à contagem física de elétrons naquele intervalo.
• Síntese Observacional: As distribuições de elétrons simuladas são convertidas em mapas sintéticos de raios gama (banda de 25–40 TeV) via uma matriz de espalhamento IC. Esses mapas são convoluídos com uma Função de Espalhamento de Ponto (PSF) Gaussiana de $\sigma_{PSF} = 0,43^\circ$ (mimando o LHAASO-KM2A) e adicionados a um modelo de ruído de fundo Poissoniano.
• Identificação de Fontes: Um framework de comparação de modelos Bayesiano usando Cadeias de Markov Monte Carlo (MCMC) é aplicado. Modelos com números crescentes de fontes Gaussianas são ajustados aos dados sintéticos. Uma fonte é considerada "identificada" se o Estatístico de Teste (TS) aumentar mais de 25 ($\Delta \log L > 12,5$) ao adicionar um componente.

Contribuições Principais e Resultados
O estudo explora sistematicamente o espaço de parâmetros, variando a distância do pulsar ($d$), o índice espectral de injeção de elétrons ($\alpha$), a força do campo magnético total ($B$), o comprimento de coerência da turbulência ($L_c$) e a razão entre a força dos campos regular e turbulento ($B_r/B_t$).

1. Formação de Fontes Miragem: As simulações confirmam que a propagação assimétrica pode criar múltiplas fontes de raios gama distintas ("fontes miragem") a partir de um único acelerador. Isso ocorre quando os elétrons propagam-se ao longo das linhas de campo magnético que, devido a efeitos de projeção, parecem originar-se de um local distante do ponto de injeção.
2. Multiplicidade de Fontes ($N$):
   • O número de fontes identificadas é fortemente influenciado pela força do campo magnético ($B$). Um $B$ mais alto leva a um resfriamento mais rápido, confinando os elétrons a escalas menores e reduzindo a probabilidade de formar múltiplas fontes distintas (menor $N$).
   • O número de fontes é negativamente correlacionado com a distância ($d$) e o comprimento de coerência ($L_c$), e positivamente correlacionado com $B^2$.
3. Separação de Fontes e Razão de Fluxo:
   • Separação: A separação máxima entre fontes identificadas é positivamente correlacionada com $L_c$ e negativamente correlacionada com $B^2$.
   • Razão de Fluxo ($\chi$): A razão entre o fluxo da fonte principal e o fluxo da fonte miragem diminui com o aumento de $B^2$. Em campos magnéticos fortes, a fonte miragem pode, às vezes, ser mais brilhante que a fonte principal.
   • Campos Regulares vs. Turbulentos: A distância de separação exibe uma dependência não monotônica de $B_r/B_t$. A separação máxima ocorre em $B_r/B_t \approx 1$, onde o campo regular fornece orientação direcional enquanto a turbulência permite o espalhamento suficiente para criar deslocamentos. Campos puramente turbulentos ou campos regulares dominantes resultam em menores separações.
4. Deslocamentos Espaciais:
   • Os deslocamentos entre o pulsar e a fonte identificada mais próxima são maximizados quando $B_r/B_t \approx 1$ e $L_c$ é grande.
   • Em regiões com grandes comprimentos de coerência (ex: regiões entre braços ou o halo Galáctico), os deslocamentos podem atingir até 50 pc.

Significância e Alegações
O artigo afirma fornecer uma explicação quantitativa para fontes de LHAASO intrigantes onde os halos de TeV parecem espacialmente desconectados de seus pulsares associados. Especificamente, os autores sugerem que fontes como 1LHAASO J0216+4237u, J0212+4254u e J0206+4302u poderiam ser halos miragem alimentados por um único pulsar de milissegundo (PSR J0218+4232), em vez de aceleradores distintos.

O estudo enfatiza que a aproximação de difusão isotrópica é insuficiente para modelar o transporte de partículas perto de aceleradores no ISM, particularmente para léptons de alta energia onde os tempos de resfriamento são mais curtos que o tempo de coerência magnética. Ao demonstrar que efeitos de projeção em campos magnéticos turbulentos podem naturalmente produzir grandes deslocamentos e múltiplas morfologias de fontes, o artigo oferece um mecanismo para reconciliar o céu de TeV observado com os modelos padrão de nebulosas de vento de pulsar. Os autores observam que, embora sua turbulência sintética assuma isotropia (ao contrário da turbulência Alfvénica anisotrópica em MHD), ela captura qualitativamente a natureza filamentar de grande escala da turbulência do ISM responsável por essas emissões assimétricas.