Relativistic jets from millisecond proto-magnetars

Resumo Técnico: Jatos Relativísticos de Proto-Magnetares de Milissegundo

Definição do Problema
Proto-magnetares de milissegundo — estrelas de nêutrons que rotacionam rapidamente ($P \sim 1$ ms) e possuem campos magnéticos intensos ($B \gtrsim 10^{15}$ G) formados em supernovas de colapso de núcleo, fusões de estrelas de nêutrons ou colapso induzido por acreção de anã branca (AIC) — têm sido há muito tempo propostos como motores centrais para surtos de raios gama (GRBs). Um desafio teórico crítico permanece: durante os primeiros segundos após o nascimento, esses objetos emitem fluxos intensos de neutrinos que aquecem a superfície e impulsionam ventos ricos em bárions. Este "vento de neutrinos" carrega o fluxo com matéria, potencialmente suprimindo a magnetização ($\sigma$) e limitando o fator de Lorentz de volume ($\Gamma_\infty$) alcançável. Modelos semianalíticos anteriores, calibrados para simulações unidimensionais e assumindo geometrias de monopolo dividido com média angular, previram $\sigma \lesssim 1$ nos tempos iniciais, sugerindo que os proto-magnetares não poderiam lançar os jatos ultra-relativísticos ($\Gamma \gtrsim 10-100$) necessários para alimentar GRBs de curta duração. A questão de se efeitos multidimensionais poderiam superar esse carregamento de bárions para produzir jatos relativísticos dentro de segundos da formação ainda não havia sido totalmente abordada.

Metodologia
Os autores apresentam as primeiras simulações de magnetohidrodinâmica de relatividade geral (GRMHD) tridimensionais de ventos de proto-magnetares aquecidos por neutrinos. O estudo utiliza uma versão modificada do código GRHydro de código aberto (parte do Einstein Toolkit) empregando um esquema de volume finito com reconstrução WENO-Z e um solver de Riemann HLLE aproximado. Características numéricas principais incluem:

• Transporte de Neutrinos: Uma aproximação de um momento ($M_0$) da equação de Boltzmann de relatividade geral é usada para modelar de forma autossistente o aquecimento e o resfriamento por neutrinos, incluindo processos $\beta$ de corrente carregada, aniquilação de pares e espalhamento.
• Evolução do Campo Magnético: O campo magnético é evoluído via o potencial magnético de 4-vetor na gauge de Lorenz generalizada, utilizando um esquema de transporte de restrição upwind para manter a restrição de solenoidal com precisão de máquina.
• Condições Iniciais: Os modelos assumem uma estrela de nêutrons de $1.4 M_\odot$ com uma temperatura central de $20$ MeV. As simulações abrangem uma gama de períodos de rotação iniciais ($P \approx 1 - 18$ ms) e forças de campo magnético polar ($B_P \approx 10^{13} - 3 \times 10^{15}$ G).
• Configuração: As simulações utilizam simetria de reflexão através do plano equatorial e utilizam refinamento de malha adaptativo (Berger-Oliger) com a grade mais interna resolvendo a superfície do proto-magnetar. Os modelos são evoluídos primeiro para um estado estacionário hidrodinâmico impulsionado por neutrinos antes de inicializar um campo dipolar poloidal.

Principais Contribuições e Resultados
As simulações revelam que efeitos multidimensionais alteram fundamentalmente a estrutura do vento em comparação com modelos 1D ou de média angular, particularmente para modelos de rotação rápida ($P \approx 1$ ms) e alta magnetização ($B_P \approx 3 \times 10^{15}$ G):

1. Estratificação Angular do Fluxo: O fluxo não é uniforme, mas divide-se em três regiões distintas:

   • Cinturão Equatorial ($\theta \gtrsim 30^\circ$): Forças centrífugas aumentam significativamente a perda de massa perto do equador rotacional. Isso cria um vento denso e sub-relativístico ($v \sim 0.1c$) com alto carregamento de bárions.
   • Jato Polar ($\theta \lesssim 10^\circ$): As linhas de campo polar experimentam um aumento mínimo de carregamento de massa por força centrífuga. Consequentemente, esta região permanece "limpa" (baixo carregamento de bárions) e altamente magnetizada.
   • Região de Transição: Latitudes intermediárias onde o vento acelera mais eficientemente.

2. Colimação Intrínseca do Jato: O vento equatorial denso e sub-relativístico atua como um meio de confinamento, colimando naturalmente o fluxo polar com menos carga de bárions. Esta autocolimação ocorre em raios pequenos ($\lesssim 100$ km) sem exigir o confinamento por ejectas de supernova ou de fusão externa.

3. Alcance de Condições Ultra-Relativísticas: Nos modelos de rotação e magnetização mais rápidas (ex: modelo B3P1), o jato polar atinge valores de magnetização de pico de $\sigma \sim 30 - 100$. Isso é suficiente para atingir fatores de Lorentz de volume $\Gamma_\infty \sim 100$ em escalas maiores. O ângulo de abertura do jato é de aproximadamente $8^\circ$, consistente com estimativas observacionais para GRBs curtos e longos.

4. Particionamento de Energia: A energia total do vento é particionada entre um jato polar ultra-relativístico e um vento equatorial de ângulo amplo e sub-relativístico. Este último carrega a maior parte da energia cinética e da massa, potencialmente alimentando a ejecta da supernova ou kilonova, enquanto o primeiro alimenta o GRB.

Significância e Alegações
O artigo afirma demonstrar que proto-magnetares de milissegundo podem lançar jatos relativísticos dentro de segundos de sua formação, resolvendo a tensão de longa data entre o carregamento de bárions impulsionado por neutrinos e os requisitos para motores de GRB. Os autores argumentam que:

• A Multidimensionalidade é Essencial: O carregamento de bárions não impede a formação de jatos relativísticos quando a dependência latitudinal das forças centrífugas e do aquecimento por neutrinos é levada em conta.
• Cenário de Motor Unificado: Um único motor central pode alimentar simultaneamente um GRB de curta duração (via o jato polar) e a supernova/kilonova associada (via o vento equatorial), explicando a partição observada entre as energias de GRB corrigidas pelo feixe e as energias de ejecta.
• Consistência Temporal: O mecanismo opera nos primeiros segundos da vida do proto-magnetar, consistente com a fase de emissão prompt de GRBs curtos.

Os autores observam que, embora seus modelos assumam campos dipolares ordenados e luminosidades de neutrinos fixas, os resultados sugerem que o surgimento de jatos relativísticos é uma consequência genérica de rotação e magnetização extremas nos estágios iniciais da evolução da proto-estrela de nêutrons. Eles reconhecem que trabalhos futuros devem abordar campos magnéticos gerados por dínamo de forma autossistente, resfriamento/spin-down dependente do tempo e a interação com discos de acreção massivos que circundam remanescentes de fusão.