Abstract:Protoplanetary disks exhibit a vertical gradient in angular momentum, rendering them susceptible to the Vertical Shear Instability (VSI). The most important condition for the onset of this mechanism is a short timescale of thermal relaxation ($\lesssim 0.1$ orbital timescales). Simulations of fully VSI active disks are characterized by turbulent, vertically extended dust layers. This is in contradiction with recent observations of the outer regions of some protoplanetary disks, which appear highly settled. In this work, we demonstrate that the process of dust coagulation can diminish the cooling rate of the gas in the outer disk and extinct the VSI activity. Our findings indicate that the turbulence strength is especially susceptible to variations in the fragmentation velocity of the grains. A small fragmentation velocity of $\approx 100 \mathrm{\, cm \,s^{-1}}$ results in a fully turbulent simulation, whereas a value of $\approx 400 \mathrm{\, cm \,s^{-1}}$ results in a laminar outer disk, being consistent with observations. We show that VSI turbulence remains relatively unaffected by variations in the maximum particle size in the inner disk regions. However, we find that dust coagulation can significantly suppress the occurrence of VSI turbulence at larger distances from the central star.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: In the protoplanetary disk, how does the relationship between vertical shear instability (VSI) and dust coagulation affect the observed disk structure? Specifically, the author explored how the dust coagulation process suppresses VSI activity by slowing down the gas cooling rate in the outer disk region, thereby explaining why the outer regions of some protoplanetary disks can remain in a highly stable state, which is inconsistent with the turbulent state predicted in previous simulations.
### Background and Problem Description
The gas in the protoplanetary disk exhibits an angular momentum gradient in the vertical direction, making it susceptible to vertical shear instability (VSI). The occurrence of VSI requires a short thermal relaxation time (usually less than 0.1 orbital periods). However, recent observations show that the outer regions of some protoplanetary disks present highly stable dust layers, which is in contradiction to the fully active VSI simulation results. These simulations usually show turbulent and vertically extended dust layers.
### Main Contributions of the Paper
By studying the influence of the dust coagulation process on the gas cooling rate, the author found that the coagulation of dust particles can significantly reduce the cooling rate in the outer disk region, thereby suppressing VSI activity. This finding indicates that the coagulation characteristics of dust particles (such as the fragmentation velocity) have an important influence on the intensity of VSI. For example, when the fragmentation velocity of dust particles is approximately 100 cm/s, the simulation results show strong turbulence in the disk; when the fragmentation velocity is increased to approximately 400 cm/s, the outer region of the disk shows a laminar state, which is consistent with the observation results.
### Key Conclusions
1. **Dust Coagulation Suppresses VSI**: The coagulation of dust particles can significantly reduce the cooling rate in the outer disk region, thereby suppressing the activity of VSI.
2. **Influence of Fragmentation Velocity**: The fragmentation velocity of dust particles has an important influence on the intensity of VSI. A lower fragmentation velocity leads to stronger turbulence, while a higher fragmentation velocity leads to a more stable disk structure.
3. **Less Influence in the Inner Disk Region**: The turbulent intensity of VSI in the inner disk region is less affected by the change in the maximum particle size, but in the outer disk region, dust coagulation can significantly suppress the occurrence of VSI.
### Formulas and Models
Several formulas are used in the paper to describe the relationship between dust coagulation and VSI:
- **Brunt - Väisälä Frequency**:
\[
N_z^2 = -\frac{1}{\rho_g c_P} \frac{\partial P}{\partial z} \frac{\partial S}{\partial z}
\]
where \( \rho_g \) is the gas density, \( P \) is the pressure, \( S \) is the gas entropy, and \( c_P \) is the specific heat capacity at constant pressure.
- **VSI Growth Conditions**:
\[
t_c < \frac{H_g}{R} \left| \frac{\beta_T}{\gamma - 1} \right| \Omega_K^{-1}
\]
where \( R \) is the distance to the central star, \( \beta_T \) is the power - law index of the temperature distribution, \( H_g \) is the pressure scale height, \( \Omega_K \) is the local Keplerian frequency, and \( \gamma=\frac{C_P}{C_V} \) is the ratio of specific heat capacities of the gas.
- **Thermal Relaxation Time**:
\[
t_{\text{NLTE}}^{\text{thin}} = 2t_{\parallel} \left[ 1 - \sqrt{1 - \frac{4t_{\parallel}^2}{t_{\text{coll}, g} t_{\text{rad}, d}}} \right]^{-1}
\]
where \( t_{\parallel} = \left( \frac{1}{t_{\text{rad}, d}} + \frac{1}{t_{\text{coll}, d}} + \frac{1}{t_{\text{coll}, g}} \right)^{-1} \)