skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: FORMING PLANETESIMALS BY GRAVITATIONAL INSTABILITY. I. THE ROLE OF THE RICHARDSON NUMBER IN TRIGGERING THE KELVIN-HELMHOLTZ INSTABILITY

Journal Article · · Astrophysical Journal
;  [1];  [2];  [3]
  1. Department of Astronomy, University of California, Berkeley, CA 94720 (United States)
  2. Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, NM 87545 (United States)
  3. Department of Physics and Astronomy, San Francisco State University, San Francisco, CA 94132 (United States)

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 three-dimensional simulations of stratified subdisks in the limit that dust particles are small and aerodynamically well coupled to gas, thereby screening out the streaming instability and isolating 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}{sub 0} and find that the critical Richardson number dividing KH-unstable from KH-stable flows is not unique; rather, Ri{sub crit} grows nearly linearly with {mu}{sub 0} for {mu}{sub 0} = 0.3-10. Plausibly, a linear dependence arises for {mu}{sub 0} << 1 because in this regime the radial Kepler shear replaces vertical buoyancy as the dominant stabilizing influence. Why this dependence should persist at {mu}{sub 0} > 1 is a new puzzle. The bulk (height-integrated) metallicity is uniquely determined by Ri and {mu}{sub 0}. Only for disks of bulk solar metallicity is Ri{sub crit} {approx} 0.2, which is close to the classical value. Our empirical stability boundary is such that a dusty sublayer can gravitationally fragment and presumably spawn planetesimals if embedded within a solar metallicity gas disk {approx}4x more massive than the minimum-mass solar nebula; or a minimum-mass disk having {approx}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.

OSTI ID:
21455066
Journal Information:
Astrophysical Journal, Vol. 718, Issue 2; Other Information: DOI: 10.1088/0004-637X/718/2/1367; ISSN 0004-637X
Country of Publication:
United States
Language:
English