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

Resumen Técnico: Las simulaciones de dinamo confirman campos predominantemente toroidales en la región cercana al núcleo de una estrella de masa intermedia

Planteamiento del problema
El reciente análisis asterosísmico de la estrella de tipo F en la secuencia principal KIC 9244992 (Takata et al. 2025) proporcionó las primeras restricciones sobre la geometría del campo magnético interno, infiriendo un campo magnético fuertemente toroidal (componente azimutal $\sim 92$ kG) en el profundo interior radiativo, específicamente en la región cercana al núcleo definida por el pico de la frecuencia de Brunt–Väisälä. Este hallazgo presenta un desafío teórico: la misma estrella exhibe una rotación diferencial notablemente plana entre la región cercana al núcleo y la región cercana a la superficie (Saio et al. 2015), con solo una diferencia del $\sim 4\%$. La teoría de dinamo estándar sugiere que los campos toroidales fuertes requieren un cizallamiento significativo, lo que plantea la cuestión de cómo se pueden mantener campos tan fuertes en presencia de una rotación diferencial global débil. Trabajos numéricos previos (p. ej., Ratnasingam et al. 2024) sugirieron que las capas de cizallamiento local podrían amplificar los campos toroidales, pero faltaba una simulación directa que coincidiera con la masa específica y las restricciones asterosísmicas de KIC 9244992.

Metodología
Los autores realizaron simulaciones magnetohidrodinámicas (MHD) tridimensionales y aneclásticas utilizando el código SPIN (Edelmann et al. 2019) para modelar una estrella de la secuencia principal de $2 M_\odot$ con una fracción de masa de hidrógeno central de $X_c = 0.5$. Las simulaciones resuelven las ecuaciones de MHD bajo la aproximación anelástica con una resolución de 1500 puntos de malla radial y un armónico esférico $\ell_{max} = 84$.

Los parámetros clave de la configuración incluyen:

• Estado inicial: Un modelo estelar 1D generado con MESA, con una velocidad convectiva "potenciada" (aproximadamente $10\times$ la teoría de longitud de mezcla) para lograr números de Reynolds computacionalmente factibles manteniendo un número de Rossby ($Ro \sim 0.1-0.3$) comparable al de KIC 9244992.
• Rotación: La estrella fue inicializada con rotación de cuerpo rígido ($\Omega \approx 1.0 \times 10^{-5}$ rad s$^{-1}$), aproximadamente 10 veces más rápida que KIC 9244992, para compensar la convección potenciada.
• Campo magnético: Se impuso un campo semilla dipolar ($100$ kG en el límite inferior) una vez que se alcanzó un estado hidrodinámico estacionario.
• Parámetros: La simulación principal utilizó un número de Reynolds ($Re$) de 240, un número de Reynolds magnético ($Rm$) de 960 y un número de Prandtl magnético ($Pm$) de $\sim 4$.
• Análisis: Los datos fueron analizados desde el día 40 hasta el 100 del tiempo físico, una vez que la energía magnética alcanzó un estado estacionario estadístico. El estudio se centró en la región cercana al núcleo (justo fuera del límite convectivo-radiativo, CRB) donde la frecuencia de Brunt–Väisälä alcanza su máximo.

Resultos clave

1. Geometría y fuerza del campo: Las simulaciones generaron de forma autoconsistente campos magnéticos donde la componente azimutal (toroidal) domina a la componente radial (poloidal) en la región cercana al núcleo. En el pico de la frecuencia de Brunt–Väisälä, la relación entre la componente de campo azimutal y la radial de raíz cuadrática media ($B_\phi/B_r$) alcanzó un máximo de 20.

   • En el CRB: $B_r \approx 33.8$ kG, $B_\phi \approx 120.8$ kG.
   • En el pico del pico de la frecuencia de Brunt–Väisälä: $B_r \approx 0.29$ kG, $B_\phi \approx 7.65$ kG.
   • Estos valores están dentro de un factor de pocos para los valores mínimos de intensidad de campo inferidos para KIC 9244992 por Takata et al. (2025).

2. Rotación diferencial y cizallamiento: Aunque el perfil de rotación promediado esféricamente mostró solo una pequeña diferencia ($\sim 3.7\%$) entre la región cercana al núcleo y la superficie —consistente con las inferencias asterosísmicas—, las simulaciones revelaron un cizallamiento latitudinal significativo. La tasa de rotación varió hasta un $20\%$ a través de las latitudes, particularmente cerca del ecuador. Este cizallamiento localizado, en lugar de un cizallamiento radial global, impulsa el efecto $\Omega$, amplificando el campo toroidal a partir del componente poloidal.

3. Estructura latitudinal: La intensidad del campo toroidal exhibió una estructura de doble pico en latitud, con máximos en latitudes medias y mínimos en los polos. Esta geometría se alinea con la inferencia asterosísmica del signo de los coeficientes de los polinomios de Legendre ($a < 0, b > 0$) reportados por Takata et al. (2025).

4. Robustez: Los autores verificaron que los resultados no son artefactos de la intensidad del campo inicial (una ejecución con $1/5$ del campo inicial alcanzó el mismo estado estacionario) ni de elecciones específicas de parámetros, siempre que el número de Reynolds sea suficientemente alto ($Re > \sim 100$) y el número de Rossby sea apropiado. Las pruebas de sensibilidad mostraron que reducir el $Re$a$\sim 30$ resultó en campos significativamente más débiles ($1-2$ kG), lo que resalta la necesidad de un $Re$ alto para la saturación del dinamo.

Significado y afirmaciones
El artículo afirma que la presencia de un campo magnético toroidal dominante en la región cercana al núcleo es una característica genérica de las estrellas con núcleo convectivo, surgiendo naturalmente de las simulaciones de dinamo global cuando se utilizan los números adimensionales apropiados ($Re$, $Ro$).

Los autores enfatizan que sus resultados proporcionan un soporte numérico independiente para la interpretación asterosísmica de KIC 9244992. Crucialmente, el estudio demuestra que campos toroidales fuertes pueden coexistir con una rotación diferencial esféricamente promediada débil, debido a que la generación del campo es impulsada por un cizallamiento latitudinal localizado. Este hallazgo sugiere que las futuras inferencias asterosísmicas no deberían asumir una geometría de campo predominantemente radial en toda la estrella, ya que la región cercana al núcleo probablemente esté dominada por campos toroidales generados por flujos de cizallamiento tridimensionales. El trabajo cierra la brecha entre las restricciones observacionales y la teoría de dinamo numérica, confirmando que la geometría de campo específica inferida para KIC 9244992 es físicamente plausible dentro del marco de la magnetohidrodinámica estelar.