# Search for giant planets in M 67

## IV. Survey results★, ★★

A. Brucalassi<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, J. Koppenhoefer<fnref n="1" />,<fnref n="2" />, R. Saglia<fnref n="1" />,<fnref n="2" />, L. Pasquini<fnref n="3" />, M. T. Ruiz<fnref n="4" />, P. Bonifacio<fnref n="5" />, L. R. Bedin<fnref n="6" />, M. Libralato<fnref n="6" />,<fnref n="7" />, K. Biazzo<fnref n="8" />, C. Melo<fnref n="9" />, C. Lovis<fnref n="10" />, and S. Randich<fnref n="11" />

Received 15 November 2015 / Accepted 26 February 2017

## ABSTRACT

**Context.** We present the results of a seven-year-long radial velocity survey of a sample of 88 main-sequence and evolved stars to reveal signatures of Jupiter-mass planets in the solar-age and solar-metallicity open cluster M 67.

**Aims.** We aim at studying the frequency of giant planets in this cluster with respect to the field stars. In addition, our sample is also ideal to perform a long-term study to compare the chemical composition of stars with and without giant planets in detail.

**Methods.** We analyzed precise radial velocity (RV) measurements obtained with the HARPS spectrograph at the European Southern Observatory (La Silla), the SOPHIE spectrograph at the Observatoire de Haute-Provence (France), the HRS spectrograph at the Hobby Eberly Telescope (Texas), and the HARPS-N spectrograph at the Telescopio Nazionale Galileo (La Palma). Additional RV data come from the CORALIE spectrograph at the Euler Swiss Telescope (La Silla). We conducted Monte Carlo simulations to estimate the occurrence rate of giant planets in our radial velocity survey. We considered orbital periods between 1.0 day and 1000 days and planet masses between 0.2 *M*<sub>1</sub> and 10.0 *M*<sub>1</sub>. We used a measure of the observational detection efficiency to determine the frequency of planets for each star.

**Results.** All the planets previously announced in this RV campaign with their properties are summarized here: 3 hot Jupiters around the main-sequence stars YBP1194, YBP1514, and YBP401, and 1 giant planet around the evolved star S364. Two additional planet candidates around the stars YBP778 and S978 are also analyzed in the present work. We discuss stars that exhibit large RV variability or trends individually. For 2 additional stars, long-term trends are compatible with new binary candidates or substellar objects, which increases the total number of binary candidates detected in our campaign to 14. Based on the Doppler-detected planets discovered in this survey, we find an occurrence of giant planets of ~18.0±8.0% in the selected period-mass range. This frequency is slightly higher but consistent within the errors with the estimate for the field stars, which leads to the general conclusion that open cluster and field statistics agree. However, we find that the rate of hot Jupiters in the cluster (~5.7±3.0%) is substantially higher than in the field.

**Key words.** techniques: radial velocities – planets and satellites: gaseous planets

## 1. Introduction

In recent years, several observational campaigns have been dedicated to search for planets in clusters or stellar associations, where the majority of stars is considered to form. Stars in

clusters constitute a homogeneous sample in age and chemical composition that is ideal for investigating the dependence of planet formation on the mass and the properties of the central star (González Hernández et al. 2013, 2010; Baumann et al. 2010; Johnson et al. 2010; Ramírez et al. 2010; Meléndez et al. 2009), for determining the formation timescale and distinguishing different migration processes (Dong et al. 2014; Quinn et al. 2014; Dawson & Murray-Clay 2013), and finally for modeling the effects of stellar encounters on the formation and evolution of planetary systems (Li & Adams 2015; Cai et al. 2016; Shara et al. 2016; Davies et al. 2014; Malmberg et al. 2011; Spurzem et al. 2009). However, the high-precision radial velocity (RV) technique and the transit method have been successful

<footnote n="1">¹ Max-Planck für extraterrestrische Physik, 85741 Garching bei München, Germany</footnote>

<footnote n="2">² University Observatory Munich, Ludwig Maximillian Universitaet, Scheinerstrasse 1, 81679 Munich, Germany</footnote>

<footnote n="3">³ ESO-European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748 Garching bei München, Germany
e-mail: abrucala@eso.org</footnote>

<footnote n="4">⁴ Astronomy Department, Universidad de Chile, 1058 Santiago, Chile</footnote>

<footnote n="5">⁵ GEPI, Observatoire de Paris, PSL Research University, CNRS, place Jules Janssen, 92190 Meudon, France</footnote>

<footnote n="6">⁶ Istituto Nazionale di Astrofisica, Osservatorio Astronomico di Padova, 35122 Padova, Italy</footnote>

<footnote n="7">⁷ Università degli Studi di Padova, UNIPD, 35122 Padova, Italy</footnote>

<footnote n="8">⁸ Istituto Nazionale di Astrofisica, Osservatorio Astrofisico di Catania, 95125 Catania, Italy</footnote>

<footnote n="9">⁹ ESO-European Southern Observatory, 3107 Santiago, Chile</footnote>

<footnote n="10">¹⁰ Observatoire de Genève, Sauverny, 1290 Versoix, Switzerland</footnote>

<footnote n="11">¹¹ Istituto Nazionale di Astrofisica, Osservatorio Astrofisico di Arcetri, 50125 Firenze, Italy</footnote>

<footnote>* Based on observations collected at the ESO 3.6m telescope (La Silla), at the 1.93 m telescope of the Observatoire de Haute-Provence (OHP, France), at the Hobby Eberly Telescope (HET, Texas), at the Telescopio Nazionale Galileo (TNG, La Palma) and at the Euler Swiss Telescope (La Silla).</footnote>

<footnote>** Individual RV measurements are available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/603/A85</footnote>
only very recently in discovering planetary-mass companions around stars belonging to open clusters (OC). Of the most recent results in this field, we mention the detection of a hot Jupiter and a multi-planetary system in the Praesepe OC (Quinn et al. 2012; Malavolta et al. 2016), of a hot Jupiter in the Hyades (Quinn et al. 2014), of two sub-Neptune planets in NGC 6811 (Meibom et al. 2013), and the discovery of five Jupiter-mass planets in M67 (Brucalassi et al. 2014, 2016). Previous RV surveys provided detections of a long-period giant planet around one of the Hyades clump giants (Sato et al. 2007) and a substellar-mass object in NGC 2423 (Lovis & Mayor 2007). No significant exoplanet candidates have been found by either ground- or by space-based transit campaigns in globular clusters (Gilliland et al. 2000; Nascimbeni et al. 2012).

These discoveries confirm that giant planets exist in a dense cluster environment, and suggest that a complete census of planet discoveries may be biased by the detection limit of the instruments and the observations available today. Moreover, results from various simulations show that dense birth environments such as stellar clusters can significantly influence the planet formation process and the resulting orbital properties of the planetary systems. Close encounters between stars or binary companions can modify the structure of any planetary system and also subsequently generate strong interactions between planets over very long timescales (Davies et al. 2014; Malmberg et al. 2011). This leads to the ejection of some planets, but it also seems to favor the conditions for the formation of hot Jupiters (Shara et al. 2016). Studying hot Jupiters in OCs can therefore shed light onto the long-standing problem of identifying their dynamical origin.

For the past seven years we have carried out a search for massive planets around main-sequence (MS) and evolved stars in the OC M67. The scientific motivations for these studies include determining the impact of a different environment on the frequency and the evolution of planetary systems with respect to field stars. As a long-term goal, we aim to study the connection between giant planet formation and stellar mass and chemical composition. M67 is the perfect target to search for planets around OC stars. Chemical analysis from several works (Randich et al. 2006; Pace et al. 2008; Önehag et al. 2011, 2014) has shown that M67 has a chemical composition (not only Fe, but also the other elements) that is extremely similar to solar, as close as allowed by the precision of the measurements. In addition, the age of M67 (3–5 Gyr), according to numerous determinations, encompasses the accepted value of the Sun, while the age determination for field stars is always rather uncertain. Finally, M67 is a rich OC, which gives us the opportunity to find many stellar candidates that share similar properties and a large number of stars with different masses, characteristics that are essential to address the questions above.

Through RV measurements obtained with HARPS at La Silla-ESO, SOPHIE at OHP, HRS at HET, HARPS-N at the TNG, and CORALIE at the *Euler* Swiss Telescope, five new giant planets have been discovered around M67 stars (Brucalassi et al. 2014, 2016): three hot Jupiters around MS stars, and two long-period planets around evolved stars.

In this paper, we present our RV campaign around stars of the OC M67 considering data obtained until March 2015, and we provide a complete census of all the stars.

The star sample is described in Sect. 2, RV observations and analyses are reported in Sect. 3. In Sect. 4 we provide more detailed information on individual objects. In Sect. 5 we report a series of simulations based on a Monte Carlo approach that we

<table><caption>**Table 1.** List of the observations.</caption><thead><tr><th>Instrument</th><th>H</th><th>S</th><th>HET</th><th>C</th><th>HN</th></tr></thead><tbody><tr><td>N. Stars</td><td>88</td><td>70</td><td>24</td><td>14</td><td>13</td></tr><tr><td>Observations</td><td>734</td><td>168</td><td>125</td><td>99</td><td>23</td></tr><tr><td>MS</td><td>481</td><td>75</td><td>59</td><td>0</td><td>12</td></tr><tr><td>TO</td><td>63</td><td>22</td><td>23</td><td>0</td><td>7</td></tr><tr><td>G</td><td>190</td><td>78</td><td>43</td><td>99</td><td>4</td></tr></tbody></table>

**Notes.** The table lists: number of observed stars, total number of observations, number of main-sequence (MS), turn-off (TO), and giant (G) stars observed for each instrument. HARPS (H), SOPHIE (S), HET, CORALIE (C), and HARPS-N (HN).

used to estimate the occurrence rate of giant planets in our radial velocity survey. Finally, we summarize our results in Sect. 6.

## 2. Sample and observations

A complete description of the sample is reported in our previous works (Pasquini et al. 2012; Brucalassi et al. 2014). We highlight here those characteristics that are most relevant for the aim of this paper.

The original M67 sample includes a total of 88 stars with V mag between 9 and 15, and a mass range of 0.9–1.4 $M_\odot$.

Main-sequence stars have been selected following Pasquini et al. (2008). In particular, we considered those stars with a membership probability higher than 60% and a proper motion lower than 6 mas/yr with respect to the average according to Yadav et al. (2008). When we considered the selection of the giants, we referred to Sanders (1977) for the membership probability, and the RV membership was derived according to Mermilliod & Mayor (2007) and Mathieu et al. (1986). We are aware that several stars that are particularly close to the turn-off point, although fulfilling our selection criteria, have been not observed because the observation time was too short.

Five different telescopes and instrument combinations have been used to obtain the RV measurements. Table 1 summarizes the number of the observations for each instrument.

The HARPS spectrograph (Mayor et al. 2003) at the ESO 3.6 m telescope was used in high-efficiency mode (EGGS mode) since the M67 stars are quite faint for this instrument. In this configuration the fiber has a aperture on the sky of 1.2 arcsec, corresponding to $R = 90\,000$, and is not equipped with an optical scrambler. The precision is limited to $5-7\ m\text{s}^{-1}$, but the efficiency is 30–40% better than the high-resolution mode (HAM mode), which allows us to obtain the highest resolving power of 115 000 with a fiber aperture of 1 arcsec. The spectral range covered is 378–691 nm. Exposure times ranged from 10 to 45 min, yielding a typical signal-to-noise ratio (S/N) per resolution element of 10 to 15 for the faintest stars. Between January 2008 and March 2015, we gathered 734 observations of 88 stars with HARPS with an associated internal precision of ~10 $m\text{s}^{-1}$. These observations represent the majority of this survey, and we therefore consider HARPS as our reference.

The SOPHIE spectrograph (Bouchy & Sophie Team 2006) at the OHP 1.93 m telescope was used in high-efficiency mode with $R = 40\,000$ and a spectral range of 387–694 nm. We analyzed 168 SOPHIE observations of 70 M67 stars with an associated internal precision of ~12 $m\text{s}^{-1}$.

For the HRS spectrograph (Tull et al. 1998) at the Hobby Eberly Telescope we opted for a configuration with $R = 60\,000$.
<figure><img src="image_1.png" /><figcaption>**Fig. 1.** Histogram showing the number of observations per star for our total sample. All observations from HARPS, HARPS-N, SOPHIE, CORALIE, and HET HRS are included in the plot and updated to March 31, 2015. The order number of the stars corresponds to Table A.1.</figcaption></figure>

and a wavelength range of 407.6–787.5 nm. We obtained 125 observations for 24 stars with a typical error bar associated with the observations of ~25 m s⁻<fnref n="1" />. During the period 2013–2015, the HET telescope was not accessible because of a telescope upgrade.

Moreover, we gathered 99 RV data points for 14 giant stars observed between 2003 and 2005 (Lovis & Mayor 2007) with the CORALIE spectrograph at the 1.2 m Euler Swiss telescope.

In addition, 27 h in service mode have been allocated with the TNG on La Palma Canary Island and the HARPS-N spectrograph. The HARPS-N is a fiber-fed echelle spectrograph, similar to HARPS on the 3.6 m ESO telescope and covers the wavelength range between 383 to 693 nm, with a spectral resolution R = 115 000. The two HARPS fibers (object + Sky/ThAr) have an aperture on the sky of 1 arcsec and are equipped with an age scrambler to provide a uniform spectrograph pupil illumination. Unfortunately, the observing campaign during the winter 2013–2014 was not very successful. The run at the TNG was impacted by bad weather and only 10 observing hours could be used, with 23 spectra for 13 stars.

Figure 1 shows a histogram with the number of observations per star updated to March 31, 2015. We have obtained 13 observations per star on average: 10 observations per star for MS stars, and 18 observations per star for giants and turn-off stars. Table A.1 shows the main data for the observed stars. In addition to the basic stellar parameters, the number of observations per star is given for each spectrograph and as a total.

## 3. RV analysis

HARPS, SOPHIE, and HARPS-N are provided with a similar automatic pipeline. The spectra are extracted from the detector images and cross-correlated with a numerical mask. For all of our stars, irrespective of the spectral type and luminosity, we used a G2V mask obtained from Sun spectra. Radial velocities are derived by fitting each resulting cross-correlation function (CCF) with a Gaussian (Baranne et al. 1996; Pepe et al. 2002). This real-time pipeline also provides an associated internal RV error $\sigma_{\text{pn}}$ (photon noise error).

For the HRS data, the radial velocities were computed using a series of dedicated routines based on IRAF<fnref n="1" />. The procedures are described in more detail in Cappetta et al. (2012). The different steps include the wavelength calibration using a Th-Ar lamp exposure acquired before and after each stellar spectrum, the normalization of the spectra, the cleaning of cosmic rays as well as telluric and sky lines, the cross-correlation of the spectrum performed order by order with a G2V star template, and finally, the computation of the heliocentric corrections. The resulting CCFs were fitted with a Gaussian function to determine the RVs. After the orders affected by telluric lines and low S/N are rejected, the final RV value is given by the average value over the retained orders. Finally, the internal RV uncertainties are calculated by $\sigma_{\text{pn}} = \text{rms}(v)/\sqrt{N}$, where v is the RV of the individual orders and N is the number of the orders.

We used nightly observations of the RV standard star HD32923 to correct all observations for each star to the zero-point of HARPS (as explained in Pasquini et al. 2012) and to take any instrument instability or any systematic velocity shifts between runs (such as the modification of the SOPHIE fiber link in June 2011 (Perruchot et al. 2011) or technical problems during the calibration phase at HARPS) into account. An additional correction was applied to the SOPHIE data to account for the low S/N of the observations. For this, we corrected our spectra using Eq. (1) in Santerne et al. (2012). Finally, we considered the total error of the RV measurements ($\sigma_{\text{obs}}$) to be the sum in quadrature of all RV error sources described above. After all the observations for each star were corrected to the zero-point of HARPS, they were analyzed together. The combined measurement uncertainties have been compared to the observed velocity dispersions to evaluate the significance of any potential velocity variation and to highlight any possible outliers that would suggest the presence of exoplanets. In Table A.1 the mean stellar RV of each star is given, together with the average associated error ($\sigma_{\text{obs}}$) and the RV dispersion of the observations ($\sigma_{\text{RV}}$). Figure 2 shows (top panel) the histogram of the observed RV scatter ($\sigma_{\text{RV}}$) for the sample stars. Binary candidates were excluded, and the stars with significant RV trend are not shown. The bulk of our observations have an RV scatter represented by a Gaussian distribution centered on ~21 ± 7 m s⁻<fnref n="1" />. Stars with an RV variability at or above 40 m s⁻<fnref n="1" /> are considered very good candidates for low-mass companion hosts.

Figure 2 also includes a second panel (bottom) showing the histogram of the ratio between the observed RV scatter ($\sigma_{\text{RV}}$) and the average associated error ($\sigma_{\text{obs}}$) for each star. Large variations for planet candidates or long-term objects (see Sect. 4) are easy to identify. A clear peak results for most of the stars. Considering only stars with $\sigma_{\text{RV}} < 4\sigma_{\text{obs}}$ (the small variations) and with no suspected companions, we calculated a peak of ~1.3 (employing a kernel density estimator to smooth the distribution), which is a sign that we may have underestimated the measurement errors. However, halos of excess scatter (RV variability) are present, which are most likely due to a combination of effects. In evolved stars some measurable intrinsic RV variability is present, as shown also in Setiawan et al. (2004) and Hekker et al. (2008), while for the faint MS stars the uncertainty in the measurements increases because of the limited S/N. Finally, some of the stars still have only a few observational points with poorly constrained RV variability, and their scatter is very likely the result of our low data statistics. Pasquini et al. (2012) also investigated whether other instrumental effects such as the observed flux or airmass (not included in the data analysis) could

<footnote n="1">¹ http://iraf.noao.edu/</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 2.** Top: histogram showing the observed velocity dispersions ($\sigma_{RV}$) for the MS stars (pink) and giants and turn-off stars (blue) of our sample. Binary candidates and stars with high RV trend have been excluded. The black line represents the Gaussian fit centered at ~21 m s⁻¹ of width $\sigma \sim 7$ m s⁻¹. Bottom: histogram showing the ratio of the RV dispersion ($\sigma_{RV}$) and the average measurements error ($\sigma_{obs}$). The symbols are the same as in the top panel.</figcaption></figure>

affect the RV measurement precision at low count levels, but no correlation was found. In order to consider these effects during the computation of the giant planet occurrence rate (see Sect. 5), we empirically inflated the estimated errors by adding in quadrature a term ($\sigma'_{RV}$) equivalent, on average, to the excess RV scatter. This is calculated by $\sigma'_{RV} = \sqrt{\sigma^2_{RV,st} - \sigma^2_{obs,st}}$, where $\sigma^2_{RV,st}$ and $\sigma^2_{obs,st}$ are the mean RV scatter and the mean estimated error, respectively, of stars without a trend and suspected companions. In Table A.1 the corrected associated error for each star is reported as $\sigma_{cor}$.

We studied the RV variations of our target stars by computing the Lomb-Scargle periodogram (Scargle 1982; Horne & Baliunas 1986). This is a commonly used technique for searching for periodic sinusoidal signals in unevenly sampled data and allows estimates of the detection threshold to be written down for periods shorter than the duration of the observations (Horne & Baliunas 1986). The significance of the sinusoid best fit of our RV values was determined by analytically calculating the false-alarm probability (FAP) level (Horne & Baliunas 1986). Afterward, we applied Levenberg-Marquardt analysis (Wright & Howard 2009, RVLIN) to fit Keplerian orbits to the radial velocity data.

<table><caption>**Table 2.** RV measurements, RV uncertainties, and instrument.</caption><thead><tr><th>Star</th><th>BJD<br>(-2 450 000)</th><th>RV<br>(km s<sup>-1</sup>)</th><th>σ<sub>obs</sub><br>(km s<sup>-1</sup>)</th><th>Instrument</th></tr></thead><tbody><tr><td rowspan="5">YBP266</td><td>4488.509042</td><td>33.78637</td><td>0.012</td><td>SOPHIE</td></tr><tr><td>4855.591058</td><td>33.77123</td><td>0.012</td><td>HARPS</td></tr><tr><td>4859.605234</td><td>33.77318</td><td>0.014</td><td>HARPS</td></tr><tr><td>4862.702784</td><td>33.77678</td><td>0.017</td><td>HARPS</td></tr><tr><td>5189.693914</td><td>33.79851</td><td>0.017</td><td>HARPS</td></tr></tbody></table>

**Notes.** All RV data points are corrected to the zero-point of HARPS. The errors $\sigma_{obs}$ do not include the correction for the excess RV scatter (see Sect. 3).

Radial velocity periodic variation can be caused by rotational inhomogeneities related to stellar surface activity, such as plages or spots, including the one that is due to magnetic cycles of several years (Santos et al. 2010). Stellar activity can be diagnosed with spectral indicators or by monitoring the shape of the spectral lines. The low S/N of our observations does not provide sufficient signal in the region of the sensitive Ca II H and K lines. We therefore followed a method similar to the one described in Pasquini & Pallavicini (1991). We investigated the presence and variability of chromospheric active regions in these stars by measuring the variations of the core of $H\alpha$ with respect to the continuum (see Dollinger 2008, for a more detailed description of how the $H\alpha$ was measured). For each case we verified the correlation between the RVs and the bisector span of the CCF (calculated following Queloz et al. 2001) or with the FWHM of the CCF.

# 4. Results

In the following, we describe the general results of our survey and discuss our analyses and interpretation of the RV time series for the most relevant targets in more detail. In particular, two further planet candidates are reported in this work. For the other stars, we refer to Table A.1, in which the main data for the observed objects are summarized. Individual RV measurements for all the stars are provided at the CDS in the format listed in Table 2.

## 4.1. Binary candidates

Twelve binary candidates were previously published in Pasquini et al. (2012). Here we present two other MS stars in the sample (YBP1051, YBP673) that either show significant linear trends or RV variations that are too large to be produced by an exoplanet or by a non-stellar object. The RV measurements of these stars are shown in Fig. 3. They display a peak-to-peak RV amplitude of ~1.0 km s⁻¹ with an RV range spanned over more than six years. All binaries are highlighted in boldface in Table A.1.

## 4.2. Planets and planetary candidates

Four planets were announced in our previous works (Brucalassi et al. 2014, 2016) for the stars YBP401, YBP1194, YBP1514, and S364. Table 3 reports the main stellar characteristics and the resulting updated orbital parameters of the Keplerian fit for all these stars. Another two stars (S978 and YBP778) show significant indication for the presence of Jovian-mass companions.

YBP778. Twenty-one RV measurements have been obtained for the MS star YBP778 since 2009: 15 with HARPS, the others with SOPHIE and HARPS-N. The typical S/N is 10 and the
<figure><img src="image_3.png" /><figcaption>**Fig. 3.** Stars with large RV or significant linear trends (suspected binaries). Black dots: HARPS measurements, red dots: SOPHIE measurements, orange dots: HARPS-N measurements, and green dots: HRS measurements. RV error bars represent $\sigma_{cor}$.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 4.** Phased RV measurements and Keplerian best fit, best-fit residuals, and bisector variation for YBP778. Same symbols as in Fig. 3.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 5.** Phased RV measurements and Keplerian best fit, best-fit residuals, and bisector variation for S978. Same symbols as in Fig. 4. Green dots: HRS measurements, blue dots: CORALIE measurements.</figcaption></figure>

average measurement uncertainty is ~17 m s⁻¹ for HARPS and HARPS-N, and ~10 m s⁻¹ for SOPHIE. A clear periodic signal can be seen in the periodogram at ~398 days (see Fig. A.3). Thus a single-planet Keplerian model was adjusted to the data, and the best-fit solution corresponds to a signal with a period of 410.4 ± 6.2 days, a semi-amplitude 158.8 ± 21.4 m s⁻¹ (see Fig. 4), and an eccentricity of 0.27 ± 0.11. Although no peak is present in the periodograms of the activity index or the bisector span, nor in the CCF FWHM (see Fig. A.3), we detected an anticorrelation between the RVs and the bisector span with a Pearson correlation coefficient of -0.56, and a slight correlation between the radial velocities and the CCF FWHM (see Fig. A.4). The correlation disappears when we consider the RV residual.

We therefore cannot completely exclude stellar activity (magnetic cycles) or a binary companion with low sin(i) as the cause of the RV variation at this point, although in the recent study of Geller et al. (2015), YBP778 is also classified as a single member. Interesting, Pace et al. (2012) noted that this star is significantly underabundant in lithium compared to the lithium abundances as a function of stellar mass for MS stars in M 67. This could argue in favor of the possibility that the star might indeed host a planet (see, e.g., Deal et al. 2015). However, Pace et al. (2012) warned that the distance of this star to the isochrone in the color-magnitude diagram (CMD; ~0.06 mag)

could also suggest the presence of a stellar or substellar companion.

S978. The star is a K4 red giant and was observed 40 times over more than five years with HARPS, HRS, SOPHIE, and HARPS-N. The average RV uncertainty is ~3.0 m s⁻¹ for HARPS and SOPHIE, ~26.0 m s⁻¹ for HRS, and ~8.0 m s⁻¹ for HARPS-N. Ten additional RV measurements were obtained with CORALIE between 2003 and 2005, with a mean measurement uncertainty of ~12.0 m s⁻¹. The strongest peak in the periodogram of the RV time series lies at about 510.4 days (see Fig. A.3). Aliases are also present in the periodogram, but at half of the power of the main signal. The bisector inverse slope as well as the activity indicator do not present any significant variation at the period of 510.4 days. We fitted a single-planet Keplerian orbit to this signal (see Fig. 5) and found an orbital solution whose parameters are reported in Table 3. The dispersion of the residuals is $\sigma(O-C) = 12.90$ m s⁻¹ and the periodogram of the residuals reveals some structures, but no significant peaks (see Fig. A.3).

Several studies indicate that M 67 stars in general have a low level of chromospheric activity. Pace & Pasquini (2004) computed activity levels for M 67 and other clusters. The paper shows that, for instance, the young Hyades and Pre-saepe stars on average have a chromospheric emission flux
<table><caption>**Table 3.** Top: stellar parameters of the M 67 stars that are newly found to host planet candidates. Bottom: orbital parameters of the planetary companions.</caption><thead><tr><th>Parameters</th><th>YBP401</th><th>YBP1194</th><th>YBP1514</th><th>SAND364</th><th>SAND978</th></tr></thead><tbody><tr><td>α (J2000)</td><td>08:51:19.05</td><td>08:51:00.81</td><td>08:51:00.77</td><td>08:49:56.82</td><td>08:51:17.48</td></tr><tr><td>δ (J2000)</td><td>+11:40:15.80</td><td>+11:48:52.76</td><td>+11:53:11.51</td><td>+11:41:33.00</td><td>+11:45:22.69</td></tr><tr><td>Spec.type</td><td>F9V</td><td>G5V</td><td>G5V</td><td>K3III</td><td>K4III</td></tr><tr><td>m<sub>V</sub> [mag]</td><td>13.70<sup>a</sup></td><td>14.6<sup>a</sup></td><td>14.77<sup>a</sup></td><td>9.8<sup>b</sup></td><td>9.71<sup>c</sup></td></tr><tr><td>B − V [mag]</td><td>0.607<sup>a</sup></td><td>0.626<sup>a</sup></td><td>0.680<sup>a</sup></td><td>1.360<sup>b</sup></td><td>1.370<sup>c</sup></td></tr><tr><td>M★ [M<sub>☉</sub>]</td><td>1.14 ± 0.02<sup>d</sup></td><td>1.01 ± 0.02<sup>d</sup></td><td>0.96 ± 0.01<sup>d</sup></td><td>1.35 ± 0.05<sup>d</sup></td><td>1.37 ± 0.02<sup>d</sup></td></tr><tr><td>log g [cgs]</td><td>4.30 ± 0.035<sup>f</sup></td><td>4.44 ± 0.05<sup>e</sup></td><td>4.57 ± 0.06<sup>g</sup></td><td>2.20 ± 0.06<sup>h</sup></td><td>1.80 ± 0.09<sup>i</sup></td></tr><tr><td>T<sub>eff</sub> [K]</td><td>6165 ± 64<sup>f</sup></td><td>5780 ± 27<sup>e</sup></td><td>5725 ± 45<sup>g</sup></td><td>4284 ± 9<sup>h</sup></td><td>4200 ± 21<sup>i</sup></td></tr><tr><td>P [days]</td><td>4.087 ± 0.003</td><td>6.960 ± 0.001</td><td>5.118 ± 0.001</td><td>120.951 ± 0.453</td><td>511.21 ± 2.04</td></tr><tr><td>T [JD]</td><td>2 456 072.4 ± 0.6</td><td>2 455 679.9 ± 0.4</td><td>2 455 986.3 ± 0.3</td><td>2 456 231.22 ± 4.26</td><td>2 456 135.92 ± 21.23</td></tr><tr><td>e</td><td>0.16 ± 0.08</td><td>0.31 ± 0.08</td><td>0.27 ± 0.09</td><td>0.35 ± 0.10</td><td>0.16 ± 0.07</td></tr><tr><td>ω [deg]</td><td>343.33 ± 62.12</td><td>109.16 ± 20.16</td><td>328.58 ± 17.78</td><td>254.62 ± 15.91</td><td>291.68 ± 35.85</td></tr><tr><td>K [m s<sup>-1</sup>]</td><td>49.29 ± 5.50</td><td>37.35 ± 4.55</td><td>50.47 ± 3.90</td><td>56.94 ± 4.26</td><td>45.48 ± 3.65</td></tr><tr><td>m sin i [M<sub>Jup</sub>]</td><td>0.42 ± 0.05</td><td>0.33 ± 0.03</td><td>0.40 ± 0.35</td><td>1.57 ± 0.11</td><td>2.18 ± 0.17</td></tr><tr><td>γ [km s<sup>-1</sup>]</td><td>33.179 ± 0.004</td><td>34.184 ± 0.003</td><td>34.058 ± 0.004</td><td>33.188 ± 0.019</td><td>34.567 ± 0.008</td></tr><tr><td>χ<sup>2</sup><sub>red</sub></td><td>0.97</td><td>0.95</td><td>0.93</td><td>1.08</td><td>1.57</td></tr><tr><td>σ(O-C) [m s<sup>-1</sup>]</td><td>12.31</td><td>11.51</td><td>14.24</td><td>15.93</td><td>12.90</td></tr></tbody></table>

**Notes.** P: period, T: time at periastron passage, e: eccentricity, ω: argument of periastron, K: semi-amplitude of the RV curve, m sin i: planetary minimum mass, γ: average radial velocity, and σ(O-C): dispersion of the Keplerian fit residuals.

**References.** (a) Yadav et al. (2008). (b) Montgomery et al. (1993). (c) Sanders (1977). (d) Pietrinferni et al. (2004) and Girardi et al. (2000). (e) Önehag et al. (2011). (f) Pasquini et al. (2008) and Pace et al. (2012). (g) Smolinski et al. (2011) and Lee et al. (2008). (h) Wu et al. (2011). (i) Jacobson et al. (2011).

in the CaII K line of $\langle F'_{K} \rangle \sim 2.1 \times 10^6$ erg cm$^{-2}$ s$^{-1}$ and $\langle F'_{K} \rangle \sim 2.43 \times 10^6$ erg cm$^{-2}$ s$^{-1}$, respectively, while M67 stars have $\langle F'_{K} \rangle \sim 0.5 \times 10^6$ erg cm$^{-2}$ s$^{-1}$, which is not enough to explain the high RV variations we observe. In a survey of the CaII H and K core strengths of a sample of 60 solar-type stars in M67, Giampapa et al. (2006) found that the distribution of the HK index (a measure of the strength of the chromospheric H and K cores) is broader than the distribution seen in the contemporary solar cycle. Significant overlap between the HK distribution of the solar cycle and that for the Sun-like stars in M67 is seen with over 70% of the solar analogs exhibiting CaII H+K strengths within the range of the modern solar cycle. About ~10% are characterized by high activity in excess of the solar maximum values, while approximately 17% have values of the HK index lower than the solar minimum. Of these, none of the stars showing enhanced activity is present in our final sample. In a following work, the same authors reported the results of the analysis of high-resolution photospheric line spectra obtained with the UVES instrument on the VLT for a subset of 15 solar-type stars in M67 selected by Giampapa et al. (2006). They found upper limits to the projected rotation velocities that are consistent with solar-like rotation (i.e., $v \sin i \le 2-3$ km s$^{-1}$) for objects with CaII chromospheric activity within the range of the contemporary solar cycle. Two solar-type stars in their sample exhibit chromospheric emission well in excess of even the solar maximum values: Sand747 and Sand1452. In one case, Sand747, the authors found it to be a spectroscopic binary. The other star, Sand1452, was also present in our original sample, but we discovered that the object was also a binary system (Pasquini et al. 2012). Furthermore, of the 15 solar-type stars analyzed in Reiners & Giampapa (2009), seven stars are in common with our sample (in particular our planet-hosts YBP1194 and YBP1514). These objects are slow rotators ($v \sin i \le 2$ km s$^{-1}$) and have Sun-like HK values.

Melo et al. (2001) used FEROS spectrograph observations to determine accurate projected rotational velocities $v \sin i$ for a sample of 28 MS, turn-off, and giant stars belonging to M67. They found that the stars show similar values of $v \sin i$ depending on their position in the CMD. Early MS G stars have a rotational velocity twice higher than the Sun, and they show a possible trend with ($B-V$) color, with redder colors corresponding to lower $v \sin i$. The stars at the turn-off are the fastest rotators, with $v \sin i$ between 6.3 and 7.6 km s$^{-1}$, while stars just above the turn-off are already significantly slower, with values between 4.6 and 4.9 km s$^{-1}$. Along the subgiant branch rotation tends to drop down, and for stars with ($B-V$) > 1 only upper limits can be found, including for the clump stars ($v \sin i \le 1.5$ km s$^{-1}$). Most of the stars in this last group are in common with our sample.

To rule out activity-related rotational modulation as the cause of the RV variations in our data, we investigated chromospheric activity by measuring the variations of the core of $H\alpha$ with respect to the continuum. Using public data with known correlations between activity indicators and RVs, we verified that such correlations could be still recorded at low S/N level of the M67 spectra. Of our targets, S364 and S978 show a variability in $H\alpha$ of 2%, YBP1514 and YBP1194 of 3% without significant periodicity and YBP401 of 4%, exhibiting all a very low activity level, while for YBP778 $H\alpha$ variability is slightly higher, with a value of 7%. The fact that these stars are of solar age and that our research is focused on finding giant planets with an expected RV variability of tens of m s$^{-1}$ makes the contamination by activity negligible.

### 4.3. Stars with long-term RV variability

The turn-off star S815 shows a peak-to-peak RV variation of the order of ~700 m s$^{-1}$ (see Fig. 6). This star is retained in the single-star sample, although the amplitude of RV values is
<figure><img src="image_6.png" /><figcaption>Fig. 6. RV time series for S815 (left) and YBP2018 (right). For this latter a tentative Keplerian solution is overplotted in blue. Same symbols as in Fig. 3.</figcaption></figure>

possibly too high for a planet. Moreover, the RV measurements would indicate a rather high eccentricity.

The MS star YBP2018 presents peak-to-peak RV variations on the order of ~400 m s⁻¹ and appears to have almost completed one orbit. Figure 6 shows the RV measurements with a tentative Keplerian solution, resulting in an orbiting object with a minimum mass of ~11.0 M<sub>Jup</sub> and a period of ~2487 days. The residuals have an rms of ~30.0 m s⁻¹, but the periodogram of the residuals does not show any clear periodicity. Extensive follow-up over several years is required for this star to understand the nature of the companion candidate.

The MS stars YBP1062 and YBP1137 exhibit a trend in RV measurements of ~28.75 m s⁻¹/yr and ~10.57 m s⁻¹/yr. Figures 7 and 8 show the RV time series with the fitted linear trend overlaid. Although the residuals of the linear fit present an RV variability of ~33.36 m s⁻¹ for YBP1062 and ~18.00 m s⁻¹ for YBP1137, the periodogram of the residuals for the two stars does not reveal any significant peaks.

The evolved star S488 also shows a high RV variability when compared to the measurement errors but caution is necessary with this star, which is located at the tip on the red giant branch. Dupree et al. (1999) found a strong Ca II emission-line with no visible change in the line asymmetries with time for this star, suggesting the presence of outward mass motions.

When we compute a Lomb-Scargle periodogram, a peak is present at ~2257 days. At this period, the measurements baseline is not yet long enough to constrain any significant orbital solution. A tentative Keplerian curve is given in Fig. 9 and corresponds to an object with a minimum mass of ~14.0 M<sub>Jup</sub> and a period of ~2332 days. The residuals have an rms amplitude of ~55 m s⁻¹, and when the main signal is removed, the periodogram of the residual shows a peak at ~95.03 days with a lower significance close to a 0.05 FAP level. Additional observations are needed to draw further conclusions.

## 5. Planet frequency

A series of simulations based on a Monte Carlo approach have been carried out to determine detection limits for the RV data as a function of planet mass and planet period, and to derive a trustworthy estimate of the occurrence rate of giant planets for our RV survey. Understanding the frequency of different types of planets around stars of different mass can provide important clues about the processes of planet formation and evolution.

Moreover, the direct comparison with similar analyses on field star samples can reveal important indications about the influence of the stellar birth environment on the evolution of planetary systems.

In general, the determination of the detection efficiency for an RV survey is based on the ability to recover a planetary signal with a given level of statistical significance. This translates into lower limits on the detectable companion mass as a function of orbital period. Several authors (e.g., Cumming et al. 2008; Howard et al. 2011) have presented statistical analyses of their planet surveys and discussed different methods to derive limits on companion detectability and unbiased distributions for a number of RV planet surveys. They used mainly two approaches: one based on χ² and F-tests to detect excess residuals above an assumed Gaussian noise (Lagrange et al. 2009; Sozzetti et al. 2009), and the other approach is based on a periodogram analysis to identify significant periodicity (Cumming 2004; Narayan et al. 2005; Mortier et al. 2012). For our project, the second approach was selected for all the stars with enough data points for a reliable periodogram analysis. For the stars where the periodogram was not feasible because of the combined effect of poor sampling and the small number of observations per star, we evaluated the detection limit using the variability of the RV values with respect to the measurement errors.

### 5.1. Method

In order to derive the detectability of planetary signals, we computed synthetic datasets, simulating a series of mass and period data pairs (~10⁶ values) of our “potential planets”, which are uniformly distributed in a logarithmical mass range of 0.2–10.0 M<sub>Jup</sub> and in a linear period range of 1.0–1000 days. From the real measurements of each star we retain the observation dates (expressed in Baricentric Julian Date) and the corresponding measurement errors, which we inflated as explained in Sect. 3. For random choices of the planet mass-period pairs and assuming eccentric orbits, we calculated the contribution to the radial velocity amplitude K from a giant planet using the relation

$$K = \left(\frac{2\pi G}{P}\right)^{1/3} \left(\frac{m_p \sin(i)}{M_{Jup}}\right) \left(\frac{M_{\star}}{M_{\odot}}\right)^{-2/3} \left(\frac{1}{\sqrt{1 - e^2}}\right), \quad (1)$$

where M<sub>★</sub> is the mass of the host star that we obtained from isochrone fitting. A random distribution of the orbit inclination
<figure><img src="image_7.png" /><figcaption>Fig. 7. Left: RV time series with the fitted linear trend of YBP1062 overlaid. Right: periodogram of the residuals after the RV trend has been removed. Same symbols as in Fig. 3.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 8. Left: RV time series with the fitted linear trend of YBP1137 overlaid. Right: periodogram of the residuals after the RV trend has been removed. Same symbols as in Fig. 3.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 9. Phased RV measurements and Keplerian best fit (left), periodogram of the RV data (center), and periodogram of the residuals for S488.</figcaption></figure>

was also assumed, and the eccentricity $e$ was allowed to vary randomly between $0.0 < e < 0.5$. Then, we derived the synthetic planetary signals following the equation for eccentric Keplerian orbits,

$$V_r(t) = a \cos v(t) + b \sin v(t) + c, \qquad (2)$$

where $a = K \cos \omega$, $b = -K \sin \omega$ and $c = Ke \cos \omega + \gamma$. Here, $K$ is the RV amplitude, $e$ the eccentricity, $\omega$ the longitude of periastron, $\gamma$ the systemic velocity of the system, $v(t)$ the true anomaly, which is a function of $t$, $P$, and $e$, and $t_0$ is the time of periastron passage. The latter was selected randomly in the time span of the observations, and $\gamma$ was allowed to vary around the RV value

of the M67 cluster obtained in Pasquini et al. (2012). We then degraded the obtained RVs by adding noise corresponding to the actual measurement error of the real observations.

### 5.1.1. Stars with at least ten data points

A periodogram-based analysis (Scargle 1982) was applied to each synthetic dataset to verify whether the simulated planet was observable. This approach was chosen for all the stars with enough measurements for a reliable periodogram analysis. We
note that this is the same procedure as we applied for the analysis of the RV measurements in Sect. 3.

For a given mass and period values, we considered a planet to be detected if we could obtain a signal from the periodogram analysis that had a power higher than the power associated with a 0.01 FAP. For each star $10^6$ trials were used.

In the two-dimensional space of orbital period and planet mass, we consider a grid of period and log-spaced mass cells within which the planet detectability (or detection efficiency, $E_{ij}$) is individually computed as the number of the potentially detected planets with respect to the total simulated planets in each bin of mass and period. On the base of this analysis, the detection efficiency was carried out for 14 MS stars and 10 evolved stars. Figures A.1 and A.2 show the detection efficiencies for the 24 stars as color contours from 0.0 to 1.0. The red regions have a 100% detectability of planets, and blue regions have low detectability. This means that in red regions (high detectability) all potential planets at the given period and mass would be detected. In the blue regions (low detectability) the detection efficiency decreases, and we are not able to discover their presence at high confidence. In Figs. A.1 and A.2, the resulting distribution of the planet detection efficiency clearly reduced with decreasing planet mass and with increasing period. This is due to the insufficient number of measurements, to the distribution of the observations, and to the weak RV signal. Moreover, windows of poor detectability in the period-mass grid occur for two main reasons: the data structure (number of observations and phase coverage), and the one-year seasonal period that affects any ground-based observing program.

### 5.1.2. Stars with poor sampling and few observations

In order to take into account the contribution of stars with poor sampling and few observations, we derived the number of detectable planets using the measurements rms. We considered the same grid of period and log-spaced mass cells as described in Sect. 5.1. In each cell of the mass-period grid, we evaluated the average rms of the synthetic radial velocities for the simulated planets (RMS<sub>sim</sub>). We set our detection threshold by identifying for each star the mass-period cells with synthetic average RV variation more than three times higher than the estimated errors (RMS<sub>sim</sub> > 3σ<sub>corr</sub>) of our real measurements. Inflated measurements errors were used in the analysis.

All potential planets with RMS<sub>sim</sub> smaller than 3σ<sub>corr</sub> were instead considered undetectable.

Finally, we folded the results of all the stars by obtaining the mass-period cell-matrix M. The value for each $M_{ij}$ cell represents the total number of poorly sampled stars for which planets in that cell would be detected.

As an example, the detection regions for a star of the sample derived using the RV rms analysis is shown in Fig. 10. In the red area the rms of the simulated RV values (RMS<sub>sim</sub>) is more than three times higher than the real measurements error (σ<sub>corr</sub>): all hypothetical planets in this mass-period range would be detected. In the blue ruled area the rms of the simulated RV values (RMS<sub>sim</sub>) is lower than 3σ<sub>corr</sub>: potential planets in this mass-period range cannot be detected with our observations.

### 5.1.3. Planetary occurrence

We computed the planet occurrence rate γ(P, M), that is, the fraction of stars of our survey that are orbited by giant planets in the selected period-mass ranges, using the following formalism.

<figure><img src="image_10.png" /><figcaption>Fig. 10. Planet detection regions for a star of the sample derived using the RV rms analysis. Red area: the rms of the simulated RV values is more than three times higher than the real measurement error, thus all potential planets can be detected. Blue ruled area: the rms of the simulated RV values is lower than three times the real measurement errors. Potential planets in this mass-period range cannot be detected with our observations.</figcaption></figure>

Considering the 24 stars analyzed with the Lomb-Scargle method, we derived the completeness, expressed as detectable planetary number, $R_{\star k}$, for each k star. First we weighted the calculated detection efficiency with the mass-period distribution of planets derived from the relation of Cumming et al. (2008): $d f \propto M^{-0.31 \pm 0.2} P^{0.26 \pm 0.1} d \log M d \log P$. Therefore we defined

$$R_{\star k} = c \sum_{ij}^{N_{\text{bins}}} E_{ij} \cdot D_{ij}, \quad (3)$$

where the sum is evaluated over the entire mass-period bins. $E_{ij}$ corresponds to the detection efficiency for each period-mass cell, $D_{ij}$ represents the mass-period distribution of planets within each bin and is calculated by

$$D_{ij} = \int_{\text{bin}} M^{-0.31 \pm 0.2} P^{0.26 \pm 0.1} d \log M d \log P. \quad (4)$$

The normalization constant c is defined conveniently inside our mass-period domain by the relation

$$1 = c \int_{1.0}^{1000} \int_{0.2}^{10} M^{-0.31 \pm 0.2} P^{0.26 \pm 0.1} d \log M d \log P. \quad (5)$$

Finally, we summed the values obtained from the analysis of all the MS and giant stars (G) as

$$R_{\text{MStot}} = \sum_{\text{MS stars}} R_{\star k}, \quad R_{\text{Gtot}} = \sum_{\text{G stars}} R_{\star k}. \quad (6)$$

$R_{\text{MStot}}$ and $R_{\text{Gtot}}$ can be interpreted as the number of planets that would be detected if every star in the sample were hosting a planet. We imposed the result to be proportional to the number of real detected planets in the seven-year-long RV survey (three MS and two giant stars):

$$3 = \gamma_{\text{MS}} R_{\text{MStot}}, \quad 2 = \gamma_{\text{G}} R_{\text{Gtot}}. \quad (7)$$

Therefore, the proportional constants γ<sub>MS</sub> and γ<sub>G</sub> represent an estimate of the planet occurrence rate of our survey for the MS and G stars respectively. In the analysis, turn-off stars with
12.5 ≤ $M_V$ ≤ 13.5 and 0.540 ≤ ($B-V$) ≤ 0.640 were considered part of the MS sample, and the others were considere part of the G sample.

Subsequently, we decided to include a correction for the total detection efficiency to take the contribution of the stars with few observations and poor sampling into account in computing the planet occurrence rate.

Two cases were evaluated.

In the first, the planet occurrence among the MS stars ($N_{MS} = 14$) and the evolved stars ($N_G = 10$) was rescaled with the ratio of $N_{MS}$ or $N_G$ over the total number of stars considered in the analysis: 14 of 56 for the MS and 10 of 20 for the giant stars. This corresponds to the more conservative approach, where all the stars with few observations do not have a companion.

In the second case, we instead decided to rescale the detection efficiency for each cell $E_{ij}$ as

$$E'_{ij} = \left( \frac{M_{ij} + N_{\text{stars}}}{N_{\text{stars}}} \right) \cdot E_{ij}, \qquad (8)$$

where $M_{ij}$ is the number of stars with potentially detectable planets for each mass-period cell calculated using the method based on the RV value rms. In this case, only stars with few observations and poor sampling were considered. $N_{\text{stars}}$ is the number of the stars with potentially detectable planets obtained with the periodogram-based analysis.

## 5.2. Discussion

When we consider only the 14 MS and the 10 evolved stars described in the previous section, the planet occurrence in our survey results in the values of $\gamma_{MS} \sim 58.3_{-31.1}^{+56.4}\%$ for the MS and $\gamma_G \sim 43.8_{-28.5}^{+57.2}\%$ for giants. The errors on the number of found planets are calculated following the prescription of Gehrels (1986). Combining the values for MS and giants gives a total occurrence rate of $\gamma_{tot} \sim 51.5_{-22.7}^{+35.0}\%$. These 24 objects represent the fraction of the stars with a detection efficiency high enough for the number of observations to observe giant planets in the selected range of mass and period (see Figs. A.1 and A.2).

Applying then the correction for the detectability based on those stars that have fewer observations, we obtained the following values for the two considered approaches. In the more conservative case, without further detections, the planet occurrence becomes $\sim 14.6_{-8.0}^{+14.1}\%$ for the MS stars, $\sim 21.9_{-14.2}^{+28.6}\%$ for the giant stars, and $\sim 16.3_{-7.0}^{+11.0}\%$ for the whole sample.

In the more general case, the procedure previously described instead gives a giant planet occurrence of $\sim 15.8_{-8.4}^{+15.3}\%$ for the MS stars and $\sim 23.0_{-15.0}^{+29.9}\%$ for the evolved stars. When referring to the whole sample, the occurrence is $\sim 18.0_{-8.0}^{+12.0}\%$.

The values from the latter more general analysis are finally considered as the final results of our study. Interestingly, they are similar to the results of RV surveys around FGK field stars that show exoplanet rates of $\sim 13\%$ for Jupiter-mass stars in approximately the same range of periods (Mayor et al. 2011; Cumming et al. 2008).

Moreover, when we investigated only the frequency of hot Jupiters around MS stars of our survey and exactly repeated the procedure described in the previous sections, but focusing on objects with periods shorter than 10 days, we found a hot Jupiter occurrence rate of $\sim 5.1_{-3.0}^{+4.9}\%$ for the conservative approach and of $\sim 5.7_{-3.0}^{+5.5}\%$ for the more general case. These values appear to be higher than what is observed in field stars (Wright et al. 2012; Howard et al. 2011). About this argument in particular, we refer

to the discussion presented in our parallel work (Brucalassi et al. 2016). We pointed out in that paper that the high rate of hot Jupiters could be favored by a formation scenario dominated by strong encounters with other stars or binary companions and subsequent planet-planet scattering, as predicted by N-body simulations. In this context, we should consider that rescaling the actual M 67 detectability with the mass-period distribution derived by Cumming et al. (2008) for the field stars as done in the previous section could produce a biased result. However, future investigations and a larger number of planet detections in stars cluster are required to assess this hypothesis.

Finally, it is worth to point out that we did not find any clear detection of long-period giant planets around MS stars, possibly because the sensitivity of our survey decreases heavily for long-period planets. However, the recent discovery in the Praesepe cluster of a second massive outer planet around Pr0211 that hosts a close hot Jupiter on a slightly eccentric orbit (Malavolta et al. 2016) suggests that further long-term monitoring could yield indications of interesting discoveries in OCs.

# 6. Summary and conclusion

We have presented the results of a long-term search program for giant planets in the solar-age, solar-metallicity open cluster M 67.

Five different instruments were used, and after finding proper zero-point corrections to HARPS, 1145 observations for 88 stars were analyzed.

The problem of combining long-term precision RV data from different instruments complicates the analysis. Zero-points offsets were derived with a limited precision, which leads to a loss of sensitivity for trends and long periods. However, this problem quite frequently occurs in long-term surveys because spectrographs receive upgrades or survey projects are transferred to new instruments. Long-term access to the same telescope or instrument configuration is therefore quite important for this type of studies.

Five stars of our sample are found to host planets. One of them (YBP401) has recently been presented in Brucalassi et al. (2016), and in that work we have refined the parameters, based on new measurements, of other two planets announced previously around YBP1194 and YBP1514 (Brucalassi et al. 2014). The other planet-host candidate S978 was discussed in this paper. We have no clear additional planet detection in our sample, but some promising or controversial cases such as YBP778 and YBP2018 call for follow-up observations. Fourteen new binaries were identified and added to the catalog of known M 67 binaries that has been created in Paper I (Pasquini et al. 2012). Moreover, we see trends in S1062 and S1137, and large RV scattering is exhibited by S815 and S488, which might be explained by stellar companions or substellar objects. However, stellar chromospherical activity and magnetic cycles may cause RV variations that might be mistaken as companions. Activity indicators such as the Ca II H&K lines, $H\alpha$, and the bisector span are fundamental to verify claimed planets from RV variability.

We provided an estimate of the planet occurrence rate $\gamma(P, M)$, that is, the fraction of stars of our survey that is orbited by Jupiter planets in the ranges of period between 1.0 day and 1000 days and has a planet mass between $0.2 M_J$ and $10.0 M_J$. Although one of the main problems of this survey was the poor and sparse sampling and the small number of observations per star, which complicates the statistical analysis and increases the uncertainty on the frequency of Jupiters in our survey, we find
a total occurrence rate of Jupiter-mass planets in the M 67 sample of ~18.0<sup>+12.0</sup><sub>-8.0</sub>% (precisely ~15.8<sup>+15.3</sup><sub>-8.4</sub>% for the MS stars and ~23.0<sup>+29.9</sup><sub>-15.0</sub> for the evolved stars). It is worth noting that on average we need a minimum number of 20 observations per star to exclude Jupiter-mass planets at high confidence when a periodogram analysis is applied.

The Jupiter-mass planets in our sample imply that our results are consistent with the planet frequency of much larger surveys carried out for field star samples. For example, from the ELODIE survey, Naef et al. (2005) estimated that a fraction of 7.5 ± 1.5% stars host giant planets with periods shorter than 10 yr, while Cumming et al. (2008) derived a frequency of 12 ± 1.6% from the Keck survey and Mayor et al. (2011) 13.9 ± 1.7% from the HARPS/CORALIE survey.

As a general conclusion, our simulation study seems to confirm the recent finding (Meibom et al. 2013) that the frequency of massive planets around stars of open clusters agrees with the frequency of these planets around field stars.

However, we note that when we only investigate the frequency of hot Jupiters around MS stars of our survey, we find a hot Jupiter occurrence rate (~5.7<sup>+5.5</sup><sub>-3.0</sub>%) that is substantially higher than what is observed in field stars (see also Brucalassi et al. 2016).

We have shown that the search for planets in open clusters is a powerful test benchmark for the theory of planet formation and stellar evolution. The new generation of spectrographs such as ESPRESSO (Pasquini et al. 2009) will extend this search more effectively. It is clear, however, that to efficiently perform a planet search survey in star clusters, a multi-object capability, even with relatively low multiplex (10–50), is desirable. Unfortunately, no such facility is currently planned.

*Acknowledgements.* L.P. acknowledges the Visiting Researcher program of the CNPq Brazilian Agency, at the Fed. Univ. of Rio Grande do Norte, Brazil. RPS thanks ESO DGDF, the whole GAPS-TNG community, the Hobby Eberly Telescope (HET) project, the PNPS and PNP of INSU – CNRS for allocating the observations and the technical support. M.T.R. received support from PFB06 CATA (CONICYT). We are grateful to Gaspare Lo Curto and Michele Cappetta for the support in data reduction analysis and helpful discussions.

## References

Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&A, 119, 373

Baumann, P., Ramírez, I., Meléndez, J., et al. 2010, A&A, 519, A87

Bouchy, F., & Sophie Team 2006, in Tenth Anniversary of 51 Peg-b, eds. L. Arnold, F. Bouchy, & C. Moutou, 319

Brucalassi, A., Pasquini, L., Saglia, R., et al. 2014, A&A, 561, L9

Brucalassi, A., Pasquini, L., Saglia, R., et al. 2016, A&A, 592, L1

Cai, M. X., Spurzem, R., & Kouwenhoven, M. B. N. 2016, IAU Symp., 312, 235

Cappetta, M., Saglia, R. P., Birkby, J. L., et al. 2012, MNRAS, 427, 1877

Cumming, A. 2004, MNRAS, 354, 1165

Cumming, A., Butler, R. P., Marcy, G. W., et al. 2008, PASP, 120, 531

Davies, M. B., Adams, F. C., Armitage, P., et al. 2014, Protostars and Planets VI, 787

Dawson, R. I., & Murray-Clay, R. A. 2013, ApJ, 767, L24

Deal, M., Richard, O., & Vauclair, S. 2015, A&A, 584, A105

Dong, S., Katz, B., & Socrates, A. 2014, ApJ, 781, L5

Dupree, A. K., Whitney, B. A., & Pasquini, L. 1999, ApJ, 520, 751

Gehrels, N. 1986, ApJ, 303, 336

Geller, A. M., Latham, D. W., & Mathieu, R. D. 2015, AJ, 150, 97

Giampapa, M. S., Hall, J. C., Radick, R. R., & Baliunas, S. L. 2006, ApJ, 651, 444

Gilliland, R. L., Brown, T. M., Guhathakurta, P., et al. 2000, ApJ, 545, L47

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371

González Hernández, J. I., Israelian, G., Santos, N. C., et al. 2010, ApJ, 720, 1592

González Hernández, J. I., Delgado-Mena, E., Sousa, S. G., et al. 2013, in New Quests in Stellar Astrophysics III: A Panchromatic View of Solar-Like Stars, With and Without Planets, eds. M. Chavez, E. Bertone, O. Vega, & V. De la Luz, ASP Conf. Ser., 472, 97

Hekker, S., Snellen, I. A. G., Aerts, C., et al. 2008, A&A, 480, 215

Horne, J. H., & Baliunas, S. L. 1986, ApJ, 302, 757

Howard, A. W., Marcy, G. W., Bryson, S. T., & et al. 2011, ApJ, 730, 10

Jacobson, H. R., Pilachowski, C. A., & Friel, E. D. 2011, AJ, 142, 59

Johnson, J. A., Aller, K. M., et al. 2010, PASP, 122, 905

Lagrange, A.-M., Desort, M., Galland, F., Udry, S., & Mayor, M. 2009, A&A, 495, 335

Lee, Y. S., Beers, T. C., Sivarami, T., et al. 2008, AJ, 136, 2050

Li, G., & Adams, F. C. 2015, MNRAS, 448, 344

Lovis, C., & Mayor, M. 2007, A&A, 472, 657

Malvolta, L., Nascimbeni, V., Piotto, G., et al. 2016, A&A, 588, A118

Malmberg, D., Davies, M. B., & Heggie, D. C. 2011, MNRAS, 411, 859

Mathieu, R. D., Latham, D. W., Griffin, R. F., & Gunn, J. E. 1986, AJ, 92, 1100

Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20

Mayor, M., Marmier, M., Lovis, C., et al. 2011, ArXiv e-prints [arXiv:1109.2497]

Meibom, S., Torres, G., Fressin, F., et al. 2013, Nature, 499, 55

Meléndez, J., Asplund, M., Gustafsson, B., & Yong, D. 2009, ApJ, 704, L66

Melo, C. H. F., Pasquini, L., & De Medeiros, J. R. 2001, A&A, 375, 851

Mermilliod, J.-C., & Mayor, M. 2007, A&A, 470, 919

Montgomery, K. A., Marschall, L. A., & Janes, K. A. 1993, AJ, 106, 181

Mortier, A., Santos, N. C., Sozzetti, A., et al. 2012, A&A, 543, A45

Naef, D., Mayor, M., Beuzit, J.-L., et al. 2005, in 13th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun, eds. F. Favata, G. A. J. Hussain, & B. Battrick, ESA SP, 560, 833

Narayan, R., Cumming, A., & Lin, D. N. C. 2005, ApJ, 620, 1002

Nascimbeni, V., Bedin, L. R., Piotto, G., De Marchi, F., & Rich, R. M. 2012, A&A, 541, A144

Önehag, A., Korn, A., Gustafsson, B., et al. 2011, A&A, 528, A85

Önehag, A., Gustafsson, B., & Korn, A. 2014, A&A, 562, A102

Pace, G., & Pasquini, L. 2004, A&A, 426, 1021

Pace, G., Pasquini, L., & François, P. 2008, A&A, 489, 403

Pace, G., Castro, M., Meléndez, J., Théado, S., & do Nascimento, J. J.-D. 2012, A&A, 541, A150

Pasquini, L., & Pallavicini, R. 1991, A&A, 251, 199

Pasquini, L., Biazzo, K., Bonifacio, P., et al. 2008, A&A, 489, 677

Pasquini, L., Manescau, A., Avila, G., et al. 2009, in Science with the VLT in the ELT Era, ed. A. Moorwood, 395

Pasquini, L., Brucalassi, A., Ruiz, M. T., et al. 2012, A&A, 545, A139

Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632

Perruchot, S., Bouchy, F., et al. 2011, in SPIE Conf. Ser., 8151

Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168

Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279

Quinn, S. N., White, R. J., Latham, D. W., et al. 2012, ApJ, 756, L33

Quinn, S. N., White, R. J., Latham, D. W., et al. 2014, ApJ, 787, 27

Ramírez, I., Asplund, M., Baumann, P., Meléndez, J., & Bensby, T. 2010, A&A, 521, A33

Randich, S., Sestito, P., Primas, F., et al. 2006, A&A, 450, 557

Reiners, A., & Giampapa, M. S. 2009, ApJ, 707, 852

Sanders, W. L. 1977, A&AS, 27, 89

Santerne, A., Díaz, R. F., Moutou, C., et al. 2012, A&A, 545, A76

Santos, N. C., Gomes da Silva, J., Lovis, C., & Melo, C. 2010, A&A, 511, A54

Sato, B., Izumiura, H., Toyota, E., et al. 2007, ApJ, 661, 527

Scargle, J. D. 1982, ApJ, 263, 835

Setiawan, J., Pasquini, L., da Silva, L., et al. 2004, A&A, 421, 241

Shara, M. M., Hurley, J. R., & Mardling, R. A. 2016, ApJ, 816, 59

Smolinski, J. P., Lee, Y. S., Beers, T. C., et al. 2011, AJ, 141, 89

Sozzetti, A., Torres, G., Latham, D. W., et al. 2009, ApJ, 697, 544

Spurzem, R., Giersz, M., Heggie, D. C., & Lin, D. N. C. 2009, ApJ, 697, 458

Tull, R. G., MacQueen, P. J., Good, J., Epps, H. W., & HET HRS Team 1998, BAAS, 30, 1263

Wright, J. T., & Howard, A. W. 2009, ApJS, 182, 205

Wright, J. T., Marcy, G. W., Howard, A. W., et al. 2012, ApJ, 753, 160

Wu, Y., Luo, A.-L., Li, H.-N., et al. 2011, RA&A, 11, 924

Yadav, R. K. S., Bedin, L. R., Piotto, G., et al. 2008, A&A, 484, 609
Appendix A: Accompanying plots and table

<table><caption>Table A.1. Observed targets in M 67.</caption><thead><tr><th>Object</th><th>B-V</th><th>M<sub>V</sub></th><th>H</th><th>S</th><th>C</th><th>HET</th><th>H-N</th><th>TOT</th><th>RV(km s<sup>-1</sup>)</th><th>σ<sub>obs</sub>(km s<sup>-1</sup>)</th><th>σ<sub>cor</sub>(km s<sup>-1</sup>)</th><th>σ<sub>RV</sub>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td><strong>YBP219</strong></td><td><strong>0.570</strong></td><td><strong>13.6</strong></td><td><strong>6</strong></td><td><strong>5</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>11</strong></td><td><strong>33.714</strong></td><td><strong>0.016</strong></td><td><strong>0.022</strong></td><td><strong>0.323</strong></td></tr><tr><td>YBP266</td><td>0.570</td><td>13.6</td><td>11</td><td>3</td><td>0</td><td>0</td><td>1</td><td>15</td><td>33.774</td><td>0.017</td><td>0.022</td><td>0.031</td></tr><tr><td>YBP285</td><td>0.663</td><td>14.5</td><td>7</td><td>2</td><td>0</td><td>0</td><td>0</td><td>9</td><td>34.397</td><td>0.016</td><td>0.022</td><td>0.019</td></tr><tr><td><strong>YBP288</strong></td><td><strong>0.637</strong></td><td><strong>13.9</strong></td><td><strong>1</strong></td><td><strong>2</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>3</strong></td><td><strong>37.691</strong></td><td><strong>0.010</strong></td><td><strong>0.018</strong></td><td><strong>1.299</strong></td></tr><tr><td>YBP291</td><td>0.570</td><td>13.5</td><td>23</td><td>3</td><td>0</td><td>0</td><td>0</td><td>26</td><td>32.521</td><td>0.017</td><td>0.022</td><td>0.026</td></tr><tr><td>YBP349</td><td>0.636</td><td>14.3</td><td>11</td><td>0</td><td>0</td><td>0</td><td>0</td><td>11</td><td>35.048</td><td>0.014</td><td>0.020</td><td>0.023</td></tr><tr><td>YBP350</td><td>0.561</td><td>13.6</td><td>11</td><td>2</td><td>0</td><td>0</td><td>0</td><td>13</td><td>33.227</td><td>0.015</td><td>0.021</td><td>0.020</td></tr><tr><td>YBP401</td><td>0.566</td><td>13.7</td><td>19</td><td>5</td><td>0</td><td>0</td><td>2</td><td>26</td><td>33.178</td><td>0.015</td><td>0.021</td><td>0.035</td></tr><tr><td>YBP473</td><td>0.658</td><td>14.4</td><td>7</td><td>0</td><td>0</td><td>0</td><td>0</td><td>7</td><td>33.266</td><td>0.016</td><td>0.022</td><td>0.023</td></tr><tr><td>YBP587</td><td>0.605</td><td>14.1</td><td>6</td><td>2</td><td>0</td><td>0</td><td>0</td><td>8</td><td>33.188</td><td>0.013</td><td>0.020</td><td>0.028</td></tr><tr><td>YBP613</td><td>0.612</td><td>13.3</td><td>10</td><td>4</td><td>0</td><td>0</td><td>0</td><td>14</td><td>33.565</td><td>0.015</td><td>0.021</td><td>0.020</td></tr><tr><td>YBP637</td><td>0.661</td><td>14.5</td><td>5</td><td>2</td><td>0</td><td>0</td><td>0</td><td>7</td><td>34.801</td><td>0.014</td><td>0.020</td><td>0.019</td></tr><tr><td><strong>YBP673</strong></td><td><strong>0.665</strong></td><td><strong>14.4</strong></td><td><strong>14</strong></td><td><strong>1</strong></td><td><strong>0</strong></td><td><strong>14</strong></td><td><strong>2</strong></td><td><strong>31</strong></td><td><strong>33.496</strong></td><td><strong>0.026</strong></td><td><strong>0.030</strong></td><td><strong>0.249</strong></td></tr><tr><td>YBP689</td><td>0.622</td><td>13.1</td><td>10</td><td>2</td><td>0</td><td>3</td><td>0</td><td>15</td><td>33.650</td><td>0.017</td><td>0.022</td><td>0.044</td></tr><tr><td>YBP750</td><td>0.598</td><td>13.6</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td><td>34.251</td><td>0.014</td><td>0.020</td><td>0.011</td></tr><tr><td><strong>YBP769</strong></td><td><strong>0.665</strong></td><td><strong>14.4</strong></td><td><strong>2</strong></td><td><strong>1</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>3</strong></td><td>CCF Double peak</td><td></td><td></td><td></td></tr><tr><td>YBP778</td><td>0.582</td><td>13.1</td><td>15</td><td>4</td><td>0</td><td>0</td><td>2</td><td>21</td><td>34.288</td><td>0.016</td><td>0.022</td><td>0.090</td></tr><tr><td>YBP809</td><td>0.696</td><td>15.0</td><td>5</td><td>0</td><td>0</td><td>0</td><td>0</td><td>5</td><td>32.864</td><td>0.011</td><td>0.019</td><td>0.011</td></tr><tr><td><strong>YBP851</strong></td><td><strong>0.617</strong></td><td><strong>14.4</strong></td><td><strong>3</strong></td><td><strong>1</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>4</strong></td><td><strong>33.759</strong></td><td><strong>0.010</strong></td><td><strong>0.018</strong></td><td><strong>1.417</strong></td></tr><tr><td><strong>YBP911</strong></td><td><strong>0.673</strong></td><td><strong>14.6</strong></td><td><strong>2</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>2</strong></td><td><strong>33.738</strong></td><td><strong>0.020</strong></td><td><strong>0.025</strong></td><td><strong>0.703</strong></td></tr><tr><td>YBP988</td><td>0.598</td><td>14.2</td><td>5</td><td>1</td><td>0</td><td>0</td><td>0</td><td>6</td><td>32.862</td><td>0.013</td><td>0.020</td><td>0.016</td></tr><tr><td>YBP1032</td><td>0.598</td><td>14.4</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>6</td><td>34.913</td><td>0.013</td><td>0.020</td><td>0.013</td></tr><tr><td>YBP1036</td><td>0.690</td><td>15.0</td><td>15</td><td>0</td><td>0</td><td>0</td><td>0</td><td>15</td><td>34.025</td><td>0.015</td><td>0.021</td><td>0.025</td></tr><tr><td><strong>YBP1051</strong></td><td><strong>0.595</strong></td><td><strong>14.1</strong></td><td><strong>15</strong></td><td><strong>3</strong></td><td><strong>0</strong></td><td><strong>5</strong></td><td><strong>0</strong></td><td><strong>23</strong></td><td><strong>33.364</strong></td><td><strong>0.016</strong></td><td><strong>0.022</strong></td><td><strong>0.156</strong></td></tr><tr><td>YBP1062</td><td>0.626</td><td>14.5</td><td>19</td><td>3</td><td>0</td><td>0</td><td>0</td><td>22</td><td>33.462</td><td>0.017</td><td>0.022</td><td>0.067</td></tr><tr><td><strong>YBP1067</strong></td><td><strong>0.642</strong></td><td><strong>14.6</strong></td><td><strong>7</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>2</strong></td><td><strong>0</strong></td><td><strong>9</strong></td><td><strong>33.667</strong></td><td><strong>0.019</strong></td><td><strong>0.024</strong></td><td><strong>1.030</strong></td></tr><tr><td>YBP1075</td><td>0.633</td><td>13.7</td><td>11</td><td>1</td><td>0</td><td>0</td><td>1</td><td>13</td><td>33.858</td><td>0.014</td><td>0.020</td><td>0.034</td></tr><tr><td>YBP1088</td><td>0.618</td><td>14.5</td><td>8</td><td>0</td><td>0</td><td>0</td><td>0</td><td>8</td><td>33.434</td><td>0.014</td><td>0.020</td><td>0.022</td></tr><tr><td><strong>YBP1090</strong></td><td><strong>0.650</strong></td><td><strong>13.8</strong></td><td><strong>2</strong></td><td><strong>1</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>3</strong></td><td><strong>35.186</strong></td><td><strong>0.021</strong></td><td><strong>0.025</strong></td><td><strong>1.265</strong></td></tr><tr><td>YBP1101</td><td>0.661</td><td>14.7</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td><td>33.484</td><td>0.015</td><td>0.021</td><td>0.029</td></tr><tr><td>YBP1129</td><td>0.583</td><td>14.2</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td><td>34.479</td><td>0.009</td><td>0.018</td><td>0.009</td></tr><tr><td>YBP1137</td><td>0.657</td><td>14.9</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>13</td><td>34.227</td><td>0.018</td><td>0.023</td><td>0.031</td></tr><tr><td>YBP1194</td><td>0.626</td><td>14.6</td><td>16</td><td>6</td><td>0</td><td>5</td><td>2</td><td>29</td><td>34.189</td><td>0.017</td><td>0.022</td><td>0.027</td></tr><tr><td>YBP1197</td><td>0.565</td><td>13.3</td><td>10</td><td>0</td><td>0</td><td>0</td><td>0</td><td>10</td><td>34.591</td><td>0.017</td><td>0.022</td><td>0.015</td></tr><tr><td>YBP1247</td><td>0.568</td><td>14.1</td><td>5</td><td>1</td><td>0</td><td>0</td><td>0</td><td>6</td><td>32.966</td><td>0.013</td><td>0.020</td><td>0.014</td></tr><tr><td>YBP1303</td><td>0.636</td><td>14.6</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td><td>33.395</td><td>0.017</td><td>0.022</td><td>0.019</td></tr><tr><td><strong>YBP1304</strong></td><td><strong>0.723</strong></td><td><strong>14.7</strong></td><td><strong>2</strong></td><td><strong>1</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>3</strong></td><td><strong>32.512</strong></td><td><strong>0.019</strong></td><td><strong>0.024</strong></td><td><strong>2.670</strong></td></tr><tr><td><strong>YBP1315</strong></td><td><strong>0.693</strong></td><td><strong>14.3</strong></td><td><strong>5</strong></td><td><strong>2</strong></td><td><strong>0</strong></td><td><strong>6</strong></td><td><strong>0</strong></td><td><strong>13</strong></td><td><strong>34.885</strong></td><td><strong>0.022</strong></td><td><strong>0.026</strong></td><td><strong>0.801</strong></td></tr><tr><td>YBP1334</td><td>0.639</td><td>14.4</td><td>6</td><td>1</td><td>0</td><td>0</td><td>0</td><td>7</td><td>33.066</td><td>0.014</td><td>0.020</td><td>0.021</td></tr><tr><td>YBP1387</td><td>0.585</td><td>14.1</td><td>5</td><td>1</td><td>0</td><td>0</td><td>0</td><td>6</td><td>34.059</td><td>0.012</td><td>0.019</td><td>0.017</td></tr><tr><td>YBP1392</td><td>0.675</td><td>14.8</td><td>8</td><td>1</td><td>0</td><td>0</td><td>0</td><td>6</td><td>34.556</td><td>0.016</td><td>0.022</td><td>0.023</td></tr><tr><td>YBP1458</td><td>0.698</td><td>15.0</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>6</td><td>33.417</td><td>0.013</td><td>0.020</td><td>0.015</td></tr><tr><td>YBP1496</td><td>0.556</td><td>13.9</td><td>8</td><td>1</td><td>0</td><td>0</td><td>0</td><td>9</td><td>34.786</td><td>0.015</td><td>0.021</td><td>0.017</td></tr><tr><td>YBP1504</td><td>0.584</td><td>14.2</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td><td>33.761</td><td>0.013</td><td>0.020</td><td>0.018</td></tr><tr><td>YBP1514</td><td>0.680</td><td>14.8</td><td>17</td><td>5</td><td>5</td><td>0</td><td>0</td><td>27</td><td>34.048</td><td>0.018</td><td>0.023</td><td>0.036</td></tr><tr><td>YBP1587</td><td>0.600</td><td>14.2</td><td>11</td><td>1</td><td>0</td><td>0</td><td>1</td><td>13</td><td>33.446</td><td>0.011</td><td>0.019</td><td>0.035</td></tr><tr><td>YBP1622</td><td>0.591</td><td>14.2</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>6</td><td>33.951</td><td>0.012</td><td>0.019</td><td>0.015</td></tr><tr><td><strong>YBP1716</strong></td><td><strong>0.619</strong></td><td><strong>13.3</strong></td><td><strong>5</strong></td><td><strong>1</strong></td><td><strong>0</strong></td><td><strong>7</strong></td><td><strong>0</strong></td><td><strong>13</strong></td><td><strong>36.205</strong></td><td><strong>0.010</strong></td><td><strong>0.018</strong></td><td><strong>0.651</strong></td></tr><tr><td>YBP1722</td><td>0.560</td><td>14.1</td><td>10</td><td>1</td><td>0</td><td>5</td><td>1</td><td>17</td><td>34.483</td><td>0.018</td><td>0.023</td><td>0.025</td></tr><tr><td>YBP1735</td><td>0.620</td><td>14.3</td><td>5</td><td>0</td><td>0</td><td>0</td><td>0</td><td>5</td><td>33.976</td><td>0.013</td><td>0.020</td><td>0.014</td></tr><tr><td><strong>YBP1758</strong></td><td><strong>0.653</strong></td><td><strong>13.2</strong></td><td><strong>2</strong></td><td><strong>1</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>3</strong></td><td><strong>29.653</strong></td><td><strong>0.011</strong></td><td><strong>0.019</strong></td><td><strong>1.521</strong></td></tr><tr><td>YBP1768</td><td>0.615</td><td>14.4</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>34.497</td><td>0.012</td><td>0.019</td><td>0.003</td></tr><tr><td>YBP1787</td><td>0.626</td><td>14.5</td><td>8</td><td>1</td><td>0</td><td>0</td><td>0</td><td>9</td><td>34.065</td><td>0.014</td><td>0.020</td><td>0.016</td></tr><tr><td>YBP1788</td><td>0.622</td><td>14.4</td><td>8</td><td>0</td><td>0</td><td>0</td><td>0</td><td>8</td><td>34.162</td><td>0.014</td><td>0.020</td><td>0.042</td></tr></tbody></table>

Notes. Object name, basic stellar parameters, number of observations for each instrument (H: HARPS, S: SOPHIE, C: CORALE, HET: HRS at Het, H-N: HARPS-N), total number of observations, mean stellar RV, average measurements error, average corrected measurements error, observed RV dispersion. Binary candidates are indicated in bold face.
**Table A.1. continued.**

<table><thead><tr><th>Object</th><th>B-V</th><th>M<sub>V</sub></th><th>H</th><th>S</th><th>C</th><th>HET</th><th>H-N</th><th>TOT</th><th>RV(km s<sup>-1</sup>)</th><th>σ<sub>obs</sub>(km s<sup>-1</sup>)</th><th>σ<sub>cor</sub>(km s<sup>-1</sup>)</th><th>σ<sub>RV</sub>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>YBP1852</td><td>0.572</td><td>14.0</td><td>7</td><td>1</td><td>0</td><td>0</td><td>0</td><td>8</td><td>32.914</td><td>0.013</td><td>0.020</td><td>0.018</td></tr><tr><td>YBP1903</td><td>0.648</td><td>14.7</td><td>5</td><td>0</td><td>0</td><td>0</td><td>0</td><td>5</td><td>33.390</td><td>0.017</td><td>0.022</td><td>0.023</td></tr><tr><td>YBP1948</td><td>0.571</td><td>14.0</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td><td>33.333</td><td>0.012</td><td>0.019</td><td>0.011</td></tr><tr><td>YBP1955</td><td>0.589</td><td>14.2</td><td>6</td><td>1</td><td>0</td><td>0</td><td>0</td><td>7</td><td>33.212</td><td>0.013</td><td>0.020</td><td>0.021</td></tr><tr><td>YBP2018</td><td>0.631</td><td>14.6</td><td>22</td><td>2</td><td>0</td><td>7</td><td>0</td><td>31</td><td>31.991</td><td>0.020</td><td>0.025</td><td>0.089</td></tr><tr><td>S364</td><td>1.360</td><td>9.8</td><td>14</td><td>6</td><td>7</td><td>5</td><td>0</td><td>32</td><td>33.198</td><td>0.010</td><td>0.018</td><td>0.042</td></tr><tr><td>S488</td><td>1.550</td><td>8.9</td><td>15</td><td>11</td><td>0</td><td>8</td><td>3</td><td>37</td><td>32.861</td><td>0.012</td><td>0.019</td><td>0.156</td></tr><tr><td>S602</td><td>0.512</td><td>12.9</td><td>6</td><td>4</td><td>0</td><td>7</td><td>0</td><td>17</td><td>33.894</td><td>0.018</td><td>0.023</td><td>0.047</td></tr><tr><td>S610</td><td>0.493</td><td>12.9</td><td>5</td><td>1</td><td>0</td><td>5</td><td>0</td><td>11</td><td>33.371</td><td>0.024</td><td>0.027</td><td>0.056</td></tr><tr><td>S657</td><td>0.559</td><td>12.3</td><td>8</td><td>0</td><td>0</td><td>0</td><td>0</td><td>8</td><td>33.234</td><td>0.008</td><td>0.017</td><td>0.015</td></tr><tr><td>S731</td><td>0.516</td><td>13.1</td><td>10</td><td>1</td><td>0</td><td>0</td><td>0</td><td>11</td><td>33.097</td><td>0.012</td><td>0.019</td><td>0.021</td></tr><tr><td>S815</td><td>0.497</td><td>12.9</td><td>14</td><td>4</td><td>0</td><td>7</td><td>3</td><td>28</td><td>34.139</td><td>0.019</td><td>0.024</td><td>0.207</td></tr><tr><td>S978</td><td>1.332</td><td>9.7</td><td>16</td><td>9</td><td>10</td><td>3</td><td>2</td><td>40</td><td>34.577</td><td>0.010</td><td>0.018</td><td>0.043</td></tr><tr><td>S989</td><td>1.048</td><td>11.4</td><td>8</td><td>1</td><td>8</td><td>3</td><td>0</td><td>20</td><td>34.795</td><td>0.019</td><td>0.024</td><td>0.024</td></tr><tr><td>S1001</td><td>0.759</td><td>12.4</td><td>3</td><td>1</td><td>0</td><td>0</td><td>0</td><td>4</td><td>33.409</td><td>0.007</td><td>0.017</td><td>0.024</td></tr><tr><td>S1010</td><td>1.069</td><td>10.5</td><td>7</td><td>3</td><td>6</td><td>0</td><td>0</td><td>16</td><td>33.746</td><td>0.009</td><td>0.018</td><td>0.025</td></tr><tr><td>S1016</td><td>1.098</td><td>10.3</td><td>7</td><td>11</td><td>0</td><td>0</td><td>0</td><td>18</td><td>33.872</td><td>0.004</td><td>0.016</td><td>0.072</td></tr><tr><td>S1054</td><td>0.859</td><td>11.2</td><td>8</td><td>1</td><td>9</td><td>0</td><td>0</td><td>18</td><td>33.500</td><td>0.015</td><td>0.021</td><td>0.025</td></tr><tr><td>S1074</td><td>1.111</td><td>10.4</td><td>3</td><td>1</td><td>9</td><td>0</td><td>0</td><td>13</td><td>34.139</td><td>0.014</td><td>0.020</td><td>0.042</td></tr><tr><td>S1084</td><td>1.086</td><td>10.5</td><td>15</td><td>1</td><td>6</td><td>5</td><td>0</td><td>27</td><td>33.900</td><td>0.007</td><td>0.017</td><td>0.021</td></tr><tr><td>S1230</td><td>0.524</td><td>13.1</td><td>4</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td><td>33.773</td><td>0.013</td><td>0.020</td><td>0.005</td></tr><tr><td>S1254</td><td>0.999</td><td>11.5</td><td>10</td><td>3</td><td>5</td><td>0</td><td>0</td><td>18</td><td>32.867</td><td>0.013</td><td>0.020</td><td>0.058</td></tr><tr><td>S1271</td><td>0.506</td><td>12.9</td><td>11</td><td>2</td><td>0</td><td>0</td><td>0</td><td>13</td><td>33.648</td><td>0.016</td><td>0.022</td><td>0.023</td></tr><tr><td>S1279</td><td>1.081</td><td>10.6</td><td>7</td><td>2</td><td>7</td><td>0</td><td>0</td><td>16</td><td>33.384</td><td>0.010</td><td>0.018</td><td>0.021</td></tr><tr><td>S1288</td><td>1.016</td><td>11.3</td><td>5</td><td>1</td><td>6</td><td>0</td><td>0</td><td>12</td><td>33.454</td><td>0.020</td><td>0.025</td><td>0.040</td></tr><tr><td>S1293</td><td>0.565</td><td>12.1</td><td>13</td><td>10</td><td>0</td><td>7</td><td>0</td><td>30</td><td>34.094</td><td>0.010</td><td>0.018</td><td>0.042</td></tr><tr><td>S1305</td><td>0.945</td><td>12.2</td><td>14</td><td>2</td><td>0</td><td>6</td><td>0</td><td>22</td><td>33.964</td><td>0.011</td><td>0.019</td><td>0.020</td></tr><tr><td>S1316</td><td>1.077</td><td>10.6</td><td>9</td><td>2</td><td>6</td><td>4</td><td>0</td><td>21</td><td>32.860</td><td>0.012</td><td>0.019</td><td>0.028</td></tr><tr><td>S1402</td><td>1.109</td><td>10.9</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td><td>33.781</td><td>0.007</td><td>0.017</td><td>0.008</td></tr><tr><td>S1479</td><td>0.682</td><td>10.5</td><td>6</td><td>1</td><td>6</td><td>0</td><td>0</td><td>13</td><td>34.319</td><td>0.011</td><td>0.019</td><td>0.013</td></tr><tr><td>S1557</td><td>1.249</td><td>10.1</td><td>13</td><td>13</td><td>6</td><td>4</td><td>1</td><td>37</td><td>33.797</td><td>0.010</td><td>0.018</td><td>0.138</td></tr><tr><td><b>S1583</b></td><td><b>0.665</b></td><td><b>14.4</b></td><td><b>2</b></td><td><b>1</b></td><td><b>0</b></td><td><b>0</b></td><td><b>0</b></td><td><b>3</b></td><td>CCF Doublepeak</td><td></td><td></td><td></td></tr><tr><td>S1592</td><td>1.032</td><td>10.5</td><td>4</td><td>2</td><td>8</td><td>0</td><td>0</td><td>14</td><td>33.639</td><td>0.011</td><td>0.018</td><td>0.024</td></tr><tr><td>S1607</td><td>0.548</td><td>12.7</td><td>12</td><td>6</td><td>0</td><td>2</td><td>2</td><td>22</td><td>33.391</td><td>0.013</td><td>0.019</td><td>0.062</td></tr></tbody></table>
<figure><img src="image_11.png" /><figcaption>**Fig. A.1.** Contours of the planet detection efficiency for the MS stars in the domains of 1–1000 days for planet period and 0.2–10 *M*<sub>J</sub> for planet mass. A 0.01 FAP level has been used as detection planet threshold in the periodogram analysis. A filled circle indicates the planet position.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Fig. A.2.** Contours of the planet detection efficiency for the evolved stars in the domains of 1–1000 days for planet period and 0.2–10 *M*<sub>j</sub> for planet mass. A 0.01 FAP level has been used as detection planet threshold in the periodogram analysis. A filled circle indicates the planet position.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. A.3.** Top: Lomb-Scargle periodogram of the RV measurements, residuals, bisector span, and FWHM for YBP778. Bottom: same plots for S978. The dashed lines correspond to 5% and 1% FAPs, calculated according to Horne & Baliunas (1986) and white noise simulations.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. A.4.** Top: RV measurements versus bisector span, residuals versus bisector span, RV measurements versus CCF FWHM, and RV measurements versus $H\alpha$ activity indicator for YBP778. The $H\alpha$ activity indicator is computed as the area below the core of $H\alpha$ line with respect to the continuum. CCF FWHM values are calculated by subtracting in quadrature the respective instrumental FWHM. The same symbols as in Fig. 4. Bottom: same plots for S978.</figcaption></figure>