skip to main content


Title: Secular chaos in white dwarf planetary systems: origins of metal pollution and short-period planetary companions
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
NSF-PAR ID:
10367292
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
513
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
p. 4178-4195
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A planet’s orbital alignment places important constraints on how a planet formed and consequently evolved. The dominant formation pathway of ultra-short-period planets (P < 1 day) is particularly mysterious as such planets most likely formed further out, and it is not well understood what drove their migration inwards to their current positions. Measuring the orbital alignment is difficult for smaller super-Earth/sub-Neptune planets, which give rise to smaller amplitude signals. Here we present radial velocities across two transits of 55 Cancri (Cnc) e, an ultra-short-period super-Earth, observed with the Extreme Precision Spectrograph. Using the classical Rossiter–McLaughlin method, we measure 55 Cnc e’s sky-projected stellar spin–orbit alignment (that is, the projected angle between the The star 55 Cancri (Cnc) A hosts five known exoplanets with minimum mass estimates ranging from approximately 8M⊕ to 3MJup and periods less than one day to nearly 20 years1–4. Of particular interest has been 55 Cnc e, one of the most massive known ultra-short-period planets (USPs) and the only planet around 55 Cnc found to transit5,6. It has an star’s spin axis and the planet’s orbit normal—will shed light on the formation and evolution of USPs, especially in the case of compact, multiplanet systems. It has been shown that USPs form a statistically distinct popula- tion of planets9 that tend to be misaligned with other planetary orbits in their system10. This suggests that USPs experience a unique migra- tion pathway that brings them close in to their host stars. This inward migration is most likely driven by dissipation due to star–planet tidal interactions that result from either non-zero eccentricities11,12 or plan- etary spin-axis tilts13. orbital period of 0.7365474 +1.3 × 10−6 days, a mass of 7.99 ± 0.33M −1.4 × 10−6 ⊕ and a radius of 1.853 +0.026 R⊕ (refs. 7,8). A precise measure of the −0.027 stellar spin–orbit alignment of 55 Cnc e—the angle between the host planet’s orbital axis and its host star’s spin axis) to be λ = 10 +17∘ with an +14∘ −20∘ unprojected angle of ψ = 23 −12∘. The best-fit Rossiter–McLaughlin model to the Extreme Precision Spectrograph data has a radial velocity semi- amplitude of just 0.41 +0.09 m s−1. The spin–orbit alignment of 55 Cnc e −0.10 favours dynamically gentle migration theories for ultra-short-period planets, namely tidal dissipation through low-eccentricity planet–planet interactions and/or planetary obliquity tides. 
    more » « less
  2. 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
  3. ABSTRACT

    A large fraction of white dwarfs (WDs) have metal-polluted atmospheres, which are produced by accreting material from remnant planetary systems. The composition of the accreted debris broadly resembles that of rocky Solar system objects. Volatile-enriched debris with compositions similar to long-period comets (LPCs) is rarely observed. We attempt to reconcile this dearth of volatiles with the premise that exo-Oort clouds (XOCs) occur around a large fraction of planet-hosting stars. We estimate the comet accretion rate from an XOC analytically, adapting the ‘loss cone’ theory of LPC delivery in the Solar system. We investigate the dynamical evolution of an XOC during late stellar evolution. Using numerical simulations, we show that 1–30 per cent of XOC objects remain bound after anisotropic stellar mass-loss imparting a WD natal kick of ${\sim}1 \, {\rm km \, s^{-1}}$. We also characterize the surviving comets’ distribution function. Surviving planets orbiting a WD can prevent the accretion of XOC comets by the star. A planet’s ‘dynamical barrier’ is effective at preventing comet accretion if the energy kick imparted by the planet exceeds the comet’s orbital binding energy. By modifying the loss cone theory, we calculate the amount by which a planet reduces the WD’s accretion rate. We suggest that the scarcity of volatile-enriched debris in polluted WDs is caused by an unseen population of 10–$100 \, \mathrm{au}$ scale giant planets acting as barriers to incoming LPCs. Finally, we constrain the amount of volatiles delivered to a planet in the habitable zone of an old, cool WD.

     
    more » « less
  4. ABSTRACT

    Diversity in the properties of exoplanetary systems arises, in part, from dynamical evolution that occurs after planet formation. We use numerical integrations to explore the relative role of secular and resonant dynamics in the long-term evolution of model planetary systems, made up of three equal mass giant planets on initially eccentric orbits. The range of separations studied is dominated by secular processes, but intersects chains of high-order mean-motion resonances. Over time-scales of 108 orbits, the secular evolution of the simulated systems is predominantly regular. High-order resonant chains, however, can be a significant source of angular momentum deficit (AMD), leading to instability. Using a time series analysis based on a Hilbert transform, we associate instability with broad islands of chaotic evolution. Previous work has suggested that first-order resonances could modify the AMD of nominally secular systems and facilitate secular chaos. We find that higher order resonances, when present in chains, can have similar impacts.

     
    more » « less
  5. Abstract

    The dynamical evolution of the solar system is chaotic with a Lyapunov time of only ∼5 Myr for the inner planets. Due to the chaos it is fundamentally impossible to accurately predict the solar system’s orbital evolution beyond ∼50 Myr based on present astronomical observations. We have recently developed a method to overcome the problem by using the geologic record to constrain astronomical solutions in the past. Our resulting optimal astronomical solution (called ZB18a) shows exceptional agreement with the geologic record to ∼58 Ma (Myr ago) and a characteristic resonance transition around 50 Ma. Here we show that ZB18a and integration of Earth’s and Mars’ spin vector based on ZB18a yield reduced variations in Earth’s and Mars’ orbital inclination and Earth’s obliquity (axial tilt) from ∼58 to ∼48 Ma—the latter being consistent with paleoclimate records. The changes in the obliquities have important implications for the climate histories of Earth and Mars. We provide a detailed analysis of solar system frequencies (gandsmodes) and show that the shifts in the variation in Earth’s and Mars’ orbital inclination and obliquity around 48 Ma are associated with the resonance transition and caused by changes in the contributions to the superposition ofsmodes, plusgsmode interactions in the inner solar system. Thegsmode interactions and the resonance transition (consistent with geologic data) are unequivocal manifestations of chaos. Dynamical chaos in the solar system hence not only affects its orbital properties but also the long-term evolution of planetary climate through eccentricity and the link between inclination and axial tilt.

     
    more » « less