High order symplectic integrators for perturbed Hamiltonian systems

We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+\epsilon B$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps with a remainder of order $O(\tau^p\epsilon +\tau^2\epsilon^2)$, where $\tau$ 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(\tau^p\epsilon +\tau^4\epsilon^2)$. The performances of these integrators are compared for the simple pendulum and the planetary 3-Body problem of Sun-Jupiter-Saturn.Comment: 24 pages, 6 figurre

Where are the Uranus Trojans?

The area of stable motion for fictitious Trojan asteroids around Uranus' equilateral equilibrium points is investigated with respect to the inclination of the asteroid's orbit to determine the size of the regions and their shape. For this task we used the results of extensive numerical integrations of orbits for a grid of initial conditions around the points L4 and L5, and analyzed the stability of the individual orbits. Our basic dynamical model was the Outer Solar System (Jupiter, Saturn, Uranus and Neptune). We integrated the equations of motion of fictitious Trojans in the vicinity of the stable equilibrium points for selected orbits up to the age of the Solar system of 5 billion years. One experiment has been undertaken for cuts through the Lagrange points for fixed values of the inclinations, while the semimajor axes were varied. The extension of the stable region with respect to the initial semimajor axis lies between 19.05 < a < 19.3 AU but depends on the initial inclination. In another run the inclination of the asteroids' orbit was varied in the range 0 < i < 60 and the semimajor axes were fixed. It turned out that only four 'windows' of stable orbits survive: these are the orbits for the initial inclinations 0 < i < 7, 9 < i < 13, 31 < i < 36 and 38 < i < 50. We postulate the existence of at least some Trojans around the Uranus Lagrange points for the stability window at small and also high inclinations.Comment: 15 pages, 12 figures, submitted to CMD

Influence of the coorbital resonance on the rotation of the Trojan satellites of Saturn

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

On the coplanar eccentric non-restricted co-orbital dynamics

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.publishe

(1173) Anchises - Thermophysical and Dynamical Studies of a Dynamically Unstable Jovian Trojan

We have performed detailed thermophysical and dynamical modelling of Jovian Trojan (1173) Anchises. Our results reveal a most unusual object. By examining observational data taken by IRAS, Akari and WISE between 11.5 and 60 microns, along with variations in its optical lightcurve, we find Anchises is most likely an elongated body, with an axes-ratio of ~1.4. This yields calculated best-fit dimensions of 170x121x121km (an equivalent diameter of 136+18/-11km). We find the observations are best fit by Anchises having a retrograde sense of rotation, and an unusually high thermal inertia (25 to 100 Jm-2s-0.5K-1). The geometric albedo is found to be 0.027 (+0.006/-0.007). Anchises therefore has one of the highest published thermal inertias of any object larger than 100km in diameter, at such large heliocentric distances, and is one of the lowest albedo objects ever observed. More observations are needed to see if there is a link between the very shallow phase curve, with almost no opposition effect, and the derived thermal properties for this large Trojan asteroid. Our dynamical investigation of Anchises' orbit has revealed it to be dynamically unstable on timescales of hundreds of Myr, similar to the unstable Neptunian Trojans 2001 QR322 and 2008 LC18. Unlike those objects, we find that Anchises' dynamical stability is not a function of its initial orbital elements, the result of the exceptional precision with which its orbit is known. This is the first time that a Jovian Trojan has been shown to be dynamically unstable, and adds weight to the idea that planetary Trojans represent a significant ongoing contribution to the Centaur population, the parents of the short-period comets. The observed instability does not rule out a primordial origin for Anchises, but when taken in concert with the result of our thermophysical analysis, suggest that it would be a fascinating target for future study.Comment: 5 figures, 3 tables, accepted for publication in Monthly Notices of the Royal Astronomical Societ

Planetary system LHS 1140 revisited with ESPRESSO and TESS

Context. LHS 1140 is an M dwarf known to host two transiting planets at orbital periods of 3.77 and 24.7 days. They were detected with HARPS and Spitzer. 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. All these properties place this system at the forefront of the habitable exoplanet exploration, and it therefore constitutes a relevant case for further astrobiological studies, including atmospheric observations. Aims. We further characterize this system by improving the physical and orbital properties of the known planets, search for additional planetary-mass components in the system, and explore the possibility of co-orbitals. Methods. We collected 113 new high-precision radial velocity observations with ESPRESSO over a 1.5-yr time span with an average photon-noise precision of 1.07 m s-1. We performed an extensive analysis of the HARPS and ESPRESSO datasets and also analyzed them together with the new TESS photometry. We analyzed the Bayesian evidence of several models with different numbers of planets and orbital configurations. Results. We significantly improve our knowledge of the properties of the known planets LHS 1140 b (Pb ∼ 24.7 days) and LHS 1140 c (Pc ∼ 3.77 days). We determine new masses with a precision of 6% for LHS 1140 b (6.48 ± 0.46 Mpdbl) and 9% for LHS 1140 c (mc = 1.78 ± 0.17 Mpdbl). This reduces the uncertainties relative to previously published values by half. Although both planets have Earth-like bulk compositions, the internal structure analysis suggests that LHS 1140 b might be iron-enriched and LHS 1140 c might be a true Earth twin. 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 ± 650 km for the habitable-zone planet LHS 1140 b. Our results also provide evidence for a new planet candidate in the system (md = 4.8 ± 1.1Mpdbl) 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 Mpdbl for periods shorter than 10 days and above 2 Mpdbl for periods up to one year. Finally, our co-orbital analysis discards co-orbital planets in the tadpole and horseshoe configurations of LHS 1140 b down to 1 Mpdbl with a 95% confidence level (twice better than with the previous HARPS dataset). Indications for a possible co-orbital signal in LHS 1140 c are detected in both radial velocity (alternatively explained by a high eccentricity) and photometric data (alternatively explained by systematics), however. Conclusions. The new precise measurements of the planet properties of the two transiting planets in LHS 1140 as well as the detection of the planet candidate LHS 1140 d make this system a key target for atmospheric studies of rocky worlds at different stellar irradiations.With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737

The TROY project: II. Multi-technique constraints on exotrojans in nine planetary systems

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. Aims. In this paper, we present new ground-based observations searching for these bodies co-orbiting with nine close-in (P < 5 days) planets, using various observing techniques. The combination of all of these techniques allows us to restrict the parameter space of any possible trojan in the system. Methods. We used multi-Technique observations, comprised of 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. Results. 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, especially in the trojan mass or radius vs. libration amplitude plane. In particular, we can set upper mass limits in the super-Earth mass regime for six of the studied systemspublishe

A Semi-Analytic Algorithm for Constructing Lower Dimensional Elliptic Tori in Planetary Systems

We adapt the Kolmogorov's normalization algorithm (which is the key element of the original proof scheme of the KAM theorem) to the construction of a suitable normal form related to an invariant elliptic torus. As a byproduct, our procedure can also provide some analytic expansions of the motions on elliptic tori. By extensively using algebraic manipulations on a computer, we explicitly apply our method to a planar four-body model not too different with respect to the real Sun--Jupiter--Saturn--Uranus system. The frequency analysis method allows us to check that our location of the initial conditions on an invariant elliptic torus is really accurate.Comment: 31 pages, 4 figure

The TROY project: Searching for co-orbital bodies to known planets: I. Project goals and first results from archival radial velocity

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.publishe