Although several thousands of exoplanets have now been detected and characterized, observational biases have led to a paucity of long-period, low-mass exoplanets with measured masses and a corresponding lag in our understanding of such planets. In this paper we report the mass estimation and characterization of the long-period exoplanet Kepler-538b. This planet orbits a Sun-like star (V = 11.27) with M_* = 0.892 +/- (0.051, 0.035) M_sun and R_* = 0.8717 +/- (0.0064, 0.0061) R_sun. Kepler-538b is a 2.215 +/- (0.040, 0.034) R_earth sub-Neptune with a period of P = 81.73778 +/- 0.00013 d. It is the only known planet in the system. We collected radial velocity (RV) observations with HIRES on Keck I and HARPS-N on the TNG. We characterized stellar activity by a Gaussian process with a quasi-periodic kernel applied to our RV and cross correlation function full width at half maximum (FWHM) observations. By simultaneously modeling Kepler photometry, RV, and FWHM observations, we found a semi-amplitude of K = 1.68 +/- (0.39, 0.38) m s^-1 and a planet mass of M_p = 10.6 +/- (2.5, 2.4) M_earth. Kepler-538b is the smallest planet beyond P = 50 d with an RV mass measurement. The planet likely consists of a significant fraction of ices (dominated by water ice), in addition to rocks/metals, and a small amount of gas. Sophisticated modeling techniques such as those used in this paper, combined with future spectrographs with ultra high-precision and stability will be vital for yielding more mass measurements in this poorly understood exoplanet regime. This in turn will improve our understanding of the relationship between planet composition and insolation flux and how the rocky to gaseous transition depends on planetary equilibrium temperature.

We present confirmation of the planetary nature of PH-2b, as well as the first mass estimates for the two planets in the Kepler-103 system. PH-2b and Kepler-103c are both long-period and transiting, a sparsely-populated category of exoplanet. We use {\it Kepler} light-curve data to estimate a radius, and then use HARPS-N radial velocities to determine the semi-amplitude of the stellar reflex motion and, hence, the planet mass. For PH-2b we recover a 3.5-$\sigma$ mass estimate of $M_p = 109^{+30}_{-32}$ M$_\oplus$ and a radius of $R_p = 9.49\pm0.16$ R$_\oplus$. This means that PH-2b has a Saturn-like bulk density and is the only planet of this type with an orbital period $P > 200$ days that orbits a single star. We find that Kepler-103b has a mass of $M_{\text{p,b}} = 11.7^{+4.31}_{-4.72}$ M$_{\oplus}$ and Kepler-103c has a mass of $M_{\text{p,c}} = 58.5^{+11.2}_{-11.4}$ M$_{\oplus}$. These are 2.5$\sigma$ and 5$\sigma$ results, respectively. With radii of $R_{\text{p,b}} = 3.49^{+0.06}_{-0.05}$ R$_\oplus$, and $R_{\text{p,c}} = 5.45^{+0.18}_{-0.17}$ R$_\oplus$, these results suggest that Kepler-103b has a Neptune-like density, while Kepler-103c is one of the highest density planets with a period $P > 100$ days. By providing high-precision estimates for the masses of the long-period, intermediate-mass planets PH-2b and Kepler-103c, we increase the sample of long-period planets with known masses and radii, which will improve our understanding of the mass-radius relation across the full range of exoplanet masses and radii.

TESS is revolutionising the search for planets orbiting bright and nearby stars. In sectors 3 and 4, TESS observed TOI-402 (TIC-120896927), a bright V=9.1 K1 dwarf also known as HD 15337, and found two transiting signals with period of 4.76 and 17.18 days and radius of 1.90 and 2.21\,\Rearth. This star was observed as part of the radial-velocity search for planets using the HARPS spectrometer, and 85 precise radial-velocity measurements were obtained over a period of 14 years. In this paper, we analyse the HARPS radial-velocity measurements in hand to confirm the planetary nature of these two signals. By reanalysing TESS photometry and host star parameters using EXOFASTv2, we find that TOI-402.01 and TOI-402.02 have periods of 4.75642$\pm$0.00021 and 17.1784$\pm$0.0016 days and radii of 1.70$\pm$0.06 and 2.52$\pm$0.11\,\Rearth\,(precision 3.6 and 4.2\%), respectively. By analysing the HARPS radial-velocity measurements, we find that those planets are both super-Earths with masses of 7.20$\pm$0.81 and 8.79$\pm$1.67\,\Mearth\,(precision 11.3 and 19.0\%), and small eccentricities compatible with zero at 2$\sigma$. Although having rather similar masses, the radius of these two planets is really different, putting them on different sides of the radius gap. With stellar irradiation 160 times more important than Earth for TOI-402.01 and only 29 times more for TOI-402.02, it is likely that photo-evaporation is at the origin of this radius difference. Those two planets, being in the same system and therefore being in the same irradiation environment are therefore extremely important to perform comparative exoplanetology across the evaporation valley and thus bring constraints on the mechanisms responsible for the radius gap.

The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets' masses, obstructing detection of Earth twins. Since 2015 July we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 years of observation, interpreting them in the context of spatially-resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frame, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise limited radial-velocity precision better than 0.43 m s$^{-1}$ per 5-minute observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s$^{-1}$. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s$^{-1}$ and an autocorrelation half-life of 15 minutes. On longer timescales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s$^{-1}$ and is correlated with the line asymmetry, but leads it by 3 days. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity.

The study of Earth-mass extrasolar planets via the radial-velocity technique and the measurement of the potential cosmological variability of fundamental constants call for very-high-precision spectroscopy at the level of $\updelta\lambda/\lambda<10^{-9}$. Wavelength accuracy is obtained by providing two fundamental ingredients: 1) an absolute and information-rich wavelength source and 2) the ability of the spectrograph and its data reduction of transferring the reference scale (wavelengths) to a measurement scale (detector pixels) in a repeatable manner. The goal of this work is to improve the wavelength calibration accuracy of the HARPS spectrograph by combining the absolute spectral reference provided by the emission lines of a thorium-argon hollow-cathode lamp (HCL) with the spectrally rich and precise spectral information of a Fabry-P\'erot-based calibration source. On the basis of calibration frames acquired each night since the Fabry-P\'erot etalon was installed on HARPS in 2011, we construct a combined wavelength solution which fits simultaneously the thorium emission lines and the Fabry-P\'erot lines. The combined fit is anchored to the absolute thorium wavelengths, which provide the `zero-point' of the spectrograph, while the Fabry-P\'erot lines are used to improve the (spectrally) local precision. The obtained wavelength solution is verified for auto-consistency and tested against a solution obtained using the HARPS Laser-Frequency Comb (LFC). The combined thorium+Fabry-P\'erot wavelength solution shows significantly better performances compared to the thorium-only calibration. The presented techniques will therefore be used in the new HARPS and HARPS-N pipeline, and will be exported to the ESPRESSO spectrograph.

K2-291 (EPIC 247418783) is a solar-type star with a radius of R_star = 0.899 $\pm$ 0.034 R_sun and mass of M_star=0.934 $\pm$ 0.038 M_sun. From K2 C13 data, we found one super-Earth planet (R_p = 1.589+0.095-0.072 R_Earth) transiting this star on a short period orbit (P = 2.225177 +6.6e-5 -6.8e-5 days). We followed this system up with adaptive-optic imaging and spectroscopy to derive stellar parameters, search for stellar companions, and determine a planet mass. From our 75 radial velocity measurements using HIRES on Keck I and HARPS-N on Telescopio Nazionale Galileo, we constrained the mass of EPIC 247418783b to M_p = 6.49 $\pm$ 1.16 M_Earth. We found it necessary to model correlated stellar activity radial velocity signals with a Gaussian process in order to more accurately model the effect of stellar noise on our data; the addition of the Gaussian process also improved the precision of this mass measurement. With a bulk density of 8.84+2.50-2.03 g cm-3, the planet is consistent with an Earth-like rock/iron composition and no substantial gaseous envelope. Such an envelope, if it existed in the past, was likely eroded away by photo-evaporation during the first billion years of the star's lifetime.

K2-36 is a K dwarf orbited by two small ($R_{\rm b}=1.43\pm0.08$ $R_\oplus$ and $R_{\rm c}=3.2\pm0.3$ $R_\oplus$), close-in ($a_{\rm b}$=0.022 AU and $a_{\rm c}$=0.054 AU) transiting planets discovered by Kepler/K2. They are representatives of two families of small planets ($R_{\rm p}$<4 $R_\oplus$) recently emerged from the analysis of Kepler data, with likely a different structure, composition and evolutionary pathways. We revise the fundamental stellar parameters and the sizes of the planets, and provide the first measurement of their masses and bulk densities, which we use to infer their structure and composition. We observed K2-36 with the HARPS-N spectrograph over $\sim$3.5 years, collecting 81 useful radial velocity measurements. The star is active, with evidence for increasing levels of magnetic activity during the observing time span. The radial velocity scatter is $\sim$17 \ms due to the stellar activity contribution, which is much larger that the semi-amplitudes of the planetary signals. We tested different methods for mitigating the stellar activity contribution to the radial velocity time variations and measuring the planet masses with good precision. We found that K2-36 is likely a $\sim$1 Gyr old system, and by treating the stellar activity through a Gaussian process regression, we measured the planet masses $m_{\rm b}$=3.9$\pm$1.1 $M_\oplus$ and $m_{\rm c}$=7.8$\pm$2.3 $M_\oplus$. The derived planet bulk densities $\rho_{\rm b}$=7.2$^{+2.5}_{-2.1}$ $g/cm^{3}$ and $\rho_{\rm c}$=1.3$^{+0.7}_{-0.5}$ $g/cm^{3}$ point out that K2-36\,b has a rocky, Earth-like composition, and K2-36\,c is a low-density sub-Neptune. Composed of two planets with similar orbital separations but different densities, K2-36 represents an optimal laboratory for testing the role of the atmospheric escape in driving the evolution of close-in, low-mass planets after $\sim$1 Gyr from their formation.

We present radial-velocity measurement of eight stars observed with the HARPS Echelle spectrograph mounted on the 3.6-m telescope in La Silla (ESO, Chile). Data span more than ten years and highlight the long-term stability of the instrument. We search for potential planets orbiting HD20003, HD20781, HD21693, HD31527, HD45184, HD51608, HD134060 and HD136352 to increase the number of known planetary systems and thus better constrain exoplanet statistics. After a preliminary phase looking for signals using generalized Lomb-Scargle periodograms, we perform a careful analysis of all signals to separate \emph{bona-fide} planets from spurious signals induced by stellar activity and instrumental systematics. We finally secure the detection of all planets using the efficient MCMC available on the Data and Analysis Center for Exoplanets (DACE web-platform), using model comparison whenever necessary. In total, we report the detection of twenty new super-Earth to Neptune-mass planets, with minimum masses ranging from 2 to 30 M$_{\rm Earth}$, and periods ranging from 3 to 1300 days. By including CORALIE and HARPS measurements of HD20782 to the already published data, we also improve the characterization of the extremely eccentric Jupiter orbiting this host.

We investigated the discrepancy between planetary mass determination using the transit timing variations (TTVs) and radial velocities (RVs), by analysing the multi-planet system Kepler-9. Despite being the first system characterised with TTVs, there are several discrepant solutions in the literature, with those reporting lower planetary densities being apparently in disagreement with high-precision RV observations. To resolve this, we gathered HARPS-N RVs at epochs that maximised the difference between the predicted RV curves from discrepant solutions in the literature. We also re-analysed the full Kepler data-set and performed a dynamical fit, within a Bayesian framework, using the newly derived central and duration times of the transits. We compared these results with the RV data and found that our solution better describes the RV observations, despite the masses of the planets being nearly half that presented in the discovery paper. We therefore confirm that the TTV method can provide mass determinations that agree with those determined using high-precision RVs. The low densities of the planets place them in the scarcely populated region of the super-Neptunes / inflated sub-Saturns in the mass-radius diagram.

Super-Earths belong to a class of planet not found in the Solar System, but which appear common in the Galaxy. Given that some super-Earths are rocky, while others retain substantial atmospheres, their study can provide clues as to the formation of both rocky planets and gaseous planets, and - in particular - they can help to constrain the role of photo-evaporation in sculpting the exoplanet population. GJ 9827 is a system already known to host 3 super-Earths with orbital periods of 1.2, 3.6 and 6.2 days. Here we use new HARPS-N radial velocity measurements, together with previously published radial velocities, to better constrain the properties of the GJ 9827 planets. Our analysis can't place a strong constraint on the mass of GJ 9827 c, but does indicate that GJ 9827 b is rocky with a composition that is probably similar to that of the Earth, while GJ 9827 d almost certainly retains a volatile envelope. Therefore, GJ 9827 hosts planets on either side of the radius gap that appears to divide super-Earths into pre-dominantly rocky ones that have radii below $\sim 1.5 R_\oplus$, and ones that still retain a substantial atmosphere and/or volatile components, and have radii above $\sim 2 R_\oplus$. That the less heavily irradiated of the 3 planets still retains an atmosphere, may indicate that photoevaporation has played a key role in the evolution of the planets in this system.

(abridged) Even if the HARPS spectrograph has been operational for more than 15 years and it provides among the most precise Doppler measurements, improvements are still possible. One known problem, for instance, is the non-fully regular block-stitching of the CCDs, which introduces, in some cases, one-year period parasitic signals in the measured radial velocity. The aim is to improve the wavelength calibration of HARPS to push further its planet-detection capabilities. The properties of the CCD stitching-induced pixel-size anomalies are determined with LED flat-field frames, and then a physical, gap-corrected map of the CCDs is used for the fitting model of the spectral orders. We also use a new thorium line list, based on much higher-accuracy measurements than the one used up to now. We derive new wavelength solutions for the 15 years of HARPS data, both before and after the 2015 fibre upgrade. We demonstrate that we correct the gap anomalies by computing the wavelength solutions of laser frequency comb exposures, both with and without taking the gap correction into account. By comparing the rms of the most stable stars of the HARPS sample, we show that we globally decrease the radial velocity dispersion of the data, especially for the data acquired after the change of fibres. Finally, the comparative analysis of several individual systems shows that we manage to attenuate the periodogram power at one year in most cases. The analysis of the RVs derived from individual stellar lines also shows that we correct the stitching-induced RV variation. This improved calibration of the HARPS spectrograph allows to go deeper in the search for low-amplitude radial-velocity signals. It will be further improved by combining the thorium calibration spectra with laser frequency comb and Fabry-Perot calibration spectra, and not only for HARPS but notably also for HARPS-N and ESPRESSO.

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.

Ultra-short period (USP) planets are a class of low mass planets with periods shorter than one day. Their origin is still unknown, with photo-evaporation of mini-Neptunes and in-situ formation being the most credited hypotheses. Formation scenarios differ radically in the predicted composition of USP planets, it is therefore extremely important to increase the still limited sample of USP planets with precise and accurate mass and density measurements. We report here the characterization of an USP planet with a period of 0.28 days around K2-141 (EPIC 246393474), and the validation of an outer planet with a period of 7.7 days in a grazing transit configuration. We derived the radii of the planets from the K2 light curve and used high-precision radial velocities gathered with the HARPS-N spectrograph for mass measurements. For K2-141b we thus inferred a radius of $1.51\pm0.05~R_\oplus$ and a mass of $5.08\pm0.41~M_\oplus$, consistent with a rocky composition and lack of a thick atmosphere. K2-141c is likely a Neptune-like planet, although due to the grazing transits and the non-detection in the RV dataset, we were not able to put a strong constraint on its density. We also report the detection of secondary eclipses and phase curve variations for K2-141b. The phase variation can be modeled either by a planet with a geometric albedo of $0.30 \pm 0.06$ in the Kepler bandpass, or by thermal emission from the surface of the planet at $\sim$3000K. Only follow-up observations at longer wavelengths will allow us to distinguish between these two scenarios.

We present the analysis of the entire HARPS observations of three stars that host planetary systems: HD1461, HD40307, and HD204313. The data set spans eight years and contains more than 200 nightly averaged velocity measurements for each star. This means that it is sensitive to both long-period and low-mass planets and also to the effects induced by stellar activity cycles. We modelled the data using Keplerian functions that correspond to planetary candidates and included the short- and long-term effects of magnetic activity. A Bayesian approach was taken both for the data modelling, which allowed us to include information from activity proxies such as $\log{(R'_{\rm HK})}$ in the velocity modelling, and for the model selection, which permitted determining the number of significant signals in the system. The Bayesian model comparison overcomes the limitations inherent to the traditional periodogram analysis. We report an additional super-Earth planet in the HD1461 system. Four out of the six planets previously reported for HD40307 are confirmed and characterised. We discuss the remaining two proposed signals. In particular, we show that when the systematic uncertainty associated with the techniques for estimating model probabilities are taken into account, the current data are not conclusive concerning the existence of the habitable-zone candidate HD40307 g. We also fully characterise the Neptune-mass planet that orbits HD204313 in 34.9 days.