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.

One of the most powerful methods used to estimate sky-projected spin-orbit angles of exoplanetary systems is through a spectroscopic transit observation known as the Rossiter-McLaughlin (RM) effect. So far mostly single RM observations have been used to estimate the spin-orbit angle, and thus there have been no studies regarding the variation of estimated spin-orbit angle from transit to transit. Stellar activity can alter the shape of photometric transit light curves and in a similar way they can deform the RM signal. In this paper we discuss several RM observations, obtained using the HARPS spectrograph, of known transiting planets that all transit extremely active stars, and by analyzing them individually we assess the variation in the estimated spin-orbit angle. Our results reveal that the estimated spin-orbit angle can vary significantly (up to 42 degrees) from transit to transit, due to variation in the configuration of stellar active regions over different nights. This finding is almost two times larger than the expected variation predicted from simulations. We could not identify any meaningful correlation between the variation of estimated spin-orbit angles and the stellar magnetic activity indicators. We also investigated two possible approaches to mitigate the stellar activity influence on RM observations. The first strategy was based on obtaining several RM observations and folding them to reduce the stellar activity noise. Our results demonstrated that this is a feasible and robust way to overcome this issue. The second approach is based on acquiring simultaneous high-precision short-cadence photometric transit light curves using TRAPPIST/SPECULOOS telescopes, which provide more information about the stellar active region's properties and allow a better RM modeling.

Radial-velocity (RV) signals induce RV variations an order of magnitude larger than the signal created by the orbit of Earth-twins, thus preventing their detection. The goal of this paper is to compare the efficiency of the different methods used to deal with stellar signals to recover extremely low-mass planets despite. However, because observed RV variations at the m/s precision level or below is a combination of signals induced by unresolved orbiting planets, by the star, and by the instrument, performing such a comparison using real data is extremely challenging. To circumvent this problem, we generated simulated RV measurements including realistic stellar and planetary signals. Different teams analyzed blindly those simulated RV measurements, using their own method to recover planetary signals despite stellar RV signals. By comparing the results obtained by the different teams with the planetary and stellar parameters used to generate the simulated RVs, it is therefore possible to compare the efficiency of these different methods. The most efficient methods to recover planetary signals {take into account the different activity indicators,} use red-noise models to account for stellar RV signals and a Bayesian framework to provide model comparison in a robust statistical approach. Using the most efficient methodology, planets can be found down to K/N= K_pl/RV_rms*sqrt{N_obs}=5 with a threshold of K/N=7.5 at the level of 80-90% recovery rate found for a number of methods. These recovery rates drop dramatically for K/N smaller than this threshold. In addition, for the best teams, no false positives with K/N > 7.5 were detected, while a non-negligible fraction of them appear for smaller K/N. A limit of K/N = 7.5 seems therefore a safe threshold to attest the veracity of planetary signals for RV measurements with similar properties to those of the different RV fitting challenge systems.

Planetary rotation rate is a key parameter in determining atmospheric circulation and hence the spatial pattern of clouds. Since clouds can exert a dominant control on planetary radiation balance, rotation rate could be critical for determining mean planetary climate. Here we investigate this idea using a three-dimensional general circulation model with a sophisticated cloud scheme. We find that slowly rotating planets (like Venus) can maintain an Earth-like climate at nearly twice the stellar flux as rapidly rotating planets (like Earth). This suggests that many exoplanets previously believed to be too hot may actually be habitable, depending on their rotation rate. The explanation for this behavior is that slowly rotating planets have a weak Coriolis force and long daytime illumination, which promotes strong convergence and convection in the substellar region. This produces a large area of optically thick clouds, which greatly increases the planetary albedo. In contrast, on rapidly rotating planets a much narrower belt of clouds form in the deep tropics, leading to a relatively low albedo. A particularly striking example of the importance of rotation rate suggested by our simulations is that a planet with modern Earth's atmosphere, in Venus' orbit, and with modern Venus' (slow) rotation rate would be habitable. This would imply that if Venus went through a runaway greenhouse, it had a higher rotation rate at that time.

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.

The hot-Jupiter WASP-10b was reported by Maciejewski et al. (2011a,b) to show transit timing variations (TTV) with an amplitude of ~ 3.5 minutes. These authors proposed that the observed TTVs were caused by a 0.1 MJup perturbing companion with an orbital period of ~ 5.23 d, and hence, close to the outer 5:3 mean motion resonance with WASP-10b. To test this scenario, we present eight new transit light curves of WASP-10b obtained with the Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 previously published ones, were modelled with a Markov-Chain Monte-Carlo transit fitting code. (...) Our homogeneously derived transit times do not support the previous claimed TTV signal, which was strongly dependent on 2 previously published transits that have been incorrectly normalised. Nevertheless, a linear ephemeris is not a statistically good fit to the transit times of WASP-10b. We show that the observed transit time variations are due to spot occultation features or systematics. We discuss and exemplify the effects of occultation spot features in the measured transit times and show that despite spot occultation during egress and ingress being difficult to distinguish in the transit light curves, they have a significant effect in the measured transit times. We conclude that if we account for spot features, the transit times of WASP-10 are consistent with a linear ephemeris with the exception of one transit (epoch 143) which is a partial transit. Therefore, there is currently no evidence for the existence of a companion to WASP-10b. Our results support the lack of TTVs of hot-Jupiters reported for the Kepler sample.

We consider the solid-solid interactions in the two body problem. The relative equilibria have been previously studied analytically and general motions were numerically analyzed using some expansion of the gravitational potential up to the second order, but only when there are no direct interactions between the orientation of the bodies. Here we expand the potential up to the fourth order and we show that the secular problem obtained after averaging over fast angles, as for the precession model of Boue and Laskar [Boue, G., Laskar, J., 2006. Icarus 185, 312-330], is integrable, but not trivially. We describe the general features of the motions and we provide explicit analytical approximations for the solutions. We demonstrate that the general solution of the secular system can be decomposed as a uniform precession around the total angular momentum and a periodic symmetric orbit in the precessing frame. More generally, we show that for a general n-body system of rigid bodies in gravitational interaction, the regular quasiperiodic solutions can be decomposed into a uniform precession around the total angular momentum, and a quasiperiodic motion with one frequency less in the precessing frame.