We study the turbulence induced in the dust layer of a protoplanetary disk based on the energetics of dust accretion due to gas drag. We estimate turbulence strength from the energy supplied by dust accretion, using the radial drift velocity of the dust particles in a laminar disk. Our estimate of the turbulence strength agrees with previous analytical and numerical research on the turbulence induced by Kelvin-Helmholtz and/or streaming instabilities for particles whose stopping time is less than the Keplerian time. For such small particles, the strongest turbulence is expected to occur when the dust-to-gas ratio of the disk is ~C_eff^(1/2) (h_g / r) ~ 10^(-2), where C_eff ~ 0.2 represents the energy supply efficiency to turbulence and h_g / r ~ 5 x 10^(-2) is the aspect ratio of the gas disk. The maximum viscosity parameter is alpha_max ~ C_eff T_s (h_g / r)^2 ~ 10^(-4) T_s, where T_s (<1) is the non-dimensional stopping time of the dust particles. Modification in the dust-to-gas ratio from the standard value, 10^(-2), by any process, results in weaker turbulence and a thinner dust layer, and consequently may accelerate the growth process of the dust particles.

Instabilities of the dust layer in a protoplanetary disk are investigated. It is known that the streaming instability develops and dust density concentration occurs in a situation where the initial dust density is uniform. This work considers the effect of initial dust density gradient vertical to the midplane. Dust and gas are treated as different fluids. Pressure of dust fluid is assumed to be zero. The gas friction time is assumed to be constant. Axisymmetric two-dimensional numerical simulation was performed using the spectral method. We found that an instability develops with a growth rate on the order of the Keplerian angular velocity even if the gas friction time multiplied by the Keplerian angular velocity is as small as 0.001. This instability is powered by two sources: (1) the vertical shear of the azimuthal velocity, and (2) the relative motion of dust and gas coupled with the dust density fluctuation due to advection. This instability diffuses dust by turbulent advection and the maximum dust density decreases. This means that the dust concentration by the streaming instability which is seen in the case of a uniform initial dust density becomes ineffective as dust density gradient increases by the dust settling toward the midplane.

We show that the planetesimal formation due to the gravitational fragmentation of a dust layer in a protoplanetary disk is possible. The dust density distribution in the dust layer would approach the constant Richardson number distribution due to the dust stirring by the shear instability and dust settling. We perform the analysis of the shear instability of dust layer in a protoplanetary disk with the constant Richardson number density distribution. Our study revealed that this distribution is stable against the shear instability even if the dust density at the midplane reaches the critical density of the gravitational instability, and the planetesimal formation through the gravitational fragmentation of the dust layer can occur even for the dust to gas surface density ratio with the solar composition.

The linear analysis of the instability due to vertical shear in the dust layer of the solar nebula is performed. The following assumptions are adopted throughout this paper: (1) The self-gravity of the dust layer is neglected. (2) One fluid model is adopted, where the dust aggregates have the same velocity with the gas due to strong coupling by the drag force. (3) The gas is incompressible. The calculations with both the Coriolis and the tidal forces show that the tidal force has a stabilizing effect. The tidal force causes the radial shear in the disk. This radial shear changes the wave number of the mode which is at first unstable, and the mode is eventually stabilized. Thus the behavior of the mode is divided into two stages: (1) the first growth of the unstable mode which is similar to the results without the tidal force, and (2) the subsequent stabilization due to an increase of the wave number by the radial shear. If the midplane dust/gas density ratio is smaller than 2, the stabilization occurs before the unstable mode grows largely. On the other hand, the mode grows faster by one hundred orders of magnitude, if this ratio is larger than 20. Because the critical density of the gravitational instability is a few hundreds times as large as the gas density, the hydrodynamic instability investigated in this paper grows largely before the onset of the gravitational instability. It is expected that the hydrodynamic instability develops turbulence in the dust layer and the dust aggregates are stirred up to prevent from settling further. The formation of planetesimals through the gravitational instabilities is difficult to occur as long as the dust/gas surface density ratio is equal to that for the solar abundance.