Variability of the inner dead zone edge in 2D radiation hydrodynamic simulations

TL;DR

This study uses 2D radiation hydrodynamic simulations to explore how thermal instability affects the inner regions of protoplanetary disks, impacting their structure, stability, and potential planet formation processes.

Contribution

It introduces detailed 2D models that analyze the emergence and effects of thermal instability and dead zone dynamics in protoplanetary disks around T Tauri stars.

Findings

TI creates pressure maxima within 1 AU of the disk.

TI disrupts and reestablishes migration traps multiple times.

TI leads to conditions favorable for streaming instability.

Abstract

The inner regions of protoplanetary disks are prone to thermal instability (TI), which can significantly impact the thermal and dynamic evolution of planet-forming regions. Observable as episodic accretion outbursts, such periodic disturbances shape the disk's vertical and radial structure. We investigate the stability of the inner disk edge around a Class II T Tauri star and analyse the consequences of TI on the inner disk's evolution in both the vertical and radial dimensions. A particular focus is laid on the emergence and destruction of solid-trapping pressure maxima. Our 2D axisymmetric radiation hydrodynamic models include the transition to the dead zone from a highly turbulent inner disk, heating by both stellar irradiation and viscous dissipation, vertical and radial radiative transport and an adaptive dust-to-gas mass ratio. The simulated time frames include both the TI- and quiescent phases. We track the TI on S-curves of thermal stability. The TI in our models can develop in disks with moderate accretion rates and results from the activation of the magnetorotational instability (MRI) in the dead zone. The TI creates an extensive MRI active region around the midplane and disrupts the stable pebble- and migration trap at the inner edge of the dead zone. Our simulations consistently show the occurrence of TI-reflares, which, together with the initial TI, produce pressure maxima in the inner disk within 1 AU, possibly providing favourable conditions for streaming instability. On a timescale of a few thousand years, TI regularly disrupts the disk's radial and vertical structure within 1 AU. While several pressure maxima are created, stable migration traps are destroyed and reinstated after the TI phase. Our models provide a foundation for more detailed investigations into phenomena such as short-term variability of accretion rates.