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.
The dynamical interaction of minor bodies (such as comets or asteroids) with planets plays an essential role in the planetary system’s architecture and evolution. As a result of these interactions, structures such as the Kuiper belt and the Oort cloud can be created. In particular, the collision of minor bodies with planets can drastically change the planet’s internal and orbital evolution. We present an analytical formulation to determine the collision time-scale for a minor body to impact a planet for arbitrary geometry. By comparing with a suite of detailed N-body simulations and an analytical method for collision time-scales in the Solar system, we confirmed the accuracy of our analytical formulation. As a proof of concept, we focused on the collision time-scales of minor bodies similar to the Jupiter-family comets and the long-period comets with a Jupiter-like planet. We show that our analytical method yields in good agreement with the numerical simulations. The formalism presented here thus provides a succinct and accurate alternative to numerical calculations.
more » « less- NSF-PAR ID:
- 10430130
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 524
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 1025-1030
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT We present a fragmentation module and a composition tracking code for the N-body code rebound. Our fragmentation code utilizes previous semi-analytical models and follows an implementation method similar to fragmentation for the N-body code mercury. In our N-body simulations with fragmentation, we decrease the collision and planet formation time-scales by inflating the particle radii by an expansion factor f and experiment with various values of f to understand how expansion factors affect the collision history and final planetary system. As the expansion factor increases, so do the rate of mergers which produces planetary systems with more planets and planets at larger orbits. Additionally, we present a composition tracking code which follows the compositional change of homogeneous bodies as a function of mass exchange and use it to study how fragmentation and the use of an expansion factor affects volatile delivery to the inner terrestrial disc. We find that fragmentation enhances radial mixing relative to perfect merging and that on average, as f increases so does the average water mass fraction of the planets. Radial mixing decreases with increasing f as collisions happen early on, before the bodies have time to grow to excited orbits and move away from their original location.
-
null (Ed.)Context. Dynamical models of Solar System evolution have suggested that the so-called P- and D-type volatile-rich asteroids formed in the outer Solar System beyond Neptune’s orbit and may be genetically related to the Jupiter Trojans, comets, and small Kuiper belt objects (KBOs). Indeed, the spectral properties of P- and D-type asteroids resemble that of anhydrous cometary dust. Aims. We aim to gain insights into the above classes of bodies by characterizing the internal structure of a large P- and D-type asteroid. Methods. We report high-angular-resolution imaging observations of the P-type asteroid (87) Sylvia with the Very Large Telescope Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. These images were used to reconstruct the 3D shape of Sylvia. Our images together with those obtained in the past with large ground-based telescopes were used to study the dynamics of its two satellites. We also modeled Sylvia’s thermal evolution. Results. The shape of Sylvia appears flattened and elongated (a/b ~1.45; a/c ~1.84). We derive a volume-equivalent diameter of 271 ± 5 km and a low density of 1378 ± 45 kg m −3 . The two satellites orbit Sylvia on circular, equatorial orbits. The oblateness of Sylvia should imply a detectable nodal precession which contrasts with the fully-Keplerian dynamics of its two satellites. This reveals an inhomogeneous internal structure, suggesting that Sylvia is differentiated. Conclusions. Sylvia’s low density and differentiated interior can be explained by partial melting and mass redistribution through water percolation. The outer shell should be composed of material similar to interplanetary dust particles (IDPs) and the core should be similar to aqueously altered IDPs or carbonaceous chondrite meteorites such as the Tagish Lake meteorite. Numerical simulations of the thermal evolution of Sylvia show that for a body of such a size, partial melting was unavoidable due to the decay of long-lived radionuclides. In addition, we show that bodies as small as 130–150 km in diameter should have followed a similar thermal evolution, while smaller objects, such as comets and the KBO Arrokoth, must have remained pristine, which is in agreement with in situ observations of these bodies. NASA Lucy mission target (617) Patroclus (diameter ≈140 km) may, however, be differentiated.more » « less
-
more » « less
Abstract Polluted white dwarfs (WDs) offer a unique way to study the bulk compositions of exoplanetary material, but it is not always clear if this material originates from comets, asteroids, moons, or planets. We combine
N -body simulations with an analytical model to assess the prevalence of extrasolar moons as WD polluters. Using a sample of observed polluted WDs, we find that the extrapolated parent body masses of the polluters are often more consistent with those of many solar system moons, rather than solar-like asteroids. We provide a framework for estimating the fraction of WDs currently undergoing observable moon accretion based on results from simulated WD planetary and moon systems. Focusing on a three-planet WD system of super-Earth to Neptune-mass bodies, we find that we could expect about one percent of such systems to be currently undergoing moon accretions as opposed to asteroid accretion. -
more » « less
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.
