We study how the interaction between the streaming instability and intrinsic gas-phase turbulence affects planetesimal formation via gravitational collapse. Turbulence impedes the formation of dense 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 fluid turbulence. We find that planetesimal formation is strongly suppressed by turbulence once velocity fluctuations exceed a threshold value of $\delta v^2 \simeq 10^{-3.5} - 10^{-3} c_s^2$. Turbulence whose strength is just below the threshold reduces the rate of solid material being converted into bound clumps. The main effect of turbulence is to thicken the mid-plane solid layer. Our results are thus consistent with planetesimal formation requiring a mid-plane solid-to-gas ratio $\epsilon \gtrsim 0.5$. We describe a method for tracking bound clumps in our simulations, and use this method to construct a mass function of planetesimals that is measured for each clump shortly after its collapse. Adopting this definition of the initial mass, instead of measuring masses at a fixed time, reduces planetesimal masses by a factor of three. For models in which planetesimals form, we show that the initial 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. Fully in situ planetesimal and planet formation may therefore not be viable for the closest-in exoplanets.

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.

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.

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.