Abstract:Turbulence in protoplanetary disks plays an important role in dust evolution and planetesimal formation. The vertical shear instability (VSI) is one of the candidate hydrodynamic mechanisms that can generate turbulence in the outer disk regions. The VSI requires rapid gas cooling in addition to vertical shear. A linear stability analysis suggests that the VSI may not operate around the midplane where gas cooling is inefficient. In this study, we investigate the nonlinear outcome of the VSI in disks with a linearly VSI-stable midplane region. We perform two-dimensional global hydrodynamical simulations of an axisymmetric disk with vertically varying cooling times. The vertical cooling time profile determines the thicknesses of the linearly VSI-stable midplane layer and unstable layers above and below the midplane. We find that the thickness of the midplane stable layer determines the vertical structure of VSI-driven turbulence in the nonlinear saturated state. We identify two types of final saturated state: (1) T states characterized by vertical turbulent motion penetrating into the VSI-stable midplane layer and (2) pT states characterized by turbulent motion confined in the unstable layers. The pT states are realized when the midplane VSI-stable layer is thicker than two gas scale heights. We also find that the VSI-driven turbulence is largely suppressed at all heights when the VSI-unstable region lying above and below the midplane is thinner than two gas scale heights. We present empirical formulas that predict the strength of VSI-driven turbulence as a function of the thicknesses of the unstable and stable layers. These formulas will be useful for studying how VSI-driven turbulence and dust grains controlling the disk cooling efficiency evolve simultaneously.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how the turbulence driven by vertical shear instability (VSI) in a protoplanetary disk reaches a saturated state in a disk with vertically varying cooling times. Specifically, through two - dimensional global hydrodynamic simulations, the author studies the vertical structure of VSI - driven turbulence in a non - linearly saturated state when the central region of the disk is linearly stable while the upper and lower layers are unstable. The paper explores two final saturated states: one is that the turbulent motion penetrates into the VSI - stable central region (T - state), and the other is that the turbulent motion is restricted to the unstable upper and lower layers (pT - state). The author also finds that when the thickness of the unstable upper and lower layers is less than two gas scale heights, the VSI - driven turbulence will be significantly suppressed at all heights. In addition, the author provides empirical formulas for predicting the intensity of VSI - driven turbulence. These formulas, as functions of the thicknesses of the unstable and stable layers, are very useful for studying how VSI - driven turbulence and dust particles that control the disk cooling efficiency evolve simultaneously. ### Key questions: 1. **Performance of VSI in the central region of the disk**: Linear stability analysis shows that VSI may not operate in the central region of the disk where the gas cooling efficiency is low. However, when VSI operates above and below the central region, can the central region become turbulent? 2. **Final saturated state of VSI - driven turbulence**: When the central region of the disk is linearly stable, what kind of vertical structure will the VSI - driven turbulence form in a non - linearly saturated state? 3. **Influence of the thicknesses of the unstable and stable layers on turbulence**: How do the thicknesses of the unstable and stable layers affect the intensity of VSI - driven turbulence? ### Main findings: - **Two saturated states**: According to the thickness of the stable layer in the center of the disk, VSI - driven turbulence can be divided into two saturated states: - **T - state**: Vertical turbulent motion penetrates into the VSI - stable central region. - **T - state**: Vertical turbulent motion is restricted to the unstable upper and lower layers. - **Influence of the thickness of the unstable layer**: When the thickness of the unstable upper and lower layers is greater than two gas scale heights, turbulence can develop significantly; when the thickness is less than two gas scale heights, turbulence is suppressed. - **Empirical formulas**: The author provides empirical formulas for predicting the intensity of VSI - driven turbulence. These formulas can be used to study the changes of turbulence during the long - term evolution of the disk. ### Significance: These findings are helpful for understanding the turbulence mechanism in protoplanetary disks, especially the performance of VSI under different cooling conditions. This is of great significance for studying the evolution of dust particles and the formation process of kilometer - sized planetesimals.

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

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

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

The Sandwich Mode for Vertical Shear Instability in Protoplanetary Disks

Vertical shear instability with dust evolution and consistent cooling times

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

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

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

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

Vertical Shear Instability with Partially Reflecting Boundary Conditions

Global Simulations of the Vertical Shear Instability with Nonideal Magnetohydrodynamic Effects

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

Global Simulations of the Vertical Shear Instability with Non-ideal Magnetohydrodynamical Effects

A high resolution simulation of protoplanetary disk turbulence driven by the vertical shear instability

The vertical shear instability in poorly ionised, magnetized protoplanetary discs

On the Vertical Shear Instability in Magnetized Protoplanetary Disks

Turbulence in Outer Protoplanetary Disks: MRI or VSI?

Hydrodynamical simulations of protoplanetary disks including irradiation of stellar photons. I. Resolution study for Vertical Shear Instability (VSI)

Thermal instabilities in accretion disks II: Numerical Experiments for the Goldreich-Schubert-Fricke Instability and the Convective Overstability in disks around young stars