The purpose of this work is to determine the location and stability of the Cassini states of a celestial body with an inviscid fluid core surrounded by a perfectly rigid mantle. Both situations where the rotation speed is either non-resonant or trapped in a p:1 spin-orbit resonance where p is a half integer are addressed. The rotation dynamics is described by the Poincar\'e-Hough model which assumes a simple motion of the core. The problem is written in a non-canonical Hamiltonian formalism. The secular evolution is obtained without any truncation in obliquity, eccentricity nor inclination. The condition for the body to be in a Cassini state is written as a set of two equations whose unknowns are the mantle obliquity and the tilt angle of the core spin-axis. Solving the system with Mercury's physical and orbital parameters leads to a maximum of 16 different equilibrium configurations, half of them being spectrally stable. In most of these solutions the core is highly tilted with respect to the mantle. The model is also applied to Io and the Moon.

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.

Eccentricity is a parameter of particular interest as it is an informative indicator of the past of planetary systems. It is however not always clear whether the eccentricity fitted on radial velocity data is real or if it is an artefact of an inappropriate modelling. In this work, we address this question in two steps: we first assume that the model used for inference is correct and present interesting features of classical estimators. Secondly, we study whether the eccentricity estimates are to be trusted when the data contain incorrectly modelled signals, such as missed planetary companions, non Gaussian noises, correlated noises with unknown covariance, etc. Our main conclusion is that data analysis via posterior distributions, with a model including a free error term gives reliable results provided two conditions. First, convergence of the numerical methods needs to be ascertained. Secondly, the noise power spectrum should not have a particularly strong peak at the semi period of the planet of interest. As a consequence, it is difficult to determine if the signal of an apparently eccentric planet might be due to another inner companion in 2:1 mean motion resonance. We study the use of Bayes factors to disentangle these cases. Finally, we suggest methods to check if there are hints of an incorrect model in the residuals. We show on simulated data the performance of our methods and comment on the eccentricities of Proxima b and 55 Cnc f.

We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators, while also reducing the effective step size whenever two planets encounter. We use a high order MVS scheme such that it is possible to integrate with large time steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e with a mutual separation of a fraction of the physical radius) while using the same time-step as in weakly perturbed problem such as the Solar System. We demonstrate the long term behaviour on systems of six super-Earths experiencing strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost we obtain much better energy conservation.

We address the expressions for the rates of the Keplerian orbital elements within a two-body problem perturbed by the tides in both partners. The formulae for these rates have appeared in the literature in various forms, at times with errors. We reconsider, from scratch, the derivation of these rates and arrive at the Lagrange-type equations which, in some details, differ from the corresponding equations obtained previously by Kaula (1964). We also write down detailed expressions for $da/dt$, $de/dt$ and $di/dt$, to order $e^4$. They differ from Kaula's expressions which contain a redundant factor of $M/(M+M^{\prime}),$ with $M$ and $M^{\prime}$ being the masses of the primary and the secondary. As Kaula was interested in the Earth-Moon system, this redundant factor was close to unity and was unimportant in his developments. This factor, however, must be reinstated when Kaula's theory is applied to a binary composed of partners of comparable masses. We have found that, while it is legitimate to simply sum the primary's and secondary's inputs in $da/dt$ or $de/dt$, this is not the case for $di/dt$. So our expression for $di/dt$ differs from that of Kaula in two regards. First, the contribution due to the dissipation in the secondary averages out when the apsidal precession is uniform. Second, we have obtained an additional term which emerges owing to the conservation of the angular momentum: a change in the inclination of the orbit causes a change of the primary's plane of equator.

Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a white dwarf and a near-spherical solid planet. This method determines the planet's inward and/or outward drift, and whether the planet will reach the Roche radius and be destroyed. We avoid constant tidal lag formulations and instead employ the self-consistent secular Darwin-Kaula expansions from Bou\'{e} & Efroimsky (2019), which feature an arbitrary frequency dependence on the quality functions. We adopt wide ranges of dynamic viscosities and spin rates for the planet in order to straddle many possible outcomes, and provide a foundation for the future study of individual systems with known or assumed rheologies. We find that: (i) massive Super-Earths are destroyed more readily than minor planets (such as the ones orbiting WD 1145+017 and SDSS J1228+1040), (ii) low-viscosity planets are destroyed more easily than high-viscosity planets, and (iii) the boundary between survival and destruction is likely to be fractal and chaotic.

Context: Seasonal variations and climate stability of a planet are very sensitive to the planet obliquity and its evolution. This is of particular interest for the emergence and sustainability of land-based life, but orbital and rotational parameters of exoplanets are still poorly constrained. Numerical explorations usually realised in this situation are thus in heavy contrast with the uncertain nature of the available data. Aims: We aim to provide an analytical formulation of the long-term spin-axis dynamics of exoplanets, linking it directly to physical and dynamical parameters, but still giving precise quantitative results if the parameters are well known. Together with bounds for the poorly constrained parameters of exoplanets, this analysis is designed to allow a quick and straightforward exploration of the spin-axis dynamics. Methods: The long-term orbital solution is decomposed in quasi-periodic series and the spin-axis Hamiltonian is expanded in powers of eccentricity and inclination. Chaotic zones are measured by the resonance overlap criterion. Bounds for the poorly known parameters of exoplanets are obtained from physical grounds (rotational breakup) and dynamical considerations (equipartition of AMD). Results: This method gives accurate results when the orbital evolution is well known. The chaotic zones for planets of the Solar System can be retrieved in details from simple analytical formulas. For less constrained planetary systems, the maximal extent of the chaotic regions can be computed, requiring only the mass, the semi-major axis and the eccentricity of the planets present in the system. Additionally, some estimated bounds of the precession constant allow to classify which observed exoplanets are necessarily out of major spin-orbit secular resonances (unless the precession rate is affected by the presence of massive satellites).

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.

In a two-planet system, due to Sundman (1912) inequality, a topological boundary can forbid close encounters between the two planets for infinite time. A system is said Hill stable if it verifies this topological condition. Hill stability is widely used in the study of extra solar planets dynamics. However people often use the coplanar and circular orbits approximation. In this paper, we explain how the Hill stability can be understood in the framework of Angular Momentum Deficit (AMD). In the secular approximation, the AMD allows to discriminate between a priori stable systems and systems for which a more in depth dynamical analysis is required. We show that the general Hill stability criterion can be expressed as a function of only the semi major axes, the masses and the total AMD of the system. The proposed criterion is only expanded in the planets-to-star mass ratio $\epsilon$ and not in the semi-major axis ratio, in eccentricities nor in the mutual inclination. Moreover the expansion in $\epsilon$ remains excellent up to values of about $10^{-3}$ even for two planets with very different mass values. We performed numerical simulations in order to highlight the sharp change of behaviour between Hill stable and Hill unstable systems. We show that Hill stable systems tend to be very regular whereas Hill unstable ones often lead to rapid planet collisions. We also remind that Hill stability does not protect from the ejection of the outer planet.

We revisit the rotation dynamics of a rigid satellite with either a liquid core or a global sub-surface ocean. In both problems, the flow of the fluid component is assumed inviscid. The study of a hollow satellite with a liquid core is based on the Poincare-Hough model which provides exact equations of motion. We introduce an approximation when the ellipticity of the cavity is low. This simplification allows to model both types of satellite in the same manner. The analysis of their rotation is done in a non-canonical Hamiltonian formalism closely related to Poincare's "forme nouvelle des equations de la mecanique". In the case of a satellite with a global ocean, we obtain a seven-degree of freedom system. Six of them account for the motion of the two rigid components, and the last one is associated with the fluid layer. We apply our model to Titan for which the origin of the obliquity is still a debated question. We show that the observed value is compatible with Titan slightly departing from the hydrostatic equilibrium and being in a Cassini equilibrium state.

The AMD-stability criterion allows to discriminate between a-priori stable planetary systems and systems for which the stability is not granted and needs further investigations. AMD-stability is based on the conservation of the Angular Momentum Deficit (AMD) in the averaged system at all orders of averaging. While the AMD criterion is rigorous, the conservation of the AMD is only granted in absence of mean-motion resonances (MMR). Here we extend the AMD-stability criterion to take into account mean-motion resonances, and more specifically the overlap of first order MMR. If the MMR islands overlap, the system will experience generalized chaos leading to instability. The Hamiltonian of two massive planets on coplanar quasi-circular orbits can be reduced to an integrable one degree of freedom problem for period ratios close to a first order MMR. We use the reduced Hamiltonian to derive a new overlap criterion for first order MMR. This stability criterion unifies the previous criteria proposed in the literature and admits the criteria obtained for initially circular and eccentric orbits as limit cases. We then improve the definition of AMD-stability to take into account the short term chaos generated by MMR overlap. We analyze the outcome of this improved definition of AMD-stability on selected multi-planet systems from the Extrasolar Planets Encyclopeadia.

The velocity distribution function (VDF) of dark matter (DM) halos in $\Lambda$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $\Lambda$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless spherical velocity components (after subtraction of streaming motions), normalized by the respective velocity dispersions. We test our VDF on 90 cluster-mass halos of a dissipationless cosmological simulation. Beyond the virial radius, $r_{\rm vir}$, our model VDF adequately reproduces that measured in the simulated halos, but no $q$-Gaussian model can adequately represent the VDF within $r_{\rm vir}$, as the speed distribution function is then flatter-topped than any $q$-Gaussian can allow. Nevertheless, our VDF fits significantly better the simulations than the commonly used Maxwellian (Gaussian) distribution, at virtually all radii within $5\,r_{\rm vir}$. Within 0.4 (1) $r_{\rm vir}$, the non-Gaussianity index $q$ is (roughly) linearly related to the slope of the density profile and also to the velocity anisotropy profile. We provide a parametrization of the modulation of $q$ with radius for both the median fits and the fit of the stacked halo. At radii of a few percent of $r_{\rm vir}$, corresponding to the Solar position in the Milky Way, our best-fit VDF, although fitting better the simulations than the Gaussian one, overproduces significantly the fraction of high velocity objects, indicating that one should not blindly use these $q$-Gaussian fits to make predictions on the direct detection rate of DM particles.

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.

This work presents an elegant formalism to model the evolution of the full two rigid body problem. The equations of motion, given in a Cartesian coordinate system, are expressed in terms of spherical harmonics and Wigner D-matrices. The algorithm benefits from the numerous recurrence relations satisfied by these functions allowing a fast evaluation of the mutual potential. Moreover, forces and torques are straightforwardly obtained by application of ladder operators taken from the angular momentum theory and commonly used in quantum mechanics. A numerical implementation of this algorithm is made. Tests show that the present code is significantly faster than those currently available in literature.

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

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.

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.

The stellar spin orientation relative to the orbital planes of multiplanet systems are becoming accessible to observations. For example, 55 Cancri is a system composed of 5 planets orbiting a member of a stellar binary for which a projected obliquity of 72+-12 deg relative to the orbit of the innermost planet has been reported (Bourrier & Hebrard 2014). This large obliquity has been attributed to the perturbation induced by the binary. Here we describe the secular evolution of similar systems and we discuss the case of the 55 Cancri system more deeply. We provide two different orbital configurations compatible with the currently available observations.

In this work we present the solution of the stellar spot problem using the Kelvin-Stokes theorem. Our result is applicable for any given location and dimension of the spots on the stellar surface. We present explicitely the result up to the second degree in the limb darkening law. This technique can be used to calculate very efficiently mutual photometric effects produced by eclipsing bodies occulting stellar spots and to construct complex spot shapes.

The non-resonant secular dynamics of compact planetary systems are modeled by a perturbing function which is usually expanded in eccentricity and absolute inclination with respect to the invariant plane. Here, the expressions are given in a vectorial form which naturally leads to an expansion in eccentricity and mutual inclination. The two approaches are equivalent in most cases, but the vectorial one is specially designed for those where a quasi-coplanar system tilts as a whole by a large amount. Moreover, the vectorial expressions of the Hamiltonian and of the equations of motion are slightly simpler than those given in terms of the usual elliptical elements. We also provide the secular perturbing function in vectorial form expanded in semimajor axis ratio allowing for arbitrary eccentricities and inclinations. The interaction between the equatorial bulge of a central star and its planets is also provided, as is the relativistic periapse precession of any planet induced by the central star. We illustrate the use of this representation for following the secular oscillations of the terrestrial planets of the solar system, and for Kozai cycles as may take place in exoplanetary systems.

The stellar spin orientation relative to the orbital planes of multiplanet systems are becoming accessible to observations. Here, we analyze and classify different types of spin-orbit evolution in compact multiplanet systems perturbed by an inclined outer companion. Our study is based on classical secular theory, using a vectorial approach developed in a separate paper. When planet-planet perturbations are truncated at the second order in eccentricity and mutual inclination, and the planet-companion perturbations are developed at the quadrupole order, the problem becomes integrable. The motion is composed of a uniform precession of the whole system around the total angular momentum, and in the rotating frame, the evolution is periodic. Here, we focus on the relative motion associated to the oscillations of the inclination between the planet system and the outer orbit, and of the obliquities of the star with respect to the two orbital planes. The solution is obtained using a powerful geometric method. With this technique, we identify four different regimes characterized by the nutation amplitude of the stellar spin-axis relative to the orbital plane of the planets. In particular, the obliquity of the star reaches its maximum when the system is in the Cassini regime where planets have more angular momentum than the star, and where the precession rate of the star is similar to that of the planets induced by the companion. In that case, spin-orbit oscillations exceed twice the inclination between the planets and the companion. Even if mutual inclination is only ~ 20 deg, this resonant case can cause the spin-orbit angle to oscillate between perfectly aligned and retrograde values.

In general, in the studies of transit light-curves and the Rossiter-McLaughlin (RM), the contribution of the planet's gravitational microlensing is neglected. Theoretical studies, have, however shown that the planet's microlensing can affect the transit light-curve and in some extreme cases cause the transit depth to vanish. In this letter, we present the results of our quantitative analysis of microlening on the RM effect. Results indicate that for massive planets in on long period orbits, the planet's microlensing will have considerable contribution to the star's RV measurements. We present the details of our study, and discuss our analysis and results.

Stellar activity features such as spots can create complications in determining planetary parameters through spectroscopic and photometric observations. The overlap of a transiting planet and a stellar spot, for instance, can produce anomalies in the transit light curves that may lead to inaccurate estimation of the transit duration, depth and timing. For instance such inaccuracies can affect the precise derivation of the planet's radius. In this paper, we present the results of a quantitative study on the effects of stellar spots on high precision transit light curves. We show that spot anomalies can lead to the estimate of a planet radius that is 4% smaller than the real value. The effects on the transit duration can also be of the order of 4%, longer or shorter. Depending on the size and distribution of spots, anomalies can also produce transit timing variations with significant amplitudes. For instance, TTVs with signal amplitudes of 200 seconds can be produced when the spot is completely dark and as large as the largest Sun spot. Our study also indicates that the minimum size of a stellar spot with detectable affects on the high precision transit light curve is around 0.03 time the stellar radius for typical Kepler telescope precision. We also show that the strategy of including more free parameters (such as transit depth and duration) in the fitting procedure to measure the transit time of each individual transit will not produce accurate results in the case of active stars.

Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic systems. A new fast method, called MAMPOSSt, which performs a maximum likelihood fit of the distribution of observed tracers in projected phase space, is developed and thoroughly tested. MAMPOSSt assumes a shape for the gravitational potential, but instead of postulating a shape for the distribution function in terms of energy and angular momentum, or supposing Gaussian line-of-sight velocity distributions, MAMPOSSt assumes a VA profile and a shape for the 3D velocity distribution, here Gaussian. MAMPOSSt requires no binning, differentiation, nor extrapolation of the observables. Tests on cluster-mass haloes from LambdaCDM cosmological simulations show that, with 500 tracers, MAMPOSSt is able to jointly recover the virial radius, tracer scale radius, dark matter scale radius and outer or constant VA with small bias (<10% on scale radii and <2% on the two other quantities) and inefficiencies of 10%, 27%, 48% and 20%, respectively. MAMPOSSt does not perform better when some parameters are frozen, and even worse when the virial radius is set to its true value, which appears to be the consequence of halo triaxiality. The accuracy of MAMPOSSt depends weakly on the adopted interloper removal scheme, including an efficient iterative Bayesian scheme that we introduce here, which can directly obtain the virial radius with as good precision as MAMPOSSt. Our tests show that MAMPOSSt with Gaussian 3D velocities is very competitive with, and up to 1000x faster than other methods. Hence, MAMPOSSt is a very powerful and rapid tool for the mass and anisotropy modeling of systems such as clusters and groups of galaxies, elliptical and dwarf spheroidal galaxies.

Context: With the announced arrival of instruments such as ESPRESSO one can expect that several systematic noise sources on the measurement of precise radial velocity will become the limiting factor instead of photon noise. A stellar companion within the fiber is such a possible noise source. Aims: With this work we aim at characterizing the impact of a stellar companion within the fiber to radial velocity measurements made by fiber-fed spectrographs. We consider the contaminant star either to be part of a binary system whose primary star is the target star, or as a background/foreground star. Methods: To carry out our study, we used HARPS spectra, co-added the target with contaminant spectra, and then compared the resulting radial velocity with that obtained from the original target spectrum. We repeated this procedure and used different tunable knobs to reproduce the previously mentioned scenarios. Results: We find that the impact on the radial velocity calculation is a function of the difference between individual radial velocities, of the difference between target and contaminant magnitude, and also of their spectral types. For the worst-case scenario in which both target and contaminant star are well centered on the fiber, the maximum contamination for a G or K star may be higher than 10 cm/s, on average, if the difference between target and contaminant magnitude is $\Delta m$ < 10, and higher than 1 m/s if $\Delta m$ < 8. If the target star is of spectral type M, $\Delta m$ < 8 produces the same contamination of 10 cm/s, and a contamination may be higher than 1 m/s

The Rossiter-McLaughlin (hereafter RM) effect is a key tool for measuring the projected spin-orbit angle between stellar spin axes and orbits of transiting planets. However, the measured radial velocity (RV) anomalies produced by this effect are not intrinsic and depend on both instrumental resolution and data reduction routines. Using inappropriate formulas to model the RM effect introduces biases, at least in the projected velocity Vsin(i) compared to the spectroscopic value. Currently, only the iodine cell technique has been modeled, which corresponds to observations done by, e.g., the HIRES spectrograph of the Keck telescope. In this paper, we provide a simple expression of the RM effect specially designed to model observations done by the Gaussian fit of a cross-correlation function (CCF) as in the routines performed by the HARPS team. We derived also a new analytical formulation of the RV anomaly associated to the iodine cell technique. For both formulas, we modeled the subplanet mean velocity v_p and dispersion beta_p accurately taking the rotational broadening on the subplanet profile into account. We compare our formulas adapted to the CCF technique with simulated data generated with the numerical software SOAP-T and find good agreement up to Vsin(i) < 20 km/s. In contrast, the analytical models simulating the two different observation techniques can disagree by about 10 sigma in Vsin(i) for large spin-orbit misalignments. It is thus important to apply the adapted model when fitting data.

We present an improved version of SOAP (Boisse et al. 2012) named "SOAP-T", which can generate the radial velocity variations and light-curves for systems consisting of a rotating spotted star with a transiting planet. This tool can be used to study the anomalies inside transit light-curves and the Rossiter-McLaughlin effect, to better constrain the orbital configuration and properties of planetary systems and active zones of their host stars. Tests of the code are presented to illustrate its performance and to validate its capability when compared with analytical models and real data. Finally, we apply SOAP-T to the active star, HAT-P-11, observed by the NASA Kepler space telescope and use this system to discuss the capability of this tool in analyzing light-curves for the cases where the transiting planet overlaps with the star's spots.

Mercury's eccentricity is chaotic and can increase so much that collisions with Venus or the Sun become possible (Laskar, 1989, 1990, 1994, 2008, Batygin & Laughlin, 2008, Laskar & Gastineau, 2009). This chaotic behavior results from an intricate network of secular resonances, but in this paper, we show that a simple integrable model with only one degree of freedom is actually able to reproduce the large variations in Mercury's eccentricity, with the correct amplitude and timescale. We show that this behavior occurs in the vicinity of the separatrices of the resonance g1-g5 between the precession frequencies of Mercury and Jupiter. However, the main contribution does not come from the direct interaction between these two planets. It is due to the excitation of Venus' orbit at Jupiter's precession frequency g5. We use a multipolar model that is not expanded with respect to Mercury's eccentricity, but because of the proximity of Mercury and Venus, the Hamiltonian is expanded up to order 20 and more in the ratio of semimajor axis. When the effects of Venus' inclination are added, the system becomes nonintegrable and a chaotic zone appears in the vicinity of the separatrices. In that case, Mercury's eccentricity can chaotically switch between two regimes characterized by either low-amplitude circulations or high-amplitude librations.

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.

Aims. We study a subset of the planetary population characterized both by HARPS and Kepler surveys. We compare the statistical properties of planets in systems with m.sin i >5-10 M_Earth and R>2 R_Earth. If we assume that the underlying population has the same characteristics, the different detection sensitivity to the orbital inclination relative to the line of sight allows us to probe the planets' mutual inclination. Methods. We considered the frequency of systems with one, two and three planets as dictated by HARPS data. We used Kepler's planetary period and host mass and radii distributions (corrected from detection bias) to model planetary systems in a simple yet physically plausible way. We then varied the mutual inclination between planets in a system according to different prescriptions (completely aligned, Rayleigh distributions and isotropic) and compared the transit frequencies with one, two or three planets with those measured by Kepler. Results. The results show that the two datasets are compatible, a remarkable result especially because there are no tunable knobs other than the assumed inclination distribution. For m.sin i cutoffs of 7-10 M_Earth, which are those expected to correspond to the radius cutoff of 2 R_Earth, we conclude that the results are better described by a Rayleigh distribution with mode of 1 deg or smaller. We show that the best-fit scenario only becomes a Rayleigh distribution with mode of 5 deg if we assume a rather extreme mass-radius relationship for the planetary population. Conclusions. These results have important consequences for our understanding of the role of several proposed formation and evolution mechanisms. They confirm that planets are likely to have been formed in a disk and show that most planetary systems evolve quietly without strong angular momentum exchanges (abridged).

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

We present a model-independent technique for calculating the time of mid-transits. This technique, named "barycenter method", uses the light-curve's symmetry to determine the transit timing by calculating the transit light-curve barycenter. Unlike the other methods of calculating mid-transit timing, this technique does not depend on the parameters of the system and central star. We demonstrate the capabilities of the barycenter method by applying this technique to some known transiting systems including several \emph{Kepler} confirmed planets. Results indicate that for complete and symmetric transit lightcurves, the barycenter method achieves the same precision as other techniques, but with fewer assumptions and much faster. Among the transiting systems studied with the barycenter method, we focus in particular on LHS 6343C, a brown dwarf that transits a member of an M+M binary system, LHS 6343AB. We present the results of our analysis, which can be used to set an upper limit on the period and mass of a possible second small perturber.

The transit timing variation (TTV) method allows the detection of non-transiting planets through their gravitational perturbations. Since TTVs are strongly enhanced in systems close to mean-motion resonances (MMR), even a low mass planet can produce an observable signal. This technique has thus been proposed to detect terrestrial planets. In this letter, we analyse TTV signals for systems in or close to MMR in order to illustrate the difficulties arising in the determination of planetary parameters. TTVs are computed numerically with an n-body integrator for a variety of systems close to MMR. The main features of these TTVs are also derived analytically. Systems deeply inside MMR do not produce particularly strong TTVs, while those close to MMR generate quasiperiodic TTVs characterised by a dominant long period term and a low amplitude remainder. If the remainder is too weak to be detected, then the signal is strongly degenerate and this prevents the determination of the planetary parameters. Even though an Earth mass planet can be detected by the TTV method if it is close to a MMR, it may not be possible to assert that this planet is actually an Earth mass planet. On the other hand, if the system is right in the center of a MMR, the high amplitude oscillation of the TTV signal vanishes and the detection of the perturber becomes as difficult as it is far from MMR.

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

Context: Several studies have so far placed useful constraints on planetary atmospheric properties using transmission spectrsocopy, and in the case of HD209458b even the radial velocity of the planet during the transit event has been reconstructed opening a new range of possibilities. AIMS. In this contribution we highlight the importance to account for the orbital eccentricity and longitude of periastron of the planetary orbit to accurately interpret the measured planetary radial velocity during the transit. Methods: We calculate the radial velocity of a transiting planet in an eccentric orbit. Given the larger orbital speed of planets with respect to their stellar companions even small eccentricities can result in detectable blue or redshift radial velocity offsets during the transit with respect to the systemic velocity, the exact value depending also on the longitude of the periastron of the planetary orbit. For an hot-jupiter planet, an eccentricity of only e=0.01 can produce a radial velocity offset of the order of the km/s. Conclusions: We propose an alternative interpretation of the recently claimed radial velocity blueshift (~2 km/s) of the planetary spectral lines of HD209458b which implies that the orbit of this system is not exactly circular. In this case, the longitude of the periastron of the stellar orbit is most likely confined in the first quadrant (and that one of the planet in the third quadrant). We highlight that transmission spectroscopy allows not only to study the compositional properties of planetary atmospheres, but also to refine their orbital parameters and that any conclusion regarding the presence of windflows on planetary surfaces coming from transmission spectroscopy measurements requires precise known orbital parameters from RV.

Since the original work of Hansen and Tisserand in the XIXth century, there have been many variations in the analytical expansion of the three-body disturbing function in series of the semi-major axis ratio. With the increasing number of planetary systems of large eccentricity, these expansions are even more interesting as they allow us to obtain for the secular systems finite expressions that are valid for all eccentricities and inclinations. We revisited the derivation of the disturbing function in Legendre polynomial, with a special focus on the secular system. We provide here expressions of the disturbing function for the planar and spatial case at any order with respect to the ratio of the semi-major axes. Moreover, for orders in the ratio of semi-major axis up to ten in the planar case and five in the spatial case, we provide explicit expansions of the secular system, and simple algorithms with minimal computation to extend this to higher order, as well as the algorithms for the computation of non secular terms.

The origin of the high inclination of Uranus' spin-axis (Uranus' obliquity) is one of the great unanswered questions about the Solar system. Giant planets are believed to form with nearly zero obliquity, and it has been shown that the present behaviour of Uranus' spin is essentially stable. Several attempts were made in order to solve this problem. Here we report numerical simulations showing that Uranus' axis can be tilted during the planetary migration, without the need of a giant impact, provided that the planet had an additional satellite and a temporary large inclination. This might have happened during the giant planet instability phase described in the Nice model. In our scenario, the satellite is ejected after the tilt by a close encounter at the end of the migration. This model can both explain Uranus' large obliquity and bring new constraints on the planet orbital evolution.

Traditionally, the mass / velocity anisotropy degeneracy (MAD) inherent in the spherical, stationary, non-streaming Jeans equation has been handled by assuming a mass profile and fitting models to the observed kinematical data. Here, the opposite approach is considered: the equation of anisotropic kinematic projection is inverted for known arbitrary anisotropy to yield the space radial velocity dispersion profile in terms of an integral involving the radial profiles of anisotropy and isotropic dynamical pressure. Then, through the Jeans equation, the mass profile is derived in terms of double integrals of observable quantities. Single integral formulas for both deprojection and mass inversion are provided for several simple anisotropy models (isotropic, radial, circular, general constant, Osipkov-Merritt, Mamon-Lokas and Diemand-Moore-Stadel). Tests of the mass inversion on NFW models with these anisotropy models yield accurate results in the case of perfect observational data, and typically better than 70% (in 4 cases out of 5) accurate mass profiles for the sampling errors expected from current observational data on clusters of galaxies. For the NFW model with mildly increasing radial anisotropy, the mass is found to be insensitive to the adopted anisotropy profile at 7 scale radii and to the adopted anisotropy radius at 3 scale radii. This anisotropic mass inversion method is a useful complementary tool to analyze the mass and anisotropy profiles of spherical systems. It provides the practical means to lift the MAD in quasi-spherical systems such as globular clusters, round dwarf spheroidal and elliptical galaxies, as well as groups and clusters of galaxies, when the anisotropy of the tracer is expected to be linearly related to the slope of its density.

In this Letter, we give new constraints on planet migration. They were obtained under the assumption that Saturn's current obliquity is due to a capture in resonance with Neptune's ascending node. If planet migration is too fast, then Saturn crosses the resonance without being captured and it keeps a small obliquity. This scenario thus gives a lower limit on the migration time scale tau. We found that this boundary depends strongly on Neptune's initial inclination. For two different migration types, we found that tau should be at least greater than 7 Myr. This limit increases rapidly as Neptune's initial inclination decreases from 10 to 1 degree. We also give an algorithm to know if Saturn can be tilted for any migration law.

We compare the distribution of optically and Halpha (Ha) selected galaxies in the Southern half of the nearby Abell 85 (A85) cluster with the recently discovered X-ray filament (XRF). We search for galaxies where star formation (SF) may have been triggered by interactions with intracluster gas or tidal pressure due to the cluster potential when entering the cluster. Our analysis is based on images obtained with CFHT MegaPrime/MegaCam (1x1 deg2 field) in four bands (ugri) and ESO 2.2mWFI (38'x36' field) in a narrow band filter corresponding to the redshifted Halpha (Ha) line and in a broad R-band filter. The LFs are estimated by statistically subtracting a reference field. Background contamination is minimized by cutting out galaxies redder than the observed red sequence in the g-i vs. i colour-magnitude diagram. The galaxy distribution shows a significantly flattened cluster, whose principal axis is slightly offset from the XRF. The analysis of the broad band LFs shows that the filament region is well populated. The filament is also independently detected as a gravitationally bound structure by the Serna & Gerbal hierarchical method. 101 galaxies are detected in Ha, among which 23 have spectroscopic redshifts in the cluster, 2 have spectroscopic redshifts higher than the cluster and 58 have photometric redshifts that tend to indicate that they are background objects.The 23 galaxies with spectroscopic redshifts in the cluster are mostly concentrated in the South part of the cluster and along the filament. We find a number of galaxies showing evidence for SF in the XRF, and all our results are consistent with the previous hypothesis that the XRF in A85 is a gravitationally bound structure made of groups falling on to the main cluster.

We investigate the LF in the very relaxed cluster Abell 496. Our analysis is based on deep images obtained at CFHT with MegaPrime/MegaCam in four bands (ugri) covering a 1x1 deg2 region, which is centered on the cluster Abell 496 and extends to near its virial radius. The LFs are estimated by statistically subtracting a reference field taken as the mean of the 4 Deep fields of the CFHTLS survey. Background contamination is minimized by cutting out galaxies redder than the observed Red Sequence in the g-i versus i colour-magnitude diagram. In Abell 496, the global LFs show a faint-end slope alpha=-1.55+/-0.06 and vary little with observing band. Without colour cuts, the LFs are much noisier but not significantly steeper. The faint-end slopes show a statistically significant steepening from alpha=-1.4+/-0.1 in the central region (extending to half a virial radius) to -1.8+/-0.1 in the Southern envelope of the cluster. Cosmic variance and uncertain star-galaxy separation are our main limiting factors in measuring the faint-end of the LFs. The large-scale environment of Abell 496, probed with the fairly complete 6dFGS catalogue, shows a statistically significant 36 Mpc long filament at PA=137 deg, but we do not find an enhanced LF along this axis. Our LFs do not display the large number of dwarf galaxies (alpha ~ -2) inferred by several authors, whose analyses may suffer from field contamination caused by inexistent or inadequate colour cuts. Alternatively, different clusters may have different faint-end slopes, but this is hard to reconcile with the wide range of slopes found for given clusters and for wide sets of clusters.