-
Turbulence Regulates the Rate of Planetesimal Formation via Gravitational Collapse
Authors:
Daniel A. Gole,
Jacob B. Simon,
Rixin Li,
Andrew N. Youdin,
Philip J. Armitage
Abstract:
We study how the interaction between the streaming instability and intrinsic gas-phase turbulence affects planetesimal formation via gravitational collapse in protoplanetary disks. Turbulence impedes the formation of particle clumps by acting as an effective turbulent diffusivity, but it can also promote planetesimal formation by concentrating solids, for example in zonal flows. We quantify the ef…
▽ More
We study how the interaction between the streaming instability and intrinsic gas-phase turbulence affects planetesimal formation via gravitational collapse in protoplanetary disks. Turbulence impedes the formation of particle clumps by acting as an effective turbulent diffusivity, but it can also promote planetesimal formation by concentrating solids, for example in zonal flows. We quantify the effect of turbulent diffusivity using numerical simulations of the streaming instability in small local domains, forced with velocity perturbations that establish approximately Kolmogorov-like turbulence. We find that planetesimal formation is suppressed by turbulence once velocity fluctuations exceed $\langle δv^2 \rangle \simeq 10^{-3.5} - 10^{-3} c_s^2$. Turbulence whose strength is just below the threshold reduces the rate at which solids are bound into clumps. Our results suggest that the well-established turbulent thickening of the mid-plane solid layer is the primary mechanism by which turbulence influences planetesimal formation and that planetesimal formation requires a mid-plane solid-to-gas ratio $ε\gtrsim 0.5$. We also quantify the initial planetesimal mass function using a new clump-tracking method to determine each planetesimal mass shortly after collapse. For models in which planetesimals form, we show that the mass function is well-described by a broken power law, whose parameters are robust to the inclusion and strength of imposed turbulence. Turbulence in protoplanetary disks is likely to substantially exceed the threshold for planetesimal formation at radii where temperatures $T \gtrsim 10^3 \ {\rm K}$ lead to thermal ionization. Planetesimal formation may therefore be unviable in the inner disk out to 2-3 times the dust sublimation radius.
△ Less
Submitted 21 October, 2020; v1 submitted 27 January, 2020;
originally announced January 2020.
-
The Nature of Turbulence in the Outer Regions of Protoplanetary Disks
Authors:
Daniel A. Gole,
Jacob B. Simon
Abstract:
We carry out a series of local, shearing box simulations of the outer regions of protoplanetary disks, where ambipolar diffusion is important due to low ionization levels, to better characterize the nature of turbulence and angular momentum transport in these disks. These simulations are divided into two groups, one with far ultraviolet (FUV) ionization included, and one without FUV. In both cases…
▽ More
We carry out a series of local, shearing box simulations of the outer regions of protoplanetary disks, where ambipolar diffusion is important due to low ionization levels, to better characterize the nature of turbulence and angular momentum transport in these disks. These simulations are divided into two groups, one with far ultraviolet (FUV) ionization included, and one without FUV. In both cases, we explore a large range in diffusivity values. We find that in the simulations without FUV, the properties of the turbulence are similar to the unstratified simulations of Bai & Stone (2011); for a given diffusivity, the MRI can still be present so long as the magnetic field is sufficiently weak. Furthermore, the dynamics of the mid-plane in these simulations are primarily controlled by the MRI. In the FUV simulations on the other hand, the MRI-active FUV layers transport strong toroidal magnetic flux to the mid-plane, which shuts off the MRI. Instead, angular momentum transport at the mid-plane is dominated by laminar magnetic fields, resulting in lower levels of turbulent Maxwell stress compared to the no-FUV simulations. Finally, we perform a temporal correlation analysis on the FUV simulations, confirming our result that the dynamics in the mid-plane region is strongly controlled by the FUV-ionized layers.
△ Less
Submitted 24 October, 2018;
originally announced October 2018.
-
Turbulence, Transport and Waves in Ohmic Dead Zones
Authors:
Daniel Gole,
Jacob B. Simon,
Stephen H. Lubow,
Philip J. Armitage
Abstract:
We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change…
▽ More
We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change. We find that the residual Reynolds stress near the mid-plane decreases with increasing dead zone thickness, becoming negligible in cases where the active to dead mass ratio is less than a few percent. This implies that purely Ohmic dead zones would be vulnerable to episodic accretion outbursts via the mechanism of Martin & Lubow (2011). We show that even thick dead zones support a large amount of kinetic energy, but this energy is largely in fluid motions that are inefficient at angular momentum transport. Confirming results from Oishi & Mac Low (2009), the perturbed velocity field in the dead zone is dominated by an oscillatory, vertically extended circulation pattern with a low frequency compared to the orbital frequency. This disturbance has the properties predicted for the lowest order r mode in a hydrodynamic disk. We suggest that in a global disk similar excitations would lead to propagating waves, whose properties would vary with the thickness of the dead zone and the nature of the perturbations (isothermal or adiabatic). Flows with similar amplitudes would buckle settled particle layers and could reduce the efficiency of pebble accretion.
△ Less
Submitted 9 June, 2016;
originally announced June 2016.
-
Substructure in bulk velocities of Milky Way disk stars
Authors:
Jeffrey L. Carlin,
James DeLaunay,
Heidi Jo Newberg,
Licai Deng,
Daniel Gole,
Kathleen Grabowski,
Ge Jin,
Chao Liu,
Xiaowei Liu,
A-Li Luo,
Haibo Yuan,
Haotong Zhang,
Gang Zhao,
Yongheng Zhao
Abstract:
We find that Galactic disk stars near the anticenter exhibit velocity asymmetries in both the Galactocentric radial and vertical components across the mid-plane as well as azimuthally. These findings are based on LAMOST spectroscopic velocities for a sample of ~400,000 F-type stars, combined with proper motions from the PPMXL catalog for which we have derived corrections to the zero points based i…
▽ More
We find that Galactic disk stars near the anticenter exhibit velocity asymmetries in both the Galactocentric radial and vertical components across the mid-plane as well as azimuthally. These findings are based on LAMOST spectroscopic velocities for a sample of ~400,000 F-type stars, combined with proper motions from the PPMXL catalog for which we have derived corrections to the zero points based in part on spectroscopically discovered galaxies and QSOs from LAMOST. In the region within 2 kpc outside the Sun's radius and +/-2 kpc from the Galactic midplane, we show that stars above the plane exhibit net outward radial motions with downward vertical velocities, while stars below the plane have roughly the opposite behavior. We discuss this in the context of other recent findings, and conclude that we are likely seeing the signature of vertical disturbances to the disk due to an external perturbation.
△ Less
Submitted 24 September, 2013;
originally announced September 2013.