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.

Correlated noise affects most astronomical datasets and to neglect accounting for it can lead to spurious signal detections, especially in low signal-to-noise conditions, which is often the context in which new discoveries are pursued. For instance, in the realm of exoplanet detection with radial velocity time series, stellar variability can induce false detections. However, a white noise approximation is often used because accounting for correlated noise when analyzing data implies a more complex analysis. Moreover, the computational cost can be prohibitive as it typically scales as the cube of the dataset size. For some restricted classes of correlated noise models, there are specific algorithms that can be used to help bring down the computational cost. This improvement in speed is particularly useful in the context of Gaussian process regression, however, it comes at the expense of the generality of the noise model. Here, we present the S+LEAF noise model, which allows us to account for a large class of correlated noises with a linear scaling of the computational cost with respect to the size of the dataset. The S+LEAF model includes, in particular, mixtures of quasiperiodic kernels and calibration noise. This efficient modeling is made possible by a sparse representation of the covariance matrix of the noise and the use of dedicated algorithms for matrix inversion, solving, determinant computation, etc. We applied the S+LEAF model to reanalyze the HARPS radial velocity time series of HD 136352. We illustrate the flexibility of the S+LEAF model in handling various sources of noise. We demonstrate the importance of taking correlated noise into account, and especially calibration noise, to correctly assess the significance of detected signals. We provide an open-source implementation of the S+LEAF model, available at https://gitlab.unige.ch/jean-baptiste.delisle/spleaf.

Since 2011, the SOPHIE spectrograph has been used to search for Neptunes and super-Earths in the Northern Hemisphere. As part of this observational program, 290 radial velocity measurements of the 6.4 V magnitude star HD 158259 were obtained. Additionally, TESS photometric measurements of this target are available. We present an analysis of the SOPHIE data and compare our results with the output of the TESS pipeline. The radial velocity data, ancillary spectroscopic indices, and ground-based photometric measurements were analyzed with classical and $\ell_1$ periodograms. The stellar activity was modeled as a correlated Gaussian noise and its impact on the planet detection was measured with a new technique. The SOPHIE data support the detection of five planets, each with $m \sin i \approx 6 M_\oplus$, orbiting HD 158259 in 3.4, 5.2, 7.9, 12, and 17.4 days. Though a planetary origin is strongly favored, the 17.4 d signal is classified as a planet candidate due to a slightly lower statistical significance and to its proximity to the expected stellar rotation period. The data also present low frequency variations, most likely originating from a magnetic cycle and instrument systematics. Furthermore, the TESS pipeline reports a significant signal at 2.17 days corresponding to a planet of radius $\approx 1.2 R_\oplus$. A compatible signal is seen in the radial velocities, which confirms the detection of an additional planet and yields a $\approx 2 M_\oplus$ mass estimate. We find a system of five planets and a strong candidate near a 3:2 mean motion resonance chain orbiting HD 158259. The planets are found to be outside of the two and three body resonances.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

Context. HD13724 is a nearby solar-type star at 43.48 $\pm$ 0.06 pc hosting a long-period low-mass brown dwarf detected with the CORALIE echelle spectrograph as part of the historical CORALIE radial-velocity search for extra-solar planets. The companion has a minimum mass of $26.77^{+4.4}_{-2.2} M_{\mathrm{Jup}}$ and an expected semi-major axis of $\sim$ 240 mas making it a suitable target for further characterisation with high-contrast imaging, in particular to measure its inclination, mass, and spectrum and thus establish its substellar nature. Aims. Using high-contrast imaging with the SPHERE instrument on the Very Large Telescope (VLT), we are able to directly image a brown dwarf companion to HD13724 and obtain a low-resolution spectrum. Methods. We combine the radial-velocity measurements of CORALIE and HARPS taken over two decades and high contrast imaging from SPHERE to obtain a dynamical mass estimate. From the SPHERE data we obtain a low resolution spectrum of the companion from Y to J band, as well as photometric measurements from IRDIS in the J, H and K bands. Results. Using high-contrast imaging with the SPHERE instrument at the VLT, we report the first images of a brown dwarf companion to the host star HD13724. It has an angular separation of 175.6 $\pm$ 4.5 mas and H-band contrast of $10.61\pm0.16$ mag and, using the age estimate of the star to be $\sim$1 Gyr, gives an isochronal mass estimate of $\sim$44 $M_{\mathrm{Jup}}$. By combining radial-velocity and imaging data we also obtain a dynamical mass of $50.5^{+3.3}_{-3.5} M_{\mathrm{Jup}}$. Through fitting an atmospheric model, we estimate a surface gravity of $\log g = 5.5$ and an effective temperature of 1000K. A comparison of its spectrum with observed T dwarfs estimates a spectral type of T4 or T4.5, with a T4 object providing the best fit.

Periodograms are common tools used to search for periodic signals in unevenly spaced time series. The significance of periodogram peaks is often assessed using false alarm probability (FAP), which in most studies assumes uncorrelated noise and is computed using numerical methods such as bootstrapping or Monte Carlo. These methods have a high computational cost, especially for low FAP levels, which are of most interest. We present an analytical estimate of the FAP of the periodogram in the presence of correlated noise, which is fundamental to analyze astronomical time series correctly. The analytical estimate that we derive provides a very good approximation of the FAP at a much lower cost than numerical methods. We validate our analytical approach by comparing it with Monte Carlo simulations. Finally, we discuss the sensitivity of the method to different assumptions in the modeling of the noise.

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.

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.

Context. Since 1998, a planet-search around main sequence stars within 50~pc in the southern hemisphere has been carried out with the CORALIE spectrograph at La Silla Observatory. Aims. With an observing time span of more than 20 years, the CORALIE survey is able to detect long term trends in data with masses and separations large enough to select ideal targets for direct imaging. Detecting these giant companion candidates will allow us to start bridging the gap between radial velocity detected exoplanets and directly imaged planets and brown dwarfs. Methods. Long-term precise Doppler measurements with the CORALIE spectrograph reveal radial velocity signatures of massive planetary companions and brown dwarfs on long-period orbits. Results. In this paper we report the discovery of new companions orbiting HD~181234, HD~13724, HD~25015, HD~92987 and HD~50499. We also report updated orbital parameters for HD~50499b, HD~92788b and HD~98649b. In addition, we confirm the recent detection of HD~92788c. The newly reported companions span a period range of 15.6 to 40.4 years and a mass domain of 2.93 to 26.77 $M_{\mathrm{Jup}}$, the latter of which straddles the nominal boundary between planets and brown dwarfs. Conclusion. We have reported the detection of five new companions and updated parameters of four known extrasolar planets. We identify at least some of these companions to be promising candidates for imaging and further characterisation.

We report the discovery of an additional substellar companion in the CoRoT-20 system based on six years of HARPS and SOPHIE radial velocity follow-up. CoRoT-20 c has a minimum mass of 17 $\pm$ 1 $M_{Jup}$ and it orbits the host star in 4.59$\pm 0.05$ years, with an orbital eccentricity of 0.60 $\pm$ 0.03. This is the first identified system with an eccentric hot Jupiter and an eccentric massive companion. The discovery of the latter might be an indication of the migration mechanism of the hot Jupiter, via Lidov-Kozai effect. We explore the parameter space to determine which configurations would trigger this type of interactions.

We report the discovery of four super-Earth planets around HD 215152, with orbital periods of 5.76, 7.28, 10.86, and 25.2 d, and minimum masses of 1.8, 1.7, 2.8, and 2.9 M_Earth respectively. This discovery is based on 373 high-quality radial velocity measurements taken by HARPS over 13 years. Given the low masses of the planets, the signal-to-noise ratio is not sufficient to constrain the planet eccentricities. However, a preliminary dynamical analysis suggests that eccentricities should be typically lower than about 0.03 for the system to remain stable. With two pairs of planets with a period ratio lower than 1.5, with short orbital periods, low masses, and low eccentricities, HD 215152 is similar to the very compact multi-planet systems found by Kepler, which is very rare in radial-velocity surveys. This discovery proves that these systems can be reached with the radial-velocity technique, but characterizing them requires a huge amount of observations.

Using high-contrast imaging with the SPHERE instrument at the VLT, we report the first images of a cold brown dwarf companion to the exoplanet host star HD4113A. The brown dwarf HD4113C is part of a complex dynamical system consisting of a giant planet, stellar host and a known wide M-dwarf companion. Its separation of $535\pm3$mas and H-band contrast of $13.35\pm0.10$mag correspond to a projected separation of 22AU and an isochronal mass estimate of $36\pm5$M$_J$ based on COND models. The companion shows strong methane absorption, and through atmospheric model fitting we estimate a surface gravity of $\log g$=5 and an effective temperature of ~500-600K. A comparison of its spectrum with observed T dwarfs indicates a late-T spectral type, with a T9 object providing the best match. By combining the observed astrometry from the imaging data with 27 years of radial velocities, we use orbital fitting to constrain its orbital and physical parameters, as well as update those of the planet HD4113Ab, discovered by previous radial velocity measurements. The data suggest a dynamical mass of $66\pm5$M$_J$ and moderate eccentricity of $0.44\pm0.08$ for the brown dwarf. This mass estimate appears to conflict with the isochronal estimate and that of similar objects, which may be caused by the newly detected object being an unresolved binary brown dwarf system or the presence of an additional object in the system. Through dynamical simulations we show that the planet may undergo strong Lidov-Kozai cycles, raising the possibility that it formed on a quasi-circular orbit and gained its currently observed high eccentricity through interactions with the brown dwarf. Follow-up observations combining radial velocities, direct imaging and Gaia astrometry will be crucial to precisely constrain the dynamical mass of the brown dwarf and allow for in-depth comparison with evolutionary and atmospheric models.

Exoplanet surveys have shown that systems with multiple low-mass planets on compact orbits are common. Except for few cases, however, the masses of these planets are generally unknown. At the very end of the main sequence, host stars have the lowest mass and hence offer the largest reflect motion for a given planet. In that context, we monitored the low-mass (0.13Msun) M dwarf YZ Cet (GJ 54.1, HIP 5643) intensively and obtained both radial velocities and stellar-activity indicators derived from both spectroscopy and photometry. We find strong evidence that it is orbited by at least three planets in compact orbits (P=1.97, 3.06, 4.66 days), with the inner two near a 2:3 mean-motion resonance. The minimum masses are comparable to that of Earth (Msini=0.75+-0.13, 0.98+-0.14, and 1.14+-0.17 Mearth) and also the lowest masses measured by radial velocity so far. We note the possibility for an even lower-mass, fourth planet with Msini=0.472+-0.096 Mearth at P=1.04 days. An n-body dynamical model is used to put further constraints on the system parameters. At 3.6 parsecs, YZ Cet is the nearest multi-planet system detected to date.

Resonant chains are groups of planets for which each pair is in resonance, with an orbital period ratio locked at a rational value (2/1, 3/2, etc.). Such chains naturally form as a result of convergent migration of the planets in the proto-planetary disk. In this article, I present an analytical model of resonant chains of any number of planets. Using this model, I show that a system captured in a resonant chain can librate around several possible equilibrium configurations. The probability of capture around each equilibrium depends on how the chain formed, and especially on the order in which the planets have been captured in the chain. Therefore, for an observed resonant chain, knowing around which equilibrium the chain is librating allows for constraints to be put on the formation and migration scenario of the system. I apply this reasoning to the four planets orbiting Kepler-223 in a 3:4:6:8 resonant chain. I show that the system is observed around one of the six equilibria predicted by the analytical model. Using N-body integrations, I show that the most favorable scenario to reproduce the observed configuration is to first capture the two intermediate planets, then the outermost, and finally the innermost.

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.

We describe an analytical method for computing the orbital parameters of a planet from the periodogram of a radial velocity signal. The method is very efficient and provides a good approximation of the orbital parameters. The accuracy is mainly limited by the accuracy of the computation of the Fourier decomposition of the signal which is sensitive to sampling and noise. Our method is complementary with more accurate (and more expensive in computer time) numerical algorithms (e.g. Levenberg-Marquardt, Markov chain Monte Carlo, genetic algorithms). Indeed, the analytical approximation can be used as an initial condition to accelerate the convergence of these numerical methods. Our method can be applied iteratively to search for multiple planets in the same system.

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.

Multi-planetary systems detected by the Kepler mission present an excess of planets close to first-order mean-motion resonances (2:1 and 3:2) but with a period ratio slightly higher than the resonant value. Several mechanisms have been proposed to explain this observation. Here we provide some clues that these near-resonant systems were initially in resonance and reached their current configuration through tidal dissipation. The argument that has been opposed to this scenario is that it only applies to the close-in systems and not to the farthest ones for which the tidal effect is too weak. Using the catalog of KOI of the Kepler mission, we show that the distributions of period ratio among the most close-in planetary systems and the farthest ones differ significantly. This distance dependent repartition is a strong argument in favor of the tidal dissipation scenario.

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.

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.

We numerically estimate the semi-major axis chaotic diffusion of the Vesta family asteroids induced by close encounters with 11 massive main belt asteroids : (1) Ceres, (2) Pallas, (3) Juno, (4) Vesta, (7) Iris, (10) Hygiea, (15) Eunomia, (19) Fortuna, (324) Bamberga, (532) Herculina, (704) Interamnia. We find that most of the diffusion is due to Ceres and Vesta. By extrapolating our results, we are able to constrain the global effect of close encounters with all the main belt asteroids. A comparison of this drift estimate with the one expected for the Yarkovsky effect shows that for asteroids whose diameter is larger than about 40 km, close encounters dominate the Yarkovsky effect. Overall, our findings confirm the standard scenario for the history of the Vesta family.