For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to represent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy loss during real-world collisions, implement it in an $N$-body code and test it. The tests provide interesting insights into dynamics of self gravitating collisional systems. They show how particle systems evolve over several free fall timescales to form central concentrations and equilibrated outer shells. The superparticle method can be extended to account for the accretional growth of objects during inelastic mergers.

The streaming instability (SI) has been extensively studied in the linear and non-linear regimes as a mechanism to concentrate solids and trigger planetesimal formation in the midplane of protoplanetary discs. A related dust settling instability (DSI) applies to particles while settling towards the midplane. The DSI has previously been studied in the linear regime, with predictions that it could trigger particle clumping away from the midplane. This work presents a range of linear calculations and non-linear simulations, performed with FARGO3D, to assess conditions for DSI growth. We expand on previous linear analyses by including particle size distributions and performing a detailed study of the amount of background turbulence needed to stabilize the DSI. When including binned size distributions, the DSI often produces converged growth rates with fewer bins than the standard SI. With background turbulence, we find that the most favorable conditions for DSI growth are weak turbulence, characterized by $\alpha \lesssim 10^{-6}$ with intermediate-sized grains that settle from one gas scale-height. These conditions could arise during a sudden decrease in disc turbulence following an accretion outburst. Ignoring background turbulence, we performed a parameter survey of local 2D DSI simulations. Particle clumping was either weak or occurred slower than particles settle. Clumping was reduced by a factor of two in a comparison 3D simulation. Overall, our results strongly disfavor the hypothesis that the DSI significantly promotes planetesimal formation. Non-linear simulations of the DSI with different numerical methods could support or challenge these findings.

The streaming instability (SI) is a mechanism to aerodynamically concentrate solids in protoplanetary disks and facilitate the formation of planetesimals. Recent numerical modeling efforts have demonstrated the increasing complexity of the initial mass distribution of planetesimals. To better constrain this distribution, we conduct SI simulations including the self-gravity with hitherto the highest resolution. To subsequently identify all of the self-bound clumps, we develop a new clump-finding tool, PLanetesimal ANalyzer (\texttt{PLAN}). We then apply a maximum likelihood estimator to fit a suite of parameterized models with different levels of complexity to the simulated mass distribution. To determine which models are best-fitting and statistically robust, we apply three model selection criteria with different complexity penalties. We find that the initial mass distribution of clumps is not universal regarding both the functional forms and parameter values. Our model selection criteria prefer models different from those previously considered in the literature. Fits to multi-segment power law models break to a steeper distribution above masses close to 100 km collapsed planetesimals, similar to observed Kuiper Belt size distributions. We find evidence for a turnover in the low mass end of the planetesimal mass distribution in our high resolution run. Such a turnover is expected for gravitational collapse, but had not previously been reported.

The New Horizons flyby of the cold classical Kuiper Belt object MU69 showed it to be a contact binary. The existence of other contact binaries in the 1-10km range raises the question of how common these bodies are and how they evolved into contact. Here we consider that the lobes of MU69 formed as a binary in the Solar nebula, and calculate its orbital evolution in the presence of gas. We find that the sub-Keplerian wind of the disk brings the drag timescales for 10km bodies under 1 Myr for quadratic-velocity drag. In the Kuiper belt, however, the drag is linear with velocity and the effect of the wind cancels out as the angular momentum gained in half an orbit is lost in the other half: the drag timescales for 10 km bodies remain over 10 Myr. In this situation we find that a combination of drag and Kozai-Lidov oscillations is a promising channel for collapse. We analytically solve the hierarchical three-body problem with drag and implement it into a Kozai plus tidal friction model. The permanent quadrupoles of the lobes make the Kozai oscillations stochastic, and as drag shrinks the semimajor axis it more easily allows the fluctuations to bring the system into contact. Evolution to contact happens very rapidly (within $10^4$ yr) in the pure, double-average quadrupole, Kozai region between $\approx 85-95^\circ$, and within 3 Myr in the drag-assisted region beyond it. The synergy between $J_2$ and gas drag widens the window of contact to 80-100$^\circ$ initial inclination, over a larger range of semimajor axes. As such, the model predicts an initial contact binary fraction of about 10\% for the cold classicals in the Kuiper belt. The speed at contact deviates from the escape velocity only because of the oblateness. For MU69, the contact velocity should be 3.3-4.2 m/s, in line with the observational evidence from the lack of deformation features and estimate of the tensile strength.

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.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

The streaming instability concentrates solid particles in protoplanetary disks, leading to gravitational collapse into planetesimals. Despite its key role in producing particle clumping and determining critical length scales in the instability's linear regime, the influence of the disk's radial pressure gradient on planetesimal properties has not been examined in detail. Here, we use streaming instability simulations that include particle self-gravity to study how the planetesimal initial mass function depends on the radial pressure gradient. Fitting our results to a power-law, ${\rm d}N / {\rm d}M_p \propto M_p^{-p}$, we find a single value $p \approx 1.6$ describes simulations in which the pressure gradient varies by $\gtrsim 2$. An exponentially truncated power-law provides a significantly better fit, with a low mass slope of $p^\prime \approx 1.3$ that weakly depends on the pressure gradient. The characteristic truncation mass is found to be $\sim M_G = 4 \pi^5 G^2 \Sigma_p^3 / \Omega^4$. We exclude the cubic dependence of the characteristic mass with pressure gradient suggested by linear considerations, finding instead a linear scaling. These results strengthen the case for a streaming-derived initial mass function that depends at most weakly on the aerodynamic properties of the disk and participating solids. A simulation initialized with zero pressure gradient---which is {\em not} subject to the streaming instability---also yields a top-heavy mass function but with modest evidence for a different shape. We discuss the consistency of the theoretically predicted mass function with observations of Kuiper Belt planetesimals, and describe implications for models of early stage planet formation..

A critical step toward the emergence of planets in a protoplanetary disk consists in accretion of planetesimals, bodies 1-1000 km in size, from smaller disk constituents. This process is poorly understood partly because we lack good observational constraints on the complex physical processes that contribute to planetesimal formation. In the outer solar system, the best place to look for clues is the Kuiper belt, where icy planetesimals survived to this day. Here we report evidence that Kuiper belt planetesimals formed by the streaming instability, a process in which aerodynamically concentrated clumps of pebbles gravitationally collapse into 100-km-class bodies. Gravitational collapse was previously suggested to explain the ubiquity of equal-size binaries in the Kuiper belt. We analyze new hydrodynamical simulations of the streaming instability to determine the model expectations for the spatial orientation of binary orbits. The predicted broad inclination distribution with 80% of prograde binary orbits matches the observations of trans-Neptunian binaries. The formation models which imply predominantly retrograde binary orbits can be ruled out. Given its applicability over a broad range of protoplanetary disk conditions, it is expected that the streaming instability seeded planetesimal formation also elsewhere in the solar system, and beyond.

The Streaming Instability (SI) is a mechanism to concentrate solids in protoplanetary disks. Nonlinear particle clumping from the SI can trigger gravitational collapse into planetesimals. To better understand the numerical robustness of the SI, we perform a suite of vertically-stratified 3D simulations with fixed physical parameters known to produce strong clumping. We vary the numerical implementation, namely the computational domain size and the vertical boundary conditions (vBCs), comparing newly-implemented outflow vBCs to the previously-used periodic and reflecting vBCs. We find strong particle clumping by the SI is mostly independent of the vBCs. However, peak particle densities are higher in larger simulation domains due to a larger particle mass reservoir. We report SI-triggered zonal flows, i.e., azimuthally-banded radial variations of gas pressure. These structures have low amplitudes, insufficient to halt particle radial drift, confirming that particle trapping in gas pressure maxima is not the mechanism of the SI. We find that outflow vBCs produce artificially large gas outflow rates at vertical boundaries. However, the outflow vBCs reduce artificial reflections at vertical boundaries, allowing more particle sedimentation, and showing less temporal variation and better convergence with box size. The radial spacing of dense particle filaments is $\sim0.15$ gas scale heights ($H$) for all vBCs, which sets the feeding zone for planetesimal growth in self-gravitating simulations. Our results validate the use of the outflow vBCs in SI simulations, even with vertical boundaries close ($\leq 0.4H$) to the disk midplane. Overall, our study demonstrates the numerical robustness of nonlinear particle clumping by the SI.

We analyze the dust morphology of 29 transition disks (TDs) observed with ALMA at (sub-) millimeter-emission. We perform the analysis in the visibility plane to characterize the total flux, cavity size, and shape of the ring-like structure. First, we found that the $M_{\rm{dust}}-M_\star$ relation is much flatter for TDs than the observed trends from samples of class II sources in different star forming regions. This relation demonstrates that cavities open in high (dust) mass disks, independent of the stellar mass. The flatness of this relation contradicts the idea that TDs are a more evolved set of disks. Two potential reasons (not mutually exclusive) may explain this flat relation: the emission is optically thick or/and millimeter-sized particles are trapped in a pressure bump. Second, we discuss our results of the cavity size and ring width in the context of different physical processes for cavity formation. Photoevaporation is an unlikely leading mechanism for the origin of the cavity of any of the targets in the sample. Embedded giant planets or dead zones remain as potential explanations. Although both models predict correlations between the cavity size and the ring shape for different stellar and disk properties, we demonstrate that with the current resolution of the observations, it is difficult to obtain these correlations. Future observations with higher angular resolution observations of TDs with ALMA will help to discern between different potential origins of cavities in TDs.

Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

Systems with tightly-packed inner planets (STIPs) are very common. Chatterjee & Tan proposed Inside-Out Planet Formation (IOPF), an in situ formation theory, to explain these planets. IOPF involves sequential planet formation from pebble-rich rings that are fed from the outer disk and trapped at the pressure maximum associated with the dead zone inner boundary (DZIB). Planet masses are set by their ability to open a gap and cause the DZIB to retreat outwards. We present models for the disk density and temperature structures that are relevant to the conditions of IOPF. For a wide range of DZIB conditions, we evaluate the gap opening masses of planets in these disks that are expected to lead to truncation of pebble accretion onto the forming planet. We then consider the evolution of dust and pebbles in the disk, estimating that pebbles typically grow to sizes of a few cm during their radial drift from several tens of AU to the inner, $\lesssim1\:$AU-scale disk. A large fraction of the accretion flux of solids is expected to be in such pebbles. This allows us to estimate the timescales for individual planet formation and entire planetary system formation in the IOPF scenario. We find that to produce realistic STIPs within reasonable timescales similar to disk lifetimes requires disk accretion rates of $\sim10^{-9}\:M_\odot\:{\rm yr}^{-1}$ and relatively low viscosity conditions in the DZIB region, i.e., Shakura-Sunyaev parameter of $\alpha\sim10^{-4}$.

Small solids embedded in gaseous protoplanetary disks are subject to strong dust-gas friction. Consequently, tightly-coupled dust particles almost follow the gas flow. This near conservation of dust-to-gas ratio along streamlines is analogous to the near conservation of entropy along flows of (dust-free) gas with weak heating and cooling. We develop this thermodynamic analogy into a framework to study dusty gas dynamics in protoplanetary disks. We show that an isothermal dusty gas behaves like an adiabatic pure gas; and that finite dust-gas coupling may be regarded as an effective heating/cooling. We exploit this correspondence to deduce that 1) perfectly coupled, thin dust layers cannot cause axisymmetric instabilities; 2) radial dust edges are unstable if the dust is vertically well-mixed; 3) the streaming instability necessarily involves a gas pressure response that lags behind dust density; 4) dust-loading introduces buoyancy forces that generally stabilizes the vertical shear instability associated with global radial temperature gradients. We also discuss dusty analogs of other hydrodynamic processes (e.g. Rossby wave instability, convective overstability, and zombie vortices), and how to simulate dusty protoplanetary disks with minor tweaks to existing codes for pure gas dynamics.

Planetesimals may form from the gravitational collapse of dense particle clumps initiated by the streaming instability. We use simulations of aerodynamically coupled gas-particle mixtures to investigate whether the properties of planetesimals formed in this way depend upon the sizes of the particles that participate in the instability. Based on three high resolution simulations that span a range of dimensionless stopping time $6 \times 10^{-3} \leq \tau \leq 2$ no statistically significant differences in the initial planetesimal mass function are found. The mass functions are fit by a power-law, ${\rm d}N / {\rm d}M_p \propto M_p^{-p}$, with $p=1.5-1.7$ and errors of $\Delta p \approx 0.1$. Comparing the particle density fields prior to collapse, we find that the high wavenumber power spectra are similarly indistinguishable, though the large-scale geometry of structures induced via the streaming instability is significantly different between all three cases. We interpret the results as evidence for a near-universal slope to the mass function, arising from the small-scale structure of streaming-induced turbulence.

Hot Jupiters receive strong stellar irradiation, producing equilibrium temperatures of $1000 - 2500 \ \mathrm{Kelvin}$. Incoming irradiation directly heats just their thin outer layer, down to pressures of $\sim 0.1 \ \mathrm{bars}$. In standard irradiated evolution models of hot Jupiters, predicted transit radii are too small. Previous studies have shown that deeper heating -- at a small fraction of the heating rate from irradiation -- can explain observed radii. Here we present a suite of evolution models for HD 209458b where we systematically vary both the depth and intensity of internal heating, without specifying the uncertain heating mechanism(s). Our models start with a hot, high entropy planet whose radius decreases as the convective interior cools. The applied heating suppresses this cooling. We find that very shallow heating -- at pressures of $1 - 10 \ \mathrm{bars}$ -- does not significantly suppress cooling, unless the total heating rate is $\gtrsim 10\%$ of the incident stellar power. Deeper heating, at $100 \ \mathrm{bars}$, requires heating at only $1\%$ of the stellar irradiation to explain the observed transit radius of $1.4 R_{\rm Jup}$ after 5 Gyr of cooling. In general, more intense and deeper heating results in larger hot Jupiter radii. Surprisingly, we find that heat deposited at $10^4 \ \mathrm{bars}$ -- which is exterior to $\approx 99\%$ of the planet's mass -- suppresses planetary cooling as effectively as heating at the center. In summary, we find that relatively shallow heating is required to explain the radii of most hot Jupiters, provided that this heat is applied early and persists throughout their evolution.

Stars form when filaments and dense cores in molecular clouds fragment and collapse due to self-gravity. In the most basic analyses of gravitational stability, the competition between self-gravity and thermal pressure sets the critical (i.e. maximum stable) mass of spheres and the critical line density of cylinders. Previous work has considered additional support from magnetic fields and turbulence. Here, we consider the effects of non-ionizing radiation, specifically the inward radiation pressure force that acts on dense structures embedded in an isotropic radiation field. Using hydrostatic, isothermal models, we find that irradiation lowers the critical mass and line density for gravitational collapse, and can thus act as a trigger for star formation. For structures with moderate central densities, $\sim10^3$ cm$^{-3}$, the interstellar radiation field in the Solar vicinity has an order unity effect on stability thresholds. For more evolved objects with higher central densities, a significant lowering of stability thresholds requires stronger irradiation, as can be found closer to the Galactic center or near stellar associations. Even when strong sources of ionizing radiation are absent or extincted, our study shows that interstellar irradiation can significantly influence the star formation process.

We study the formation of planetesimals in protoplanetary disks from the gravitational collapse of solid over-densities generated via the streaming instability. To carry out these studies, we implement and test a particle-mesh self-gravity module for the Athena code that enables the simulation of aerodynamically coupled systems of gas and collisionless self-gravitating solid particles. Upon employment of our algorithm to planetesimal formation simulations, we find that (when a direct comparison is possible) the Athena simulations yield predicted planetesimal properties that agree well with those found in prior work using different numerical techniques. In particular, the gravitational collapse of streaming-initiated clumps leads to an initial planetesimal mass function that is well-represented by a power-law, dN/dM ~ M^{-p}, with p = 1.6 +/- 0.1, which equates to a differential size distribution dN/dR ~ R^{-q}, with q = 2.8 +/- 0.1. We find no significant trends with resolution from a convergence study of up to 512^3 grid zones and Npar = 1.5x10^8 particles. Likewise, the power-law slope appears indifferent to changes in the relative strength of self-gravity and tidal shear, and to the time when (for reasons of numerical economy) self-gravity is turned on, though the strength of these claims is limited by small number statistics. For a typically assumed radial distribution of minimum mass solar nebula solids (assumed here to have dimensionless stopping time \tau = 0.3), our results support the hypothesis that bodies on the scale of large asteroids or Kuiper Belt Objects could have formed as the high-mass tail of a primordial planetesimal population.

We quantify the thermodynamic requirement for the Vertical Shear Instability and evaluate its relevance to realistic protoplanetary disks as a potential route to hydrodynamic turbulence.

Planetesimals form in gas-rich protoplanetary disks around young stars. However, protoplanetary disks fade in about 10 Myr. The planetesimals (and also many of the planets) left behind are too dim to study directly. Fortunately, collisions between planetesimals produce dusty debris disks. These debris disks trace the processes of terrestrial planet formation for 100 Myr and of exoplanetary system evolution out to 10 Gyr. This chapter begins with a summary of planetesimal formation as a prelude to the epoch of planetesimal destruction. Our review of debris disks covers the key issues, including dust production and dynamics, needed to understand the observations. Our discussion of extrasolar debris keeps an eye on similarities to and differences from Solar System dust.

The vertical shear instability (VSI) offers a potential hydrodynamic mechanism for angular momentum transport in protoplanetary disks (PPDs). The VSI is driven by a weak vertical gradient in the disk's orbital motion, but must overcome vertical buoyancy, a strongly stabilizing influence in cold disks, where heating is dominated by external irradiation. Rapid radiative cooling reduces the effective buoyancy and allows the VSI to operate. We quantify the cooling timescale $t_c$ needed for efficient VSI growth, through a linear analysis of the VSI with cooling in vertically global, radially local disk models. We find the VSI is most vigorous for rapid cooling with $t_c<\Omega_\mathrm{K}^{-1}h|q|/(\gamma -1)$ in terms of the Keplerian orbital frequency, $\Omega_\mathrm{K}$; the disk's aspect-ratio, $h\ll1$; the radial power-law temperature gradient, $q$; and the adiabatic index, $\gamma$. For longer $t_c$, the VSI is much less effective because growth slows and shifts to smaller length scales, which are more prone to viscous or turbulent decay. We apply our results to PPD models where $t_c$ is determined by the opacity of dust grains. We find that the VSI is most effective at intermediate radii, from $\sim5$AU to $\sim50$AU with a characteristic growth time of $\sim30$ local orbital periods. Growth is suppressed by long cooling times both in the opaque inner disk and the optically thin outer disk. Reducing the dust opacity by a factor of 10 increases cooling times enough to quench the VSI at all disk radii. Thus the formation of solid protoplanets, a sink for dust grains, can impede the VSI.

Giant planet formation by core accretion requires a core that is sufficiently massive to trigger runaway gas accretion in less that the typical lifetime of protoplanetary disks. We explore how the minimum required core mass, M_crit, depends on a non-ideal equation of state and on opacity changes due to grain growth, across a range of stellocentric distances from 5-100 AU. This minimum M_crit applies when planetesimal accretion does not substantially heat the atmosphere. Compared to an ideal gas polytrope, the inclusion of molecular hydrogen (H_2) dissociation and variable occupation of H_2 rotational states increases M_crit. Specifically, M_crit increases by a factor of ~2 if the H_2 spin isomers, ortho- and parahydrogen, are in thermal equilibrium, and by a factor of ~2-4 if the ortho-to-para ratio is fixed at 3:1. Lower opacities due to grain growth reduce M_crit. For a standard disk model around a Solar mass star, we calculate M_crit ~ 8 M_Earth at 5 AU, decreasing to ~5 M_Earth at 100 AU, for a realistic EOS with an equilibrium ortho-to-para ratio and for grain growth to cm-sizes. If grain coagulation is taken into account, M_crit may further reduce by up to one order of magnitude. These results for the minimum critical core mass are useful for the interpretation of surveys that find exoplanets at a range of orbital distances.

In the core accretion hypothesis, giant planets form by gas accretion onto solid protoplanetary cores. The minimum (or critical) core mass to form a gas giant is typically quoted as 10 Earth masses. The actual value depends on several factors: the location in the protoplanetary disk, atmospheric opacity, and the accretion rate of solids. Motivated by ongoing direct imaging searches for giant planets, this study investigates core mass requirements in the outer disk. To determine the fastest allowed rates of gas accretion, we consider solid cores that no longer accrete planetesimals, as this would heat the gaseous envelope. Our spherical, two-layer atmospheric cooling model includes an inner convective region and an outer radiative zone that matches onto the disk. We determine the minimum core mass for a giant planet to form within a typical disk lifetime of 3 Myr. The minimum core mass declines with disk radius, from ~8.5 Earth masses at 5 AU to ~3.5 Earth masses at 100 AU, with standard interstellar grain opacities. Lower temperatures in the outer disk explain this trend, while variations in disk density are less influential. At all distances, a lower dust opacity or higher mean molecular weight reduces the critical core mass. Our non-self-gravitating, analytic cooling model reveals that self-gravity significantly affects early atmospheric evolution, starting when the atmosphere is only ~10% as massive as the core.

We present a brief overview of a splinter session on M dwarf stars as planet hosts that was organized as part of the Cool Stars 17 conference. The session was devoted to reviewing our current knowledge of M dwarf stars and exoplanets in order to prepare for current and future exoplanet searches focusing in low mass stars. We review the observational and theoretical challenges to characterize M dwarf stars and the importance of accurate fundamental parameters for the proper characterization of their exoplanets and our understanding on planet formation.

This pedagogical chapter covers the theory of planet formation, with an emphasis on the physical processes relevant to current research. After summarizing empirical constraints from astronomical and geophysical data, we describe the structure and evolution of protoplanetary disks. We consider the growth of planetesimals and of larger solid protoplanets, followed by the accretion of planetary atmospheres, including the core accretion instability. We also examine the possibility that gas disks fragment directly into giant planets and/or brown dwarfs. We defer a detailed description of planet migration and dynamical evolution to other work, such as the complementary chapter in this series by Morbidelli (available at arXiv:1106.4114).

The Pluto system provides a unique local laboratory for the study of binaries with multiple low mass companions. In this paper, we study the orbital stability of P4, the most recently discovered moon in the Pluto system. This newfound companion orbits near the plane of the Pluto-Charon binary, roughly halfway between the two minor moons Nix and Hydra. We use a suite of few body integrations to constrain the masses of Nix and Hydra, and the orbital parameters of P4. For the system to remain stable over the age of the Solar System, the masses of Nix and Hydra likely do not exceed 5e16 kg and 9e16 kg, respectively. These upper limits assume a fixed mass ratio between Nix and Hydra at the value implied by their median optical brightness. Our study finds that stability is more sensitive to their total mass and that a downward revision of Charon's eccentricity (from our adopted value of 0.0035) is unlikely to significantly affect our conclusions. Our upper limits are an order of magnitude below existing astrometric limits on the masses of Nix and Hydra. For a density at least that of ice, the albedos of Nix and Hydra would exceed 0.3. This constraint implies they are icy, as predicted by giant impact models. Even with these low masses, P4 only remains stable if its eccentricity e < 0.02. The 5:1 commensurability with Charon is particularly unstable, Combining stability constraints with the observed mean motion places the preferred orbit for P4 just exterior to the 5:1 resonance. These predictions will be tested when the New Horizons satellite visits Pluto. Based on the results for the Pluto-Charon system, we expect that circumbinary, multi-planet systems will be more widely spaced than their singleton counterparts. Further, circumbinary exoplanets close to the three-body stability boundary, such as those found by Kepler, are less likely to have other companions nearby.

Modelling the formation of super-km-sized planetesimals by gravitational collapse of regions overdense in small particles requires numerical algorithms capable of handling simultaneously hydrodynamics, particle dynamics and particle collisions. While the initial phases of radial contraction are dictated by drag forces and gravity, particle collisions become gradually more significant as filaments contract beyond Roche density. Here we present a new numerical algorithm for treating momentum and energy exchange in collisions between numerical superparticles representing a high number of physical particles. We adopt a Monte Carlo approach where superparticle pairs in a grid cell collide statistically on the physical collision time-scale. Collisions occur by enlarging particles until they touch and solving for the collision outcome, accounting for energy dissipation in inelastic collisions. We demonstrate that superparticle collisions can be consistently implemented at a modest computational cost. In protoplanetary disc turbulence driven by the streaming instability, we argue that the relative Keplerian shear velocity should be subtracted during the collision calculation. If it is not subtracted, density inhomogeneities are too rapidly diffused away, as bloated particles exaggerate collision speeds. Local particle densities reach several thousand times the mid-plane gas density. We find efficient formation of gravitationally bound clumps, with a range of masses corresponding to contracted radii from 100 to 400 km when applied to the asteroid belt and 150 to 730 km when applied to the Kuiper belt, extrapolated using a constant self-gravity parameter. The smaller planetesimals are not observed at low resolution, but the masses of the largest planetesimals are relatively independent of resolution and treatment of collisions.

We develop a general method to fit the planetary distribution function (PLDF) to exoplanet survey data. This maximum likelihood method accommodates more than one planet per star and any number of planet or target star properties. Application to \Kepler data relies on estimates of the efficiency of discovering transits around Solar type stars by Howard et al. (2011). These estimates are shown to agree with theoretical predictions for an ideal transit survey. Using announced \Kepler planet candidates, we fit the PLDF as a joint powerlaw in planet radius, down to 0.5 R_Eart, and orbital period, up to 50 days. The estimated number of planets per star in this sample is ~ 0.7 --- 1.4, where the broad range covers systematic uncertainties in the detection efficiency. To analyze trends in the PLDF we consider four planet samples, divided between shorter and longer periods at 7 days and between large and small radii at 3 R_Earth. At longer periods, the size distribution of the small planets, with index \alpha = -1.2 \pm 0.2 steepens to \alpha = -2.0 \pm 0.2 for the larger planet sample. For shorter periods, the opposite is seen: smaller planets follow a steep powerlaw, \alpha = -1.9 \pm 0.2 that is much shallower, \alpha = -0.7 \pm 0.2 at large radii. The observed deficit of intermediate-sized planets at the shortest periods may arise from the evaporation and sublimation of Neptune and Saturn-like planets. If the trend and explanation hold, it would be spectacular observational confirmation of the core accretion and migration hypotheses, and allow refinement of these theories.

This document covers five broad topics from theory of small bodies: planetesimal formation, cosmochemistry, thermal evolution, collisions and dynamics. Each of these topics is described in a separate section, where we prioritize selected main issues over completeness. The text points out the principal unresolved problems in each area, and suggests ways how progress can be made. This includes support for new code development, observations and experimental work that can be used to constrain theory, new research directions, and studies of cross-over regimes where the system's behavior is determined by several competing processes. The suggested development areas are placed in the context of NASA space exploration.

We study the gravitational instability (GI) of small solids in a gas disk as a mechanism to form planetesimals. Dissipation from gas drag introduces secular GI, which proceeds even when standard GI criteria for a critical density or Toomre's $Q$ predict stability. We include the stabilizing effects of turbulent diffusion, which suppresses small scale GI. The radially wide rings that do collapse contain up to $\sim 0.1$ Earth masses of solids. Subsequent fragmentation of the ring (not modeled here) would produce a clan of chemically homogenous planetesimals. Particle radial drift time scales (and, to a lesser extent, disk lifetimes and sizes) restrict the viability of secular GI to disks with weak turbulent diffusion, characterized by $\alpha \lesssim 10^{-4}$. Thus midplane dead zones are a preferred environment. Large solids with radii $\gtrsim 10$ cm collapse most rapidly because they partially decouple from the gas disk. Smaller solids, even below $\sim$ mm-sizes could collapse if particle-driven turbulence is weakened by either localized pressure maxima or super-Solar metallicity. Comparison with simulations that include particle clumping by the streaming instability shows that our linear model underpredicts rapid, small scale gravitational collapse. Thus the inclusion of more detailed gas dynamics promotes the formation of planetesimals. We discuss relevant constraints from Solar System and accretion disk observations.

The intense irradiation received by hot Jupiters suppresses convection in the outer layers of their atmospheres and lowers their cooling rates. "Inflated" hot Jupiters, i.e., those with anomalously large transit radii, require additional sources of heat or suppressed cooling. We consider the effect of forced turbulent mixing in the radiative layer, which could be driven by atmospheric circulation or by another mechanism. Due to stable stratification in the atmosphere, forced turbulence drives a downward flux of heat. Weak turbulent mixing slows the cooling rate by this process, as if the planet was irradiated more intensely. Stronger turbulent mixing buries heat into the convective interior, provided the turbulence extends to the radiative-convective boundary. This inflates the planet until a balance is reached between the heat buried into and radiated from the interior. We also include the direct injection of heat due to the dissipation of turbulence or other effects. Such heating is already known to slow planetary cooling. We find that dissipation also enhances heat burial from mixing by lowering the threshold for turbulent mixing to drive heat into the interior. Strong turbulent mixing of heavy molecular species such as TiO may be necessary to explain stratospheric thermal inversions. We show that the amount of mixing required to loft TiO may overinflate the planet by our mechanism. This possible refutation of the TiO hypothesis deserves further study. Our inflation mechanism requires a deep stratified layer that only exists when the absorbed stellar flux greatly exceeds the intrinsic emitted flux. Thus it would be less effective for more luminous brown dwarfs and for longer period gas giants, including Jupiter and Saturn.

A large fraction of 100-km-class low-inclination objects in the classical Kuiper Belt (KB) are binaries with comparable mass and wide separation of components. A favored model for their formation was capture during the coagulation growth of bodies in the early KB. Instead, recent studies suggested that large objects can rapidly form in the protoplanetary disks when swarms of locally concentrated solids collapse under their own gravity. Here we examine the possibility that KB binaries formed during gravitational collapse when the excess of angular momentum prevented the agglomeration of available mass into a solitary object. We find that this new mechanism provides a robust path toward the formation of KB binaries with observed properties, and can explain wide systems such as 2001 QW322 and multiples such as (47171) 1999 TC36. Notably, the gravitational collapse is capable of producing 100% binary fraction for a wide range of the swarm's initial angular momentum values. The binary components have similar masses (80% have the secondary-over-primary radius ratio >0.7) and their separation ranges from ~1,000 to ~100,000 km. The binary orbits have eccentricities from e=0 to ~1, with the majority having e<0.6. The binary orbit inclinations with respect to the initial angular momentum of the swarm range from i=0 to ~90 deg, with most cases having i<50 deg. Our binary formation mechanism implies that the primary and secondary components in each binary pair should have identical bulk composition, which is consistent with the current photometric data. We discuss the applicability of our results to the Pluto-Charon, Orcus-Vanth, (617) Patroclus-Menoetius and (90) Antiope binary systems.

Recent direct imaging discoveries suggest a new class of massive, distant planets around A stars. These widely separated giants have been interpreted as signs of planet formation driven by gravitational instability, but the viability of this mechanism is not clear cut. In this paper, we first discuss the local requirements for fragmentation and the initial fragment mass scales. We then consider whether the fragment's subsequent growth can be terminated within the planetary mass regime. Finally, we place disks in the larger context of star formation and disk evolution models. We find that in order for gravitational instability to produce planets, disks must be atypically cold in order to reduce the initial fragment mass. In addition, fragmentation must occur during a narrow window of disk evolution, after infall has mostly ceased, but while the disk is still sufficiently massive to undergo gravitational instability. Under more typical conditions, disk-born objects will likely grow well above the deuterium burning planetary mass limit. We conclude that if planets are formed by gravitational instability, they must be the low mass tail of the distribution of disk-born companions. To validate this theory, on-going direct imaging surveys must find a greater abundance of brown dwarf and M-star companions to A-stars. Their absence would suggest planet formation by a different mechanism such as core accretion, which is consistent with the debris disks detected in these systems.

Planets are built from planetesimals: solids larger than a kilometer which grow by colliding pairwise. Planetesimals themselves are unlikely to form by two-body collisions; sub-km objects have gravitational fields individually too weak, and electrostatic attraction is too feeble for growth beyond a few cm. We review the possibility that planetesimals form when self-gravity brings together vast ensembles of small particles. Even when self-gravity is weak, aerodynamic processes can accumulate solids relative to gas, paving the way for gravitational collapse. Particles pile up as they drift radially inward. Gas turbulence stirs particles, but can also seed collapse by clumping them. While the feedback of solids on gas triggers vertical shear instabilities that obstruct self-gravity, this same feedback triggers streaming instabilities that strongly concentrate particles. Numerical simulations find that solids 10-100 cm in size gravitationally collapse in turbulent disks. We outline areas for progress, including the possibility that still smaller objects self-gravitate.

We present three-dimensional numerical simulations of particle clumping and planetesimal formation in protoplanetary disks with varying amounts of solid material. As centimeter-size pebbles settle to the mid-plane, turbulence develops through vertical shearing and streaming instabilities. We find that when the pebble-to-gas column density ratio is 0.01, corresponding roughly to solar metallicity, clumping is weak, so the pebble density rarely exceeds the gas density. Doubling the column density ratio leads to a dramatic increase in clumping, with characteristic particle densities more than ten times the gas density and maximum densities reaching several thousand times the gas density. This is consistent with unstratified simulations of the streaming instability that show strong clumping in particle dominated flows. The clumps readily contract gravitationally into interacting planetesimals of order 100 km in radius. Our results suggest that the correlation between host star metallicity and exoplanets may reflect the early stages of planet formation. We further speculate that initially low metallicity disks can be particle enriched during the gas dispersal phase, leading to a late burst of planetesimal formation.

We study the behavior of magnetorotational turbulence in shearing box simulations with a radial and azimuthal extent up to ten scale heights. Maxwell and Reynolds stresses are found to increase by more than a factor two when increasing the box size beyond two scale heights in the radial direction. Further increase of the box size has little or no effect on the statistical properties of the turbulence. An inverse cascade excites magnetic field structures at the largest scales of the box. The corresponding 10% variation in the Maxwell stress launches a zonal flow of alternating sub- and super-Keplerian velocity. This in turn generates a banded density structure in geostrophic balance between pressure and Coriolis forces. We present a simplified model for the appearance of zonal flows, in which stochastic forcing by the magnetic tension on short time-scales creates zonal flow structures with life-times of several tens of orbits. We experiment with various improved shearing box algorithms to reduce the numerical diffusivity introduced by the supersonic shear flow. While a standard finite difference advection scheme shows signs of a suppression of turbulent activity near the edges of the box, this problem is eliminated by a new method where the Keplerian shear advection is advanced in time by interpolation in Fourier space.

This pedagogical review covers an unsolved problem in the theory of protoplanetary disks: the growth of dust grains into planetesimals, solids at least a kilometer in size. I summarize timescale constraints imposed on planetesimal formation by circumstellar disk observations, analysis of meteorites, and aerodynamic radial migration. The infall of ~meter-sized solids in a hundred years is the most stringent constraint. I review proposed mechanisms for planetesimal formation. Collisional coagulation models are informed by laboratory studies of microgravity collisions. The gravitational collapse (or Safronov-Goldreich-Ward) hypothesis involves detailed study of the interaction between solid particles and turbulent gas. I cover the basics of aerodynamic drag in protoplanetary disks, including radial drift and vertical sedimentation. I describe various mechanisms for particle concentration in gas disks -- including turbulent pressure maxima, drag instabilities and long-lived anticylonic vortices. I derive a general result for the minimum size for a vortex to trap particles in a sub-Keplerian disk. Recent numerical simulations demonstrate that particle clumping in turbulent protoplanetary disks can trigger gravitational collapse. I discuss several outstanding issues in the field.

Proposed mechanisms for the formation of km-sized solid planetesimals face long-standing difficulties. Robust sticking mechanisms that would produce planetesimals by coagulation alone remain elusive. The gravitational collapse of smaller solids into planetesimals is opposed by stirring from turbulent gas. This proceeding describes recent works showing that "particle feedback," the back-reaction of drag forces on the gas in protoplanetary disks, promotes particle clumping as seeds for gravitational collapse. The idealized streaming instability demonstrates the basic ability of feedback to generate particle overdensities. More detailed numerical simulations show that the particle overdensities produced in turbulent flows trigger gravitational collapse to planetesimals. We discuss surprising aspects of this work, including the large (super-Ceres) mass of the collapsing bound cluster, and the finding that MHD turbulence aids gravitational collapse.

This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs''. It contains 15 sections (\S1.1 -- \S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (\S1.3). Dissipative collisions remove energy from the motion of the particles by collisional cooling (\S1.4), an effect that allows gravitational collapse to occur in somewhat less massive discs (\S1.7). A resolution study of the gravitational collapse of the boulders is presented in \S1.6. We find that gravitational collapse can occur in progressively less massive discs as the grid resolution is increased, likely due to the decreased smoothing of the particle-mesh self-gravity solver with increasing resolution. In \S1.10 we show that it is in good agreement with the Goldreich & Ward (1973) stability analysis to form several-hundred-km-sized bodies, when the analysis is applied to 5 AU and to regions of increased boulder column density. \S11 is devoted to the measurement of random speeds and collision speeds between boulders. We find good agreement between our measurements and analytical theory for the random speeds, but the measured collision speeds are 3 times lower than expected from analytical theory. Higher resolution studies, and an improved analytical theory of collision speeds that takes into account epicyclic motion, will be needed to determine whether collision speeds have converged. In \S1.12 we present models with no magnetic fields. The boulder layer still exhibits strong clumping, due to the streaming instability, if the global solids-to-gas ratio is increased by a factor 3. Gravitational collapse occurs as readily as in magnetised discs.

The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.

We describe the diffusion and random velocities of solid particles due to stochastic forcing by turbulent gas. We include the orbital dynamics of Keplerian disks, both in-plane epicycles and vertical oscillations. We obtain a new result for the diffusion of solids. The Schmidt number (ratio of gas to particle diffusivity) is Sc = 1 + (Omega t_stop)^2, in terms of the particle stopping time, t_stop, and the orbital frequency, Omega. The standard result, Sc = 1 + t_stop/t_eddy, in terms of the eddy turnover time, t_eddy, is shown to be incorrect. The main difference is that Sc rises quadratically, not linearly, with stopping time. Consequently, particles larger than ~ 10 cm in protoplanetary disks will suffer less radial diffusion and will settle closer to the midplane. Such a layer of boulders would be more prone to gravitational collapse. Our predictions of RMS speeds, vertical scale height and diffusion coefficients will help interpret numerical simulations. We confirm previous results for the vertical stirring of particles (scale heights and random velocities), and add a correction for arbitrary ratios of eddy to orbital times. The particle layer becomes thinner for t_eddy > 1/Omega, with the strength of turbulent diffusion held fixed. We use two analytic techniques -- the Hinze-Tchen formalism and the Fokker-Planck equation with velocity diffusion -- with identical results when the regimes of validity overlap. We include simple physical arguments for the scaling of our results.

We present local simulations that verify the linear streaming instability that arises from aerodynamic coupling between solids and gas in protoplanetary disks. This robust instability creates enhancements in the particle density in order to tap the free energy of the relative drift between solids and gas, generated by the radial pressure gradient of the disk. We confirm the analytic growth rates found by Youdin & Goodman (2005) using grid hydrodynamics to simulate the gas and, alternatively, particle and grid representations of the solids. Since the analytic derivation approximates particles as a fluid, this work corroborates the streaming instability when solids are treated as particles. The idealized physical conditions -- axisymmetry, uniform particle size, and the neglect of vertical stratification and collisions -- provide a rigorous, well-defined test of any numerical algorithm for coupled particle-gas dynamics in protoplanetary disks. We describe a numerical particle-mesh implementation of the drag force, which is crucial for resolving the coupled oscillations. Finally we comment on the balance of energy and angular momentum in two-component disks with frictional coupling. A companion paper details the non-linear evolution of the streaming instability into saturated turbulence with dense particle clumps.

We present simulations of the non-linear evolution of streaming instabilities in protoplanetary disks. The two components of the disk, gas treated with grid hydrodynamics and solids treated as superparticles, are mutually coupled by drag forces. We find that the initially laminar equilibrium flow spontaneously develops into turbulence in our unstratified local model. Marginally coupled solids (that couple to the gas on a Keplerian time-scale) trigger an upward cascade to large particle clumps with peak overdensities above 100. The clumps evolve dynamically by losing material downstream to the radial drift flow while receiving recycled material from upstream. Smaller, more tightly coupled solids produce weaker turbulence with more transient overdensities on smaller length scales. The net inward radial drift is decreased for marginally coupled particles, whereas the tightly coupled particles migrate faster in the saturated turbulent state. The turbulent diffusion of solid particles, measured by their random walk, depends strongly on their stopping time and on the solids-to-gas ratio of the background state, but diffusion is generally modest, particularly for tightly coupled solids. Angular momentum transport is too weak and of the wrong sign to influence stellar accretion. Self-gravity and collisions will be needed to determine the relevance of particle overdensities for planetesimal formation.

We analyze the gravitational collapse of solids subject to gas drag in a protoplanetary disk. We also study the stirring of solids by turbulent fluctuations to determine the velocity dispersion and thickness of the midplane particle layer. The usual thresholds for determining gravitational instability in disks, Toomre's criterion and/or the Roche density, do not apply. Dissipation of angular momentum allows instability at longer wavelengths, lower densities, and higher velocity dispersions than without drag. Small solids will slowly leak into axisymmetric rings since initial collapse occurs over many orbits. Growth is fastest when particle stopping times are comparable to orbital times. Our analysis of particle stirring by turbulence is consistent with previous results for tightly coupled particles, but is generalized to loose coupling where epicyclic motions contribute to random velocities. A companion paper applies these results to turbulent protoplanetary disks.

We show that small solids in low mass, turbulent protoplanetary disks collect into self-gravitating rings. Growth is faster than disk lifetimes and radial drift times for moderately strong turbulence, characterized by dimensionless diffusivities, $\alpha_g < 10^{-6} -- 10^{-3}$ when particles are mm-sized. This range reflects a strong dependance on disk models. Growth is faster for higher particle surface densities. Lower gas densities and larger solids also give faster growth, as long as aerodynamic coupling is tight. In simple power law models, growth is slowest around ~0.3 AU, where drag coupling is strongest for mm-sized solids. Growth is much faster close to the star where orbital times are short, with implications for in situ formation of short period extrasolar planets. Growth times also decrease toward the outer disk where lower gas densities allow greater particle settling. Beyond roughly Kuiper Belt distances however, solids are sufficiently decoupled from gas that dissipative gravitational instabilities are less effective. Turbulence not only slows growth, but also increases radial wavelengths. The initial solid mass in an unstable ring can be ~0.01 M_Earth or greater, huge compared to km-sized planetesimals. Nonlinear fragmentation, which has not been studied in detail, will lower the final planetesimal mass. We consider applications to the asteroid belt and discuss the alternate hypothesis of collisional agglomeration.

Interpenetrating streams of solids and gas in a Keplerian disk produce a local, linear instability. The two components mutually interact via aerodynamic drag, which generates radial drift and triggers unstable modes. The secular instability does not require self-gravity, yet it generates growing particle density perturbations that could seed planetesimal formation. Growth rates are slower than dynamical, but faster than radial drift, timescales. Growth rates, like streaming velocities, are maximized for marginal coupling (stopping times comparable dynamical times). Fastest growth occurs when the solid to gas density ratio is order unity and feedback is strongest. Curiously, growth is strongly suppressed when the densities are too nearly equal. The relation between background drift and wave properties is explained by analogy with Howard's semicircle theorem. The three-dimensional, two-fluid equations describe a sixth order (in the complex frequency) dispersion relation. A terminal velocity approximation allows simplification to an approximate cubic dispersion relation. To describe the simplest manifestation of this instability, we ignore complicating (but possibly relevant) factors like vertical stratification, dispersion of particle sizes, turbulence, and self-gravity. We consider applications to planetesimal formation and compare our work to other studies of particle-gas dynamics.

Solid particles in protoplanetary disks that are sufficiently super-solar in metallicity overcome turbulence generated by vertical shear to gravitationally condense into planetesimals. Super-solar metallicities result if solid particles pile up as they migrate starward as a result of aerodynamic drag. Previous analyses of aerodynamic drift rates that account for mean flow differences between gas and particles yield particle pile-ups. We improve on these studies not only by accounting for the collective inertia of solids relative to that of gas, but also by including the transport of angular momentum by turbulent stresses within the particle layer. These turbulent stresses are derived in a physically self-consistent manner from the structure of marginally Kelvin-Helmholtz turbulent flows. They are not calculated using the usual plate drag formulae, whose use we explain is inappropriate. Accounting for the relative inertia of solids to gas retards, but does not prevent, particle pile-ups, and generates more spatially extended regions of metal enrichment. Turbulent transport hastens pile-ups. We conclude that particle pile-up is a robust outcome in sufficiently passive protoplanetary disks. Connections to observations of circumstellar disks, including the Kuiper Belt, and the architectures of planetary systems are made.

The debate over whether kilometer-sized solids, or planetesimals, assemble by collision-induced chemical sticking or by gravity-driven unstable modes remains unsettled. In light of recent work showing that gravitational growth can occur despite turbulent stirring, we critically evaluate the collisional hypothesis. Impact speeds in protoplanetary disks reach 50 m/s in a laminar disk and may be larger in turbulent disks. We consider the role of elastic and plastic deformations, restructuring of "rubble piles," and sticky organic matter. Coagulation is possible for dust grains with large surface area-to-volume ratios. Simple energetic arguments show that bouncing, cratering, and fragmentation should dominate collisional dynamics between millimeter and kilometer sizes.

We investigate the formation of planetesimals via the gravitational instability of solids that have settled to the midplane of a circumstellar disk. Vertical shear between the gas and a subdisk of solids induces turbulent mixing which inhibits gravitational instability. Working in the limit of small, well-coupled particles, we find that the mixing becomes ineffective when the surface density ratio of solids to gas exceeds a critical value. Solids in excess of this precipitation limit can undergo midplane gravitational instability and form planetesimals. However, this saturation effect typically requires increasing the local ratio of solid to gaseous surface density by factors of two to ten times cosmic abundances, depending on the exact properties of the gas disk. We discuss existing astrophysical mechanisms for augmenting the ratio of solids to gas in protoplanetary disks by such factors, and investigate a particular process that depends on the radial variations of orbital drift speeds induced by gas drag. This mechanism can concentrate millimeter sized chondrules to the supercritical surface density in < few x 10^6 years, a suggestive timescale for the disappearance of dusty disks in T Tauri stars. We discuss the relevance of our results to some outstanding puzzles in planet formation theory: the size of the observed solar system, and the rapid type I migration of Earth mass bodies.