The RoPES project with HARPS and HARPS-N. I. A system of super-Earths orbiting the moderately active K-dwarf HD 176986

We report the discovery of a system of two super-Earths orbiting the moderately active K-dwarf HD 176986. This work is part of the RoPES RV program of G- and K-type stars, which combines radial velocities (RVs) from the HARPS and HARPS-N spectrographs to search for short-period terrestrial planets. HD 176986 b and c are super-Earth planets with masses of 5.74 and 9.18 M$_{\oplus}$, orbital periods of 6.49 and 16.82 days, and distances of 0.063 and 0.119 AU in orbits that are consistent with circular. The host star is a K2.5 dwarf, and despite its modest level of chromospheric activity (log(R'hk) = - 4.90 +- 0.04), it shows a complex activity pattern. Along with the discovery of the planets, we study the magnetic cycle and rotation of the star. HD 176986 proves to be suitable for testing the available RV analysis technique and further our understanding of stellar activity.Comment: 21 pages, 24 figures, 7 table

The HADES RV Programme with HARPS-N at TNG XI. GJ 685 b: a warm super-Earth around an active M dwarf

Small rocky planets seem to be very abundant around low-mass M-type stars. Their actual planetary population is however not yet precisely understood. Currently several surveys aim to expand the statistics with intensive detection campaigns, both photometric and spectroscopic. We analyse 106 spectroscopic HARPS-N observations of the active M0-type star GJ 685 taken over the past five years. We combine these data with photometric measurements from different observatories to accurately model the stellar rotation and disentangle its signals from genuine Doppler planetary signals in the RV data. We run an MCMC analysis on the RV and activity indexes time series to model the planetary and stellar signals present in the data, applying Gaussian Process regression technique to deal with the stellar activity signals. We identify three periodic signals in the RV time series, with periods of 9, 24, and 18 d. Combining the analyses of the photometry of the star with the activity indexes derived from the HARPS-N spectra, we identify the 18 d and 9 d signals as activity-related, corresponding to the stellar rotation period and its first harmonic respectively. The 24 d signals shows no relations with any activity proxy, so we identify it as a genuine planetary signal. We find the best-fit model describing the Doppler signal of the newly-found planet, GJ 685\,b, corresponding to an orbital period $P_b = 24.160^{+0.061}_{-0.047}$ d and a minimum mass $M_P \sin i = 9.0^{+1.7}_{-1.8}$ M$_\oplus$. We also study a sample of 70 RV-detected M-dwarf planets, and present new statistical evidence of a difference in mass distribution between the populations of single- and multi-planet systems, which can shed new light on the formation mechanisms of low-mass planets around late-type stars.Comment: 18 pages, 13 figures, accepted for publication in A&

HADES RV Programme with HARPS-N at TNG. VII. Rotation and activity of M-Dwarfs from time-series high-resolution spectroscopy of chromospheric indicators

We aim to investigate the presence of signatures of magnetic cycles and rotation on a sample of 71 early M-dwarfs from the HADES RV programme using high-resolution time-series spectroscopy of the Ca II H & K and Halpha chromospheric activity indicators, the radial velocity series, the parameters of the cross correlation function and the V-band photometry. We used mainly HARPS-N spectra, acquired over four years, and add HARPS spectra from the public ESO database and ASAS photometry light-curves as support data, extending the baseline of the observations of some stars up to 12 years. We provide log(R'hk) measurements for all the stars in the sample, cycle length measurements for 13 stars, rotation periods for 33 stars and we are able to measure the semi-amplitude of the radial velocity signal induced by rotation in 16 stars. We complement our work with previous results and confirm and refine the previously reported relationships between the mean level of chromospheric emission, measured by the log(R'hk), with the rotation period, and with the measured semi-amplitude of the activity induced radial velocity signal for early M-dwarfs. We searched for a possible relation between the measured rotation periods and the lengths of the magnetic cycle, finding a weak correlation between both quantities. Using previous v sin i measurements we estimated the inclinations of the star's poles to the line of sight for all the stars in the sample, and estimate the range of masses of the planets GJ 3998 b and c (2.5 - 4.9 Mearth and 6.3 - 12.5 Mearth), GJ 625 b (2.82 Mearth), GJ 3942 b (7.1 - 10.0 Mearth) and GJ 15A b (3.1 - 3.3 Mearth), assuming their orbits are coplanar with the stellar rotation.Comment: 19 pages, 16 figures, 10 table

HADES RV Programme with HARPS-N at TNG VI. GJ 3942 b behind dominant activity signals

Short- to mid-term magnetic phenomena on the stellar surface of M-type stars cannot only resemble the effects of planets in radial velocity data, but also may hide them. We analyze 145 spectroscopic HARPS-N observations of GJ 3942 taken over the past five years and additional photometry to disentangle stellar activity effects from genuine Doppler signals as a result of the orbital motion of the star around the common barycenter with its planet. To achieve this, we use the common methods of pre-whitening, and treat the correlated red noise by a first-order moving average term and by Gaussian-process regression following an MCMC analysis. We identify the rotational period of the star at 16.3 days and discover a new super-Earth, GJ 3942 b, with an orbital period of 6.9 days and a minimum mass of 7.1 Me. An additional signal in the periodogram of the residuals is present but we cannot claim it to be related to a second planet with sufficient significance at this point. If confirmed, such planet candidate would have a minimum mass of 6.3 Me and a period of 10.4 days, which might indicate a 3:2 mean-motion resonance with the inner planet

Characterization of the K2-38 planetary system: Unraveling one of the densest planets known to date

.-- Toledo-Padrón, B. et al.Context. An accurate characterization of the known exoplanet population is key to understanding the origin and evolution of planetary systems. Determining true planetary masses through the radial velocity (RV) method is expected to experience a great improvement thanks to the availability of ultra-stable echelle spectrographs. Aims. We took advantage of the extreme precision of the new-generation echelle spectrograph ESPRESSO to characterize the transiting planetary system orbiting the G2V star K2-38 located at 194 pc from the Sun with V 11.4. This system is particularly interesting because it could contain the densest planet detected to date. Methods. We carried out a photometric analysis of the available K2 photometric light curve of this star to measure the radius of its two known planets, K2-38b and K2-38c, with Pb = 4.01593 ± 0.00050 d and Pc = 10.56103 ± 0.00090 d, respectively. Using 43 ESPRESSO high-precision RV measurements taken over the course of 8 months along with the 14 previously published HIRES RV measurements, we modeled the orbits of the two planets through a Markov chain Monte Carlo analysis, significantly improving their mass measurements. Results. Using ESPRESSO spectra, we derived the stellar parameters, Teff = 5731 ± 66, log g = 4.38 ± 0.11 dex, and [Fe/H] = 0.26 ± 0.05 dex, and thus the mass and radius of K2-38, Ma = 1.03-0.02+0.04 MaS and Ra = 1.06-0.06+0.09 RaS. We determine new values for the planetary properties of both planets. We characterize K2-38b as a super-Earth with RP = 1.54 ± 0.14 RaS and Mp = 7.3-1.0+1.1 MaS, and K2-38c as a sub-Neptune with RP = 2.29 ± 0.26 RaS and Mp = 8.3-1.3+1.3 MaS. Combining the radius and mass measurements, we derived a mean density of ρp = 11.0-2.8+4.1 g cm-3 for K2-38b and ρp = 3.8-1.1+1.8 g cm-3 for K2-38c, confirming K2-38b as one of the densest planets known to date. Conclusions. The best description for the composition of K2-38b comes from an iron-rich Mercury-like model, while K2-38c is better described by a rocky-model with H2 envelope. The maximum collision stripping boundary shows how giant impacts could be the cause for the high density of K2-38b. The irradiation received by each planet places them on opposite sides of the radius valley. We find evidence of a long-period signal in the RV time-series whose origin could be linked to a 0.25-3 MJ planet or stellar activity.With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737

HADES RV Programme with HARPS-N at TNG XII. The abundance signature of M dwarf stars with planets

Most of our current knowledge on planet formation is still based on the analysis of main-sequence, solar-type stars. Conversely, detailed chemical studies of large samples of M-dwarf planet hosts are still missing. We develop for the first time a methodology to determine stellar abundances of elements others than iron for M dwarf stars from high-resolution, optical spectra. Our methodology is based on the use of principal component analysis and sparse Bayesian's methods. We made use of a set of M dwarfs orbiting around an FGK primary with known abundances to train our methods. We applied our methods to derive stellar metalliticies and abundances of a large sample of M dwarfs observed within the framework of current radial velocity surveys. We then used a sample of nearby FGK stars to cross-validate our technique by comparing the derived abundance trends in the M dwarf sample with those found on the FGK stars. The metallicity distribution of the different subsamples shows that M dwarfs hosting giant planets show a planet-metallicity correlation as well as a correlation with the stellar mass. M dwarfs hosting low-mass planets do not seem to follow the planet-metallicity correlation. We also found that the frequency of low-mass planets does not depend on the mass of the stellar host. These results seem in agreement with previous works. However, we note that for giant planet hosts our metallicities predict a weaker planet metallicity correlation but a stronger mass-dependency than photometric values. We show, for the first time, that there seems to be no differences in the abundance distribution of elements different from iron between M dwarfs with and without known planets. Our data shows that low-mass stars with planets follow the same metallicity, mass, and abundance trends than their FGK counterparts.Comment: Accepted for publication by Astronomy & Astrophysic

The HADES RV Programme with HARPS-N@TNG VIII. Gl15A: A multiple wide planetary system sculpted by binary interaction

We present 20 years of radial velocity (RV) measurements of the M1 dwarf Gl15A, combining 5 years of intensive RV monitoring with the HARPS-N spectrograph with 15 years of archival HIRES/Keck RV data. We carry out an MCMC-based analysis of the RV time series, inclusive of Gaussian Process (GP) approach to the description of stellar activity induced RV variations. Our analysis confirms the Keplerian nature and refines the orbital solution for the 11.44-day period super Earth, Gl15A\,b, reducing its amplitude to $1.68^{+0.17}_{-0.18}$ m s$^{-1}$ ($M \sin i = 3.03^{+0.46}_{-0.44}$ M$_\oplus$), and successfully models a long-term trend in the combined RV dataset in terms of a Keplerian orbit with a period around 7600 days and an amplitude of $2.5^{+1.3}_{-1.0}$ m s$^{-1}$, corresponding to a super-Neptune mass ($M \sin i = 36^{+25}_{-18}$ M$_\oplus$) planetary companion. We also discuss the present orbital configuration of Gl15A planetary system in terms of the possible outcomes of Lidov-Kozai interactions with the wide-separation companion Gl15B in a suite of detailed numerical simulations. In order to improve the results of the dynamical analysis, we derive a new orbital solution for the binary system, combining our RV measurements with astrometric data from the WDS catalogue. The eccentric Lidov-Kozai analysis shows the strong influence of Gl15B on the Gl15A planetary system, which can produce orbits compatible with the observed configuration for initial inclinations of the planetary system between $75^\circ$ and $90^\circ$, and can also enhance the eccentricity of the outer planet well above the observed value, even resulting in orbital instability, for inclinations around $0^\circ$ and $15^\circ - 30^\circ$. The Gl15A system is the multi-planet system closest to Earth, at $3.57$ pc, and hosts the longest period RV sub-jovian mass planet discovered so far.Comment: 19 pages, 14 figures, accepted for publication in A\&

HADES RV Programme with HARPS-N at TNG XIII. A sub-Neptune around the M dwarf GJ 720 A

Context. The high number of super-Earth and Earth-like planets in the habitable zone (HZ) detected around M-dwarf stars in the last years has revealed these stellar objects to be the key for planetary radial velocity (RV) searches. Aims. Using the HARPS-N spectrograph within The HArps-n red Dwarf Exoplanet Survey (HADES) we reach the precision needed to detect small planets with a few Earth masses using the RV technique. Methods. We obtained 138 HARPS-N RV measurements between 2013 May and 2020 September of GJ 720 A, classified as an M0.5V star located at a distance of 15.56 pc. To characterize the stellar variability and to discern the periodic variation due to the Keplerian signals from those related to stellar activity, the HARPS-N spectroscopic activity indicators and the simultaneous photometric observations were analyzed. The combined analysis of HARPS-N RVs and activity indicators let us to address the nature of the periodic signals. The final model and the orbital planetary parameters were obtained by fitting simultaneously the stellar variability and the Keplerian signal using a Gaussian process regression and following a Bayesian criterion. Results. The HARPS-N RV periodic signals around 40 d and 100 d have counterparts at the same frequencies in HARPS-N activity indicators and photometric light curves. Then we attribute these periodicities to stellar activity the former period being likely associated with the stellar rotation. GJ 720 A shows the most significant signal at 19.466$\pm$0.005 d with no counterparts in any stellar activity indices. We hence ascribe this RV signal, having a semiamplitude of 4.72$\pm$0.27 m/s , to the presence of a sub-Neptune mass planet. The planet GJ 720 Ab has a minimum mass of 13.64$\pm$0.79 M$_{\oplus}$, it is in circular orbit at 0.119$\pm$0.002 AU from its parent star, and lies inside the inner boundary of the HZ around its parent star