Gravitational instability (GI) of a dust-rich layer at the midplane of a gaseous circumstellar disk is one proposed mechanism to form planetesimals, the building blocks of rocky planets and gas giant cores. Self-gravity competes against the Kelvin-Helmholtz instability (KHI): gradients in dust content drive a vertical shear which risks overturning the dusty subdisk and forestalling GI. To understand the conditions under which the disk can resist the KHI, we perform 3D simulations of stratified subdisks in the limit that dust particles are small and aerodynamically well coupled to gas. This limit screens out the streaming instability and isolates the KHI. Each subdisk is assumed to have a vertical density profile given by a spatially constant Richardson number Ri. We vary Ri and the midplane dust-to-gas ratio mu and find that the critical Richardson number dividing KH-unstable from KH-stable flows is not unique; rather Ri_crit grows nearly linearly with mu for mu=0.3-10. Only for disks of bulk solar metallicity is Ri_crit ~ 0.2, close to the classical value. Our results suggest that a dusty sublayer can gravitationally fragment and presumably spawn planetesimals if embedded within a solar metallicity gas disk ~4x more massive than the minimum-mass solar nebula; or a minimum-mass disk having ~3x solar metallicity; or some intermediate combination of these two possibilities. Gravitational instability seems possible without resorting to the streaming instability or to turbulent concentration of particles.

Dust at the midplane of a circumstellar disk can become gravitationally unstable and fragment into planetesimals if the local dust-to-gas density ratio mu is sufficiently high. We simulate how dust settles in passive disks and ask how high mu can become. We settle the dust using a 1D code and test for dynamical stability using a 3D shearing box code. This scheme allows us to explore the behavior of small particles having short but non-zero stopping times in gas: 0 < t_stop << the orbital period. The streaming instability is thereby filtered out. Dust settles until shearing instabilities in the edges of the dust layer threaten to overturn the entire layer. In this state of marginal stability, mu=2.9 for a disk whose bulk (height-integrated) metallicity is solar. For a disk whose bulk metallicity is 4x solar, mu reaches 26.4. These maximum values of mu, which depend on the background radial pressure gradient, are so large that gravitational instability of small particles is viable in disks whose bulk metallicities are just a few (<4) times solar. Earlier studies assumed that dust settles until the Richardson number Ri is spatially constant. Our simulations are free of this assumption but provide support for it within the dust layer's edges, with the proviso that Ri increases with bulk metallicity in the same way that we found in Paper I. Only modest enhancements in bulk metallicity are needed to spawn planetesimals directly from small particles.