A Parameterized YSO Accretion Disk Model with Increasing Accretion Rate: Predicted Outburst Lightcurves

Resumen Técnico: Un Modelo de Disco de Acreción de Objetos Estelares Jóvenes Parametrizado con Tasa de Acreción Creciente

Enunciado del Problema
Los Objetos Estelares Jóvenes (YSO), particularmente la subclase FU Orionis (FU Ori), experimentan erupciones de acreción dramáticas caracterizadas por aumentos de luminosidad de varios órdenes de magnitud a lo largo de meses a años, seguidos de decaimientos durante décadas. Si bien la transición de la acreción magnetosférica en estado bajo a la acreción en erupción de estado alto es un fenómeno reconocido, comprender la evolución específica de los componentes del sistema —fotosfera estelar, choques magnetosféricos, discos de gas viscosos y discos de polvo pasivos— durante estos eventos sigue siendo un desafío. Los modelos existentes a menudo dependen de simulaciones numéricas complejas y costosas computacionalmente (por ejemplo, MHD) o de suposiciones de estado estacionario que podrían no capturar los cambios morfológicos dependientes del tiempo observados en las curvas de luz multibanda. Existe la necesidad de un marco parametrizado que acople estos componentes físicos para simular la evolución temporal de las distribuciones de energía espectral (SED) y las curvas de luz en respuesta a tasas de acreción variables.

Metodología
Los autores presentan un pipeline de modelado parametrizado construido sobre el código YSOpy, diseñado para simular fotometría de series temporales y espectros para estrellas jóvenes en acreción. A diferencia de la transferencia radiativa completa o las simulaciones numéricas autoconsistentes, este enfoque utiliza una infraestructura modular para calcular espectros para cuatro componentes distintos del sistema por separado antes de sumarlos:

1. Fotosfera Estelar: Modelada como un cuerpo negro o una capa atmosférica.
2. Choques de Acreción Magnetosféricos: Modelados utilizando una aproximación de capa de hidrógeno donde la temperatura ($T_{hotspot}$), la densidad electrónica ($n_e$) y la profundidad óptica ($\tau$) se derivan de manera autoconsistente a partir de la tasa de acreción ($\dot{M}$) y la velocidad de caída libre ($v_{ff}$).
3. Disco de Gas Calentado Viscamente: Modelado utilizando un perfil de temperatura modificado de Shakura & Sunyaev (1973) para el disco interno, utilizando rejillas atmosféricas BT-Settl para temperaturas $>1400$ K y cuerpos negros para regiones más frías.
4. Disco de Polvo Calentado Pasivamente: Modelado utilizando un perfil de temperatura de Chiang & Goldreich (1997) para polvo irradiado, con un límite interior determinado por el radio de sublimación del polvo ($R_{sub}$).

El pipeline itera a través de un perfil de tasa de acreción dependiente del tiempo ($\dot{M}(t)$). En cada paso de tiempo, el modelo calcula dinámicamente radios críticos, incluido el radio de truncamiento magnetosférico ($R_{trunc}$), el radio de corrotación ($R_{co-rot}$) y el radio del disco de polvo interior ($R_{dust,inner}$). El radio del disco de gas interior ($R_{in}$) se determina por el equilibrio de las presiones magnética y de acreción, transitando desde $R_{co-rot}$ a bajas $\dot{M}$ hasta $R_{trunc}$ (o $R_*$) a altas $\dot{M}$. El modelo considera la "desactivación" de la acreción magnetosférica poloidal cuando el disco alcanza la superficie estelar ($R_{in} \approx R_*$). Finalmente, los espectros de alta resolución con resolución temporal se convolucionan con funciones de transmisión de filtros estándar (por ejemplo, Gaia, 2MASS, WISE) para generar curvas de luz y curvas de color multibanda.

Contribuciones Clave

• Marco Parametrizado Integrado: El artículo introduce un pipeline optimizado que acopla componentes estelares, de choque, de gas y de polvo en un único modelo evolutivo temporal, permitiendo la simulación de ciclos completos de erupción sin el costo computacional de las simulaciones MHD 3D.
• Física de Choques Autoconsistente: El modelo deriva temperaturas de puntos calientes y densidades electrónicas directamente de la tasa de acreción y la velocidad de caída libre, asegurando la conservación de la energía ($L_{slab} = L_{blackbody}$) mediante un factor de escala $\chi$.
• Límites de Disco Dinámicos: El modelo rastrea explícitamente la evolución del límite del disco de gas interior y el frente de sublimación del polvo como funciones de $\dot{M}$, capturando la transición de un disco limitado por corrotación a un disco truncado magnéticamente, y finalmente a un disco que toca la superficie estelar.
• Herramientas de Diagnóstico Multibanda: La infraestructura genera curvas de luz y curvas de color sintéticas en bandas ópticas e infrarrojas, proporcionando una herramienta de diagnóstico para interpretar las morfologías de erupción observadas.

Resultados
Los autores aplicaron el modelo a una tasa de acreción que aumenta linealmente y a cuatro perfiles de acreción específicos que imitan objetos FU Ori observados: HBC 722 (lineal de tres partes), V890 Aur (tangente hiperbólica), V960 Mon (exponencial bilateral) y Gaia 17bpi (logístico $\times$ Gaussiano). Los hallazgos clave incluyen:

• Evolución de la Dominancia de Componentes:

  • Baja $\dot{M}$: El sistema está dominado por la fotosfera estelar y el disco de polvo pasivo; los choques magnetosféricos son insignificantes o más fríos que la fotosfera.
  • $\dot{M}$ Creciente: A medida que aumenta $\dot{M}$, los puntos calientes magnetosféricos se iluminan rápidamente, dominando primero las bandas UV y ópticas azules (por ejemplo, banda u).
  • Pico/Alta $\dot{M}$: El calentamiento viscoso del disco de gas interior se convierte en la fuente de flujo dominante, particularmente en el infrarrojo cercano (bandas Gaia, 2MASS). Si $\dot{M}$ supera $\sim 10^{-5} - 10^{-4} M_\odot$ yr$^{-1}$, el radio de truncamiento del disco se mueve hacia adentro hasta la superficie estelar, apagando la acreción poloidal y la emisión de puntos calientes.
  • Comportamiento en Infrarrojo Medio: Las bandas de infrarrojo medio (WISE) responden principalmente al calentamiento y expansión del frente de sublimación del polvo, a menudo brillando antes que las bandas ópticas a medida que el radio del disco de polvo interior se mueve hacia afuera.

• Morfología de la Curva de Luz:

  • Las curvas de luz ópticas rojas y de infrarrojo cercano generalmente siguen el perfil de la tasa de acreción de entrada.
  • Las curvas de luz de infrarrojo medio son más sensibles a la geometría y el calentamiento del disco de polvo interior, mostrando a veces un brillo "pre-cursor" o una temporización distinta en relación con las bandas ópticas.
  • La evolución del color varía: los sistemas a menudo se vuelven "más azules cuando son más brillantes" durante el ascenso debido al calentamiento de los puntos calientes y del disco interior, pero pueden enrojecerse durante el decaimiento a medida que el disco viscoso se enfría y la fotosfera vuelve a emerger.

• Sensibilidad a Parámetros:

  • Intensidad del Campo Magnético: Campos más fuertes retrasan el acercamiento del disco a la superficie estelar, prolongando la fase de acreción magnetosférica y creando curvas de luz con múltiples picos en el UV/óptico.
  • Periodo de Rotación: Las variaciones en el periodo de corrotación afectan principalmente al radio del disco interior en estado bajo, pero tienen un impacto negligible en las curvas de luz durante la fase de erupción de alta acreción.

Significado y Afirmaciones
Los autores posicionan este trabajo como una herramienta ilustrativa para interpretar los mecanismos físicos detrás de las erupciones de FU Ori, en lugar de un modelo predictivo para objetos individuales específicos. Afirman que su enfoque parametrizado demuestra exitosamente cómo diferentes componentes del sistema dominan el flujo en diversas longitudes de onda a lo largo de un ciclo de erupción.

El artículo destaca que la temporización del brillo a través de las longitudes de onda ofrece poder de diagnóstico: el brillo en UV/óptico azul indica el inicio del calentamiento por choque magnetosférico, mientras que el brillo en infrarrojo medio señala la expansión del frente de sublimación del polvo. Los autores sugieren que futuras instalaciones como LSST y UVEX podrían utilizar estas firmas predichas para identificar erupciones en etapas tempranas. Además, el modelo proporciona un marco para inferir perfiles de acreción subyacentes y tasas de acreción en estado bajo a partir de curvas de luz y curvas de color observadas, cerrando la brecha entre descripciones fenomenológicas simples y simulaciones físicas complejas. Los autores reconocen que, aunque su modelo simplifica la reestructuración del campo magnético y las capas límite de acreción ecuatorial, captura la física esencial de la transición de estados de acreción similares a T-Tauri a estados similares a FU Ori.