The Role of Thermal Instability in Accretion Outbursts in High-Mass Stars

TL;DR

This study investigates the role of thermal instability in episodic accretion outbursts of high-mass young stellar objects, comparing simulation results with observations and highlighting the need for alternative mechanisms for certain burst types.

Contribution

It demonstrates that thermal instability models can replicate long outburst durations in high-mass stars but are insufficient for short or multiple outbursts, suggesting other processes are involved.

Findings

Thermal instability models match long burst durations in HMYSOs.

Models struggle to reproduce short, rapid outbursts and multiple events.

Alternative mechanisms like gravitational instabilities may explain observed phenomena.

Abstract

High-mass young stellar objects (HMYSOs) can exhibit episodic bursts of accretion, accompanied by intense outflows and luminosity variations. Thermal Instability (TI) due to Hydrogen ionisation is among the most promising mechanisms of episodic accretion in low mass ($M_*\lesssim 1M_{\odot}$) protostars. Its role in HMYSOs has not yet been elucidated. Here, we investigate the properties of TI outbursts in young, massive ($M_*\gtrsim 5M_{\odot}$) stars, and compare them to those observed so far. Our simulations show that modelled TI bursts can replicate the durations and peak accretion rates of long (a few years to decades) outbursts observed in HMYSOs with similar mass characteristics. However, they struggle with short-duration (less than a year) bursts with short (a few weeks or months) rise times, suggesting the need for alternative mechanisms. Moreover, while our models match the durations of longer bursts, they fail to reproduce the multiple outbursts seen in some HMYSOs, regardless of model parameters. We also emphasise the significance of not just evaluating model accretion rates and durations, but also performing photometric analysis to thoroughly evaluate the consistency between model predictions and observational data. Our findings suggest that some other plausible mechanisms, such as gravitational instabilities and disc fragmentation can be responsible for generating the observed outburst phenomena in HMYSOs and underscore the need for further investigation into alternative mechanisms driving short outbursts. However, the physics of TI is crucial in sculpting the inner disc physics in the early bright epoch of massive star formation, and comprehensive parameter space exploration and the use of 2D modeling are essential for obtaining a more detailed understanding of the underlying physical processes.