Free Access
Issue
A&A
Volume 625, May 2019
Article Number A126
Number of page(s) 17
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/201834969
Published online 24 May 2019

© ESO 2019

1 Introduction

Most of the early surveys hunting for exoplanets, which employed the radial-velocity (RV) method, directed their efforts towards dwarf stars of spectral type G or K, usually in a range around the mass of the Sun (e.g. Queloz et al. 2001; Valenti & Fischer 2005; Tamuz et al. 2008). Instead, in recent years M dwarfs have become the most promising targets for the hunt for low-mass rocky planets (e.g. Bonfils et al. 2013; Dressing & Charbonneau 2013; Sozzetti et al. 2013; Astudillo-Defru et al. 2017a) due to their more advantageous mass and radius ratios compared to solar-type stars. Moreover, with the availability of high-precision spectrographs mounted on 4 m-class telescopes, led by HARPS at La Silla (Northern High Accuracy Radial velocity Planet Searcher, Mayor et al. 2003) and recentlybacked up by its younger twin HARPS-N at TNG (Cosentino et al. 2012), it was possible to reach 1 m s−1 precision allowing the detection of small Earth-like rocky planets (e.g. Anglada-Escudé et al. 2016a; Astudillo-Defru et al. 2017b).

It is becoming clear that giant gas planets are less frequent around low-mass stars than around those of solar type, as expected from theoretical studies (e.g. Laughlin et al. 2004; Mordasini et al. 2009), while low-mass rocky planets are proving to be much more common around low-mass stars, both in RVs (e.g. Tuomi et al. 2014, and references therein) and transit (e.g. Gaidos et al. 2016, and references therein) observations. Nevertheless, the complete characterization of this abundant population of rocky planets around M dwarfs is hindered by the strong effects of the photospheric and magnetic activity of the latter, which can produce RV signals as large as tens of meters per second. This can result in stellar signals being mistaken for planetary signals or otherwise uncertain results for systems around active stars (e.g. Bonfils et al. 2007; Baluev 2013a; Robertson et al. 2014), and also around quieter targets, which can still present periodic signals of unclear nature (e.g. Robertson et al. 2015; Anglada-Escudé et al. 2016b).

The Harps-n red Dwarf Exoplanet Survey (HADES) programme is a collaboration between the Italian Global Architecture of Planetary Systems (GAPS, Covino et al. 2013; Desidera et al. 2013; Poretti et al. 2016) Consortium1, the Institut de Ciències de l’Espai de Catalunya (ICE), and the Instituto de Astrofísica de Canarias (IAC). The aim of the survey is to characterise the exoplanetary systems around a well-defined sample of M dwarfs, with spectral type between dM0 and dM3. High-precision RVs of the sample have been collected over the course of five years with the HARPS-N at TNG spectrograph. Several planets have been discovered as part of the survey (e.g. Affer et al. 2016; Suárez Mascareño et al. 2017; Perger et al. 2017a; Pinamonti et al. 2018), and also some more general studies have already been performed on the samples, studying both the stellar properties (Maldonado et al. 2017; Scandariato et al. 2017; Suárez Mascareño et al. 2018; González-Álvarez et al. 2019) and the preliminary planetary population statistics (Perger et al. 2017b). Moreover, a collaboration between the HADES and CARMENES programs resulted in the detection of the super-Earth Gl 49 b (Perger et al. 2019).

In this work we present the search for planetary companions around the M dwarf GJ 685, based on high-precision spectroscopic observations carried out with HARPS-N as part of the HADES programme. We also take advantage of ancillary photometricobservations of the target to better constrain the stellar activity signal in the RVs.

In Sect. 2 we describe the Doppler measurements of GJ 685 collected for this analysis, and in Sect. 3 we briefly discuss the physical properties of the host star. The independent analyses of two photometric datasets are presented in Sect. 4. We describe our periodogram analyses of the RV data and stellar activity indexes in Sect. 5. We proceed to find the best-fit parameters for the models describing the activity and RV time series via a Monte Carlo Markov chain (MCMC) analysis in Sect. 6. Finally, we summarise and discuss our findings in the more general context of the current population of M-dwarf RV-detected planetary systems in Sect. 7.

thumbnail Fig. 1

HARPS-N RV time series of GJ 685.

2 Spectroscopic observations

As part of the HADES RV programme, GJ 685 has been observed from BJD = 2 456 439.6 (27 May 2013) to BJD = 2 458 044.4 (17 October 2017) with the HARPS-N spectrograph, connected by fibres to the Nasmyth B focus through a Front End Unit of the 3.58 m Telescopio Nazionale Galileo (TNG) in La Palma, Spain. HARPS-N is a fibre-fed, cross-dispersed echelle spectrograph with a spectral resolution of 115 000, covering a wavelength range from 3830 to 6900 Å. We observed with fixed integration times of 900 s to obtain data of sufficient signal-to-noise ratio (SN > 20) and to average out short-term periodic oscillations of the star, such as p-modes (Dumusque et al. 2011).

The total number of data points acquired was 106 over a time span of 1605 days. The time series is shown in Fig. 1. The observations were gathered without the simultaneous Th-Ar calibration, which could contaminate the Ca II H and K lines due to the long exposure times and the relative faintness of M-dwarf targets in the blue part of the spectra. Moreover, these lines are crucial in the analysis of stellar activity (e.g. Giampapa et al. 1989; Forveille et al. 2009), which is particularly important for active late-type stars. A precise acquisition of the Ca II H and K lines was therefore preferred over a better correction of possible instrumental drifts. Perger et al. (2017b) used other GAPS target spectra gathered by the Italian team on the same nights using the simultaneous Th-Ar calibration to quantify this effect, and found a mean inter-night instrumental drift of the order of 1 m s−1 over the whole HADES sample. Nevertheless, for some of the brightest targets of the sample, some spectra could be safely taken with the simultaneous drift calibrations without problems for the Ca II H and K observations (Perger et al. 2017a). This helped us to monitor the instrumental drift over the single target: for GJ 685, four spectra were collected with the simultaneous Th-Ar calibrations, measuring a mean inter-night instrumental drift of ~ 0.5 m s−1. Any residual drift in the time series was taken into account in our final model as discussed in Sect. 6.2.

The data reduction and the RV extraction were performed using the HARPS-N Data Reduction Software (DRS, Lovis & Pepe 2007) and the TERRA pipeline (Template-Enhanced Radial velocity Re-analysis Application, Anglada-Escudé & Butler 2012), respectively. The latter is considered to be more accurate when applied to M-dwarfs, with respect to the DRS. For a more thorough discussion of the DRS and TERRA performances on the HADES targets see Perger et al. (2017b). The mean internal error of the TERRA data is 1.02 m s−1, with a few low-S/N data with σRV > 2.0 m s−1. With an rms of 6.16 m s−1, GJ 685 is one of the HADES targets with the largest RV dispersion. The TERRA pipeline also corrected the RV data for the perspective acceleration of GJ 685, d vr ∕dt = 0.11 m s−1 yr−1.

Table 1

Stellar parameters for the target GJ 685.

3 Stellar properties of GJ 685

The star GJ 685 is a high proper motion nearby (π = 69.825 ± 0.039 mas) M0-type dwarf. We used the stellar parameters published by Maldonado et al. (2017), which were calculated applying the empirical relations by Maldonado et al. (2015) to the same HARPS-N spectra from which we derived the RV time series. The parallax and proper motions of the star were taken from the Gaia Data Release 2 (Gaia Collaboration 2018). All the stellar parameters of GJ 685 are listed in Table 1.

Suárez Mascareño et al. (2018) studied the presence of signatures of magnetic cycles and rotation on the stars of the HADES sample, measuring the rotation periods and logRHK$\log R'_{\text{HK}}$ for several stars of the sample2. For GJ 685, Suárez Mascareño et al. (2018) derived, from the variability in the S-index and Hα activity indexes and RV time series, a rotation period of 16.3 ± 4.2 d finding no evidence for the presence of a magnetic cycle. The rotation period value is listed in Table 1 along with the measured value of logRHK$\log R'_{\text{HK}}$.

4 Photometric monitoring

As with most of the targets of the HADES sample, GJ 685 has been monitored photometrically by means of the APACHE (Sozzetti et al. 2013) and EXORAP (EXOplanetary systems Robotic APT2 Photometry) surveys. The two surveys perform regular follow-up observations of HADES M-dwarf targets to constrain the stellar rotation periods by analysing the photometricvariability. We briefly discuss the analyses of the datasets collected by the two surveys in the following sections.

4.1 APACHE photometry

GJ 685 was monitored for 64 nights between BJD = 2 456 456.4 (12 Jun 2013) and BJD = 2 456 793.6 (15 May 2014) with one of the five 40 cm telescopes composing the APACHE array located at the Astronomical Observatory of the Autonomous Region of the Aosta Valley (OAVdA, + 45.7895 N, + 7.478 E, 1650 m.a.s.l.). The observations were collected following the standard APACHE procedure, and the images were reduced with the standard pipeline TEEPEE by the APACHE team (Giacobbe et al. 2012).

Since Suárez Mascareño et al. (2018) calculated the rotation period of GJ 685 to be Prot = 16.3 ± 4.2 d, we performed Generalized Lomb Scargle periodogram analysis (GLS, Zechmeister & Kürster 2009) of the photometric data, looking for similar periodicities. To do so we binned the data over each night to average out the short-period noise due to the high number of very close data points. This resulted in a time series of 64 data points over a time span of 337 d, with an rms of 0.012 mag.

The results of the GLS analysis, which covered periods between 1 d and the time span of the time series, are shown in Fig. 2, with the highest peak at a period of Prot, AP = 16.85 ± 0.12 d with a theoretical false alarm probability (FAP) of 1.9%, in good agreement with the rotation period found by Suárez Mascareño et al. (2018). It is also worth noticing that there is no significant peak at longer periods.

4.2 EXORAP photometry

In the framework of the EXORAP project, we observed GJ 685 using an 80 cm f∕8 Ritchey-Chretien robotic telescope (APT2) located at Serra la Nave on the Mt. Etna and operated by the INAF-Catania Astrophysical Observatory. We collected approximately 200 measurements in each band between 5 May 2014 and 6 Sept 2017. Data reduction was preformed with overscan, bias, dark subtraction, and flat fielding with IRAF procedures and visually inspected tocheck the data quality (see Affer et al. 2016 for details).

The scatter of the B and V photometry is slightly larger than 0.01 mag, which corresponds to the intra-night sensitivity of the survey. Thissuggests that there is some jitter of stellar origin in the collected data. The pooled variance (PV) analysis (Scandariato et al. 2017, and references therein) identifies a significant timescale around 20 d, but the precision of the PV technique is not sufficient to distinguish the 18 and 24 d periodicities that we discuss in Sect. 5.1. The GLS periodogram analysis shows a low-significance peak again around 20 d, while the GP analysis converges directly to an 18 d period, excluding longer periods closer to 24 d.

The analysis of the RI photometry does not lead to any significant result, as the scatter of the data of ~ 0.01 mag is dominated by the inter-night sensitivity of the survey. This is consistent with a scenario where the photometric scatter is dominated by cool photospheric spots, whose contrast against the unspotted photosphere is larger in the bluer bands than in R and I. This is also consistent with other similar analyses we published in other papers of the HADES series.

thumbnail Fig. 2

Upperpanel: GLS periodograms of the APACHE photometric data. The dotted and dashed horizontal lines indicate the 1 and 10% FAP levels, respectively. Lower panel: APACHE light curve phase-folded over the 16.85 d period found by the GLS periodogram. The red line represents the best-fit model.

5 Periodogram analyses

We started by analysing our spectroscopic data by means of GLS periodograms in order to identify significant periodicities in our time series, and compare the signals identified in the RV and activity indexes time series to pinpoint activity-related signals in the RV time series.

5.1 Radial velocity periodograms

First we analysed the RV time series, identifying additional significant periodicities by pre-whitening until no significant signal was found in the GLS periodogram below the FAP = 10% level. We computed the GLS FAPs via bootstrap randomisation with 10 000 iterations (Endl et al. 2001). The resulting periodograms are shown in Fig. 3. The strongest peak in the first periodogram is at P = 9 d, which is probably due to stellar activity, being roughly half the rotation period of the star Prot = 16.3 ± 4.2 d, as derived by Suárez Mascareño et al. (2018) and confirmed by our photometric analysis in the previous section. After subtracting this signal we see that a peak around P = 23.66 d rises from the periodogram, along with a clustered peak around P = 18 d. The latter can be related again to the activity signal of the star, proving that a simple sinusoidal fit with P = P1∕2 = 9 d is not sufficient to model the influence of active regions on the RVs. We thus identify Prot = 18 d as the rotation period of the star. It is worth noting that this rotation period is longer than the one derived from the analysis of the APACHE photometry in Sect. 4.1. This can be explained as a consequence of two effects: firstly, activity signals in the photometry can differ from those present in the RV time series, depending on the nature of the active regions present on the star (e.g. Kürster et al. 2003; Dumusque et al. 2014); secondly, the rotation signal is relatively weak in the photometry, with a theoretical FAP > 1% in the APACHE analysis and not appearing clearly in the EXORAP datasets, and thus a relatively large uncertainty in the rotation period is not surprising. It is also worth noting that a peak at P = 18 d is indeed present in the periodogram in Fig. 2, even if weaker than the main peak. The P = 23.7 d signal instead does not appear to be easily related to the assumed rotation period of the star or its harmonics, and it is also relatively strong in the periodogram (FAP = 0.47%). We therefore suspect it to be a genuine Keplerian Doppler shift due to an orbiting planet, hereafter GJ 685 b. The semi-amplitude of the signal at P = 23.661 ± 0.037 d identified in the periodogram is K = 3.11 ± 0.53 m s−1.

For a more comprehensive view of the periodicities present in the RV data, we also studied the time series with the Bayesian Generalized Lomb-Scargle (BGLS, Mortier et al. 2015) and FREquency DEComposer (FREDEC, Baluev 2013b) to compare their different results as suggested in Pinamonti et al. (2017). The BGLS periodogram results are very similar to those from GLS, with the 9 d peak dominating the first periodogram, and the P = 23.7 d signal emerging in the residuals analysis. A difference arises in the second residual analysis compared to the bottom panel of Fig. 3, since the periodic signals are much weaker, with a theoretical BGLS FAP (as defined in Pinamonti et al. 2017) greater than 10%3. The FREDEC analysis instead produces as best fit solution the five-signal solution P1 = 18.46 d, P2 = 18.10 d, P3 = 9.66 d, P4 = 9.05 d, P5 = 8.99 d, with a FAP = 1%. As an alternative but slightly less significant (FAP ≃ 1.5%) solution, the algorithm proposes a four-signal solution P1 = 23.71 d, P2 = 17.03 d, P3 = 9.65 d, P4 = 9.04 d, which includes also the 23.7 d signal found by the other techniques. It is worth noting the presence of multiple periodic signals near the stellar rotation period and its first harmonic, which proves that the RV stellar activity signal is strongly quasi-periodic.

thumbnail Fig. 3

GLS periodograms of the RV time series and residuals after consecutive signal fits. The red vertical lines indicate the potential rotation period (dashed) and its first harmonic (dot-dashed), while the dashed blue vertical lines marks the orbital period of the planetary candidate. The dotted, dash-dotted, and dashed horizontal lines indicate the 0.1, 1, and 10% FAP levels, respectively.

5.2 Stellar activity analysis

To expand the analysis of the stellar activity signals of GJ 685 performed by Suárez Mascareño et al. (2018), we derived the complete time series for all the line profile indicators evaluated with the method from Lanza et al. (2018). These latter authors derived several line profile asymmetry indicators by computing the cross-correlation function (CCF) between a mask and the stellar spectra. For our study, we selected three asymmetry indicators: the bisector inverse span (BIS), ΔV (which compute thedifference between the RV values computed with a Gaussian and bi-Gaussian best fit of the CCF, respectively), and Vasy (which quantifies the asymmetry in the RV spectral line information content). In addition to this asymmetry analysis, we also derived the activity indexes based on the stellar Ca II H and K, Hα, Na I D1 D2, and He I D3 spectral lines following the procedure described in Gomes da Silva et al. (2011).

As a first-order test of the effect of stellar chromospheric activity on the RV time series, we checked for correlations between the asymmetry and activity indexes, andthe RV datasets. We computed the Pearson correlation coefficients for the different combinations of RV and activity indicators, and no significant correlation was identified (|ρ|0.3$\left|\rho\right|\,{\lesssim}\,0.3$ for all the indexes).

The GLS periodograms of these asymmetry and activity indexes are shown in Figs. 4 and 5. It is worth noting that for the asymmetry indicator Vasy we show the results obtained from the new definition by Lanza et al. (2018) (Vasy(mod)), since, as they stated, the original definition from Figueira et al. (2013) is sensitive to genuine Doppler shifts of the star, and thus presents misleading signals at the periods of actual planets, leading to erroneous rejections.

We can see in Fig. 4 that the asymmetry indicators show only weak periodic signals, which for two of them, BIS and ΔV, correspond to the first harmonic of the stellar rotation period. In Fig. 5 we see instead that the Ca II H and K and Hα present strong signals corresponding to the stellar rotation, while the He I D3 index shows only a long-term signal around 400 d4. The GLS analysis of the Na I D1 D2 time series produces no significant signals.

Our GLS analysis of different asymmetry and activity indicators confirms the results of Suárez Mascareño et al. (2018), providing a slightly longer stellar rotation period of 18 d, with clear signals at both the rotation period and its first harmonic. This confirms the rotation period, Prot = 18 d, identified inthe RV time series in Sect. 5.1, with respect to the shorter Prot, AP = 16.85 d found in the APACHE photometric data in Sect. 4.1. Moreover it is worth noting how no peak was identified in the indicators time series at periods corresponding or close to the Pb = 23 d period of the new planet candidate GJ 685 b.

thumbnail Fig. 4

GLS periodograms of the asymmetry indicators. The red vertical lines indicate the potential rotation period (dashed) and its first harmonic (dot-dashed), while the dashed blue vertical lines mark the orbital period of the planetary candidate discussed in Sect. 5.1. The dotted horizontal lines indicate the 1% FAP level.

thumbnail Fig. 5

GLS periodograms of the activity indexes. The red vertical lines indicate the potential rotation period (dashed) and its first harmonic (dot-dashed), while the dashed blue vertical lines mark the orbital period of the planetary candidate discussed in Sect. 5.1. The dotted horizontal lines indicate the 1% FAP level.

6 MCMC analysis

We then proceeded to expand our analysis of the RV and activity indexes time series with a combined fit of the Keplerian and stellar activity signals. A very common method to model and subtract the stellar activity correlated “noise” from RV time series is the Gaussian Process (GP) regression (e.g. Haywood et al. 2014; Grunblatt et al. 2015; Dumusque et al. 2017; Pinamonti et al. 2018). This technique has proven to be especially effective when adopting a quasi-periodic covariance function described by four parameters, called hyper-parameters: K(t,t)=h2×exp[(tt)22λ2sin2(π(tt)θ)2w2]+(σdata2(t)+σjit2)δt,t, \begin{eqnarray*}K(t, t^{\prime})\,&{=}&\,h^2\times\exp\bigg[-\frac{(t-t^{\prime})^2}{2\lambda^2} - \frac{\textrm{sin}^{2}\left(\dfrac{\pi(t-t^{\prime})}{\theta}\right)}{2w^2}\bigg] \nonumber \\ && +\, (\sigma^{2}_{\textrm{data}}(t)\,+\,\sigma_{\textrm{jit}}^{2})\cdot\delta_{\textrm{t, t}^{\prime}}, \end{eqnarray*}(1)

where t and t′ indicate two different epochs; h is the amplitude of the correlations; θ represents the period of the correlated signal (and corresponds to the rotation period of the star in our model); w is the length scale of the periodic component; and λ is the correlation decay timescale (which can be related to the decay time of the active regions); σdata(t) is the data internal error at time t for each instrument; σjit is the additional uncorrelated “jitter” term used in the analysis of the RVs; and δt,t$\delta_{\textrm{t, t}^{\prime}}$ is the Kronecker delta function.

For a more thorough description of the GP kernel and hyper-parameters, see Pinamonti et al. (2018).

We applied the GP regression as part of an MCMC analysis, performed via the publicly available emcee algorithm (Foreman-Mackey et al. 2013) and GEORGE Python library (Ambikasaran et al. 2015). We used 150 random walkers to sample the parameter space. The posterior distributions were derived after applying a burn-in as explained in Eastman et al. (2013, and references therein). To evaluate the convergence of the different MCMC analyses, we calculated the integrated correlation time for each of the parameters, and stopped the code after a number of steps equal to 150 times the largest autocorrelation times of all the parameters (Foreman-Mackey et al. 2013).

Table 2

Priors and best-fit results for the Gaussian process regression analysis of the Hα activity indicator.

6.1 Activity indexes GP analysis

First we performed a GP emcee fit of the activity index time series, as a reference for the following GP+planetary signal analysis of the RV time series. We show only the Hα time series because, as shown in Fig. 5, the stellar rotation signal of the star is clearly present in the time series, and some studies also suggest the Hα to be the best indicators of the activity of early- and mid-M type stars (e.g. Robertson et al. 2013). We nonetheless also performed an emcee analysis of the Ca II H and K time series (not shown), which produced analogous results.

We chose uniform priors for all the GP hyper-parameters. For the λ hyper-parameter, which represents the correlation decay timescale, and can range over several orders of magnitude, we adopted a uniform prior in logarithmic scale to avoid oversampling the long scales. The adopted priors and best-fit results for the GP hyper-parameters of the analysis are listed in Table 2, while the a posteriori distributions are shown in Fig. 6.

The GP fit confirms that the rotation period of the star, corresponding to the hyper-parameter θ, is close to ~18 d, as inferred from the GLS analysis. It is also worth noting that the median value of the λ hyper-parameter, which is related to the evolution timescale of the active regions, is of the order of three to four times the rotation period. This is consistent with an independent analysis of a few M dwarfs in the HADES sample reported in Scandariato et al. (2017). The uncertainties of this hyper-parameter are very large due to the tail of the distribution at short timescales, which can be observed in Fig. 6. However, this uncertainty does not affect the determination of the rotation period θ, which is well constrained within a 1 d interval. In addition, this analysis does not present any evidence of longer-period signals in the activity indexes of GJ 685, in particular around the candidate planet period Pb = 24 d.

6.2 emcee analysis of the RV time series

To better understand the structure of the stellar activity signal in the RV time series we performed a GP analysis of the dataset, adopting very similar priors to those adopted in the analysis of the Hα activity index time series (see Table 2), and also adding an uncorrelated jitter term to the model. In Table 3 we can see the chosen priors and the best fit values. Since from the analysis of the Hα time series no evidence emerged of long correlation decay timescales, we restricted the prior of λ in the interval [1, 500] d, again uniform in logarithmic space. Even though the analysis of the activity and asymmetry indexes in the previous sections pointed out the absence of stellar activity signals at periods larger than 20 d and in particular close to the orbital period of the planet candidate, we decided to keep the prior of the rotation period θ over the interval [5, 30] d, as in the analysis described in the previous section. In this analysis, we add an uncorrelated jitter term σjit, which takes into account any additional uncorrelated stellar noise that is not corrected by the GP model. This can also be used to take into account possible residual errors from the instrumental drift correction (see Sect. 2) as done in Affer et al. (2019).

In Fig. 7 the posterior distributions of the fit parameters are shown. We can see that, even if the GLS periodogram identifies P = 9 d as the strongest period, the GP correctly identifies the 18 d rotation period of the star, taking into account the complex nature of the stellar activity RV signal due to its quasi-periodic nature. It is also worth noting from Table 3 that the value found in this analysis for the hyper-parameter λ is smaller than the one previously found in the analysis of the Hα index time series, even though it remains consistent within the large error bars.

We then repeated the emcee analysis adding a Keplerian planetary signal to the model, to recover the best-fit parameters of the candidate planet GJ 685 b. We used wide priors for the planetary parameters, so as not to force the solution on the value found with the GLS periodogram, with the prior of the orbital period ranging up to 100 d. For the same reason, the priors on the planetary parameters were all chosen to be uninformative uniform priors. We can see the best-fit solution in the right column of Table 3, while the posterior distributions are shown in Fig. 8. It is worth noting that, while the amplitude of the signal is quite similar to the value recovered from the periodogram analysis, the period is slightly longer, due to the simultaneous fitting of the stellar activity signal with the GP. A series of aliases of the orbital period P appear to be present, slightly correlated with T0, in the bottom right contour plot in Fig. 8. However, as we can see in the histograms of the posterior distributions of the orbital parameters, these aliases are not significant and do not change the best-fit value of the orbital period. Moreover, we can see in the right column of Table 3 that the best-fit value of λ is larger than in the pure GP analysis (left column), and much more similar to that found in the analysis of the Hα: this can be explained by the fact that in the pure-GP model the presence of the un-modelled planetary signal interferes with the fit of the stellar activity, and a shorter decay timescale is needed to compensate this effect. Once the planetary signal is taken into account in the analysis, the decay timescale returns to the more accurate value found in the activity index analysis. It is also worth noting that the jitter term retrieved by the emcee analysis is small, σjit=1.460.32+0.33$\sigma_{\text{jit}}\,{=}\,1.46^{+0.33}_{-0.32}$ m s−1, suggesting low levels of uncorrelated stellar noise, as well as confirming the low levels of residual instrumental drift as discussed in Sect. 2.

The quasi-periodic stellar model obtained from the simultaneous GP + 1 planet fit is shown in Fig. 9 compared to the RV time series residuals after the subtraction of the 24 d planetary signal. We can see both the fine correspondence between the data and our stellar model, and the strength and variability of the stellar signals throughout the four years of HADES observations.

As additional evidence of the presence of the candidate planet GJ 685 b, we computed the Bayesian information criterion (BIC, Schwarz et al. 1978) for the two models, obtaining BIC = 635 and BIC = 619 for the pureGP and GP + 1 planet models, respectively. There is a very strong statistical evidence, ΔBIC = 16, in favour of the presence of the planetary signal. Figure 10 shows the phase-folded RV time series, closely following the planetary model after subtraction of the stellar activity signal.

In Table 3 are also shown the values of minimum mass, MP sin i, and semi-major axis, aP, derived from the best-fit orbital parameters of the GP + planet model. The derived minimum mass is 9.01.8+1.7$9.0^{+1.7}_{-1.8}$ M, placing GJ 685 b within the Super-Earth regime. Since single-planet systems are known to show a wide range of eccentricities (e.g. Rodigas & Hinz 2009; Limbach & Turner 2015), we also tested an eccentric model for GJ 685 b in order to constrain the possible eccentricity of the orbit. This additional analysis resulted in a best-fit value of e=0.140.10+0.18$e\,{=}\,0.14 ^{+0.18} _{- 0.10}$, consistent with zero within 1.5-σ, and also with a higher value of BIC = 628 with respectto the circular model.

We found no evidence for the presence of additional short-period signals in the RV time series: we computed the GLS periodogram of the RV residuals after the subtraction of the GP + 1 planet model, and no significant signal below the 10% FAP level was found. Additionally, we tested a GP + 2 planets MCMC model, with the parameters of the second planet left free to explore a wide parameter space; we did not find any dominant signal and no statistical improvement of its BIC over the GP + 1 planet model. Moreover, as reported in Table 3, the best-fit value of the acceleration in our final model is 0.0010 ± 0.0024 m s−1 d−1, largely consistent with zero, suggesting the absence of long-period signals.

Table 3

Priors and best-fit results for the emcee analysis of GJ 685 RV time series.

thumbnail Fig. 6

Posterior distributions of the fitted (hyper-)parameters of the GP quasi-periodic model applied to the time series of Hα activity index. The vertical dashed lines denote the median and the 16th and 84th percentiles.

thumbnail Fig. 7

Posterior distributions of the fitted (hyper-)parameters of the pure GP model applied to the RV time series. The vertical dashed lines denote the median and the 16 th and 84 th percentiles.

7 Summary and discussion

We investigated 106 spectroscopic observations of GJ 685 obtained over 4.4 yr with HARPS-N at the TNG in La Palma and additional photometry from the APACHE and EXORAP programs. We used RVs derived from the TERRA pipeline along with activity and asymmetry indexes derived from the same HARPS-N spectra and used them to monitor the stellar chromospheric activity of the target.

The RV time series of GJ 685 is dominated by three peaks at 9, 18, and 24 d. Our spectroscopic analysis, strengthened by the analyses of two independent photometric light curves of the target, confirms that the 18 and 9d period signals are related to the stellar activity, corresponding respectively to the stellar rotation period and its first harmonic. On the other hand, the 24 d period signal seems not to be related to any stellar effects, and is best described as a Keplerian signal caused by an orbiting planet, GJ 685 b.

To derive the minimum mass and orbital parameters of GJ 685 b we fitted the RV time series with a Keplerian model combined with a GP quasi-periodic model to take into account the stellar activity signal. We obtained a period Pb = 24.160 d, a semi-major axis a = 0.1344 AU, and a minimum mass Mb sin i = 9.0 M. The GP quasi-periodic model improves the precision of the rotation period of the stars computed by Suárez Mascareño et al. (2018), finding a best-fit value of 18.150.16+0.15$18.15^{+0.15}_{-0.16}$ d. The amplitude of the stellar activity signal is h=6.050.72+0.94$h\,{=}\,6.05^{+0.94}_{-0.72}$ m s−1, more than twice the amplitude of the Keplerian signal of GJ 685 b, similar to the case of GJ 3942 (Perger et al. 2017a). It is also worth noting that GJ 685 presents the largest stellar RV signal of the HADES targets with planetary companions detected to date. Moreover, even if the strongest periodic signal present in the RV time series is the rotation period first harmonic, Prot∕2 = 9 d, the GP model easily identifies the actual rotation period of the star, Prot = 18 d, as the source of the RV modulation. Even though the prior adopted for the θ hyper-parameter was very wide, including also the first harmonic value of 9 d, the model converged naturally on the rotation period. This again demonstrates the effectiveness of GP quasi-periodic models when dealing with complex stellar activity signals producing several different peaks in the periodogram, even if the strongest peak does not correspond to the stellar rotation period.

We tested an eccentric model for the orbit of GJ 685 b. A single-planet system could be expected to show significant eccentricity, since the orbit of this planet should not be circularized by tides owing to the relatively large separation from the star. The analysis resulted in a best-fit value of eccentricity consistent with zero within 1.5-σ, and a posterior distribution (not shown) strongly peaked at zero. Moreover, the BIC presented strong statistical evidence in favour of the circular-orbit model. We therefore adopted the null-eccentricity model as the best representation of the orbit of GJ 685 b.

Due to its close orbit to the host star, GJ 685 b is unlikely to host an atmosphere capable of maintaining liquid water on its surface. Following the definition of habitable zone (HZ) from Kopparapu et al. (2013), we computed the inner edge of the HZ for GJ 685, with the most optimistic limits (“recent Venus”), which correspond to a semi-major axis of aHZ = 0.190 AU, significantly larger than the orbit of GJ 685 b.

GJ 685 b is the seventh extrasolar planet discovered within the HADES program (Affer et al. 2016, 2019; Suárez Mascareño et al. 2017; Perger et al. 2017a; Pinamonti et al. 2018). Figures 11 and 12 show GJ 685 b and the other HADES planets compared to the current population of RV detected planets orbiting M dwarfs selected from the NASA Exoplanet Archive5 : the sample is composed of70 extrasolar planets discovered with the RV method and orbiting stars of spectral type M0 and M9. Even if the sample is relatively small, some information can be gleaned concerning the planetary and stellar parameters of these systems. First of all, it is worth noting that there seems to be no high-mass planets orbiting low-mass late-type M dwarfs (bottom-left panel Fig. 11). Similarly, there seems to be a scarcity of long-period planets orbiting stars with M* < 0.3 M, which is thesame threshold below which only low-mass planets have been detected. However, these effects are strongly affected by observational bias, since late-type M dwarfs are very faint and difficult to observe, and therefore are usually excluded from RV exoplanet surveys; for example the HADES sample has a median stellar mass of 0.5 M, with no target below0.3 M (Perger et al. 2017b). This bias should be solved in the near future, since RV surveys such as CARMENES (Quirrenbach et al. 2014) and SPIRou (Moutou et al. 2017) are intensively monitoring large samples of nearby M dwarfs in search of extrasolar planets; the CARMENES sample for example includes ≃ 120 stars with M* < 0.3 M (Reiners et al. 2018). Nevertheless, to date no such high-mass planet has been announced.

It is also worth noting from the bottom-right panel of Fig. 11 that there appears to be a dependence of the planetary minimum mass on the stellar metallicity: we computed the Pearson correlation coefficient, which found a weak-to-moderate correlation, ρ = 0.30 p-value = 1.5%. A similar dependence is expected from theoretical models (e.g. Mordasini et al. 2012) and was observed also for solar-mass dwarf stars (e.g. Mortier et al. 2012; Wang & Fischer 2015), while not for planet-hosting giants (Mortier et al. 2013). Moreover, although this effect was observed to be strong for giant planets, which are much more abundant around high-metallicity stars (e.g. Santos et al. 2005; Gaidos & Mann 2014), it is still debated for low-mass Neptune-like planets and super-Earths: some studies found no evidence of a correlation between metallicity and planetary occurrence rates for such small planets (e.g. Sousa et al. 2008; Gaidos et al. 2016), while Wang & Fischer (2015), analysing a large sample of Kepler-candidates, pointed out that a similar correlation should be present, though weakened, down to terrestrial planets. Courcol et al. (2016) studied a sample of mass-measured exoplanets, and found evidence for the frequency of exo-Neptunes (M ∈ [10, 40] M) to be correlated with stellar metallicity, while this was not the case for super-Earths (M < 10 M). Similarly, we could divide our sample of RV-detected planets between masses higher and lower than 10 M (green dashed line in the bottom-right panel of Fig. 11): the two sub-samples have both weaker correlations, ρ = 0.23 p-value = 30% for M > 10 M, and ρ = 0.22 p-value = 14% for M < 10 M. However, the high-mass subsample contains only 22 planets, which could be the cause of the non-detection of the expected correlation. Therefore, based on the selected sample, we are not able to confirm whether or not the trend observed by Courcol et al. (2016) affects M dwarfs in the same way as solar-mass stars.

Analysing only the subsample of exoplanets discovered by the HADES program it is difficult to confirm the properties discussed above. The HADES planets are mostly found in a small region of the parameter space, since they usually have low masses and relatively short periods. Also, by construction of the survey, HADES targets are in narrow ranges of stellar parameters (Perger et al. 2017b; Maldonado et al. 2017). This will allow a focused analysis of the characteristics of planetary systems around specific stellar hosts, and more results will become available as the analysis of the complete survey sample draws near. It is worth noting that the distribution of HADES planets is not surprising, since they are mostly found near the medians of the overall period and minimum mass distributions, P˜=19$\tilde{P}\,{=}\,19$ d M˜sini=7.6$\tilde{M} \sin i\,{=}\,7.6$ M.

Recently, Luque et al. (2018) suggested that masses of single and multiple systems around M dwarfs should follow different distributions. To test this hypothesis, Fig. 12 shows the minimum masses and orbital periods of our sample of M-dwarfs planets, distinguishing between single and multiple planetary systems. We can observe that the two populations appear to have similar distributions, with the single planets being on average slightly more massive than those found in multiple systems6. To asses the statistical significance of this difference, we performed a two-sided Kolmogorov-Smirnov (K-S) test of the minimum mass distributions of the two populations. We obtained a p-value p = 1.3%, thus reinforcing our hypothesis that planets found in multi- and single-planet systems tend to have different masses. This seems to confirm the effect found by Luque et al. (2018), who studied a different sample of M-dwarf extrasolar planets with masses measured from RVs and TTVs. These latter authors suggested two possible explanations, connected to the formation of low-mass planets: that either (i) the formation of super-Earth impedes the formation of smaller Earth-like planets in the same system, or (ii) super-Earth planets around M dwarfs are formed by a pile-up of several low-mass planets. If the former were true, the single more-massive planet populations would have a mass distribution lower than that of the summed mass of the planets in multiple systems, due to the absence of additional smaller-mass planets not able to form. We tested this hypothesis on our dataset, performing another K-S test on the two distributions, and obtained a very high p-value p = 47%. We found no evidence in favour of the first formation mechanism proposed by Luque et al. (2018), suggesting that M-dwarf super-Earths are formed by aggregation of smaller planets.

To further compare the properties of the single- and multi-planet populations around M dwarfs, Fig. 13 shows the eccentricities and orbital periods of the planets in the considered sample. Several studies pointed out a wide distribution of orbital eccentricities of exoplanets, as well as significant correlation with orbital periods (e.g. Stepinski & Black 2000; Kipping 2013) and multiplicities (e.g. Rodigas & Hinz 2009; Limbach & Turner 2015). We can see in Fig. 13 that, even if the observed distribution of exoplanets orbiting solar-type stars extend up to e ~ 0.9, for M-dwarf RV planets the eccentricities are usually e < 0.5, even for long-period giant planets. Moreover, comparing the single- and multi-planet populations, we find no difference in their eccentricity distribution; they are very similar, with a K-S test suggesting the two distributions to be identical (p-value p = 99%). This is further evidence of the difference between the populations of extrasolar planets orbiting M dwarfs and solar-type stars. Focusing on the sub-sample of low-mass planets (Mp sin i < 30 M), it is worth noting that eccentricities are usually poorly constrained: out of 54 low-mass planets, only 7 have eccentricities with a significance higher than 3-σ. This is not surprising, since eccentricities are often difficult to precisely constrain, in particular for single-planet systems (e.g. Wittenmyer et al. 2013).

In Fig. 12, there also appears to be a correlation between the minimum mass and the orbital period, as was suggested by early formation models (Zucker & Mazeh 2002, and references therein). However, the Pearson correlation coefficient does not favour such a correlation in our sample: ρ = 0.22 p-value = 5.0%. It is also important to note that observational biases have a strong influence on the mass and period distributions, since smaller-mass and longer-period planets are more difficult to detect. Recent studies on the occurrence rates of extrasolar planets around M dwarfs, taking into account detection probabilities and detection biases, suggest that low-mass long-period planets could in fact be abundant around M dwarfs (Bonfils et al. 2013; Tuomi et al. 2014), but such an extended analysis is beyond the scope of this paper. Moreover, it is worth noting that M stars are a very heterogeneous group, with early- and late-M dwarfs differing in both their internal structure and magnetic activity. Therefore, a more accurate study of the properties of planetary systems around low-mass stars should consider different early- and late-M stars separately. For this reason, the target of the HADES programme is precisely to study the population of extrasolar planets over a well-defined sample of targets with similar spectral type and stellar properties.

A thorough and unbiased analysis of the detection efficiencies and planetary occurrence rates in the HADES sample, completing and expanding the preliminary statistical analysis from Perger et al. (2017b), will be the object of a future publication (Pinamonti et al., in prep.).

thumbnail Fig. 8

Posterior distributions of the fitted (hyper-)parameters of the GP + 1 planet model applied to the RV time series. The verticaldashed lines denote the median and the 16th and 84th percentiles.

thumbnail Fig. 9

Upper panel: best-fit stellar quasi-periodic signal obtained from the GP + 1 planet model (orange line) compared to the RV residuals (blue points). Lower panel: magnification of the last HARPS-N observing season. The grey area represents the ± 1σ uncertainties of the stellar activity model.

thumbnail Fig. 10

Phase-folded for the RV signal of GJ 685 b, after the subtraction of the stellar correlated signal.

thumbnail Fig. 11

Properties of RV exoplanets orbiting M dwarfs. Top-left panel: orbital period as a function of the mass of the host star. Top-right panel: orbital period as a function of the metallicity of the host star. Bottom-left panel: minimum mass as a function of the mass of the host star. Bottom right panel: minimum mass as a function of the metallicity of the host star. The red dots represent the M-dwarf hosted RV extrasolar planets, the blue dots represent the HADES planets, and the blue star represents the newly detected planet GJ 685 b.

thumbnail Fig. 12

Properties of RV exoplanets orbiting M dwarfs. Minimum mass as a function of the orbital period of the planet. The red dots represent planets orbiting M dwarfs in multiple systems, while the blue dots represent single planets.

thumbnail Fig. 13

Properties of RV exoplanets orbiting M dwarfs. Orbital eccentricity as a function of the orbital period of the planet. The red dots represent planets orbiting M dwarfs in multiple systems, while the blue dots represent single planets. The marker size represents the minimum mass of the planet.

Acknowledgements

GAPS acknowledges support from INAF through the Progetti Premiali funding scheme of the Italian Ministry of Education, University, and Research. M.Pi. gratefully acknowledges the support from the European Union Seventh Framework Programme (FP7/2007-2013) under Grant Agreement No. 313014 (ETAEARTH). L.A., A.F.L., G.M., J.M., and M.D. acknowledge financial support from Progetto Premiale 2015 FRONTIERA funding scheme of the Italian Ministry of Education, University, and Research. G.S. acknowledges financial support from “Accordo ASI-INAF” No. 2013-016-R.0 July 9, 2013 and July 9, 2015. M.Pe. and I.R. acknowledge support from the Spanish Ministry of Economy and Competitiveness (MINECO) and the Fondo Europeo de Desarrollo Regional (FEDER) through grant ESP2016-80435-C2-1-R, as well as the support of the Generalitat de Catalunya/CERCA program. J.I.G.H., R.R. and B.T.P. acknowledge financial support from the Spanish Ministry project MINECO AYA2017-86389-P, and J.I.G.H. from the Spanish MINECO under the 2013 Ramón y Cajal programme MINECO RYC-2013-14875. This work is based on observations made with the Italian Telescopio Nazionale Galileo (TNG), operated on the island of La Palma by the Fundación Galielo Galilei of the Istituto Nazionale di Astrofisica (INAF) at the Spanish Observatorio del Roque de los Muchachos (ORM) of the Instituto de Astrofísica de Canarias (IAC). This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

Appendix A Observation log for GJ 685

In this section we report the observational data collected with the HARPS-N spectrograph as part the HADES project and used in the present study. We list in Table A.1 the observation dates (barycentric Julian date or BJD), the RVs calculated by the TERRA pipeline (Anglada-Escudé & Butler 2012), and asymmetry indexes BIS, ΔV, and Vasy derived as described in Lanza et al. (2018). The RV errors reported are the formal ones, not including the jitter term. In Table A.2 we list the observation dates and the activity indexes Ca II H and K, Hα, Na I D1 D2, and He I D3, derived following the procedure by Gomes da Silva et al. (2011).

Table A.1

RV and asymmetry indexes data of the 106 observed HARPS-N spectra of GJ 685.

Table A.2

Activity indexes data of the 106 observed HARPS-N spectra of GJ 685.

References

  1. Affer, L., Micela, G., Damasso, M., et al. 2016, A&A, 593, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Affer, L., Damasso, M., Micela, G., et al. 2019, A&A, 622, A193 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Ambikasaran, S., Foreman-Mackey, D., Greengard, L., Hogg, D. W., & O’Neil, M. 2015, IEEE Trans. Pattern Anal. Mach. Intell., 38, 2 [Google Scholar]
  4. Anglada-Escudé, G., & Butler, R. P. 2012, ApJS, 200, 15 [NASA ADS] [CrossRef] [Google Scholar]
  5. Anglada-Escudé, G., Amado, P. J., Barnes, J., et al. 2016a, Nature, 536, 437 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  6. Anglada-Escudé, G., Tuomi, M., Arriagada, P., et al. 2016b, ApJ, 830, 74 [NASA ADS] [CrossRef] [Google Scholar]
  7. Astudillo-Defru, N., Forveille, T., Bonfils, X., et al. 2017a, A&A, 602, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Astudillo-Defru, N., Díaz, R. F., Bonfils, X., et al. 2017b, A&A, 605, L11 [Google Scholar]
  9. Astudillo-Defru, N., Delfosse, X., Bonfils, X., et al. 2017c, A&A, 600, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Baluev, R. V. 2013a, MNRAS, 429, 2052 [NASA ADS] [CrossRef] [Google Scholar]
  11. Baluev, R. V. 2013b, Astron. Comput., 3, 50 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bonfils, X., Mayor, M., Delfosse, X., et al. 2007, A&A, 474, 293 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Bonfils, X., Delfosse, X., Udry, S., et al. 2013, A&A, 549, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Cosentino, R., Lovis, C., Pepe, F., et al. 2012, in Ground-based and Airborne Instrumentation for Astronomy IV, Proc. SPIE, 8446, 84461V [Google Scholar]
  15. Courcol, B., Bouchy, F., & Deleuil, M. 2016, MNRAS, 461, 1841 [NASA ADS] [CrossRef] [Google Scholar]
  16. Covino, E., Esposito, M., Barbieri, M., et al. 2013, A&A, 554, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  18. Desidera, S., Sozzetti, A., Bonomo, A. S., et al. 2013, A&A, 554, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Dressing, C. D., & Charbonneau, D. 2013, ApJ, 767, 95 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dumusque, X., Lovis, C., Udry, S., & Santos, N. C. 2011, in The Astrophysics of Planetary Systems: Formation, Structure, and Dynamical Evolution, eds. A. Sozzetti, M. G. Lattanzi, & A. P. Boss, IAU Symp., 276, 530 [NASA ADS] [Google Scholar]
  21. Dumusque, X., Boisse, I., & Santos, N. C. 2014, ApJ, 796, 132 [NASA ADS] [CrossRef] [Google Scholar]
  22. Dumusque, X., Borsa, F., Damasso, M., et al. 2017, A&A, 598, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Eastman, J., Gaudi, B. S., & Agol, E. 2013, PASP, 125, 83 [NASA ADS] [CrossRef] [Google Scholar]
  24. Endl, M., Kürster, M., Els, S., Hatzes, A. P., & Cochran, W. D. 2001, A&A, 374, 675 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Figueira, P., Santos, N. C., Pepe, F., Lovis, C., & Nardetto, N. 2013, A&A, 557, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
  27. Forveille, T., Bonfils, X., Delfosse, X., et al. 2009, A&A, 493, 645 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Gaidos, E., & Mann, A. W. 2014, ApJ, 791, 54 [NASA ADS] [CrossRef] [Google Scholar]
  30. Gaidos, E., Mann, A. W., Kraus, A. L., & Ireland, M. 2016, MNRAS, 457, 2877 [NASA ADS] [CrossRef] [Google Scholar]
  31. Giacobbe, P., Damasso, M., Sozzetti, A., et al. 2012, MNRAS, 424, 3101 [NASA ADS] [CrossRef] [Google Scholar]
  32. Giampapa, M. S., Cram, L. E., & Wild, W. J. 1989, ApJ, 345, 536 [NASA ADS] [CrossRef] [Google Scholar]
  33. Gomes da Silva, J., Santos, N. C., Bonfils, X., et al. 2011, A&A, 534, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. González-Álvarez, E., Micela, G., Maldonado, J., et al. 2019, A&A, 624, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Grunblatt, S. K., Howard, A. W., & Haywood, R. D. 2015, ApJ, 808, 127 [NASA ADS] [CrossRef] [Google Scholar]
  36. Haywood, R. D., Collier Cameron, A., Queloz, D., et al. 2014, MNRAS, 443, 2517 [NASA ADS] [CrossRef] [Google Scholar]
  37. Kipping, D. M. 2013, MNRAS, 434, L51 [NASA ADS] [CrossRef] [Google Scholar]
  38. Kopparapu, R. K., Ramirez, R., Kasting, J. F., et al. 2013, ApJ, 765, 131 [NASA ADS] [CrossRef] [Google Scholar]
  39. Kürster, M., Endl, M., Rouesnel, F., et al. 2003, A&A, 403, 1077 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Lanza, A. F., Malavolta, L., Benatti, S., et al. 2018, A&A, 616, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Laughlin, G., Bodenheimer, P., & Adams, F. C. 2004, ApJ, 612, L73 [NASA ADS] [CrossRef] [Google Scholar]
  42. Limbach, M. A., & Turner, E. L. 2015, Proc. Natl. Acad. Sci., 112, 20 [NASA ADS] [CrossRef] [Google Scholar]
  43. Lovis, C., & Pepe, F. 2007, A&A, 468, 1115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Luque, R., Nowak, G., Pallé, E., et al. 2018, A&A, 620, A171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Maldonado, J., Affer, L., Micela, G., et al. 2015, A&A, 577, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Maldonado, J., Scandariato, G., Stelzer, B., et al. 2017, A&A, 598, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
  48. Mordasini, C., Alibert, Y., & Benz, W. 2009, A&A, 501, 1139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Mordasini, C., Alibert, Y., Benz, W., Klahr, H., & Henning, T. 2012, A&A, 541, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Mortier, A., Santos, N. C., Sozzetti, A., et al. 2012, A&A, 543, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Mortier, A., Santos, N. C., Sousa, S. G., et al. 2013, A&A, 557, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Mortier, A., Faria, J. P., Correia, C. M., Santerne, A., & Santos, N. C. 2015, A&A, 573, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Moutou, C., Vigan, A., Mesa, D., et al. 2017, A&A, 602, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Perger, M., Ribas, I., Damasso, M., et al. 2017a, A&A, 608, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Perger, M., García-Piquer, A., Ribas, I., et al. 2017b, A&A, 598, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Perger, M., Scandariato, G., Ribas, I., et al. 2019, A&A, 624, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Pinamonti, M., Sozzetti, A., Bonomo, A. S., & Damasso, M. 2017, MNRAS, 468, 3775 [NASA ADS] [CrossRef] [Google Scholar]
  58. Pinamonti, M., Damasso, M., Marzari, F., et al. 2018, A&A, 617, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Poretti, E., Boccato, C., Claudi, R., et al. 2016, Mem. Soc. Astron. It., 87, 141 [NASA ADS] [Google Scholar]
  60. Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Quirrenbach, A., Amado, P. J., Caballero, J. A., et al. 2014, in Ground-based and Airborne Instrumentation for Astronomy V, Proc. SPIE, 9147, 91471F [Google Scholar]
  62. Reiners, A., Zechmeister, M., Caballero, J. A., et al. 2018, A&A, 612, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Robertson, P., Endl, M., Cochran, W. D., & Dodson-Robinson, S. E. 2013, ApJ, 764, 3 [NASA ADS] [CrossRef] [Google Scholar]
  64. Robertson, P., Mahadevan, S., Endl, M., & Roy, A. 2014, Science, 345, 440 [NASA ADS] [CrossRef] [Google Scholar]
  65. Robertson, P., Roy, A., & Mahadevan, S. 2015, ApJ, 805, L22 [NASA ADS] [CrossRef] [Google Scholar]
  66. Rodigas, T. J., & Hinz, P. M. 2009, ApJ, 702, 716 [NASA ADS] [CrossRef] [Google Scholar]
  67. Santos, N. C., Israelian, G., Mayor, M., et al. 2005, A&A, 437, 1127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Scandariato, G., Maldonado, J., Affer, L., et al. 2017, A&A, 598, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Schwarz, G. 1978, Ann. Stat., 6, 461 [Google Scholar]
  70. Sousa, S. G., Santos, N. C., Mayor, M., et al. 2008, A&A, 487, 373 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
  71. Sozzetti, A., Bernagozzi, A., Bertolini, E., et al. 2013, Eur. Phys. J. Web Conf., 47, 03006 [CrossRef] [Google Scholar]
  72. Stepinski, T. F., & Black, D. C. 2000, A&A, 356, 903 [NASA ADS] [Google Scholar]
  73. Suárez Mascareño, A., González Hernández, J. I., Rebolo, R., et al. 2017, A&A, 605, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Suárez Mascareño, A., Rebolo, R., González Hernández, J. I., et al. 2018, A&A, 612, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Tamuz, O., Ségransan, D., Udry, S., et al. 2008, A&A, 480, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Tuomi, M., Jones, H. R. A., Barnes, J. R., Anglada-Escudé, G., & Jenkins, J. S. 2014, MNRAS, 441, 1545 [NASA ADS] [CrossRef] [Google Scholar]
  77. Valenti, J. A., & Fischer, D. A. 2005, VizieR Online Data Catalog: J/ApJS/159/141 [Google Scholar]
  78. Wang, J., & Fischer, D. A. 2015, AJ, 149, 14 [NASA ADS] [CrossRef] [Google Scholar]
  79. Wittenmyer, R. A., Wang, S., Horner, J., et al. 2013, ApJS, 208, 2 [NASA ADS] [CrossRef] [Google Scholar]
  80. Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Zucker, S., & Mazeh, T. 2002, ApJ, 568, L113 [NASA ADS] [CrossRef] [Google Scholar]

2

Suárez Mascareño et al. (2018) extended the definition of logRHK$\log R&#x0027;_{\text{HK}}$ for application on M-dwarfs spectra, following a procedure very similar to the one used by Astudillo-Defru et al. (2017c).

3

For the comparative analysis of the different periodogram algorithms we computed only the theoretical FAPs, since the bootstrap simulation, in particular applied to FREDEC, is very time-consuming.

4

In the analysis of the residuals of the He I D3 time series (not shown), only an additional signal at P = 9 d emerges, corresponding to the rotation period first harmonic.

6

The minimum mass distribution of detected exoplanets could be affected by observation bias, hiding the presence of additional undetected companions.

All Tables

Table 1

Stellar parameters for the target GJ 685.

Table 2

Priors and best-fit results for the Gaussian process regression analysis of the Hα activity indicator.

Table 3

Priors and best-fit results for the emcee analysis of GJ 685 RV time series.

Table A.1

RV and asymmetry indexes data of the 106 observed HARPS-N spectra of GJ 685.

Table A.2

Activity indexes data of the 106 observed HARPS-N spectra of GJ 685.

All Figures

thumbnail Fig. 1

HARPS-N RV time series of GJ 685.

In the text
thumbnail Fig. 2

Upperpanel: GLS periodograms of the APACHE photometric data. The dotted and dashed horizontal lines indicate the 1 and 10% FAP levels, respectively. Lower panel: APACHE light curve phase-folded over the 16.85 d period found by the GLS periodogram. The red line represents the best-fit model.

In the text
thumbnail Fig. 3

GLS periodograms of the RV time series and residuals after consecutive signal fits. The red vertical lines indicate the potential rotation period (dashed) and its first harmonic (dot-dashed), while the dashed blue vertical lines marks the orbital period of the planetary candidate. The dotted, dash-dotted, and dashed horizontal lines indicate the 0.1, 1, and 10% FAP levels, respectively.

In the text
thumbnail Fig. 4

GLS periodograms of the asymmetry indicators. The red vertical lines indicate the potential rotation period (dashed) and its first harmonic (dot-dashed), while the dashed blue vertical lines mark the orbital period of the planetary candidate discussed in Sect. 5.1. The dotted horizontal lines indicate the 1% FAP level.

In the text
thumbnail Fig. 5

GLS periodograms of the activity indexes. The red vertical lines indicate the potential rotation period (dashed) and its first harmonic (dot-dashed), while the dashed blue vertical lines mark the orbital period of the planetary candidate discussed in Sect. 5.1. The dotted horizontal lines indicate the 1% FAP level.

In the text
thumbnail Fig. 6

Posterior distributions of the fitted (hyper-)parameters of the GP quasi-periodic model applied to the time series of Hα activity index. The vertical dashed lines denote the median and the 16th and 84th percentiles.

In the text
thumbnail Fig. 7

Posterior distributions of the fitted (hyper-)parameters of the pure GP model applied to the RV time series. The vertical dashed lines denote the median and the 16 th and 84 th percentiles.

In the text
thumbnail Fig. 8

Posterior distributions of the fitted (hyper-)parameters of the GP + 1 planet model applied to the RV time series. The verticaldashed lines denote the median and the 16th and 84th percentiles.

In the text
thumbnail Fig. 9

Upper panel: best-fit stellar quasi-periodic signal obtained from the GP + 1 planet model (orange line) compared to the RV residuals (blue points). Lower panel: magnification of the last HARPS-N observing season. The grey area represents the ± 1σ uncertainties of the stellar activity model.

In the text
thumbnail Fig. 10

Phase-folded for the RV signal of GJ 685 b, after the subtraction of the stellar correlated signal.

In the text
thumbnail Fig. 11

Properties of RV exoplanets orbiting M dwarfs. Top-left panel: orbital period as a function of the mass of the host star. Top-right panel: orbital period as a function of the metallicity of the host star. Bottom-left panel: minimum mass as a function of the mass of the host star. Bottom right panel: minimum mass as a function of the metallicity of the host star. The red dots represent the M-dwarf hosted RV extrasolar planets, the blue dots represent the HADES planets, and the blue star represents the newly detected planet GJ 685 b.

In the text
thumbnail Fig. 12

Properties of RV exoplanets orbiting M dwarfs. Minimum mass as a function of the orbital period of the planet. The red dots represent planets orbiting M dwarfs in multiple systems, while the blue dots represent single planets.

In the text
thumbnail Fig. 13

Properties of RV exoplanets orbiting M dwarfs. Orbital eccentricity as a function of the orbital period of the planet. The red dots represent planets orbiting M dwarfs in multiple systems, while the blue dots represent single planets. The marker size represents the minimum mass of the planet.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.