Resolution Dependence in Magnetohydrodynamic Simulations of Neutrino-Driven Core-Collapse Supernovae

Título: Dependência de Resolução em Simulações Magnetohidrodinâmicas de Supernovas de Colapso do Núcleo Acionadas por Neutrinos

1. Problema e Contexto

As supernovas de colapso do núcleo (CCSNe) são os eventos finais de estrelas massivas. O mecanismo de explosão mais estudado é o acionado por neutrinos, onde a reabsorção de neutrinos aquece o material atrás da onda de choque estagnada, permitindo sua reativação. Embora simulações hidrodinâmicas puras tenham avançado, a inclusão de campos magnéticos em simulações 3D de neutrinos permanece um desafio.

O problema central abordado neste estudo é a dependência da resolução numérica em simulações de supernovas magnetohidrodinâmicas (MHD) 3D. Especificamente, os autores investigam como a resolução da grade computacional e a força inicial do campo magnético afetam:

• A reativação da onda de choque (shock revival).
• A dinâmica da explosão e a energia final.
• As propriedades do remanescente compacto (estrela de nêutrons proto-neutrônica - PNS).
• A emissão de neutrinos e ondas gravitacionais.

Existe a necessidade de entender se campos magnéticos fortes (da ordem de $10^{12}$ G), possivelmente gerados por dinamos de pequena escala ou fossilizados, podem influenciar explosões em progenitores não rotativos ou de rotação lenta, e se as simulações atuais estão suficientemente resolvidas para capturar esses efeitos.

2. Metodologia

Os autores realizaram uma série de simulações 3D utilizando o código CoCoNuT-FMT (uma versão Newtoniana MHD com transporte de neutrinos de múltiplos grupos).

• Progenitor: Uma estrela de $13 M_\odot$ não rotativa e magnetizada, gerada pelo código MESA.
• Configurações de Campo Magnético: Dois cenários de campo dipolar inicial no centro da estrela:
  • Campo forte: $B \approx 10^{12}$ G.
  • Campo fraco: $B \approx 10^{10}$ G.
• Variação de Resolução: Foram executados seis modelos com diferentes resoluções de grade ($N_r \times N_\theta \times N_\phi$):
  • vLRes: Baixa resolução (275 x 64 x 128).
  • LRes: Média resolução (550 x 64 x 128).
  • MRes: Alta resolução (550 x 128 x 256).
  • HRes: Muito alta resolução (550 x 256 x 512).
  • Modelos com campo fraco (vLRes10 e MRes10) foram comparados com seus equivalentes de campo forte.
• Física: Inclui equações de estado SFHo, transporte de neutrinos FMT (Fast Multi-group Transport) e tratamento de campos magnéticos com limpeza hiperbólica de divergência.

3. Contribuições Principais

• Primeiro Estudo de Resolução MHD 3D: Este é o primeiro estudo sistemático de resolução focado em supernovas acionadas por neutrinos com magnetohidrodinâmica completa em 3D.
• Análise do Dynamo de Pequena Escala: Demonstra como a resolução afeta a eficiência do dynamo de pequena escala (SSD) na região de ganho e na PNS, crucial para a amplificação do campo magnético em estrelas não rotativas.
• Acoplamento Campo Magnético-Neutrino: Revela um mecanismo indireto onde campos magnéticos fortes, amplificados em altas resoluções, deformam a PNS, alterando a emissão de neutrinos e, consequentemente, a dinâmica da explosão.

4. Resultados Chave

A. Dinâmica da Explosão e Reativação do Choque:

• O tempo de reativação do choque é independente da resolução e da força inicial do campo magnético para todos os modelos (ocorrendo entre 0,15s e 0,20s após o bounce).
• No entanto, a energia de explosão e o raio do choque pós-reativação são fortemente dependentes da resolução. Modelos de maior resolução (HRes) expandem mais rápido e atingem energias maiores.
• Campos magnéticos fortes não impulsionam a explosão diretamente (não causam a reativação), mas têm um impacto subsidiário: modelos com campos fortes e alta resolução exibem energias de explosão maiores do que seus equivalentes com campos fracos.
• A energia turbulenta cinética é convertida eficientemente em energia magnética na região de ganho, em vez de aumentar a pressão turbulenta hidrodinâmica (diferente do que ocorre em simulações puramente hidrodinâmicas).

B. Emissão de Neutrinos:

• Há uma correlação inversa entre a resolução e a luminosidade/energia média dos neutrinos. Modelos de maior resolução (com campos magnéticos mais fortes devido ao dynamo mais eficiente) exibem luminosidades de neutrinos ~5-8% menores.
• Mecanismo: Campos magnéticos fortes deformam a PNS, tornando-a mais alongada nos polos. Isso faz com que parte dos neutrinos seja emitida de regiões de raio maior e temperatura mais baixa, reduzindo a luminosidade total.
• Modelos de baixa resolução (vLRes) sofrem de alta difusão numérica, impedindo a sustentação de campos magnéticos fortes, resultando em propriedades de neutrinos semelhantes independentemente da força do campo inicial.

C. Propriedades do Remanescente (PNS):

• Convecção e Tensões: Campos magnéticos fortes e alta resolução levam a uma zona de convecção da PNS mais extensa e vigorosa. As tensões magnéticas (Maxwell) tornam-se comparáveis às tensões cinéticas (Reynolds) nas regiões internas da PNS.
• Rotação e Kick:
  • Torques magnéticos dominam o spin-up (aceleração da rotação) da estrela de nêutrons em modelos com campos fortes.
  • Períodos de rotação estimados variam de 0,024 s (MRes, campo forte) a 0,178 s (vLRes, campo forte).
  • O modelo MRes10 (campo inicial fraco, alta resolução) também experimenta spin-up significativo devido à amplificação do campo, atingindo um período de 8,12 s, enquanto o vLRes10 (baixa resolução) não.
  • As velocidades de kick (recuo) são baixas (10-71 km/s) e dependem da assimetria global da explosão, não apenas da energia total.

D. Ondas Gravitacionais:

• A inclusão de campos fortes e alta resolução altera qualitativamente o espectro de ondas gravitacionais, aumentando a potência e a faixa de frequências.
• Um sinal de modo g do núcleo (identificado anteriormente em 2D) foi observado no modelo MRes, sugerindo que efeitos MHD podem excitar esse modo em 3D.

5. Significado e Conclusões

O estudo conclui que as simulações de supernovas MHD não atingiram a convergência mesmo na resolução mais alta testada (0,7 graus angulares), diferentemente de simulações puramente hidrodinâmicas que convergem em ~1 grau.

• Sensibilidade à Resolução: A sub-resolução de modelos MHD leva a resultados quantitativamente diferentes, pois a difusão numérica impede a amplificação correta dos campos magnéticos via dynamo de pequena escala.
• Impacto Físico: Campos magnéticos fortes, mesmo em progenitores não rotativos, podem ser amplificados para níveis dinamicamente relevantes ($>10^{14}$ G na superfície da PNS), transformando o remanescente em um magnetar. Eles alteram a estrutura da PNS, a emissão de neutrinos e a redistribuição de momento angular.
• Implicações Futuras: É necessário realizar estudos de resolução de longo prazo e comparações de códigos independentes para quantificar as incertezas e entender plenamente a conexão entre o progenitor e o remanescente em cenários magnetizados.

Em resumo, a resolução numérica é um fator crítico que determina a força do campo magnético amplificado, o que, por sua vez, retroage na estrutura da estrela de nêutrons, na emissão de neutrinos e na dinâmica final da explosão.