Source author record

Andrew N. Youdin

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

19 published item(s)

preprint2022arXiv

The 3D dust and opacity distribution of protoplanets in multi-fluid global simulations

The abundance and distribution of solids inside the Hill sphere are central to our understanding of the giant planet dichotomy. Here, we present a three-dimensional characterization of the dust density, mass flux, and mean opacities in the envelope of sub-thermal and super-thermal mass planets. We simulate the dynamics of multiple dust species in a global protoplanetary disk model accounting for dust feedback. We find that the meridional flows do not effectively stir dust grains at scales of the Bondi sphere. Thus the dust-settling driven by the stellar gravitational potential sets the latitudinal dust density gradient within the planet envelope. Not only does the planet's potential enhance this gradient, but also the spiral wakes serve as another source of asymmetry. These asymmetries substantially alter the inferred mean Rosseland and Planck opacities. In cases with the moderate-to-strong dust settling, the opacity gradient can range from a few percent to more than two orders of magnitude between the mid-plane and the polar regions of the Bondi sphere. Finally, we show that this strong latitudinal opacity gradient can introduce a transition between optically thick and thin regimes at the scales of the planet envelope. We suggest that this transition is likely to occur when the equilibrium scale height of hundred-micron-sized particles is smaller than the Hill radius of the forming planet. This work calls into question the adoption of a constant opacity derived from well-mixed distributions and demonstrates the need for global radiation hydrodynamics models of giant planet formation which account for dust dynamics.

preprint2021arXiv

Planetesimal formation around the snow line. II. Dust or pebbles?

Around the snow line, icy pebbles and silicate dust may locally pile-up and form icy and rocky planetesimals via streaming instability and/or gravitational instability. We perform 1D diffusion-advection simulations that include the back-reaction to radial drift and diffusion of icy pebbles and silicate dust, ice sublimation, release of silicate dust, and their recycling through recondensation and sticking onto pebbles outside the snow line. We use a realistic description of the scale height of silicate dust obtained from Ida et al. and that of pebbles including the effects of a Kelvin-Helmholtz instability. We study the dependence of solid pile-up on distinct effective viscous parameters for turbulent diffusions in the radial and vertical directions ($α_{\rm Dr}$ and $α_{\rm Dz}$) and for the gas accretion to the star ($α_{\rm acc}$) as well as that on the pebble-to-gas mass flux ($F_{\rm p/g}$). We derive the sublimation width of drifting icy pebbles which is a critical parameter to characterize the pile-up of silicate dust and pebbles around the snow line. We identify a parameter space (in the $F_{\rm p/g}-α_{\rm acc}-α_{\rm Dz}(=α_{\rm Dr})$ space) where pebbles no longer drift inward to reach the snow line due to the back-reaction that slows down radial velocity of pebbles. We show that the pile-up of solids around the snow line occurs in a broader range of parameters for $α_{\rm acc}=10^{-3}$ than for $α_{\rm acc}=10^{-2}$. Above a critical $F_{\rm p/g}$ value, the runaway pile-up of silicate dust inside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \ll 1$, while that of pebbles outside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \sim 1$. Our results imply that a distinct evolutionary path could produce a diversity of outcomes in terms of planetesimal formation around the snow line.

preprint2020arXiv

Dust Settling Instability in Protoplanetary Discs

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 $α\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.

preprint2020arXiv

Planetesimal formation around the snow line: I. Monte Carlo simulations of silicate dust pile-up in a turbulent disk

Context: The formation of rocky planetesimals is a long-standing problem in planet formation theory. One of the possibilities is that it results from gravitational instability as a result of pile-up of small silicate dust particles released from sublimating icy pebbles that pass the snow line. Aims: We want to understand and quantify the role of the water snow line for the formation of rock-rich and ice-rich planetesimals. In this paper, we focus on the formation of rock-rich planetesimals. A companion paper examines the combined formation of both rock-rich and ice-rich planetesimals. Methods: We develop a new Monte Carlo code to calculate the radial evolution of silicate particles in a turbulent accretion disk, accounting for the back-reaction (i.e., inertia) of the particles on their radial drift velocity and diffusion. Results depend in particular on the particle injection width (determined from the radial sublimation width of icy pebbles), the pebble scale height and the pebble mass flux through the disk. The scale height evolution of the silicate particles, which is the most important factor for the runaway pile-up, is automatically calculated in this Lagrange method. Results: From the numerical results, we derive semi-analytical relations for the scale height of the silicate dust particles and the particles-to-gas density ratio at the midplane, as functions of a pebble-to-gas mass flux ratio and the $α$ parameters for disk gas accretion and vertical/radial diffusion. We find that the runaway pile-up of the silicate particles (formation of rocky planetesimals) occurs if the pebble-to-gas mass flux ratio is $> [(α_{Dz}/α_{acc})/3 \times 10^{-2}]^{1/2}$ where $α_{Dz}$ and $α_{acc}$ are the $α$ parameters for vertical turbulent diffusion and disk gas accretion.

preprint2020arXiv

Superparticle Method for Simulating Collisions

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.

preprint2017arXiv

Structure and Evolution of Internally Heated Hot Jupiters

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.

preprint2016arXiv

The Mass and Size Distribution of Planetesimals Formed by the Streaming Instability. I. The Role of Self-Gravity

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 τ= 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.

preprint2016arXiv

The Role of Non-ionizing Radiation Pressure in Star Formation: The Stability of Cores and Filaments

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.

preprint2015arXiv

Cooling Requirements for the Vertical Shear Instability in Protoplanetary Disks

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<Ω_\mathrm{K}^{-1}h|q|/(γ-1)$ in terms of the Keplerian orbital frequency, $Ω_\mathrm{K}$; the disk's aspect-ratio, $h\ll1$; the radial power-law temperature gradient, $q$; and the adiabatic index, $γ$. 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.

preprint2015arXiv

Planetesimals in Debris Disks

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.

preprint2014arXiv

Minimum Core Masses for Giant Planet Formation With Realistic Equations of State and Opacities

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.

preprint2014arXiv

On the Minimum Core Mass for Giant Planet Formation at Wide Separations

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.

preprint2012arXiv

Circumbinary Chaos: Using Pluto's Newest Moon to Constrain the Masses of Nix & Hydra

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.

preprint2012arXiv

From Disks to Planets

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).

preprint2011arXiv

On the Formation of Planetesimals via Secular Gravitational Instabilities with Turbulent Stirring

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 $α\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.

preprint2011arXiv

The Exoplanet Census: A General Method, Applied to Kepler

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 α= -1.2 \pm 0.2 steepens to α= -2.0 \pm 0.2 for the larger planet sample. For shorter periods, the opposite is seen: smaller planets follow a steep powerlaw, α= -1.9 \pm 0.2 that is much shallower, α= -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.

preprint2010arXiv

The Mechanical Greenhouse: Burial of Heat by Turbulence in Hot Jupiter Atmospheres

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.

preprint2010arXiv

The Runts of the Litter: Why planets formed through gravitational instability can only be failed binary stars

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.