The HARPS-N Rocky Planet Search

The HARPS-N Rocky Planet Search

I. HD 219134 b: A transiting rocky planet in a multi-planet system at 6.5 pc from the Sun thanks: The photometric time series and radial velocities used in this work are available in electronic form at the CDS via anonymous ftp to ( or via
F. Motalebi 1 Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 1    S. Udry 1 Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 1    M. Gillon 2Institut d’Astrophysique et de Géophysique, Université de Liège, Allée du 6 Août 17, Bat. B5C, 4000 Liège, Belgium 2    C. Lovis 1 Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 1    D. Ségransan 1 Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 1    L. A. Buchhave 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 34Centre for Stars and Planet Formation, Natural History Museum of Denmark, University of Copenhagen, DK-1350 Copenhagen, Denmark 4    B. O. Demory 5Cavendish Laboratory, J J Thomson Avenue, Cambridge CB3 0HE, UK 5    L. Malavolta 6Dipartimento di Fisica e Astronomia ”Galileo Galilei” , Universita’ di Padova, Vicolo dell’Osservatorio 3, 35122 Padova, Italy 6    C. D. Dressing 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    D. Sasselov 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    K. Rice 7SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH93HJ, UK 7    D. Charbonneau 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    A. Collier Cameron 8SUPA, School of Physics & Astronomy, University of St. Andrews, North Haugh, St. Andrews Fife, KY16 9SS, UK 8    D. Latham 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    E. Molinari 9INAF - Fundación Galileo Galilei, Rambla José Ana Fernandez Pérez 7, 38712 Berña Baja, Spain 910INAF - IASF Milano, via Bassini 15, 20133, Milano, Italy 10    F. Pepe 1 Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 1    L. Affer 11INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy 11    A. S. Bonomo 12INAF - Osservatorio Astrofisico di Torino, via Osservatorio 20, 10025 Pino Torinese, Italy 12    R. Cosentino 9INAF - Fundación Galileo Galilei, Rambla José Ana Fernandez Pérez 7, 38712 Berña Baja, Spain 9    X. Dumusque 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    P. Figueira 13Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, PT4150-762 Porto, Portugal 13    A. F. M. Fiorenzano 9INAF - Fundación Galileo Galilei, Rambla José Ana Fernandez Pérez 7, 38712 Berña Baja, Spain 9    S. Gettel 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    A. Harutyunyan 9INAF - Fundación Galileo Galilei, Rambla José Ana Fernandez Pérez 7, 38712 Berña Baja, Spain 9    R. D. Haywood 8SUPA, School of Physics & Astronomy, University of St. Andrews, North Haugh, St. Andrews Fife, KY16 9SS, UK 8    J. Johnson 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    E. Lopez 7SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH93HJ, UK 7    M. Lopez-Morales 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    M. Mayor 1 Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 1    G. Micela 11INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy 11    A. Mortier 8SUPA, School of Physics & Astronomy, University of St. Andrews, North Haugh, St. Andrews Fife, KY16 9SS, UK 8    V. Nascimbeni 6Dipartimento di Fisica e Astronomia ”Galileo Galilei” , Universita’ di Padova, Vicolo dell’Osservatorio 3, 35122 Padova, Italy 6    D. Philips 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    G. Piotto 6Dipartimento di Fisica e Astronomia ”Galileo Galilei” , Universita’ di Padova, Vicolo dell’Osservatorio 3, 35122 Padova, Italy 6    D. Pollacco 14Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK 14    D. Queloz 1 Observatoire de Genève, Université de Genève, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 15Cavendish Laboratory, J J Thomson Avenue, Cambridge CB3 0HE, UK 5    A. Sozzetti 12INAF - Osservatorio Astrofisico di Torino, via Osservatorio 20, 10025 Pino Torinese, Italy 12    A. Vanderburg 3Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 3    C. A. Watson 15Astrophysics Research Centre, School of Mathematics and Physics, Queen’s University Belfast, Belfast, BT7 1NN, UK 15
Received on 24th of June 2015; Accepted on 30th of July 2015
Key Words.:
planetary systems: super-Earth – techniques: radial velocity – techniques: photometry – stars:individual: HD 219134 – binaries: eclipsing – instrument: HARPS-N

We know now from radial-velocity surveys and transit space missions that planets only a few times more massive than our Earth are frequent around solar-type stars. Fundamental questions about their formation history, physical properties, internal structure, and atmosphere composition are, however, still to be solved. We present here the detection of a system of four low-mass planets around the bright (V=5.5) and close-by (6.5 pc) star HD 219134. This is the first result of the Rocky Planet Search program with HARPS-N on the Telescopio Nazionale Galileo in La Palma. The inner planet orbits the star in  days, on a quasi-circular orbit with a semi-major axis of 0.0382  AU. Spitzer observations allowed us to detect the transit of the planet in front of the star making HD 219134 b the nearest known transiting planet to date. From the amplitude of the radial-velocity variation ( ms) and observed depth of the transit ( ppm), the planet mass and radius are estimated to be  M and  R, leading to a mean density of  g cm, suggesting a rocky composition. One additional planet with minimum mass of  M moves on a close-in, quasi-circular orbit with a period of 6.7650.005 days. The third planet in the system has a period of days and a minimum mass of  M, at 0.2340.002 AU from the star. Its eccentricity is . The period of this planet is close to the rotational period of the star estimated from variations of activity indicators ( days). The planetary origin of the signal is, however, the preferred solution as no indication of variation at the corresponding frequency is observed for activity-sensitive parameters. Finally, a fourth additional longer-period planet of mass of  M orbits the star in 1190 days, on an eccentric orbit () at a distance of 2.140.27 AU.

1 Introduction

The statistical analysis and occurrence rate of the small-mass planets in the range of super-Earth to Neptune are discussed in several studies with the main motivation of better understanding the orbital and physical characteristics of this population of planets, in order to constrain their formation processes (Mordasini et al. 2012; Benz et al. 2014). The radial-velocity planet search program with the HARPS spectrograph on the ESO 3.6-m telescope (Mayor et al. 2003, 2011), a similar survey with the Keck telescope (Howard et al. 2010), and the NASA Kepler transit space mission (Borucki et al. 2011) have in particular contributed in a tremendous way to our knowledge of the population of small-mass/size planets around solar-type stars.

The most common planets detected by the Kepler mission peak around 2 R (Howard et al. 2012; Fressin et al. 2013; Marcy et al. 2014). Even if such planets do not exist in our Solar System, they are found around more than 30 % of solar-type hosts. Furthermore, many of them are found in very coplanar multiple systems (Fabrycky et al. 2014; Howard et al. 2012; Figueira et al. 2012), tightly packed close to the central star, a new challenge to explain for planet formation models (Ogihara et al. 2015). Over the last 12 years the extra-solar planet zoo has also been supplied with low-mass planets detected by the HARPS GTO planet search program and successive subsequent ESO Large programs. Notable examples include HD 40307 (Mayor et al. 2009), HD 10180 (Lovis et al. 2011b), HD 20794 and HD 85512 (Pepe et al. 2011), or Alpha Cen B (Dumusque et al. 2012). Analysis of the results by Mayor et al. (2011) provided a list of additional super-Earth and mini-Neptunes unveiled by the survey, as well as a first statistical analysis of the properties and occurrence rate of the super-Earths and Neptune-mass planets around solar-type stars. These preliminary findings were confirmed by the fantastic statistics and precision of the Kepler detections (Fressin et al. 2013).

After a decade focusing mainly on the detection and the determination of occurrence rate of low-mass planets, a significant observational effort is now dedicated to the planet characterization. Transit results provide the planet radius and in combination with radial velocities, the mean density of the planet can be derived. Spectral features of exoplanet atmospheres may also be revealed by space and ground-based high-resolution transmission spectroscopy in the visible and near infrared. One of the main difficulties of such characterization is the availability of bright targets favourable for follow-up observations. In particular exquisite planetary physical parameter determination will be key in lifting the intrinsic degeneracy in the determination of the planet composition when several components (gas, silicates, metals) and chemical species are mixed in the planet interior. This was a strong driver for transit search follow-up of super-Earths detected with radial velocities around very bright stars, from space with the MOST, Hubble or Spitzer space telescopes (Gillon et al. 2012, 2010; Demory et al. 2015a), and for the development of a new generation of space transit missions to be launched within the coming decade: CHEOPS/ESA (Fortier et al. 2014), TESS/NASA (Ricker et al. 2014), and PLATO/ESA (Rauer et al. 2014).

In order to start to fulfill the need for good estimates of planet physical parameters and considering the successful achievements with HARPS, a similar design has been implemented for the HARPS-N spectrograph. The instrument was built to become an efficient exoplanet hunter in the northern hemisphere, specifically aimed at the follow-up and characterization of the brightest transit candidates from the Kepler mission. HARPS-N is a high-precision echelle spectrograph located at the 3.6 m Telescopio Nazionale Galileo (TNG) at the Roque de los Muchachos Observatory on La Palma, Spain. It was built by an international consortium of Swiss, Italian, UK and US partners, led by Geneva university. It began operations in August 2012. Eighty nights per year for five years were granted to the consortium for a science program including the follow-up of Kepler candidates at high precision, and an additional survey, the Rocky Planet Search (RPS), to search for small planets around bright and quiet nearby stars in the northern hemisphere when Kepler field is down.

While waiting for the new generation of ground- and space-based transit searches targeting bright stars, planet density estimates from the follow-up of the Kepler candidates with HARPS-N have already been presented in several papers (Pepe et al. 2013; Dumusque et al. 2014; Bonomo et al. 2014; Dressing et al. 2015). They are populating the low-radius regime of the mass-radius relation for small exoplanets with candidates for which precise mass and radius could be determined.

We present in this paper the first results from the Rocky Planet Search: a planetary system around HD 219134, composed of 3 inner super-Earths and an outer sub-Saturn mass planet, with the 4 planets spread between 0.038 and 2 AU from the central star. Thanks to the high precision photometric observations from the Spitzer space telescope, the inner planet is furthermore observed to transit in front of the star. After a short presentation of the Rocky Planet Search program in Sect. 2, and the stellar properties of HD 219134 in Sect. 3, the spectroscopic and photometric observations of the new system are described in Sect. 4 and 5. Sect. 6 provides a discussion of the stability of the system and of the composition of the inner planet. We finally draw some conclusions in Sect. 7.

2 The HARPS-N Rocky Planet Search program

The GTO (Guaranteed Time Observation) granted to the HARPS-N consortium is dedicated to two programs: i) the confirmation of the planetary nature and the characterization of Kepler candidates and ii) a Rocky Planet Search (RPS). The aim of the RPS program is to perform a systematic search for low-mass planets around nearby quiet stars visible from the northern hemisphere, through an intensive monitoring of the radial velocity of quiet stars at very high precision. HARPS-N is stabilized and well-controlled in pressure and temperature to minimize instrumental drifts, and so to secure sub-m/s radial velocities over long periods of time. More details on technical aspects are presented in Cosentino et al. (2012, 2014).

Target ID Distance [pc] V Spec.Type
HD38 11.750.38 8.20 K2
HIP1368 14.990.25 8.99 M0
HD3651 11.110.09 5.88 K0V
HD4628 7.460.05 5.74 K2V
HD10476 7.470.05 5.24 K1V
HD10436 13.430.20 8.42 K5V
HD16160 7.210.05 5.79 K3V
HD19305 14.770.39 9.07 M0V
HD232979 10.190.12 8.62 K8V
HD32147 8.810.06 6.22 K3V
HIP27188 13.690.21 9.02 M0
HD41593 15.450.16 6.76 K0
HD47752 17.520.35 8.08 K2
HD48948 16.400.38 8.59 M0
HIP36357 17.550.29 7.73 K2V
HD62613 17.040.11 6.55 G8V
HD65277 17.460.39 8.05 K5V
HD65583 16.800.16 6.97 G8V
HIP39826 17.120.54 9.41 M0
HIP42220 13.900.24 9.28 M2
HIP43534 16.490.37 9.26 K5
HD79210 6.190.20 7.64 M0V
HD79211 6.270.26 7.70 K2
HD79969 17.530.23 7.20 K3V
HD84035 17.780.31 8.13 K5V
HD88230 4.870.20 6.60 K8V
HIP51525 15.670.27 8.85 K7V
HIP51547 17.490.42 9.63 M0
HD94765 17.550.26 7.37 K0
HD97101 11.930.15 8.31 K8V
HD97503 17.950.41 8.70 K5V
HD99492 17.990.47 7.58 K2V
HD103095 9.160.07 6.42 G8V
HD110315 14.300.22 7.91 K2
HD111631 10.780.11 8.49 M0.5V
HD122064 10.100.06 6.49 K3V
HD128165 13.420.12 7.24 K3V
HD144579 14.370.12 6.66 G8V
HD147379 10.660.14 8.61 M0V
HD151288 9.770.08 8.10 K7V
HD158633 12.800.08 6.44 K0V
HD157881 7.720.06 7.54 K7V
HD166620 11.100.07 6.38 K2V
HD173818 14.120.25 8.81 K5
HD185144 5.770.02 4.67 K0V
HD184489 14.470.35 9.35 K5
HD190007 13.110.18 7.46 K4V
HD200779 14.840.30 8.27 K5
HD201091 3.480.02 5.20 K5V
HD201092 3.500.01 6.05 K7V
HD219134 6.530.03 5.57 K3V
Table 1: List of the RPS targets in the GTO program on HARPS-N drawn from the Hipparcos catalogue and van Leeuwen (2007).

The sample: The first step in establishing the RPS program was to select a sample of stars best suited for long-term precise radial-velocity measurements. Uncertainties in such observations are mainly linked to noise from three different origins: photon noise, instrumental noise, and stellar intrinsic signals. The first important criterion for this program was thus to concentrate on bright stars, in our closest neighbourhood. This is also in the interest of potential follow-up studies for characterization of the planet properties. The sample was thus selected from the Hipparcos catalogue within a distance limit set to 18 pc111The limit was chosen iteratively in order to have a sample large enough and covering a full range of right ascension, with an over density during winter when the Kepler field in not visible. The second step was to focus on quiet stars. Based on the CORAVEL data and investigation in the literature, we rejected known spectroscopic binaries as well as stars with  kms and stars with from our sample. In addition, K dwarfs are favoured in our sample compared to G dwarfs because of their lower level of stellar ”noise” (p-mode, granulation and activity jitter; Dumusque et al. 2011b). Their habitable zone is also closer in. Adding all the criteria together, we ended up with a sample of 51 quiet stars with a range of spectral types from G8 to M0. They are listed in Table 1.

Observational strategy: To minimize the effect of stellar noise with short typical time scales (p-modes and granulations), we applied the observational strategy implemented for the HARPS very high-precision observations (Pepe et al. 2011). With this strategy, we observe each target with 15 minutes exposure time to damp the p-modes effect and make a 2nd (and possibly a 3rd) measurement of the target well-spread in time during the night to damp the effect of granulation (for more details see Dumusque et al. 2011b). For bright stars, in order to avoid saturation the 15 minutes on target are split into sub-observations. Typical average signal-to-noise ratio (SNR) of the spectra obtained are between 200 and 400, measured at  nm.

With continuous monitoring during the past two years of operation of HARPS-N we have acquired hundreds of precise radial velocities of the stars in the RPS sample. The typical photon-noise precision per [sub-]observation is about 0.3 ms. In Fig. 1, we show the radial-velocity rms for all the stars in the RPS program. This rms naturally includes photon noise, instrumental effects (telescope, spectrograph, detector), stellar intrinsic ”noise”, and of course signatures of still undetected planets. With a mode around 1.48 ms, the distribution is very similar to the one obtained for the sample of solar-type stars followed at high precision with HARPS in the southern hemisphere.

Figure 1: Histogram of the dispersion of the radial-velocity measurements of stars in the RPS program on HARPS-N.

We present here the first result of the HARPS-N RPS program: the discovery of a planetary system with a transiting planet around HD 219134.

3 Stellar characteristics of HD 219134

The very bright star HD 219134 has been extensively studied in the literature. Basic photometric properties of this star were taken from the Hipparcos catalogue (van Leeuwen 2007). A precise estimate of the radius of the star,  R, is available from interferometric measurements (Boyajian et al. 2012). Combined with a luminosity of  L, the Stefan-Boltzmann law gives  K. We will adopt this value for our photometric analysis of the transiting planet in the system in Sect. 5.

Technique [K]
Photometry 0.26
EWs 0.35

Valenti & Fischer (2005); Mishenina et al. (2008); Ramírez et al. (2013); Prugniel et al. (2011); This paper

Table 2: Comparative summary table for the atmospheric parameters derived for HD 219134.

Over the past few years, the star has been the object of several studies aiming at determining photospheric parameters and chemical abundance analyses (Valenti & Fischer 2005; Mishenina et al. 2008, 2012; Kovtyukh et al. 2003; Ramírez et al. 2013; Prugniel et al. 2011). Results of these studies are reported in Table 2 for comparison.

We also derived the atmospheric stellar parameters directly from HARPS-N spectra using two approaches: one based on Equivalent Width (EW) determination and one using the Stellar Parameter Classification (SPC) tool.

For the Equivalent Width (EW) approach, we followed the procedure described in Sect. 3 of Dumusque et al. (2014). We used the 2014 version of the local thermodynamic equilibrium code MOOG (Sneden 1973) together with the Kurucz grid of atmosphere models (Castelli & Kurucz 2004; Kurucz 1992), while oscillator strength values in the line list from Sousa et al. (2011) were updated in accordance with the solar iron abundance from Asplund et al. (2009). In order to construct high-qualitity spectra, spectra within the same visit (sub-observations) were co-added, obtaining at the time of the analysis 88 spectra with SNRs between 400 and 700. EWs were then measured for each of them and their mean and rms were used for the atmospheric parameter determination. We only retained the lines within the range and with a dispersion lower than either 1 m or 5 % of the mean EW. Despite the quality of our spectra, the preliminary analysis resulted in a poor determination of the microtubulent velocity . Following the calibration from Ramírez et al. (2013) and Tsantaki et al. (2013), we decided to fix the microturbulent velocity to . The resulting atmospheric parameters are reported in Table  2, with the gravity from FeII lines already increased by according to the calibration in Sect. 3.1 of Mortier et al. (2014). Note that the derived parameters are dominated by systematic errors (e.g. choice of the oscillator strength, continuum placement) rather than random errors on EW measurements.

Finally, we used the SPC tool, an alternative approach to derive atmospheric parameters by cross-correlating an observed spectrum with a library of synthetic spectra (Buchhave et al. 2012, 2014). With SPC we obtain an effective temperature K, a surface gravity of and kms. The metallicity derived by SPC is , from a mix of metallic absorption lines in the wavelength range between 5050 to 5360 . The derived values are reported in Table 2 as well.

The mass of HD 219134, M M, was estimated through the Synthetic Clusters Isochrones & Stellar Tracks tool (SYCLIST)222 SYCLIST allows the user to determine stellar parameters through a Bayesian-based interpolation of the grid of Geneva stellar evolution models. For HD 219134, we used the taken from Boyajian et al. (2012), the apparent magnitude and parallax from the Hipparcos catalogue (van Leeuwen 2007), and the metallicity derived from our spectral analysis as input parameters. For the latter we used the average of the values obtained from the EWs and SPC approaches.

The mean activity index log R is estimated from the HARPS-N spectra. We derive an average value of with a dispersion of 0.06 dex. To estimate the rotational period of the star we used a periodogram analysis of the activity indicators (, CCF FWHM and CCF bisector span time series) that yields a most significant peak at 42.3 days (see Sect. 4.2.3).

Table 3 summarises the stellar parameters obtained from the various catalogues and analyses mentioned above, with the final values selected for the star when several estimates are available.

Parameter Ref.
V 5.57 1
B-V 0.99 1
SpTyp K3V 1
[mas] 1
L [L] 2
R [R] 2
[K] 2
[dex] this paper
[kms] this paper
[km s] this paper
[dex] this paper
M [M] this paper, 2
¡¿ this paper
P [days] 42.3 0.1 this paper
Table 3: Stellar parameters taken from the Hipparcos catalogue (van Leeuwen 2007), interferometric measurement (Boyajian et al. 2012) and estimated from HARPS-N spectra adopting EWs and SPC and stellar evolution models SYCLIST for M (see text).

4 Spectro-velocimetric observations

4.1 The HARPS-N data

HD 219134 has been monitored by the HARPS-N spectrograph for close to three years, from BJD= 2456148.7 (9th of Aug. 2012) to BJD = 2457195.7 (22nd of June 2015). To mitigate the effects of stellar oscillations, the strategy was to stay on target for 10 to 15 minutes. Because of the brightness of the star, the observations were split in several individual sub-observations to avoid saturation. We obtained a total of 481 data points spread over 99 epochs (nightly averaged values). The nightly averaged values are displayed in Fig. 2. Observations were performed using the simultaneous thorium calibration technique. The mean uncertainty on the individual RVs due to photon noise and known calibration noise is 0.4 ms. This corresponds to an average SNR of 268 per pixel at = 550 nm. The data reduction was carried out with the latest version of the HARPS-N pipeline (DRS 3.7) using the cross-correlation technique (Pepe et al. 2002; Baranne et al. 1996). On top of precise radial velocities, the pipeline directly provides parameters estimated from the cross-correlation function (CCF) of the spectrum: CCF full width at half maximum (FWHM), CCF contrast (minimum vs continuum), CCF bisector span inverse slope, and Ca II activity index S and log R (Fig. 2).

The raw rms dispersion of the radial velocities is 3.57 ms. Once de-trended from an obvious drift, the observed radial velocities still show a dispersion of 2.77 ms, which is significantly above the typical dispersion of such quiet stars, calling for a search for additional coherent signals in the data.

Figure 2: Radial-velocity (top) and activity index (bottom) time series for HD 219134.

4.2 Data analysis

4.2.1 Data Modelling

The first step of the radial-velocity data analysis consists in identifying significant periodic signals in the data. This was done using the General Lomb-Scargle periodogram algorithm (GLS, Zechmeister & Kürster 2009) applied to the nightly averaged radial-velocity measurements to which a systematic error of 1 ms was quadratically added. False alarm probabilities were estimated through a bootstrap approach by permuting the nightly averaged data. Once a significant peak was located at a given period, the corresponding Keplerian was adjusted and removed. The process was repeated several times until no significant peak remained. For multiple Keplerians, all parameters were re-adjusted at each step of the analysis. During the analysis, we identified one radial velocity outlier (5  at jdb=56829.7) and decided to remove it to allow for a robust frequency analysis.

Figure 3: Generalized Lomb Scargle periodogram of the radial ve- locities of HD 219134, after removing - from top to bottom - a long period Keplerian, and successively the 3.09 d signal, the 46.78 d signal and finally the 6.75 day Keplerian. Dashed lines indicate 0.1%, 1%, and 10% level of false alarm probabilities.

Following this procedure, a first low frequency oscillation at 1190 day was identified with a semi-amplitude of  4.5 ms and a 0.8 % false alarm probability (FAP). After removing the corresponding Keplerians, two highly significant peaks were seen in the periodogram, at 3.09 and 46.78 days, with semi-amplitudes of 2.33 and 1.94 ms, and FAPs smaller than 0.01 % and 1 %, respectively (Fig. 3). One additional peak remains at 6.76 d with a significant FAP smaller than 1 % and a semi-amplitude of 1.1 ms. No more periodic signals remain in the data after subtraction of the corresponding 4 Keplerian model. This purely frequentist approach, using the nightly average data, has the advantage of being ”simple and fast” and is also conservative in terms of detection limits.

To double check the purely frequentist approach, we also analysed models of different complexity ranging from 1 to 4 Keplerian. The best model was chosen by comparing the Bayesian Information Criterion (BIC, Kass & Raftery 1995) between the different models:

where is the maximum of the log likelihood, the number of free parameters in the model, and is the number of measurements. In model comparisons based on the Bayes factor, when the prior information is smaller than the information provided by the data, the BIC (a.k.a Schwarz criterion) indicates that the model with the highest probability is the one with the minimum BIC. A summary of the comparison between a few models is provided in Table 4. The 4 Keplerian model has the lowest BIC and differs from the 3 Keplerian model by BIC=34. It is usually assumed that a BIC of 20 between two models is considered as strong evidence in favour of the most complex one which leads us to adopt the 4 Keplerians as our best solution.

Model (Period [days]) BIC BIC
K1 (1190) 611 6.35 2.63
K2 (1190, 3.09) 389 221 3.90 2.00
K3 (1190, 3.09, 46.8) 275 114 2.47 1.54
K4 (1190, 3.09, 46.8, 6.76) 241 34 1.89 1.31
Table 4: Models tested on the 98 nightly averaged radial velocities with an additional noise of 1 ms quadratically added to the radial-velocity error bars. The 4 Keplerian model is preferred since it has the lowest BIC, and residual dispersion. These values are estimated based on the maximum likelihood.

To obtain robust confidence intervals for the Keplerian parameters as well as an estimate of the additional noise present in the data (nuisance parameter below), we further probe the parameter space with a Markov Chain Monte Carlo algorithm (MCMC) with Metropolis-Hasting. An abundant literature discusses in much detail the implementation of MCMC posterior sampling (Andrieu & Thoms (2008) for a pure statistical approach; Gregory (2005a, b), Collier Cameron et al. (2007), and Pollacco et al. (2008) for exoplanet searches). Our MCMC probes the following set of parameters: , , , (the mean longitude at a given epoch) while the noise model follows a simple normal law with standard deviation derived from the observation errors and a nuisance parameter (). Jeffrey’s priors are used for the period, the radial-velocity semi-amplitude, and the nuisance parameter while uniform priors are used for the other parameters.

Figure 4: Phase-folded radial-velocity measurements of HD 219134 with the corresponding Keplerian model (solid line) for each of the 3 inner super-Earths, after removing the contribution of all the other planets in the system. From top to bottom, we have the 3.09-, 6.76-, and 46.78-day periods.

4.2.2 Orbital solution

The orbital elements corresponding to the four Keplerian model and transit timing constraint are listed in table 5 while the phase folded radial velocities are displayed in figures 4 and 5. The shortest period signal, with a radial velocity semi-amplitude of 2.33 ms at 3.09 days, is clearly identified as a super-Earth (4.46 M) in a circular orbit. One additional low amplitude signal (K=1.1 ms) is present at 6.76 days, corresponding to a planet with minimum mass of 2.67 M. The existence of this signal is corroborated by both the FAP and the BIC estimators but its strength might be affected by the un-perfect modelling of the long-period signal (see below).

The third signal with a period of 46.78 days and a semi-amplitude of 1.94 ms corresponds to a super-Earth of 8.67  M with an eccentric orbit of . We show, in the following section, that this period, which is close to the rotational period of the star (42.3 days), is not of stellar origin. Its relatively high eccentricity may also originate from an un-perfect modelling of the outer signal.

The long period radial velocity trend clearly shows two extrema, which favour a Keplerian instead of a polynomial drift, with a period of 1190 days and a semi-amplitude of 4.5 ms. As explained below, its origin cannot be linked to any long term stellar activity fluctuation and is undoubtedly of planetary origin. Due to the incomplete coverage of the orbit, both the eccentricity and the period of the planet remain moderately constrained but still leads to a good estimate of its minimum-mass, i.e. 62 M.

Model K4+
HD 219134 b HD 219134 c HD 219134 d HD 219134 e

undefined undefined
  kms   ms
s [m/s]
Table 5: Orbital solution for the 4 Keplerian model (K4) and planet inferred parameters for the system around HD 219134. is the expected date of transit. is the mean longitude at the time 2457126.7001 day, corresponding to the transit timing (see Sect. 5)
Figure 5: Radial-velocity measurements as a function of time of HD 219134 with the corresponding Keplerian model (solid line) for the outer planet in the system, after removing the contribution of the 3 inner ones.
Figure 6: Plot of the star activity index as a function of the radial-velocity residuals after removing the 3 shorter-period planet contributions for HD 219134. No clear positive correlation is observed, contrary to what is expected in the case of magnetic cycles (Lovis et al. 2011a).

4.2.3 Periodic signals: planet vs activity-related origin

Long-period magnetic cycle: The radial-velocity and activity index measurements of the star HD 219134 are displayed in Fig. 2. Although both exhibit long-term variations, the time series is very stable over the first 750 days of the survey and slightly increases over the last 300 days whereas the radial velocity constantly decrease over the  1050 day observation time span. No correlation is observed between RVs and log R as shown on Fig. 6. The long term trend of the is most probably related to the magnetic cycle of the star (Lovis et al. 2011a; Dumusque et al. 2011a; Meunier & Lagrange 2013) while the observed 1190 day period corresponds to a long period companion. Moreover, at the observed level of activity, we expect a low impact of stellar activity on the observed radial velocities.

Rotation period: Although different, the period of the third planet (46.8 days) is not very far from the rotational period of the star estimated from coherent variations of the activity indicators. A periodic variation around 42.3 days is indeed observed in the log R as well as in the CCF FWHM and in the bisector span of the radial velocities, whereas no signal appears around 46.8 days, as illustrated in Fig. 7 showing the GLS periodograms of these parameters zoomed around the periods of interest. We can therefore safely consider the 42.3 days as a valid estimation of the rotational period of the star (P). Conversely, no signal at 42.3 days appears in the radial-velocity data, supporting the planetary solution.

Figure 7: From Top to Bottom the panels show the GLS periodograms of radial velocities, , CCF bisector span, and CCF FWHM of the 2-year data set, zoomed around the periods of interest. The thick red line in each panel, indicates the planetary period (46.8 days) and the thick green line indicates the stellar rotation period (42.3 days). The thinner lines are the corresponding yearly aliases.

In order to avoid any misinterpretation of the stellar activity as a planetary signal, we pushed our investigations a bit further. We first examined the yearly aliases (Dawson & Fabrycky 2010) of the 46.8 and 42.3 days signals (Fig. 7) and confirmed that one period is not an alias of the other. In a second step, we considered several subsets of the data to check the persistence of the planetary signal over time and to mitigate the possible effects of discontinuities in the data sampling. The same features are observed. This confirms that the planetary signal at 46.8 days is present in the radial velocities at any time, and that in the same data no signal appears at the estimated (42.3 days).

A final argument in favour of the planetary interpretation of the 46.8 day signal is provided in Fig. 8 by the absence of a correlation between the activity index and the residuals around a 3-planet solution (leaving the 46.8 day period out). Such a correlation would be expected if the radial-velocity variation is induced by activity-related spots or plages on the star surface.

Figure 8: Radial velocities after subtracting the contributions of the planets at 3.09, 6.76 and 1190 days, displayed as a function of the activity index detrended with a polynomial of degree 3. No correlation is observed supporting so the planetary interpretation of the 46.8-day signal still present in the data.

5 Transit detection with Spitzer space telescope

5.1 Spitzer observations

Thanks to its short-period of 3.09 d, HD 219134 b has an interestingly high geometric transit probability of 9%. In Feb 2015, we requested 9.5hr of Director’s Discretionary Time (DDT) on the Spitzer Space Telescope to continuously monitor the 2- transit window of the planet as derived from our analysis of the HARPS-N RVs. As demonstrated by its co-detection and subsequent studies of the transits and occultations of 55 Cnc e (Demory et al. 2011, 2012, 2015b; Gillon et al. 2012) and its confirmation of the transiting nature of HD 97658 b (Van Grootel et al. 2014), Spitzer is indeed a very powerful facility to search for and measure with high-precision the transits of low-mass RV planets, thanks to its trailing orbit that allows monitoring the same star continuously for a complete transit window, and to its ultra-high photometric precision (a few dozens of ppm per time interval of 30 min for 55 Cnc). Our program was approved by the Spitzer Science Center (SSC) Director under the ID 11180 (PI: M. Gillon), and the SSC managed to schedule it for 2015 Apr 14th, a few days before the end of the current visibility window of the star.

We observed HD 219134 at 4.5 m with the Spitzer/IRAC detector (Fazio et al. 2004) in subarray mode (32x32 pixels windowing of the detector), the extremely fast Fowler sampling (0.01s) maximizing the duty cycle and SNR. No dithering pattern was applied to the telescope (continuous staring). HD 219434 being an extremely bright star for Spitzer, we used the shortest available integration time of 0.01s, allowing the counts to remain in the linear regime of the detector. The observations were performed from 2015 Apr 14th 01h58 to 11h18 UT. We used the recently introduced PCRS peak-up mode (Grillmair et al. 2012; Ingalls et al. 2014) which was previously used by Ballard et al. (2014) to estimate the infrared transit depth of Kepler-93 b. This mode provides enhanced accuracy in the position of the target on the detector, leading to a significant decrease of the so-called ‘pixel phase effect’ that is the most important source of correlated noise in high-SNR staring mode observation with IRAC InSb arrays (e.g. Knutson et al. 2008). The run consisted of a 9hr-long science Astronomical Observational Requests (AOR) preceded by a short (30 min) AOR to allow the spacecraft to stabilize.

5.2 Data reduction

After download to Earth and basic calibration with the Spitzer pipeline S19.1.0, the images were made available to us by SSC through the Spitzer Heritage Archive (SHA) web interface333 under the form of Basic Calibrated Data (BCD). Each subarray mode BCD is composed of a cube of 64 subarray images of 3232 pixels (pixel scale = 1.2 arc second). We used the following strategy to reduce these BCDs. We first converted fluxes from the Spitzer units of specific intensity (MJy/sr) to photon counts, then aperture photometry was performed on each subarray image with the IRAF/DAOPHOT444IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. software (Stetson 1987). We tested different aperture radii, and selected 2.3 pixels as the radius minimizing at best the white and red noises in the residuals of a short data fitting analysis. The centre and width of the Point-Spread Functions (PSF) were measured by fitting a 2D-Gaussian profile on each image. The distribution of the measurements was then looked at, and measurements having a visually discrepant position relative to the bulk of the data were then discarded. For each block of 64 subarray images, we then discarded the discrepant values for the measurements of flux, background, - and -positions using a 10- median clipping for the four parameters, and the resulting values were averaged, the photometric errors being taken as the errors on the average flux measurements. Finally, a 50- slipping median clipping was used on the resulting light curves to discard outliers (due, e.g., to cosmic hits).

Our resulting light curve counted 9396 measurements. It is shown in Fig. 9, with the evolution of relevant external parameters (PSF - and -center and PSF width, background). Its time sampling being much shorter than the structures of the expected transit and of the typical timescale of the Spitzer systematics, we binned the light curve to time intervals of 30s for the sake of computational speed of the data analysis. Nevertheless, we verified with a shorter version of the data analysis procedure described below that our results are insensitive to the binning of the photometry.

5.3 Data analysis

We analysed the Spitzer photometric time-series with our adaptative MCMC code (see Gillon et al. 2012 and references therein). The assumed photometric model consisted of the eclipse model of Mandel & Agol (2002) to represent the possible transit of HD 219134 b, multiplied by a baseline model aiming to represent the other astrophysical and instrumental effects at the source of photometric variations. We assumed a quadratic limb-darkening law for the star. We based the selection of the baseline model on the minimization of the Bayesian Information Criterion (BIC, Schwarz 1978).

Figure 9: Evolution of the following measured parameters in the Spitzer images for HD 219134: PSF centre - and -positions (), PSF - and -width (), background counts (bottom left), and stellar fluxes (bottom right). The times of the start and end of the detected transit are shown as red vertical lines.

Following Gillon et al. (2014), the instrumental models included three types of low-order polynomials. The first one had as variables the - and -positions of the centre of the PSF to represent the ‘pixel phase’ effect of the IRAC InSb arrays (e.g. Knutson et al. 2008). The second one had as variables the PSF widths in the - and/or the -direction, its inclusion in the baseline model strongly increasing the quality of the fit for Warm Spitzer photometry (see also Lanotte et al. 2014). The third, optional, function was a polynomial of the logarithm of time + a slope to represent a sharp decrease of the detector response at the start of the AOR (‘ramp’ effect, Knutson et al. 2008). To improve the quality of the modelling of the pixel phase effect, especially the fitting of its highest frequency components, we supplemented the - and -polynomial with the Bi-Linearly-Interpolated Sub-pixel Sensitivity (BLISS) mapping method (Stevenson et al. 2012). The sampling of the positions space was selected so that at least five measurements fall within the same sub-pixel. We refer the reader to Gillon et al. (2014) for more details.

Figure 10: Spitzer photometry divided by two different best-fit baseline models. The first one () aims to remove position- and PSF effects, and consists in a 4th order polynomial of the PSF -and centres and widths, supplemented with the BLISS mapping method (Stevenson et al. 2012) to remove high-frequency position effects. The second baseline model () adds to the first one a slope, and a quadratic function of the logarithm of time to model the sharp decrease of the counts at the beginning of the run (negative ramp). Both light curves are binned per 0.005d = 7.2min intervals.
Figure 11: Spitzer photometry divided by the best-fit baseline model and binned per 0.005d = 7.2min, with the best-fit transit model over-imposed in red. Below are shown the -shifted residuals of the fit binned per 7.2 min and 30 min intervals. Their standard deviations are, respectively, 57 ppm and 25 ppm. Over the light curve is shown the prior probability distribution function derived for the transit timing of HD 219134 b from our analysis of the HARPS-N RVs.

Assuming no transit of HD 219134 b, the light curve corrected for the instrumental model described above showed a transit-like structure of ppm depth and lasting 50 min (see Fig. 10), this structure remaining if higher order terms are used in the polynomial functions. Our first reflex was to check that this transit-like structure did not correspond to any odd behaviour of the external parameters, which was not the case (Fig. 9). Before identifying the structure to the searched transit, we nevertheless performed a large set of short MCMC analyses assuming different baseline models, and assuming or not that the structure was a transit of the planet. For each baseline model, we computed the BIC difference between the best-fit models with and without transit to estimate the Bayes factor in favour of the transit hypothesis, multiplying the likelihood term by 9/91, the prior transit/no-transit probability ratio. In these tests, we multiplied the photometric errors by a Correction Factor (CF, see below) of 2.72, the highest value that we ever encountered in our past experiences of high-precision photometry with Warm Spitzer. It corresponds to an AOR targeting 55 Cnc (Demory et al., in prep.) for which the PCRS peak-up mode did not work properly. Doing so, we thus assumed that the HD 219134 light curve was strongly affected by some correlated noise poorly reproduced by our instrumental model. At the end, the lowest value that we obtained for the Bayes factor under this extreme hypothesis was 1000 in favour of the transit hypothesis. We could thus conclude with certainty the transiting nature of the planet.

We then performed a longer MCMC analysis to probe the posterior probability distribution of the transit parameters. The orbit of HD 219134 b was assumed to be circular in this MCMC analysis. The jump parameters of the MCMC, i.e. the parameters randomly perturbed at each step of the Markov Chains, were the following.

  • The stellar mass , radius , effective temperature , and metallicity [Fe/H]. For these four parameters, normal prior probability distribution functions based on the values given in Table 3 were assumed.

  • The planet/star area ratio .

  • The impact parameter of the planet, where is the orbital semi-major axis and is the orbital inclination. A full-transit configuration corresponds to .

  • The time of inferior conjunction for the planet, corresponding to the mid-time of the transit.

The limb-darkening of the star was modeled by a quadratic law (Claret 2000). Values for the two quadratic limb-darkening coefficients and were drawn at all steps of the MCMC from normal distributions with expectations and standard deviations drawn from the tables of Claret & Bloemen (2011) for the Spitzer 4.5 m bandpass and for the stellar atmospheric parameters drawn at the same step.

Five chains of 100,000 steps were performed for each analysis, their convergence being checked using the statistical test of Gelman & Rubin (1992). They followed a preliminary chain of 100,000 steps performed to estimate the need to rescale the photometric errors, at the end of which the standard deviation of the residuals was compared to the mean photometric errors, and the resulting ratios were stored. represents the under- or overestimation of the white noise of each measurement. On its side, the red noise present in the light curve (i.e. the inability of our model to represent perfectly the data) was taken into account as described in Gillon et al. (2010), i.e. a scaling factor was determined from the standard deviations of the binned and unbind residuals for different binning intervals ranging from 5 to 120 minutes, the largest values being kept as . At the end, the error bars were multiplied by the correction factor . The derived values for and were 1.30 and 1.01, resulting in .

Table 6 presents the resulting values plus error bars for the transit and planet’s parameters, while Fig 11 shows the light-curve corrected for the systematics, the best-fit transit model, and the residuals.

Transit parameters Value
Depth ppm
Impact parameter
Timing BJD
Duration min
Physical parameters
Planet radius
Orbital inclination deg
Table 6: Transit and physical parameters deduced from the MCMC analysis of the Spitzer photometry.

6 Discussion

6.1 Dynamical stability

A very important and necessary a posteriori consistency check of a planetary system characterization is needed to verify that the dynamical evolution of the system is viable on the long term, assuring the persistence of the system from the end stage of its formation (when the protoplanetary disk disappears) till its observation today. Pure n-body integrations of the 4-planet system, using both the GENGA Code (Grimm & Stadel 2014) launched through the DACE interface 555DACE is a platform of the National Centre for Competence in Research ‘PlanetS’, that can be accessed at, and a fourth-order Hermite scheme (Makino 1991), with initial conditions derived from parameters in Table 5, and assuming coplanarity and an inclination of 5 degrees from the transit observation, shows that the system is stable for more than orbits of the outermost planet.

A more complete analysis, taking into account general relativity and tides as well as longer-term secular effects, is beyond the scope of this paper. As expected, a preliminary check indicates, however, that general relativity and tides will act to limit the secular growth of the eccentricities of the inner planets, favouring the long-term stability of the system.

6.2 Bulk composition of the planet HD  219134 b

HD  219134 b is the nearest transiting super-Earth known today. The radius and mass of the planet have been characterized to 6% and 9% accuracy, allowing us to place tight constraints on the bulk density. It thus will become one of the best targets for internal structure and atmosphere characterization with the Spitzer and HST and the up-coming space follow-up missions e.g. TESS, JWST, CHEOPS and PLATO. Its position in a stellar distance vs planetary mass diagram is provided in Fig. 12, in comparison with the other most favourable known cases for characterization within 40 pc from the sun and up to 35 M: GJ 436 b (Butler et al. 2004; Gillon et al. 2007), 55 Cnc e (Gillon et al. 2012; de Mooij et al. 2014), GJ 1214 b (Charbonneau et al. 2009), HD 97658 b (Howard et al. 2011; Dragomir et al. 2012; Van Grootel et al. 2014), GJ 3470 b (Bonfils et al. 2012; Demory et al. 2013) and HAT-P-11 b (Dittmann et al. 2009).

Figure 12: Planet mass vs distance to the host stars for planets in close neighborhood. Transiting planets with good mass and radius (density) determinations are shown as in this diagram. The planets detected by radial velocity are also shown using the minimum-mass as a proxy for the mass (data taken from
Figure 13: Mass-radius relation for planets with radii smaller than 2.7 R and with masses determined to a precision better than 20 % (updated from Dressing et al. (2015)). The shaded gray region in the lower right indicates planets with iron content exceeding the maximum value predicted from models of collisional stripping (Marcus et al. 2010) . The solid lines are theoretical mass-radius curves (Zeng & Sasselov 2013) for planets with compositions of 100 % HO (blue), 25 % MgSiO - 75 % HO (purple), 50 % MgSiO - 50 % HO (green), 100 % MgSiO (black), 50 % MgSiO - 50 % Fe (red), and 100 % Fe (orange). In this diagram, the position of HD 219134 b is almost overlapping the point for CoRoT-7 b. It belongs to a group of planets including Kepler-36 b, Kepler-93 b, and Kepler-10 b.

Combining the spectroscopic and space-based photometric data, we estimate a density of gcm for the planet. This density is consistent with the value of 6.90 gcm that would be predicted for a  R planet obeying the Earth-like compositional model presented in Dressing et al. (2015).

The compositional tracks employed in Dressing et al. (2015) (reproduced and updated in Fig. 13) are based on interior structure models by Zeng & Sasselov (2013) that represent small planets as fully differentiated iron cores surrounded by lower density magnesium silicate mantles. These models provide a useful framework for comparing relative planet compositions, but the absolute core mass fractions are underestimated slightly because the Zeng & Sasselov (2013) models do not incorporate the presence of lighter elements in the core and the inclusion of water in the mantle. Accordingly, the Earth-like compositional track presented in Dressing et al. (2015) corresponds to a model composition of 83% MgSiO and 17% Fe whereas the actual core mass fraction of the Earth is closer to 30%. Zeng et al. (2015) have recently updated their model framework to account for the presence of lighter elements in the core and the presence of water in the mantle. Employing the new models, we find that the population of highly irradiated dense planets (HD 219134b, CoRoT-7b, Kepler-10b, Kepler-36b, Kepler-78b, and Kepler-93b) are best described by a two-component iron-magnesium silicate model with a core mass fraction of approximately 22-23%.

7 Conclusion

We have presented in this paper the first result from the Rocky Planet Search (RPS) program conducted with HARPS-N, as a planetary system of three inner super-Earths and one outer sub-Saturn planet, hosted by the bright and nearby quiet K dwarf HD 219134. The planet separations between 0.039 and 0.23 AU called for a search of potential transits of the inner planet(s) with the Spitzer space telescope. The successful detection of the transit of HD 219134 b makes the star the closest (6.5 pc) and brightest (V=5.5) star known to date with a transiting planet (super-Earth). This system is thus becoming one of the most favourable ones for follow-up observations, for better constraining the system architecture or aiming at the characterization of the planet physical properties. First, it provides an exquisite constraining point in the mass-radius diagram (Fig. 13) for the composition of the planet, found to be of terrestrial-equivalent composition with a core mass fraction of the order of 22 %. The quality of the measurements of the radius, mass and then mean density actually foreshadows what can be expected from the future transit missions in preparation that will target bright stars (CHEOPS, TESS, PLATO). We also know from Kepler results that multi-transiting systems of small-size planets are numerous. It is then now highly suitable to search for traces of transits of the other planets in the systems. Finally, even if a potential atmosphere around the planet is expected a priori to be tiny, the brightness of the system makes it worth trying to detect features of this atmosphere in the UV, visible and NIR, from space and from the ground, especially in preparation for future measurements with larger facilities (JWST, TMT).

The HARPS-N project was funded by the Prodex program of the Swiss Space Office (SSO), the Harvard University Origin of Life Initiative (HUOLI), the Scottish Universities Physics Alliance (SUPA), the University of Geneva, the Smithsonian Astrophysical Observatory (SAO), and the Italian National Astrophysical Institute (INAF), University of St. Andrews, Queen’s University Belfast, and University of Edinburgh. The research leading to these results has received funding from the European Union Seventh Framework program (FP7/2007- 2013) under grant agreement No. 313014 (ETAEARTH). C.D. is supported by a National Science Foundation Graduate Research Fellowship. P.F. acknowledges support by Fundaçao para a Ciência e a Tecnologia (FCT) through Investigador FCT contracts of reference IF/01037/2013 and POPH/FSE (EC) by FEDER funding through the program ”Programa Operacional de Factores de Competitividade - COMPETE”. This work has been carried out in the frame of the National Centre for Competence in Research ‘PlanetS’ supported by the Swiss National Science Foundation (SNSF). S.U., C.L., D.S. and F.P. acknowledge the financial support of the SNSF. This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was provided by NASA. M. Gillon is Research Associate at the Belgian Scientific Research Fund (F.R.S-FNRS), and he is extremely grateful to NASA and SSC Director for having supported his searches for RV planets with Spitzer. PF further acknowledges support from Fundação para a Ciência e a Tecnologia (FCT) in the form of an exploratory project of reference IF/01037/2013CP1191/CT0001. RDH was supported by STFC studentship grant ST/J500744/1 during the course of this work. CAW acknowledges support from STFC grant ST/L000709/1. This publication was made possible by a grant from the John Templeton Foundation. The opinions expressed in this publication are those of the authors and do not necessarily reflect the views of the John Templeton Foundation. This material is based upon work supported by the National Aeronautics and Space Administration under Grant No. NNX15AC90G issued through the Exoplanets Research Program.


  • Andrieu & Thoms (2008) Andrieu, C. & Thoms, J. 2008, Statistics and Computing, 18, 343
  • Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, Annual Review of Astronomy and Astrophysics, 47, 481
  • Ballard et al. (2014) Ballard, S., Chaplin, W. J., Charbonneau, D., et al. 2014, ApJ, 790, 12
  • Baranne et al. (1996) Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373
  • Benz et al. (2014) Benz, W., Ida, S., Alibert, Y., Lin, D., & Mordasini, C. 2014, Protostars and Planets VI, 691
  • Bonfils et al. (2012) Bonfils, X., Gillon, M., Udry, S., et al. 2012, A&A, 546, A27
  • Bonomo et al. (2014) Bonomo, A. S., Sozzetti, A., Lovis, C., et al. 2014, A&A 1409,4592
  • Borucki et al. (2011) Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 736, 19
  • Boyajian et al. (2012) Boyajian, T. S., von Braun, K., van Belle, G., et al. 2012, ApJ, 757, 112
  • Buchhave et al. (2014) Buchhave, L. A., Bizzarro, M., Latham, D. W., et al. 2014, Nature, 509, 593
  • Buchhave et al. (2012) Buchhave, L. A., Latham, D. W., Johansen, A., et al. 2012, Nature, 486, 375
  • Butler et al. (2004) Butler, R. P., Vogt, S. S., Marcy, G. W., et al. 2004, ApJ, 617, 580
  • Castelli & Kurucz (2004) Castelli, F. & Kurucz, R. L. 2004, IAU Symp.210
  • Charbonneau et al. (2009) Charbonneau, D., Berta, Z. K., Irwin, J., et al. 2009, Nature, 462, 891
  • Claret (2000) Claret, A. 2000, A&A, 363, 1081
  • Claret & Bloemen (2011) Claret, A. & Bloemen, S. 2011, A&A, 529, A75
  • Collier Cameron et al. (2007) Collier Cameron, A., Wilson, D. M., West, R. G., et al. 2007, MNRAS, 380, 1230
  • Cosentino et al. (2014) Cosentino, R., Lovis, C., Pepe, F., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9147, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 8
  • Cosentino et al. (2012) Cosentino, R., Lovis, C., Pepe, F., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8446, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 1
  • Dawson & Fabrycky (2010) Dawson, R. I. & Fabrycky, D. C. 2010, ApJ, 722, 937
  • de Mooij et al. (2014) de Mooij, E. J. W., López-Morales, M., Karjalainen, R., Hrudkova, M., & Jayawardhana, R. 2014, ApJ, 797, L21
  • Demory et al. (2015a) Demory, B.-O., Ehrenreich, D., Queloz, D., et al. 2015a, MNRAS, 450, 2043
  • Demory et al. (2011) Demory, B.-O., Gillon, M., Deming, D., et al. 2011, A&A, 533, A114
  • Demory et al. (2015b) Demory, B.-O., Gillon, M., Madhusudhan, N., & Queloz, D. 2015b, ArXiv e-prints
  • Demory et al. (2012) Demory, B.-O., Gillon, M., Seager, S., et al. 2012, ApJ, 751, L28
  • Demory et al. (2013) Demory, B.-O., Torres, G., Neves, V., et al. 2013, ApJ, 768, 154
  • Dittmann et al. (2009) Dittmann, J. A., Close, L. M., Green, E. M., Scuderi, L. J., & Males, J. R. 2009, ApJ, 699, L48
  • Dragomir et al. (2012) Dragomir, D., Matthews, J. M., Howard, A. W., et al. 2012, ApJ, 759, L41
  • Dressing et al. (2015) Dressing, C. D., Charbonneau, D., Dumusque, X., et al. 2015, ApJ, 800, 135
  • Dumusque et al. (2014) Dumusque, X., Bonomo, A. S., Haywood, R. D., et al. 2014, ApJ, 789, 154
  • Dumusque et al. (2011a) Dumusque, X., Lovis, C., Ségransan, D., et al. 2011a, A&A, 535, A55
  • Dumusque et al. (2012) Dumusque, X., Pepe, F., Lovis, C., et al. 2012, Nature, 491, 207
  • Dumusque et al. (2011b) Dumusque, X., Udry, S., Lovis, C., Santos, N. C., & Monteiro, M. J. P. F. G. 2011b, A&A, 525, A140
  • Fabrycky et al. (2014) Fabrycky, D. C., Lissauer, J. J., Ragozzine, D., et al. 2014, ApJ, 790, 146
  • Fazio et al. (2004) Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10
  • Figueira et al. (2012) Figueira, P., Marmier, M., Boué, G., et al. 2012, A&A, 541, A139
  • Fortier et al. (2014) Fortier, A., Beck, T., Benz, W., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9143, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 2
  • Fressin et al. (2013) Fressin, F., Torres, G., Charbonneau, D., et al. 2013, ApJ, 766, 81
  • Gelman & Rubin (1992) Gelman, A. & Rubin, D. B. 1992, Statist. Sci., 7, 457
  • Gillon et al. (2010) Gillon, M., Deming, D., Demory, B.-O., et al. 2010, A&A, 518, A25
  • Gillon et al. (2007) Gillon, M., Demory, B.-O., Barman, T., et al. 2007, A&A, 471, L51
  • Gillon et al. (2012) Gillon, M., Demory, B.-O., Benneke, B., et al. 2012, A&A, 539, A28
  • Gillon et al. (2014) Gillon, M., Demory, B.-O., Madhusudhan, N., et al. 2014, A&A, 563, A21
  • Gregory (2005a) Gregory, P. C. 2005a, ApJ, 631, 1198
  • Gregory (2005b) Gregory, P. C. 2005b, Bayesian Logical Data Analysis for the Physical Sciences: A Comparative Approach with ‘Mathematica’ Support (Cambridge University Press)
  • Grillmair et al. (2012) Grillmair, C. J., Carey, S. J., Stauffer, J. R., et al. 2012, Pointing effects and their consequences for Spitzer IRAC exoplanet observations
  • Grimm & Stadel (2014) Grimm, S. L. & Stadel, J. G. 2014, ApJ, 796, 23
  • Howard et al. (2011) Howard, A. W., Johnson, J. A., Marcy, G. W., et al. 2011, ApJ, 730, 10
  • Howard et al. (2012) Howard, A. W., Marcy, G. W., Bryson, S. T., et al. 2012, ApJS, 201, 15
  • Howard et al. (2010) Howard, A. W., Marcy, G. W., Johnson, J. A., et al. 2010, Science, 330, 653
  • Ingalls et al. (2014) Ingalls, J. G., Carey, S. J., Lowrance, P. J., Grillmair, C. J., & Stauffer, J. R. 2014, Using drift scans to improve astrometry with Spitzer
  • Kass & Raftery (1995) Kass, E. R. & Raftery, E. A. 1995, Journal of the American Statistical Association, 90, 773
  • Knutson et al. (2008) Knutson, H. A., Charbonneau, D., Allen, L. E., Burrows, A., & Megeath, S. T. 2008, ApJ, 673, 526
  • Kovtyukh et al. (2003) Kovtyukh, V. V., Soubiran, C., Belik, S. I., & Gorlova, N. I. 2003, A&A, 411, 559
  • Kurucz (1992) Kurucz, R. L. 1992, in IAU Symposium, Vol. 149, The Stellar Populations of Galaxies, ed. B. Barbuy & A. Renzini, 225
  • Lanotte et al. (2014) Lanotte, A. A., Gillon, M., Demory, B.-O., et al. 2014, A&A, 572, A73
  • Lovis et al. (2011a) Lovis, C., Dumusque, X., Santos, N. C., et al. 2011a, ArXiv e-prints
  • Lovis et al. (2011b) Lovis, C., Ségransan, D., Mayor, M., et al. 2011b, A&A, 528, A112
  • Makino (1991) Makino, J. 1991, ApJ, 369, 200
  • Mandel & Agol (2002) Mandel, K. & Agol, E. 2002, ApJ, 580, L171
  • Marcus et al. (2010) Marcus, R. A., Sasselov, D., Stewart, S. T., & Hernquist, L. 2010, ApJ, 719, L45
  • Marcy et al. (2014) Marcy, G. W., Weiss, L. M., Petigura, E. A., et al. 2014, Proceedings of the National Academy of Science, 111, 12655
  • Mayor et al. (2009) Mayor, M., Bonfils, X., Forveille, T., et al. 2009, A&A, 507, 487
  • Mayor et al. (2011) Mayor, M., Marmier, M., Lovis, C., et al. 2011, A&A Submitted
  • Mayor et al. (2003) Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20
  • Meunier & Lagrange (2013) Meunier, N. & Lagrange, A.-M. 2013, A&A, 551, A101
  • Mishenina et al. (2008) Mishenina, T. V., Soubiran, C., Bienaymé, O., et al. 2008, A&A, 489, 923
  • Mishenina et al. (2012) Mishenina, T. V., Soubiran, C., Kovtyukh, V. V., Katsova, M. M., & Livshits, M. A. 2012, A&A, 547, A106
  • Mordasini et al. (2012) Mordasini, C., Alibert, Y., Georgy, C., et al. 2012, A&A, 547, A112
  • Mortier et al. (2014) Mortier, A., Sousa, S. G., Adibekyan, V. Z., Brandão, I. M., & Santos, N. C. 2014, A&A, 572, A95
  • Ogihara et al. (2015) Ogihara, M., Morbidelli, A., & Guillot, T. 2015, ArXiv e-prints
  • Pepe et al. (2013) Pepe, F., Cameron, A. C., Latham, D. W., et al. 2013, Nature, 503, 377
  • Pepe et al. (2011) Pepe, F., Lovis, C., Ségransan, D., et al. 2011, A&A, 534, A58
  • Pepe et al. (2002) Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632
  • Pollacco et al. (2008) Pollacco, D., Skillen, I., Collier Cameron, A., et al. 2008, MNRAS, 385, 1576
  • Prugniel et al. (2011) Prugniel, P., Vauglin, I., & Koleva, M. 2011, A&A, 531, A165
  • Ramírez et al. (2013) Ramírez, I., Allende Prieto, C., & Lambert, D. L. 2013, ApJ, 764, 78
  • Rauer et al. (2014) Rauer, H., Catala, C., Aerts, C., et al. 2014, Experimental Astronomy, 38, 249
  • Ricker et al. (2014) Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9143, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 20
  • Schwarz (1978) Schwarz, G. 1978, Ann. Statist., 6, 461
  • Sneden (1973) Sneden, C. 1973, ApJ, 184, 839
  • Sousa et al. (2011) Sousa, S. G., Santos, N. C., Israelian, G., Mayor, M., & Udry, S. 2011, A&A, 533, A141
  • Stetson (1987) Stetson, P. B. 1987, PASP, 99, 191
  • Stevenson et al. (2012) Stevenson, K. B., Harrington, J., Fortney, J. J., et al. 2012, ApJ, 754, 136
  • Tsantaki et al. (2013) Tsantaki, M., Sousa, S. G., Adibekyan, V. Z., et al. 2013, A&A, 555, A150
  • Valenti & Fischer (2005) Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141
  • Van Grootel et al. (2014) Van Grootel, V., Gillon, M., Valencia, D., et al. 2014, ApJ, 786, 2
  • van Leeuwen (2007) van Leeuwen, F. 2007, A&A, 474, 653
  • Zechmeister & Kürster (2009) Zechmeister, M. & Kürster, M. 2009, A&A, 496, 577
  • Zeng et al. (2015) Zeng, L., Jacobsen, S., & Sasselov, D. D. 2015, in American Astronomical Society Meeting Abstracts, Vol. 225, American Astronomical Society Meeting Abstracts, 406.02
  • Zeng & Sasselov (2013) Zeng, L. & Sasselov, D. 2013, PASP, 125, 227

Note added in proof During the refereeing process, we learned about an independent detection by Vogt et al. (Laughlin, private communication) reporting additional planets in the system, based on long-term radial velocities obtained with the Keck and APF telescopes.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description