Pebble Accretion Scenario Research Articles

The bouncing barrier revisited: Impact on key planet formation processes and observational signatures

Context. A leading paradigm in planet formation is currently the streaming instability and pebble accretion scenario. Notably, dust must grow into sizes in a specific regime of Stokes numbers in order to make the processes in the scenario viable and sufficiently effective. The dust growth models currently in use do not implement some of the growth barriers suggested to be relevant in the literature. Aims. We investigate if the bouncing barrier, when effective, has an impact on the timescales and efficiencies of processes such as the streaming instability and pebble accretion as well as on the observational appearance of planet-forming disks. Methods. We implemented a formalism for the bouncing barrier into the publicly available dust growth model DustPy and ran a series of models to understand the impact. Results. We found that the bouncing barrier has a significant effect on the dust evolution in planet-forming disks. In many cases, it reduces the size of the typical or largest particles available in the disk; it produces a very narrow, almost monodisperse, size distribution; and it removes most μm-sized grains in the process, with an impact on scattered light images. It modifies the settling and therefore the effectiveness of and timescales for the streaming instability and for pebble accretion. An active bouncing barrier may well have observational consequences: It may reduce the strength of the signatures of small particles (e.g., the 10 μm silicate feature), and it may create additional shadowed regions visible in scattered light images. Conclusions. Modeling of planet formation that leans heavily on the streaming instability and on pebble accretion should take the bouncing barrier into account. The complete removal of small grains in our model is not consistent with observations. However, this could be resolved by incomplete vertical mixing or some level of erosion in collisions.

Composition of giant planets: The roles of pebbles and planetesimals

One of the current challenges of planet formation theory is to explain the enrichment of observed exoplanetary atmospheres. While past studies have focused on scenarios where either pebbles or planetesimals are the main drivers of heavy element enrichment, here we combine the two approaches to understand whether the composition of a planet can constrain its formation pathway. We study three different formation scenarios: pebble accretion, pebble accretion with planetesimal formation inside the disc, and combined pebble and planetesimal accretion. We used the CHEMCOMP code to perform semi-analytical 1D simulations of protoplanetary discs, including viscous evolution, pebble drift, and simple chemistry to simulate the growth of planets from planetary embryos to gas giants as they migrate through the disc, while simultaneously tracking their composition. Our simulations confirm that the composition of the planetary atmosphere is dominated by the accretion of gas vapour enriched by inward-drifting and evaporating pebbles. Including planetesimal formation hinders this enrichment because many pebbles are locked into planetesimals and cannot evaporate and enrich the disc. This results in a dramatic drop in accreted heavy elements in the cases of planetesimal formation and accretion, demonstrating that planetesimal formation needs to be inefficient in order to explain planets with high heavy element content. On the other hand, accretion of planetesimals enhances the refractory component of the atmosphere, leading to low volatile-to-refractory ratios in the atmosphere, in contrast to the majority of pure pebble simulations. However, low volatile-to-refractory ratios can also be achieved in the pure pebble accretion scenario if the planet migrates all the way into the inner disc and accretes gas that is enriched in evaporated refractories. Distinguishing these two scenarios requires knowledge about the planet’s atmospheric C/H and O/H ratios, which are much higher in the pure pebble scenario compared to the planetesimal formation and accretion scenario. This implies that a detailed knowledge of the composition of planetary atmospheres could help to distinguish between the different formation scenarios.

Global N-body simulations of circumbinary planet formation around Kepler-16 and -34 analogues I: Exploring the pebble accretion scenario

ABSTRACT Numerous circumbinary planets have been discovered in surveys of transiting planets. Often, these planets are found to orbit near to the zone of dynamical instability, close to the central binary. The existence of these planets has been explained by hydrodynamical simulations that show that migrating circumbinary planets, embedded in circumbinary discs, halt at the central cavity that is formed by the central binary. Transit surveys are naturally most sensitive to finding circumbinary planets with the shortest orbital periods. The future promise of detecting longer period systems using radial-velocity searches, combined with the anticipated detection of numerous circumbinary planets by ESA’s PLATO mission, points to the need to model and understand the formation and evolution of circumbinary planets in a more general sense than has been considered before. With this goal in mind, we present a newly developed global model of circumbinary planet formation that is based on the mercury6 symplectic N-body integrator, combined with a model for the circumbinary disc and prescriptions for a range of processes involved in planet formation such as pebble accretion, gas envelope accretion, and migration. Our results show that under reasonable assumptions, the pebble accretion scenario can produce circumbinary systems that are similar to those observed, and in particular is able to produce planets akin to Kepler-16b and Kepler-34b. Comparing our results to other systems, we find that our models also adequately reproduce such systems, including multiplanet systems. Resonances between neighbouring planets are frequently obtained, whilst ejections of planets by the central binary acts as an effective source of free floating planets.

Dry or water world? How the water contents of inner sub-Neptunes constrain giant planet formation and the location of the water ice line

In the pebble accretion scenario, the pebbles that form planets drift inward from the outer disk regions, carrying water ice with them. At the water ice line, the water ice on the inward drifting pebbles evaporates and is released into the gas phase, resulting in water-rich gas and dry pebbles that move into the inner disk regions. Large planetary cores can block the inward drifting pebbles by forming a pressure bump outside their orbit in the protoplanetary disk. Depending on the relative position of a growing planetary core relative to the water ice line, water-rich pebbles might be blocked outside or inside the water ice line. Pebbles blocked outside the water ice line do not evaporate and thus do not release their water vapor into the gas phase, resulting in a dry inner disk, while pebbles blocked inside the water ice line release their water vapor into the gas phase, resulting in water vapor diffusing into the inner disk. As a consequence, close-in sub-Neptunes that accrete some gas from the disk should be dry or wet, respectively, if outer gas giants are outside or inside the water ice line, assuming that giant planets form fast, as has been suggested for Jupiter in our Solar System. Alternatively, a sub-Neptune could form outside the water ice line, accreting a large amount of icy pebbles and then migrating inward as a very wet sub-Neptune. We suggest that the water content of inner sub-Neptunes in systems with giant planets that can efficiently block the inward drifting pebbles could constrain the formation conditions of these systems, thus making these sub-Neptunes exciting targets for detailed characterization (e.g., with JWST, ELT, or ARIEL). In addition, the search for giant planets in systems with already characterized sub-Neptunes can be used to constrain the formation conditions of giant planets as well.

How planets grow by pebble accretion

During their formation, planets form large, hot atmospheres due to the ongoing accretion of solids. It has been customary to assume that all solids end up at the center, constituting a “core” of refractory materials, whereas the envelope remains metal-free. However, recent work, as well as observations by the Juno mission, indicate that the distinction may not be so clear cut. Indeed, small silicate, pebble-sized particles will sublimate in the atmosphere when they hit the sublimation temperature (T ~ 2000 K). In this paper we extend previous analytical work to compute the properties of planets within such a pebble accretion scenario. We conduct 1D numerical calculations of the atmosphere of an accreting planet, solving the stellar structure equations, augmented by a nonideal equation of state that describes a hydrogen and helium-silicate vapor mixture. Calculations terminate at the point where the total mass in metal is equal to that of the H+He gas, which we numerically confirm as the onset of runaway gas accretion. When pebbles sublimate before reaching the core, insufficient (accretion) energy is available to mix dense, vapor-rich lower layers with the higher layers of lower metallicity. A gradual structure in which Z decreases with radius is therefore a natural outcome of planet formation by pebble accretion. We highlight, furthermore, that (small) pebbles can act as the dominant source of opacity, preventing rapid cooling and presenting a channel for (mini-)Neptunes to survive in gas-rich disks. Nevertheless, once pebble accretion subsides, the atmosphere rapidly clears followed by runaway gas accretion. We consider atmospheric recycling to be the most probable mechanism to have stalled the growth of the envelopes of these planets.

Probing the impact of varied migration and gas accretion rates for the formation of giant planets in the pebble accretion scenario

ABSTRACT The final orbital position of growing planets is determined by their migration speed, which is essentially set by the planetary mass. Small mass planets migrate in type-I migration, while more massive planets migrate in type-II migration, which is thought to depend mostly on the viscous evolution rate of the disc. A planet is most vulnerable to inward migration before it reaches type-II migration and can lose a significant fraction of its semimajor axis at this stage. We investigated the influence of different disc viscosities, the dynamical torque, and gas accretion from within the horseshoe region as mechanisms for slowing down planet migration. Our study confirms that planets growing in low viscosity environments migrate less, due to the earlier gap opening and slower type-II migration rate. We find that taking the gas accretion from the horseshoe region into account allows an earlier gap opening and this results in less inward migration of growing planets. Furthermore, this effect increases the planetary mass compared to simulations that do not take the effect of gas accretion from the horseshoe region. Moreover, combining the effect of the dynamical torque with the effect of gas accretion from the horseshoe region, significantly slows down inward migration. Taking these effects into account could allow the formation of cold Jupiters (a > 1 au) closer to the water ice line region compared to previous simulations that did not take these effects into account. We, thus, conclude that gas accretion from within the horseshoe region and the dynamical torque play crucial roles in shaping planetary systems.

The eccentricity distribution of giant planets and their relation to super-Earths in the pebble accretion scenario

Observations of the population of cold Jupiter planets (r >1 AU) show that nearly all of these planets orbit their host star on eccentric orbits. For planets up to a few Jupiter masses, eccentric orbits are thought to be the outcome of planet–planet scattering events taking place after gas dispersal. We simulated the growth of planets via pebble and gas accretion as well as the migration of multiple planetary embryos in their gas disc. We then followed the long-term dynamical evolution of our formed planetary system up to 100 Myr after gas disc dispersal. We investigated the importance of the initial number of protoplanetary embryos and different damping rates of eccentricity and inclination during the gas phase for the final configuration of our planetary systems. We constrained our model by comparing the final dynamical structure of our simulated planetary systems to that of observed exoplanet systems. Our results show that the initial number of planetary embryos has only a minor impact on the final orbital eccentricity distribution of the giant planets, as long as the damping of eccentricity and inclination is efficient. If the damping is inefficient (slow), systems with a larger initial number of embryos harbour larger average eccentricities. In addition, for slow damping rates, we observe that scattering events are already common during the gas disc phase and that the giant planets that formed in these simulations match the observed giant planet eccentricity distribution best. These simulations also show that massive giant planets (above Jupiter mass) on eccentric orbits are less likely to host inner super-Earths as they get lost during the scattering phase, while systems with less massive giant planets on nearly circular orbits should harbour systems of inner super-Earths. Finally, our simulations predict that giant planets are not single, on average, but they live in multi-planet systems.

Planet formation and disc mass dependence in a pebble-driven scenario for low-mass stars

ABSTRACT Measured disc masses seem to be too low to form the observed population of planetary systems. In this context, we develop a population synthesis code in the pebble accretion scenario, to analyse the disc mass dependence on planet formation around low-mass stars. We base our model on the analytical sequential model presented by Ormel, Liu, and Schoonenberg and analyse the populations resulting from varying initial disc mass distributions. Starting out with seeds the mass of Ceres formed by streaming instability inside the ice-line, we grow the planets using the pebble accretion process and migrate them inwards using type I migration. The next planets are formed sequentially after the previous planet crosses the ice line. We explore different initial distributions of disc masses to show the dependence of this parameter with the final planetary population. Our results show that compact close-in resonant systems can be pretty common around M dwarfs between 0.09 and 0.2 M⊙ only when the discs considered are more massive than what is being observed by sub-mm disc surveys. The minimum disc mass to form a Mars-like planet is found to be about 2 × 10−3 M⊙. Small variations in the disc mass distribution also manifest in the simulated planet distribution. The paradox of disc masses might be caused by an underestimation of the disc masses in observations, by a rapid depletion of mass in discs by planets growing within 1 million years, or by deficiencies in our current planet formation picture.

Survival of Primordial Planetary Atmospheres: Photodissociation-driven Mass Loss

Abstract The most widely studied mechanism of mass loss from extrasolar planets is photoevaporation via X-ray and ultraviolet ionization, primarily in the context of highly irradiated planets. However, the extreme ultraviolet dissociation of hydrogen molecules can also theoretically drive atmospheric evaporation on low-mass planets. For temperate planets such as the early Earth, impact erosion is expected to dominate in the traditional planetesimal accretion model, but it would be greatly reduced in pebble accretion scenarios, allowing other mass-loss processes to be major contributors. We apply the same prescription for photoionization to this photodissociation mechanism and compare it to an analysis of other possible sources of mass loss in pebble accretion scenarios. We find that there is no clear path to evaporating the primordial atmosphere accreted by an early Earth analog in a pebble accretion scenario. Impact erosion could remove ∼2300 bars of hydrogen if 1% of the planet’s mass is accreted as planetesimals, while the combined photoevaporation processes could evaporate ∼750 bars of hydrogen. Photodissociation is likely a subdominant but significant component of mass loss. Similar results apply to super-Earths and mini-Neptunes. This mechanism could also preferentially remove hydrogen from a planet’s primordial atmosphere, thereby leaving a larger abundance of primordial water compared to standard dry formation models. We discuss the implications of these results for models of rocky planet formation, including Earth’s formation, and the possible application of this analysis to mass loss from observed exoplanets.

Influence of migration models and thermal torque on planetary growth in the pebble accretion scenario

Low-mass planets that are in the process of growing larger within protoplanetary disks exchange torques with the disk and change their semi-major axis accordingly. This process is called type I migration and is strongly dependent on the underlying disk structure. As a result, there are many uncertainties about planetary migration in general. In a number of simulations, the current type I migration rates lead to planets reaching the inner edge of the disk within the disk lifetime. A new kind of torque exchange between planet and disk, the thermal torque, aims to slow down inward migration via the heating torque. The heating torque may even cause planets to migrate outwards, if the planetary luminosity is large enough. Here, we study the influence on planetary migration of the thermal torque on top of previous type I models. We find that the formula of Paardekooper et al. (2011, MNRAS, 410, 293) allows for more outward migration than that of Jiménez & Masset (2017, MNRAS, 471, 4917) in most configurations, but we also find that planets evolve to very similar mass and final orbital radius using both formulae in a single planet-formation scenario, including pebble and gas accretion. Adding the thermal torque can introduce new, but small, regions of outwards migration if the accretion rates onto the planet correspond to typical solid accretion rates following the pebble accretion scenario. If the accretion rates onto the planets become very large, as could be the case in environments with large pebble fluxes (e.g., high-metallicity environments), the thermal torque can allow more efficient outward migration. However, even then, the changes for the final mass and orbital positions in our planet formation scenario are quite small. This implies that for single planet evolution scenarios, the influence of the heating torque is probably negligible.

The Galilean Satellites Formed Slowly from Pebbles

Abstract It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter’s formation. However, understanding the specifics of the formation process is challenging, as both small particles (pebbles) and the satellites are subject to fast migration processes. Here we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede’s mass is determined by pebble isolation. Our slow pebble accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; and (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

Inner rocky super-Earth formation: distinguishing the formation pathways in viscously heated and passive discs

Observations have revealed that super-Earths (planets up to 10 Earth masses) are the most abundant type of planets in the inner systems. Their formation is strongly linked to the structure of the protoplanetary disc, which determines growth and migration. In the pebble accretion scenario, planets grow to the pebble isolation mass, at which the planet carves a small gap in the gas disc halting the pebble flux and thus its growth. The pebble isolation mass scales with the disc’s aspect ratio, which directly depends on the heating source of the protoplanetary disc. I compare the growth of super-Earths in viscously heated discs, where viscous heating dissipates within the first million years, and discs purely heated by the central star with super-Earth observations from the Kepler mission. This allows two formation pathways of super-Earths to be distinguished in the inner systems within this framework. Planets growing within 1 Myr in the viscously heated inner disc reach pebble isolation masses that correspond directly to the inferred masses of the Kepler observations for systems that feature planets in resonance or not in resonance. However, to explain the period ratio distribution of Kepler planets – where most Kepler planet pairs are not in mean motion resonance configurations – a fraction of these resonant chains has to be broken. In case the planets are born early in a viscously heated disc, these resonant chains thus have to be broken without planetary mergers, for example through the magnetic rebound effect, and the final system architecture should feature low mutual inclinations. If super-Earths form either late or in purely passive discs, the pebble isolation mass is too small (around 2–3 Earth masses) to explain the Kepler observations, implying that planetary mergers have to play a significant role in determining the final system architecture. Resonant planetary systems thus have to experience mergers already during the gas disc phase, so the planets can get trapped in resonance after reaching 5–10 Earth masses. In case instabilities are dominating the system architecture, the systems should not be flat, but feature mutually inclined orbits. This implies that future observations of planetary systems with radial velocities (RV) and transits (for example through the Transiting Exoplanet Survey Satellite (TESS) and its follow up RV surveys) could distinguish between these two formation channels of super-Earth and thus constrain planet formation theories.

Atmospheric signatures of giant exoplanet formation by pebble accretion

Atmospheric chemical abundances of giant planets lead to important constraints on planetary formation and migration. Studies have shown that giant planets that migrate through the protoplanetary disk can accrete substantial amounts of oxygen-rich planetesimals, leading to super-solar metallicities in the envelope and solar or sub-solar C/O ratios. Pebble accretion has been demonstrated to play an important role in core accretion and to have growth rates that are consistent with planetary migration. The high pebble accretion rates allow planetary cores to start their growth beyond 10 AU and subsequently migrate to cold (>~ 1 AU), warm (~0.1 AU- 1AU) or hot (<~ 0.1 AU) orbits. In this work we investigate how the formation of giant planets via pebble accretion influences their atmospheric chemical compositions. We find that under the standard pebble accretion scenario, where the core is isolated from the envelope, the resulting metallicities (O/H and C/H ratios) are sub-solar, while the C/O ratios are super-solar. Planets that migrate through the disk to become hot Jupiters accrete substantial amounts of water vapour, but still acquire slightly sub-solar O/H and super-solar C/O of 0.7-0.8. The metallicity can be substantially sub-solar (~0.2-0.5x solar) and the C/O can even approach 1.0 if the planet accretes its envelope mostly beyond the CO2 ice line, i.e. cold Jupiters or hot Jupiters that form far out and migrate in by scattering. Allowing for core erosion yields significantly super-solar metallicities and solar or sub-solar C/O, which can also be achieved by other means, e.g. photo-evaporation and late-stage planetesimal accretion.