We present the discovery and characterization of WASP-148, a new extrasolar system that includes at least two giant planets. The host star is a slowly rotating inactive late-G dwarf with a V=12 magnitude. The planet WASP-148b is a hot Jupiter of 0.72 R_Jup and 0.29 M_Jup that transits its host with an orbital period of 8.80 days. We found the planetary candidate with the SuperWASP photometric survey, then characterized it with the SOPHIE spectrograph. Our radial velocity measurements subsequently revealed a second planet in the system, WASP-148c, with an orbital period of 34.5 days and a minimum mass of 0.40 M_Jup. No transits of this outer planet were detected. The orbits of both planets are eccentric and fall near the 4:1 mean-motion resonances. This configuration is stable on long timescales, but induces dynamical interactions so that the orbits differ slightly from purely Keplerian orbits. In particular, WASP-148b shows transit-timing variations of typically 15 minutes, making it the first interacting system with transit-timing variations that is detected on ground-based light curves. We establish that the mutual inclination of the orbital plane of the two planets cannot be higher than 35 degrees, and the true mass of WASP-148c is below 0.60 M_Jup. We present photometric and spectroscopic observations of this system that cover a time span of ten years. We also provide their Keplerian and Newtonian analyses; these analyses should be significantly improved through future TESS~observations.

Most Neptune-mass planets in close-in orbits (orbital periods less than a few days) present nonzero eccentricity, typically around 0.15. This is somehow unexpected, as these planets undergo strong tidal dissipation that should circularize their orbits in a time-scale shorter than the age of the system. In this paper we discuss some mechanisms that can oppose to bodily tides, namely, thermal atmospheric tides, evaporation of the atmosphere, and excitation from a distant companion. In the first two cases, the eccentricity can increase consistently, while in the last one, the eccentricity can only be excited for a limited amount of time (that may nevertheless exceed the age of the system). We show the limitations of these different mechanisms and how some of them could, depending on specific properties of the observed planetary systems, account for their presently observed eccentricities.

This white book stems from the contributions presented at the Portuguese SKA Days, held on the 6th and 7th February 2018 with the presence of the SKA Deputy Director General Alistair McPherson and the SKA Science Director Robert Braun. This initiative was held to promote the Square Kilometer Array (SKA) - the world's largest radio telescope - among the Portuguese scientific and business communities with support from the Portuguese Science and Technology Foundation (FCT) with the contribution of Portuguese policy makers and researchers. The meeting was very successful in providing a detailed overview of the SKA status, vision and goals and describes most of the Portuguese contributions to science, technology and the related industry aspirations

We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

The main purpose of this work is to present a time-domain implementation of the Andrade rheology, instead of the traditional expansion in terms of a Fourier series of the tidal potential. This approach can be used in any fully three dimensional numerical simulation of the dynamics of a system of many deformable bodies. In particular, it allows large eccentricities, large mutual inclinations, and it is not limited to quasi-periodic perturbations. It can take into account an extended class of perturbations, such as chaotic motions, transient events, and resonant librations. The results are presented by means of a concrete application: the analysis of the libration of Enceladus. This is done by means of both analytic formulas in the frequency domain and direct numerical simulations. We do not a priori assume that Enceladus has a triaxial shape, the eventual triaxiality is a consequence of the satellite motion and its rheology. As a result we obtain an analytic formula for the amplitude of libration that incorporates a new correction due to the rheology. Our results provide an estimation of the amplitude of libration of the core of Enceladus as 0.6% of that of the shell. They also reproduce the observed 10 GW of tidal heat generated by Enceladus with a value of $0.17\times 10^{14}$Pa$\cdot$s for the global effective viscosity under both Maxwell and Andrade rheology.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

In some planetary systems, the orbital periods of two of its members present a commensurability, usually known by mean-motion resonance. These resonances greatly enhance the mutual gravitational influence of the planets. As a consequence, these systems present uncommon behaviors, and their motions need to be studied with specific methods. Some features are unique and allow us a better understanding and characterization of these systems. Moreover, mean-motion resonances are a result of an early migration of the orbits in an accretion disk, so it is possible to derive constraints on their formation. Here we review the dynamics of a pair of resonant planets and explain how their orbits evolve in time. We apply our results to the HD 45365 planetary system.

The Carrington event is considered to be one of the most extreme space weather events in observational history within a series of magnetic storms caused by extreme interplanetary coronal mass ejections (ICMEs) from a large and complex active region (AR) emerged on the solar disk. In this article, we study the temporal and spatial evolutions of the source sunspot active region and visual aurorae, and compare this storm with other extreme space weather events on the basis of their spatial evolution. Sunspot drawings by Schwabe, Secchi, and Carrington describe the position and morphology of the source AR at that time. Visual auroral reports from the Russian Empire, Iberia, Ireland, Oceania, and Japan fill the spatial gap of auroral visibility and revise the time series of auroral visibility in mid to low magnetic latitudes (MLATs). The reconstructed time series is compared with magnetic measurements and shows the correspondence between low to mid latitude aurorae and the phase of magnetic storms. The spatial evolution of the auroral oval is compared with those of other extreme space weather events in 1872, 1909, 1921, and 1989 as well as their storm intensity, and contextualizes the Carrington event, as one of the most extreme space weather events, but likely not unique.

We present the obliquity measurement, that is, the angle between the normal angle of the orbital plane and the stellar spin axis, of the sub-Neptune planet HD3167 c, which transits a bright nearby K0 star. We study the orbital architecture of this multi-planet system to understand its dynamical history. We also place constraints on the obliquity of planet d based on the geometry of the planetary system and the dynamical study of the system. New observations obtained with HARPS-N at the Telescopio Nazionale Galileo (TNG) were employed for our analysis. The sky-projected obliquity was measured using three different methods: the Rossiter-McLaughlin anomaly, Doppler tomography, and reloaded Rossiter-McLaughlin techniques. We performed the stability analysis of the system and investigated the dynamical interactions between the planets and the star. HD3167 c is found to be nearly polar with sky-projected obliquity, $\lambda$ = -97$^{\circ} \pm$ 23$^{\circ}$. This misalignment of the orbit of planet c with the spin axis of the host star is detected with 97\% confidence. The analysis of the dynamics of this system yields coplanar orbits of planets c and d. It also shows that it is unlikely that the currently observed system can generate this high obliquity for planets c and d by itself. However, the polar orbits of planets c and d could be explained by the presence of an outer companion in the system. Follow-up observations of the system are required to confirm such a long-period companion.

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 $\pi/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.

Short-period planets suffer from extreme tidal forces from their parent stars causing them to deform and attain non-spherical shapes. The planet shapes, modeled here as triaxial ellipsoids, can have an impact on the observed transit light curves and the planetary parameters derived. We investigate the detectability of tidal deformation in short-period planets from their transit light curves and the instrumental precision needed. We show how the detection of deformation from the light curve allows us to obtain an observational estimate of the second fluid Love number which gives valuable insight about the planet's internal structure. We adopted a model to calculate the shape of a planet due to the external potentials acting on it and used this model to modify the ellc transit tool. Our model is parameterized by the Love number, hence for a given light curve we can derive the value of the Love number that best matches the observations. We simulated the known cases of WASP-103b and WASP-121b expected to be highly deformed. Our analyses showed that instrumental precision $\leq$50ppm/min is needed to reliably estimate the Love number and detect tidal deformation. This precision can be achieved for WASP-103b in ~40 transits using HST and in ~300 transits using the forthcoming CHEOPS. However, fewer transits will be required for short-period planets that may be found around bright stars in the TESS and PLATO survey missions. The unprecedented precisions expected from PLATO and JWST can permit the detection of deformation with a single transit. However, the effects of instrumental and astrophysical noise must be well-considered as they can increase the number of transits required to reach the 50 ppm/min detection limit. We also show that improper modeling of limb darkening can act to bury signals related to the planet's shape thereby leading us to infer sphericity for a deformed planet.

Many triple-star systems have an inner pair with an orbital period of a few days only. A common mechanism to explain the short-period pile-up present in the observations is the migration through Lidov-Kozai cycles combined with tidal friction. Here, we revisit this mechanism and aim to determine the initial orbital configurations leading to this process. We show that the mutual inclination of the triple-star system is not the only critical parameter, since the eccentricity as well as the argument of the pericenter of the inner orbit also play an important role in the establishment of the Lidov-Kozai migration. Our framework is the secular hierarchical three-body problem (octupole order approximation) with general relativity corrections, including the effects of tides, stellar oblateness and magnetic spin-down braking. Both the orbital and the spin evolutions are considered. Extensive numerical simulations with uniform and non-uniform distributions of the initial orbital parameters are carried out, and unbiased initial conditions leading to Lidov-Kozai migration are revealed. Finally, we highlight the importance of the initial "Kozai constant" $h=\sqrt{(1-e^2)}\cos{i}$ in the dynamical evolution of triple-star systems, by showing that phase portraits at given $h$-values unveil different evolution paths.

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.

Stellar systems consisting of multiple stars tend to undergo tidal interactions when the separations between the stars are short. While tidal phenomena have been extensively studied, a certain tidal effect exclusive to hierarchical triples (triples in which one component star has a much wider orbit than the others) has hardly received any attention, mainly due to its complexity and consequent resistance to being modelled. This tidal effect is the tidal perturbation of the tertiary by the inner binary, which in turn depletes orbital energy from the inner binary, causing the inner binary separation to shrink. In this paper, we develop a fully numerical simulation of these "tertiary tides" by modifying established tidal models. We also provide general insight as to how close a hierarchical triple needs to be in order for such an effect to take place, and demonstrate that our simulations can effectively retrieve the orbital evolution for such systems. We conclude that tertiary tides are a significant factor in the evolution of close hierarchical triples, and strongly influence at least ~ 1% of all multiple star systems.

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.

We investigate the dynamics of the (47171) Lempo triple system, also known by 1999TC$_{36}$. We derive a full 3D $N$-body model that takes into account the orbital and spin evolution of all bodies, which are assumed triaxial ellipsoids. We show that, for reasonable values of the shapes and rotational periods, the present best fitted orbital solution for the Lempo system is chaotic and unstable in short time-scales. The formation mechanism of this system is unknown, but the orbits can be stabilised when tidal dissipation is taken into account. The dynamics of the Lempo system is very rich, but depends on many parameters that are presently unknown. A better understanding of this systems thus requires more observations, which also need to be fitted with a complete model like the one presented here.

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.

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.

The stability of planets in the alpha-Centauri AB stellar system has been studied extensively. However, most studies either focus on the orbital plane of the binary or consider inclined circular orbits. Here, we numerically investigate the stability of a possible planet in the alpha-Centauri AB binary system for S-type orbits in an arbitrary spatial configuration. In particular, we focus on inclined orbits and explore the stability for different eccentricities and orientation angles. We show that large stable and regular regions are present for very eccentric and inclined orbits, corresponding to libration in the Lidov-Kozai resonance. We additionally show that these extreme orbits can survive over the age of the system, despite the effect of tides. Our results remain qualitatively the same for any compact binary system.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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.

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 $\alpha$, which is proportional to the mass ratio $m_2/m_1$. Therefore, when $\alpha$ 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.

Aims:We aim to significantly increase the number of detected extra-solar planets in a magnitude-limited sample to improve our knowledge of their orbital element distributions and thus obtain better constraints for planet-formation models. Methods: Radial-velocity data were taken at Haute-Provence Observatory (OHP, France) with the ELODIE echelle spectrograph. Results: We report the presence of a planet orbiting HD 196885 A, with an orbital period of 1349 days. This star was previously suggested to host a 386 -day planet, but we cannot confirm its existence. We also detect the presence of a stellar companion, HD 196885 B, and give some constraints on its orbit.

Cassini states correspond to equilibria of the spin axis of a body when its orbit is perturbed. They were initially described for satellites, but the spin axis of stars and planets undergoing strong dissipation can also evolve into some equilibria. For small satellites, the rotational angular momentum is usually much smaller than the total angular momentum, so classical methods for finding Cassini states rely on this approximation. Here we present a more general approach, which is valid for the secular quadrupolar non-restricted problem with spin. Our method is still valid when the precession rate and the mutual inclination of the orbits are not constant. Therefore, it can be used to study stars with close-in companions, or planets with heavy satellites, like the Earth-Moon system.

In this paper, we present a formalism designed to model tidal interaction with a viscoelastic body made of Maxwell material. Our approach remains regular for any spin rate and orientation, and for any orbital configuration including high eccentricities and close encounters. The method is to integrate simultaneously the rotation and the position of the planet as well as its deformation. We provide the equations of motion both in the body frame and in the inertial frame. With this study, we generalize preexisting models to the spatial case and to arbitrary multipole orders using a formalism taken from quantum theory. We also provide the vectorial expression of the secular tidal torque expanded in Fourier series. Applying this model to close-in exoplanets, we observe that if the relaxation time is longer than the revolution period, the phase space of the system is characterized by the presence of several spin-orbit resonances, even in the circular case. As the system evolves, the planet spin can visit different spin-orbit configurations. The obliquity is decreasing along most of these resonances, but we observe a case where the planet tilt is instead growing. These conclusions derived from the secular torque are successfully tested with numerical integrations of the instantaneous equations of motion on HD 80606b. Our formalism is also well adapted to close-in super-Earths in multiplanet systems which are known to have non-zero mutual inclinations.

The competition between the torques induced by solid and thermal tides drives the rotational dynamics of Venus-like planets and super-Earths orbiting in the habitable zone of low-mass stars. The tidal responses of the atmosphere and telluric core are related to their respective physical properties and strongly depend on the tidal frequency. The resulting torque determines the possible equilibrium states of the planet's spin. We compute here an analytic expression for the total tidal torque exerted on a Venus-like planet. This expression is used to characterize the equilibrium rotation of the body. Close to the star, the solid tide dominates. Far from it, the thermal tide drives the rotational dynamics of the planet. The transition regime corresponds to the habitable zone, where prograde and retrograde equilibrium states appear. We demonstrate the strong impact of the atmospheric properties and of the rheology of the solid part on the rotational dynamics of Venus-like planets, highlighting the key role played by dissipative mechanisms in the stability of equilibrium configurations.

We present a novel approach for analysing radial velocity data that combines two features: all the planets are searched at once and the algorithm is fast. This is achieved by utilizing compressed sensing techniques, which are modified to be compatible with the Gaussian processes framework. The resulting tool can be used like a Lomb-Scargle periodogram and has the same aspect but with much fewer peaks due to aliasing. The method is applied to five systems with published radial velocity data sets: HD 69830, HD 10180, 55 Cnc, GJ 876 and a simulated very active star. The results are fully compatible with previous analysis, though obtained more straightforwardly. We further show that 55 Cnc e and f could have been respectively detected and suspected in early measurements from the Lick observatory and Hobby-Eberly Telescope available in 2004, and that frequencies due to dynamical interactions in GJ 876 can be seen.

We investigate the orbital and rotational evolution of the CoRoT-7 two-planet system, assuming that the innermost planet behaves like a Maxwell body. We numerically resolve the coupled differential equations governing the instantaneous deformation of the inner planet together with the orbital motion of the system. We show that, depending on the relaxation time for the deformation of the planet, the orbital evolution has two distinct behaviours: for relaxation times shorter than the orbital period, we reproduce the results from classic tidal theories, for which the eccentricity is always damped. However, for longer relaxation times, the eccentricity of the inner orbit is secularly excited and can grow to high values. This mechanism provides an explanation for the present high eccentricity observed for CoRoT-7 b, as well as for other close-in super-Earths in multiple planetary systems.

We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin-orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star-planet-satellite systems and triple stellar systems.

Aims: Several studies suggest that the activity level of a planet-host star can be influenced by the presence of a close-by orbiting planet. Moreover, the interaction mechanisms that have been proposed, magnetic interaction and tidal interaction, exhibit a very different dependence on orbital separation between the star and the planet. A detection of activity enhancement and characterization of its dependence on planetary orbital distance can, in principle, allow us to characterize the physical mechanism behind the activity enhancement. Methods: We used the HARPS-N spectrograph to measure the stellar activity level of HD 80606 during the planetary periastron passage and compared the activity measured to that close to apastron. Being characterized by an eccentricity of 0.93 and an orbital period of 111 days, the system's extreme variation in orbital separation makes it a perfect target to test our hypothesis. Results: We find no evidence for a variation in the activity level of the star as a function of planetary orbital distance, as measured by all activity indicators employed log($R'_{HK}$), H$_\alpha$, NaI, and HeI. None of the models employed, whether magnetic interaction or tidal interaction, provides a good description of the data. Conclusions: We find no evidence for star-planet interaction in HD\,80606 at the moment of the periastron passage of its very eccentric planet. The straightforward explanation for the non-detection is the absence of interaction as a result of a low magnetic field strength on either the planet or the star and of the low level of tidal interaction between the two. However, we cannot exclude two scenarios: i) the interaction can be instantaneous and of magnetic origin, being concentrated on the substellar point and its surrounding area, and ii) the interaction can lead to a delayed activity enhancement. (abridged)

With about 2000 extrasolar planets confirmed, the results show that planetary systems have a whole range of unexpected properties. We present a full investigation of the HD219828 system, a bright metal-rich star for which a hot neptune has previously been detected. We used a set of HARPS, SOPHIE, and ELODIE radial velocities to search for the existence of orbiting companions to HD219828. A dynamical analysis is also performed to study the stability of the system and to constrain the orbital parameters and planet masses. We announce the discovery of a long period (P=13.1years) massive (msini=15.1MJup) companion (HD219828c) in a very eccentric orbit (e=0.81). The same data confirms the existence of a hot-neptune, HD219828b, with a minimum mass of 21 MEarth and a period of 3.83days. The dynamical analysis shows that the system is stable. The HD219828 system is extreme and unique in several aspects. First, among all known exoplanet systems it presents an unusually high mass ratio. We also show that systems like HD219828, with a hot neptune and a long-period massive companion are more frequent than similar systems with a hot jupiter instead. This suggests that the formation of hot neptunes follows a different path than the formation of their hot jovian counterparts. The high mass, long period, and eccentricity of HD219828c also make it a good target for Gaia astrometry as well as a potential target for atmospheric characterisation, using direct imaging or high-resolution spectroscopy. Astrometric observations will allow us to derive its real mass and orbital configuration. If a transit of HD219828b is detected, we will be able to fully characterise the system, including the relative orbital inclinations. With a clearly known mass, HD219828c may become a benchmark object for the range in between giant planets and brown dwarfs.

Extrasolar planets abound in almost any possible configuration. However, until five years ago, there was a lack of planets orbiting closer than 0.5 au to giant or subgiant stars. Since then, recent detections have started to populated this regime by confirming 13 planetary systems. We discuss the properties of these systems in terms of their formation and evolution off the main sequence. Interestingly, we find that $70.0\pm6.6$ % of the planets in this regime are inner components of multiplanetary systems. This value is 4.2$\sigma$ higher than for main-sequence hosts, which we find to be $42.4\pm0.1$ %. The properties of the known planets seem to indicate that the closest-in planets (a < 0.06 au) to main-sequence stars are massive (i.e., hot Jupiters) and isolated and that they are subsequently engulfed by their host as it evolves to the red giant branch, leaving only the predominant population of multiplanetary systems in orbits 0.06 < a < 0.5 au. We discuss the implications of this emerging observational trend in the context of formation and evolution of hot Jupiters.

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.

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.

We present new radial velocity measurements of eight stars secured with the spectrograph SOPHIE at the 193-cm telescope of the Haute-Provence Observatory allowing the detection and characterization of new giant extrasolar planets. The host stars are dwarfs of spectral types between F5 and K0 and magnitudes between 6.7 and 9.6; the planets have minimum masses M_p sin i between 0.4 to 3.8 M_Jup and orbital periods of several days to several months. The data allow only single planets to be discovered around the first six stars (HD143105, HIP109600, HD35759, HIP109384, HD220842, and HD12484), but one of them shows the signature of an additional substellar companion in the system. The seventh star, HIP65407, allows the discovery of two giant planets, just outside the 12:5 resonance in weak mutual interaction. The last star, HD141399, was already known to host a four-planetary system; our additional data and analyses allow new constraints to be put on it. We present Keplerian orbits of all systems, together with dynamical analyses of the two multi-planetary systems. HD143105 is one of the brightest stars known to host a hot Jupiter, which could allow numerous follow-up studies to be conducted despite this is not a transiting system. The giant planets HIP109600b, HIP109384b, and HD141399c are located in the habitable zone of their host star.

Cassini states correspond to the equilibria of the spin axis of a body when its orbit is perturbed. They were initially described for planetary satellites, but the spin axes of black hole binaries also present this kind of equilibria. In previous works, Cassini states were reported as spin-orbit resonances, but actually the spin of black hole binaries is in circulation and there is no resonant motion. Here we provide a general description of the spin dynamics of black hole binary systems based on a Hamiltonian formalism. In absence of dissipation the problem is integrable and it is easy to identify all possible trajectories for the spin for a given value of the total angular momentum. As the system collapses due to radiation reaction, the Cassini states are shifted to different positions, which modifies the dynamics around them. This is why the final spin distribution may differ from the initial one. Our method provides a simple way of predicting the distribution of the spin of black hole binaries at the end of the inspiral phase.

Here we show how to determine the orbital parameters of a system composed of a star and N companions (that can be planets, brown-dwarfs or other stars), using a simple Fourier analysis of the radial velocity data of the star. This method supposes that all objects in the system follow keplerian orbits around the star and gives better results for a large number of observational points. The orbital parameters may present some errors, but they are an excellent starting point for the traditional minimization methods such as the Levenberg-Marquardt algorithms.

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass $0.1~M_\odot$, Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water.

In this paper we explore the possibility that the recently detected reflected light signal of 51\,Peg\,b could be caused by a ring system around the planet. We use a simple model to compare the observed signal with the expected signal from a short-period giant planet with rings. We also use simple dynamical arguments to understand the possible geometry of such a system. We provide evidence that, to a good approximation, the observations are compatible with the signal expected from a ringed planet, assuming that the rings are non-coplanar with the orbital plane. However, based on dynamical arguments, we also show that this configuration is unlikely. In the case of coplanar rings we then demonstrate that the incident flux on the ring surface is about 2\% the value received by the planet, a value that renders the ring explanation unlikely. The results suggest that the signal observed cannot in principle be explained by a planet+ring system. We discuss, however, the possibility of using reflected light spectra to detect and characterize the presence of rings around short-period planets. Finally, we show that ring systems could have already been detected by photometric transit campaigns, but their signal could have been easily misinterpreted by the expected light curve of an eclipsing binary.

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\nu/2$, where $\nu = n_b - n$, and $k$ is an integer. Moreover, when the libration at natural frequency has the same magnitude as $\nu$, 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.

We study the stability of mean-motion resonances (MMR) between two planets during their migration in a protoplanetary disk. We use an analytical model of resonances, and describe the effect of the disk by a migration timescale (T_{m,i}) and an eccentricity damping timescale (T_{e,i}) for each planet (i=1,2 respectively for the inner and outer planet). We show that the resonant configuration is stable if T_{e,1}/T_{e,2} > (e_1/e_2)^2. This general result can be used to put constraints on specific models of disk-planet interactions. For instance, using classical prescriptions for type I migration, we show that when the angular momentum deficit (AMD) of the inner orbit is larger than the outer's orbit AMD, resonant systems must have a locally inverted disk density profile to stay locked in resonance during the migration. This inversion is very untypical of type I migration and our criterion can thus provide an evidence against classical type I migration. That is indeed the case for the Jupiter-mass resonant systems HD 60532b, c (3:1 MMR), GJ 876b, c (2:1 MMR), and HD 45364b, c (3:2 MMR). This result may be an evidence for type II migration (gap opening planets), which is compatible with the large masses of these planets.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

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.

In this paper we present a new approach to tidal theory. Assuming a Maxwell viscoelastic rheology, we compute the instantaneous deformation of celestial bodies using a differential equation for the gravity field coefficients. This method allows large eccentricities and it is not limited to quasi-periodic perturbations. It can take into account an extended class of perturbations, including chaotic motions and transient events. We apply our model to some already detected eccentric hot Jupiters and super-Earths in planar configurations. We show that when the relaxation time of the deformation is larger than the orbital period, spin-orbit equilibria arise naturally at half-integers of the mean motion, even for gaseous planets. In the case of super-Earths, these equilibria can be maintained for very low values of eccentricity. Our method can also be used to study planets with complex internal structures and other rheologies.

A large number of observed exoplanets are part of multiple planet systems. Most of these systems are sufficiently close-in to be tidally evolving. In such systems, there is a competition between the excitation caused by planet-planet interactions and tidal damping. Using as an example two multiple planet systems, which host planets in the surface liquid water habitable zone (HZ): Kepler-62 and Kepler-186, we show the importance and effect of both planetary and stellar tides on the dynamical evolution of planets and on the climate of the HZ planets.

Planets orbiting very close to their host stars have been found, some of them on the verge of tidal disruption. The ellipsoidal shape of these planets can significantly differ from a sphere, which modifies the transit light curves. Here we present an easy method for taking the effect of the tidal bulge into account in the transit photometric observations. We show that the differences in the light curve are greater than previously thought. When detectable, these differences provide us an estimation of the fluid Love number, which is invaluable information on the internal structure of close-in planets. We also derive a simple analytical expression to correct the bulk density of these bodies, that can be 20% smaller than current estimates obtained assuming a spherical radius.

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.

In this paper we report a new transiting warm giant planet: KOI-1257 b. It was first detected in photometry as a planet-candidate by the ${\it Kepler}$ space telescope and then validated thanks to a radial velocity follow-up with the SOPHIE spectrograph. It orbits its host star with a period of 86.647661 d $\pm$ 3 s and a high eccentricity of 0.772 $\pm$ 0.045. The planet transits the main star of a metal-rich, relatively old binary system with stars of mass of 0.99 $\pm$ 0.05 Msun and 0.70 $ \pm $ 0.07 Msun for the primary and secondary, respectively. This binary system is constrained thanks to a self-consistent modelling of the ${\it Kepler}$ transit light curve, the SOPHIE radial velocities, line bisector and full-width half maximum (FWHM) variations, and the spectral energy distribution. However, future observations are needed to confirm it. The PASTIS fully-Bayesian software was used to validate the nature of the planet and to determine which star of the binary system is the transit host. By accounting for the dilution from the binary both in photometry and in radial velocity, we find that the planet has a mass of 1.45 $ \pm $ 0.35 Mjup, and a radius of 0.94 $ \pm $ 0.12 Rjup, and thus a bulk density of 2.1 $ \pm $ 1.2 g.cm$^{-3}$. The planet has an equilibrium temperature of 511 $\pm$ 50 K, making it one of the few known members of the warm-jupiter population. The HARPS-N spectrograph was also used to observe a transit of KOI-1257 b, simultaneously with a joint amateur and professional photometric follow-up, with the aim of constraining the orbital obliquity of the planet. However, the Rossiter-McLaughlin effect was not clearly detected, resulting in poor constraints on the orbital obliquity of the planet.

Planets with masses between 0.1 - 10 M_earth are believed to host dense atmospheres. These atmospheres can play an important role on the planet's spin evolution, since thermal atmospheric tides, driven by the host star, may counterbalance gravitational tides. In this work we study the long-term spin evolution of Earth-sized exoplanets. We generalize previous works by including the effect of eccentric orbits and obliquity. We show that under the effect of tides and core-mantle friction, the obliquity of the planets evolve either to 0 or 180 degrees. The rotation of these planets is also expected to evolve into a very restricted number of equilibrium configurations. In general, none of this equilibria is synchronous with the orbital mean motion. The role of thermal atmospheric tides becomes more important for Earth-sized planets in the habitable zones of their systems, so they cannot be neglected when we search for their potential habitability.

We study the evolution of two planets around a star, in mean-motion resonance and undergoing tidal effect. We derive an integrable analytical model of mean-motion resonances of any order which reproduce the main features of the resonant dynamics. Using this simplified model, we obtain a criterion showing that depending on the balance of the tidal dissipation in both planets, their final period ratio may stay at the resonant value, increase above, or decrease below the resonant value. Applying this criterion to the two inner planets orbiting GJ163, we deduce that the current period ratio (2.97) could be the outcome of dissipation in the 3:1 MMR provided that the innermost planet is gaseous (slow dissipation) while the second one is rocky (faster dissipation). We perform N-body simulations with tidal dissipation to confirm the results of our analytical model. We also apply our criterion on GJ581b, c (5:2 MMR) and reproduce the current period ratio (2.4) if the inner planet is gaseous and the outer is rocky (as for GJ163). Finally, we apply our model to the Kepler mission's statistics. We show that the excess of planets pairs close to first order MMR but in external circulation, i.e., with period ratios P_out/P_in > (p+1)/p for the resonance (p+1):p, can be reproduced by tidal dissipation in the inner planet. There is no need for any other dissipative mechanism, provided that these systems left the resonance with non-negligible eccentricities.

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\nu/2$, where $n$ is the orbital mean motion, $\nu$ the orbital libration frequency, and $k$ an integer. In addition, when the natural rotational libration frequency due to the axial asymmetry, $\sigma$, has the same magnitude as $\nu$, the rotation becomes chaotic. Saturn coorbital satellites are synchronous since $\nu\ll\sigma$, 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.

We explore the relations between physical and orbital properties of planets and properties of their host stars to identify the main observable signatures of the formation and evolution processes of planetary systems. We use a large sample of FGK dwarf planet hosts with stellar parameters derived in a homogeneous way from the SWEET-Cat database to study the relation between stellar metallicity and position of planets in the period-mass diagram. In the second part we use all the RV-detected planets orbiting FGK stars to explore the role of planet-disk and planet-planet interaction on the evolution of orbital properties of planets with masses above 1MJup. We show that planets orbiting metal-poor stars have longer periods than those in metal-rich systems. This trend is valid for masses at least from 10MEarth to 4MJup. Earth-like planets orbiting metal-rich stars always show shorter periods (fewer than 20 days) than those orbiting metal-poor stars. We also found statistically significant evidence that very high mass giants have on average more eccentric orbits than giant planets with lower mass.Finally, we show that the eccentricity of planets with masses higher than 4MJup tends to be lower for planets with shorter periods. Our results suggest that the planets in the P-MP diagram are evolving differently because of a mechanism that operates over a wide range of planetary masses. This mechanism is stronger or weaker depending on the metallicity of the respective system. One possibility is that planets in metal-poor disks form farther out from their central star and/or they form later and do not have time to migrate as far as the planets in metal-rich systems. The trends and dependencies obtained for very high mass planetary systems suggest that planet-disk interaction is a very important and orbit-shaping mechanism for planets in the high-mass domain. Shortened.

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.

The meter-per-second precision achieved by today velocimeters enables the search for 1-10 M_Earth planets in the habitable zone of cool stars. This paper reports on the detection of 3 planets orbiting GJ163 (HIP19394), a M3 dwarf monitored by our ESO/HARPS search for planets. We made use of the HARPS spectrograph to collect 150 radial velocities of GJ163 over a period of 8 years. We searched the RV time series for coherent signals and found 5 distinct periodic variabilities. We investigated the stellar activity and casted doubts on the planetary interpretation for 2 signals. Before more data can be acquired we concluded that at least 3 planets are orbiting GJ163. They have orbital periods of P_b=8.632+-0.002, P_c=25.63+-0.03 and P_d=604+-8 days and minimum masses msini = 10.6+-0.6, 6.8+-0.9, and 29+-3 M_Earth, respectively. We hold our interpretations for the 2 additional signals with periods P_(e)=19.4 and P_(f)=108 days. The inner pair presents an orbital period ratio of 2.97, but a dynamical analysis of the system shows that it lays outside the 3:1 mean motion resonance. GJ163c, in particular, is a super-Earth with an equilibrium temperature of T_eq = (302+-10) (1-A)^(1/4) K and may lie in the so called habitable zone for albedo values (A=0.34-0.89) moderately higher than that of Earth (A_Earth=0.2-0.3).

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.

Many exoplanets have been observed close to their parent stars with orbital periods of a few days. As for the major satellites of the Jovian planets, the figure of these planets is expected to be strongly shaped by tidal forces. However, contrarily to Solar System satellites, exoplanets may present high values for the obliquity and eccentricity due to planetary perturbations, and may also be captured in spin-orbit resonances different from the synchronous one. Here we give a general formulation of the equilibrium figure of those bodies, that makes no particular assumption on the spin and/or orbital configurations. The gravity field coefficients computed here are well suited for describing the figure evolution of a body whose spin and orbit undergo substantial variations in time.

Close-in planetary systems detected by the Kepler mission present an excess of periods ratio that are just slightly larger than some low order resonant values. This feature occurs naturally when resonant couples undergo dissipation that damps the eccentricities. However, the resonant angles appear to librate at the end of the migration process, which is often believed to be an evidence that the systems remain in resonance. Here we provide an analytical model for the dissipation in resonant planetary systems valid for low eccentricities. We confirm that dissipation accounts for an excess of pairs that lie just aside from the nominal periods ratios, as observed by the Kepler mission. In addition, by a global analysis of the phase space of the problem, we demonstrate that these final pairs are non-resonant. Indeed, the separatrices that exist in the resonant systems disappear with the dissipation, and remains only a circulation of the orbits around a single elliptical fixed point. Furthermore, the apparent libration of the resonant angles can be explained using the classical secular averaging method. We show that this artifact is only due to the severe damping of the amplitudes of the eigenmodes in the secular motion.

Impact basins identified by Mariner 10 and Messenger flyby images provide us a fossilized record of the impactor flux of asteroids on Mercury during the last stages of the early Solar System. The distribution of these basins is not uniform across the surface, and is consistent with a primordial synchronous rotation (Wieczorek et al. 2012). By analyzing the size of the impacts, we show that the distribution for asteroid diameters D < 110 km is compatible with an index power law of 1.2, a value that matches the predicted primordial distribution of the main-belt. We then derive a simple collisional model coherent with the observations, and when combining it with the secular evolution of the spin of Mercury, we are able to reproduce the present 3/2 spin-orbit resonance (about 50% of chances), as well as a primordial synchronous rotation. This result is robust with respect to variations in the dissipation and collisional models, or in the initial spin state of the planet.

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.

The planet Mercury is locked in a spin-orbit resonance where it rotates three times about its spin axis for every two orbits about the Sun. The current explanation for this unique state assumes that the initial rotation of this planet was prograde and rapid, and that tidal torques decelerated the planetary spin to this resonance. When core-mantle boundary friction is accounted for, capture into the 3/2 resonance occurs with a 26% probability, but the most probable outcome is capture into one of the higher-order resonances. Here we show that if the initial rotation of Mercury were retrograde, this planet would be captured into synchronous rotation with a 68% probability. Strong spatial variations of the impact cratering rate would have existed at this time, and these are shown to be consistent with the distribution of pre-Calorian impact basins observed by Mariner 10 and MESSENGER. Escape from this highly stable resonance is made possible by the momentum imparted by large basin-forming impact events, and capture into the 3/2 resonance occurs subsequently under favourable conditions.

Short period planets are subject to intense energetic irradiations from their stars. It has been shown that this can lead to significant atmospheric mass-loss and create smaller mass planets. Here, we analyse whether the evaporation mechanism can affect the orbit of planets. The orbital evolution of a planet undergoing evaporation is derived analytically in a very general way. Analytical results are then compared with the period distribution of two classes of inner exoplanets: Jupiter-mass planets and Neptune-mass planets. These two populations have a very distinct period distribution, with a probability lower than 10^-4 that they were derived from the same parent distribution. We show that mass ejection can generate significant migration with an increase of orbital period that matches very well the difference of distribution of the two populations. This would happen if the evaporation emanates from above the hottest region of planet surface. Thus, migration induced by evaporation is an important mechanism that cannot be neglected.

Planets close to their host stars are believed to undergo significant tidal interactions, leading to a progressive damping of the orbital eccentricity. Here we show that, when the orbit of the planet is excited by an outer companion, tidal effects combined with gravitational interactions may give rise to a secular increasing drift on the eccentricity. As long as this secular drift counterbalances the damping effect, the eccentricity can increase to high values. This mechanism may explain why some of the moderate close-in exoplanets are observed with substantial eccentricity values.

We present here in full details the linear secular theory with tidal damping that was used to constraint the fit of the HD10180 planetary system in (Lovis et al. 2011). The theory is very general and can provide some intuitive understanding of the final state of a planetary system when one or more planets are close to their central star. We globally recover the results of (Mardling 2007), but we show that in the HD209458 planetary system, the consideration of the tides raised by the central star on the planet lead to believe that the eccentricity of HD209458b is most probably much smaller than 0.01.

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.

We investigate the dynamical evolution of hierarchical three-body systems under the effect of tides, when the ratio of the orbital semi-major axes is small and the mutual inclination is relatively large (greater than 20 degrees). Using the quadrupolar non-restricted approximation for the gravitational interactions and the viscous linear model for tides, we derive the averaged equations of motion in a vectorial formalism which is suitable to model the long-term evolution of a large variety of exoplanetary systems in very eccentric and inclined orbits. In particular, it can be used to derive constraints for stellar spin-orbit misalignment, capture in Cassini states, tidal-Kozai migration, or damping of the mutual inclination. Because our model is valid for the non-restricted problem, it can be used to study systems of identical mass or for the outer restricted problem, such as the evolution of a planet around a binary of stars. Here, we apply our model to three distinct situations: 1) the HD80606 planetary system, for which we obtain the probability density function distribution for the misalignment angle, with two pronounced peaks of higher probability around 53 and 109 degrees; 2) the HD98800 binary system, for which we show that initial prograde orbits inside the observed disc may become retrograde and vice-versa, only because of tidal migration within the binary stars; 3) the HD11964 planetary system, for which we show that tidal dissipation combined with gravitational perturbations may lead to a decrease in the mutual inclination, and a fast circularization of the inner orbit.

Understanding the consequences of the gravitational interaction between a star and a planet is fundamental to the study of exoplanets. The solution of the two-body problem shows that the planet moves in an elliptical path around the star and that each body moves in an ellipse about the common center of mass. The basic properties of such a system are derived from first principles and described in the context of detecting exoplanets.

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.

Context. Low-mass extrasolar planets are presently being discovered at an increased pace by radial velocity and transit surveys, opening a new window on planetary systems. Aims. We are conducting a high-precision radial velocity survey with the HARPS spectrograph which aims at characterizing the population of ice giants and super-Earths around nearby solar-type stars. This will lead to a better understanding of their formation and evolution, and yield a global picture of planetary systems from gas giants down to telluric planets. Methods. Progress has been possible in this field thanks in particular to the sub-m/s radial velocity precision achieved by HARPS. We present here new high-quality measurements from this instrument. Results. We report the discovery of a planetary system comprising at least five Neptune-like planets with minimum masses ranging from 12 to 25 M_Earth, orbiting the solar-type star HD 10180 at separations between 0.06 and 1.4 AU. A sixth radial velocity signal is present at a longer period, probably due to a 65-M_Earth object. Moreover, another body with a minimum mass as low as 1.4 M_Earth may be present at 0.02 AU from the star. This is the most populated exoplanetary system known to date. The planets are in a dense but still well-separated configuration, with significant secular interactions. Some of the orbital period ratios are fairly close to integer or half-integer values, but the system does not exhibit any mean-motion resonances. General relativity effects and tidal dissipation play an important role to stabilize the innermost planet and the system as a whole. Numerical integrations show long-term dynamical stability provided true masses are within a factor ~3 from minimum masses. We further note that several low-mass planetary systems exhibit a rather "packed" orbital architecture with little or no space left for additional planets. (Abridged)

Tidal effects arise from differential and inelastic deformation of a planet by a perturbing body. The continuous action of tides modify the rotation of the planet together with its orbit until an equilibrium situation is reached. It is often believed that synchronous motion is the most probable outcome of the tidal evolution process, since synchronous rotation is observed for the majority of the satellites in the Solar System. However, in the 19th century, Schiaparelli also assumed synchronous motion for the rotations of Mercury and Venus, and was later shown to be wrong. Rather, for planets in eccentric orbits synchronous rotation is very unlikely. The rotation period and axial tilt of exoplanets is still unknown, but a large number of planets have been detected close to the parent star and should have evolved to a final equilibrium situation. Therefore, based on the Solar System well studied cases, we can make some predictions for exoplanets. Here we describe in detail the main tidal effects that modify the secular evolution of the spin and the orbit of a planet. We then apply our knowledge acquired from Solar System situations to exoplanet cases. In particular, we will focus on two classes of planets, "Hot-Jupiters" (fluid) and "Super-Earths" (rocky with atmosphere).

Precise radial-velocity measurements for data acquired with the HARPS spectrograph infer that three planets orbit the M4 dwarf star GJ876. In particular, we confirm the existence of planet "d", which orbits every 1.93785 days. We find that its orbit may have significant eccentricity (e=0.14), and deduce a more accurate estimate of its minimum mass of 6.3 Earth masses. Dynamical modeling of the HARPS measurements combined with literature velocities from the Keck Observatory strongly constrain the orbital inclinations of the "b" and "c" planets. We find that i_b = 48.9 degrees and i_c = 48.1 degrees, which infers the true planet masses of M_b = 2.64 Jupiter masses and M_c = 0.83 Jupiter masses, respectively. Radial velocities alone, in this favorable case, can therefore fully determine the orbital architecture of a multi-planet system, without the input from astrometry or transits. The orbits of the two giant planets are nearly coplanar, and their 2:1 mean motion resonance ensures stability over at least 5 Gyr. The libration amplitude is smaller than 2 degrees, suggesting that it was damped by some dissipative process during planet formation. The system has space for a stable fourth planet in a 4:1 mean motion resonance with planet "b", with a period around 15 days. The radial velocity measurements constrain the mass of this possible additional planet to be at most that of the Earth.

Neptune Trojans and Plutinos are two subpopulations of trans-Neptunian objects located in the 1:1 and the 3:2 mean motion resonances with Neptune, respectively, and therefore protected from close encounters with the planet. However, the orbits of these two kinds of objects may cross very often, allowing a higher collisional rate between them than with other kinds of trans-Neptunian objects, and a consequent size distribution modification of the two subpopulations. Observational colors and absolute magnitudes of Neptune Trojans and Plutinos show that i) there are no intrinsically bright (large) Plutinos at small inclinations, ii) there is an apparent excess of blue and intrinsically faint (small) Plutinos, and iii) Neptune Trojans possess the same blue colors as Plutinos within the same (estimated) size range do. For the present subpopulations we analyzed the most favorable conditions for close encounters/collisions and address any link there could be between those encounters and the sizes and/or colors of Plutinos and Neptune Trojans. We also performed a simultaneous numerical simulation of the outer Solar System over 1 Gyr for all these bodies in order to estimate their collisional rate. We conclude that orbital overlap between Neptune Trojans and Plutinos is favored for Plutinos with large libration amplitudes, high eccentricities, and small inclinations. Additionally, with the assumption that the collisions can be disruptive creating smaller objects not necessarily with similar colors, the present high concentration of small Plutinos with small inclinations can thus be a consequence of a collisional interaction with Neptune Trojans and such hypothesis should be further analyzed.

Long-term precise Doppler measurements with the CORALIE spectrograph revealed the presence of two massive companions to the solar-type star HD202206. Although the three-body fit of the system is unstable, it was shown that a 5:1 mean motion resonance exists close to the best fit, where the system is stable. We present here an extensive dynamical study of the HD202206 system aiming at constraining the inclinations of the two known companions, from which we derive possible ranges of value for the companion masses. We study the long term stability of the system in a small neighborhood of the best fit using Laskar's frequency map analysis. We also introduce a numerical method based on frequency analysis to determine the center of libration mode inside a mean motion resonance. We find that acceptable coplanar configurations are limited to inclinations to the line of sight between 30 and 90 degrees. This limits the masses of both companions to roughly twice the minimum. Non coplanar configurations are possible for a wide range of mutual inclinations from 0 to 90 degrees, although $\Delta\Omega = 0 [\pi]$ configurations seem to be favored. We also confirm the 5:1 mean motion resonance to be most likely. In the coplanar edge-on case, we provide a very good stable solution in the resonance, whose $\chi^2$ does not differ significantly from the best fit. Using our method to determine the center of libration, we further refine this solution to obtain an orbit with a very low amplitude of libration, as we expect dissipative effects to have dampened the libration.

Some of the satellites in the Solar System, including the Moon, appear to have been captured from heliocentric orbits at some point in their past, and then have evolved to the present configurations. The exact process of how this trapping occurred is unknown, but the dissociation of a planetesimal binary in the gravitational field of the planet, gas drag, or a massive collision seem to be the best candidates. However, all these mechanisms leave the satellites in elliptical orbits that need to be damped to the present almost circular ones. Here we give a complete description of the secular tidal evolution of a satellite just after entering a bounding state with the planet. In particular, we take into account the spin evolution of the satellite, which has often been assumed synchronous in previous studies. We apply our model to Triton and successfully explain some geophysical properties of this satellite, as well as the main dynamical features observed for the Neptunian system.

The present obliquity of Mercury is very low (less than 0.1 degree), which led previous studies to always adopt a nearly zero obliquity during the planet's past evolution. However, the initial orientation of Mercury's rotation axis is unknown and probably much different than today. As a consequence, we believe that the obliquity could have been significant when the rotation rate of the planet first encountered spin-orbit resonances. In order to compute the capture probabilities in resonance for any evolutionary scenario, we present in full detail the dynamical equations governing the long term evolution of the spin, including the obliquity contribution. The secular spin evolution of Mercury results from tidal interactions with the Sun, but also from viscous friction at the core-mantle boundary. Here, this effect is also regarded with particular attention. Previous studies show that a liquid core enhances drastically the chances of capture in spin-orbit resonances. We confirm these results for null obliquity, but we find that the capture probability generally decreases as the obliquity increases. We finally show that, when core-mantle friction is combined with obliquity evolution, the spin can evolve into some unexpected configurations as the synchronous or the 1/2 spin-orbit resonance.

We present a uniform study of the chemical abundances of 12 elements (Si, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Na, Mg, and Al) derived from the spectra of 451 stars observed as part of one of the HARPS GTO planet search programs. Sixty eight of these are planet-bearing stars. The main goals of our work are: i) the investigation of possible differences between the abundances of stars with and without planets; ii) the study of the possible differences in the abundances of stars in the thin and the thick disc. We confirm that there is a systematically higher metallicity in planet host stars, when compared to non planet-hosts, common to all studied species. We also found that there is no difference in the galactic chemical evolution trends of the stars with and without planets. Stars that harbour planetary companions simply appear to be in the high metallicity tail of the distribution. We also confirm that Neptunian and super-Earth class planets may be easier to find at lower metallicities. A statistically significative abundance difference between stars of the thin and the thick disc was found for [Fe/H] $<$ 0. However, the populations from the thick and the thin disc cannot be clearly separated.

In a recent paper it was reported a planetary system around the star HD60532, composed by two giant planets in a possible 3:1 mean motion resonance, that should be confirmed within the next decade. Here we show that the analysis of the global dynamics of the system allows to confirm this resonance. The present best fit to data already corresponds to this resonant configuration and the system is stable for at least 5Gry. The 3:1 resonance is so robust that stability is still possible for a wide variety of orbital parameters around the best fit solution and also if the inclination of the system orbital plane with respect to the plane of the sky is as small as 15 deg. Moreover, if the inclination is taken as a free parameter in the adjustment to the observations, we find an inclination ~ 20 deg, which corresponds to M_b =3.1 M_Jup and M_c = 7.4 M_Jup for the planetary companions.

Precise radial-velocity measurements with the HARPS spectrograph reveal the presence of two planets orbiting the solar-type star HD45364. The companion masses are 0.187 Mjup and 0.658 Mjup, with semi-major axes of 0.681 AU and 0.897 AU, and eccentricities of 0.168 and 0.097, respectively. A dynamical analysis of the system further shows a 3:2 mean motion resonance between the two planets, which prevents close encounters and ensures the stability of the system over 5 Gyr. This is the first time that such a resonant configuration has been observed for extra-solar planets, although there is an analogue in our Solar System formed by Neptune and Pluto. This singular planetary system may provide important constraints on planetary formation and migration scenarios.

The rotation of Mercury is presently captured in a 3/2 spin-orbit resonance with the orbital mean motion. The capture mechanism is well understood as the result of tidal interactions with the Sun combined with planetary perturbations. However, it is now almost certain that Mercury has a liquid core, which should induce a contribution of viscous friction at the core-mantle boundary to the spin evolution. This last effect greatly increases the chances of capture in all spin-orbit resonances, being 100% for the 2/1 resonance, and thus preventing the planet from evolving to the presently observed configuration. Here we show that for a given resonance, as the chaotic evolution of Mercury's orbit can drive its eccentricity to very low values during the planet's history, any previous capture can be destabilized whenever the eccentricity becomes lower than a critical value. In our numerical integrations of 1000 orbits of Mercury over 4 Gyr, the spin ends 99.8% of the time captured in a spin-orbit resonance, in particular in one of the following three configurations: 5/2 (22%), 2/1 (32%) and 3/2 (26%). Although the present 3/2 spin-orbit resonance is not the most probable outcome, we also show that the capture probability in this resonance can be increased up to 55% or 73%, if the eccentricity of Mercury in the past has descended below the critical values 0.025 or 0.005, respectively.

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.

The equilibrium rotation of tidally evolved "Earth-like" extra-solar planets is often assumed to be synchronous with their orbital mean motion. The same assumption persisted for Mercury and Venus until radar observations revealed their true spin rates. As many of these planets follow eccentric orbits and are believed to host dense atmospheres, we expect the equilibrium rotation to differ from the synchronous motion. Here we provide a general description of the allowed final equilibrium rotation states of these planets, and apply this to already discovered cases in which the mass is lower than twelve Earth-masses. At low obliquity and moderate eccentricity, it is shown that there are at most four distinct equilibrium possibilities, one of which can be retrograde. Because most presently known "Earth-like" planets present eccentric orbits, their equilibrium rotation is unlikely to be synchronous.

Over the past two years, the search for low-mass extrasolar planets has led to the detection of seven so-called 'hot Neptunes' or 'super-Earths' around Sun-like stars. These planets have masses 5-20 times larger than the Earth and are mainly found on close-in orbits with periods of 2-15 days. Here we report a system of three Neptune-mass planets with periods of 8.67, 31.6 and 197 days, orbiting the nearby star HD 69830. This star was already known to show an infrared excess possibly caused by an asteroid belt within 1 AU (the Sun-Earth distance). Simulations show that the system is in a dynamically stable configuration. Theoretical calculations favour a mainly rocky composition for both inner planets, while the outer planet probably has a significant gaseous envelope surrounding its rocky/icy core; the outer planet orbits within the habitable zone of this star.

Eccentricity or obliquity tides have been proposed as the missing energy source that may explain the anomalously large radius of some transiting ``hot Jupiters''. To maintain a non-zero and large obliquity, it was argued that the planets can be locked in a Cassini state, i.e. a resonance between spin and orbital precessions. We compute the tidal heating within ``inflated'' close-in giant planets with a non-zero eccentricity or obliquity. We further inspect whether the spin of a ``hot Jupiter'' could have been trapped and maintained in a Cassini state during its early despinning and migration. We estimate the capture probability in a spin-orbit resonance between $\sim$ 0.5 AU (a distance where tidal effects become significant) and 0.05 AU for a wide range of secular orbital frequencies and amplitudes of gravitational perturbations. Numerical simulations of the spin evolution are performed to explore the influence of tidal despinning and migration processes on the resonance stability. We find that tidal heating within a non-synchronous giant planet is about twice larger than previous estimates based on the hypothesis of synchronization. Chances of capture in a spin-orbit resonance are very good around 0.5 AU but they decrease dramatically with the semi-major axis. Furthermore, even if captured, both tidal despinning and migration processes cause the tidal torque to become large enough that the obliquity ultimately leaves the resonance and switches to near $0^{\circ}$. Locking a ``hot Jupiter'' in an isolated spin-orbit resonance is unlikely at 0.05 AU but could be possible at larger distances. Another mechanism is then required to maintain a large obliquity and create internal heating through obliquity tides

The mu Ara planetary system is rather complex: It contains two already known planets, mu Ara b with P=640 days and mu Ara c with P=9.64 days, and a third companion on a wide but still poorly defined orbit. Even with three planets in the system, the data points keep anomalously high dispersion around the fitted solution. The high residuals are only partially due to the strong p-mode oscillations of the host star. We have therefore studied in this paper the possible presence of a fourth planet in the system. During the past years we have carried out additional and extremely precise radial-velocity measurements with the HARPS spectrograph. We provide in this paper a full orbital solution of the planetary system around mu Ara. It turns out to be the second system known to harbor 4 planetary companions. Thanks to the new data points acquired with HARPS we can confirm the presence of mu Ara c at P=9.64 days, which produces a coherent RV signal over more than two years. The new orbital fit sets the mass of mu Ara c to 10.5 M_Earth. Furthermore, we present the discovery of mu Ara d, a new planet on an almost circular 310 days-period and with a mass of 0.52 M_Jup. Finally, we give completely new orbital parameters for the longest-period planet, mu Ara e. It is the first time that this companion is constrained by radial-velocity data into a dynamical stable orbit, which leaves no doubt about its planetary nature. (Abridged).

Long-term precise Doppler measurements with the CORALIE spectrograph reveal the presence of a second planet orbiting the solar-type star HD202206. The radial-velocity combined fit yields companion masses of m_2\sini = 17.4 M_Jup and 2.44 M_Jup, semi-major axes of a = 0.83 AU and 2.55 AU, and eccentricities of e = 0.43 and 0.27, respectively. A dynamical analysis of the system further shows a 5/1 mean motion resonance between the two planets. This system is of particular interest since the inner planet is within the brown-dwarf limits while the outer one is much less massive. Therefore, either the inner planet formed simultaneously in the protoplanetary disk as a superplanet, or the outer Jupiter-like planet formed in a circumbinary disk. We believe this singular planetary system will provide important constraints on planetary formation and migration scenarios.