Search NASA⌕ Search

SEARCH · Search NASA

Results for “Planetary systems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

A Machine Learns to Predict the Stability of Highly Diverse Multi-planetary Systems [Slides]

We’ve discovered over 4,000 exoplanets in more than 500 systems. Why do none of these systems resemble our Solar System? How do planetary systems organize themselves? How do planets’ orbits evolve over time? What causes orbital instability? What are indicators of instability? Evolution of high multiplicity planetary systems is an analytically unsolved problem. Assessing the stability of these systems is usually done through computationally expensive N-body simulations. The report provides further details about the author's research, methods and conclusions.

79 ASTRONOMY AND ASTROPHYSICS↗

The Dark Planets of the WASP-47 Planetary System

Exoplanet discoveries have demonstrated a vast range of planetary system architectures. The demographic of compact planetary systems are especially interesting from the perspective of planetary formation and the evolution of orbital dynamics. Another interesting demographic is that of giant planets in eccentric orbits, since these planets have likely had a dynamical history involving planet–planet scattering events. The WASP-47 system is particularly fascinating since it combines these two demographics, having both compact planetary orbits and a giant planet on an eccentric orbit within the system the habitable zone. Here we provide an analysis of the WASP-47 system from the perspective of atmospheric detection and characterization. We discuss the system architecture and the potential for additional long-period planets. We simulate expected phase variations as a function of planet orbital phase for the system due to the combined effect of the planets. We present an analysis of precision photometry of WASP-47 from the K2 mission, phased on each of the planets. The analysis rules out the detection of phase signatures for the two inner-most planets, enabling constraints upon their albedos and atmospheric properties. Our study concludes that WASP-47b is an example of a “dark” planet with a tentative geometric albedo of 0.016 and a 1σ upper limit of 0.17. The WASP-47e data are consistent with a broad range of albedos, but also show early evidence of having a relatively low albedo. The growing number of dark, short-period giant planets provide the framework of an ideal sample for studying low albedo dependence on atmospheric composition.

79 ASTRONOMY AND ASTROPHYSICS↗

Two Young Planetary Systems around Field Stars with Ages between 20 and 320 Myr from TESS

Planets around young stars trace the early evolution of planetary systems. We report the discovery and validation of two planetary systems with ages ≲300 Myr from observations by the Transiting Exoplanet Survey Satellite (TESS). The 40--320 Myr old G star TOI-251 hosts a 2.74{sub −0.18}{sup +0.18} R{sub ⊕} mini-Neptune with a 4.94 day period. The 20--160 Myr old K star TOI-942 hosts a system of inflated Neptune-sized planets, with TOI-942b orbiting in a period of 4.32 days with a radius of 4.81{sub −0.20}{sup +0.20} R{sub ⊕} and TOI-942c orbiting in a period of 10.16 days with a radius of 5.79{sub −0.18}{sup +0.19} R{sub ⊕}. Though we cannot place either host star into a known stellar association or cluster, we can estimate their ages via their photometric and spectroscopic properties. Both stars exhibit significant photometric variability due to spot modulation, with measured rotation periods of ∼3.5 days. These stars also exhibit significant chromospheric activity, with age estimates from the chromospheric calcium emission lines and X-ray fluxes matching that estimated from gyrochronology. Both stars also exhibit significant lithium absorption, similar in equivalent width to well-characterized young cluster members. TESS has the potential to deliver a population of young planet-bearing field stars, contributing significantly to tracing the properties of planets as a function of their age.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Machine Learns to Predict the Stability of Highly Diverse Multi-planetary Systems

Machine learning has proven to be an invaluable tool for characterizing the stability of planets in simplified planetary systems. In this work, we investigate the performance of a machine learning classifier on tightlypacked systems containing a rich diversity of planets, from Earths to Jupiters. Using information derived from short numerical simulations about a planet’s early orbital evolution and its relationship with the most massive planets in the system, we train a random forest classifier to predict instability with a > 88 percent accuracy. Our classifier relies on relative planet masses and the standard deviation of eccentricity for much of its predictive power. Most misclassified planets lie along a multi-dimensional boundary between stable and unstable planets, indicating that their early orbital evolution is ambiguous. The major reason for misclassification in this work is timescale: because our classifier uses information from only the first 137 years of simulation data, it is blind to late time interactions that cause or prevent instability. Machine learning methods like those utilized in this work provide powerful tools to complement numerical simulations across a wide range of planetary architectures.

79 ASTRONOMY AND ASTROPHYSICS↗

TOI-2076 and TOI-1807: Two Young, Comoving Planetary Systems within 50 pc Identified by TESS that are Ideal Candidates for Further Follow Up

We report the discovery of two planetary systems around comoving stars: TOI-2076 (TIC 27491137) and TOI-1807 (TIC 180695581). TOI-2076 is a nearby (41.9 pc) multiplanetary system orbiting a young (204 ± 50 Myr), bright (K = 7.115 in TIC v8.1) start. TOI-1807 hosts a single transiting planet and is similarly nearby (42.58 pc), similarly young (180 ± 40 Myr ), and bright. Both targets exhibit significant, periodic variability due to starspots, characteristic of their young ages. Using photometric data collected by TESS we identify three transiting planets around TOI-2076 with radii of R {sub b} = 3.3 ± 0.04 R {sub ⊕}, R {sub c} = 4.4 ± 0.05 R {sub ⊕}, and R {sub d} = 4.1 ± 0.07 R {sub ⊕}. Planet TOI-2076b has a period of P {sub b} = 10.356 days. For both TOI-2076c and d, TESS observed only two transits, separated by a 2 yr interval in which no data were collected, preventing a unique period determination. A range of long periods (<17 days) are consistent with the data. We identify a short-period planet around TOI-1807 with a radius of R {sub b} = 1.8 ± 0.04 R {sub ⊕} and a period of P {sub b} = 0.549 days. Their close proximity, and bright, cool host stars, and young ages make these planets excellent candidates for follow up. TOI-1807b is one of the best-known small (R < 2 R{sub ⊕}) planets for characterization via eclipse spectroscopy and phase curves with JWST. TOI-1807b is the youngest ultra-short-period planet discovered to date, providing valuable constraints on formation timescales of short-period planets. Given the rarity of young planets, particularly in multiple-planet systems, these planets present an unprecedented opportunity to study and compare exoplanet formation, and young planet atmospheres, at a crucial transition age for formation theory.

47 OTHER INSTRUMENTATION↗

Architecture of TOI-561 planetary system

ABSTRACT We present new observations from CHEOPS (CHaracterising ExOPlanet Satellite) and TESS (Transiting Exoplanet Survey Satellite) to clarify the architecture of the planetary system hosted by the old Galactic thick disc star TOI-561. Our global analysis, which also includes previously published photometric and radial velocity data, incontrovertibly proves that TOI-561 is hosting at least four transiting planets with periods of 0.44 d (TOI-561 b), 10.8 d (TOI-561 c), 25.7 d (TOI-561 d), and 77.1 d (TOI-561 e) and a fifth non-transiting candidate, TOI-561f with a period of 433 d. The precise characterization of TOI-561’s orbital architecture is interesting since old and metal-poor thick disc stars are less likely to host ultrashort-period super-Earths like TOI-561 b. The new period of planet -e is consistent with the value obtained using radial velocity alone and is now known to be $77.14399\pm 0.00025$ d, thanks to the new CHEOPS and TESS transits. The new data allowed us to improve its radius ($R_p = 2.517 \pm 0.045\,\mathrm{ R}_{\rm{\oplus }}$ from 5 per cent to 2 per cent precision) and mass ($M_p = 12.4 \pm 1.4\, \mathrm{ M}_{\rm{\oplus }}$) estimates, implying a density of $\rho _p = 0.778 \pm 0.097\, \rho _{\rm{\oplus }}$. Thanks to recent TESS observations and the focused CHEOPS visit of the transit of TOI-561 e, a good candidate for exomoon searches, the planet’s period is finally constrained, allowing us to predict transit times through 2030 with 20-min accuracy. We present an updated version of the internal structure of the four transiting planets. We finally performed a detailed stability analysis, which confirmed the long-term stability of the outer planet TOI-561 f.

Piotto, G. (ORCID:0000000299376387)↗

Dynamical Mass Estimates of the β Pictoris Planetary System through Gaussian Process Stellar Activity Modeling

Nearly 15 yr of radial velocity (RV) monitoring and direct imaging enabled the detection of two giant planets orbiting the young, nearby star β Pictoris. The δ Scuti pulsations of the star, which overwhelm planetary signals, need to be carefully suppressed. In this work, we independently revisit the analysis of the RV data following a different approach than available in the literature to model the activity of the star. We show that a Gaussian process (GP) with a stochastically driven damped harmonic oscillator kernel can model the δ Scuti pulsations. It provides similar results to parametric models but with a simpler framework, using only three hyperparameters. It also enables us to model poorly sampled RV data that were excluded from previous analyses, hence extending the RV baseline by nearly five years. Altogether, the orbit and mass of both planets can be constrained from RV only, which was not possible with the parametric modeling. To characterize the system more accurately, we also perform a joint fit of all available relative astrometry and RV data. Our orbital solutions for β Pic b favor a low eccentricity of 0.029$_{−0.024}^{+0.061}$ and a relatively short period of 21.1$_{−0.8}^{+2.0}$ yr. The orbit of β Pic c is eccentric with 0.206$_{−0.063}^{+0.074}$ with a period of 3.36 ± 0.03 yr. We find model-independent masses of 11.7 ± 1.4 and 8.5 ± 0.5 M Jup for β Pic b and c, respectively, assuming coplanarity. The mass of β Pic b is consistent with the hottest start evolutionary models, at an age of 25 ± 3 Myr. A direct detection of β Pic c would provide a second calibration measurement in a coeval system.

79 ASTRONOMY AND ASTROPHYSICS↗

SOAR TESS Survey. I. Sculpting of TESS Planetary Systems by Stellar Companions

The Transiting Exoplanet Survey Satellite (TESS) is finding transiting planet candidates around bright, nearby stars across the entire sky. The large field of view, however, results in low spatial resolution; therefore, multiple stars contribute to almost every TESS light curve. High angular resolution imaging can detect the previously unknown companions to planetary candidate hosts that dilute the transit depths, lead to host star ambiguity, and, in some cases, are the source of false-positive transit signals. We use speckle imaging on the Southern Astrophysical Research (SOAR) telescope to search for companions to 542 TESS planet candidate hosts in the southern sky. We provide correction factors for the 117 systems with resolved companions due to photometric contamination. The contamination in TESS due to close binaries is similar to that found in surveys of Kepler planet candidates. For the solar-type population, we find a deep deficit of close binary systems with projected stellar separations less than 100 au among planet candidate hosts (44 observed binaries compared to 124 expected based on field binary statistics). The close binary suppression among TESS planet candidate hosts is similar to that seen for the more distant Kepler population. We also find a large surplus of TESS planet candidates in wide binary systems detected in both SOAR and Gaia DR2 (119 observed binaries compared to 77 expected). These wide binaries almost exclusively host giant planets, however, suggesting that orbital migration caused by perturbations from the stellar companion may lead to planet–planet scattering and suppress the population of small planets in wide binaries. Both trends are also apparent in the M dwarf planet candidate hosts.

79 ASTRONOMY AND ASTROPHYSICS↗

Second-generation dust in planetary systems: the case of HD 163296

Observations indicate that large, dust-laden protoplanetary discs are common. Some features, like gaps, rings, and spirals, suggest they may host young planets, which can excite the orbits of nearby leftover planetesimals. Energetic collisions among these bodies can lead to the production of second-generation dust. Grains produced by collisions may have a dynamical behaviour different from that of first-generation, primordial dust out of which planetesimals and planets formed. Here, we aim to study these differences for the HD 163296 system and determine whether dynamical signatures in the mixture of the two dust populations can help separate their contributions. We use three-dimensional (3D) hydrodynamic models to describe the gaseous disc with three, Saturn- to Jupiter-mass, embedded planets. Dust grains, of sizes $1\, \mu \mathrm{m}$–$1\, \mathrm{mm}$, are treated as Lagrangian particles with resolved thermodynamics and mass-loss. Initial disc and planet configurations are derived from observation-based work, which indicates low gas viscosity. The 3D approach also allows us to detect the formation of vortices induced by Rossby waves, where dust becomes concentrated and may contribute to planetesimal formation. We find that the main differences in the dynamical behaviour of first- and second-generation dust occur in the vertical distribution. The two populations have similar distributions around the disc mid-plane, although second-generation dust shows longer residence times close to the radial locations of the planets’ gas gaps. Sedimentation rates of $\mu$m-sized grains are comparable to or lower than the production rates by planetesimals’ collisions, making this population potentially observable. These outcomes can be extended to similar systems harbouring giant planets.

79 ASTRONOMY AND ASTROPHYSICS↗

A search for transit timing variations in the HATS-18 planetary system

ABSTRACT HATS-18 b is a transiting planet with a large mass and a short orbital period, and is one of the best candidates for the detection of orbital decay induced by tidal effects. We present extensive photometry of HATS-18 from which we measure 27 times of mid-transit. Two further transit times were measured from data from the Transiting Exoplanet Survey Satellite (TESS) and three more taken from the literature. The transit timings were fitted with linear and quadratic ephemerides and an upper limit on orbital decay was determined. This corresponds to a lower limit on the modified stellar tidal quality factor of $Q_\star ^{\, \prime } \gt 10^{5.11 \pm 0.04}$. This is at the cusp of constraining the presence of enhanced tidal dissipation due to internal gravity waves. We also refine the measured physical properties of the HATS-18 system, place upper limits on the masses of third bodies, and compare the relative performance of TESS and the 1.54 m Danish Telescope in measuring transit times for this system.

Southworth, John (ORCID:0000000238073198)↗

Orbital Refinement and Stellar Properties for the HD 9446, HD 43691, and HD 179079 Planetary Systems

The Transit Ephemeris Refinement and Monitoring Survey is a project that aims to detect transits of intermediate-long period planets by refining orbital parameters of the known radial velocity planets using additional data from ground-based telescopes, calculating a revised transit ephemeris for the planet, then monitoring the planet host star during the predicted transit window. Here we present the results from three systems that had high probabilities of transiting planets: HD 9446 b and c, HD 43691 b, and HD 179079 b. We provide new radial velocity (RV) measurements that are then used to improve the orbital solution for the known planets. We search the RV data for indications of additional planets in orbit and find that HD 9446 shows a strong linear trend of 4.8σ. Using the newly refined planet orbital solutions, which include a new best-fit solution for the orbital period of HD 9446 c, and an improved transit ephemerides, we found no evidence of transiting planets in the photometry for each system. Transits of HD 9446 b can be ruled out completely and transits HD 9446 c and HD 43691 b can be ruled out for impact parameters up to b = 0.5778 and b = 0.898, respectively, due to gaps in the photometry. A transit of HD 179079 b cannot be ruled out, however, due to the relatively small size of this planet compared to the large star and thus low signal to noise. We determine properties of the three host stars through spectroscopic analysis and find through photometric analysis that HD 9446 exhibits periodic variability.

79 ASTRONOMY AND ASTROPHYSICS↗

Tycho: Realistically Simulating Exoplanets within Stellar Clusters. I. Improving the Monte Carlo Approach

To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses >0.5M {sub ⊙} (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additional bodies (as in the case of HD 285507b). However, it has proven difficult to investigate these effects, as many current numerical solvers for the multi-scale N-body problem are simplified and limited in scope. To remedy this, we aim to create a physically complete computational solution to explore the role of stellar close encounters and interplanetary interactions in producing the observed exoplanet populations for both open cluster stars and field stars. We present a new code, Tycho, which employs a variety of different computational techniques, including multiple N-body integration methods, close encounter handling, modified Monte Carlo scattering experiments, and a variety of empirically informed initial conditions. We discuss the methodology in detail, and its implementation within the AMUSE software framework. Approximately 1% of our systems are promptly disrupted by star-star encounters contributing to the rogue planets occurrence rate. Additionally, we find that close encounters which that perturb long-period planets lead to 38.3% of solar-system-like planetary systems becoming long-term unstable.

79 ASTRONOMY AND ASTROPHYSICS↗

A Planetary Defense System For Responsive Space Missions

As King et al succinctly write, “Collision between the Earth and small objects in the solar system are known to occur regularly, and can result in extinction-level events at their most energetic or more commonly widespread loss of life or economic destruction even for relatively low-energy impacts.” Currently, the JPL Sentry System and the NASA NEO Surveyor mission track, identify, and monitor potential asteroid threats. The JPL Sentry System4 is an automated collision monitoring system which scans the current asteroid catalog for potential impact trajectories with Earth. A 2017 study led by P. Tricarico titled “The near-Earth asteroid population from two decades of observations” analyzes the combined observations of nine leading asteroid surveys. The NASA NEO Surveyor mission seeks to identify 90% of near-Earth objects (NEO) with diameters greater than 140 meters (460 feet) within 10 years. The NEO Surveyor combines data from two infrared imaging detectors to measure NEO size, composition, shape, rotational state, and orbits. Recently, the NASA Double Asteroid Redirection Test (DART) mission successfully altered the orbit of asteroid moon Dimorphos in 2022. This demonstrated that orbit alteration is possible. Thus, collisions between the Earth and asteroids are both likely and demonstrate potential for great harm. There is ongoing work to identify, catalog, and track potential threats. However, when a threat is identified there is not yet a system in place to target the threatening asteroid and defend the planet from potentially existential threats. Building on the DART Mission, we explore a planetary defense system to mitigate asteroid threats.

42 ENGINEERING↗

Giant Planet Scatterings and Collisions: Hydrodynamics, Merger-Ejection Branching Ratio, and Properties of the Remnants

Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple “sticky-sphere” prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and N-body integrations to investigate the fluid effects in close encounters between gas giants and the long-term evolution of closely-packed planetary systems. Focusing on parabolic encounters between Jupiter-like planets with MJ and 2MJ, we find that quick mergers occur when the impact parameter rp (the pericenter separation between the planets) is less than 2RJ, and the merger conserved at least 97% of the initial mass. Strong tidal effects can affect the “binary-planet” orbit when rp is between 2RJ and 4RJ. We quantify these effects using a set of fitting formulae that can be implemented in N-body codes. We run a suite of N-body simulations with and without the formulae for systems of two giant planets initially in unstable, nearly circular and nearly coplanar orbits. The fluid (tidal) effects significantly increase the branching ratio of planetary mergers relative to ejections by doubling the effective collision radius. While the fluid effects do not change the distributions of semi-major axis and eccentricity of each type of remnant planets (mergers vs surviving planets in ejections), the overall orbital properties of planet scattering remnants are strongly affected due to the branching ratio change. We also find that the merger products have broad distributions of spin magnitudes and obliquities.

79 ASTRONOMY AND ASTROPHYSICS↗

A Numerical Method for Determining the Elements of Circumbinary Orbits and Its Application to Circumbinary Planets and the Satellites of Pluto-Charon

Planets and satellites orbiting a binary system exist in the solar system and extrasolar planetary systems. Their orbits can be significantly different from Keplerian orbits, if they are close to the binary and the secondary-to-primary mass ratio is high. A proper description of a circumbinary orbit is in terms of the free eccentricity e {sub free} at the epicyclic frequency κ {sub 0}, forced eccentricity e {sub forced} at the mean motion n {sub 0}, and oscillations at higher frequencies forced by the non-axisymmetric components of the binary’s potential. We show that accurate numerical values for the amplitudes and frequencies of these terms can be extracted from numerical orbit integrations by applying fast Fourier transformation (FFT) to the cylindrical distance between the circumbinary object and the center of mass of the binary as a function of time. We apply this method to three Kepler circumbinary planets and the satellites of Pluto-Charon. For the satellite Styx of Pluto-Charon, the FFT results for κ {sub 0} and e {sub free} differ significantly from the first-order analytic value and the value reported by Showalter & Hamilton, respectively. We show that the deviation in κ {sub 0} is likely due to the effect of the 3:1 mean-motion resonance and discuss the implications of the lower value for e {sub free}.

79 ASTRONOMY AND ASTROPHYSICS↗

Planetary Spin and Obliquity from Mergers

In planetary systems with sufficiently small inter-planet spacing, close encounters can lead to planetary collisions/mergers or ejections. Here, we study the spin property of the merger products of two giant planets in a statistical manner using numerical simulations and analytical modeling. Planetary collisions lead to rapidly rotating objects and a broad range of obliquities. We find that, under typical conditions for two-planet scatterings, the distributions of spin magnitude S and obliquity ${\theta }_{\mathrm{SL}}$ of the merger products have simple analytical forms: f S ∝ S and ${f}_{\cos {\theta }_{\mathrm{SL}}}\propto {(1-{\cos }^{2}{\theta }_{\mathrm{SL}})}^{-1/2}$. Through parameter studies, we determine the regime of validity for the analytical distributions of spin and obliquity. Since planetary mergers are a major outcome of planet–planet scatterings, observational search for the spin/obliquity signatures of exoplanets would provide important constraints on the dynamical history of planetary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

A possible misaligned orbit for the young planet AU Mic c

ABSTRACT The AU Microscopii planetary system is only 24 Myr old, and its geometry may provide clues about the early dynamical history of planetary systems. Here, we present the first measurement of the Rossiter–McLaughlin effect for the warm sub-Neptune AU Mic c, using two transits observed simultaneously with the European Southern Observatory's (ESO's) Very Large Telescope (VLT)/Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO), CHaracterising ExOPlanet Satellite (CHEOPS), and Next-Generation Transit Survey (NGTS). After correcting for flares and for the magnetic activity of the host star, and accounting for transit-timing variations, we find the sky-projected spin–orbit angle of planet c to be in the range $\lambda _{\mathrm{c}}=67.8_{-49.0}^{+31.7}$ degrees (1$\sigma$). We examine the possibility that planet c is misaligned with respect to the orbit of the inner planet b ($\lambda _{\mathrm{b}}=-2.96_{-10.30}^{+10.44}$), and the equatorial plane of the host star, and discuss scenarios that could explain both this and the planet’s high density, including secular interactions with other bodies in the system or a giant impact. We note that a significantly misaligned orbit for planet c is in some degree of tension with the dynamical stability of the system, and with the fact that we see both planets in transit, though these arguments alone do not preclude such an orbit. Further observations would be highly desirable to constrain the spin–orbit angle of planet c more precisely.

Yu, H. (ORCID:0000000209716078)↗

Dynamical Instabilities in Systems of Multiple Short-period Planets Are Likely Driven by Secular Chaos: A Case Study of Kepler-102

We investigated the dynamical stability of high-multiplicity Kepler and K2 planetary systems. Our numerical simulations find instabilities in ∼20% of the cases on a wide range of timescales (up to 5 × 10{sup 9} orbits) and over an unexpectedly wide range of initial dynamical spacings. To identify the triggers of long-term instability in multiplanet systems, we investigated in detail the five-planet Kepler-102 system. Despite having several near-resonant period ratios, we find that mean-motion resonances are unlikely to directly cause instability for plausible planet masses in this system. Instead, we find strong evidence that slow inward transfer of angular momentum deficit (AMD) via secular chaos excites the eccentricity of the innermost planet, Kepler-102 b, eventually leading to planet–planet collisions in ∼80% of Kepler-102 simulations. Kepler-102 b likely needs a mass ≳0.1 M {sub ⊕}, hence a bulk density exceeding about half Earth’s, in order to avoid dynamical instability. To investigate the role of secular chaos in our wider set of simulations, we characterize each planetary system’s AMD evolution with a “spectral fraction” calculated from the power spectrum of short integrations (∼5 × 10{sup 6} orbits). We find that small spectral fractions (≲0.01) are strongly associated with dynamical stability on long timescales (5 × 10{sup 9} orbits) and that the median time to instability decreases with increasing spectral fraction. Our results support the hypothesis that secular chaos is the driver of instabilities in many nonresonant multiplanet systems and also demonstrate that the spectral analysis method is an efficient numerical tool to diagnose long-term (in)stability of multiplanet systems from short simulations.

79 ASTRONOMY AND ASTROPHYSICS↗