# The HARPS search for southern extrasolar planets.
## XXIII. 8 planetary companions to low-activity solar-type stars

Dominique Naef, Michel Mayor, Gaspare Lo Curto, Francois Bouchy, Christophe Lovis, Claire Moutou, Willy Benz, Francesco Pepe, Didier Queloz, Nuno C. Santos, et al.

► To cite this version:

Dominique Naef, Michel Mayor, Gaspare Lo Curto, Francois Bouchy, Christophe Lovis, et al. The HARPS search for southern extrasolar planets. XXIII. 8 planetary companions to low-activity solar-type stars. Astronomy & Astrophysics - A&A, 2010, 523, pp.A15. 10.1051/0004-6361/200913616. hal-00662790

HAL Id: hal-00662790

https://hal.science/hal-00662790v1

Submitted on 30 Nov 2021

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# The HARPS search for southern extrasolar planets*,**

## XXIII. 8 planetary companions to low-activity solar-type stars

D. Naef<fnref n="1" />, M. Mayor<fnref n="1" />, G. Lo Curto<fnref n="2" />, F. Bouchy<fnref n="3" />,<fnref n="4" />, C. Lovis<fnref n="1" />, C. Moutou<fnref n="5" />, W. Benz<fnref n="6" />, F. Pepe<fnref n="1" />, D. Queloz<fnref n="1" />, N.C. Santos<fnref n="7" />, D. Ségransan<fnref n="1" />, S. Udry<fnref n="1" />, X. Bonfils<fnref n="1" />,<fnref n="8" />, X. Delfosse<fnref n="8" />, T. Forveille<fnref n="8" />, G. Hébrard<fnref n="3" />, C. Mordasini<fnref n="6" />,<fnref n="9" />, C. Perrier<fnref n="8" />, I. Boisse<fnref n="3" />, and D. Sosnowska<fnref n="1" />

Received 6 November 2009 / Accepted 10 August 2010

## ABSTRACT

In this paper, we present our HARPS radial-velocity data for eight low-activity solar-type stars belonging to the HARPS volume-limited sample: HD 6718, HD 8535, HD 28254, HD 290327, HD 43197, HD 44219, HD 148156, and HD 156411. Keplerian fits to these data reveal the presence of low-mass companions around these targets. With minimum masses ranging from 0.58 to 2.54 $M_{Jup}$, these companions are in the planetary mass domain. The orbital periods of these planets range from slightly less than one to almost seven years. The eight orbits presented in this paper exhibit a wide variety of eccentricities: from 0.08 to above 0.8.

**Key words.** planetary systems – techniques: radial velocities – techniques: spectroscopic

## 1. Introduction

In this paper, we present new detections from the HARPS planet search programme. It has been ongoing since 2003. Its detections are based on radial velocities obtained with the *High Accuracy Radial velocity Planet Searcher* (HARPS, Mayor et al. 2003). HARPS is a high-resolution echelle spectrograph mounted on the ESO-3.6-m telescope in La Silla (Chile). The instrument is installed in a vacuum vessel, and it is thermally controlled. These characteristics result in a very high stability of the intrumental radial-velocity zero point (measured drifts smaller than 1 m s⁻<fnref n="1" /> over one night). HARPS possesses two optical fibres. The first one collects the target light, whereas the second one can be used for recording a reference Thorium-Argon spectrum. The second spectrum has been used for measuring the residual drift of the instrument since the time of the latest wavelength calibration, allowing us to apply a correction to the velocities (Baranne et al. 1996). The second fibre can also be used for measuring the spectrum of the sky background. With more than 70 detections

obtained so far, including the lightest planet (Gl 581 e, Mayor et al. 2009), the HARPS contribution to the exoplanet research field is a major one.

Since its beginning, the HARPS planet search programme has been divided into several sub-programmes related to specific scientific questions. One of these sub-programmes is the HARPS volume-limited programme to which about 20% of the HARPS Guaranteed Time Observation (GTO) programme allocated time was dedicated. The sample of the HARPS volume-limited programme is an extension up to 57.5 pc of the CORALIE planet search sample (volume-limited sample up to 50 pc, see Udry et al. 2000). It contains about 850 F8-M0 dwarfs. As in the case of the CORALIE sample, the maximum distance decreases with the *B* − *V* colour index for targets with late spectral types (later than K0). Our target list will be published in another paper (Lo Curto et al., in prep.).

Targets in this sample are observed at lower signal-to-noise ratios (typically 40–50) than the ones in the high-precision programme and without using the simultaneous calibration. The typical photon-noise error obtained in this case is about 1–2 m s⁻<fnref n="1" /> depending on the target spectral type and projected rotational velocity. This strategy, optimized for the follow-up of a large stellar sample, is a good compromise between precision and observing efficiency.

Rather short exposure times (less than 5 min) are usually enough for reaching the targetted signal-to-noise ratio. An

<footnote n="1">¹ Observatoire Astronomique de l'Université de Genève, 51 Ch. des Maillettes, 1290 Sauverny, Switzerland
e-mail: dominique.naef@unige.ch</footnote>

<footnote n="2">² European Southern Observatory, Karl Schwarzschild Strasse 2, 85748 Garching bei München, Germany</footnote>

<footnote n="3">³ Institut d'Astrophysique de Paris, UMR7095 CNRS, Université Pierre & Marie Curie, 98bis boulevard Arago, 75014 Paris, France</footnote>

<footnote n="4">⁴ Observatoire de Haute-Provence, 04870 St-Michel L'Observatoire, France</footnote>

<footnote n="5">⁵ Laboratoire d'Astrophysique de Marseille, OAMP, Université Aix-Marseille and CNRS, 38 rue Frédéric Jolliot-Curie, 13388 Marseille Cedex 13, France</footnote>

<footnote n="6">⁶ Physikalisches Institut Universität Bern, Sidlerstrasse 5, 3012 Bern, Switzerland</footnote>

<footnote n="7">⁷ Centro de Astrofísica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal</footnote>

<footnote n="8">⁸ Laboratoire d'Astrophysique de Grenoble, UMR5571 CNRS, Université Joseph Fourier, BP 53, 38041 Grenoble Cedex 9, France</footnote>

<footnote n="9">⁹ Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany</footnote>

<footnote>* Based on observations made with the HARPS instrument on the ESO 3.6-m telescope at the La Silla Observatory (Chile) under the GTO programme ID 072.C-0488.</footnote>

<footnote>** Radial-velocity tables are only available in electronic form 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/523/A15</footnote>
exposure time of 5 min is too short for correctly averaging out
intrinsic stellar signals due for example to pulsation or granula-
tion. Our velocities can thus be affected by those effects that in
some cases can reach a peak-to-peak amplitude of several m s⁻¹
(see e.g. the case of µ Arae, Santos et al. 2004a; Bouchy et al.
2005). The amplitudes of these signals depend on the target spec-
tral types and on the averaging strategy. Expected radial-velocity
errors due to intrinsic signals are greater for earlier spectral types
and for evolved targets.

With such an observing strategy, our programme is mostly
sensitive to Saturn or Jupiter-like planets but is also sensitive
to lower-mass planets around stars with a low intrinsic radial-
velocity variability. Also, we generally avoid taking exposures
shorter than 1 min in order to reduce the impact of observing
overheads. As a consequence, the brightest targets in our sam-
ple frequently have photon-noise errors significantly lower than
the targetted one. For illustrating these facts, we mention the re-
cent detections of a 5.5 M⊕ planet around HD 215497 (Lo Curto
et al. 2010), of the 12.6 M⊕ planet HD 125595 b (Ségransan et al.
2010) and of the 14.4 M⊕ planet around BD −08 2823 b (Hébrard
et al. 2010).

The main aim of this observational effort is a better char-
acterization of orbital elements distributions of gaseous planets.
With the increasing size of the sample of known planets, a quan-
titative comparison between the observed sample and the predic-
tions of planet formation models becomes possible and reliable
(see for example Mordasini et al. 2009). With the eight planetary
companions presented in this paper, the total number of planets
detected in the volume-limited sample now reaches 39 (see Pepe
et al. 2004a; Moutou et al. 2005, 2009b; Lo Curto et al. 2006,
2010; Naef et al. 2007; Mordasini et al. 2010; Ségransan et al.
2010).

In this paper, we present our HARPS radial-velocity
data for 8 solar-type stars hosting low-mass compan-
ions: HD 6718, HD 8535, HD 28254, HD 290327, HD 43197,
HD 44219, HD 148156, and HD 156411. This paper is organized
as follows. We present the main characteristics of the 8 host
stars in Sect. 2. Section 3 contains the description of the radial-
velocity sets and orbital solutions obtained for these targets.
Finally, we provide concluding remarks in Sect. 4.

In this paper, we present our HARPS radial-velocity
data for 8 solar-type stars hosting low-mass compan-
ions: HD 6718, HD 8535, HD 28254, HD 290327, HD 43197,
HD 44219, HD 148156, and HD 156411. This paper is organized
as follows. We present the main characteristics of the 8 host
stars in Sect. 2. Section 3 contains the description of the radial-
velocity sets and orbital solutions obtained for these targets.
Finally, we provide concluding remarks in Sect. 4.

## 2. Characteristics of the host stars

We list the main characteristics of the nine host stars in Table 1.
Apparent V-band magnitudes *m*<sub>V</sub> and colour indexes *B* − *V* are
from the HIPPARCOS Catalogue (ESA 1997). Distances *d* are
based on the improved HIPPARCOS astrometric parallaxes *π* de-
rived by van Leeuwen (2007). Absolute magnitudes *M*<sub>V</sub> were
computed using *m*<sub>V</sub> and *d*.

For each star, we derived the atmospheric parameters T<sub>eff</sub>,
log g, and [Fe/H] by performing a detailed local thermodynam-
ical equilibrium spectroscopic analysis. A high signal-to-noise
co-added HARPS spectrum was built for each target. The set of
Fe I/Fe II equivalent widths used for the analysis were measured
on the co-added spectra with the ARES code (Sousa et al. 2007).
The list of spectral lines for the analysis, as well as more details
on the method can be found in Santos et al. (2004b).

Stellar luminosities *L* were computed using HIPPARCOS
apparent magnitudes and parallaxes and the bolometric correc-
tion obtained via the *T*<sub>eff</sub>-based calibration in Flower (1996).
Stellar masses *M*, and radii *R*<sub>*</sub> were derived from *T*<sub>eff</sub>, [Fe/H],
and *M*<sub>V</sub> by interpolating the Girardi et al. (2000) theoretical
isochrones. Following Fernandes & Santos (2004), a relative

uncertainty of ≈10% can be assumed for *M*<sub>*</sub>. We also com-
puted surface gravities, log *g*, using *M*<sub>*</sub> and *R*<sub>*</sub>. These values
are in fair agreement with the spectroscopic ones: a mean dif-
ference (spectroscopy–isochrones) of 0.06 dex and an rms of
the differences equal to 0.11 dex. Spectral types listed in Table 1
are based on the stellar parameters we determined in this pa-
per. Jenkins et al. (2008) have published iron abundances for
HD 43197 and HD 148156. Their values, [Fe/H] = 0.37 ± 0.09
and 0.22 ± 0.15, respectively, are in good agreement with ours.
The atmospheric parameters for HD 44219 have also been deter-
mined by Robinson et al. (2007). They found T<sub>eff</sub> = 5723 ± 82,
log *g* = 4.37 ± 0.13 and [Fe/H] = −0.05 ± 0.07. These values
agree with ours within uncertainties.

Activity indexes log R'_{HK} were extracted from individual HARPS spectra using a method similar to the one described in Santos et al. (2000). Values indicated in Table 1 are average activity indexes obtained for each target. All the targets exhibit a very low level of activity. Our log R'_{HK} values for HD 43197 and HD 148156 are in good agreement with those published in Jenkins et al. (2008): -5.12 and -5.08, respectively. We did not try to estimate chromospheric ages for our targets because these two quantities are no longer correlated for weakly active stars (Pace & Pasquini 2004). However, we can infer for most of these stars (if not all) a safe age constraint based on their low activity indexes and rather high luminosities: age > 3 Gyr. The projected rotational velocities v sin i were derived from the widths of the HARPS cross-correlation functions (CCF) following the method presented in Santos et al. (2002).

Finally, HD 28254 is the bright component of a visual bi-
nary. According to the Catalogue of Components of Double and
Multiple stars (Dommanget & Nys 2002), the visual companion,
CCDMJ04248-5037B, is located at a separation of 4.3" and at
a position angle of 254 degrees. With an apparent magnitude of
*m*<sub>V</sub> = 13.8, this object is 6.11 mag fainter (i.e. a flux ratio less
than 0.4%). It was thus clearly outside the detection capabilities
of HIPPARCOS. Assuming the two stars form a bound system
and using the distance of HD 28254 (54.7 pc), we computed the
projected separation between them: 235 AU. We also computed
the absolute magnitude of the B component: *M*<sub>V</sub> = 10.11 corre-
sponding to an M0V-M2V star. With this absolute magniitude
and the mass-luminosity relations in Delfosse et al. (2000), we
computed a mass of 0.48 *M*<sub>☉</sub> for this object. Its minimum period
is about 1000 yr. Assuming that this companion is on a circular
edge-on orbit, the maximal radial-velocity amplitude it would
induce on the primary star is of the order of 1.1 km s<sup>-1</sup>.

3. HARPS radial-velocity data and orbital solutions

As explained before, stars belonging to the HARPS volume-
limited sample are generally observed without simultaneous
thorium-argon reference. Possible instrumental drifts therefore
remain uncorrected. This has only a limited impact on our re-
sults because: 1) the instrument is very stable: on average, the
radial-velocity drifts of HARPS remain below 0.5 m s⁻¹ over one
night and 2) the targetted photon-noise error for this programme
is 1–2 m s⁻¹. We stress again the fact that our observing strat-
egy is not optimally designed for averaging out stellar intrinsic
signals. We list in Table 2 the median exposure times for all the
targets. This median is below 3 min in all cases. We thus expect
our velocities to be affected at some level by stellar intrinsic jit-
ter in particular the velocities of our early type and/or slightly
evolved targets, namely HD 8535 (G0V), HD 28254 (G1IV/V),
HD 148156 (F8V), and HD 156411 (F8IV/V).
<table><caption>**Table 1.** Observed and inferred stellar characteristics of the host stars.</caption><thead><tr><th>Parameter</th><th>Unit</th><th>HD 6718</th><th>HD 8535</th><th>HD 28254</th><th>HD 290327</th></tr></thead><tbody><tr><td>Sp. Type</td><td></td><td>G5V</td><td>G0V</td><td>G1IV/V</td><td>G8V</td></tr><tr><td>m<sub>V</sub></td><td></td><td>8.45</td><td>7.70</td><td>7.69</td><td>8.99</td></tr><tr><td>B-V</td><td></td><td>0.662</td><td>0.553</td><td>0.722</td><td>0.761</td></tr><tr><td>π</td><td>(mas)</td><td>18.23 ± 0.76</td><td>19.03 ± 0.60</td><td>18.29 ± 0.53</td><td>17.65 ± 1.57</td></tr><tr><td>d</td><td>(pc)</td><td>54.9<sup>+2.4</sup><sub>-2.2</sub></td><td>52.5<sup>+1.7</sup><sub>-1.6</sub></td><td>54.7<sup>+1.6</sup><sub>-1.5</sub></td><td>56.7<sup>+5.5</sup><sub>-4.6</sub></td></tr><tr><td>M<sub>V</sub></td><td></td><td>4.754</td><td>4.097</td><td>4.001</td><td>5.224</td></tr><tr><td>B.C.</td><td></td><td>-0.145</td><td>-0.028</td><td>-0.102</td><td>-0.126</td></tr><tr><td>L</td><td>(L<sub>☉</sub>)</td><td>1.13</td><td>1.86</td><td>2.18</td><td>0.72</td></tr><tr><td>T<sub>eff</sub></td><td>(K)</td><td>5746 ± 19</td><td>6136 ± 18</td><td>5664 ± 35</td><td>5552 ± 21</td></tr><tr><td>log g<sup>a</sup></td><td>(cgs)</td><td>4.48 ± 0.03</td><td>4.46 ± 0.06</td><td>4.12 ± 0.05</td><td>4.42 ± 0.04</td></tr><tr><td>log g<sup>b</sup></td><td>(cgs)</td><td>4.40</td><td>4.35</td><td>4.13</td><td>4.39</td></tr><tr><td>[Fe/H]</td><td></td><td>-0.06 ± 0.02</td><td>0.06 ± 0.02</td><td>0.36 ± 0.03</td><td>-0.11 ± 0.02</td></tr><tr><td>M<sub>*</sub></td><td>(M<sub>☉</sub>)</td><td>0.96</td><td>1.13</td><td>1.06</td><td>0.90</td></tr><tr><td>R<sub>*</sub></td><td>(R<sub>☉</sub>)</td><td>1.02 ± 0.03</td><td>1.19 ± 0.04</td><td>1.48 ± 0.06</td><td>1.00 ± 0.01</td></tr><tr><td>log R'<sub>HK</sub></td><td></td><td>-4.97</td><td>-4.95</td><td>-5.10</td><td>-4.96</td></tr><tr><td>v sin i</td><td>(km s<sup>-1</sup>)</td><td>1.76 ± 1.0</td><td>1.41 ± 1.0</td><td>2.50 ± 1.0</td><td>1.44 ± 1.0</td></tr><tr><td>Parameter</td><td>Unit</td><td>HD 43197</td><td>HD 44219</td><td>HD 148156</td><td>HD 156411</td></tr><tr><td>Sp. Type</td><td></td><td>G8V</td><td>G2V</td><td>F8V</td><td>F8IV/V</td></tr><tr><td>m<sub>V</sub></td><td></td><td>8.98</td><td>7.69</td><td>7.69</td><td>6.67</td></tr><tr><td>B-V</td><td></td><td>0.817</td><td>0.687</td><td>0.560</td><td>0.614</td></tr><tr><td>π</td><td>(mas)</td><td>17.76 ± 1.22</td><td>19.83 ± 0.78</td><td>19.38 ± 0.75</td><td>18.25 ± 0.49</td></tr><tr><td>d</td><td>(pc)</td><td>56.3<sup>+4.2</sup><sub>-3.6</sub></td><td>50.4<sup>+2.1</sup><sub>-1.9</sub></td><td>51.6<sup>+2.1</sup><sub>-1.9</sub></td><td>54.8<sup>+1.5</sup><sub>-1.4</sub></td></tr><tr><td>M<sub>V</sub></td><td></td><td>5.227</td><td>4.177</td><td>4.127</td><td>2.976</td></tr><tr><td>B.C.</td><td></td><td>-0.137</td><td>-0.085</td><td>-0.009</td><td>-0.060</td></tr><tr><td>L</td><td>(L<sub>☉</sub>)</td><td>0.73</td><td>1.82</td><td>1.78</td><td>5.38</td></tr><tr><td>T<sub>eff</sub></td><td>(K)</td><td>5508 ± 46</td><td>5752 ± 16</td><td>6308 ± 28</td><td>5900 ± 15</td></tr><tr><td>log g<sup>a</sup></td><td>(cgs)</td><td>4.31 ± 0.08</td><td>4.21 ± 0.03</td><td>4.56 ± 0.07</td><td>4.07 ± 0.04</td></tr><tr><td>log g<sup>b</sup></td><td>(cgs)</td><td>4.42</td><td>4.20</td><td>4.36</td><td>3.87</td></tr><tr><td>[Fe/H]</td><td></td><td>0.40 ± 0.04</td><td>0.03 ± 0.01</td><td>0.29 ± 0.02</td><td>-0.12 ± 0.01</td></tr><tr><td>M<sub>*</sub></td><td>(M<sub>☉</sub>)</td><td>0.96</td><td>1.00</td><td>1.22</td><td>1.25</td></tr><tr><td>R<sub>*</sub></td><td>(R<sub>☉</sub>)</td><td>1.00 ± 0.01</td><td>1.32 ± 0.07</td><td>1.21 ± 0.03</td><td>2.16 ± 0.09</td></tr><tr><td>log R'<sub>HK</sub></td><td></td><td>-5.06</td><td>-5.03</td><td>-4.94</td><td>-5.05</td></tr><tr><td>v sin i</td><td>(km s<sup>-1</sup>)</td><td>2.18 ± 1.0</td><td>2.22 ± 1.0</td><td>5.70 ± 1.0</td><td>3.30 ± 1.0</td></tr></tbody></table>

**Notes.** (a) From our spectroscopic LTE analysis; (b) using $M_*$ and $R_*$.

Our radial velocities can be obtained from the extracted spectra by cross-correlating them with a numerical template. Numerical templates are chosen to minimize the mismatch with the target spectral type. Our radial-velocity uncertainties include two well-quantified contributions: photon noise and calibration error. An additional error of 0.5 m s⁻¹ is added quadratically to account for the absence of simultaneous calibration and thus the possibility of a small uncorrected radial-velocity drift of the HARPS zero point.

All the orbital solutions presented in the following sections were obtained with software that includes a genetic algorithm. In general, the only a priori information used by this algorithm is the type of model we want to use (i.e., 1-Keplerian, 2-Keplerian, 1-Keplerian plus a linear/quadratic drift, etc.). Once the global parameter space minimum is identified, we use it as the starting point for a standard Levenberg-Marquardt minimization out of which the final solution is obtained. In all the cases, we have verified the absence of correlation between the post-fit residuals and the measured cross-correlation function bisectors.

Finally, uncertainties on the fitted parameters are determined using Monte-Carlo simulations.

### 3.1. HD 6718

The 22 HARPS radial-velocity measurements were gathered for HD 6718 (HIP 5301) between BJD = 2452991 (December 18, 2003) and BJD = 2455019 (July 6, 2009). These velocities have a mean uncertainty of 1.59 m s⁻¹. We have fitted a long-period Keplerian orbit to these data. The orbital parameters resulting from this fit are listed in Table 2. With these parameters and assuming a mass of 0.96 $M_☉$ for the host star, we computed a minimum mass of 1.56 $M_{Jup}$ for the companion and a semi-major axis of 3.56 AU. Residuals to the fitted orbit have a very low dispersion ($\sigma_{O-C} = 1.79$ m s⁻¹), and a 1-planet model is thus sufficient for explaining the detected signal. We display our data and the fitted orbit in panel a of Fig. 1.
<table><caption>**Table 2.** Orbital solutions for the eight host stars.</caption><thead><tr><th>Parameter</th><th>Unit</th><th>HD 6718</th><th>HD 8535</th><th>HD 28254</th><th>HD 290327</th></tr></thead><tbody><tr><td>P</td><td>(d)</td><td>2496 ± 176</td><td>1313 ± 28</td><td>1116 ± 26</td><td>2443<sup>+205</sup><sub>-117</sub></td></tr><tr><td>T</td><td>(JD<sup>a</sup>)</td><td>357 ± 251</td><td>537 ± 92</td><td>49<sup>+43</sup><sub>-55</sub></td><td>1326<sup>+230</sup><sub>-412</sub></td></tr><tr><td>e</td><td>(km s<sup>-1</sup>)</td><td>0.10<sup>+0.11</sup><sub>-0.04</sub></td><td>0.15<sup>+0.09</sup><sub>-0.05</sub></td><td>0.81<sup>+0.05</sup><sub>-0.02</sub></td><td>0.08<sup>+0.08</sup><sub>-0.03</sub></td></tr><tr><td>γ</td><td>(km s<sup>-1</sup>)</td><td>34.7509 ± 0.0013</td><td>2.4588 ± 0.0008</td><td>-9.315<sup>b</sup> ± 0.002</td><td>29.559 ± 0.003</td></tr><tr><td>ω</td><td>(deg)</td><td>286<sup>+64</sup><sub>-35</sub></td><td>61<sup>+18</sup><sub>-43</sub></td><td>301<sup>+3</sup><sub>-4</sub></td><td>268<sup>+94</sup><sub>-19</sub></td></tr><tr><td>K<sub>1</sub></td><td>(m s<sup>-1</sup>)</td><td>24.1 ± 1.5</td><td>11.8 ± 0.8</td><td>37.3<sup>+5.1</sup><sub>-0.9</sub></td><td>41.3 ± 2.9</td></tr><tr><td>Slope</td><td>(m s<sup>-1</sup> yr<sup>-1</sup>)</td><td>...</td><td>...</td><td>2.629 ± 0.015</td><td>...</td></tr><tr><td>Curvature</td><td>(m s<sup>-1</sup> yr<sup>-2</sup>)</td><td>...</td><td>...</td><td>-0.9290 ± 0.00012</td><td>...</td></tr><tr><td>f(m)</td><td>(10<sup>-9</sup> M<sub>☉</sub>)</td><td>3.56<sup>+0.71</sup><sub>-0.52</sub></td><td>0.22<sup>+0.06</sup><sub>-0.04</sub></td><td>1.22<sup>+0.30</sup><sub>-0.17</sub></td><td>17.6<sup>+3.2</sup><sub>-2.3</sub></td></tr><tr><td>a<sub>1</sub> sin i</td><td>(10<sup>-3</sup> AU)</td><td>5.50<sup>+0.45</sup><sub>-0.24</sub></td><td>1.41<sup>+0.11</sup><sub>-0.09</sub></td><td>2.25<sup>+0.19</sup><sub>-0.13</sub></td><td>9.24<sup>+0.61</sup><sub>-0.29</sub></td></tr><tr><td>m<sub>2</sub> sin i</td><td>(M<sub>Jup</sub>)</td><td>1.56<sup>+0.11</sup><sub>-0.10</sub></td><td>0.68<sup>+0.07</sup><sub>-0.04</sub></td><td>1.16<sup>+0.10</sup><sub>-0.06</sub></td><td>2.54<sup>+0.17</sup><sub>-0.14</sub></td></tr><tr><td>a</td><td>(AU)</td><td>3.56<sup>+0.24</sup><sub>-0.15</sub></td><td>2.45<sup>+0.04</sup><sub>-0.06</sub></td><td>2.15<sup>+0.04</sup><sub>-0.05</sub></td><td>3.43<sup>+0.20</sup><sub>-0.12</sub></td></tr><tr><td>N<sub>mes</sub></td><td></td><td>22</td><td>45</td><td>32</td><td>18</td></tr><tr><td>Med(t<sub>exp</sub>)</td><td>(s)</td><td>158</td><td>106</td><td>60</td><td>173</td></tr><tr><td>⟨S/N⟩</td><td></td><td>52</td><td>58</td><td>51</td><td>45</td></tr><tr><td>⟨ε<sub>RV</sub>⟩</td><td>(m s<sup>-1</sup>)</td><td>1.59</td><td>1.57</td><td>1.34</td><td>1.73</td></tr><tr><td>Span</td><td>(d)</td><td>2028</td><td>2220</td><td>1989</td><td>1986</td></tr><tr><td>σ<sub>O-C</sub></td><td>(m s<sup>-1</sup>)</td><td>1.79</td><td>2.49</td><td>2.19</td><td>1.60</td></tr><tr><td>χ<sup>2</sup><sub>red</sub></td><td></td><td>1.50</td><td>2.48</td><td>3.08</td><td>1.22</td></tr><tr><th>Parameter</th><th>Unit</th><th>HD 43197</th><th>HD 44219</th><th>HD 148156</th><th>HD 156411</th></tr><tr><td>P</td><td>(d)</td><td>327.8 ± 1.2</td><td>472.3<sup>+6.3</sup><sub>-5.0</sub></td><td>1027 ± 28</td><td>842.2 ± 14.5</td></tr><tr><td>T</td><td>(JD<sup>a</sup>)</td><td>713.8 ± 1.7</td><td>585.6<sup>+7.2</sup><sub>-3.1</sub></td><td>707 ± 20</td><td>356<sup>+46</sup><sub>-33</sub></td></tr><tr><td>e</td><td>(km s<sup>-1</sup>)</td><td>0.83<sup>+0.05</sup><sub>-0.01</sub></td><td>0.61<sup>+0.07</sup><sub>-0.09</sub></td><td>0.52<sup>+0.04</sup><sub>-0.09</sub></td><td>0.22 ± 0.08</td></tr><tr><td>γ</td><td>(km s<sup>-1</sup>)</td><td>72.512<sup>+4</sup><sub>-11</sub></td><td>-12.0732<sup>c</sup> ± 0.0026</td><td>-1.7480<sup>+0.0013</sup><sub>-0.0009</sub></td><td>-38.7512 ± 0.0011</td></tr><tr><td>ω</td><td>(deg)</td><td>251<sup>+7</sup><sub>-12</sub></td><td>147.4<sup>+6.2</sup><sub>-3.2</sub></td><td>35<sup>+17</sup><sub>-12</sub></td><td>45<sup>+22</sup><sub>-14</sub></td></tr><tr><td>K<sub>1</sub></td><td>(m s<sup>-1</sup>)</td><td>32.4<sup>+10.9</sup><sub>-3.3</sub></td><td>19.4 ± 3.0</td><td>17.5<sup>+1.5</sup><sub>-0.5</sub></td><td>14.0 ± 0.8</td></tr><tr><td>Slope</td><td>(m s<sup>-1</sup> yr<sup>-1</sup>)</td><td>...</td><td>-0.637 ± 0.016</td><td>...</td><td>...</td></tr><tr><td>Curvature</td><td>(m s<sup>-1</sup> yr<sup>-2</sup>)</td><td>...</td><td>-0.76390 ± 0.000</td><td>...</td><td>...</td></tr><tr><td>f(m)</td><td>(10<sup>-9</sup> M<sub>☉</sub>)</td><td>0.21<sup>+0.13</sup><sub>-0.04</sub></td><td>0.17<sup>+0.05</sup><sub>-0.03</sub></td><td>0.35<sup>+0.08</sup><sub>-0.04</sub></td><td>0.23<sup>+0.04</sup><sub>-0.03</sub></td></tr><tr><td>a<sub>1</sub> sin i</td><td>(10<sup>-3</sup> AU)</td><td>0.55<sup>+0.10</sup><sub>-0.04</sub></td><td>0.66<sup>+0.06</sup><sub>-0.04</sub></td><td>1.41<sup>+0.09</sup><sub>-0.06</sub></td><td>1.06<sup>+0.08</sup><sub>-0.05</sub></td></tr><tr><td>m<sub>2</sub> sin i</td><td>(M<sub>Jup</sub>)</td><td>0.60<sup>+0.12</sup><sub>-0.04</sub></td><td>0.58<sup>+0.06</sup><sub>-0.04</sub></td><td>0.85<sup>+0.06</sup><sub>-0.05</sub></td><td>0.74<sup>+0.05</sup><sub>-0.04</sub></td></tr><tr><td>a</td><td>(AU)</td><td>0.92<sup>+0.01</sup><sub>-0.02</sub></td><td>1.19 ± 0.02</td><td>2.45<sup>+0.04</sup><sub>-0.05</sub></td><td>1.88<sup>+0.03</sup><sub>-0.04</sub></td></tr><tr><td>N<sub>mes</sub></td><td></td><td>33</td><td>46</td><td>42</td><td>50</td></tr><tr><td>Med(t<sub>exp</sub>)</td><td>(s)</td><td>171</td><td>63</td><td>92</td><td>90</td></tr><tr><td>⟨S/N⟩</td><td></td><td>48</td><td>51</td><td>56</td><td>73</td></tr><tr><td>⟨ε<sub>RV</sub>⟩</td><td>(m s<sup>-1</sup>)</td><td>1.42</td><td>1.55</td><td>2.24</td><td>1.26</td></tr><tr><td>Span</td><td>(d)</td><td>1943</td><td>1988</td><td>2168</td><td>2231</td></tr><tr><td>σ<sub>O-C</sub></td><td>(m s<sup>-1</sup>)</td><td>1.44</td><td>2.39</td><td>3.69</td><td>2.94</td></tr><tr><td>χ<sup>2</sup><sub>red</sub></td><td></td><td>1.19</td><td>2.62</td><td>2.89</td><td>5.21</td></tr></tbody></table>

**Notes.** Errobars are Monte-Carlo based 1σ uncertainties. ⟨ε<sub>RV</sub>⟩ is the radial-velocity uncertainty weighted average. σ<sub>O-C</sub> is the weighted rms of the residuals to the fitted orbits. Med(t<sub>exp</sub>) is the median exposure time and ⟨S/N⟩ is the average signal-to-noise ratio at 550 nm.

## 3.2. HD 8535

A set of 45 HARPS radial velocities was obtained for HD 8535 (HIP 6511) between BJD = 2 452 850 (July 29, 2003) and BJD = 2 455 070 (August 26, 2009). The mean radial-velocity uncertainty of these measurements is 1.57 m s<sup>-1</sup>. The presence of a long-period signal is obvious (see panel b of Fig. 1). The results of the Keplerian fit to our data are presented in Table 2. The

dispersion of the residuals to the fitted orbit are somewhat larger than expected (χ<sup>2</sup><sub>red</sub> = 2.48). We have tried to fit orbits including additional signals (linear or quadratic trend; 2-Keplerian model) but we were unable to significantly improve the quality of the solution. We have also verified the absence of any correlation between the residuals and the CCF bisectors.

Assuming a primary mass of 1.13 M<sub>☉</sub>, we computed the minimum mass of the companion and the semi-major axis of its

<footnote>(a) JD = BJD − 2 454 000; (b) at BJD = 2 454 106.6; (c) at BJD = 2 454 512.6.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** HARPS radial-velocity data for the eight stars presented in this paper. Temporal velocities (filled dots) and fitted orbital solutions (solid lines) are displayed in the *upper panels*. Residuals to the fitted orbits are displayed in the *bottom panels*. For HD 28254 and HD 44219, the dashed lines in the *upper panels* show additional quadratic radial-velocity drifts that were included to the orbital fits. For all the other targets, the *upper panels* dashed lines represent the systemic velocity $v$.</figcaption></figure>

orbit: $m_2 \sin i = 0.68 M_{Jup}$ and $a = 2.45$ AU. The fitted orbit and the fit residuals are displayed in panel **b** of Fig. 1. The phase-folded orbit and velocities are displayed in panel **a** of Fig. 2.

### 3.3. HD 28254

We collected 32 HARPS radial velocities for HD 28254 (HIP 20606) between BJD = 2452941 (October 28, 2003) and BJD = 2454931 (April 9, 2009). The average velocity uncertainty of these measurements is $1.34 \text{ m s}^{-1}$. As indicated in Sect. 2, HD 28254 is the main component of a visual binary system. The presence of a nearby visual companion can induce radial-velocity errors via the seeing-dependent contamination of the spectrum of the main component. Pepe et al. (2004b) have

simulated this effect for HARPS and using the La Silla seeing statistics. We have performed the same simulations using the characteristics of the HD 28254 system ($\rho = 4.3''$, $\Delta m_v = 6.11$). The expected effect is negligible in this case ($\ll 0.1 \text{ m s}^{-1}$).

A long-period signal is visible on the data (see panel **c** of Fig. 1). From the visual inspection of these data, it is also obvious that the Keplerian signal is highly eccentric. We first attempted to fit a single orbit to our velocities. With $e > 0.99$, the obtained orbit was not physical so we rejected it. After this initial unsuccessful attempt, we tried to fit orbits with additional signals. The following two models did not give satisfactory results: the 1-Keplerian plus linear drift model and the 2-Keplerian one. The best fit we obtained and finally adopted is based on a 1-Keplerian plus quadratic drift model. Despite the
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** HARPS phase-folded velocities and fitted orbits for the 6 targets observed for at least two orbital periods. The radial-velocities of HD 28254 and HD 44219 have been corrected for the fitted quadratic drift.</figcaption></figure>

poor coverage of the second radial-velocity maximum, hence the poorly constrained velocity semi-amplitude, the planetary nature of HD 23254 b is well established. The orbital period is very well constrained. Solutions with significantly larger amplitudes are possible, but they also have much larger eccentricities. The effect of the increasing amplitude on the minimum mass of the companion is largely compensated for by the increasing eccentricity as the mass function scales with $K_1^3(1 - e^2)^{3/2}$.

More data are necessary for better constraining the shape of the orbit of HD 28254 b. In particular, a better coverage of the periastron passage will allow us to constrain the radial-velocity semi-amplitude and orbital eccentricity. These future observations will also help us characterizing the secondary signal better. The orbital solution we present here is thus a preliminary one. The observed additional quadratic signal remains compatible with the signal that the visual companion described in Sect. 2 could induce.

Using our fitted parameters, we computed the minimum mass of the planetary companion and its orbital semi-major axis: $m_2 \sin i = 1.16 M_{\text{Jup}}$ and $a = 2.15$ AU. Our radial-velocity data and the fitted orbit are displayed in panel **c** of Fig. 1. The phase-folded orbit and velocities are displayed in panel **b** of Fig. 2.

### 3.4. HD 290327

Between BJD = 2452945 (November 1, 2003) and BJD = 2454931 (April 9, 2009), we have observed HD 290327 (HIP 25191) 18 times with HARPS. The obtained radial velocities have a mean uncertainty of 1.73 m s⁻¹. A long-period signal is well visible on the data. We list in Table 2 the parameters resulting from a Keplerian fit to the data. With these parameters and using the primary mass of Table 1 (0.90 $M_\odot$), we computed the minimum mass of the companion that induces the detected radial-velocity signal and the semi-major axis of its orbit: $m_2 \sin i = 2.54 M_{\text{Jup}}$ and $a = 3.43$ AU. The dispersion of the residuals to the fitted orbit is very low: 1.6 m s⁻¹. Our velocity data and the fitted orbit are displayed in panel **d** of Fig. 1.

### 3.5. HD 43197

We have gathered a set of 33 radial velocities for HD 43197 (HIP 29550). These data were obtained between BJD = 2452989 (December 15, 2003) and BJD = 2454932 (April 10, 2009), and they have a mean velocity uncertainty of 1.42 m s⁻¹. A visual inspection of the data displayed in panel **e** of Fig. 1 shows a periodic highly eccentric signal. The parameters resulting from a Keplerian fit to the data are listed in Table 2. The fitted eccentricity is indeed very high and very well constrained in spite of our poor coverage of the radial-velocity minimum. The rms of the residuals to the fitted orbit is very low: 1.42 m s⁻¹. This 1-planet fit is thus satisfactory. With a primary mass of 0.96 $M_\odot$, the computed minimum mass of the companion is 0.6 $M_{\text{Jup}}$. The orbital semi-major axis is 0.92 AU. Our data, the fitted orbit, and the residuals are displayed in panel **e** of Fig. 1. Panel **c** of Fig. 2 shows the phase-folded orbit.

### 3.6. HD 44219

The HARPS data set we have in hand for HD 44219 (HIP 30114) consists of 46 radial velocities measured between BJD = 2452944 (October 31, 2003) and BJD = 2454932 (April 10, 2009). These velocities have an average uncertainty of 1.55 m s⁻¹. We first fitted a 1-Keplerian model to these data. The residuals to this fitted orbit were a bit large ($\sigma_{O-C} = 2.95$ m s⁻¹). We then tried to include an additional signal in the fitted model. Adding a linear drift improved the quality of the fit in a marginally significant way ($\chi^2_{\text{red}}$ decreasing from 3.79 to 3.30; f-test probability = 71%). The improvement resulting from including a quadratic trend in the fitted model was more significant: $\chi^2_{\text{red}} = 2.62$ and f-test probability equals 90%. We also tried to fit a 2-planet model but we were unable to find any convincing solution. We finally decided to adopt the Keplerian + quadratic drift model. The parameters resulting from this fit are listed in Table 2. The orbital parameters obtained with our initial 1-planet fit did not differ much from the ones in the adopted solution ($P = 447$ d; $e = 0.57$; $K_1 = 17$ m s⁻¹ and $m_2 \sin i = 0.52 M_{\text{Jup}}$).
The dispersion of the residuals to the finally selected solution is still a bit large ($\sigma_{O-C} = 2.39 \text{ m s}^{-1}$), but this is probably due to the still not very well constrained additional signal. With this solution and assuming $M_* = 1 M_\odot$, the computed minimum mass for HD 44219 b is $0.58 M_{\text{Jup}}$ and its orbital semi-major axis is 1.19 AU. We display the fitted model and the temporal residuals in panel *f* of Fig. 1. The phase-folded orbit and velocities are displayed in panel *d* of Fig. 2.

### 3.7. HD 148156

We have collected 42 radial velocities with HARPS for HD 148156 (HIP 80680) between BJD = 2 452 852 (July 31, 2003) and BJD = 2 455 020 (July 7, 2009). The mean radial-velocity uncertainty of this data set is $2.24 \text{ m s}^{-1}$. A long-period signal is visible in the data displayed in panel *g* of Fig. 1. The results of the Keplerian fit performed on our velocity set are presented in Table 2. The rms of the residuals is a bit large ($\sigma_{O-C} = 3.69 \text{ m s}^{-1}$; $\chi^2_{\text{red}} = 2.89$). We were unable to significantly improve the fit by adding extra signals. No obvious correlation between CCF bisectors and residuals could be detected. The origin of these abnormal residuals remains unknown, but they may be due to intrinsic stellar jitter.

The minimum mass of HD 148156 b is $0.85 M_{\text{Jup}}$ and the semi-major axis of its orbit is 2.45 AU (assuming $M_* = 1.22 M_\odot$). The fitted orbit and the fit residuals are displayed in panel *g* of Fig. 1. The phase-folded orbit and velocities are displayed in panel *e* of Fig. 2.

### 3.8. HD 156411

We obtained 50 HARPS radial-velocity measurements of HD 156411 (HIP 84787). These measurements were taken between BJD = 2 452 850 (July 29, 2003) and BJD = 2 455 081 (September 6, 2009) and they have a mean velocity uncertainty of $1.26 \text{ m s}^{-1}$. A long-period signal is visible on the data displayed in panel *h* of Fig. 1. We fitted a Keplerian orbit to these data. The resulting parameters are listed in Table 2. The dispersion of the residuals to the fitted orbit is too large ($\sigma_{O-C} = 2.94 \text{ m s}^{-1}$), and the fit quality quite poor ($\chi^2_{\text{red}} = 5.21$).

We tried to fit models including additional signals. First, we tried to add linear and quadratic drifts. In both cases, the fit quality improved a bit ($\chi^2_{\text{red}} = 4.80$ and 4.70, respectively) but the corresponding f-test probabilities indicated that the achieved improvements were only marginally significant (at the $1\sigma$ level). A 2-planet fit was also tried, but we were unable to find any convincing solution. We also verified the absence of correlation between CCF bisectors and residuals. HD 156411 is slightly evolved so we speculate that intrinsic stellar signals are responsible for the abnormal residuals.

For this paper, we finally decided to adopt our initial 1-planet orbital solution. With this solution and assuming a mass of $1.25 M_\odot$ for the host star, we computed a minimum mass of $0.74 M_{\text{Jup}}$ and an orbital semi-major axis of $1.88 \text{ AU}$ for HD 156411 b. The fitted orbit and the temporal residuals are displayed in panel *h* of Fig. 1. The phase-folded orbit and velocities are displayed in panel *f* of Fig. 2.

## 4. Conclusions

In this paper, we have presented our HARPS radial-velocity data for 8 low-activity solar-type stars orbited by at least one long-period planetary companion. HD 6718 has a $1.56 M_{\text{Jup}}$ companion on a 2496-d low-eccentricity orbit ($e=0.1$).

The companion orbiting HD 8535 is a $0.68 M_{\text{Jup}}$ planet with a period of 1313 days and an eccentricity of 0.15. HD 28254 b is a $1.16 M_{\text{Jup}}$ planet on a 1116-d highly eccentric orbit ($e=0.81$). The companion to HD 290327 has a minimum mass of $2.54 M_{\text{Jup}}$. Its 2442-d orbit is nearly circular ($e=0.08$). As for HD 28254 b, the $0.6 M_{\text{Jup}}$ companion to HD 43197 is on a highly eccentric orbit: $e=0.83$. Its orbital period is 327 d. HD 44219 is orbited by a $0.58 M_{\text{Jup}}$ planet on a 472-d 0.61 eccentricity orbit. HD 148156 b is a $0.85 M_{\text{Jup}}$ planet orbiting its parent star on a 1027-d eccentric orbit ($e=0.52$) and the companion to HD 156411 is a $0.74 \text{ Jupiter-mass planet with a mildly eccentric}$ ($e=0.22$) $842$-d orbit.

We find some evidence of additional longer-period companions around HD 28254 and HD 44219. In both cases, we are unable to fit full double-Keplerian orbits because of our short observing time span. Future observations of these two systems will probably allow us to better characterize the companions responsible for these signals.

Finally, with the addition of HD 28254 b and HD 43197 b, the population of exoplanets with extremely high eccentricities now amounts to 6. The other planets with $e \ge 0.8$ are HD 4113 ($e=0.903$, Tamuz et al. 2008), 156846 ($e=0.8472$, Tamuz et al. 2008), HD 20782 b ($e=0.92$, Jones et al. 2006), and HD 80606 b ($e=0.934$, Moutou et al. 2009a). Five of these stars hosting highly eccentric planets are members of wide visual binaries (see Tamuz et al. 2008 for HD 4113 and HD 156846, Desidera & Barbieri 2007 for HD 20792, and Naef et al. 2001 for HD 80606). Only HD 43197 is not known to be in a multiple stellar system. This large number of eccentric-planet hosts in binaries tend to favour scenarios in which the extreme orbital eccentricity results from the influence of a third body via Kozai oscillations (Kozai 1962) or Kozai migration (Wu & Murray 2003).

*Acknowledgements.* We thank the Swiss National Science Foundation (SNSF) and Geneva University for their continuous support of our planet search programmes. N.C.S. would like to acknowledge the support from the European Research Council/European Community under the FP7 through a Starting Grant, as well from Fundação para a Ciência e a Tecnologia (FCT), Portugal, through programme Ciência 2007, and in the form of grants PTDC/CTE-AST/098528/2008 and PTDC/CTE-AST/098604/2008. This research has made use of the Simbad database, operated at the CDS in Strasbourg, France.

## References

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

* Bouchy, F., Bazot, M., Santos, N. C., Vauclair, S., & Sosnowska, D. 2005, A&A, 440, 609

* Delfosse, X., Forveille, T., Ségransan, D., et al. 2000, A&A, 364, 217

* Desidera, S., & Barbieri, M. 2007, A&A, 462, 345

* Dommanget, J., & Nys, O. 2002, IAU Commission on Double Stars, 148, 4

* ESA. 1997, The HIPPARCOS and TYCHO catalogue, ESA-SP, 1200

* Fernandes, J., & Santos, N. C. 2004, A&A, 427, 607

* Flower, P. J. 1996, ApJ, 469, 355

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

* Hébrard, G., Udry, S., Lo Curto, G., et al. 2010, A&A, 512, A46

* Jenkins, J. S., Jones, H. R. A., Pavlenko, Y., et al. 2008, A&A, 485, 571

* Jones, H. R. A., Butler, R. P., Tinney, C. G., et al. 2006, MNRAS, 369, 249

* Kozai, Y. 1962, AJ, 67, 591

* Lo Curto, G., Mayor, M., Clausen, J. V., et al. 2006, A&A, 451, 345

* Lo Curto, G., Mayor, M., Benz, W., et al. 2010, A&A, 512, A48

* Mayor, M., Bonfils, X., Forveille, T., et al. 2009, A&A, 507, 487

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

* Mordasini, C., Alibert, Y., Benz, W., & Naef, D. 2009, A&A, 501, 1161

* Mordasini, C., Mayor, M., Udry, S., et al. 2010, A&A, accepted

* Moutou, C., Mayor, M., Bouchy, F., et al. 2005, A&A, 439, 367

* Moutou, C., Hébrard, G., Bouchy, F., et al. 2009a, A&A, 498, L5
Moutou, C., Mayor, M., Lo Curto, G., et al. 2009b, A&A, 496, 513

Moutou, C., Mayor, M., Lo Curto, G., et al. 2010, A&A, submitted

Naef, D., Latham, D. W., Mayor, M., et al. 2001, A&A, 375, L27

Naef, D., Mayor, M., Benz, W., et al. 2007, A&A, 470, 721

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

Pepe, F., Mayor, M., Queloz, D., et al. 2004a, A&A, 423, 385

Pepe, F., Mayor, M., Queloz, D., & Udry, S. 2004b, in *Planetary Systems in the Universe*, ed. A. Penny, IAU Symp., 202, 103

Robinson, S. E., Ammons, S. M., Kretke, K. A., et al. 2007, ApJS, 169, 430

Santos, N. C., Mayor, M., Naef, D., et al. 2000, A&A, 361, 265

Santos, N. C., Mayor, M., Naef, D., et al. 2002, A&A, 392, 215

Santos, N. C., Bouchy, F., Mayor, M., et al. 2004a, A&A, 426, L19

Santos, N. C., Israelian, G., & Mayor, M. 2004b, A&A, 415, 1153

Ségransan, D., Mayor, M., Udry, S., et al. 2010, A&A, submitted

Sousa, S. G., Santos, N. C., Israelian, G., Mayor, M., & Monteiro, M. J. P. F. G. 2007, A&A, 469, 783

Tamuz, O., Ségransan, D., Udry, S., et al. 2008, A&A, 480, L33

Udry, S., Mayor, M., Queloz, D., Naef, D., & Santos, N. 2000, in *From Extrasolar Planets to Cosmology: The VLT Opening Symposium*, ed. J. Bergeron, & A. Renzini, 571

van Leeuwen, F. 2007, A&A, 474, 653

Wu, Y., & Murray, N. 2003, ApJ, 589, 605