Kick bimodality of neutron stars and mode dependence of their parameters

Enunciado do Problema
Estrelas de nêutrons (ENs) são observadas possuindo velocidades espaciais que excedem significativamente as de suas estrelas progenitoras, um fenômeno atribuído a chutes natais recebidos durante explosões de supernovas. Embora uma distribuição bimodal de velocidades de chute (compreendendo modos de baixa e alta velocidade) tenha sido sustentada por modelagem teórica e síntese de população (ex., Igoshev 2020), a confirmação observacional permanece desafiadora. Derivar a distribuição de chute natal a partir de velocidades transversais 2D observadas é não trivial devido à velocidade radial desconhecida, incertezas sobre os locais de nascimento e evolução cinemática dentro do potencial Galáctico. Além disso, estudos existentes não investigaram minuciosamente se parâmetros físicos específicos de pulsares de rádio correlacionam-se com esses distintos modos de velocidade de chute, particularmente distinguindo entre efeitos causados por viés de seleção e diferenças físicas intrínsecas.

Metodologia
Os autores analisaram uma amostra de 202 pulsares de rádio normais e isolados selecionados do Catálogo de Pulsares ATNF (v2.7.0). Os critérios de seleção exigiram medições robustas de paralaxe ou distâncias de medida de dispersão < 10 kpc, movimentos próprios disponíveis, parâmetros de spin ($P, \dot{P}$) conhecidos e a exclusão de pulsares reciclados, binários e fontes exóticas.

Para classificar esses pulsares em modos de baixa ou alta velocidade, os autores empregaram um modelo de integração de trajetória:

1. Integração de Trajetória: Para cada pulsar, 100 valores de distância foram amostrados (baseados em paralaxe ou medida de dispersão). Usando coordenadas galácticas e movimentos próprios, 100 posições atuais foram geradas.
2. Reconstrução de Velocidade: Como as velocidades radiais ($v_r$) não são observadas, os autores assumiram uma orientação de velocidade isotrópica relativa ao Padrão Local de Repouso (LSR). Eles amostraram o ângulo $\theta$ entre a linha de visada e o vetor de velocidade 3D a partir de uma distribuição uniforme truncada (Uniform(0.05$\pi$, 0.95$\pi$)) para evitar valores extremos, gerando $10^3$ vetores de velocidade por distância. Isso resultou em $10^5$ combinações de posição-velocidade por pulsar.
3. Integração Reversa: Usando o potencial Galáctico (McMillan 2017), as trajetórias foram integradas reversamente no tempo até $1,5 \times \tau_c$ (idade característica). Os locais de nascimento foram identificados ou por interseções com o plano Galáctico ou em $t = -\tau_c$ caso nenhuma interseção ocorresse.
4. Cálculo da Velocidade de Chute: Em cada potencial local de nascimento, a velocidade de chute ($v_{kick} = |\vec{v} - \vec{v}_{LSR}|$) foi calculada.
5. Categorização: Usando os parâmetros da distribuição de chute bimodal de Igoshev (2020) (uma soma ponderada de duas Maxwellianas: $\sigma_1 = 56$ km s$^{-1}$, $\sigma_2 = 336$ km s$^{-1}$, com pesos $w_1=0,2, w_2=0,8$), os autores calcularam a probabilidade de cada pulsar pertencer ao modo de baixa ou alta velocidade. Os pulsares foram categorizados como pertencentes "inequivocamente" a um modo se a probabilidade excedesse 0,5 em todos os possíveis locais de nascimento; caso contrário, foram categorizados com base na maioria dos locais de nascimento.

Principais Resultados

• Distribuição de Modos: Aproximadamente 30% (60 de 202) da amostra pertence ao modo de baixa velocidade, consistente com estimativas anteriores para pulsares jovens. Os ~70% restantes pertencem ao modo de alta velocidade.
• Efeitos de Viés de Seleção: Como esperado, pulsares de alta velocidade foram encontrados para serem, em média, mais jovens (idades de spin-down mais curtas) e mais distantes do que os de baixa velocidade. Essas diferenças foram confirmadas por testes de Kolmogorov-Smirnov (K-S) ($p < 0,01$) e são atribuídas a efeitos de seleção observacional (dificuldade em medir pequenas velocidades a grandes distâncias e objetos de alta velocidade deixando o volume de pesquisa).
• Largura de Pulso: Nenhuma diferença significativa foi encontrada na largura de pulso ($W_{10}$ e $W_{50}$) entre os dois modos ($p \approx 0,5$), sugerindo que este parâmetro não está sujeito a forte viés de seleção ou dependência intrínseca do modo.
• Correlação de Campo Magnético: Uma diferença estatisticamente significativa ($p \approx 0,02$) foi detectada nas distribuições de campo magnético superficial. Pulsares no modo de baixa velocidade apresentam uma abundância excessiva de campos magnéticos mais baixos ($B \lesssim 10^{12}$ G). Notavelmente, nenhum pulsar com $B \lesssim 10^{11}$ G foi identificado no modo de alta velocidade. Esta discrepância persiste mesmo ao considerar distância e idade, embora os autores notem que pode ser parcialmente influenciada pelo decaimento do campo magnético ao longo do tempo.

Significância e Alegações
O artigo afirma fornecer uma classificação detalhada de pulsares de rádio normais em grupos de velocidade de chute bimodal e analisar as distribuições de parâmetros resultantes. A principal significância reside na identificação de uma correlação entre os modos de velocidade de chute e a intensidade do campo magnético. Os autores afirmam que, embora o viés de seleção e o decaimento do campo magnético possam contribuir para a diferença observada, a origem da discrepância não é totalmente clara. Eles sugerem que, se a estrutura de chute bimodal decorre de diferenças na física da explosão de supernovas, a diferença de campo magnético pode compartilhar uma origem semelhante, potencialmente ligando os mecanismos de chute de ENs ao campo magnético dipolar externo. No entanto, os autores mantêm-se modestos, reconhecendo que a modelagem detalhada de síntese de população é necessária para confirmar essas origens e que os efeitos de seleção não podem ser inteiramente descartados. O estudo não propõe novas aplicações experimentais, mas oferece uma análise estatística de dados existentes para refinar a compreensão da cinemática e propriedades magnéticas de ENs.