BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew Youdin

Andrew Youdin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-phastro-ph.EP

Catalog footprint

What is connected

8works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2020arXiv

Demographics of Planetesimals Formed by the Streaming Instability

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Adding particle collisions to the formation of asteroids and Kuiper belt objects via streaming instabilities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Particle Clumping and Planetesimal Formation Depend Strongly on Metallicity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Zonal Flows and Long-Lived Axisymmetric Pressure Bumps in Magnetorotational Turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

From Grains to Planetesimals: Les Houches Lecture

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Protoplanetary Disk Turbulence Driven by the Streaming Instability: Linear Evolution and Numerical Methods

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Protoplanetary Disk Turbulence Driven by the Streaming Instability: Non-Linear Saturation and Particle Concentration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Rapid planetesimal formation in turbulent circumstellar 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.

0 citations0 reviewsNo code linkedOpen access