Small planets ($\sim$1--3.9 $\Rearth$) constitute more than half of the inventory of the 4000-plus exoplanets discovered so far. Smaller planets are sufficiently dense to be rocky, but those with radii larger than $\sim$1.6 $\Rearth$ are thought to display in many cases hydrogen/helium gaseous envelopes up to $\sim$30\% of the planetary mass. These low-mass planets are highly irradiated and the question of their origin, evolution, and possible links remains open. Here we show that close-in ocean planets affected by greenhouse effect display hydrospheres in supercritical state, which generate inflated atmospheres without invoking the presence of large hydrogen/helium gaseous envelopes. We present a new set of mass-radius relationships for ocean planets with different compositions and different equilibrium temperatures, which are found to be well adapted to low-density sub-Neptune planets. Our model suggests that super-Earths and water-rich sub-Neptunes could belong to the same family of planets, i.e. hydrogen/helium-free planets, with differences between their interiors simply resulting from the variation in the water content.

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.

The detection of low-mass transiting exoplanets in multiple systems brings new constraints to planetary formation and evolution processes and challenges the current planet formation theories. Nevertheless, only a mere fraction of the small planets detected by Kepler and K2 have precise mass measurements, which are mandatory to constrain their composition. We aim to characterise the planets that orbit the relatively bright star K2-138. This system is dynamically particular as it presents the longest chain known to date of planets close to the 3:2 resonance. We obtained 215 HARPS spectra from which we derived the radial-velocity variations of K2-138. Via a joint Bayesian analysis of both the K2 photometry and HARPS radial-velocities (RVs), we constrained the parameters of the six planets in orbit. The masses of the four inner planets, from b to e, are 3.1, 6.3, 7.9, and 13.0 $\mathrm{M}_{\oplus}$ with a precision of 34%, 20%, 18%, and 15%, respectively. The bulk densities are 4.9, 2.8, 3.2, and 1.8 g cm$^{-3}$, ranging from Earth to Neptune-like values. For planets f and g, we report upper limits. Finally, we predict transit timing variations of the order two to six minutes from the masses derived. Given its peculiar dynamics, K2-138 is an ideal target for transit timing variation (TTV) measurements from space with the upcoming CHaracterizing ExOPlanet Satellite (CHEOPS) to study this highly-packed system and compare TTV and RV masses.

Periodic radial velocity variations in the nearby M-dwarf star Gl411 are reported, based on measurements with the SOPHIE spectrograph. Current data do not allow us to distinguish between a 12.95-day period and its one-day alias at 1.08 days, but favour the former slightly. The velocity variation has an amplitude of 1.6 m/s, making this the lowest-amplitude signal detected with SOPHIE up to now. We have performed a detailed analysis of the significance of the signal and its origin, including extensive simulations with both uncorrelated and correlated noise, representing the signal induced by stellar activity. The signal is significantly detected, and the results from all tests point to its planetary origin. Additionally, the presence of an additional acceleration in the velocity time series is suggested by the current data. On the other hand, a previously reported signal with a period of 9.9 days, detected in HIRES velocities of this star, is not recovered in the SOPHIE data. An independent analysis of the HIRES dataset also fails to unveil the 9.9-day signal. If the 12.95-day period is the real one, the amplitude of the signal detected with SOPHIE implies the presence of a planet, called Gl411 b, with a minimum mass of around three Earth masses, orbiting its star at a distance of 0.079 AU. The planet receives about 3.5 times the insolation received by Earth, which implies an equilibrium temperature between 255 K and 350 K, and makes it too hot to be in the habitable zone. At a distance of only 2.5 pc, Gl411 b, is the third closest low-mass planet detected to date. Its proximity to Earth will permit probing its atmosphere with a combination of high-contrast imaging and high-dispersion spectroscopy in the next decade.

We present the detection of a Warm Neptune orbiting the M-dwarf Gl378, using radial velocity measurements obtained with the SOPHIE spectrograph at the Observatoire de Haute-Provence. The star was observed in the context of the SOPHIE exoplanets consortium's subprogramme dedicated to finding planets around M-dwarfs. Gl378 is an M1 star, of solar metallicity, at a distance of 14.96 pc. The single planet detected, Gl378 b, has a minimum mass of 13.02 $\rm M_{Earth}$ and an orbital period of 3.82 days, which place it at the lower boundary of the Hot Neptune desert. As one of only a few such planets around M-dwarfs, Gl378 b provides important clues to the evolutionary history of these close-in planets. In particular, the eccentricity of 0.1 may point to a high-eccentricity migration. The planet may also have lost part of its envelope due to irradiation.

We report the discovery of the super-Earth K2-265 b detected with K2 photometry. The planet orbits a bright (V_mag = 11.1) star of spectral type G8V with a period of 2.37 days. We obtained high-precision follow-up radial velocity measurements from HARPS, and the joint Bayesian analysis showed that K2-265 b has a radius of 1.71 +/- 0.11 R_earth and a mass of 6.54 +/- 0.84 M_earth, corresponding to a bulk density of 7.1 +/- 1.8 g/cm^3 . Composition analysis of the planet reveals an Earth-like, rocky interior, with a rock mass fraction of 80%. The short orbital period and small radius of the planet puts it below the lower limit of the photoevaporation gap, where the envelope of the planet could have eroded due to strong stellar irradiation, leaving behind an exposed core. Knowledge of the planet core composition allows us to infer the possible formation and evolution mechanism responsible for its current physical parameters.

We report the detection of two exoplanets and a further tentative candidate around the M-dwarf stars Gl96 and Gl617A, based on radial velocity measurements obtained with the SOPHIE spectrograph at the Observatoire de Haute-Provence. Both stars were observed in the context of the SOPHIE exoplanet consortium's dedicated M-dwarf subprogramme, which aims to detect exoplanets around nearby M-dwarf stars through a systematic survey. For Gl96, we present the discovery of a new exoplanet at 73.9 d with a minimum mass of 19.66 earth masses. Gl96 b has an eccentricity of 0.44, placing it among the most eccentric planets orbiting M stars. For Gl617A we independently confirm a recently reported exoplanet at 86.7 d with a minimum mass of 31.29 earth masses. Both Gl96 b and Gl617A b are potentially within the habitable zone, though Gl96 b's high eccentricity may take it too close to the star at periapsis.

The Earth, Venus, Mars, and some extrasolar terrestrial planets have a mass and radius that is consistent with a mass fraction of about 30% metallic core and 70% silicate mantle. At the inner frontier of the solar system, Mercury has a completely different composition, with a mass fraction of about 70% metallic core and 30% silicate mantle. Several formation or evolution scenarios are proposed to explain this metal-rich composition, such as a giant impact, mantle evaporation, or the depletion of silicate at the inner-edge of the proto-planetary disk. These scenarios are still strongly debated. Here we report the discovery of a multiple transiting planetary system (K2-229), in which the inner planet has a radius of 1.165+/-0.066 Rearth and a mass of 2.59+/-0.43 Mearth. This Earth-sized planet thus has a core-mass fraction that is compatible with that of Mercury, while it was expected to be similar to that of the Earth based on host-star chemistry. This larger Mercury analogue either formed with a very peculiar composition or it has evolved since, e.g. by losing part of its mantle. Further characterisation of Mercury-like exoplanets like K2-229 b will help putting the detailed in-situ observations of Mercury (with Messenger and BepiColombo) into the global context of the formation and evolution of solar and extrasolar terrestrial planets.

Modeling the interior of exoplanets is essential to go further than the conclusions provided by mean density measurements. In addition to the still limited precision on the planets' fundamental parameters, models are limited by the existence of degeneracies on their compositions. Here we present a model of internal structure dedicated to the study of solid planets up to ~10 Earth masses, i.e. Super-Earths. When the measurement is available, the assumption that the bulk Fe/Si ratio of a planet is similar to that of its host star allows us to significantly reduce the existing degeneracy and more precisely constrain the planet's composition. Based on our model, we provide an update of the mass-radius relationships used to provide a first estimate of a planet's composition from density measurements. Our model is also applied to the cases of two well-known exoplanets, CoRoT-7b and Kepler-10b, using their recently updated parameters. The core mass fractions of CoRoT-7b and Kepler-10b are found to lie within the 10-37% and 10-33% ranges, respectively, allowing both planets to be compatible with an Earth-like composition. We also extend the recent study of Proxima Centauri b, and show that its radius may reach 1.94 Earth radii in the case of a 5 Earth masses planet, as there is a 96.7% probability that the real mass of Proxima Centauri b is below this value.

The multi-planetary system HD 106315 was recently found in K2 data . The planets have periods of $P_b \sim9.55$ and $P_c \sim 21.06\,$days, and radii of $ r_b = 2.44 \pm 0.17\, $ and $r_c = 4.35 \pm 0.23\, $ $R_{\oplus}$. The brightness of the host star (V=9.0 mag) makes it an excellent target for transmission spectroscopy. However, to interpret transmission spectra it is crucial to measure the planetary masses. We obtained high precision radial velocities for HD~106315 to determine the mass of the two transiting planets discovered with Kepler K2. Our successful observation strategy was carefully tailored to mitigate the effect of stellar variability. We modelled the new radial velocity data together with the K2 transit photometry and a new ground-based partial transit of HD 106315c to derive system parameters. We estimate the mass of HD 106315b to be 12.6 $\pm$ 3.2 $M_{\oplus}$ and the density to be $4.7 \pm 1.7\, g\,cm^{-3}$, while for HD 106315c we estimate a mass of 15.2 $\pm$ 3.7 $M_{\oplus}$ and a density of $1.01 \pm 0.29\, $g\,cm$^{-3}$. Hence, despite planet c having a radius almost twice as large as planet b, their masses are consistent with one another. We conclude that HD 106315c has a thick hydrogen-helium gaseous envelope. A detailed investigation of HD 106315b using a planetary interior model constrains the core mass fraction to be 5-29\%, and the water mass fraction to be 10-50\%. An alternative, not considered by our model, is that HD 106315b is composed of a large rocky core with a thick H-He envelope. Transmission spectroscopy of these planets will give insight into their atmospheric compositions and also help constrain their core compositions.

The ice giants Uranus and Neptune are the least understood class of planets in our solar system but the most frequently observed type of exoplanets. Presumed to have a small rocky core, a deep interior comprising ~70% heavy elements surrounded by a more dilute outer envelope of H2 and He, Uranus and Neptune are fundamentally different from the better-explored gas giants Jupiter and Saturn. Because of the lack of dedicated exploration missions, our knowledge of the composition and atmospheric processes of these distant worlds is primarily derived from remote sensing from Earth-based observatories and space telescopes. As a result, Uranus's and Neptune's physical and atmospheric properties remain poorly constrained and their roles in the evolution of the Solar System not well understood. Exploration of an ice giant system is therefore a high-priority science objective as these systems (including the magnetosphere, satellites, rings, atmosphere, and interior) challenge our understanding of planetary formation and evolution. Here we describe the main scientific goals to be addressed by a future in situ exploration of an ice giant. An atmospheric entry probe targeting the 10-bar level, about 5 scale heights beneath the tropopause, would yield insight into two broad themes: i) the formation history of the ice giants and, in a broader extent, that of the Solar System, and ii) the processes at play in planetary atmospheres. The probe would descend under parachute to measure composition, structure, and dynamics, with data returned to Earth using a Carrier Relay Spacecraft as a relay station. In addition, possible mission concepts and partnerships are presented, and a strawman ice-giant probe payload is described. An ice-giant atmospheric probe could represent a significant ESA contribution to a future NASA ice-giant flagship mission.

We explore the possible Proxima Centauri b's interiors assuming the planet belongs to the class of dense solid planets (rocky with possible addition of water) and derive the corresponding radii. To do so, we use an internal structure model that computes the radius of the planet along with the locations of the different layers of materials, assuming that its mass and bulk composition are known. Lacking detailed elementary abundances of the host star to constrain the planet's composition, we base our model on solar system values. We restrained the simulations to the case of solid planets without massive atmospheres. With these assumptions, the possible radius of Proxima Centauri b spans the 0.94--1.40 $R_\oplus$ range. The minimum value is obtained considering a 1.10 $M_\oplus$ Mercury-like planet with a 65% core mass fraction, whereas the highest radius is reached for 1.46 $M_\oplus$ with 50% water in mass, constituting an ocean planet. Although this range of radii still allows very different planet compositions, it helps characterizing many aspects of Proxima Centauri b, such as the formation conditions of the system or the current amount of water on the planet. This work can also help ruling out future measurements of the planet's radius that would be physically incompatible with a solid planetary body.

The origin of Phobos and Deimos is still an open question. Currently, none of the three proposed scenarios for their origin (intact capture of two distinct outer solar system small bodies, co-accretion with Mars, and accretion within an impact-generated disk) is able to reconcile their orbital and physical properties. Here, we investigate the expected mineralogical composition and size of the grains from which the moons once accreted assuming they formed within an impact-generated accretion disk. A comparison of our results with the present day spectral properties of the moons allows us to conclude that their building blocks cannot originate from a magma phase, thus preventing their formation in the innermost part of the disk. Instead, gas-to-solid condensation of the building blocks in the outer part of an extended gaseous disk is found as a possible formation mechanism as it does allow reproducing both the spectral and physical properties of the moons. Such a scenario may finally reconcile their orbital and physical properties alleviating the need to invoke an unlikely capture scenario to explain their physical properties.

Molecular oxygen has been detected in the coma of comet 67P/Churyumov-Gerasimenko with abundances in the 1-10% range by the ROSINA-DFMS instrument on board the Rosetta spacecraft. Here we find that the radiolysis of icy grains in low-density environments such as the presolar cloud may induce the production of large amounts of molecular oxygen. We also show that molecular oxygen can be efficiently trapped in clathrates formed in the protosolar nebula, and that its incorporation as crystalline ice is highly implausible because this would imply much larger abundances of Ar and N2 than those observed in the coma. Assuming that radiolysis has been the only O2 production mechanism at work, we conclude that the formation of comet 67P/Churyumov-Gerasimenko is possible in a dense and early protosolar nebula in the framework of two extreme scenarios: (1) agglomeration from pristine amorphous icy grains/particles formed in ISM and (2) agglomeration from clathrates that formed during the disk's cooling. The former scenario is found consistent with the strong correlation between O2 and H2O observed in 67P/C-G's coma while the latter scenario requires that clathrates formed from ISM icy grains that crystallized when entering the protosolar nebula.

On November 12, 2014, the ESA/Rosetta descent module Philae landed on the Abydos site of comet 67P/Churyumov-Gerasimenko. Aboard this module, the Ptolemy mass spectrometer measured a CO/CO2 ratio of 0.07 +/- 0.04 which differs substantially from the value obtained in the coma by the Rosetta/ROSINA instrument, suggesting a heterogeneity in the comet nucleus. To understand this difference, we investigated the physico-chemical properties of the Abydos subsurface leading to CO/CO2 ratios close to that observed by Ptolemy at the surface of this region. We used a comet nucleus model that takes into account different water ice phase changes (amorphous ice, crystalline ice and clathrates), as well as diffusion of molecules throughout the pores of the matrix. The input parameters of the model were optimized for the Abydos site and the ROSINA CO/CO2 measured ratio is assumed to correspond to the bulk value in the nucleus. We find that all considered structures of water ice are able to reproduce the Ptolemy observation with a time difference not exceeding ~50 days, i.e. lower than ~2% on 67P/Churyumov-Gerasimenko's orbital period. The suspected heterogeneity of 67P/Churyumov-Gerasimenko's nucleus is also found possible only if it is constituted of crystalline ices. If the icy phase is made of amorphous ice or clathrates, the difference between Ptolemy and ROSINA's measurements would rather originate from the spatial variations in illumination on the nucleus surface. An eventual new measurement of the CO/CO2 ratio at Abydos by Ptolemy could be decisive to distinguish between the three water ice structures.

We investigate the thermal evolution of comet 67P/Churyumov-Gerasimenko's subsurface in the Seth_01 region, where active pits have been observed by the ESA/Rosetta mission. Our simulations show that clathrate destabilization and amorphous ice crystallization can occur at depths corresponding to those of the observed pits in a timescale shorter than 67P/Churyumov-Gerasimenko's lifetime in the comet's activity zone in the inner solar system. Sublimation of crystalline ice down to such depths is possible only in the absence of a dust mantle, which requires the presence of dust grains in the matrix small enough to be dragged out by gas from the pores. Our results are consistent with both pits formation via sinkholes or subsequent to outbursts, the dominant process depending on the status of the subsurface porosity. A sealed dust mantle would favor episodic and disruptive outgassing as a result of an increasing gas pressure in the pores, while a high porosity should allow the formation of large voids in the subsurface due to the continuous escape of volatiles. We finally conclude that the subsurface of 67P/Churyumov-Gerasimenko is not uniform at a spatial scale of 100-200~m.