Dynamo Simulations Confirm Predominantly Toroidal Fields in Near-Core Region of an Intermediate-Mass Star

Resumo Técnico: Simulações de Dínamo Confirmam Campos Predominantemente Toroidais na Região Próxima ao Núcleo de uma Estrela de Massa Intermediária

Problema
A recente análise asterosísmica da estrela de sequência principal do tipo F, KIC 9244992 (Takata et al. 2025), forneceu as primeiras restrições sobre a geometria do campo magnético interno, inferindo um campo magnético predominantemente toroidal forte (componente azimutal $\sim 92$ kG) no interior radiativo profundo, especificamente na região próxima ao núcleo definida pelo pico da frequência de Brunt–Väisälä. Este achado apresenta um desafio teórico: a mesma estrela exibe uma rotação diferencial notavelmente plana entre a região próxima ao núcleo e a região próxima à superfície (Saio et al. 2015), com apenas uma diferença de $\sim 4\%$. A teoria de dínamo padrão sugere que campos toroidais fortes requerem cisalhamento significativo, levantando a questão de como campos tão fortes podem ser mantidos na presença de uma rotação diferencial global fraca. Trabalhos numéricos anteriores (ex., Ratnasingam et al. 2024) sugeriram que camadas de cisalhamento local poderiam amplificar campos toroidais, mas faltava uma simulação direta que correspondesse à massa específica e às restrições asterossísmicas de KIC 9244992.

Metodologia
Os autores realizaram simulações magnetohidrodinâmicas (MHD) tridimensionais e anelásticas utilizando o código SPIN (Edelmann et al. 2019) para modelar uma estrela de sequência principal de $2 M_\odot$ com uma fração de massa de hidrogênio central de $X_c = 0,5$. As simulações resolvem as equações MHD na aproximação anelástica com uma resolução de 1500 pontos de grade radial e harmônicos esféricos $\ell_{max} = 84$.

Os principais parâmetros de configuração incluem:

• Estado Inicial: Um modelo estelar 1D gerado pelo MESA, com uma velocidade convectiva "potencializada" (aproximadamente $10\times$ a teoria de comprimento de mistura) para alcançar números de Reynolds computacionalmente viáveis, mantendo um número de Rossby ($Ro \sim 0,1-0,3$) comparável ao de KIC 9244992.
• Rotação: A estrela foi inicializada com rotação de corpo rígido ($\Omega \approx 1,0 \times 10^{-5}$ rad s$^{-1}$), aproximadamente 10 vezes mais rápida que KIC 9244992, para compensar a convecção potencializada.
• Campo Magnético: Um campo semente dipolo ($100$ kG na fronteira inferior) foi imposto uma única vez após o alcance de um estado hidrodinâmico estacionário.
• Parâmetros: A simulação principal utilizou um número de Reynolds ($Re$) de 240, número de Reynolds magnético ($Rm$) de 960 e número de Prandtl magnético ($Pm$) de $\sim 4$.
• Análise: Os dados foram analisados do dia 40 ao 100 do tempo físico, uma vez que a energia magnética atingiu um estado estacionário estatístico. O estudo focou na região próxima ao núcleo (logo fora da fronteira convectivo-radiativa, CRB), onde a frequência de Brunt–Väisälä atinge o pico.

Principais Resultados

1. Geometria e Força do Campo: As simulações geraram de forma autossuficiente campos magnéticos onde a componente azimutal (toroidal) domina a componente radial (poloidal) na região próxima ao núcleo. No pico da frequência de Brunt–Väisälä, a razão entre o campo azimutal e o radial de raiz quadrada média ($B_\phi/B_r$) atingiu um máximo de 20.

   • Na CRB: $B_r \approx 33,8$ kG, $B_\phi \approx 120,8$ kG.
   • No pico do pico de Brunt–Väisälä: $B_r \approx 0,29$ kG, $B_\phi \approx 7,65$ kG.
   • Estes valores estão dentro de um fator de poucos para os valores de força de campo mínima inferidos para KIC 9244992 por Takata et al. (2025).

2. Rotação Diferencial e Cisalhamento: Embora o perfil de rotação com média esférica tenha mostrado apenas uma pequena diferença ($\sim 3,7\%$) entre a região próxima ao núcleo e a superfície — consistente com as inferências asterossísmicas — as simulações revelaram um diferencial de rotação latitudinal significativo. A taxa de rotação variou em até $20\%$ através das latitudes, particularmente perto do equador. Este cisalhamento localizado, em vez do cisalhamento radial global, impulsiona o efeito $\Omega$, amplificando o campo toroidal a partir do componente poloidal.

3. Estrutura Latitudinal: A força do campo toroidal exibiu uma estrutura de pico duplo em latitude, com máximos em latitudes médias e mínimos nos polos. Esta geometria alinha-se com a inferência asterossísmica do sinal dos coeficientes de polinômio de Legendre ($a < 0, b > 0$) relatado por Takata et al. (2025).

4. Robustez: Os autores verificaram que os resultados não são artefatos da força do campo inicial (uma execução com $1/5$ do campo inicial atingiu o mesmo estado estacionário) ou de escolhas específicas de parâmetros, desde que o número de Reynolds seja suficientemente alto ($Re > \sim 100$) e o número de Rossby seja apropriado. Testes de sensibilidade mostraram que reduzir o $Re$ para $\sim 30$ resultou em campos significativamente mais fracos ($1-2$ kG), destacando a necessidade de um $Re$ elevado para a saturação do dínamo.

Significado e Alegações
O artigo afirma que a presença de um campo magnético toroidal dominante na região próxima ao núcleo é uma característica genérica de estrelas com convecção central, surgindo naturalmente de simulações de dínamo globais quando os números adimensionais apropriados ($Re$, $Ro$) são utilizados.

Os autores enfatizam que seus resultados fornecem suporte numérico independente para a interpretação asterossísmica de KIC 9244992. Crucialmente, o estudo demonstra que campos toroidais fortes podem coexistir com uma rotação diferencial de média esférica fraca, porque a geração do campo é impulsionada por cisalhamento latitudinal localizado. Esta descoberta sugere que futuras inferências asterossísmicas não devem assumir uma geometria de campo predominantemente radial em toda a estrela, pois a região próxima ao núcleo é provavelmente dominada por campos toroidais gerados por fluxos de cisalhamento tridimensionais. Este trabalho preenche a lacuna entre as restrições observacionais e a teoria de dínamo numérica, confirmando que a geometria de campo específica inferida para KIC 9244992 é fisicamente plausível dentro do quadro da magnetohidrodinâmica estelar.