skip to main content


Title: Planetesimal Accretion at Short Orbital Periods
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
Award ID(s):
2006752
NSF-PAR ID:
10445107
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
954
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 61
Size(s):
["Article No. 61"]
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Using the tree-based N-body code ChaNGa, we explore the consequences of planetesimal-planetesimal accretion in the context of systems of tightly-packed inner planets (STIPs). Formation models for STIPs often begin with the assumption that the building blocks for the planets previously underwent a phase of oligarchic growth, in which a handful of large bodies develop and leave behind a dynamically hot population of residual planetesimals. We show that the conditions required for oligarchic growth become difficult to achieve close to the star, where the dynamical and orbital timescales are short. In this case, planetary embryos can form much larger than the typical isolation mass. We explore a range of physically motivated initial conditions for the starting distribution of planetesimals and find that a number of different outcomes for the resulting embryo and residual planetesimal populations are possible. In some cases, a transition region develops partway out in the disk, the location of which serves as an artifact of the original planetesimal distribution. 
    more » « less
  2. Using high-resolution N-body simulations, we investigate the outcome of terrestrial planet formation at short (< 100 day) orbital periods under a migration-free model. The collisional and dynamical evolution of systems of nearly 106 self-interacting planetesimals are directly followed through the final planet assembly phase. This is done by first modeling the planetesimal evolution with the tree-based N-body code ChaNGa, and then passing the results to the hybrid-symplectic N-body code genga, once the particle count has dropped sufficiently. Previously, we showed that oligarchic growth fails to operate at arbitrarily short orbital periods. This leaves a distinct feature in the mass and orbital distribution of the planetary embryos. In this most recent work, we explore whether this boundary between oligarchic and non-oligarchic growth leaves any kind of imprint on the terrestrial planets that form. If so, this would provide an important clue to evaluate whether migration played a significant role in shaping the architecture systems of tightly-packed inner planets. 
    more » « less
  3. null (Ed.)
    During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modeling. Calculated solid/molten alloy partitioning of C increases greatly with liquid S concentration, and inferred parent body C concentrations range from 0.0004 to 0.11 wt%. Parent bodies fall into two compositional clusters characterized by cores with medium and low C/S. Both of these require significant planetesimal degassing, as metamorphic devolatilization on chondrite-like precursors is insufficient to account for their C depletions. Planetesimal core formation models, ranging from closed-system extraction to degassing of a wholly molten body, show that significant open-system silicate melting and volatile loss are required to match medium and low C/S parent body core compositions. Greater depletion in C relative to S is the hallmark of silicate degassing, indicating that parent body core compositions record processes that affect composite silicate/iron planetesimals. Degassing of bare cores stripped of their silicate mantles would deplete S with negligible C loss and could not account for inferred parent body core compositions. Devolatilization during small-body differentiation is thus a key process in shaping the volatile inventory of terrestrial planets derived from planetesimals and planetary embryos. 
    more » « less
  4. ABSTRACT

    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
  5. ABSTRACT

    A planetary system can undergo multiple episodes of intense dynamical activities throughout its life, resulting in the production of star-grazing planetesimals (or exocomets) and pollution of the host star. Such activity is especially pronounced when giant planets interact with other small bodies during the system’s evolution. However, due to the chaotic nature of the dynamics, it is difficult to determine the properties of the perturbing planet(s) from the observed planetesimal-disruption activities. In this study, we examine the outcomes of planetesimal-planet scatterings in a general setting. We focus on one-planet systems, and determine the likelihood and time-scale of planetesimal disruption by the host star as a function of the planet properties. We obtain a new analytical expression for the minimum distance a scattering body can reach, extending previous results by considering finite planet eccentricity and non-zero planetesimal mass. Through N-body simulations, we derive the distribution of minimum distances and the likelihood and time-scales of three possible outcomes of planetesimal-planet scatterings: collision with the planet, ejection, and disruption by the star. For planetesimals with negligible mass, we identify four defining dimensionless parameters (the planet eccentricity, planet-to-star mass ratio, planet radius to semimajor axis ratio, and the stellar disruption radius to planet semimajor axis ratio) that enable us to scale the problem and generalize our findings to a wide range of orbital configurations. Using these results, we explore three applications: falling evaporating bodies in the β Pictoris system, white dwarf pollution due to planetesimal disruption and planet engulfment by main-sequence stars.

     
    more » « less