Born to be recycled: a comprehensive population synthesis of the Galactic millisecond pulsars

1. Problema e Contexto

Os pulsares de milissegundo (MSPs) são os objetos mais antigos, mas também os mais rápidos, conhecidos na população de pulsares. A hipótese predominante para sua formação é o mecanismo de "reciclagem", onde um pulsar normal em um sistema binário acreta matéria de uma companheira, ganhando momento angular e acelerando sua rotação.

Apesar de ser um cenário amplamente aceito, os detalhes físicos da reciclagem (como o decaimento do campo magnético, a evolução do binário e a geometria de emissão) ainda não são completamente compreendidos. Estudos anteriores de síntese populacional frequentemente simplificaram a evolução, ignorando a evolução completa do binário desde o nascimento ou não considerando a detecção simultânea em rádio e raios gama. Além disso, questões como a origem do excesso de raios gama (GeV) no centro galáctico e a contribuição de canais de formação alternativos (como o colapso induzido por acreção de anãs brancas) permanecem abertas.

2. Metodologia

Os autores desenvolveram um algoritmo de síntese populacional abrangente para pulsares em sistemas binários, integrando evolução estelar, evolução do pulsar e modelos de detecção observacional.

• Evolução do Binário: Utilizaram o código SEVN (Stellar EVolution for N-body) para simular a evolução estelar e as interações binárias (como overflow de lóbulo de Roche e transferência de massa). As estrelas secundárias começam como estrelas da sequência principal.
• Evolução do Pulsar: Diferente do SEVN padrão, os autores implementaram um modelo autoconsistente para a evolução do pulsar que inclui:
  • Magnetosfera force-free (sem forças).
  • Decaimento do campo magnético (seguindo uma lei de potência com três escalas de tempo características: 0,05, 1,4 e 2,3 Myr).
  • Frenagem gravitacional (ondas gravitacionais) para rotações muito rápidas.
  • Torques de acreção, frenagem magnética e perda de radiação.
• Geometria de Emissão:
  • Rádio: Modelado pela geometria do polar cap (tampa polar).
  • Raios Gama: Modelado pelo cenário de striped wind (vento listrado), onde a emissão ocorre na folha de corrente do vento do pulsar.
• Modelo de Detecção:
  • Rádio: Simulação baseada nos levantamentos FAST GPPS e PMPS (Parkes), calculando fluxo, largura de pulso e relação sinal-ruído (S/N).
  • Raios Gama: Comparação com o catálogo 3PC do telescópio Fermi-LAT, utilizando mapas de sensibilidade all-sky.
• Parâmetros Iniciais: Incluem distribuições de massa (IMF de Kroupa), massa da companheira, excentricidade, ângulo de inclinação e velocidade de kick (baixa dispersão de 70 km/s para manter o sistema ligado).

3. Principais Contribuições

• Simulação Completa: É o primeiro estudo a seguir a evolução completa de um binário desde o nascimento da estrela de nêutrons até o estado de MSP, incluindo a fase de acreção e o subsequente resfriamento, aplicando critérios de detecção realistas em rádio e raios gama simultaneamente.
• Modelo de Campo Magnético Variável: Incorpora um decaimento do campo magnético não uniforme, permitindo que diferentes pulsares sigam trajetórias distintas no diagrama $P-\dot{P}$.
• Análise de Geometria: Estuda detalhadamente como os ângulos de visão ($\zeta$) e inclinação magnética ($\chi$) afetam a detecção, especialmente a dificuldade de observar pulsares com eixos de rotação e magnéticos alinhados.
• Investigação de Canais Alternativos: Avaliou a contribuição do Colapso Induzido por Acreção (AIC) de anãs brancas de oxigênio-neônio (ONeWDs) para a população de pulsares levemente reciclados.

4. Resultados Chave

• Reprodução da População: O modelo reproduz bem a distribuição observada de pulsares de rádio e raios gama no diagrama $P-\dot{P}$, validado por testes de Kolmogorov-Smirnov (KS).
• Pulsares Não Identificados: Estima-se que existam menos de 220 pulsares não identificados no catálogo 4FGL do Fermi-LAT, uma redução significativa em relação a previsões anteriores.
• Duração da Acreção: A maioria dos pulsares acreta matéria por durações típicas de $10^8$ a $10^9$ anos.
• Alinhamento de Eixos: O tempo de alinhamento entre o eixo de rotação e a normal ao plano orbital é curto. Consequentemente, cerca de 80% da população de pulsares binários tem o ângulo $\alpha$ próximo de 0° (alinhados).
• Massa dos Pulsares: A massa média dos MSPs reciclados é de aproximadamente 1,8 $M_\odot$, com alguns atingindo até 2,7 $M_\odot$. Isso corrobora observações recentes de sistemas "spider" (pulsares que ablam suas companheiras).
• Origem do Excesso GeV: Pulsares nascidos nos braços espirais da Galáxia contribuem com apenas ~5% do excesso de raios gama no centro galáctico (cerca de 190 MSPs). Isso reforça a hipótese de que, se o excesso for inteiramente de origem de pulsares, a maioria deve vir de aglomerados globulares que migraram para o centro ou de nascimentos no bojo galáctico.
• Geometria de Detecção: É difícil observar pulsares reciclados com eixos de rotação e magnéticos alinhados ($\chi \le 10^\circ$). Além disso, altos ângulos de visão ($\zeta$) parecem necessários para a detecção.
• Contribuição de ONeWDs: Apenas uma fração minúscula (~0,0076%) das anãs brancas ONe em sistemas binários contribui para a população de pulsares levemente reciclados via AIC. No entanto, esse canal pode explicar parte dos pulsares com períodos mais longos que não são totalmente explicados pela acreção em estrelas de nêutrons nascidas de supernovas.
• População de "Spiders": O modelo consegue reproduzir a existência de sistemas tipo "black widow" e "redback", mas subestima o número total observado, sugerindo a necessidade de outros canais de formação ou ajustes na eficiência de acreção/evaporação.

5. Significância e Conclusões

Este trabalho fornece uma das simulações mais realistas da evolução de pulsares reciclados até o momento. Ao integrar a evolução binária detalhada com modelos de emissão e detecção, os autores validam o cenário de reciclagem como o principal mecanismo de formação de MSPs.

Os resultados têm implicações importantes para:

1. Física de Estrelas de Nêutrons: Suportam a existência de estrelas de nêutrons com massas superiores a 2 $M_\odot$, restringindo a equação de estado da matéria nuclear.
2. Astrofísica Galáctica: Clarificam a origem do excesso de raios gama no centro da Via Láctea, apontando para uma origem em aglomerados globulares ou no bojo, e não apenas em estrelas do campo galáctico.
3. Futuros Levantamentos: Preveem que o Square Kilometre Array (SKA) aumentará a detecção de pulsares reciclados em rádio em cerca de 3 vezes, e um instrumento de raios gama 10 vezes mais sensível que o Fermi/LAT aumentaria as detecções em raios gama em 5 vezes.

Em suma, o estudo confirma a eficácia do mecanismo de reciclagem, mas destaca que a compreensão completa da população de MSPs, especialmente os sistemas "spider" e os pulsares de raios gama silenciosos, requer a consideração de canais de formação alternativos e refinamentos nos parâmetros de evolução binária.