Abstract:Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust. Aims. We aim to study the effects of the vertical shear instability (VSI) on the thermal relaxation times, and vice versa. We are particularly interested in the influence of the initial dust grain size on the VSI and whether the emerging turbulence is sustained over long timescales. Results. We find that the emergence of the VSI is strongly dependent on the initial dust grain size. Coagulation also counteracts the emergence of hydrodynamic turbulence in our simulations, as shown by others before. Starting a simulation with larger grains (100 $\mu$m) generally leads to a less turbulent outcome. While the inner disk regions (within $\sim$ 70 au) develop turbulence in all three simulations, we find that the simulations with larger particles do not develop VSI in the outer disk. Conclusions. Our simulations with dynamically calculated thermal accommodation times based on the drifting and settling dust distribution show that the VSI, once developed in a disk, can be sustained over long timescales, even if grain growth is occurring. The VSI corrugates the dust layer and even diffuses the smaller grains into the upper atmosphere, where they can cool the gas. Whether the instability can emerge for a specific stratification depends on the initial dust grain sizes and the initial dust scale height. If the grains are initially $\gtrsim$ 100 $\mu$m and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the mutual influence between vertical shear instability (VSI) and thermal relaxation time, especially the role in the consistency of dust evolution and cooling time in protoplanetary disks. Specifically, the research focuses on the following aspects: 1. **The influence of VSI on thermal relaxation time and vice versa**: The research aims to explore how VSI affects the gas thermal relaxation time and how the thermal relaxation time in turn affects the development of VSI. Of particular interest is the influence of the initial dust particle size on VSI and whether the resulting turbulence can persist on long - time scales. 2. **The influence of dust particle size on VSI**: The research examines the influence of different initial dust particle sizes (1 micron, 10 microns and 100 microns) on the development of VSI. In particular, the researchers want to know whether larger initial dust particles will suppress the occurrence and development of VSI. 3. **The influence of dust distribution on VSI**: The research also explores the influence of the initial dust distribution (including the height of the dust layer and the size distribution of dust particles) on VSI. These factors determine the cooling time, which in turn affects the occurrence and development of VSI. 4. **The performance of VSI in different regions**: The research analyzes the different performances of VSI in the inner region (30 - 55 astronomical units) and the outer region (115 - 150 astronomical units) of the protoplanetary disk. The results show that VSI is more likely to develop in the inner region, while it is more restricted in the outer region. ### Formula Explanation The key formula involved in the paper is the time scale for the rapid cooling required for VSI to occur: \[ t_c \lesssim \frac{H_g}{R} \left| \beta_T \right| \gamma \frac{2}{(\gamma - 1)} \left( \frac{z}{H_g} \right)^{-1} \Omega_K^{-1} \] where: - $ H_g $ is the gas pressure scale height of the disk, - $ R $ is the distance to the central star, - $ \beta_T $ is the index of the radial temperature structure $ T \propto R^{\beta_T} $ - $ \gamma $ is the specific heat ratio, - $ z $ is the distance from the disk mid - plane, - $ \Omega_K $ is the Keplerian angular frequency. In addition, the research also involves the calculation formula of thermal relaxation time: \[ t_{\text{NLTE}}^{\text{thin}} \approx t_{\text{gas}}^{\text{coll}} = \frac{\gamma}{\gamma - 1} \frac{1}{n_s \bar{v}_{\text{gas}} \sigma_s} \] where: - $ n_s=\frac{\rho_d}{\frac{4}{3} \pi \rho_m a_s^3} $ is the dust number density, - $ \sigma_s = \pi a_s^2 $ is the collision cross - section, - $ \rho_m $ is the material density of the dust, - $ \rho_d $ is the dust volume density. Through these formulas, researchers can simulate and evaluate the development of VSI under different conditions and its impact on protoplanetary disks.

Vertical shear instability with dust evolution and consistent cooling times. On the importance of the initial dust distribution

Vertical shear instability with dust evolution and consistent cooling times

Dust Coagulation Reconciles Protoplanetary Disk Observations with the Vertical Shear Instability. I. Dust Coagulation and the VSI Dead Zone

The Sandwich Mode for Vertical Shear Instability in Protoplanetary Disks

Effects of Dust Evolution on the Vertical Shear Instability in the Outer Regions of Protoplanetary Disks

A self-consistent model for dust settling and the vertical shear instability in protoplanetary disks

Two saturated states of the vertical shear instability in protoplanetary disks with vertically varying cooling times

Dust dynamics and vertical settling in gravitoturbulent protoplanetary discs

Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating

Thermal Structure Determines Kinematics: Vertical Shear Instability in Stellar Irradiated Protoplanetary Disks

Dynamics of Dusty Vortices II: Stability of 2D dust laden vortices

The vertical shear instability in poorly ionised, magnetized protoplanetary discs

Instabilities in dusty non-isothermal proto-planetary discs

Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks II. Secondary instabilities and stability regions

Dust-vortex instability in the regime of well-coupled grains

Dust Settling Instability in Protoplanetary Discs

THREE-DIMENSIONAL SIMULATIONS OF KELVIN-HELMHOLTZ INSTABILITY IN SETTLED DUST LAYERS IN PROTOPLANETARY DISKS

Vertical Shearing Instabilities in Radially Shearing Disks: The Dustiest Layers of the Protoplanetary Nebula

Particle Dynamics in 3D Self-gravitating Disks II: Strong Gas Accretion and Thin Dust Disks

Dust settling against hydrodynamic turbulence in protoplanetary discs