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. Methods. We ran three axisymmetric hydrodynamic simulations of a protoplanetary disk including four dust fluids that initially resemble MRN size distributions of different initial grain sizes. From the local dust densities, we calculated the thermal accommodation timescale of dust and gas and used the result as the thermal relaxation time of the gas in our simulation. We included the effect of dust growth by applying the monodisperse dust growth rate and the typical growth limits. 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 μm) generally leads to a less turbulent outcome. While the inner disk regions (within ∼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 ≳100 μm and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.

Azimuthal and Vertical Streaming Instability at High Dust-to-gas Ratios and on the Scales of Planetesimal Formation

Streaming Instability and Turbulence: Conditions for Planetesimal Formation

Streaming instabilities in accreting protoplanetary disks: A parameter study

Dust Settling Instability in Protoplanetary Discs

Thresholds for Particle Clumping by the Streaming Instability

Observation of Aerodynamic Instability in the Flow of a Particle Stream in a Dilute Gas

Effects of Dust Coagulation on Streaming Instability

Probing Conditions for Strong Clumping by the Streaming Instability: Small Dust Grains and Low Dust-to-gas Density Ratio

Nonlinear Evolution of Streaming Instabilities in Accreting Protoplanetary Disks

Dust dynamics and vertical settling in gravitoturbulent protoplanetary discs

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

On Secular Gravitational Instability in Vertically Stratified Disks

Gravoturbulent Planetesimal Formation: The Positive Effect of long-lived Zonal Flows

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

On the Formation of Planetesimals via Secular Gravitational Instabilities with Turbulent Stirring

Thousands of planetesimals: Simulating the streaming Instability in very large computational domains

Dynamics of Solids in the Midplane of Protoplanetary Disks: Implications for Planetesimal Formation

Planetesimal formation around the snow line: I. Monte Carlo simulations of silicate dust pile-up in a turbulent disk

A Panoptic Model for Planetesimal Formation and Pebble Delivery

Streaming instability in a global patch simulation of protoplanetary disks