skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: An Analytical Theory for the Growth from Planetesimals to Planets by Polydisperse Pebble Accretion

Pebble accretion is recognized as a significant accelerator of planet formation. Yet only formulae for single-sized (monodisperse) distribution have been derived in the literature. These can lead to significant underestimates for Bondi accretion, for which the best accreted pebble size may not be the one that dominates the mass distribution. We derive in this paper the polydisperse theory of pebble accretion. We consider a power-law distribution in pebble radius, and we find the resulting surface and volume number density distribution functions. We derive also the exact monodisperse analytical pebble accretion rate for which 3D accretion and 2D accretion are limits. In addition, we find analytical solutions to the polydisperse 2D Hill and 3D Bondi limits. We integrate the polydisperse pebble accretion numerically for the MRN distribution, finding a slight decrease (by an exact factor 3/7) in the Hill regime compared to the monodisperse case. In contrast, in the Bondi regime, we find accretion rates 1–2 orders of magnitude higher compared to monodisperse, also extending the onset of pebble accretion to 1–2 orders of magnitude lower in mass. We find megayear timescales, within the disk lifetime, for Bondi accretion on top of planetary seeds of masses 10−6to 10−4M, over a significant range of the parameter space. This mass range overlaps with the high-mass end of the planetesimal initial mass function, and thus pebble accretion is possible directly following formation by streaming instability. This alleviates the need for mutual planetesimal collisions as a major contribution to planetary growth.

more » « less
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Page Range / eLocation ID:
Article No. 60
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    Despite many methods developed to find young massive planets in protoplanetary discs, it is challenging to directly detect low-mass planets that are embedded in discs. On the other hand, the core-accretion theory suggests that there could be a large population of embedded low-mass young planets at the Kelvin-Helmholtz (KH) contraction phase. We adopt both 1D models and 3D simulations to calculate the envelopes around low-mass cores (several to tens of M⊕) with different luminosities, and derive their thermal fluxes at radio wavelengths. We find that, when the background disc is optically thin at radio wavelengths, radio observations can see through the disc and probe the denser envelope within the planet’s Hill sphere. When the optically thin disc is observed with the resolution reaching one disc scale height, the radio thermal flux from the planetary envelope around a 10 M⊕ core is more than 10 per cent higher than the flux from the background disc. The emitting region can be extended and elongated. Finally, our model suggests that the au-scale clump at 52 au in the TW Hydrae disc revealed by ALMA is consistent with the envelope of an embedded 10–20 M⊕ planet, which can explain the detected flux, the spectral index dip, and the tentative spirals. The observation is also consistent with the planet undergoing pebble accretion. Future ALMA and ngVLA observations may directly reveal more such low-mass planets, enabling us to study core growth and even reconstruct the planet formation history using the embedded ‘protoplanet’ population.

    more » « less

    The origin of Uranus and Neptune has long been challenging to explain, due to the large orbital distances from the Sun. After a planetary embryo has been formed, the main accretion processes are likely pebble, gas, and planetesimal accretion. Previous studies of Uranus and Neptune formation typically do not consider all three processes; and furthermore, do not investigate how the formation of the outer planet impacts the inner planet. In this paper, we study the concurrent formation of Uranus and Neptune via pebble, gas, and planetesimal accretion. We use a dust-evolution model to predict the size and mass flux of pebbles, and derive our own fit for gas accretion. We do not include migration, but consider a wide range of formation locations between 12 and $40\, \textrm {au}$. If the planetary embryos form at the same time and with the same mass, our formation model with an evolving dust population is unable to produce Uranus and Neptune analogues. This is because the mass difference between the planets and the H–He mass fractions become too high. However, if the outer planetary embryo forms earlier and/or more massive than the inner embryo, the two planets do form in a few instances when the disc is metal-rich and dissipates after a few Myr. Furthermore, our study suggests that in situ formation is rather unlikely. Nevertheless, giant impacts and/or migration could potentially aid in the formation, and future studies including these processes could bring us one step closer to understanding how Uranus and Neptune formed.

    more » « less
  3. Abstract

    Formation models in which terrestrial bodies grow via the pairwise accretion of planetesimals have been reasonably successful at reproducing the general properties of the Solar System, including small-body populations. However, planetesimal accretion has not yet been fully explored in the context of the wide variety of recently discovered extrasolar planetary systems, in particular those that host short-period terrestrial planets. In this work, we use directN-body simulations to explore and understand the growth of planetary embryos from planetesimals in disks extending down to ≃1 day orbital periods. We show that planetesimal accretion becomes nearly 100% efficient at short orbital periods, leading to embryo masses that are much larger than the classical isolation mass. For rocky bodies, the physical size of the object begins to occupy a significant fraction of its Hill sphere toward the inner edge of the disk. In this regime, most close encounters result in collisions, rather than scattering, and the system does not develop a bimodal population of dynamically hot planetesimals and dynamically cold oligarchs, as is seen in previous studies. The highly efficient accretion seen at short orbital periods implies that systems of tightly packed inner planets should be almost completely devoid of any residual small bodies. We demonstrate the robustness of our results to assumptions about the initial disk model, and we also investigate the effects that our simplified collision model has on the emergence of this non-oligarchic growth mode in a planet-forming disk.

    more » « less

    Gaps imaged in protoplanetary discs are suspected to be opened by planets. We compute the present-day mass accretion rates $\dot{M}_{\rm p}$ of seven hypothesized gap-embedded planets, plus the two confirmed planets in the PDS 70 disc. The accretion rates are based on disc gas surface densities Σgas from C18O observations, and planet masses Mp from simulations fitted to observed gaps. Assuming accretion is Bondi-like, we find in eight out of nine cases that $\dot{M}_{\rm p}$ is consistent with the time-averaged value given by the current planet mass and system age, Mp/tage. As system ages are comparable to circumstellar disc lifetimes, these gap-opening planets may be undergoing their last mass doublings, reaching final masses of $M_{\rm p} \sim 10\rm{\!-\!}10^2 \, M_\oplus$ for the non-PDS 70 planets, and $M_{\rm p} \sim 1\!-\!10 \, M_{\rm J}$ for the PDS 70 planets. For another 15 gaps without C18O data, we predict Σgas by assuming their planets are accreting at their time-averaged $\dot{M}_{\rm p}$. Bondi accretion rates for PDS 70b and c are orders of magnitude higher than accretion rates implied by measured U-band and H α fluxes, suggesting most of the accretion shock luminosity emerges in as yet unobserved wavebands, or that the planets are surrounded by dusty, highly extincting, quasi-spherical circumplanetary envelopes. Thermal emission from such envelopes or from circumplanetary discs, on Hill sphere scales, peaks at wavelengths in the mid-to-far-infrared and can reproduce observed mm-wave excesses.

    more » « less

    Secular oscillations in multiplanet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This ‘secular chaos’ readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging method, with emphasis on the relationship between the planets’ properties and the time-scale and efficiency of chaotic diffusion. We find that secular chaos can excite extreme eccentricities on time-scales spanning several orders of magnitude in a given system. We apply our results to the evolution of planetary systems around white dwarfs (WDs), specifically the tidal disruption and high-eccentricity migration of planetesimals and planets. We find that secular chaos in a planetesimal belt driven by large (≳10 M⊕), distant ($\gtrsim 10 \, \mathrm{au}$) planets can sustain metal accretion on to a WD over Gyr time-scales. We constrain the total mass of planetesimals initially present within the chaotic zone by requiring that the predicted mass delivery rate to the Roche limit be consistent with the observed metal accretion rates of WDs with atmospheric pollution throughout the cooling sequence. Based on the occurrence of long-period exoplanets and exo-asteroid belts, we conclude that secular chaos can be a significant (perhaps dominant) channel for polluting solitary WDs. Secular chaos can also produce short-period planets and planetesimals around WDs in concert with various circularization mechanisms. We discuss prospects for detecting exoplanets driving secular chaos around WDs using direct imaging and microlensing.

    more » « less