No Arabic abstract
We derive a semi-analytic criterion for the presence of chaos in compact, eccentric multiplanet systems. Beyond a minimum semimajor-axis separation, below which the dynamics are chaotic at all eccentricities, we show that (i) the onset of chaos is determined by the overlap of two-body mean motion resonances (MMRs), like it is in two-planet systems; (ii) secular evolution causes the MMR widths to expand and contract adiabatically, so that the chaotic boundary is established where MMRs overlap at their greatest width. For closely spaced two-planet systems, a near-symmetry strongly suppresses this secular modulation, explaining why the chaotic boundaries for two-planet systems are qualitatively different from cases with more than two planets. We use these results to derive an improved angular-momentum-deficit (AMD) stability criterion, i.e., the critical system AMD below which stability should be guaranteed. This introduces an additional factor to the expression from Laskar and Petit (2017) that is exponential in the interplanetary separations, which corrects the AMD threshold toward lower eccentricities by a factor of several for tightly packed configurations. We make routines for evaluating the chaotic boundary available to the community through the open-source SPOCK package.
We derive a criterion for the onset of chaos in systems consisting of two massive, eccentric, coplanar planets. Given the planets masses and separation, the criterion predicts the critical eccentricity above which chaos is triggered. Chaos occurs where mean motion resonances overlap, as in Wisdom (1980)s pioneering work. But whereas Wisdom considered only nearly circular planets, and hence examined only first order resonances, we extend his results to arbitrarily eccentric planets (up to crossing orbits) by examining resonances of all orders. We thereby arrive at a simple expression for the critical eccentricity. We do this first for a test particle in the presence of a planet, and then generalize to the case of two massive planets, based on a new approximation to the Hamiltonian (Hadden, in prep). We then confirm our results with detailed numerical simulations. Finally, we explore the extent to which chaotic two-planet systems eventually result in planetary collisions.
We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.
We review the occurrence of the patterns of the onset of chaos in low-dimensional nonlinear dissipative systems in leading topics of condensed matter physics and complex systems of various disciplines. We consider the dynamics associated with the attractors at period-doubling accumulation points and at tangent bifurcations to describe features of glassy dynamics, critical fluctuations and localization transitions. We recall that trajectories pertaining to the routes to chaos form families of time series that are readily transformed into networks via the Horizontal Visibility algorithm, and this in turn facilitates establish connections between entropy and Renormalization Group properties. We discretize the replicator equation of game theory to observe the onset of chaos in familiar social dilemmas, and also to mimic the evolution of high-dimensional ecological models. We describe an analytical framework of nonlinear mappings that reproduce rank distributions of large classes of data (including Zipfs law). We extend the discussion to point out a common circumstance of drastic contraction of configuration space driven by the attractors of these mappings. We mention the relation of generalized entropy expressions with the dynamics along and at the period doubling, intermittency and quasi-periodic routes to chaos. Finally, we refer to additional natural phenomena in complex systems where these conditions may manifest.
Understanding the relationship between long-period giant planets and multiple smaller short-period planets is critical for formulating a complete picture of planet formation. This work characterizes three such systems. We present Kepler-65, a system with an eccentric (e=0.28+/-0.07) giant planet companion discovered via radial velocities (RVs) exterior to a compact, multiply-transiting system of sub-Neptune planets. We also use precision RVs to improve mass and radius constraints on two other systems with similar architectures, Kepler-25 and Kepler-68. In Kepler-68 we propose a second exterior giant planet candidate. Finally, we consider the implications of these systems for planet formation models, particularly that the moderate eccentricity in Kepler-65s exterior giant planet did not disrupt its inner system.
We study the stability of a family of spherical equilibrium models of self-gravitating systems, the so-called $gamma-$models with Osipkov-Merritt velocity anisotropy, by means of $N-$body simulations. In particular, we analyze the effect of self-consistent $N-$body chaos on the onset of radial-orbit instability (ROI). We find that degree of chaoticity of the system associated to its largest Lyapunov exponent $Lambda_{rm max}$ has no appreciable relation with the stability of the model for fixed density profile and different values of radial velocity anisotropy. However, by studying the distribution of the Lyapunov exponents $lambda_{rm m}$ of the individual particles in the single-particle phase space, we find that more anisotropic systems have a larger fraction of orbits with larger $lambda_{rm m}$.