-
The TROY project III. Exploring co-orbitals around low-mass stars
Authors:
O. Balsalobre-Ruza,
J. Lillo-Box,
D. Barrado,
A. Correia,
J. P. Faria,
P. Figueira,
A. Leleu,
P. Robutel,
N. Santos,
E. Herrero-Cisneros
Abstract:
Co-orbital objects, also known as trojans, are frequently found in simulations of planetary system formation. In these configurations, a planet shares its orbit with other massive bodies. It is still unclear why there have not been any co-orbitals discovered thus far in exoplanetary systems or even pairs of planets found in such a 1:1 mean motion resonance. Reconciling observations and theory is a…
▽ More
Co-orbital objects, also known as trojans, are frequently found in simulations of planetary system formation. In these configurations, a planet shares its orbit with other massive bodies. It is still unclear why there have not been any co-orbitals discovered thus far in exoplanetary systems or even pairs of planets found in such a 1:1 mean motion resonance. Reconciling observations and theory is an open subject in the field. The main objective of the TROY project is to conduct an exhaustive search for exotrojans using diverse observational techniques. In this work, we analyze the radial velocity time series informed by transits, focusing the search around low-mass stars. We employed the alpha-test method on confirmed planets searching for shifts between spectral and photometric mid-transit times. This technique is sensitive to mass imbalances within the planetary orbit, allowing us to identify non-negligible co-orbital masses. Among the 95 transiting planets examined, we find one robust exotrojan candidate with a significant 3-sigma detection. Additionally, 25 exoplanets show compatibility with the presence of exotrojan companions at a 1-sigma level, requiring further observations to better constrain their presence. For two of those weak candidates, we find dimmings in their light curves within the predicted Lagrangian region. We established upper limits on the co-orbital masses for either the candidates and null detections. Our analysis reveals that current high-resolution spectrographs effectively rule out co-orbitals more massive than Saturn around low-mass stars. This work points out to dozens of targets that have the potential to better constraint their exotrojan upper mass limit with dedicated radial velocity observations. We also explored the potential of observing the secondary eclipses of the confirmed exoplanets to enhance the exotrojan search.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Trojan Exoplanets
Authors:
Philippe Robutel,
Adrien Leleu
Abstract:
Co-orbital exoplanets are a by-product of the models of formation of planetary systems. However, none have been detected in nature thus far. Although challenging, the observation of co-orbital exoplanets would provide valuable information on the formation of planetary systems as well as on the interactions between planets and their host star. After a brief review of the stability and formation iss…
▽ More
Co-orbital exoplanets are a by-product of the models of formation of planetary systems. However, none have been detected in nature thus far. Although challenging, the observation of co-orbital exoplanets would provide valuable information on the formation of planetary systems as well as on the interactions between planets and their host star. After a brief review of the stability and formation issues of co-orbital systems, some observational methods dedicated to their detection are presented.
△ Less
Submitted 23 February, 2024;
originally announced February 2024.
-
Long-term dust dynamics in Didymos and Dimorphos system: production, stability, and transport
Authors:
Gustavo Madeira,
Sebastien Charnoz,
Nicolas Rambaux,
Philippe Robutel
Abstract:
Target of NASA's DART mission, the system of Didymos and Dimorphos will once again be visited by a space mission -- ESA's Hera mission, scheduled to be launch in 2024. Hera will arrive in the system approximately 4 years after the DART impact, a long period compared to Dimorphos' orbital period (about 12 hours). It is therefore imperative to understand the dynamics of material in this environment…
▽ More
Target of NASA's DART mission, the system of Didymos and Dimorphos will once again be visited by a space mission -- ESA's Hera mission, scheduled to be launch in 2024. Hera will arrive in the system approximately 4 years after the DART impact, a long period compared to Dimorphos' orbital period (about 12 hours). It is therefore imperative to understand the dynamics of material in this environment on a long timescale. Here, we explore the long-term dynamics of the binary system (65038) Didymos, in the context of the perturbed, planar, circular and restricted 3-body problem. We design an analytical description for a symmetrical top-shaped object, the shape assumed for the Didymos, while the Dimorphos is considered an ellipsoid. In the absence of external effects, we identify seven stable equatorial regions where particles persist for more than a decade. However, in the presence of the solar radiation effect, the lifetime of small particles (<mm) is in the order of days, being unlikely that Hera spacecraft will encounter clusters of millimetre and sub-millimetre particles in stable equatorial orbits. Nonetheless, large objects may reside in the region for some years, particularly in quasi-satellite orbits, the most stable orbits in the system. Additionally, interplanetary dust impacts onto Didymos populate the region, extending up to a distance of approximately 1500 meters from the primary center, with young dust. These impacts are responsible for a transfer of dust mainly from Didymos to Dimorphos. If the interplanetary dust impacts generate metric-sized boulders, they may persist in the system for years, in first sort orbits around Didymos.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Dynamics of co-orbital exoplanets in a first order resonance chain with tidal dissipation
Authors:
Jérémy Couturier,
Philippe Robutel,
Alexandre C. M. Correia
Abstract:
Co-orbital planets (in a $1:1$ mean motion resonance) can be formed within a Laplace resonance chain. Here, we develop a secular model to study the dynamics of the resonance chain $p:p:p+1$, where the co-orbital pair is in a first-order mean motion resonance with the outermost third planet. Our model takes into account tidal dissipation through the use of a Hamiltonian version of the constant time…
▽ More
Co-orbital planets (in a $1:1$ mean motion resonance) can be formed within a Laplace resonance chain. Here, we develop a secular model to study the dynamics of the resonance chain $p:p:p+1$, where the co-orbital pair is in a first-order mean motion resonance with the outermost third planet. Our model takes into account tidal dissipation through the use of a Hamiltonian version of the constant time-lag model, which extends the Hamiltonian formalism of the point-mass case. We show the existence of several families of equilibria, and how these equilibria extend to the complete system. In one family, which we call the main branch, a secular resonance between the libration frequency of the co-orbitals and the precession frequency of the pericentres has unexpected dynamical consequences when tidal dissipation is added. We report the existence of two distinct mechanisms that make co-orbital planets much more stable within the $p:p:p+1$ resonance chain rather than outside it. The first one is due to negative real parts of the eigenvalues of the linearised system with tides, in the region of the secular resonance mentioned above. The second one comes from non-linear contributions of the vector field and it is due to eccentricity damping. These two stabilising mechanisms increase the chances of a still-to-come detection of exoplanets in the co-orbital configuration.
△ Less
Submitted 17 April, 2022;
originally announced April 2022.
-
An analytical model for tidal evolution in co-orbital systems $\,$I. Application to exoplanets
Authors:
Jérémy Couturier,
Philippe Robutel,
Alexandre C. M. Correia
Abstract:
Close-in co-orbital planets (in a 1:1 mean motion resonance) can experience strong tidal interactions with the central star. Here, we develop an analytical model adapted to the study of the tidal evolution of those systems. We use a Hamiltonian version of the constant time-lag tidal model, which extends the Hamiltonian formalism developed for the point-mass case. We show that co-orbital systems un…
▽ More
Close-in co-orbital planets (in a 1:1 mean motion resonance) can experience strong tidal interactions with the central star. Here, we develop an analytical model adapted to the study of the tidal evolution of those systems. We use a Hamiltonian version of the constant time-lag tidal model, which extends the Hamiltonian formalism developed for the point-mass case. We show that co-orbital systems undergoing tidal dissipation either favour the Lagrange or the anti-Lagrange configurations, depending on the system parameters. However, for all range of parameters and initial conditions, both configurations become unstable, although the timescale for the destruction of the system can be larger than the lifetime of the star. We provide an easy-to-use criterion to determine if an already known close-in exoplanet may have an undetected co-orbital companion.
△ Less
Submitted 18 June, 2021;
originally announced June 2021.
-
Six transiting planets and a chain of Laplace resonances in TOI-178
Authors:
A. Leleu,
Y. Alibert,
N. C. Hara,
M. J. Hooton,
T. G. Wilson,
P. Robutel,
J. -B. Delisle,
J. Laskar,
S. Hoyer,
C. Lovis,
E. M. Bryant,
E. Ducrot,
J. Cabrera,
L. Delrez,
J. S. Acton,
V. Adibekyan,
R. Allart,
C. Allende Prieto,
R. Alonso,
D. Alves,
D. R. Anderson,
D. Angerhausen,
G. Anglada Escudé,
J. Asquier,
D. Barrado
, et al. (130 additional authors not shown)
Abstract:
Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this cont…
▽ More
Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth's density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes.
△ Less
Submitted 22 January, 2021;
originally announced January 2021.
-
Planetary system LHS 1140 revisited with ESPRESSO and TESS
Authors:
J. Lillo-Box,
P. Figueira,
A. Leleu,
L. Acuña,
J. P. Faria,
N. Hara,
N. C. Santos,
A. C. M. Correia,
P. Robutel,
M. Deleuil,
D. Barrado,
S. Sousa,
X. Bonfils,
O. Mousis,
J. M. Almenara,
N. Astudillo-Defru,
E. Marcq,
S. Udry,
C. Lovis,
F. Pepe
Abstract:
LHS 1140 is an M dwarf known to host two known transiting planets at orbital periods of 3.77 and 24.7 days. The external planet (LHS 1140 b) is a rocky super-Earth that is located in the middle of the habitable zone of this low-mass star, placing this system at the forefront of the habitable exoplanet exploration. We further characterize this system by improving the physical and orbital properties…
▽ More
LHS 1140 is an M dwarf known to host two known transiting planets at orbital periods of 3.77 and 24.7 days. The external planet (LHS 1140 b) is a rocky super-Earth that is located in the middle of the habitable zone of this low-mass star, placing this system at the forefront of the habitable exoplanet exploration. We further characterize this system by improving the physical and orbital properties and search for additional planetary-mass components in the system, also exploring the possibility of co-orbitals. We collected 113 new radial velocity observations with ESPRESSO over a 1.5-year time span with an average photon-noise precision of 1.07 m/s. We determine new masses with a precision of 6% for LHS 1140 b ($6.48 \pm 0.46~M_{\oplus}$) and 9% for LHS 1140 c ($m_c=1.78 \pm 0.17~M_{\oplus}$), reducing by half the previously published uncertainties. Although both planets have Earth-like bulk compositions, the internal structure analysis suggests that LHS 1140 b might be iron-enriched. In both cases, the water content is compatible to a maximum fraction of 10-12% in mass, which is equivalent to a deep ocean layer of $779 \pm 650$ km for the habitable-zone planet LHS 1140 b. Our results also provide evidence for a new planet candidate in the system ($m_d= 4.8\pm1.1~M_{\oplus}$) on a ~78.9-day orbital period, which is detected through three independent methods. The analysis also allows us to discard other planets above 0.5 $M_{\oplus}$ for periods shorter than 10 days and above 2 $M_{\oplus}$ for periods up to one year. Finally, our analysis discards co-orbital planets of LHS 1140 b down to 1 $M_{\oplus}$. Indications for a possible co-orbital signal in LHS 1140 c are detected in both radial velocity and photometric data, however. The new characterization of the system make it a key target for atmospheric studies of rocky worlds at different stellar irradiations
△ Less
Submitted 14 October, 2020;
originally announced October 2020.
-
Co-orbital exoplanets from close period candidates: The TOI-178 case
Authors:
Adrien Leleu,
Jorge Lillo-Box,
Marko Sestovic,
Philippe Robutel,
Alexandre Correia,
Nathan Hara,
Daniel Angerhausen,
Simon Grimm,
Jean Schneider
Abstract:
Despite the existence of co-orbital bodies in the solar system, and the prediction of the formation of co-orbital planets by planetary system formation models, no co-orbital exoplanets (also called trojans) have been detected thus far. Here we study the signature of co-orbital exoplanets in transit surveys when two planet candidates in the system orbit the star with similar periods. Such pair of c…
▽ More
Despite the existence of co-orbital bodies in the solar system, and the prediction of the formation of co-orbital planets by planetary system formation models, no co-orbital exoplanets (also called trojans) have been detected thus far. Here we study the signature of co-orbital exoplanets in transit surveys when two planet candidates in the system orbit the star with similar periods. Such pair of candidates could be discarded as false positives because they are not Hill-stable. However, horseshoe or long libration period tadpole co-orbital configurations can explain such period similarity. This degeneracy can be solved by considering the Transit Timing Variations (TTVs) of each planet. We then focus on the three planet candidates system TOI-178: the two outer candidates of that system have similar orbital period and had an angular separation near $π/3$ during the TESS observation of sector 2. Based on the announced orbits, the long-term stability of the system requires the two close-period planets to be co-orbitals. Our independent detrending and transit search recover and slightly favour the three orbits close to a 3:2:2 resonant chain found by the TESS pipeline, although we cannot exclude an alias that would put the system close to a 4:3:2 configuration. We then analyse in more detail the co-orbital scenario. We show that despite the influence of an inner planet just outside the 2:3 mean-motion resonance, this potential co-orbital system can be stable on the Giga-year time-scale for a variety of planetary masses, either on a trojan or a horseshoe orbit. We predict that large TTVs should arise in such configuration with a period of several hundred days. We then show how the mass of each planet can be retrieved from these TTVs.
△ Less
Submitted 22 January, 2019;
originally announced January 2019.
-
La stabilité des lunes de Saturne, Janus et Épiméthée: de l'observation astronomique à la théorie KAM
Authors:
Alexandre Pousse,
Laurent Niederman,
Philippe Robutel
Abstract:
Popular science article associated with the work `On the co-orbital motion in the three-body problem: existence of quasi-periodic horseshoe-shaped orbits" (arXiv:1806.07262) from the same authors. Janus and Epimetheus are two moons of Saturn which exhibit a really peculiar dynamics. As they orbit on circular trajectories whose radii are only 50 km apart (less than their respective diameters), ever…
▽ More
Popular science article associated with the work `On the co-orbital motion in the three-body problem: existence of quasi-periodic horseshoe-shaped orbits" (arXiv:1806.07262) from the same authors. Janus and Epimetheus are two moons of Saturn which exhibit a really peculiar dynamics. As they orbit on circular trajectories whose radii are only 50 km apart (less than their respective diameters), every four (terrestrial) years the bodies are getting closer and their mutual gravitational influence leads to a swapping of the orbits: the outer moon becoming the inner one and vice-versa. In this article, we describe how, from this specific astronomical problem to the KAM theory, we came to prove the existence of perpetually stable trajectories associated with the Janus and Epimetheus orbits.
△ Less
Submitted 25 July, 2018;
originally announced July 2018.
-
The TROY project: II. Multi-technique constraints on exotrojans in nine planetary systems
Authors:
J. Lillo-Box,
A. Leleu,
H. Parviainen,
P. Figueira,
M. Mallonn,
A. C. M. Correia,
N. C. Santos,
P. Robutel,
M. Lendl,
H. M. J. Boffin,
J. P. Faria,
D. Barrado,
J. Neal
Abstract:
Co-orbital bodies are the byproduct of planet formation and evolution, as we know from the Solar System. Although planet-size co-orbitals do not exists in our planetary system, dynamical studies show that they can remain stable for long periods of time in the gravitational well of massive planets. Should they exist, their detection is feasible with the current instrumentation. In this paper, we pr…
▽ More
Co-orbital bodies are the byproduct of planet formation and evolution, as we know from the Solar System. Although planet-size co-orbitals do not exists in our planetary system, dynamical studies show that they can remain stable for long periods of time in the gravitational well of massive planets. Should they exist, their detection is feasible with the current instrumentation. In this paper, we present new ground-based observations searching for these bodies co-orbiting with nine close-in (P<5 days) planets, using different observing techniques. The combination of all of them allows us to restrict the parameter space of any possible trojan in the system. We use multi-technique observations (radial velocity, precision photometry and transit timing variations), both newly acquired in the context of the TROY project and publicly available, to constrain the presence of planet-size trojans in the Lagrangian points of nine known exoplanets. We find no clear evidence of trojans in these nine systems through any of the techniques used down to the precision of the observations. However, this allows us to constrain the presence of any potential trojan in the system, specially in the trojan mass/radius versus libration amplitude plane. In particular, we can set upper mass limits in the super-Earth mass regime for six of the studied systems.
△ Less
Submitted 2 July, 2018;
originally announced July 2018.
-
On the co-orbital motion in the three-body problem: existence of quasi-periodic horseshoe-shaped orbits
Authors:
Laurent Niederman,
Alexandre Pousse,
Philippe Robutel
Abstract:
Janus and Epimetheus are two moons of Saturn with very peculiar motions. As they orbit around Saturn on quasi-coplanar and quasi-circular trajectories whose radii are only 50 km apart (less than their respective diameters), every four (terrestrial) years the bodies approach each other and their mutual gravitational influence lead to a swapping of the orbits: the outer moon becomes the inner one an…
▽ More
Janus and Epimetheus are two moons of Saturn with very peculiar motions. As they orbit around Saturn on quasi-coplanar and quasi-circular trajectories whose radii are only 50 km apart (less than their respective diameters), every four (terrestrial) years the bodies approach each other and their mutual gravitational influence lead to a swapping of the orbits: the outer moon becomes the inner one and vice-versa. This behavior generates horseshoe-shaped trajectories depicted in an appropriate rotating frame. In spite of analytical theories and numerical investigations developed to describe their long-term dynamics, so far very few rigorous long-time stability results on the "horseshoe motion" have been obtained even in the restricted three-body problem. Adapting the idea of Arnol'd (1963) to a resonant case (the co-orbital motion is associated with trajectories in 1:1 mean motion resonance), we provide a rigorous proof of existence of 2-dimensional elliptic invariant tori on which the trajectories are similar to those followed by Janus and Epimetheus. For this purpose, we apply KAM theory to the planar three-body problem.
△ Less
Submitted 14 February, 2020; v1 submitted 18 June, 2018;
originally announced June 2018.
-
Towards completing Planetary Systems: The role of minor bodies on life growth and survival
Authors:
Jorge Lillo-Box,
David Kipping,
Isabel Rebollido,
Pedro Figueira,
Adrien Leleu,
Alexandre Correia,
Philippe Robutel,
Nuno C. Santos,
David Barrado,
Benjamín Montesinos,
Tjarda Boekholt
Abstract:
The search for extrasolar planets in the past decades has shown that planets abound in the Solar neighborhood. While we are still missing an Earth twin, the forthcoming space missions and ground-based instrumentation are already driven to achieve this goal. But, in order to fully understand the conditions for life appearing in the Solar System, we still miss some pieces of the planetary system jig…
▽ More
The search for extrasolar planets in the past decades has shown that planets abound in the Solar neighborhood. While we are still missing an Earth twin, the forthcoming space missions and ground-based instrumentation are already driven to achieve this goal. But, in order to fully understand the conditions for life appearing in the Solar System, we still miss some pieces of the planetary system jigsaw puzzle, namely a deeper understanding of the minor bodies. Trojans, moons, and comets are tracers of the formation and evolution processes of planetary systems. These missing pieces are also critical to understand the emergence and evolution of life over millions of years. With the large crop of planetary systems discovered so far and yet to be detected with the forthcoming missions, the hunt for minor bodies in extrasolar systems is a natural continuation of our search for real Solar System- and, in particular, Earth- analogs. This white paper is focused on detection of these minor components and their relevance in the emergence, evolution and survival of life.
△ Less
Submitted 11 March, 2018;
originally announced March 2018.
-
On the coplanar eccentric non restricted co-orbital dynamics
Authors:
A. Leleu,
P. Robutel,
A. C. M. Correia
Abstract:
We study the phase space of eccentric coplanar co-orbitals in the non-restricted case. Departing from the quasi-circular case, we describe the evolution of the phase space as the eccentricities increase. We find that over a given value of the eccentricity, around $0.5$ for equal mass co-orbitals, important topological changes occur in the phase space. These changes lead to the emergence of new co-…
▽ More
We study the phase space of eccentric coplanar co-orbitals in the non-restricted case. Departing from the quasi-circular case, we describe the evolution of the phase space as the eccentricities increase. We find that over a given value of the eccentricity, around $0.5$ for equal mass co-orbitals, important topological changes occur in the phase space. These changes lead to the emergence of new co-orbital configurations and open a continuous path between the previously distinct trojan domains near the $L_4$ and $L_5$ eccentric Lagrangian equilibria. These topological changes are shown to be linked with the reconnection of families of quasi-periodic orbits of non-maximal dimension.
△ Less
Submitted 20 November, 2017;
originally announced November 2017.
-
Secular Dynamics of Multiplanetary Circumbinary Systems
Authors:
Eduardo Andrade-Ines,
Philippe Robutel
Abstract:
We present an analytical formalism to study the secular dynamics of a system consisting of N-2 planets orbiting a binary star in outer orbits. We introduce a canonical coordinate system and expand the disturbing function in terms of canonical elliptic elements, combining both Legendre polynomials and Laplace coefficients, to obtain a general formalism for the secular description of this type of co…
▽ More
We present an analytical formalism to study the secular dynamics of a system consisting of N-2 planets orbiting a binary star in outer orbits. We introduce a canonical coordinate system and expand the disturbing function in terms of canonical elliptic elements, combining both Legendre polynomials and Laplace coefficients, to obtain a general formalism for the secular description of this type of configuration. With a quadractic approximation of the development, we present a simplified analytical solution for the planetary orbits for both the single planet and the two-planet cases. From the two-planet model, we show that the inner planet accelerates the precession rate of the binary pericenter, which, in turn, may enter in resonance with the secular frequency of the outer planet, characterizing a secular resonance. We calculate an analytical expression for the approximate location of this resonance and apply it to known circumbinary systems, where we show that it can occur at relatively close orbits, for example at 2.4 au for the Kepler-38 system. With a more refined model, we analyse the dynamics of this secular resonance and we show that a bifurcation of the corresponding fixed points can affect the long term evolution and stability of planetary systems. By comparing our results with complete integrations of the exact equations of motion, we verified the accuracy of our analytical model.
△ Less
Submitted 6 November, 2017;
originally announced November 2017.
-
The TROY project: Searching for co-orbital bodies to known planets. I. Project goals and first results from archival radial velocity
Authors:
J. Lillo-Box,
D. Barrado,
P. Figueira,
A. Leleu,
N. C. Santos,
A. C. M. Correia,
P. Robutel,
J. P. Faria
Abstract:
The detection of Earth-like planets, exocomets or Kuiper belts show that the different components found in the solar system should also be present in other planetary systems. Trojans are one of these components and can be considered fossils of the first stages in the life of planetary systems. Their detection in extrasolar systems would open a new scientific window to investigate formation and mig…
▽ More
The detection of Earth-like planets, exocomets or Kuiper belts show that the different components found in the solar system should also be present in other planetary systems. Trojans are one of these components and can be considered fossils of the first stages in the life of planetary systems. Their detection in extrasolar systems would open a new scientific window to investigate formation and migration processes. In this context, the main goal of the TROY project is to detect exotrojans for the first time and to measure their occurrence rate (eta-Trojan). In this first paper, we describe the goals and methodology of the project. Additionally, we used archival radial velocity data of 46 planetary systems to place upper limits on the mass of possible trojans and investigate the presence of co-orbital planets down to several tens of Earth masses. We used archival radial velocity data of 46 close-in (P<5 days) transiting planets (without detected companions) with information from high-precision radial velocity instruments. We took advantage of the time of mid-transit and secondary eclipses (when available) to constrain the possible presence of additional objects co-orbiting the star along with the planet. This, together with a good phase coverage, breaks the degeneracy between a trojan planet signature and signals coming from additional planets or underestimated eccentricity. We identify nine systems for which the archival data provide 1-sigma evidence for a mass imbalance between L4 and L5. Two of these systems provide 2-sigma detection, but no significant detection is found among our sample. We also report upper limits to the masses at L4/L5 in all studied systems and discuss the results in the context of previous findings.
△ Less
Submitted 3 October, 2017;
originally announced October 2017.
-
Spin dynamics of close-in planets exhibiting large TTVs
Authors:
J. -B. Delisle,
A. C. M. Correia,
A. Leleu,
P. Robutel
Abstract:
We study the spin evolution of close-in planets in compact multi-planetary systems. The rotation period of these planets is often assumed to be synchronous with the orbital period due to tidal dissipation. Here we show that planet-planet perturbations can drive the spin of these planets into non-synchronous or even chaotic states. In particular, we show that the transit timing variation (TTV) is a…
▽ More
We study the spin evolution of close-in planets in compact multi-planetary systems. The rotation period of these planets is often assumed to be synchronous with the orbital period due to tidal dissipation. Here we show that planet-planet perturbations can drive the spin of these planets into non-synchronous or even chaotic states. In particular, we show that the transit timing variation (TTV) is a very good probe to study the spin dynamics, since both are dominated by the perturbations of the mean longitude of the planet. We apply our model to KOI-227b and Kepler-88b, which are both observed undergoing strong TTVs. We also perform numerical simulations of the spin evolution of these two planets. We show that for KOI-227b non-synchronous rotation is possible, while for Kepler-88b the rotation can be chaotic.
△ Less
Submitted 12 May, 2017;
originally announced May 2017.
-
Detection of co-orbital planets by combining transit and radial-velocity measurements
Authors:
Adrien Leleu,
Philippe Robutel,
Alexandre C M Correia,
Jorge Lillo-Box
Abstract:
Co-orbital planets have not yet been discovered, although they constitute a frequent by-product of planetary formation and evolution models. This lack may be due to observational biases, since the main detection methods are unable to spot co-orbital companions when they are small or near the Lagrangian equilibrium points. However, for a system with one known transiting planet (with mass $m_1$), we…
▽ More
Co-orbital planets have not yet been discovered, although they constitute a frequent by-product of planetary formation and evolution models. This lack may be due to observational biases, since the main detection methods are unable to spot co-orbital companions when they are small or near the Lagrangian equilibrium points. However, for a system with one known transiting planet (with mass $m_1$), we can detect a co-orbital companion (with mass $m_2$) by combining the time of mid-transit with the radial-velocity data of the star. Here, we propose a simple method that allows the detection of co-orbital companions, valid for eccentric orbits, that relies on a single parameter $α$, which is proportional to the mass ratio $m_2/m_1$. Therefore, when $α$ is statistically different from zero, we have a strong candidate to harbour a co-orbital companion. We also discuss the relevance of false positives generated by different planetary configurations.
△ Less
Submitted 28 February, 2017;
originally announced February 2017.
-
On the co-orbital motion in the planar restricted three-body problem: the quasi-satellite motion revisited
Authors:
Alexandre Pousse,
Philippe Robutel,
Alain Vienne
Abstract:
In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a non zero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet it descri…
▽ More
In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a non zero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet it describes the same trajectory as the one of a retrograde satellite even though the planet acts as a perturbator. In the last few years, the discoveries of asteroids in this type of motion made the term "quasi-satellite" more and more present in the literature. However, some authors rather use the term "retrograde satellite" when referring to this kind of motion in the studies of the restricted problem in the rotating frame. In this paper we intend to clarify the terminology to use, in order to bridge the gap between the perturbative co-orbital point of view and the more general approach in the rotating frame. Through a numerical exploration of the co-orbital phase space, we describe the quasi-satellite domain and highlight that it is not reachable by low eccentricities by averaging process. We will show that the quasi-satellite domain is effectively included in the domain of the retrograde satellites and neatly defined in terms of frequencies. Eventually, we highlight a remarkable high eccentric quasi-satellite orbit corresponding to a frozen ellipse in the heliocentric frame. We extend this result to the eccentric case (planet on an eccentric motion) and show that two families of frozen ellipses originate from this remarkable orbit.
△ Less
Submitted 5 November, 2018; v1 submitted 21 March, 2016;
originally announced March 2016.
-
Secular Dynamics of S-type Planetary Orbits in Binary Star Systems: Applicability Domains of First- and Second-Order Theories
Authors:
Eduardo Andrade-Ines,
Cristian Beaugé,
Tatiana Michtchenko,
Philippe Robutel
Abstract:
We analyse the secular dynamics of planets on S-type coplanar orbits in tight binary systems, based on first- and second-order analytical models, and compare their predictions with full N-body simulations. The perturbation parameter adopted for the development of these models depends on the masses of the stars and on the semimajor axis ratio between the planet and the binary.
We show that each m…
▽ More
We analyse the secular dynamics of planets on S-type coplanar orbits in tight binary systems, based on first- and second-order analytical models, and compare their predictions with full N-body simulations. The perturbation parameter adopted for the development of these models depends on the masses of the stars and on the semimajor axis ratio between the planet and the binary.
We show that each model has both advantages and limitations. While the first-order analytical model is algebraically simple and easy to implement, it is only applicable in regions of the parameter space where the perturbations are sufficiently small. The second-order model, although more complex, has a larger range of validity and must be taken into account for dynamical studies of some real exoplanetary systems such as $γ$-Cephei and HD 41004A. However, in some extreme cases, neither of these analytical models yields quantitatively correct results, requiring either higher-order theories or direct numerical simulations.
Finally, we determine the limits of applicability of each analytical model in the parameter space of the system, giving an important visual aid to decode which secular theory should be adopted for any given planetary system in a close binary.
△ Less
Submitted 11 December, 2015;
originally announced December 2015.
-
On the rotation of co-orbital bodies in eccentric orbits
Authors:
A. Leleu,
P. Robutel,
A. C. M. Correia
Abstract:
We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin-orbit resonances, one de…
▽ More
We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin-orbit resonances, one depending on the eccentricity, one depending on the orbital libration frequency, and another depending on the pericenter's dynamics. We can estimate the width and the location of the different resonant islands in the phase space, predicting which are the more likely to capture the spin of the rotating body. In some regions of the phase space the resonant islands may overlap, giving rise to chaotic rotation.
△ Less
Submitted 25 March, 2016; v1 submitted 30 October, 2015;
originally announced October 2015.
-
Detectability of quasi-circular co-orbital planets. Application to the radial velocity technique
Authors:
Adrien Leleu,
Philippe Robutel,
Alexandre C. M. Correia
Abstract:
Several celestial bodies in co-orbital configurations exist in the solar system. However, co-orbital exoplanets have not yet been discovered. This lack may result from a degeneracy between the signal induced by co-orbital planets and other orbital configurations. Here we determine a criterion for the detectability of quasi-circular co-orbital planets and develop a demodulation method to bring out…
▽ More
Several celestial bodies in co-orbital configurations exist in the solar system. However, co-orbital exoplanets have not yet been discovered. This lack may result from a degeneracy between the signal induced by co-orbital planets and other orbital configurations. Here we determine a criterion for the detectability of quasi-circular co-orbital planets and develop a demodulation method to bring out their signature from the observational data. We show that the precision required to identify a pair of co-orbital planets depends only on the libration amplitude and on the planet's mass ratio. We apply our method to synthetic radial velocity data, and show that for tadpole orbits we are able to determine the inclination of the system to the line of sight. Our method is also valid for planets detected through the transit and astrometry techniques.
△ Less
Submitted 25 March, 2016; v1 submitted 8 September, 2015;
originally announced September 2015.
-
Spin-orbit coupling and chaotic rotation for circumbinary bodies. Application to the small satellites of the Pluto-Charon system
Authors:
Alexandre C. M. Correia,
Adrien Leleu,
Nicolas Rambaux,
Philippe Robutel
Abstract:
We investigate the resonant rotation of circumbinary bodies in planar quasi-circular orbits. Denoting $n_b$ and $n$ the orbital mean motion of the inner binary and of the circumbinary body, respectively, we show that spin-orbit resonances exist at the frequencies $n\pm kν/2$, where $ν= n_b - n$, and $k$ is an integer. Moreover, when the libration at natural frequency has the same magnitude as $ν$,…
▽ More
We investigate the resonant rotation of circumbinary bodies in planar quasi-circular orbits. Denoting $n_b$ and $n$ the orbital mean motion of the inner binary and of the circumbinary body, respectively, we show that spin-orbit resonances exist at the frequencies $n\pm kν/2$, where $ν= n_b - n$, and $k$ is an integer. Moreover, when the libration at natural frequency has the same magnitude as $ν$, the resonances overlap and the rotation becomes chaotic. We apply these results to the small satellites in the Pluto-Charon system, and conclude that their rotations are likely chaotic. However, the rotation can also be stable and not synchronous for small axial asymmetries.
△ Less
Submitted 20 August, 2015; v1 submitted 22 June, 2015;
originally announced June 2015.
-
Rigorous treatment of the averaging process for co-orbital motions in the planetary problem
Authors:
Philippe Robutel,
Laurent Niederman,
Alexandre Pousse
Abstract:
We develop a rigorous analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. By constructing a complex domain of holomorphy for the planetary Hamilto-nian, we estimate the size of the transformation that maps this Hamil-tonian to its first order averaged over one of the fast angles. After having derived an integrable approximation of the average…
▽ More
We develop a rigorous analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. By constructing a complex domain of holomorphy for the planetary Hamilto-nian, we estimate the size of the transformation that maps this Hamil-tonian to its first order averaged over one of the fast angles. After having derived an integrable approximation of the averaged problem, we bound the distance between this integrable approximation and the averaged Hamiltonian. This finally allows to prove rigorous theorems on the behavior of co-orbital motions over a finite but large timescale.
△ Less
Submitted 1 September, 2015; v1 submitted 9 June, 2015;
originally announced June 2015.
-
The family of Quasi-satellite periodic orbits in the circular co-planar RTBP
Authors:
Alexandre Pousse,
Philippe Robutel,
Alain Vienne
Abstract:
In the circular case of the coplanar Restricted Three-body Problem, we studied how the family of quasi-satellite (QS) periodic orbits allows to define an associated libration center. Using the averaged problem, we highlighted a validity limit of this one: for QS orbits with low eccentricities, the averaged problem does not correspond to the real problem. We do the same procedure to L 3 , L 4 and L…
▽ More
In the circular case of the coplanar Restricted Three-body Problem, we studied how the family of quasi-satellite (QS) periodic orbits allows to define an associated libration center. Using the averaged problem, we highlighted a validity limit of this one: for QS orbits with low eccentricities, the averaged problem does not correspond to the real problem. We do the same procedure to L 3 , L 4 and L 5 emerging periodic orbits families and remarked that for very high eccentricities F L4 and F L5 merge with F L3 which bifurcates to a stable family.
△ Less
Submitted 11 December, 2014;
originally announced December 2014.
-
Spin-orbit coupling and chaotic rotation for eccentric coorbital bodies
Authors:
Adrien Leleu,
Philippe Robutel,
A. C. M. Correia
Abstract:
The presence of a co-orbital companion induces the splitting of the well known Keplerian spin-orbit resonances. It leads to chaotic rotation when those resonances overlap.
The presence of a co-orbital companion induces the splitting of the well known Keplerian spin-orbit resonances. It leads to chaotic rotation when those resonances overlap.
△ Less
Submitted 14 November, 2014;
originally announced November 2014.
-
Spin-orbit resonances and rotation of coorbital bodies in quasi-circular orbits
Authors:
Philippe Robutel,
A. C. M. Correia,
Adrien Leleu
Abstract:
The rotation of asymmetric bodies in eccentric Keplerian orbits can be chaotic when there is some overlap of spin-orbit resonances. Here we show that the rotation of two coorbital bodies (two planets orbiting a star or two satellites of a planet) can also be chaotic even for quasi-circular orbits around the central body. When dissipation is present, the rotation period of a body on a nearly circul…
▽ More
The rotation of asymmetric bodies in eccentric Keplerian orbits can be chaotic when there is some overlap of spin-orbit resonances. Here we show that the rotation of two coorbital bodies (two planets orbiting a star or two satellites of a planet) can also be chaotic even for quasi-circular orbits around the central body. When dissipation is present, the rotation period of a body on a nearly circular orbit is believed to always end synchronous with the orbital period. Here we demonstrate that for coorbital bodies in quasi-circular orbits, stable non-synchronous rotation is possible for a wide range of mass ratios and body shapes. We further show that the rotation becomes chaotic when the natural rotational libration frequency, due to the axial asymmetry, is of the same order of magnitude as the orbital libration frequency.
△ Less
Submitted 12 October, 2014;
originally announced October 2014.
-
Spin-orbit coupling and chaotic rotation for coorbital bodies in quasi-circular orbits
Authors:
Alexandre C. M. Correia,
Philippe Robutel
Abstract:
Coorbital bodies are observed around the Sun sharing their orbits with the planets, but also in some pairs of satellites around Saturn. The existence of coorbital planets around other stars has also been proposed. For close-in planets and satellites, the rotation slowly evolves due to dissipative tidal effects until some kind of equilibrium is reached. When the orbits are nearly circular, the rota…
▽ More
Coorbital bodies are observed around the Sun sharing their orbits with the planets, but also in some pairs of satellites around Saturn. The existence of coorbital planets around other stars has also been proposed. For close-in planets and satellites, the rotation slowly evolves due to dissipative tidal effects until some kind of equilibrium is reached. When the orbits are nearly circular, the rotation period is believed to always end synchronous with the orbital period. Here we demonstrate that for coorbital bodies in quasi-circular orbits, stable non-synchronous rotation is possible for a wide range of mass ratios and body shapes. We show the existence of an entirely new family of spin-orbit resonances at the frequencies $n\pm kν/2$, where $n$ is the orbital mean motion, $ν$ the orbital libration frequency, and $k$ an integer. In addition, when the natural rotational libration frequency due to the axial asymmetry, $σ$, has the same magnitude as $ν$, the rotation becomes chaotic. Saturn coorbital satellites are synchronous since $ν\llσ$, but coorbital exoplanets may present non-synchronous or chaotic rotation. Our results prove that the spin dynamics of a body cannot be dissociated from its orbital environment. We further anticipate that a similar mechanism may affect the rotation of bodies in any mean-motion resonance.
△ Less
Submitted 3 December, 2013; v1 submitted 26 November, 2013;
originally announced November 2013.
-
On the co-orbital motion of two planets in quasi-circular orbits
Authors:
Philippe Robutel,
Alexandre Pousse
Abstract:
We develop an analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. The Hamiltonian, averaged over one of the planetary mean longitude, is expanded in power series of eccentricities and inclinations. The model, which is valid in the entire co-orbital region, possesses an integrable approximation modeling the planar and quasi-circular motions. F…
▽ More
We develop an analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. The Hamiltonian, averaged over one of the planetary mean longitude, is expanded in power series of eccentricities and inclinations. The model, which is valid in the entire co-orbital region, possesses an integrable approximation modeling the planar and quasi-circular motions. First, focusing on the fixed points of this approximation, we highlight relations linking the eigenvectors of the associated linearized differential system and the existence of certain remarkable orbits like the elliptic Eulerian Lagrangian configurations, the Anti-Lagrange (Giuppone et al., 2010) orbits and some second sort orbits discovered by Poincaré. Then, the variational equation is studied in the vicinity of any quasi-circular periodic solution. The fundamental frequencies of the trajectory are deduced and possible occurrence of low order resonances are discussed. Finally, with the help of the construction of a Birkhoff normal form, we prove that the elliptic Lagrangian equilateral configurations and the Anti-Lagrange orbits bifurcate from the same fixed point L4.
△ Less
Submitted 26 April, 2013; v1 submitted 3 April, 2013;
originally announced April 2013.
-
Influence of the coorbital resonance on the rotation of the Trojan satellites of Saturn
Authors:
Philippe Robutel,
Nicolas Rambaux,
Maryame El Moutamid
Abstract:
The Cassini spacecraft collects high resolution images of the saturnian satellites and reveals the surface of these new worlds. The shape and rotation of the satellites can be determined from the Cassini Imaging Science Subsystem data, employing limb coordinates and stereogrammetric control points. This is the case for Epimetheus (Tiscareno et al. 2009) that opens elaboration of new rotational mod…
▽ More
The Cassini spacecraft collects high resolution images of the saturnian satellites and reveals the surface of these new worlds. The shape and rotation of the satellites can be determined from the Cassini Imaging Science Subsystem data, employing limb coordinates and stereogrammetric control points. This is the case for Epimetheus (Tiscareno et al. 2009) that opens elaboration of new rotational models (Tiscareno et al. 2009; Noyelles 2010; Robutel et al. 2011). Especially, Epimetheus is characterized by its horseshoe shape orbit and the presence of the swap is essential to introduce explicitly into rotational models. During its journey in the saturnian system, Cassini spacecraft accumulates the observational data of the other satellites and it will be possible to determine the rotational parameters of several of them. To prepare these future observations, we built rotational models of the coorbital (also called Trojan) satellites Telesto, Calypso, Helene, and Polydeuces, in addition to Janus and Epimetheus. Indeed, Telesto and Calypso orbit around the L_4 and L_5 Lagrange points of Saturn-Tethys while Helene and Polydeuces are coorbital of Dione. The goal of this study is to understand how the departure from the Keplerian motion induced by the perturbations of the coorbital body, influences the rotation of these satellites. To this aim, we introduce explicitly the perturbation in the rotational equations by using the formalism developed by Erdi (1977) to represent the coorbital motions, and so we describe the rotational motion of the coorbitals, Janus and Epimetheus included, in compact form.
△ Less
Submitted 13 February, 2012; v1 submitted 4 November, 2011;
originally announced November 2011.
-
Analytical description of physical librations of Saturnian coorbital satellites Janus and Epimetheus
Authors:
Philippe Robutel,
Nicolas Rambaux,
Julie Castillo-Rogez
Abstract:
Janus and Epimetheus are famously known for their distinctive horseshoe-shaped orbits resulting from a 1:1 orbital resonance. Every four years these two satellites swap their orbits by a few tens of kilometers as a result of their close encounter. Recently Tiscareno et al. (2009) have proposed a model of rotation based on images from the Cassini orbiter. These authors inferred the amplitude of rot…
▽ More
Janus and Epimetheus are famously known for their distinctive horseshoe-shaped orbits resulting from a 1:1 orbital resonance. Every four years these two satellites swap their orbits by a few tens of kilometers as a result of their close encounter. Recently Tiscareno et al. (2009) have proposed a model of rotation based on images from the Cassini orbiter. These authors inferred the amplitude of rotational librational motion in longitude at the orbital period by fitting a shape model to the recent Cassini ISS images. By a quasiperiodic approximation of the orbital motion, we describe how the orbital swap impacts the rotation of the satellites. To that purpose, we have developed a formalism based on quasi-periodic series with long and short-period librations. In this framework, the amplitude of the libration at the orbital period is found proportional to a term accounting for the orbital swap. We checked the analytical quasi-periodic development by performing a numerical simulation and find both results in good agreement. To complete this study, the results regarding the short-period librations are studied with the help of an adiabatic-like approach.
△ Less
Submitted 10 September, 2010; v1 submitted 2 March, 2010;
originally announced March 2010.
-
Frozen Orbits at high eccentricity and inclination: Application to Mercury orbiter
Authors:
N. Delsate,
P. Robutel,
A. Lemaitre,
T. Carletti
Abstract:
We hereby study the stability of a massless probe orbiting around an oblate central body (planet or planetary satellite) perturbed by a third body, assumed to lie in the equatorial plane (Sun or Jupiter for example) using an Hamiltonian formalism.
We are able to determine, in the parameters space, the location of the frozen orbits, namely orbits whose orbital elements remain constant on averag…
▽ More
We hereby study the stability of a massless probe orbiting around an oblate central body (planet or planetary satellite) perturbed by a third body, assumed to lie in the equatorial plane (Sun or Jupiter for example) using an Hamiltonian formalism.
We are able to determine, in the parameters space, the location of the frozen orbits, namely orbits whose orbital elements remain constant on average, to characterize their stability/unstability and to compute the periods of the equilibria.
The proposed theory is general enough, to be applied to a wide range of probes around planet or natural planetary satellites.
The BepiColombo mission is used to motivate our analysis and to provide specific numerical data to check our analytical results.
Finally, we also bring to the light that the coefficient J_2 is able to protect against the increasing of the eccentricity due to the Kozai-Lidov effect.
△ Less
Submitted 1 March, 2010;
originally announced March 2010.
-
The resonant structure of Jupiter's trojan asteroids-II. What happens for different configurations of the planetary system
Authors:
Philippe Robutel,
Julien Bodossian
Abstract:
In a previous paper, we have found that the resonance structure of the present Jupiter Trojan swarms could be split up into four different families of resonances. Here, in a first step, we generalize these families in order to describe the resonances occurring in Trojan swarms embedded in a generic planetary system. The location of these families changes under a modification of the fundamental f…
▽ More
In a previous paper, we have found that the resonance structure of the present Jupiter Trojan swarms could be split up into four different families of resonances. Here, in a first step, we generalize these families in order to describe the resonances occurring in Trojan swarms embedded in a generic planetary system. The location of these families changes under a modification of the fundamental frequencies of the planets and we show how the resonant structure would evolve during a planetary migration. We present a general method, based on the knowledge of the fundamental frequencies of the planets and on those that can be reached by the Trojans, which makes it possible to predict and localize the main events arising in the swarms during migration. In particular, we show how the size and stability of the Trojan swarms are affected by the modification of the frequencies of the planets. Finally, we use this method to study the global dynamics of the Jovian Trojan swarms when Saturn migrates outwards. Besides the two resonances found by Morbidelli et al (2005) which could have led to the capture of the current population just after the crossing of the 2:1 orbital resonance, we also point out several sequences of chaotic events that can influence the Trojan population.
△ Less
Submitted 12 June, 2009; v1 submitted 20 September, 2008;
originally announced September 2008.
-
High order symplectic integrators for perturbed Hamiltonian systems
Authors:
J. Laskar,
P. Robutel
Abstract:
We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+εB$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps with a remainder of order $O(τ^pε+τ^2ε^2)$, where $τ$ is the stepsize of the integrator. The analytical expressions of the leading terms of the remainders are given at…
▽ More
We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+εB$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps with a remainder of order $O(τ^pε+τ^2ε^2)$, where $τ$ is the stepsize of the integrator. The analytical expressions of the leading terms of the remainders are given at all orders. In many cases, a corrector step can be performed such that the remainder becomes $O(τ^pε+τ^4ε^2)$. The performances of these integrators are compared for the simple pendulum and the planetary 3-Body problem of Sun-Jupiter-Saturn.
△ Less
Submitted 4 May, 2000;
originally announced May 2000.
