Neutron Star Mass across Binary Pulsar Subpopulations: Mass-Spin Correlation, Mass Distributions, and Moment of Inertia Effects

Resumo Técnico: Massa de Estrelas de Neutrons em Subpopulações de Pulsares Binários

Definição do Problema
A formação e evolução de estrelas de nêutrons (EN) em sistemas binários estão impressas em suas massas, períodos de rotação e parâmetros orbitais. Enquanto o paradigma de "reciclagem" postula que a acreção de uma companheira faz o spin-up de ENs para períodos de milissegundos, a origem física da distribuição de massa observada permanece em debate. Especificamente, não está claro se os pulsares reciclados de alta massa ($\sim 1,8 M_\odot$) nasceram massivos ou cresceram via acreção. Além disso, dois mecanismos físicos concorrentes preveem correlações opostas entre a massa da EN ($M$) e o período de rotação ($P$) em populações recicladas:

1. Reciclagem impulsionada por acreção: O ganho de momento angular é acoplado ao ganho de massa, prevendo uma correlação negativa ($\rho < 0$; rotações mais rápidas associadas a massas maiores).
2. Mecanismo do Momento de Inércia (MoI): ENs de menor massa possuem menores momentos de inércia e, portanto, giram mais rapidamente para um dado momento angular acretado, prevendo uma correlação positiva ($\rho > 0$).
   As amostras atuais são limitadas por incertezas de medição e efeitos de seleção, tornando difícil distinguir esses mecanismos ou determinar se a distribuição de massa reciclada é intrinsecamente bimodal ou uma população unimodal assimétrica.

Metodologia
Os autores apresentam uma análise bayesiana hierárquica de $\sim 50$ pulsares de rádio binários galácticos com massas de EN medidas, classificados em sistemas Pulsar–Anã Branca (PSR–WD) e Estrela de Nêutrons Dupla (DNS).

• Curadoria de Dados: A amostra é compilada de catálogos estabelecidos de massa de pulsar. Medições repetidas são combinadas via ponderação de variância inversa, com cautela aplicada a sistemas com medições discrepantes (ex: J0348+0432, B1855+09).
• Classificação: Os sistemas são classificados pelo tipo de binária (primárias PSR–WD, primárias DNS, secundárias DNS) em vez do período de rotação para evitar circularidade. Companheiras de anã branca são subclassificadas como Hélio (He) ou Carbono-Oxigênio/Oxigênio-Neon (CO/ONe) com base em um limiar de massa de $0,45 M_\odot$.
• Estrutura Estatística: Um modelo gaussiano bivariado hierárquico é empregado para estimar correlações intrínsecas entre massa e proxies de reciclagem (período de rotação, período orbital, excentricidade, inclinação). Esta abordagem modela explicitamente os erros de medição como variáveis latentes para corrigir o viés de diluição de regressão (atenuação).
• Correções de Seleção: A análise incorpora correções para detectabilidade Doppler-orbital (usando o fator de eficiência $\gamma_1$) e considera vieses de inclinação inerentes às medições de massa por atraso de Shapiro.
• Comparação de Modelos: Para testar o mecanismo de MoI, os autores comparam modelos de regressão do período de rotação contra a massa diretamente versus contra o momento de inércia (derivado de relações universais), utilizando o Critério de Informação Bayesiano (BIC).

Contribuições Principais e Resultados

1. Anti-correlação Massa–Período de Rotação:
   O estudo encontra evidência moderada de uma anti-correlação entre a massa e o período de rotação na população reciclada agrupada ($N=47$ primárias).

   • Resultado: $\rho = -0,26$ com um intervalo de credibilidade de 90% de $[-0,47, -0,01]$ e uma probabilidade posterior unilateral $P(\rho < 0) = 0,96$ (aprox. $1,7\sigma$).
   • Robustez: O resultado se mantém ao incluir DNS candidatas e aplicar correções de seleção.
   • Nuance de Subpopulação: A amostra PSR–WD sozinha mostra uma tendência sugestiva ($\rho = -0,19$), enquanto as primárias DNS sozinhas pendem levemente para o positivo ($\rho = +0,13$), embora a amostra DNS seja pequena demais ($N=10$) para ser conclusiva. O significado do agrupamento é impulsionado principalmente pela inclusão de primárias DNS, que ocupam uma região de baixa massa e período longo, estendendo a diagonal negativa no plano massa–rotação.

2. Degenerescência de Mecanismos:
   Embora a correlação negativa seja consistente com a reciclagem impulsionada por acreção, os autores não conseguem distinguir estatisticamente isso do mecanismo de MoI.

   • Razão: Na faixa de massa observada, o momento de inércia $I(M)$ da EN é quase linear. Consequentemente, uma regressão do período de rotação sobre a massa é estatisticamente indistinguível de uma regressão sobre o $I(M)$ (diferença de BIC $\Delta \text{BIC} = -0,2$).

3. Distribuições de Massa e Bimodalidade:

   • Reciclados vs. Não Reciclados: Primárias recicladas são significativamente mais massivas ($1,50 M_\odot$) do que secundárias DNS não recicladas ($1,30 M_\odot$), consistente com achados anteriores.
   • Bimodalidade: Dentro da população reciclada, uma mistura de duas gaussianas é preferida em relação a uma gaussiana simétrica única, mas uma distribuição única com assimetria à direita (skew-normal) fornece a descrição mais parcimoniosa (menor BIC). Isso sugere que a aparente bimodalidade pode ser um artefato de uma população única com uma cauda de alta massa, em vez de duas populações de nascimento distintas.

4. Tipo e Massa da Companheira:
   Estrelas de nêutrons com companheiras de anã branca de Hélio são marginalmente mais massivas ($1,547 M_\odot$) do que aquelas com companheiras CO/ONe ($1,484 M_\odot$), uma diferença de $\sim 0,06 M_\odot$. Esta diferença é de significância estatística marginal ($p=0,11$), mas alinha-se com a expectativa de que os canais de WD de He envolvem transferência de massa mais extensa.

5. Correlação Massa da Companheira–Excentricidade em DNS:
   Em um framework hierárquico, os autores confirmam uma forte correlação positiva entre a massa da companheira (segunda nascida) e a excentricidade orbital em sistemas DNS ($\rho = +0,82$). Isso apoia o elo teórico entre supernovas de baixo chute (produzindo sistemas de baixa massa e baixa excentricidade) e a formação de DNSs.

Significância e Alegações
O artigo enquadra seus achados como evidência moderada ou sugestiva em vez de uma detecção definitiva, primariamente devido ao tamanho modesto da amostra ($N \approx 50$).

• Interpretação: Os resultados apoiam um quadro de "dois canais": sistemas PSR–WD são "crescidos por acreção", exibindo a anti-correlação massa–rotação e desvios de massa, enquanto sistemas DNS são "dominados pela massa de nascimento", onde as massas refletem condições de formação em vez de história de reciclagem.
• Limitações: A incapacidade de distinguir entre os mecanismos de acreção e MoI é uma consequência direta da relação quase linear entre massa e momento de inércia na faixa observada.
• Perspectivas Futuras: Os autores preveem que alcançar uma detecção confiável de $3\sigma$ da correlação massa–rotação exigirá uma amostra de $N \gtrsim 125–150$ sistemas reciclados de massa precisa, provavelmente um feito da era SKA. Eles observam que quebrar a degenerescência acreção/MoI pode exigir a medição do spin e da massa de ENs não recicladas (ex: segundos pulsares nascidos em sistemas DNS) para isolar propriedades de nascimento de efeitos de reciclagem.

O trabalho reformula o debate sobre as distribuições de massa recicladas, sugerindo que a cauda de alta massa observada é melhor descrita por uma distribuição unimodal assimétrica do que por uma população de nascimento distinta e bimodal.