Evolution of proto-neutron stars to pulsars, magnetars and central compact objects

Resumo Técnico: Evolução de Proto-Estrelas de Nêutrons para Pulsares, Magnetares e Objetos Compactos Centrais

Definição do Problema
A origem dos campos magnéticos ultrafortes observados em magnetares ($B \sim 10^{14}–10^{15}$ G) permanece um objeto de debate. Embora a hipótese do "campo fóssil" (conservação de fluxo das estrelas progenitoras) possa explicar os campos de estrelas de nêutrons convencionais e de Objetos Compactos Centrais (CCOs), ela tem dificuldade em explicar os campos de magnetares sem invocar condições progenitoras extremas e raras. Por outro lado, a hipótese do dínamo, que requer rotação rápida e convecção vigorosa durante a fase de proto-estrela de nêutrons (PNS), é amplamente aceita para magnetares, mas enfrenta um desafio: a taxa de nascimento inferida de magnetares é comparável à taxa de supernovas de colapso de núcleo, sugerindo que a formação de magnetares não é tão rara quanto os requisitos de parâmetros extremos implicariam. Além disso, cálculos padrão de conservação de fluxo sugerem que, mesmo com o colapso de um núcleo de ferro, os campos resultantes ($B \sim 10^{10}$ G) são insuficientes para explicar os campos típicos de pulsares ($B \sim 10^{12}$ G), implicando a necessidade de um mecanismo de dínamo mesmo para pulsares comuns. O problema central é reconciliar as altas taxas de nascimento de vários tipos de estrelas de nêutrons com suas diversas intensidades de campo magnético usando um arcabouço evolutivo unificado.

Metodologia
Os autores empregam um modelo de dínamo $\alpha-\Omega$ reduzido, unidimensional e impulsionado por cisalhamento, adaptado de Wickramasinghe et al. (2014), que originalmente o aplicaram à fusão de anãs brancas. O modelo rastreia a evolução dos campos magnéticos poloidal ($B_p$) e toroidal ($B_\phi$) durante a fase de PNS. As características principais da metodologia incluem:

• Equações Governantes: O modelo resolve equações diferenciais acopladas para $B_p$ e $B_\phi$. O campo toroidal é gerado via o efeito $\Omega$ (rotação diferencial atuando sobre o campo poloidal), enquanto o campo poloidal é regenerado via o efeito $\alpha$ (convecção turbulenta atuando sobre o campo toroidal).
• Evolução da PNS: Diferente de modelos estáticos, este arcabouço incorpora a contração dinâmica da PNS de um raio inicial $R_0 \sim 40$ km para um raio final $R_\infty \approx 12$ km. A evolução do raio é parametrizada com base em Camelio et al. (2017) para várias massas e equações de estado (EoS).
• Massa e Momento Angular: O modelo contabiliza a diminuição da massa gravitacional devido ao resfriamento por neutrinos e a conservação do momento angular, o que impulsiona o aumento da velocidade de rotação (spin-up) da PNS.
• Viscosidade e Amortecimento de Cisalhamento: O amortecimento da rotação diferencial (taxa de cisalhamento $q$) é modelado via uma combinação de viscosidade impulsionada pela instabilidade magnetorotacional (MRI) e viscosidade convectiva. A taxa de cisalhamento decai ao longo do tempo, terminando eventualmente o dínamo.
• Perda de Fluxo e Estabilidade: O modelo inclui um termo de reconexão para contabilizar a perda de fluxo magnético (definida inicialmente em 20%) e incorpora critérios de estabilidade (Braithwaite 2009), onde instabilidades de campo determinam as escalas de tempo de decaimento de $B_p$ e $B_\phi$.
• Varredura do Espaço de Parâmetros: Os autores integram numericamente as equações usando o método Runge-Kutta, variando sistematicamente os parâmetros iniciais, incluindo a taxa de cisalhamento ($q_0$), massa inicial ($M_0$), raio inicial ($R_0$), campos semente iniciais e o período rotacional final ($P_\infty$).

Resultos Principais
As simulações revelam que a configuração final do campo magnético é altamente sensível à taxa de rotação inicial e ao cisalhamento, permitendo que um único mecanismo de dínamo produza todo o espectro de tipos observados de estrelas de nêutrons:

1. Limiar de Taxa de Cisalhamento: O valor de saturação do campo toroidal ($B_{\phi, \infty}$) depende fortemente da taxa de cisalhamento inicial $q_0$. Existe um limiar crítico em $q_0 \simeq 0,23$.

   • Abaixo do Limiar ($q_0 < 0,23$): O efeito $\alpha$ não é suficientemente desencadeado. O campo poloidal permanece próximo ao valor ditado pela conservação de fluxo ($B_p \approx 2,5 \times 10^{10}$ G), e o campo toroidal não cresce significativamente. Este regime corresponde aos Objetos Compactos Centrais (CCOs).
   • Acima do Limiar: O dínamo torna-se ativo. Taxas de cisalhamento mais altas levam a campos toroidais mais fortes.

2. Formação de Diferentes Populações:

   • Magnetares: Para os períodos iniciais mais curtos (rotação rápida, ex: $P_\infty \sim 2,5$ ms), o modelo produz intensidades do tipo magnetar com $B_p \simeq 10^{15}$ G e $B_\phi \simeq 10^{16}$ G.
   • Pulsares Ordinários: Para rotadores iniciais mais lentos, o modelo produz campos típicos de pulsares ($B_p \simeq 10^{12}$ G, $B_\phi \simeq 10^{14}$ G).
   • Magnetares de Baixo-B / Pulsares de Alto-B: Taxas de rotação intermediárias produzem campos na faixa de magnetares de baixo-B ou pulsares de alto campo.

3. Papel da Conservação de Fluxo: A intensidade do campo poloidal inicial é determinada em grande parte pela conservação de fluxo durante o colapso. No entanto, o processo de dínamo (especificamente o efeito $\alpha$) é necessário para amplificar esses campos aos níveis observados de pulsares e magnetares. O campo toroidal é genericamente dominado pelo efeito $\Omega$.

4. Objetos Compactos Centrais (CCOs): O modelo sugere que os CCOs não passam por um processo de dínamo vigoroso devido à sua rotação inicial lenta ($P_0 \sim 100$ ms). Seus campos poloidais são definidos pela conservação de fluxo, mas eles ainda podem adquirir campos toroidais amplificados apenas pelo efeito $\Omega$, explicando seus pequenos campos dipolares, mas potencialmente significativos componentes toroidais internos.

Significância e Alegações
O artigo afirma reconciliar a aparente discrepância entre as altas taxas de nascimento de magnetares e as condições extremas frequentemente exigidas pelos cenários de formação. Ao demonstrar que um processo de dínamo operando durante a fase de PNS pode naturalmente produzir um contínuo de intensidades de campo magnético dependendo da taxa de rotação inicial, os autores propõem que:

• Tanto os pulsares comuns quanto os magnetares provavelmente passam por uma fase de dínamo; a diferença reside na eficiência deste processo impulsionado pela rotação inicial.
• A hipótese do "campo fóssil" é insuficiente para explicar os campos de pulsares típicos, necessitando de um dínamo mesmo para estrelas de nêutrons não-magnetares.
• Os CCOs representam o final de baixa rotação deste espectro, onde o dínamo é suprimido, deixando campos dominados pela conservação de fluxo.

Os autores enfatizam que seu modelo unidimensional, embora simplificado, captura as dependências essenciais de rotação, cisalhamento e flutuabilidade, e se alinha qualitativamente com simulações tridimensionais mais complexas (ex: Raynaud et al. 2020). Eles concluem que o dínamo $\alpha-\Omega$ fornece um arcabouço unificado para compreender a diversidade dos campos magnéticos de estrelas de nêutrons jovens sem exigir condições de progenitores raras e extremas para magnetares.