Asteroid (514107) Kaepaokaawela is the first example of an object in the 1/1 mean motion resonance with Jupiter with retrograde motion around the Sun. Its orbit was shown to be stable over the age of the Solar System which implies that it must have been captured from another star when the Sun was still in its birth cluster. Kaepaokaawela orbit is also located at the peak of the capture probability in the coorbital resonance. Identifying the periodic orbits that Kaepaokaawela and similar asteroids followed during their evolution is an important step towards precisely understanding their capture mechanism. Here, we find the families of periodic orbits in the two-dimensional retrograde coorbital problem and analyze their stability and bifurcations into three-dimensional periodic orbits. Our results explain the radical differences observed in 2D and 3D coorbital capture simulations. In particular, we find that analytical and numerical results obtained for planar motion are not always valid at infinitesimal deviations from the plane.

Asteroid (514107) 2015 BZ509 was discovered recently in Jupiter's co-orbital region with a retrograde motion around the Sun. The known chaotic dynamics of the outer Solar System have so far precluded the identification of its origin. Here, we perform a high-resolution statistical search for stable orbits and show that asteroid (514107) 2015 BZ509 has been in its current orbital state since the formation of the Solar System. This result indicates that (514107) 2015 BZ509 was captured from the interstellar medium 4.5 billion years in the past as planet formation models cannot produce such a primordial large-inclination orbit with the planets on nearly-coplanar orbits interacting with a coplanar debris disk that must produce the low-inclination small-body reservoirs of the Solar System such as the asteroid and Kuiper belts. This result also implies that more extrasolar asteroids are currently present in the Solar System on nearly-polar orbits.

The classical disturbing function of the three-body problem widely used in planetary dynamics studies is an expansion of the gravitational interaction of the three-body problem with respect to zero eccentricity and zero inclination. This restricts its validity to nearly coplanar orbits. Motivated by the dynamical study of asteroids, Centaurs and transneptunian objects with arbitrary inclinations, we derive a series expansion of the gravitational interaction with respect to an arbitrary reference inclination that generalises our work on the polar and retrograde disturbing functions. The new disturbing function, like the polar one, may model any resonance as expansion order is unrelated to resonance order. The powers of eccentricity and inclination of the force amplitude of a $p$:$q$ resonance depend only on the parity of the resonance order $|p-q|$. Disturbing functions with non zero reference inclinations are thus physically different from the classical disturbing function as the former are based on the three-dimensional three-body problem and the latter on the two-dimensional one. We illustrate the use of the new disturbing function by showing that what is known as pure eccentricity resonances are intrinsically dependent on inclination contrary to the prediction of the classical disturbing function. We determine the inclination dependence of the resonance widths of the 2:1 and 3:1 prograde and retrograde inner resonances with Jupiter as well as those of the asymmetric librations of the 1:2 and 1:3 prograde outer resonances with Neptune.

Capture in mean motion resonance has been observed in the Solar System for small objects with prograde as well as retrograde orbits of moderate inclinations. However, no example of an object with a nearly polar orbit was known to be in resonance with a planet. In this Letter, we report that the nearly-polar transneptunian object (471325), nicknamed Niku, is in a 7:9 resonance with Neptune, with a mean lifetime in resonance of 16 +- 11 million years. While entrance and exit in the 7:9 resonance is caused by close encounters with Neptune the resonant configuration provides a temporary protection mechanism against disruptive close encounters with this planet. The other nearly polar transneptunian objects do not seem to be in resonance with the planets with the possible exception of 2008 KV42, also known as Drac, that has a small chance of being in the 8:13 resonance with Neptune.

The classical disturbing function of the three-body problem is based on an expansion of the gravitational interaction in the vicinity of nearly coplanar orbits. Consequently, it is not suitable for the identification and study of resonances of the Centaurs and transneptunian objects on nearly polar orbits with the solar system planets. Here, we provide a series expansion algorithm of the gravitational interaction in the vicinity of polar orbits and produce explicitly the disturbing function to fourth order in eccentricity and inclination cosine. The properties of the polar series differ significantly from those of the classical disturbing function: the polar series can model any resonance as the expansion order is not related to the resonance order. The powers of eccentricity and inclination of the force amplitude of a $p$:$q$ resonance do not depend on the value of the resonance order $|p-q|$ but only on its parity. Thus all even resonance order eccentricity amplitudes are $\propto e^2$ and odd ones $\propto e$ to lowest order in eccentricity $e$. With the new findings on the structure of the polar disturbing function and the possible resonant critical arguments, we illustrate the dynamics of the polar resonances 1:3, 3:1, 2:9 and 7:9 where transneptunian object 471325 could currently be locked.

The effect of radial drift rate on mean motion resonance capture is studied for prograde, polar and retrograde orbits. We employ the numerical framework of our earlier exploration of resonance capture at arbitrary inclination. Randomly constructed samples of massless particles are set to migrate radially from outside the orbit of a Jupiter-mass planet at different drift rates totalling more than $1.6\times 10^6$ numerical simulations. Slower drift rates reduce overall capture probability especially for prograde orbits and enhance capture at specific initial inclinations of high order resonances such as the outer 1:5, 1:4, 1:3, 2:5, 3:7 and 5:7. Global capture is reduced with increasing eccentricity at all inclinations as high order resonances capture more particles that are subsequently lost by disruptive close encounters with the planet. The relative efficiency of retrograde resonances at long-lived capture with respect to prograde resonances is explained by the reduced effect of planet encounters as such events occur with a shorter duration and a higher relative velocity for retrograde motion. Capture in the coorbital 1:1 resonance is marginally affected by the radial drift rate except for nearly coplanar retrograde eccentric orbits whose capture likelihood is increased significantly with slower drift rates. An unexpected finding is the presence of a dynamical corridor for capture in high order inner prograde resonances with initial inclinations in the range [50$^\circ$,80$^\circ$] especially at the inner 5:2 resonance whose capture likelihood peaks at 80% to 90% depending on initial eccentricity.

Motivated by the dynamics of resonance capture, we study numerically the coorbital resonance for inclination180 >=I>=0 in the circular restricted three-body problem. We examine the similarities and differences between planar and three dimensional coorbital resonance capture and seek their origin in the stability of coorbital motion at arbitrary inclination. After we present stability maps of the planar prograde and retrograde coorbital resonances, we characterize the new coorbital modes in three dimensions. We see that retrograde mode I (R1) and mode II (R2) persist as we change the relative inclination, while retrograde mode III (R3) seems to exist only in the planar problem. A new coorbital mode (R4) appears in 3D which is a retrograde analogue to an horseshoe-orbit. The Kozai-Lidov resonance is active for retrograde orbits as well as prograde orbits and plays a key role in coorbital resonance capture. Stable coorbital modes exist at all inclinations, including retrograde and polar obits. This result confirms the robustness the coorbital resonance at large inclination and encourages the search for retrograde coorbital companions of the solar system's planets.

Resonance capture is studied numerically in the three-body problem for arbitrary inclinations. Massless particles are set to drift from outside the 1:5 resonance with a Jupiter-mass planet thereby encountering the web of the planet's diverse mean motion resonances. Randomly constructed samples explore parameter space for inclinations from 0 to 180 deg with 5deg increments totalling nearly 6x10^5 numerical simulations. Thirty resonances internal and external to the planet's location are monitored. We find that retrograde resonances are unexpectedly more efficient at capture than prograde resonances and that resonance order is not necessarily a good indicator of capture efficiency at arbitrary inclination. Capture probability drops significantly at moderate sample eccentricity for initial inclinations in the range [10deg,110deg]. Orbit inversion is possible for initially circular orbits with inclinations in the range [60deg,130deg]. Capture in the 1:1 coorbital resonance occurs with great likelihood at large retrograde inclinations. The planet's orbital eccentricity, if larger than 0.1, reduces the capture probabilities through the action of the eccentric Kozai-Lidov mechanism. A capture asymmetry appears between inner and outer resonances as prograde orbits are preferentially trapped in inner resonances. The relative capture efficiency of retrograde resonance suggests that the dynamical lifetimes of Damocloids and Centaurs on retrograde orbits must be significantly larger than those on prograde orbits implying that the recently identified asteroids in retrograde resonance, 2006 BZ8, 2008 SO218, 2009 QY6 and 1999 LE31(Morais and Namouni, 2013, MNRAS 436, L30) may be among the oldest small bodies that wander between the outer giant planets.

We test a crossing orbit stability criterion for eccentric planetary systems, based on Wisdom's criterion of first order mean motion resonance overlap (Wisdom, 1980). We show that this criterion fits the stability regions in real exoplanet systems quite well. In addition, we show that elliptical orbits can remain stable even for regions where the apocenter distance of the inner orbit is larger than the pericenter distance of the outer orbit, as long as the initial orbits are aligned. The analytical expressions provided here can be used to put rapid constraints on the stability zones of multi-planetary systems. As a byproduct of this research, we further show that the amplitude variations of the eccentricity can be used as a fast-computing stability indicator.

We continue the investigation of the dynamics of retrograde resonances initiated in Morais & Giuppone (2012). After deriving a procedure to deduce the retrograde resonance terms from the standard expansion of the three-dimensional disturbing function, we concentrate on the planar problem and construct surfaces of section that explore phase-space in the vicinity of the main retrograde resonances (2/-1, 1/-1 and 1/-2). In the case of the 1/-1 resonance for which the standard expansion is not adequate to describe the dynamics, we develop a semi-analytic model based on numerical averaging of the unexpanded disturbing function, and show that the predicted libration modes are in agreement with the behavior seen in the surfaces of section.

We identify a set of asteroids among Centaurs and Damocloids, that orbit contrary to the common direction of motion in the Solar System and that enter into resonance with Jupiter and Saturn. Their orbits have inclinations I >= 140 deg and semi-major axes a < 15 AU. Two objects are currently in retrograde resonance with Jupiter: 2006 BZ8 in the 2/-5 resonance and 2008 SO218 in the 1/-2 resonance. One object, 2009 QY6, is currently in the 2/-3 retrograde resonance with Saturn. These are the first examples of Solar System objects in retrograde resonance. The present resonant configurations last for several thousand years. Brief captures in retrograde resonance with Saturn are also possible during the 20,000 years integration timespan, particularly in the 1/-1 resonance (2006 BZ8) and the 9/-7 resonance (1999 LE31).

Hierarchical two-planet systems, in which the inner body's semi-major axis is between 0.1 and 0.5 AU, usually present high eccentricity values, at least for one of the orbits. As a result of the formation process, one may expect that planetary systems with high eccentricities also have high mutual inclinations. However, here we show that tidal effects combined with gravitational interactions damp the initial mutual inclination to modest values in timescales that are shorter than the age of the system. This effect is not a direct consequence of tides on the orbits, but it results from a secular forcing of the inner planet's flattening. We then conclude that these hierarchical planetary systems are unlikely to present very high mutual inclinations, at least as long as the orbits remain outside the Lidov-Kozai libration areas. The present study can also be extended to systems of binary stars and to planet-satellite systems.

We present a new observational campaign, DWARF, aimed at detection of circumbinary extrasolar planets using the timing of the minima of low-mass eclipsing binaries. The observations will be performed within an extensive network of relatively small to medium-size telescopes with apertures of ~20-200 cm. The starting sample of the objects to be monitored contains (i) low-mass eclipsing binaries with M and K components, (ii) short-period binaries with sdB or sdO component, and (iii) post-common-envelope systems containing a WD, which enable to determine minima with high precision. Since the amplitude of the timing signal increases with the orbital period of an invisible third component, the timescale of project is long, at least 5-10 years. The paper gives simple formulas to estimate suitability of individual eclipsing binaries for the circumbinary planet detection. Intrinsic variability of the binaries (photospheric spots, flares, pulsation etc.) limiting the accuracy of the minima timing is also discussed. The manuscript also describes the best observing strategy and methods to detect cyclic timing variability in the minima times indicating presence of circumbinary planets. First test observation of the selected targets are presented.

We investigate the stability of prograde versus retrograde planets in circular binary systems using numerical simulations. We show that retrograde planets are stable up to distances closer to the perturber than prograde planets. We develop an analytical model to compute the prograde and retrograde mean motion resonances' locations and separatrices. We show that instability is due to single resonance forcing, or caused by nearby resonances' overlap. We validate our results regarding the role of single resonances and resonances' overlap on orbit stability, by computing surfaces of section of the CR3BP. We conclude that the observed enhanced stability of retrograde planets with respect to prograde planets is due to essential differences between the phase-space topology of retrograde versus prograde resonances (at p/q mean motion ratio, prograde resonance is of order p - q while retrograde resonance is of order p + q).

The radial velocity method for detecting extra-solar planets relies on measuring the star's wobble around the system's center of mass. Since this is an indirect method, we may ask if there are other dynamical effects that can mimic such wobble. In recent articles\cite{1,2,3}, we modeled the effect of a nearby binary system on a star's radial velocity. We showed that, if we are unaware of this nearby binary, for instance because one component is unresolved or both components are faint stars, the binary's effect may mimic a planet. Here, we review this work, explaining in which circumstances the binary's effect may mimic a planet and we discuss what can be done in practice in order to distinguish between these two scenarios (planet or nearby binary).

The HD\,196885 system is composed of a binary star and a planet orbiting the primary. The orbit of the binary is fully constrained by astrometry, but for the planet the inclination with respect to the plane of the sky and the longitude of the node are unknown. Here we perform a full analysis of the HD\,196885 system by exploring the two free parameters of the planet and choosing different sets of angular variables. We find that the most likely configurations for the planet is either nearly coplanar orbits (prograde and retrograde), or highly inclined orbits near the Lidov-Kozai equilibrium points, i = 44^{\circ} or i = 137^{\circ} . Among coplanar orbits, the retrograde ones appear to be less chaotic, while for the orbits near the Lidov-Kozai equilibria, those around \omega= 270^{\circ} are more reliable, where \omega_k is the argument of pericenter of the planet's orbit with respect to the binary's orbit. From the observer's point of view (plane of the sky) stable areas are restricted to (I1, \Omega_1) \sim (65^{\circ}, 80^{\circ}), (65^{\circ},260^{\circ}), (115^{\circ},80^{\circ}), and (115^{\circ},260^{\circ}), where I1 is the inclination of the planet and \Omega_1 is the longitude of ascending node.

We model the secular evolution of a star's orbit when it has a nearby binary system. We assume a hierarchical triple system where the inter-binary distance is small in comparison with the distance to the star. We show that the major secular effect is precession of the star's orbit around the binary system's centre of mass. We explain how we can obtain this precession rate from the star's radial velocity data, and thus infer the binary system's parameters. We show that the secular effect of a nearby binary system on the star's radial velocity can sometimes mimic a planet. We analyze the radial velocity data for {\nu}-octantis A which has a nearby companion ({\nu}-octantis B) and we obtain retrograde precession of (-0.86 \pm 0.02)\degree/yr. We show that if {\nu}-octantis B was itself a double star, it could mimic a signal with similarities to that previously identified as a planet of {\nu}-octantis A. Nevertheless, we need more observations in order to decide in favor of the double star hypothesis.

Most extrasolar planets currently known were discovered by means of an indirect method that measures the stellar wobble caused by the planet. We previously studied a triple system composed of a star and a nearby binary on circular coplanar orbits. We showed that although the effect of the binary on the star can be differentiated from the stellar wobble caused by a planet, because of observational limitations the two effects may often remain indistinguishable. Here, we develop a model that applies to eccentric and inclined orbits. We show that the binary's effect is more likely to be mistaken by planet(s) in the case of coplanar motion observed equator-on. Moreover, when the orbits are eccentric, the magnitude of the binary's effect may be larger than in the circular case. Additionally, an eccentric binary can mimic two planets with orbital periods in the ratio 2/1. However, when the star's orbit around the binary's center of mass has a high eccentricity and a reasonably well-constrained period, it should be easier to distinguish the binary's effect from a planet.

The traditional method for detecting extra-solar planets relies on measuring a small stellar wobble which is assumed to be caused by a planet orbiting the star. Recently, it was suggested that a similar stellar wobble could be caused by a close binary system (Schneider and Cabrera, 2006). Here we show that, although the effect of a close binary system can at first sight be mistaken as a planetary companion to the star, more careful analysis of the observational data should allow us to distinguish between the two effects.