The effect of variable stellar magnetic fields on the spin state of T Tauri stars

Resumen Técnico: El Efecto de los Campos Magnéticos Estelares Variables en el Estado de Rotación de las Estrellas T Tauri

Planteamiento del Problema
La evolución de la rotación estelar de las Estrellas T Tauri Clásicas (CTTSs) presenta un desafío significativo para la comprensión de la evolución de los discos protoplanetarios (PPD) y la formación de exoplanetas. Las observaciones revelan periodos de rotación que oscilan entre $\sim$1 y 1 los días, los cuales permanecen relativamente constantes durante la vida esperada del PPD ($\lesssim$10 Myr), lo que contradice la expectativa de que la contracción estelar debería conducir a un aumento significativo de la velocidad de rotación (spin-up). Aunque se presume que las interacciones magnéticas entre la estrella y el disco eliminan el momento angular (AM), observaciones recientes indican que un número considerable de CTTSs existen en un estado de aceleración con parámetros de rapidez ($\omega_s = (r_t/R_{co})^{3/2}$) significativamente inferiores a los predichos por el "estado de torque cero" (ZTS) o los valores de equilibrio de rotación. Específicamente, muchas estrellas exhiben $\omega_s < 0.3$, lo que desafía la suposición de que las estrellas evolucionan hacia un equilibrio de torque. Este artículo investiga si las variaciones temporales en la fuerza del campo magnético dipolar estelar ($B_\star$) pueden explicar estas desviaciones del equilibrio.

Metodología
Los autores utilizaron el código hidrodinámico implícito TAPIR, que modela la estructura y evolución de los PPD combinando ecuaciones hidrodinámicas con viscosidad turbulenta (Shakura & Sunyaev 1973) y un modelo de viscosidad por capas (zona muerta frente a capas turbulentas). El modelo incorpora de forma autoconsistente:

• Interacción Magnética Estrella-Disco: El radio de truncamiento ($r_t$) se calcula equilibrando la presión magnética contra las presiones de ram y de gas.
• Evolución de la Rotación Estelar: El modelo rastrea la derivada temporal de la velocidad angular estelar, contabilizando los torques externos (acreción, vientos estelares impulsados por acreción [APSW] y eyecciones magnetosféricas [MEs]) y la evolución interna (contracción a lo largo de la trayectoria de Hayashi).
• Fotoevaporación: Un modelo de fotoevaporación por rayos X controla las tasas de pérdida de masa.

Se simularon dos escenarios de $B_\star$ variable para masas estelares de 0.3, 0.7 y 1.0 $M_\odot$, inicializados en el ZTS ($\omega_s \approx 0.72$):

1. Variación a Largo Plazo ("Modelo Largo"): Simulando la disminución de la fuerza del campo dipolar ($B_\star$ de 1 kG a 0.1 kG) sobre un periodo de $\tau_B \approx 0.5$ Myr, correspondiente al desarrollo de un núcleo radiativo.
2. Variación a Corto Plazo ("Modelo Corto"): Simulando ciclos magnéticos estelares donde $B_\star$ oscila sinusoidalmente sobre $\tau_{cyc} = 10$ años, imitando posibles reversiones de polaridad o ciclos de actividad.

Contribuciones Clave y Resultados
El estudio demuestra que las variaciones temporales en $B_\star$ alteran significativamente el estado de rotación estelar, siendo la magnitud del efecto dependiente de la relación entre la escala de tiempo de la variación magnética, la escala de tiempo de aceleración ($\tau_{spin}$) y la escala de tiempo viscosa del disco ($\tau_\nu$).

• Desarrollo del Núcleo Radiativo a Largo Plazo: Cuando $B_\star$ disminuye en una escala de tiempo más corta que $\tau_{spin}$ ($\tau_B < \tau_{spin}$), el radio de truncamiento ($r_t$) se ajusta rápidamente (dinámicamente), mientras que el periodo de rotación estelar ($P_\star$) y el radio de corrotación ($R_{co}$) evolucionan lentamente. Este desacoplamiento causa una disminución rápida de $\omega_s$ (por ejemplo, de 0.72 a $\sim$0.17). A medida que la tasa de acreción disminuye posteriormente debido a la evolución viscosa y la fotoevaporación, $r_t$ aumenta, devolviendo eventualmente a $\omega_s$ hacia la unidad. Los autores encuentran que para estrellas con masas $M_\star \gtrsim 0.7 M_\odot$, el desarrollo de un núcleo radiativo puede desplazar el estado de rotación fuera del ZTS, explicando potencialmente los bajos valores observados de $\omega_s$ en CTTSs más antiguas.
• Ciclos Magnéticos a Corto Plazo: Para ciclos con $\tau_{cyc} \ll \tau_{\nu,in}$ (escala de tiempo viscosa del disco interno $\sim$1000 años), el material del disco no puede ajustarse rápidamente al cambio en $r_t$. Esto resulta en una "reacción inversa" donde la tasa de acreción varía inversamente a los cambios en el campo magnético (por ejemplo, un $B_\star$ más alto empuja $r_t$ hacia afuera, pero el disco no puede evacuarse lo suficientemente rápido, lo que conduce a una mayor acreción). Esta reacción inversa debilita el efecto neto sobre $\omega_s$, causando que este oscile (por ejemplo, entre 0.50 y 0.87) pero evitando una desviación grande y sostenida del ZTS. Sin embargo, si $\tau_{cyc} \gtrsim \tau_{\nu,in}$, el disco puede ajustarse y el efecto sobre $\omega_s$ se vuelve más pronunciado.
• Comparación con Observaciones: Los autores compararon los resultados de su "Modelo Largo" con los datos observacionales de Pittman et al. (2025). Categorizaron las estrellas como "jóvenes" (edad $< \tau_{conv} + \tau_B$) y "viejas" (edad $\ge \tau_{conv} + \tau_B$). Las estrellas "viejas" exhibieron un $\omega_s$ medio más bajo ($0.24 \pm 0.16$) en comparación con las estrellas "jóvenes" ($0.55 \pm 0.40$), lo cual es consistente con la predicción del modelo de que el desarrollo de un núcleo radiativo impulsa a $\omega_s$ hacia valores más bajos.

Significancia y Reivindicaciones
El artículo sostiene que las variaciones en los campos magnéticos estelares proporcionan un mecanismo plausible para explicar al menos una parte de las estrellas T Tauri observadas que existen fuera del equilibrio de rotación predicho o del ZTS. Específicamente, la transición de un estado completamente convectivo a uno parcialmente convectivo (formación de un núcleo radiativo) puede impulsar un cambio significativo y temporal en el estado de rotación hacia parámetros de rapidez más bajos.

Los autores mantienen la modestia respecto a la completitud de su explicación. Señalan que:

• La dispersión observada en $\omega_s$ también puede estar influenciada por incertidumbres en las estimaciones de la edad y la masa estelar, así como por la variabilidad a corto plazo en las tasas de acreción.
• Los modelos de rotación actuales pueden ser incompletos, pudiendo omitir otros mecanismos de frenado como los vientos de disco en las regiones más internas.
• El escenario de los ciclos magnéticos a corto plazo requiere mayor confirmación observacional de ciclos regulares en estrellas T Tauri y mediciones simultáneas de $B_\star$, $r_t$ y $P_\star$.

El estudio concluye que conectar la acreción, la rotación estelar y las propiedades magnéticas requiere de más trabajo teórico y observacional, particularmente involucrando misiones futuras como el Early eVolution Explorer (EVE) y sondeos como SPIRou para constreñir los ciclos magnéticos y su impacto en la evolución de la rotación.