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

Technical Summary: A Parameterized YSO Accretion Disk Model with Increasing Accretion Rate

Problem Statement
Young Stellar Objects (YSOs), particularly the FU Orionis (FU Ori) subclass, undergo dramatic accretion outbursts characterized by luminosity increases of several orders of magnitude over months to years, followed by decays over decades. While the transition from low-state magnetospheric accretion to high-state outburst accretion is a recognized phenomenon, understanding the specific evolution of system components—stellar photosphere, magnetospheric shocks, viscous gas disks, and passive dust disks—during these events remains challenging. Existing models often rely on complex, computationally expensive numerical simulations (e.g., MHD) or steady-state assumptions that may not capture the time-dependent morphological changes observed in multi-band lightcurves. There is a need for a parameterized framework that couples these physical components to simulate the temporal evolution of spectral energy distributions (SEDs) and lightcurves in response to varying accretion rates.

Methodology
The authors present a parameterized modeling pipeline built upon the YSOpy code, designed to simulate time-series photometry and spectra for accreting young stars. Unlike full radiative transfer or self-consistent numerical simulations, this approach uses a modular infrastructure to calculate spectra for four distinct system components separately before summing them:

1. Stellar Photosphere: Modeled as a blackbody or atmospheric slab.
2. Magnetospheric Accretion Shocks: Modeled using a hydrogen slab approximation where the temperature ($T_{hotspot}$), electron density ($n_e$), and optical depth ($\tau$) are derived self-consistently from the accretion rate ($\dot{M}$) and free-fall velocity ($v_{ff}$).
3. Viscously Heated Gas Disk: Modeled using a modified Shakura & Sunyaev (1973) temperature profile for the inner disk, utilizing BT-Settl atmospheric grids for temperatures $>1400$ K and blackbodies for cooler regions.
4. Passively Heated Dust Disk: Modeled using a Chiang & Goldreich (1997) temperature profile for irradiated dust, with an inner boundary determined by the dust sublimation radius ($R_{sub}$).

The pipeline iterates through a time-dependent accretion rate profile ($\dot{M}(t)$). At each time step, the model dynamically calculates critical radii, including the magnetospheric truncation radius ($R_{trunc}$), the co-rotation radius ($R_{co-rot}$), and the inner dust disk radius ($R_{dust,inner}$). The inner gas disk radius ($R_{in}$) is determined by the balance of magnetic and accretion pressures, transitioning from $R_{co-rot}$ at low $\dot{M}$ to $R_{trunc}$ (or $R_*$) at high $\dot{M}$. The model accounts for the "turning off" of poloidal magnetospheric accretion when the disk reaches the stellar surface ($R_{in} \approx R_*$). Finally, the time-resolved high-resolution spectra are convolved with standard filter transmission functions (e.g., Gaia, 2MASS, WISE) to generate multi-band lightcurves and color curves.

Key Contributions

• Integrated Parameterized Framework: The paper introduces a streamlined pipeline that couples stellar, shock, gas, and dust components into a single time-evolving model, allowing for the simulation of full outburst cycles without the computational cost of 3D MHD simulations.
• Self-Consistent Shock Physics: The model derives hotspot temperatures and electron densities directly from the accretion rate and free-fall velocity, ensuring energy conservation ($L_{slab} = L_{blackbody}$) via a scaling factor $\chi$.
• Dynamic Disk Boundaries: The model explicitly tracks the evolution of the inner gas disk boundary and the dust sublimation front as functions of $\dot{M}$, capturing the transition from a co-rotation-limited disk to a magnetically truncated disk, and finally to a disk touching the stellar surface.
• Multi-Band Diagnostic Tools: The infrastructure generates synthetic lightcurves and color curves across optical and infrared bands, providing a diagnostic tool to interpret observed outburst morphologies.

Results
The authors applied the model to a linearly increasing accretion rate and four specific accretion profiles mimicking observed FU Ori objects: HBC 722 (three-piece linear), V890 Aur (hyperbolic tangent), V960 Mon (two-sided exponential), and Gaia 17bpi (logistic $\times$ Gaussian). Key findings include:

• Component Dominance Evolution:

  • Low $\dot{M}$: The system is dominated by the stellar photosphere and passive dust disk; magnetospheric shocks are negligible or cooler than the photosphere.
  • Rising $\dot{M}$: As $\dot{M}$ increases, magnetospheric hotspots brighten rapidly, dominating the UV and blue optical bands (e.g., u-band) first.
  • Peak/High $\dot{M}$: Viscous heating of the inner gas disk becomes the dominant flux source, particularly in the near-infrared (Gaia, 2MASS bands). If $\dot{M}$ exceeds $\sim 10^{-5} - 10^{-4} M_\odot$ yr$^{-1}$, the disk truncation radius moves inward to the stellar surface, shutting off poloidal accretion and hotspot emission.
  • Mid-Infrared Behavior: Mid-infrared bands (WISE) respond primarily to the heating and expansion of the dust sublimation front, often brightening before the optical bands as the inner dust disk radius moves outward.

• Lightcurve Morphology:

  • Red optical and near-infrared lightcurves generally track the input accretion rate profile.
  • Mid-infrared lightcurves are more sensitive to the geometry and heating of the inner dust disk, sometimes exhibiting a "pre-cursor" brightening or distinct timing relative to optical bands.
  • Color evolution varies: systems often become "bluer-when-brighter" during the rise due to hotspot and inner disk heating, but may redden during the decay as the viscous disk cools and the photosphere re-emerges.

• Parameter Sensitivity:

  • Magnetic Field Strength: Stronger fields delay the disk's approach to the stellar surface, prolonging the magnetospheric accretion phase and creating multi-peaked lightcurves in the UV/optical.
  • Rotation Period: Variations in the co-rotation period primarily affect the low-state inner disk radius but have negligible impact on the lightcurves during the high-accretion outburst phase.

Significance and Claims
The authors position this work as an illustrative tool for interpreting the physical mechanisms behind FU Ori outbursts rather than a predictive model for specific individual objects. They claim that their parameterized approach successfully demonstrates how different system components dominate the flux at various wavelengths throughout an outburst cycle.

The paper highlights that the timing of brightening across wavelengths offers diagnostic power: UV/blue optical brightening indicates the onset of magnetospheric shock heating, while mid-infrared brightening signals the expansion of the dust sublimation front. The authors suggest that future facilities like LSST and UVEX could utilize these predicted signatures to identify early-stage outbursts. Furthermore, the model provides a framework for inferring underlying accretion profiles and low-state accretion rates from observed lightcurves and color curves, bridging the gap between simple phenomenological descriptions and complex physical simulations. The authors acknowledge that while their model simplifies magnetic field restructuring and equatorial accretion boundary layers, it captures the essential physics of the transition from T-Tauri-like to FU Ori-like accretion states.