No Arabic abstract
The dynamical stability of tightly packed exoplanetary systems remains poorly understood. While for a two-planet system a sharp stability boundary exists, numerical simulations of three and more planet systems show that they can experience instability on timescales up to billions of years. Moreover, an exponential trend between the planet orbital separation measured in units of Hill radii and the survival time has been reported. While these findings have been observed in numerous numerical simulations, little is known of the actual mechanism leading to instability. Contrary to a constant diffusion process, planetary systems seem to remain dynamically quiescent for most of their lifetime before a very short unstable phase. In this work, we show how the slow chaotic diffusion due to the overlap of three-body resonances dominates the timescale leading to the instability for initially coplanar and circular orbits. While the last instability phase is related to scattering due to two-planet mean motion resonances (MMR), for circular orbits the two-planets MMR are too far separated to destabilize systems initially away from them. We develop an analytical model to generalize the empirical trend obtained for equal mass and equally-spaced planets to general systems. We obtain an analytical estimate of the survival time consistent with simulations over four orders of magnitude for the planet to star mass ratio $epsilon$, and 6 to 8 orders of magnitude for the instability time. We also confirm that measuring the orbital spacing in terms of Hill radii is not adapted and that the right spacing unit scales as $epsilon^{1/4}$. We predict that beyond a certain spacing, the three-planet resonances are not overlapped, which results in an increase of the survival time. We finally discuss the extension of our result to more general systems, containing more planets on initially non circular orbits.
Transit surveys have revealed a significant population of compact multi-planet systems, containing several sub-Neptune-mass planets on close-in, tightly-packed orbits. These systems are thought to have formed through a final phase of giant impacts, which would tend to leave systems close to the edge of stability. Here, we assess this hypothesis, comparing observed eccentricities in systems exhibiting transit-timing variations (TTVs), with the maximum eccentricities compatible with long-term stability. We use the machine-learning classifier SPOCK (Tamayo et al. 2020) to rapidly classify the stability of numerous initial configurations and hence determine these stability limits. While previous studies have argued that multi-planet systems are often maximally packed, in the sense that they could not host any additional planets, we find that the existing planets in these systems have measured eccentricities below the limits allowed by stability by a factor of 2--10. We compare these results against predictions from the giant impact theory of planet formation, derived from both $N$-body integrations and theoretical expectations that in the absence of dissipation, the orbits of such planets should be distributed uniformly throughout the phase space volume allowed by stability. We find that the observed systems have systematically lower eccentricities than this scenario predicts, with a median eccentricity about 4 times lower than predicted. These findings suggest that if such systems formed through giant impacts, then some dissipation must occur to damp their eccentricities. This may take place during formation, perhaps through interactions with the natal gas disk or a leftover population of planetesimals, or over longer timescales through the coupling of tidal and secular processes.
Despite over three hundred years of effort, no solutions exist for predicting when a general planetary configuration will become unstable. We introduce a deep learning architecture to push forward this problem for compact systems. While current machine learning algorithms in this area rely on scientist-derived instability metrics, our new technique learns its own metrics from scratch, enabled by a novel internal structure inspired from dynamics theory. Our Bayesian neural network model can accurately predict not only if, but also when a compact planetary system with three or more planets will go unstable. Our model, trained directly from short N-body time series of raw orbital elements, is more than two orders of magnitude more accurate at predicting instability times than analytical estimators, while also reducing the bias of existing machine learning algorithms by nearly a factor of three. Despite being trained on compact resonant and near-resonant three-planet configurations, the model demonstrates robust generalization to both non-resonant and higher multiplicity configurations, in the latter case outperforming models fit to that specific set of integrations. The model computes instability estimates up to five orders of magnitude faster than a numerical integrator, and unlike previous efforts provides confidence intervals on its predictions. Our inference model is publicly available in the SPOCK package, with training code open-sourced.
Excess emission, associated with warm, dust belts, commonly known as exozodis, has been observed around a third of nearby stars. The high levels of dust required to explain the observations are not generally consistent with steady-state evolution. A common suggestion is that the dust results from the aftermath of a dynamical instability, an event akin to the Solar Systems Late Heavy Bombardment. In this work, we use a database of N-body simulations to investigate the aftermath of dynamical instabilities between giant planets in systems with outer planetesimal belts. We find that, whilst there is a significant increase in the mass of material scattered into the inner regions of the planetary system following an instability, this is a short-lived effect. Using the maximum lifetime of this material, we determine that even if every star has a planetary system that goes unstable, there is a very low probability that we observe more than a maximum of 1% of sun-like stars in the aftermath of an instability, and that the fraction of planetary systems currently in the aftermath of an instability is more likely to be limited to <0.06. This probability increases marginally for younger or higher mass stars. We conclude that the production of warm dust in the aftermath of dynamical instabilities is too short-lived to be the dominant source of the abundantly observed exozodiacal dust.
We revisit the discovery and implications of the first candidate systems to contain multiple transiting exoplanets. These systems were discovered using data from the Kepler space telescope. The initial paper, presenting five systems (Steffen et al. 2010), was posted online at the time the project released the first catalog of Kepler planet candidates. The first extensive analysis of the observed population of multis was presented in a follow-up paper published the following year (Lissauer et al. 2011a). Multiply-transiting systems allow us to answer a variety of important questions related to the formation and dynamical evolution of planetary systems. These two papers addressed a wide array of topics including: the distribution of orbital period ratios, planet size ratios, system architectures, mean-motion resonance, orbital eccentricities, planet validation and confirmation, and the identification of different planet populations. They set the stage for many subsequent, detailed studies by other groups. Intensive studies of individual multiplanet systems provided some of Keplers most important exoplanet discoveries. As we examine the scientific impact of the first of these systems, we also present some history of the people and circumstances surrounding their discoveries.
We announce the discovery of a planetary system with 7 transiting planets around a Kepler target, a current record for transiting systems. Planets b, c, e and f are reported for the first time in this work. Planets d, g and h were previously reported in the literature (Batalha et al. 2013), although here we revise their orbital parameters and validate their planetary nature. Planets h and g are gas giants and show strong dynamical interactions. The orbit of planet g is perturbed in such way that its orbital period changes by 25.7h between two consecutive transits during the length of the observations, which is the largest such perturbation found so far. The rest of the planets also show mutual interactions: planets d, e and f are super-Earths close to a mean motion resonance chain (2:3:4), and planets b and c, with sizes below 2 Earth radii, are within 0.5% of the 4:5 mean motion resonance. This complex system presents some similarities to our Solar System, with small planets in inner orbits and gas giants in outer orbits. It is, however, more compact. The outer planet has an orbital distance around 1 AU, and the relative position of the gas giants is opposite to that of Jupiter and Saturn, which is closer to the expected result of planet formation theories. The dynamical interactions between planets are also much richer.