Issue 
A&A
Volume 525, January 2011



Article Number  A95  
Number of page(s)  24  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201015427  
Published online  03 December 2010 
Search for browndwarf companions of stars^{⋆,}^{⋆⋆}
^{1}
Observatoire de Genève, Université de Genève,
51 Chemin des Maillettes,
1290
Sauverny,
Switzerland
email: johannes.sahlmann@unige.ch
^{2}
European Southern Observatory, KarlSchwarzschildStr. 2, 85748
Garching bei München,
Germany
^{3}
Centro de Astrofísica, Universidade do Porto,
Rua das Estrelas,
4150762
Porto,
Portugal
^{4} Departamento de Física e Astronomia, Faculdade de Ciências,
Universidade do Porto, Portugal
^{5}
Department of Geophysics and Planetary Sciences, Tel Aviv
University, Tel Aviv
69978,
Israel
Received:
19
July
2010
Accepted:
23
September
2010
Context. The frequency of browndwarf companions in close orbit around Sunlike stars is low compared to the frequency of planetary and stellar companions. There is presently no comprehensive explanation of this lack of browndwarf companions.
Aims. By combining the orbital solutions obtained from stellar radialvelocity curves and Hipparcos astrometric measurements, we attempt to determine the orbit inclinations and therefore the masses of the orbiting companions. By determining the masses of potential browndwarf companions, we improve our knowledge of the companion massfunction.
Methods. The radialvelocity solutions revealing potential browndwarf companions are obtained for stars from the CORALIE and HARPS planetsearch surveys or from the literature. The best Keplerian fit to our radialvelocity measurements is found using the LevenbergMarquardt method. The spectroscopic elements of the radialvelocity solution constrain the fit to the intermediate astrometric data of the new Hipparcos reduction. The astrometric solution and the orbit inclination are found using nonlinear χ^{2}minimisation on a twoparameter search grid. The statistical confidence of the adopted orbital solution is evaluated based on the distributionfree permutation test.
Results. The discovery of nine new browndwarf candidates orbiting stars in the CORALIE and HARPS radialvelocity surveys is reported. New CORALIE radial velocities yielding accurate orbits of six previouslyknown hosts of potential browndwarf companions are presented. Including the literature targets, 33 hosts of potential browndwarf companions are examined. Employing innovative methods, we use the new reduction of the Hipparcos data to fully characterise the astrometric orbits of six objects, revealing Mdwarf companions of masses between 90 M_{J} and 0.52 M_{⊙}. In addition, the masses of two companions can be restricted to the stellar domain. The companion to HD 137510 is found to be a brown dwarf. At 95% confidence, the companion of HD 190228 is also a brown dwarf. Twentythree companions remain browndwarf candidates. On the basis of the CORALIE planetsearch sample, we obtain an upper limit of 0.6% for the frequency of browndwarf companions around Sunlike stars. We find that the companionmass distribution function increases toward the lower end of the browndwarf mass range, suggesting that we detect the highmass tail of the planetary distribution.
Conclusions. Our findings agree with the results of previous similar studies and confirm the pronounced paucity of browndwarf companions around Sunlike stars. They are affected by the Hipparcos astrometric precision and mission duration, which limits the minimum detectable companion mass, and some of the remaining candidates are probably browndwarf companions.
Key words: binaries: spectroscopic / brown dwarfs / stars: lowmass / techniques: radial velocities / astrometry
Based on observations made with the CORALIE spectrograph on the Swiss telescope at the European Southern Observatory in La Silla, Chile, with the ESA Hipparcos astrometry satellite, and with the HARPS instrument on the ESO 3.6m telescope (GTO programme 072.C0488).
The CORALIE and HARPS radialvelocity measurements discussed in this paper are only available in electronic form at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/vol/pge
© ESO, 2010
1. Introduction
The mass distribution of close companions to Sunlike stars exhibits a clear deficit of brown dwarfs compared to planets and stellar binaries (Grether & Lineweaver 2006). With a frequency of less than 1% (Marcy & Butler 2000), browndwarf companions to solartype stars at separations below 10 AU are much less common than planetary companions with a frequency of ~7% (e.g. Udry & Santos 2007) and stellar binaries with a frequency of ~13% (Duquennoy & Mayor 1991; Halbwachs et al. 2003) at similar separations. Brown dwarfs are substellar objects, massive enough to burn deuterium but too light to permit hydrogen burning, which places them in the mass range of approximately 13–80 Jupiter masses (Burrows et al. 1997; Chabrier & Baraffe 2000; Spiegel et al. 2010). They therefore constitute intermediate objects between planets and stars, though the lower end of the browndwarf mass range overlaps with the one of massive planets (e.g., Luhman et al. 2009), and distinguishing planets and brown dwarfs may require us to trace the individual formation processes. A comprehensive description of the formation of lowmass objects, covering the range of planets to very lowmass stars, remains undefined, although the progress in this direction is formidable as reported in the reviews of Whitworth et al. (2007), Luhman et al. (2007), and Burgasser et al. (2007).
A few companions with masses of 18–60 Jupiter masses (M_{J}) have been identified by astrometry (e.g., Martioli et al. 2010; Benedict et al. 2010) or discovered in transiting systems (e.g., Deleuil et al. 2008; Bouchy et al. 2010; Anderson et al., in preparation). Yet, most potential brown dwarf companions are discovered in radial velocity surveys (e.g. Nidever et al. 2002). In these highprecision programs, designed to achieve planet detection, they are easily spotted and characterised because of their large signatures, by far exceeding the radial velocity measurement precision (e.g., Patel et al. 2007; Bouchy et al. 2009). However, radialvelocity measurements alone do not constrain the orbit inclination. Therefore, they cannot reveal the companion mass, but determine its lower limit. Observations of complementary effects such as the transit lightcurve (Deleuil et al. 2008) or the astrometric motion of the host star are required to solve this ambiguity and determine the companion mass.
Highprecision astrometry with the Hubble Space Telescope (HST) has been used to find the masses of planetary companions discovered by radial velocity (Benedict et al. 2010, and references therein) and has demonstrated the capability of astrometry to distinguish between planetary, browndwarf, or stellar companions. Groundbased optical interferometers achieve the necessary precision to detect orbital motions of binary and multiple systems and to determine the orbital parameters in combination with radialvelocity measurements (Hummel et al. 1993; Muterspaugh et al. 2006; Lane et al. 2007). The astrometric detection of planetary orbits using IRinterferometry is the aim of the extrasolarplanet search with PRIMA project (ESPRI, Launhardt et al. 2008). On a larger scale, the GAIA satellite (e.g. Lindegren 2010) will advance the field of astrometric detection and characterisation of planetary systems, because of its outstanding measurement precision (Casertano et al. 2008).
Hipparcos astrometry has been extensively used to constrain the masses of substellar and stellar companions in multiple systems (Perryman et al. 1996; Pourbaix & Jorissen 2000; Zucker & Mazeh 2001; Torres 2007). Han et al. (2001) used Hipparcos astrometry to analyse 30 extrasolarplanet candidates and concluded that half of these companions were rather brown or M dwarfs. With the consequent rectification by Pourbaix & Arenou (2001), showing that the results of Han et al. (2001) had been biased by the fitting procedure, it became clear that the fitting of astrometric data has to be performed very carefully, especially when the size of the orbital signature approaches the measurement precision of the instrument (Pourbaix 2001). All of these past works used the original Hipparcos data (Perryman et al. 1997).
In 2007, the new reduction of the raw Hipparcos data (van Leeuwen 2007) was released, which represents a significant improvement in quality compared to the original data (van Leeuwen 2007). It also includes a revised version of the intermediate astrometric data in a new format that facilitates the search for the signature of orbital motion. We use the new Hipparcos reduction for our study, which aims to determine the masses of companions to stars detected using highprecision radialvelocity measurements and found to have minimum masses in the browndwarf domain. Our target sample is composed of stars from the CORALIE and HARPS planet surveys. In addition, we included a list of targets selected from the literature. To validate our method of astrometric orbit determination, we also analysed a few comparison targets with published astrometric orbits. This study explores the brown dwarf desert and is in line with the foregoing works by Halbwachs et al. (2000) and Zucker & Mazeh (2001).
The paper is organised as follows. The browndwarf candidates from the CORALIE and HARPS surveys are presented in Sect. 2 with an additional 18 candidates selected from the literature. The method of combining radialvelocity orbits and Hipparcos astrometry is described in Sect. 3. The results are summarised in Sect. 4 and discussed in Sect. 5. We conclude in Sect. 6.
2. Brown dwarf candidates
Our target sample consists of 33 stars that exhibit radial velocity variations caused by a companion with minimum mass M_{2} sin i in the browndwarf mass range of 13–80 M_{J}. It contains 14 stars from the CORALIE survey (Udry et al. 2000) and one star from the HARPS planet search program (Mayor et al. 2003). To date, the CORALIE planet survey has contributed to the discovery of more than 50 extrasolar planets (e.g. Naef et al. 2001; Mayor et al. 2004; Tamuz et al. 2008). CORALIE is an optical echelle spectrograph mounted on the 1.2 m Swiss Telescope located at the European Southern Observatory in La Silla, Chile. A description of the instrument can be found in Queloz et al. (2000) and the references therein. With its unmatched precision, the HARPS instrument permitted the discovery of more than 70 planetary companions to date, including earthmass planets (e.g. Pepe et al. 2007; Santos et al. 2007; Mayor et al. 2009; Forveille et al. 2009). The characteristics, radial velocities, and orbital solutions of the 15 CORALIE and HARPS stars are presented and discussed in Sects. 2.1–2.3. In addition, 18 potential browndwarf host stars with radialvelocity orbits were selected from the literature and are listed in Sect. 2.4. These include HD 190228, initially a planethost star (Sivan et al. 2004), whose importance we recognised during the analysis.
Hipparcos parameters of the surveyed stars.
Derived stellar parameters of the surveyed stars.
2.1. Stellar characteristics
The identifiers and basic stellar characteristics of the CORALIE stars and HIP 103019 from the HARPS survey are listed in Table 1. The apparent visual magnitude V, the colour B − V , and the parallax ϖ are from the new Hipparcos reduction (van Leeuwen 2007), whereas the spectral type is from the original Hipparcos catalogue (Perryman et al. 1997). Table 2 displays the derived stellar parameters. The uncertainty in V is smaller than 0.002 mag. The absolute magnitude M_{V} is derived from the Hipparcos magnitude and parallax. The effective temperature T_{eff}, the surface gravity log g, and the metallicity [Fe/H] are derived from the spectroscopic analysis of high signaltonoise ratio (S/N) spectra with the method presented in Santos et al. (2004) using the FeI and FeII lines listed in Sousa et al. (2008). The stellar rotation parameter ν sin i is derived from the calibration of the CORALIE or HARPS crosscorrelation function by Santos et al. (2002). Finally, the stellar mass of the primary M_{1} and the age are estimated from the theoretical isochrones of Girardi et al. (2000) and a Bayesian estimation method described in da Silva et al. (2006)^{1}. Stellar ages of mainsequence dwarfs are usually not well constrained and we quote the 1σ confidence interval. The spectroscopic analysis was problematic for HIP 103019, possibly because of its faintness and late spectral type. Errors in the parameters of this star are possibly underestimated. The stellar characteristics of the targets from the literature can be found in the respective references given in Sect. 2.4.
2.2. Radialvelocity measurements and orbital solutions
Optical highresolution spectra of the stars in Table 2 were collected with the CORALIE and HARPS spectrographs over time intervals extending to 11 years. Radial velocities were estimated by means of the crosscorrelations of the extracted stellar spectra with numerical templates, which depend on the target’s spectral type. During CORALIE observations, a reference thoriumargon spectrum is recorded simultaneously with the stellar spectrum and used to measure and correct for residual zeropoint drifts (Baranne et al. 1996). Photon noise limits the obtained precision to typically 2–4 m s^{1} per epoch. To account for systematic drifts, an external error of 5 m s^{1} is quadratically added to CORALIE radialvelocity uncertainties before performing the period search and the model adjustment. The abbreviations C98 and C07 refer to the CORALIE instrument before and after its upgrade in 2007, respectively (Ségransan et al. 2010). The two parameters γ_{C98} and γ_{C07} are introduced to account for the possibly different velocity offsets of these instrument configurations.
Because of the stability of the HARPS instrument and the targeted precision of 2 m s^{1} per epoch for the volumelimited programme, HIP 103019 was observed without thoriumargon reference (cf. Naef et al. 2010). To account for possible drifts, an additional error of 0.5 m s^{1} is quadratically added to the HARPS radialvelocity uncertainties.
The Keplerianorbit solution is found by imposing the model function Eq. (11) with 6 parameters (K_{1}, e, ω, T, P, and γ) to the radial velocities of a given star. The fit determines the values of the period P, the eccentricity e, the time of periastron passage T_{0}, the longitude of periastron ω, the radialvelocity semiamplitude K_{1}, and the systemic velocity γ. In general, measurements with uncertainties larger than 15 m s^{1} are not used in the model adjustments.
A singlecompanion model closely reproduces the data of 11 stars with reducedchi values of χ_{r} = 0.7–1.3. For 4 stars, HD 17289, HD 89707, HD 211847, and HIP 103019, the fit is not as good with χ_{r} = 1.7–2.0 and produces residuals of unusually large amplitude, without contesting the presence of a large orbital motion. Introducing additional model parameters, corresponding to a linear or quadratic drift or a second companion does not significantly improve the fit. Therefore, the moderate fit quality has to be attributed to either underestimated radial velocity errors or to additional noise sources, such as stellar activity. This is discussed for the individual targets in the next section.
The orbital periods of HD 4747 and HD 211847 are very long and not covered by the measurement timespan. Hence, the bestfit period is poorly constrained and realistic confidence intervals have to be derived by a different method. We chose a criterion based on the fitquality metric χ_{r} and its error σ_{χ} obtained for the best fit. During a second fit, we fix the period during the adjustment, such that the resulting χvalue increases to χ_{r} + 3σ_{χ}. In this way, we obtain upper and lower limits to the orbital parameters, which are listed at the bottom of Tables 3 and 4. The method that searches for the astrometric signature in the Hipparcos data presented in Sect. 3 relies on the precise knowledge of the radialvelocity parameters. Because the periods of HD 4747 and HD 211847 are poorly constrained and the orbital coverage of Hipparcos is very short, we do not perform the combined astrometric analysis for these targets. They are presented as part of the browndwarf candidates characterised with CORALIE.
Table 3 displays the orbital elements of the bestfit solutions and Table 4 lists the basic statistics and the derived orbit parameters (see also Fig. 1). The parameter T_{0} is given as reduced barycentric Julian date (BJD), that is JD^{⋆} = BJD − 2 400 000. In addition, N_{RV} is the total number of measurements, ΔT is the observation timespan, S/N is the median S/N per pixel at 550 nm, ⟨ σ_{RV} ⟩ is the mean measurement uncertainty, χ_{r} is the reducedchi value of the fit, G.o.F. is the goodness of fit, σ_{(O−C)} gives the rootmeansquare dispersion of the fit residuals, a_{1} sin i is the minimum astrometric semimajor axis expressed in unit of length, f(m) is the mass function, a_{rel} is the semimajor axis of the relative orbit, and finally M_{2} sin i is the minimum companion mass. The quoted errors correspond to 1σ confidence intervals derived from 5000 Monte Carlo simulations. Figures 4 and 5 show the phasefolded radial velocities and the adopted orbital solution. All measurements are available in electronic form at the CDS as indicated on the title page.
Fig. 1 Orbital characteristics of the 15 browndwarf candidates observed with CORALIE and HARPS. Filled symbols indicate the companions with stellar mass identified by this study. The error bars are often smaller than the symbol size, though the orbits of HD 4747 and HD 211847 are not covered and therefore less constrained, see the text for details. Top: minimum companion mass and separation. Bottom: eccentricity and orbital period. 
Fig. 2 Radial velocities (top) and residuals (bottom) to the best fit, represented by the solid line, of HD 4747 measured with HIRES (orange diamonds), CORALIE C98 (red circles), and CORALIE C07 (blue squares). 
Singlecompanion Keplerian orbital solutions of the surveyed stars.
Companion parameters for the surveyed stars.
2.3. Details of the surveyed stars
We now present the 15 stars observed with CORALIE and HARPS, and discuss details about their respective radialvelocity solutions:

HD 3277 (HIP 2790) is found tohost a potential browndwarf companion. The star was listed as anuncertain spectroscopic binary in Tokovininet al. (2006), who did not havesufficient radialvelocity measurements to constrain the orbit.

HD 4747 (HIP 3850) has a potential browndwarf companion discovered with radial velocities by Nidever et al. (2002). We combine CORALIE observations with the KECK HIRES radial velocities of Nidever et al. (2002) retrieved from the SB9catalogue (Pourbaix et al. 2004). However, the orbital period is very long and not covered by the measurements (see Fig. 2). Tables 3 and 4 show the bestfit solution and the solutions with fixed period such that the χvalue increases to χ_{r} + 3σ_{χ} = 1.56. The resulting period range is P = 7900–29 000 days with corresponding minimum companionmass M_{2} sin i = 39.6–58.1 M_{J}.

HD 17289 (HIP 12726) is found to host a potential browndwarf companion. Goldin & Makarov (2007) independently found an orbital signature from Hipparcos astrometry alone with a period of 536 ± 12 days, which is close to our solution (562.1 ± 0.3 days). However, the orbital parameters derived by Goldin & Makarov (2007) are less precise and less accurate than the radialvelocity solution. The dispersion in the fit residuals is abnormally large for CORALIE. The detailed analysis of the crosscorrelation function (CCF) indicates that its bisector velocity span and fullwidthathalfmaximum (FWHM) are correlated with the orbital phase (Fig. 6). As we show in Sect. 4, the companion of HD 17289 is a star with mass of 0.52 M_{⊙}. The two stars are very close (<50 mas, Sect. 4.1) and their estimated intensity ratio in the visible is about 1:180. Therefore, the companion’s light is also injected in the science fibre of CORALIE and its signature appears in the collected spectra. The spectral lines of the companion distort the CCF depending on the orbital phase. As expected, the CCF width is minimal when the radialvelocity curves of both components cross and equal the systemic velocity. The detailed modelling of this doublelined spectroscopic binary is beyond the scope of this paper, but we have identified the binary as the cause of a phasedependent radialvelocity bias, which explains the excess noise in our measurements.

HD 30501 (HIP 22122) is found to host a potential browndwarf companion.

HD 43848 (HIP 29804) has a potential browndwarf companion discovered with radial velocity by Minniti et al. (2009). An astrometric analysis using Hipparcos and the solution of Minniti et al. (2009) was performed by Sozzetti & Desidera (2010). We combined the published measurements^{2} with the CORALIE velocities to adjust the Keplerian model. While the orbital period is well defined, the distributions of e, K_{1}, T_{0}, and ω obtained from the Monte Carlo simulations are bimodal. This is because the lower velocity turnover is not well sampled by the measurements. The fitquality in terms of χ^{2} is comparable for all Monte Carlo simulations and cannot be used to identify the correct solution. In Fig. 7, we show the eccentricity and semiamplitude distributions. Because the realisation of the lowereccentricity orbit is more probable (77%) than the higheccentricity orbit (23%), we consider only the orbits with e ≲ 0.78 to find the final solution and perform the astrometric analysis. For completeness, we also list the less probable and highereccentricity solution in Tables 3 and 4.

HD 52756 (HIP 33736) is found to host a potential browndwarf companion.

HD 53680 (HIP 34052) is found to host a potential browndwarf companion.

HD 74014 (HIP 42634) has a potential browndwarf companion discovered with radial velocities by Patel et al. (2007). The original orbit is compatible with our solution obtained from 118 CORALIE velocities covering the complete orbit.

HD 89707 (HIP 50671) has a potential browndwarf companion announced by Duquennoy & Mayor (1991) and analysed by Halbwachs et al. (2000) based on radial velocities collected with the CORAVEL and ELODIE instruments. The original orbit is compatible with our solution obtained from 37 CORALIE velocities. A combined astrometric analysis using Hipparcos was performed by Zucker & Mazeh (2001). We find slight evidence of reemission at the bottom of the CaII H absorption line (Fig. 8) and derive an activity index of from a high S/N CORALIE spectrum. Using the relation of Santos et al. (2000), we expect an additional activityinduced radialvelocity noise of 12 m s^{1}, which explains the unusuallylarge fit residuals encountered for HD 89707. Accounting for the additional noise during the adjustment yields an acceptable fit quality with χ_{r} = 0.91 ± 0.13.

HD 154697 (HIP 83770) is found to host a potential browndwarf companion.

HD 164427A (HIP 88531) has a potential browndwarf companion announced by Tinney et al. (2001), whose original orbit is compatible with our solution. A combined astrometric analysis using Hipparcos was performed by Zucker & Mazeh (2001).

HD 189310 (HIP 99634) is found to host a potential browndwarf companion in a shortperiod orbit of 14.2 days.

HD 167665 (HIP 89620) has a potential browndwarf companion announced by Patel et al. (2007), whose original orbit is compatible with our solution.
Fig. 3 Radial velocities (top) and residuals (bottom) to the best fit, represented by the solid line, of HD 211847 measured with CORALIE C98 (red circles) and CORALIE C07 (blue squares).
Fig. 4 Phasefolded radial velocities of 8 stars with potential browndwarf companions. Red circles and blue squares indicate measurements with CORALIE C98 and C07, respectively. For HD 43848, the orange diamonds show the MIKE measurements. The solid lines correspond to the bestfit solutions. The error bars are smaller than the symbol size.
Fig. 5 Phasefolded radial velocities of 5 stars with potential browndwarf companions. Error bars are smaller than the symbol size.

HD 211847 (HIP 110340) is found to host a potential browndwarf companion. The longperiod orbit is not completely covered by our measurements (Fig. 3) and the CORALIE upgrade occurred at minimum radial velocity. To obtain a reasonable solution, we therefore fix the offset between the parameters γ_{C98} and γ_{C07} to zero, which is the mean offset value found for three other stars of the same spectral type. As for HD 4747, we derive realistic confidence intervals by fixing the period during the adjustment, such that the resulting χvalue increases to χ_{r} + 3σ_{χ} = 2.32. The period lower limit is P = 3750 days and we set the upper limit to P = 100 000 days, because we cannot derive an upper limit using the defined criterion. The corresponding minimum companion mass is M_{2} = 17.0–23.3 M_{J}. The dispersion in the bestfit residuals is abnormally large for CORALIE. We find evidence of reemission at the bottom of CaII H and CaII K (see Fig. 9) and derive an activity index of from a high S/N CORALIE spectrum. Using the relation of Santos et al. (2000), we expect an additional activityinduced radialvelocity noise of 11 m s^{1}, which explains the fit residuals encountered for HD 211847. Accounting for the additional noise during the adjustment yields an acceptable fit quality with χ_{r} = 1.15 ± 0.15.

HIP 103019 is discovered to host a potential browndwarf companion using HARPS. The star is included in the search for planets in a volumelimited sample (e.g. Naef et al. 2010). Twentysix radial velocities were collected with HARPS with a mean uncertainty of 2.3 m s^{1}. The dispersion in the fit residuals is abnormally large for HARPS (σ_{O−C} = 3.57 m s^{1}). Because of its faintness, we were not able to measure the activity level of HIP 103019. We therefore cannot exclude an additional radialvelocity jitter induced by the presence of spots on the surface of the star. The fit is satisfactory (χ_{r} = 1.04) if we include additional noise of 3.5 m s^{1}.
2.4. Browndwarf candidates from the literature
In addition to the stars from the CORALIE and HARPS surveys, we selected 18 stars from the literature using the following criteria. They have announced radialvelocity companions with minimum mass of 13–80 M_{J}, were published after 2002, and are not flagged as binary stars in the new Hipparcos reduction. The selected stars are listed in Table 5 and details are given below:
List of targets selected from the literature.

HD 13189 is a K2 II star with a very uncertain massof 2–7 M_{⊙} and a companion announced by Hatzes et al. (2005) with M_{2} sin i = 8–20 M_{J}.

HD 30339 (F8, 1.1 M_{⊙}), HD 65430 (K0V, 0.78 M_{⊙}), HD 140913 (G0V, 0.98 M_{⊙}), and HD 174457 (F8, 1.19 M_{⊙}) have companions announced by Nidever et al. (2002) with M_{2} sin i ranging between 42–78 M_{J}.

HD 38529 is a G4 IV star of mass 1.48 M_{⊙}. Companion announced by Fischer et al. (2003) with M_{2} sin i = 12.7 M_{J}. Refined orbits are given by Wright et al. (2009) and Benedict et al. (2010). We use the latest orbital elements for analysis.

HD 91669 is a K0 dwarf of mass 0.88 ± 0.04 M_{⊙}. Companion announced by Wittenmyer et al. (2009) with M_{2} sin i = 30.6 M_{J}.

HD 107383 is a G8 III giant star of mass 2.7 ± 0.3 M_{⊙} and a companion announced by Liu et al. (2008) with M_{2} sin i = 19.4 M_{J}.

HD 119445 is a G6 III giant star of mass 3.9 ± 0.4 M_{⊙} and a companion announced by Omiya et al. (2009) with M_{2} sin i = 37.6 M_{J}.

HD 131664 is a G3V star of mass 1.10 ± 0.03 M_{⊙}. Companion announced by Moutou et al. (2009) with M_{2} sin i = 18.15 M_{J}.

HD 137510 is a G0 IV star of mass 1.42 M_{⊙} and a companion announced by Endl et al. (2004). The updated data by Butler et al. (2006) lists M_{2} sin i = 22.7 M_{J}.

GJ 595 is a M3 dwarf of mass 0.28 M_{⊙}. Companion announced by Nidever et al. (2002) with M_{2} sin i = 60 M_{J}.

HD 162020 is a K3 V star of mass 0.75 M_{⊙} and a companion announced by Udry et al. (2002) with M_{2} sin i = 14.4 M_{J}.

HD 168443 has two companions with minimum masses in the planetary (HD 168443b) and brown dwarf domain (HD 168443c, M_{2} sin i = 18 M_{J}) discovered by Marcy et al. (2001) and Udry et al. (2004; see also Udry et al. 2002). Updated orbital elements are given by Wright et al. (2009), which we use for the astrometric analysis for the outer companion HD 168443c.

HD 180777 is a A9V star of mass 1.7 ± 0.1 M_{⊙} and a companion announced by Galland et al. (2006) with M_{2} sin i = 25 M_{J}.

HD 190228 is a G5 IV star of mass 0.83 M_{⊙} and a planetarymass companion with M_{2} sin i = 3.58 M_{J} (Perrier et al. 2003). Zucker & Mazeh (2001) performed an astrometric analysis using Hipparcos and derived M_{2} = 67 ± 29 M_{J} at 95% confidence based on the spectroscopic elements of Sivan et al. (2004; see also Sect. 3.8.2).

HD 191760 is a G3 IV star of mass 1.26 ± 0.06 M_{⊙}. Companion announced by Jenkins et al. (2009) with M_{2} sin i = 38.17 M_{J}.

HD 202206 is a G6 V star of mass 1.15 M_{⊙} and has two companions with minimum masses in the planetary (HD 202206c) and brown dwarf domain (HD 202206b, M_{2} sin i = 17.4 M_{J}) discovered by Udry et al. (2002) and Correia et al. (2005).
We identified our target sample of stars with potential browndwarf companions and obtained the orbital parameters accessible from radialvelocity measurements. As the next step, we search the intermediate astrometric data of the new Hipparcos reduction for signatures of the orbital motion corresponding to the radialvelocity orbit.
Fig. 7 Eccentricities (left) and velocity semiamplitudes (right) of 10 000 Monte Carlo simulations performed to find the radialvelocity orbit of HD 43848. The bimodal structure is clearly seen in both histograms. The final solution is found by considering only the lower eccentricity realisations, which is done by separating the distributions at m s^{1} and correspondingly at ê ~ 0.78. The probability of obtaining a high or loweccentricity solution is calculated by counting the number of realisations above and below this threshold value and equals 77% and 23%, respectively. 
3. Combined astrometric analysis
Fig. 8 CaII H emission region for HD 89707 showing slight reemission at the bottom of this absorption line. 
Fig. 9 CaII H emission region for HD 211847 showing reemission at the bottom of the absorption line. 
We describe a method to search the intermediate astrometric data (IAD) of the new Hipparcos reduction for the orbital signatures of stars, whose spectroscopic elements are known from a reliable radialvelocity solution. The significance of the derived orbit is inferred from the distributionfree permutation test (Good 1994; Zucker & Mazeh 2001).
Techniques to perform the simultaneous fitting to radialvelocity and Hipparcos astrometry data for substellar companions and verify the statistical confidence of the solution have been developed by Zucker & Mazeh (2000), Pourbaix (2001), Halbwachs et al. (2000), and Zucker & Mazeh (2001). Reffert & Quirrenbach (2006) and Sozzetti & Desidera (2010) used similar approaches to claim the detection of three browndwarf companions, although with a significance not exceeding 2σ. Independently, successful detections of the astrometric signature of exoplanet candidates using the HST fine guidance sensor (HST FGS) optical interferometer are reported for a few stars with spectroscopic elements known from radial velocity (e.g. Bean et al. 2007; Martioli et al. 2010). Even though the orbital phase is not always fully covered by HST observations, for instance the orbital coverage for HD 33636 is 20% (Bean et al. 2007), their superior precision makes detection possible.
Our analysis is applied to stars with an orbital solution obtained from radial velocities that provides the spectroscopic elements P, e, T_{0}, ω, and K_{1}. The system is modelled as a primary mass M_{1} orbited by one invisible companion of mass M_{2}. Using Hipparcos astrometric measurements, we seek to determine the two remaining unknown parameters that characterise the orbit, which are the inclination i and the longitude of the ascending node Ω. The analysis is performed in three steps. First, the Hipparcos abscissa, possibly containing the companion signature, is constructed from the catalogue data. This is important, because the IAD contains only abscissa residuals, which are the astrometric residuals after subtraction of the model adopted by the new Hipparcos reduction. Second, the bestfit values for i and Ω are found by determining the global χ^{2} minimum on a twodimensional search grid. Finally, the significance of the obtained solution is evaluated.
3.1. Constructing the Hipparcos abscissa
On the basis of the Hipparcos identifier of a given star, the astrometric data is retrieved from the new reduction catalogue. It includes the astrometric parameters of the solution, the solution type, and the goodness of fit. The IAD is read from the resrec folder in the catalogue DVD of van Leeuwen (2007) and contains the satellite orbit number, the epoch t, the parallax factor Π, the scan angle orientation ψ, the abscissa residual δΛ, and the abscissa error σ_{Λ} for every satellite scan.
For the stars in our sample, we encounter the solution types “1”, “5”, “7”, and “9”. The solution type “5” indicates that the five standard astrometric parameters fit the data reasonably well, whereas a solution type “1” is termed a stochastic solution and is adopted when the fiveparameter fit is not satisfactory and neither orbital nor acceleration models improve the solution in terms of χ^{2}. The solution types “7” and “9” are given, when the solution has to include acceleration parameters (propermotion derivatives of first and second order) to obtain a reasonable fit. The solution types “1”, “7”, and “9” may be the first indications of a possible astrometric perturbation by an unseen companion. The Hipparcos abscissa Λ_{HIP} is then given by: (1)where the modified right ascension is given by α^{⋆} = αcosδ, the declination δ, the associated proper motions μ_{α⋆} and μ_{δ}, and the parallax ϖ. The epoch t is given as time in years since t_{HIP} = 1991.25 A.D. = 2 448 348.8125 JD and the correction factors for solution types “7” and “9” are given by (2)(3)where p is either α^{⋆} or δ (van Leeuwen 2007), and and are the acceleration and the change in acceleration in the coordinate direction p, respectively. For solution types “1” and “5”, Δ_{7}p and Δ_{9}p vanish and the five parameter model holds.
Because of the linear nature of Eq. (1), it is usually convenient to decompose each of the five standard parameters into a constant term and an offset, e.g. , and to set the constant terms to zero to compute the abscissa. As explained in Sect. 3.4, however, our formalism requires the value of the parallax ϖ to be present in Eq. (1), whereas we can safely set for the calculation of the abscissa.
3.2. Model function
The model function describes the astrometric signature along the scan axis ψ of an isolated star, orbited by one invisible companion. To develop our approach, we list the standard formulae (e.g. Hilditch 2001) applicable to this system. The astrometric motion of the star is characterised by P, e, T_{0}, i, ω, Ω, and the astrometric semimajor axis a, expressed in angular units. We write the ThieleInnes constants as (4)where c_{A,B,F,G} are functions of i, Ω, and ω and a is expressed in milliarcseconds (mas). The elliptical rectangular coordinates X and Y are functions of eccentric anomaly E and eccentricity If the singlecompanion assumption is true, the Hipparcos abscissa is described by a linear equation with 9 parameters, which are the 5 standard astrometric parameters plus the 4 ThieleInnes constants (8)The relationship between astrometrically and spectroscopically derived parameters is (Pourbaix 2001) (9)where K_{1}, P, and ϖ are expressed in m s^{1}, years, and mas, respectively. Hence, the astrometric model function can also be written as a nonlinear function of 12 independent parameters (10)The observed radial velocity v_{rad} of this orbit is given by (11)where γ is the systemic velocity. The true anomaly θ and the semiamplitude K_{1} are constrained by (12)(13)
3.3. Definition of the search grid
One technique to find a possible orbital signature consists of defining a twodimensional search grid in inclination i and ascending node Ω to solve for the remaining parameters of the model function (Zucker & Mazeh 2000; Reffert & Quirrenbach 2006). In this way, the goodnessoffit is obtained for each point on the grid and the optimal solution can be identified as the global χ^{2} minimum. The search grid is defined by i = 0–180° and Ω = 0–360°, to avoid negative semimajor axes and sample the complete range of values of the cosine function (cf. Eqs. (4) and (9)). We use a grid size of 30 × 30, which provides a reasonable compromise between computation time and grid resolution.
3.4. Joint confidence intervals
For each point on the iΩgrid, we find the least squares solution by solving a linear equation using matrix inversion, where e, P, T_{0}, ω, K_{1}, i, and Ω are fixed. The χ^{2} value hence obtained for each gridpoint allows us to derive joint confidence intervals on the iΩgrid. The parallax dependence of the ThieleInnes constants is removed by writing (14)and defining the new constant (15)We note that a_{1} is expressed in AU and corresponds to the angular semimajor axis a of the astrometric orbit. We can now rewrite Eq. (8) as (16)This relation is linear in the five remaining free parameters (α^{⋆}, δ, ϖ, μ_{α⋆}, μ_{δ}) and can easily be solved analytically. The companion signature is solely contained in an additive modulation Υ(i,Ω,e,P,T_{0},ω,K_{1}) of the parallax factor. It is important to realise that the presence of Υ forces the parallax (instead of a parallax offset) to be present in the model function, which originates in the abscissa reconstruction from Eq. (1). For the remaining parameters, it is sufficient to consider the offsets Δα^{⋆}, Δδ, Δμ_{α⋆}, and Δμ_{δ} to the catalogue values because the abscissa was constructed with .
For every combination of i and Ω, the linear equation Eq. (16) is solved, yielding the five astrometric parameters and the corresponding χ^{2}. The bestfit parameters are identified by the minimum χ^{2} value on the iΩgrid. The results of the linear adjustment are not directly used in the quoted final solution, but they serve as a consistency check of the results from the nonlinear fitting and are used for graphical illustration of the joint confidence intervals shown in Fig. 12. Lowsignificance orbits can show several local χ^{2}minima corresponding to approximately opposite orbit orientations (see also e.g. Zucker & Mazeh 2001; Reffert & Quirrenbach 2006). We find that moderate and highsignificance orbits have one global χ^{2}minimum and the confidence contours cover a small area of the iΩspace (similarly to HD 53680 in Fig. 12). The solution parameters corresponding to an oppositeorientation orbit are always beyond the 4σ contour, i.e. all orbital parameters of significant solutions are unambiguously determined.
3.5. Determining the final solution
While the solution method discussed in the previous section is favourable because the model function is linear, its accuracy is limited by the resolution of the iΩgrid. For instance, the resolution of a square grid with 900 points is 6° and 12° in i and Ω, respectively. To avoid this limitation, we determine the optimal solution by means of χ^{2}minimisation of a model function, where i and Ω are free parameters and thus can adopt continuous values. For this purpose, the iΩgrid is used to define the starting values of a nonlinear least squares fit by the LevenbergMarquardt method. The model function is given by Eq. (10), but is limited to seven free parameters: the five astrometric parameters α^{⋆}, δ, ϖ, μ_{α⋆}, μ_{δ}, plus i and Ω. We select the nonlinear solution yielding the smallest χ^{2} and perform 1000 Monte Carlo simulations. Every Monte Carlo realisation consists of generating a set of Hipparcos abscissa measurements and consequent χ^{2}minimisation of the nonlinear model and therefore provides 1000 sets of solution elements.
To incorporate the uncertainties in the spectroscopic parameters, we perform the complete analysis of 100 sets of spectroscopic parameters, where each set is randomly drawn from a Gaussian distributions with mean and standard deviation given by the radialvelocity solution and its error, respectively (as in Sozzetti & Desidera 2010). The final solution accounts for all 1000 Monte Carlo solutions, which are obtained for each of the 100 sets of spectroscopic parameters. To achieve this, we combine all solutions to yield distributions of 100 000 values for each parameter in Eq. (10). The final parameter value and its error is derived from the mean and the 1σ confidence interval of the associated distribution, respectively. The values and errors of derived quantities, such as the companion mass, are obtained in the same way, thus take into account possible correlations between individual parameters.
3.6. Outliers in the astrometric data
Outlying astrometric measurements are examined if they deviate by more than 4σ from the orbital solution. Here, σ is calculated as the rootmeansquare of the fit residuals. The datapoint is discarded if there are at least three measurements of the same satellite orbit number and if it is more than 3σ away from these datapoints. If fewer than three datapoints of the same orbit number are present, they are all removed. The analysis is then iterated. This procedure is similar to the one applied by Pourbaix & Jorissen (2000) and in some cases also agrees with the outlier rejection, which is applied by the new Hipparcos reduction itself. Outliers are removed for the objects GJ 595, HD 3277, HD 43848, HD 74842, HD 154697, HD 164427A, HD 167665, HD 191760, and HIP 103019. This never requires more than one iteration.
3.7. Statistical significance of the orbit
Orbital solutions can be found for any set of astrometric data and the crucial step in the analysis is to determine the credibility of the derived orbits. We accomplish this by applying two statistical tests: the Ftest for its simplicity and the permutation test because of its lack of underlying assumptions.
The Ftest is extensively used for the statistical analysis of astrometric orbit signatures in Hipparcos data (Pourbaix & Jorissen 2000; Pourbaix & Arenou 2001; Pourbaix 2001; Reffert & Quirrenbach 2006), although it relies on the assumption that the measurement errors are Gaussian. That this assumption is not necessarily fulfilled for Hipparcos data is mentioned by Zucker & Mazeh (2001), who do not apply the Ftest to their analysis. The Ftest yields the probability that the null assumption that no orbital motion present is true, by comparing the χ^{2}values of two models with differing number of parameters. In our case, we consider the value of the 7parameter solution by nonlinear minimisation for each of the 100 draws from the spectroscopic elements. The Ftest is applied to the median value with respect to the value, obtained from the standard 5parameter solution, and results in the null probability, which is tabulated in Table 8 and gives the probability that the simpler model is valid. We use the Ftest as an additional indicator of the statistical confidence of our solution, but the conclusive argument is derived from the permutation test.
The permutation test has the important advantage of being part of the distributionfree tests (Zucker & Mazeh 2001) and therefore does not rely on the assumption of Gaussian error distribution. For our purpose, it can be summarised by the idea that a periodic signal present in a given data set will be destroyed by a random permutation of the individual datapoints in all but a few cases. In the present analysis, the periodicity is defined by the radial velocity solution. To derive a significance, the semimajor axis of the solution orbit, derived with the original astrometric data, is compared to the semimajor axes of pseudoorbits, which are derived from randomly permuted astrometric data. If a signal is present in the astrometric data, it will be destroyed by the random permutations and the solution orbit will be easily distinguishable in the pseudoorbits because of its large amplitude.
Because of the prereduction present in the original Hipparcos IAD (ESA 1997), its permutation had to be performed with great care (Zucker & Mazeh 2001). Fortunately, this became considerably easier with the availability of the new Hipparcos reduction IAD. The permutation is applied only to the list of abscissa residuals δΛ, while taking the possible contribution of the 7 and 9parameter solution into account. The remaining data, i.e. epochs, parallax factors, and scan angle orientations are left unchanged. To obtain the pseudoorbits, 1000 permutations are performed and the analysis is run on the synthesised data, using the nominal spectroscopic elements. The statistical significance of the solution orbit is equal to the percentage of pseudoorbits, which have a smaller semimajor axis than the solution orbit (Zucker & Mazeh 2001). Examples of high and lowsignificance orbits are shown in Fig. 13. In contrast to Halbwachs et al. (2000), the analysis of the pseudoorbits is performed with spectroscopic parameters identical to the nominal ones, hence the method of comparing the distribution of a/σ_{a} to the expected RayleighRice law is not applicable.
3.8. Validation of the method
The combined analysis of radialvelocity and astrometric data can lead to erroneous conclusions, if the results are not carefully checked (Han et al. 2001; Pourbaix & Arenou 2001). In particular, the relationship between the inclination and the semimajor axis has to be investigated during the interpretation, since over and underestimation of these parameters during the fitting process affects the derived companion masses. To permit us to have high confidence in our results, we validate our method by both simulation and comparison with published work.
3.8.1. Simulation
Simulations are employed to investigate the effect of Hipparcos measurement precision on the derived orbital parameters and the orbit significance. Of particular interest is the behaviour of the occurring biases as a function of S/N given by the ratio of orbital semimajor axis to measurement precision.
Instead of modelling a complete dataset, we chose the case of HD 17289, whose orbit is discussed in Sect. 4 and for which the Hipparcos observations cover twice the orbital period. To inject the astrometric signature of a companion, we use the satellite configuration for this star, i.e. the scan angles, the parallax factors, and the measurement precision. Noise is generated by adding a random bias term to each measurement, which is drawn from a Gaussian distribution with a standard deviation given by the satellite measurement precision. Configurations with different S/N are generated by altering the radialvelocity semiamplitude K_{1} of the simulated orbit, while keeping all other parameters fixed. The emulated Hipparcos errors in this simulation are consequently Gaussian, which is a simplifying assumption and may not be fulfilled in reality. However, the simulation is valid for investigating the biases of the method.
The results are illustrated in Figs. 10 and 11, where the accuracy of the determined values for inclination and semimajor axis are displayed as a function of statistical significance, derived from the permutation test. The accuracy is expressed in terms of the difference between determined and true value, i.e. for the inclination it is computed as i_{derived} − i_{true}. In general, the accuracy of the derived a and i increases with S/N and significance. For significances above 2σ, the accuracy is comparable to the precision. However, the semimajor axis retains a small bias even at high S/N, whose magnitude and sign is influenced by the satellite configuration concerning orbit orientation and period. The determination of the orbit inclination greatly improves with significance. To investigate whether our analysis is biased even when exceeding the 3σ significance, we repeated the simulation for 3 different satellite configuration, i.e. three different stars. We find no sign of a systematic bias, favouring large or small inclinations and semimajor axes. Above a 3σ significance, derived and true values always agree within the error bars. This simulation shows that a large bias can be introduced in the solution, if an orbit significance of 1σ is used to qualify a detection. In contrast, the biases are small if higher significance is requested.
We find that an orbit is detected at greater than 3σ significance, when its semimajor axis amounts to ~70% of the measurement uncertainty. This result may depend on the assumption of Gaussian errors and the specific satellite configuration, but it confirms the warning of Pourbaix (2001) that derived orbits with sizes comparable to the instrument precision have to be evaluated very carefully.
Fig. 10 Simulated difference between the derived semimajor axis a_{derived} and its true value a_{true}. Top: a_{derived} as function of a_{true}. Filled circles, open circles, and diamonds designate orbits of 1–2σ, 2–3σ, and >3σ significance, respectively. The solid line has unity slope and the dashed line shows the singlemeasurement precision of σ_{Λ} = 7.2 mas. Bottom: the relative deviation calculated as (a_{derived} − a_{true})/a_{true} as a function of significance. The accuracy and relative precision of the solution improve with significance. At the detection limit of 3σ, the deviation from the true value is smaller than 0.5 mas. 
Fig. 11 Simulated difference between the derived inclination i and its true value. The accuracy and the precision of the solution improve with significance. At the detection limit of 3σ, the deviation from the true value is smaller than 0.02°. 
Fig. 12 Joint confidence contours on the iΩgrid for the lowsignificance orbit of HD 167665 (left) and the highsignificance orbit of HD 53680 (right). The contour lines correspond to confidences at the 1σ (solid), 2σ (dashed), 3σ (dotted), and 4σ (dashdotted) levels. The crosses indicate the position of the best nonlinear adjustment solution for each of the 100 draws and the star corresponds to the adopted orbit. 
Fig. 13 Histogram showing the semimajor axes of the 1000 pseudoorbits for the lowsignificance orbit of HD 167665 (left) and the highsignificance orbit of HD 53680 (right). The vertical solid line indicates the semimajor axis of the nonpermuted, bestfit solution and the dotted line indicates the median Hipparcos singlemeasurement precision σ_{Λ} for the respective star. 
Results for the comparison sample.
3.8.2. Comparison sample
For independent validation of our method and to investigate the impact of the new reduction, we screened the literature for objects for which a similar analysis was performed. We chose the three stars HR 6046, GJ 1069, and HD 190228 with companions found by Torres (2007), Halbwachs et al. (2000), and Zucker & Mazeh (2001)^{3}, respectively, that exhibit diverse orbit inclinations and companion masses. All these works use the original IAD, whereas we use the new reduction IAD.
For all three stars, good agreement is found between our results and the published works, both in terms of value and 1σ uncertainty in semimajor axis and companion mass, as shown in Table 6. The orbit significances derived from pseudoorbits given by Zucker & Mazeh (2001) for two objects match our results. For HD 190228, we find the same solution as Zucker & Mazeh (2001) and a 2σ significance (see Sect. 4.2 for further discussion of this object). We confirm the detection of the GJ 1069 orbit, caused by a low mass star (Zucker & Mazeh 2001; Halbwachs et al. 2000), and the binary orbit of HR 6046 (Torres 2007) at greater than 3σ significance. In particular, the case of GJ 1069, which is now analysed by three independent teams finding essentially identical results, attests that our analysis is robust.Both the simulation results and the analysis of a sample of comparison stars confirm that our method for the combined analysis of radialvelocity orbits and Hipparcos astrometry is correct. The consideration of the orbit significance ensures that the derived result is biasfree to a reasonable extent.
3.9. Companion mass limits
Astrometric solution elements for targets with orbits that have significances above 1σ.
The combination of astrometric and spectroscopic data makes it possible to find a complete orbital solution. However, the result can be misleading or wrong, if the significance of the solution is not verified. For significant orbits, the derived parameters, e.g., the companion mass, and their errors are valid. For lowsignificance orbits, the derived parameters and their errors are very probably false. Therefore, the distribution function for the companion mass of a lowsignificance orbit cannot be used to derive upper and lower mass limits.
The question is then turned to the definition of a significant orbit. In the literature, various significance estimators (e.g. Ftest or permutation test) and detection limits are employed. We use a 3σ detection limit based on the permutation test and deem orbits at less than 2σ unreliable. For a 2σ orbit, we adopt a reduced validity and expect that the mass limits do hold.
In principle, it is possible to derive the upper mass limit for any companion by employing simulations similar to the one described in Sect. 3.8.1. This simulation is extremely elaborate and its application to the list of host stars with undetected astrometric motion is beyond the scope of this paper. However, we have demonstrated that an orbit is detected at greater than 3σ significance when its semimajor axis amounts to ~70% of the measurement precision, if the Hipparcos observations cover one complete orbit. Including an additional confidence margin, we thus claim that the signature of a fully covered orbit would have been detected, if its semimajor axis equalled the singlemeasurement precision. Hence, an upper limit to the companion mass can be set by enforcing the equation a = σ_{Λ} for orbits with significances smaller than 2σ and N_{orb} > 1.
4. Results
For all objects in our sample, we are able to obtain an orbital solution from the combination of the radialvelocity orbit and the Hipparcos astrometric data, but not all of these orbits are credible. To distinguish a real orbit from meaningless output of the fitting procedure, we use the permutation test to evaluate the orbit significance. This results in the detection of 6 orbits having high significance above 3σ (Tables 7 and 8). The orbits of 5 objects have moderate significance between 2σ and 3σ. Three orbits have low significance between 1σ and 2σ and the orbits of 17 objects have verylow significance (<1σ, Table 10)^{4}.
The Tables 3, 9, and 10 list the spectroscopic and astrometric elements of the adopted orbital solution together with the orbit significance and additional information: N_{orb} is the number of orbits covered by the Hipparcos observations, σ_{Λ} is the median Hipparcos singlemeasurement precision, a sin i is the minimum astrometric semimajor axis, M_{2} (3σ) are upper and lower companionmass limits at 3σ, a_{rel} is the semimajor axis of the relative orbit, is the reduced chisquare value of the adopted 7parameter solution, and the null probability gives the result of the Ftest. The quoted errors correspond to Monte Carlo based 1σ confidence intervals.
For orbits with significances below 2σ and where the orbital period is fully covered by Hipparcos observations, we derived the upper limit to the companion mass M_{2,up−lim} by enforcing the equation a = σ_{Λ} (cf. Sect. 3.9).
Solution parameters and companion masses for stars with orbits that have significances above 1σ.
Spectroscopic elements of literature targets with orbits of greater than 1σ significance.
Parameters for targets with verylow significance orbits (<1σ).
4.1. Orbits with high significance exceeding 3σ
Clear orbital signatures are detected for six target stars. The orbit significances exceed 99.7%, which we adopt as a criterion for undoubtful detection. The orbits show a large variety in size (a = 2.3–35 mas) and inclination (i = 12–173°) and reveal companion masses ranging from 90 M_{J} to 0.52 M_{⊙}. The obtained relative precision of the companion mass is typically of the order of 10% and does not include any contribution by the primarymass uncertainty. The visualisation of Hipparcos astrometric data is not very intuitive, because positions are measured in one dimension only, i.e. along the scanangle orientation ψ. We choose a representation similar to Torres (2007) and the resulting stellar astrometric orbits are shown in Fig. 20. For a clearer representation, we compute normal points for each satellite orbit number. Their values and errors are given by the mean and standard deviation of the Hipparcos abscissae obtained during one given satellite orbit, respectively. Normal points are used only for the visualisation and not during the orbit adjustment. We prove that orbit detections are also possible in apparently difficult conditions, such as incomplete orbit coverage by the satellite (N_{orb} = 0.5 for HD 43848) and an instrument precision comparable to the the orbit size (HD 164427A). The final 7parameter fit for the detected orbits has an acceptable reduced chisquare value (–1.2) in all but one case. For HD 17289, this value is very small at 0.3 and the Hipparcos precision is unusually large at 7.2 mas, which may indicate that the astrometric errors for this star are overestimated. As for the radialvelocity measurements, the light emitted by the stellar companion of this star may have disturbed the Hipparcos observations.
Three out of these six stars have nonstandard Hipparcos solutions. Thus, first indications for the astrometric disturbance induced by the companion were already detected during the standard Hipparcos data analysis. An additional outcome of the companion solution are refined positions, distances, and proper motions of the objects (Table 7). For five stars, the parallax estimate becomes more precise, whereas its value is compatible within the error bars. Only the parallax of HIP 103019 is 1.9 mas larger than the value given by the new reduction. Changes in position and proper motion are of the order of 1 mas and 1 mas yr^{1}, respectively, with extrema of 41 mas and 22 mas yr^{1}, respectively. Our results for the individual objects are described in the following:

HD 17289: The orbit inclination is173.2 ± 0.5° and the companion mass is 0.52 ± 0.05 M_{⊙}, which makes it the most massive companion detected in this study. HD17289 is thus identified as a binary star with mass ratio q = 2. Goldin & Makarov (2007) used solely Hipparcos astrometry to find an astrometric orbit with semimajor axis , which is compatible with our more precise value of a = 10.8 ± 0.7. Because the orbit is completely characterised, we can calculate the projected angular separation between the two components. Along the orbit, their separation on the sky ranges between 15–49 mas. The massluminosity relations of Delfosse et al. (2000) estimate an apparent Vmagnitude of the secondary of V_{2} = 13.18. The primary magnitude obtained from Hipparcos is V = 7.56 (Table 1) and the expected intensity contrast between the secondary and primary is ~1:180 in the visible. The combination of close separation and moderate intensity contrast causes a small and periodic bias in the radialvelocity measurement (see Sect. 2.3).

HD 30501: The orbit inclination is deg and the companion mass is 90 ± 12 M_{J}, which is the lightest companion detected in this study. This companion is just at the threshold between a very lowmass star and a brown dwarf.

HD 43848: The orbit inclination is deg and the companion is a very low mass star with mass 0.10 ± 0.01 M_{⊙}. This value is slightly lower and more precise than the one derived by Sozzetti & Desidera (2010), who obtain based on the radial velocities of Minniti et al. (2009). This solution is found with the lowereccentricity orbit (cf. Sect. 2.2). For completeness, we include the formal astrometric solution for the less probable and highereccentricity radialvelocity orbit in Tables 7 and 8. The resulting companion mass obtained in this case would be higher at M_{2} = 0.11 ± 0.02 M_{⊙}. Hipparcos astrometry cannot be used to determine the true eccentricity, because the fit quality in terms of χ^{2} is comparable in both cases.

HD 53680: The orbit inclination is deg and the companion is a low mass star with mass 0.22 ± 0.02 M_{⊙}.

HD 164427A: The orbit inclination is deg and the companion is a low mass star with mass 0.26 ± 0.06 M_{⊙}. The stellar companion has also been detected by Zucker & Mazeh (2001), who derived a slightly less inclined orbit and therefore obtained the higher mass of 0.35 ± 0.09 M_{⊙}.

HIP 103019: The orbit inclination is deg and the companion is a very low mass star with mass 0.18 ± 0.03 M_{⊙}.
4.2. Orbits with moderate significances between 2σ and 3σ
The orbits of 5 stars are found to have significances of 95.4–99.7%. Although we do not consider that these solutions accurately characterise the stellar orbits, the moderate level of significance indicates that the orbital motion is present in the astrometric data. The formal orbits are shown in Fig. 21 and allow us to obtain a visual impression of the fit quality. Four out of five orbits have sizes comparable to or smaller than the instrument precision and the astrometric signal is therefore hidden. The permutation test allows us to assess the credibility of the derived orbit especially in this difficult regime. Since there is signal in the astrometric data, the orbit inclination is constrained and we expect the companion mass to be confined within the 3σ limits given in Table 8. The individual objects are discussed below:

HD 3277: The companion is of stellar nature andits mass is confined to within 0.14–0.57 M_{⊙}. The formal companion mass at moderate significance is 0.33 ± 0.07 M_{⊙}. The upper masslimit that we independently derived from a = σ_{Λ} is M_{2,up−lim} = 0.51 M_{⊙}, which is compatible with our solution.

HD 131664: The companion is a very lowmass star or a brown dwarf with a mass of 44–153 M_{J}. The available data is not sufficient to distinguish between these two cases. The formal companion mass at moderate significance is 90 ± 16 M_{J}. Although the orbit significance agrees with the result of Sozzetti & Desidera (2010), who obtain a companion mass of with 95% confidence, the mass range derived from our analysis is significantly higher.

HD 154697: The companion is a lowmass star with a mass confined within 0.08–0.22 M_{⊙}. The formal companion mass at moderate significance is 0.15 ± 0.02 M_{⊙}. The longitude of the ascending node Ω of this orbit is hardly constrained and is correlated with the position offsets Δα^{⋆} and Δδ. The Ω − Δα^{⋆}correlation is shown in Fig. 14. Because the adopted solution is derived from the mean of the complete distribution, it does not represent a valid solution with a valid orbit, as is also visible in Fig. 21, where the fit residuals are large. However, this does not affect the validity of the derived companion mass, because that is independent of Ω.

HD 174457: The companion is a very lowmass star or a brown dwarf. Its mass is confined to within 0.06–0.18 M_{⊙}. The formal companion mass at moderate significance is 0.10 ± 0.02 M_{⊙}. The upper masslimit that we independently derived from a = σ_{Λ} is M_{2,up−lim} = 0.12 M_{⊙} and is compatible with our solution.

HD 190228: The companion is probably a brown dwarf, although we cannot exclude a very lowmass star. Assuming a primary mass of 0.83 M_{⊙} (Perrier et al. 2003), we find the formal companion mass of 49 ± 15 M_{J} with lower and upper limits of 10–99 M_{J} (cf. Table 6). The upper masslimit derived from a = σ_{Λ} independently reproduces this result and yields M_{2,up−lim} = 99 M_{J}.
Fig. 14 Correlation between the adjusted position and longitude Ω present in the distribution of the combined 100 000 Monte Carlo realisations for HD 154697. The adopted solution is marked with the diamond. 
4.3. Orbits with low significances between 1σ and 2σ
The orbits of three stars are found at low significance of 68.3–95.4%, i.e. the astrometric data contains very little or no orbital signal. As we have shown by simulation (cf. Sect. 3.8.1), orbital solutions at this significance level are prone to large biases. Therefore, the solution parameters shall not be considered valid in these cases. The formal solution of the fitting procedure is given for completeness in Tables 7 and 8. We discuss the individual objects below:

HD 74014: The poor Hipparcos coverageof 10% of this longperiod orbit impedes the detection ofthe astrometric orbit. With a minimum semimajor axis ofa sin i = 9.8 mas, the orbit would have been detected on a longer timespan, given the instrument precision of 2.5 mas. Because the orbit is not fully covered, we cannot set an upper companionmass limit M_{2,up−lim}.

HD 168443: The formal solution for the companion mass of 34 ± 9 M_{J} is in agreement with the result of Reffert & Quirrenbach (2006), who derive 34 ± 12 M_{J} using the radial velocities of Marcy et al. (2001) and the original Hipparcos data. From our analysis, using the CORALIE radial velocities and the new Hipparcos reduction IAD, we conclude that the mass obtained for the outer companion HD 168443c is of low confidence. Because the orbit is not fully covered, we cannot set an upper companionmass limit M_{2,up−lim}.

HD 191760: The minimum orbit size of this star is very small (a sin i = 0.46 mas) and Hipparcos astrometry does not reveal the signature of the companion discovered by Jenkins et al. (2009). However, we can set an upper masslimit for the companion of M_{2,up−lim} = 0.28 M_{⊙}.
4.4. Orbits with verylow significances and the detection of one browndwarf companion
The analysis of 17 objects yields orbits with significances below 1σ. Because the astrometric data does not contain orbital information, the derived result is physically meaningless and only a mathematical solution. Therefore we do not list the solution parameters. Table 10 shows the targets and their basic parameters. Two characteristics can be observed that impede the astrometric orbit detection. In the one case, applicable to both HD 65430 and HD 167665, the orbital coverage is poor and the orbit is not detected, although the expected signal is large compared to the instrument precision. In the other case, for instance for GJ 595, HD 13189, and HD 140913, the minimum orbital signature given by a sin i is very small compared to the instrument precision.
The companions of six stars have maximum masses M_{2,up−lim} above 0.63 M_{⊙}, which does not represent a considerable constraint of the system. Nine companions have upper masslimits of M_{2,up−lim} = 0.11–0.63 M_{⊙}, thus are in the Mdwarf mass range. Finally, the mass of the companion of HD 137510 is confined to be between M_{2} sin i = 22.7 M_{J} and M_{2,up−lim} = 64.4 M_{J}. Therefore, we find that HD 137510b has to be a brown dwarf. This detection was possible because of the confidence we gained in the capability of Hipparcos to detect astrometric orbital motion. Using a simpler argument, Endl et al. (2004) derived an upper masslimit of 94 M_{J}, which was not stringent enough to prove the substellar nature of the companion.
This sample also contains HD 38529, for which Benedict et al. (2010) found the orbit inclination of the outer companion HD 38529c and derived its mass of M_{2} = 17.6 M_{J} using HST FGS astrometry. Our analysis made use of the radial velocity orbit given by Benedict et al. (2010), but failed to detect an orbit reliably, although the Hipparcos precision of 1.6 mas is close to the HST orbit size (a = 1.05 ± 0.06 mas). However, the Hipparcos orbit coverage is poor at 0.4, which may explain the nondetection.
5. Discussion
The compilation of close companions to Sunlike stars with minimum masses of 10–80 M_{J} of Sozzetti & Desidera (2010) lists 39 objects. With the discovery of 9 new companions, we have increased the known number of these objects by about 20%. Ten companions in the Sozzetti & Desidera (2010) list are probably stars. We confirm the stellar nature in two cases (HD 43848B and HD 164427B) and detect 6 new stellar companions (HD 3277B, HD 17289B, HD 30501B, HD 53680B, HD 154697B, and HIP 103019B) at the 3σ significance level. All of these companions are M dwarfs with masses between 90 M_{J} and 0.52 M_{⊙}. For 3 companions, we derive mass limits that enclose the boundary between lowmass stars and brown dwarfs: HD 131664 (M_{2} = 44–153 M_{J}), HD 174457 (M_{2} = 65–190 M_{J}), and HD 190228 (M_{2} = 10–99 M_{J}). HD 190228 was initially a planethost star discovered by Sivan et al. (2004). We analysed it to validate our method and consequently recognised its importance as a promising brown dwarf candidate.
Furthermore, we set upper mass limits below 0.63 M_{⊙} on the companions of nine stars. These companions are therefore Mdwarfs or brown dwarfs. The upper masslimits of 12 companions within our sample remain unconstrained because either the derived limit is very high (>0.63 M_{⊙}) or the orbital period is not covered by Hipparcos measurements. Finally, we found that the companion of HD 137510 is a brown dwarf with a mass range of 22.7–64.4 M_{J}.
Our results are summarised in Fig. 15, where we show the companionmass constraints obtained from Hipparcos astrometry for all targets in our analysis sample.
Fig. 15 Constraints on the companion masses of the target sample. Open squares show the minimum companion mass M_{2} sin i derived from radial velocities. Blue filled circles indicate the formal companion masses M_{2} derived from astrometry for orbits with >2σsignificance. Greyshaded areas show the companion mass interval constrained by astrometry. The vertical dashed line indicates the adopted threshold mass of 80 M_{J} between brown dwarfs and stars. Twelve systems could not be constrained by Hipparcos astrometry, either because the derived mass limit was very high (>0.63 M_{⊙}) or because the orbital period was not covered by Hipparcos measurements. 
Fig. 16 Comparison of the orbit significance indicators derived from the permutation test and from the Ftest. 
Fig. 17 Inclinations and astrometric semimajor axes of the orbital solutions presented in Tables 7 and 8. The inclination is displayed as angular deviation from a faceon orbit (corresponding to i = 0° or i = 180°). Filled symbols indicate the orbits with highsignificance above 3σ, while open circles stand for orbits with moderate significance between 2σ and 3σ. 
5.1. Combining radial velocities and Hipparcos astrometry
We have used the astrometric data of the new Hipparcos reduction to place constraints on the masses of potential browndwarf companions detected in radialvelocity surveys. Radialvelocity measurements of a stellar orbit yield five of the seven Keplerian elements but do not constrain the orbit’s inclination and longitude of the ascending node. We have demonstrated how these two remaining elements can be determined, if the orbital signature is present in the Hipparcos astrometry, and how to distinguish a significant orbit.
We used the permutation test to define the criterion for orbit detection. The implementation of this powerful test is simplified considerably by the use of the new Hipparcos reduction, because the individual measurements of the IAD are independent. A source of possible errors and complications in the original Hipparcos data is therefore avoided. We emphasize the importance of eliminating outliers of the IAD, because otherwise the final result can be falsified. On a couple of comparison targets drawn from the literature, we have verified that our analysis method is robust. For these targets, the performance of the new and the original Hipparcos IAD is comparable and no major difference in the results was observed.
In addition to the permutation test, we simultaneously performed the Ftest and can thus compare the performance of both significance indicators, keeping in mind that the Ftest is only valid if Gaussian errors are assumed. Figure 16 displays the respective values listed in Tables 8 and 10. We find a general agreement between the two indicators, i.e. the null probability decreases with increasing permutation significance. However, the permutation test is more stringent. If we had chosen to rely on the Ftest and imposed a null probability of 0.3% as detection limit, we would have detected 10 orbits. These orbits overlap with the 6 orbits that we did detect with the criterion based on the permutation test (cf. Table 8).
A problem of similar studies in the past was the appearance of biases towards small orbit inclinations, thus towards high companion masses (Han et al. 2001; Pourbaix 2001). Using simulations, we showed that our method is not affected by this bias when using an appropriate detection limit in terms of orbit significance. Figure 17 shows the distributions of orbital inclination and semimajor axis for all orbits with greater than 2σ significance. No preference for very small inclinations can be observed, as it would be expected from biased solutions (Pourbaix 2001). However, we note that no orbit is detected at inclinations above 50°. This is explicable because in this study, we have selected stars with minimum companion masses of M_{2} sin i = 13–80 M_{J}. As shown by Grether & Lineweaver (2006), the frequency of browndwarf companions is very low compared to the stellar companions and therefore small inclinations are more likely to be detected.
We showed that our analysis is ultimately limited by the astrometric precision of the satellite and therefore cannot constrain orbits of size typically smaller than a = 2 mas, which in most cases will impede the detection of a true brown dwarf companion. In addition, the Hipparcos mission duration of about 1000 days limits the periods of detectable orbits to ~2000 days.
In summary, we have developed a new method to combine a radial velocity solution with Hipparcos astrometric data that has several advantages compared to previous similar approaches: it uses the new Hipparcos reduction, which simplifies both the mathematical formulation and the performance of the permutation test in contrast to the implementation of e.g. Zucker & Mazeh (2001) for the original Hipparcos data. The new method defines rigorous detection limits in terms of orbit significance and yields companion mass limits for lowsignificance orbits, provided they are covered by astrometric measurements. It is stringent enough to avoid the output of biased results (an example is the case of HD 38529, see Reffert & Quirrenbach 2006; and Benedict et al. 2010). Because of the general format of the new reduction IAD, our method is universal and applicable to any astrometric data^{5}.
Fig. 18 Cumulative mass distribution of potential browndwarf companions in the CORALIE survey. The blue dashed line shows the distribution of all 21 candidates. The black solid line shows the distribution for the 11 remaining candidates after removal of the 10 stellar companions. The companion of HD 38529 with M_{2} = 17.6 M_{J} (Benedict et al. 2010) is included. For comparison, the red dashdotted line shows the expected cumulative distribution if all companions had masses of 80 M_{J} based on the assumption that the orbits are randomly oriented in space. 
Fig. 19 Histogram of minimum masses of potential browndwarf companions in the CORALIE survey. The lightgrey histogram shows the distribution of all 21 candidates. The darkgrey histogram shows the distribution of the 11 remaining candidates after removal of the 10 stellar companions. 
5.2. Browndwarf companions in the CORALIE survey
Fig. 20 Visualisation of the highsignificance orbits. Top panels: astrometric stellar orbits projected on the sky. North is up and east is left. The solid red line shows the orbital solution and open circles mark the individual Hipparcos measurements. Dashed lines with orientation along the scan angle ψ and length given by the O–C residual of the orbital solution connect the measurements with the predicted location from our model. The blue solid circles show the normal points for each satellite orbit number. The curl at the lower left corner indicates the orientation of orbital motion. Bottom panels: O–C residuals for the normal points of the orbital solution (filled blue circles) and of the standard 5parameter model without companion (black crosses). The error bars of the normal points correspond to the dispersion of Hipparcos measurements if there are several per satellite orbit and to the individual Hipparcos abscissa error if there is only one measurement. 
Fig. 21 Visualisation of the formal solutions for the moderatesignificance orbits. The panels are analogous to Fig. 20. The formal orbit of HD 154697 is not valid and exhibits large residuals (see the text for explanation). It is shown here to illustrate the possible problems that can occur when visualising astrometric orbits. Despite their moderate significance, the four remaining orbits visually appear well constrained. 
The CORALIE planetsearch sample contains about 1600 Sunlike stars within 50 pc. After a duration of 12 years, 21 potential brown companions have been observed in this survey. Eighteen have been discussed in this work and three additional hosts of potential browndwarf companions are listed by Sozzetti & Desidera (2010)^{6} and contained in the CORALIE sample. Assuming that most candidates contained in this survey have already been found, we can derive a typical frequency of 1.3% for potential browndwarf candidates in close orbits (<10 AU) around Sunlike stars. Seven of these candidates are identified as stellar companions by this study. In addition, Zucker & Mazeh (2001) showed that the companions of HD 18445, HD 112758, and HD 217580 are of stellar nature. Thus, at least 10 out of 21 browndwarf candidates are lowmass stars, which corresponds to a rate of 48%. Based on the uniform stellar sample surveyed by the CORALIE planet search, we thus obtain an upper limit of 0.6% for the frequency of closein browndwarf companions around Sunlike stars within 50 pc, which is close to the frequency of <0.5% initially quoted by Marcy & Butler (2000).
Figure 18 shows the distribution of M_{2} sin i in the browndwarf mass range for the CORALIE sample. While the cumulative distribution (blue dashed line) of all 21 candidates is reasonably compatible with a linear increase, i.e. a flat distribution function, the situation changes significantly after removal of the 10 stellar companions. As illustrated in Fig. 19, most stellar companions have M_{2} sin i > 45 M_{J}, and after their removal, a cumulative distribution with particular shape emerges (black solid line in Fig. 18). It exhibits a steep ascent in the range of 13–25 M_{J}, which contains one half of the objects, followed by a slower increase up to 60 M_{J}. Surprisingly, there is no companion left of minimum mass higher than 60 M_{J}.
In summary, the emerging cumulative distribution does not support a uniform distribution function of companion masses. The possibility of all companions in the 13–25 M_{J} range being stellar companions seen at low inclination, is both statistically unlikely and also inplausible because their astrometric signature would be large and probably detectable. Instead, these companions supposedly have M_{2} ≈ M_{2} sin i and may represent the highmass tail of the planetary distribution.
We note that the distribution found here, which is based on radial velocities, is compatible with the first browndwarf companions found in transiting systems, such as CoRoT3b (M_{2} = 21.7 M_{J}, Deleuil et al. 2008) and CoRoT15b (M_{2} = 63.3 M_{J}, Bouchy et al. 2010).
5.3. Metallicities of the browndwarf host stars
We showed that the companion of HD 137510 is a brown dwarf. In addition, the companion of HD 190228 can be considered to be a brown dwarf at 95% confidence. The host stars have metallicities of [Fe/H] = 0.373 dex (Butler et al. 2006) and [Fe/H] = −0.24 dex (Perrier et al. 2003), respectively. Compared to the mean metallicity ⟨ [Fe/H] ⟩ = −0.09 dex with rootmeansquare σ_{⟨ [Fe/H] ⟩} = 0.24 dex derived for the uniform sample of 451 nearby Sunlike stars by Sousa et al. (2008), the two browndwarf hosts are therefore metalrich in one case and unremarkable in the other. The question of whether the hosts of brown dwarf companions have a particular metallicity distribution, as the planet hosts do (Santos et al. 2001), cannot be answered with this small sample, thus remains open.
6. Conclusions
We have presented the radialvelocity solutions of 15 stars with potential brown dwarf companions, of which 9 are new discoveries. At this point, the CORALIE and HARPS surveys have increased the number of known potential browndwarf companions of Sunlike stars by 20%.
Eight companions in our search sample were identified as being of stellar nature with masses of above 80 M_{J}. These include one companion with a mass of 90 M_{J} and thus very close to the browndwarf mass boundary. The upper and lower mass limits of three companions were derived, but enclose the starbrown dwarf mass boundary and can therefore not be used to identify the object’s nature. Yet, the companion of HD 190228 is a brown dwarf at 95% confidence. Upper mass limits were set for the companions of nine stars and revealed a brown dwarf orbiting HD 137510. The masses of 12 companions could not be constrained.
Our findings are in accordance with the presence of the browndwarf desert, because only three brown dwarfs are identified in our search sample of 33 companions: HD 190228b (at 95% confidence), HD 137510b, and HD 38529b (Benedict et al. 2010). In contrast, eight browndwarf candidates were identified as M dwarfs. On the basis of the CORALIE planetsearch sample, we obtain an upper limit of 0.6% for the frequency of browndwarf companions around Sunlike stars and confirm their pronounced paucity. We find that the companionmass distribution function rises at the lower end of the browndwarf mass range, suggesting that we detect the highmass tail of the planetary distribution function.
We find good agreement with the previous similar studies of Halbwachs et al. (2000) and Zucker & Mazeh (2001). Our results are affected by the Hipparcos astrometric precision and mission duration, which limits the minimum detectable mass of the combined analysis. There are indications, that some of the remaining potential browndwarf companions are indeed brown dwarfs. Those could be identified in the future, when astrometry of higher precision is available. An example is the case of HD 38529c, whose mass we were unable to constrain with Hipparcos data, but that was identified as a brown dwarf by Benedict et al. (2010) using HST astrometry.
Higherprecision astrometry will be required to obtain an accurate census of lowmass companions in close orbit around Sunlike stars and derive the relative frequencies of planets, brown dwarfs, and lowmass stars. The next major step in this direction is the GAIA astrometry satellite, which is the successor of Hipparcos and will make it possible to determine the masses of many radialvelocity companions.
The web interface is http://stev.oapd.inaf.it/cgibin/param
By comparing the published radial velocities to our CORALIE measurements, we discovered that the published Julian dates corresponding to the radial velocities of Minniti et al. (2009) are misprinted and offset by –2000 days.
HD 190228 has also been analysed by Pourbaix & Arenou (2001).
As motivated in Sect. 2.3, we do not present an astrometric analysis for HD 4747 and HD 211847 and therefore the number of objects listed here adds up to 31 instead of 33.
This is true provided that the implicit approximations are valid, i.e. the astrometric precision is larger than about 1 mas, see e.g. Konacki et al. (2002).
Acknowledgments
We thank the numerous CORALIE observers and the Geneva observatory technical staff at the Swiss telescope for assuring highquality radialvelocity measurements over the last decade. We thank the Swiss National Science Foundation and Geneva University for continuous support of our planet search programmes. N.C.S. would like to thank the support by 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 a Ciência 2007 contract funded by FCT/MCTES (Portugal) and POPH/FSE (EC), and in the form of grant reference PTDC/CTEAST/098528/2008. S.Z. acknowledges financial support by the Israel Science Foundation (grant No. 119/07). This research has made use of the Smithsonian/NASA Astrophysics Data System (ADS) and of the Centre de Données astronomiques de Strasbourg (CDS). All calculations were performed with the freeware Yorick.
References
 Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bean, J. L., McArthur, B. E., Benedict, G. F., et al. 2007, AJ, 134, 749 [Google Scholar]
 Benedict, G. F., McArthur, B. E., Bean, J. L., et al. 2010, AJ, 139, 1844 [NASA ADS] [CrossRef] [Google Scholar]
 Bouchy, F., Hébrard, G., Udry, S., et al. 2009, A&A, 505, 853 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bouchy, F., Deleuil, M., Guillot, T., et al. 2010, A&A, 525, A68 [Google Scholar]
 Burgasser, A. J., Reid, I. N., Siegler, N., et al. 2007, Protostars and Planets V, 427 [Google Scholar]
 Burrows, A., Marley, M., Hubbard, W. B., et al. 1997, ApJ, 491, 856 [NASA ADS] [CrossRef] [Google Scholar]
 Butler, R. P., Wright, J. T., Marcy, G. W., et al. 2006, ApJ, 646, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Casertano, S., Lattanzi, M. G., Sozzetti, A., et al. 2008, A&A, 482, 699 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chabrier, G., & Baraffe, I. 2000, ARA&A, 38, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Correia, A. C. M., Udry, S., Mayor, M., et al. 2005, A&A, 440, 751 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 da Silva, L., Girardi, L., Pasquini, L., et al. 2006, A&A, 458, 609 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deleuil, M., Deeg, H. J., Alonso, R., et al. 2008, A&A, 491, 889 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Delfosse, X., Forveille, T., Ségransan, D., et al. 2000, A&A, 364, 217 [NASA ADS] [Google Scholar]
 Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 [NASA ADS] [Google Scholar]
 Endl, M., Hatzes, A. P., Cochran, W. D., et al. 2004, ApJ, 611, 1121 [NASA ADS] [CrossRef] [Google Scholar]
 ESA 1997, VizieR Online Data Catalog, 1239, 0 [Google Scholar]
 van Leeuwen, F. 2007, Astrophysics and Space Science Library, Vol. 350, Hipparcos, the New Reduction of the Raw Data [Google Scholar]
 Fischer, D. A., Marcy, G. W., Butler, R. P., et al. 2003, ApJ, 586, 1394 [NASA ADS] [CrossRef] [Google Scholar]
 Forveille, T., Bonfils, X., Delfosse, X., et al. 2009, A&A, 493, 645 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Galland, F., Lagrange, A., Udry, S., et al. 2006, A&A, 452, 709 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 [Google Scholar]
 Goldin, A., & Makarov, V. V. 2007, ApJS, 173, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Good, P. 1994, Permutation Tests – A Practical Guide to Resampling Methods for Testing Hypotheses (NewYork: SpringerVerlag) [Google Scholar]
 Grether, D., & Lineweaver, C. H. 2006, ApJ, 640, 1051 [NASA ADS] [CrossRef] [Google Scholar]
 Halbwachs, J. L., Arenou, F., Mayor, M., Udry, S., & Queloz, D. 2000, A&A, 355, 581 [NASA ADS] [Google Scholar]
 Halbwachs, J. L., Mayor, M., Udry, S., & Arenou, F. 2003, A&A, 397, 159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Han, I., Black, D. C., & Gatewood, G. 2001, ApJ, 548, L57 [NASA ADS] [CrossRef] [Google Scholar]
 Hatzes, A. P., Guenther, E. W., Endl, M., et al. 2005, A&A, 437, 743 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hilditch, R. W. 2001, An Introduction to Close Binary Stars (Cambridge University Press) [Google Scholar]
 Hummel, C. A., Armstrong, J. T., Quirrenbach, A., et al. 1993, AJ, 106, 2486 [NASA ADS] [CrossRef] [Google Scholar]
 Jenkins, J. S., Jones, H. R. A., Goździewski, K., et al. 2009, MNRAS, 398, 911 [NASA ADS] [CrossRef] [Google Scholar]
 Konacki, M., Maciejewski, A. J., & Wolszczan, A. 2002, ApJ, 567, 566 [NASA ADS] [CrossRef] [Google Scholar]
 Lane, B. F., Muterspaugh, M. W., Fekel, F. C., et al. 2007, ApJ, 669, 1209 [NASA ADS] [CrossRef] [Google Scholar]
 Launhardt, R., Queloz, D., Henning, T., et al. 2008, in Proc. SPIE, 7013 [Google Scholar]
 Lindegren, L. 2010, in IAU Symp., ed. S. A. Klioner, P. K. Seidelmann, & M. H. Soffel, 261, 296 [Google Scholar]
 Liu, Y., Sato, B., Zhao, G., et al. 2008, ApJ, 672, 553 [NASA ADS] [CrossRef] [Google Scholar]
 Luhman, K. L., Joergens, V., Lada, C., et al. 2007, Protostars and Planets V, 443 [Google Scholar]
 Luhman, K. L., Mamajek, E. E., Allen, P. R., & Cruz, K. L. 2009, ApJ, 703, 399 [NASA ADS] [CrossRef] [Google Scholar]
 Marcy, G. W., & Butler, R. P. 2000, PASP, 112, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Marcy, G. W., Butler, R. P., Vogt, S. S., et al. 2001, ApJ, 555, 418 [NASA ADS] [CrossRef] [Google Scholar]
 Martioli, E., McArthur, B. E., Benedict, G. F., et al. 2010, ApJ, 708, 625 [NASA ADS] [CrossRef] [Google Scholar]
 Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
 Mayor, M., Udry, S., Naef, D., et al. 2004, A&A, 415, 391 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mayor, M., Bonfils, X., Forveille, T., et al. 2009, A&A, 507, 487 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Minniti, D., Butler, R. P., LópezMorales, M., et al. 2009, ApJ, 693, 1424 [NASA ADS] [CrossRef] [Google Scholar]
 Moutou, C., Mayor, M., Lo Curto, G., et al. 2009, A&A, 496, 513 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Muterspaugh, M. W., Lane, B. F., Konacki, M., et al. 2006, ApJ, 636, 1020 [NASA ADS] [CrossRef] [Google Scholar]
 Naef, D., Mayor, M., Pepe, F., et al. 2001, A&A, 375, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Naef, D., Mayor, M., Lo Curto, G., et al. 2010, A&A, 523, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nidever, D. L., Marcy, G. W., Butler, R. P., Fischer, D. A., & Vogt, S. S. 2002, ApJS, 141, 503 [NASA ADS] [CrossRef] [Google Scholar]
 Omiya, M., Izumiura, H., Han, I., et al. 2009, PASJ, 61, 825 [NASA ADS] [Google Scholar]
 Patel, S. G., Vogt, S. S., Marcy, G. W., et al. 2007, ApJ, 665, 744 [Google Scholar]
 Pepe, F., Correia, A. C. M., Mayor, M., et al. 2007, A&A, 462, 769 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Perrier, C., Sivan, J., Naef, D., et al. 2003, A&A, 410, 1039 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Perryman, M. A. C., Lindegren, L., Arenou, F., et al. 1996, A&A, 310, L21 [NASA ADS] [Google Scholar]
 Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 [NASA ADS] [Google Scholar]
 Pourbaix, D. 2001, A&A, 369, L22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pourbaix, D., & Arenou, F. 2001, A&A, 372, 935 [Google Scholar]
 Pourbaix, D., & Jorissen, A. 2000, A&AS, 145, 161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pourbaix, D., Tokovinin, A. A., Batten, A. H., et al. 2004, A&A, 424, 727 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Queloz, D., Mayor, M., Weber, L., et al. 2000, A&A, 354, 99 [NASA ADS] [Google Scholar]
 Reffert, S., & Quirrenbach, A. 2006, A&A, 449, 699 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Israelian, G., & Mayor, M. 2001, A&A, 373, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Israelian, G., & Mayor, M. 2004, A&A, 415, 1153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Israelian, G., Mayor, M., et al. 2005, A&A, 437, 1127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Mayor, M., Bouchy, F., et al. 2007, A&A, 474, 647 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, N. C., Mayor, M., Naef, D., et al. 2000, A&A, 361, 265 [NASA ADS] [Google Scholar]
 Santos, N. C., Mayor, M., Naef, D., et al. 2002, A&A, 392, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ségransan, D., Udry, S., Mayor, M., et al. 2010, A&A, 511, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sivan, J., Mayor, M., Naef, D., et al. 2004, in Planetary Systems in the Universe, ed. A. Penny, IAU Symp., 202, 124 [Google Scholar]
 Sousa, S. G., Santos, N. C., Mayor, M., et al. 2008, A&A, 487, 373 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Sozzetti, A., & Desidera, S. 2010, A&A, 509, A260000 [Google Scholar]
 Spiegel, D. S., Burrows, A., & Milsom, J. A. 2010, ApJ, accepted [arXiv:1008.5150] [Google Scholar]
 Tamuz, O., Ségransan, D., Udry, S., et al. 2008, A&A, 480, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tinney, C. G., Butler, R. P., Marcy, G. W., et al. 2001, ApJ, 551, 507 [NASA ADS] [CrossRef] [Google Scholar]
 Tokovinin, A., Thomas, S., Sterzik, M., & Udry, S. 2006, A&A, 450, 681 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Torres, G. 2007, AJ, 134, 1916 [NASA ADS] [CrossRef] [Google Scholar]
 Udry, S., Mayor, M., Naef, D., et al. 2000, A&A, 356, 590 [NASA ADS] [Google Scholar]
 Udry, S., Mayor, M., Naef, D., et al. 2002, A&A, 390, 267 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Udry, S., Mayor, M., & Queloz, D. 2004, in Planetary Systems in the Universe, ed. A. Penny, IAU Symp., 202, 12 [Google Scholar]
 Udry, S., & Santos, N. C. 2007, ARA&A, 45, 397 [NASA ADS] [CrossRef] [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Whitworth, A., Bate, M. R., Nordlund, Å., Reipurth, B., & Zinnecker, H. 2007, Protostars and Planets V, 459 [Google Scholar]
 Wittenmyer, R. A., Endl, M., Cochran, W. D., et al. 2009, AJ, 137, 3529 [Google Scholar]
 Wright, J. T., Upadhyay, S., Marcy, G. W., et al. 2009, ApJ, 693, 1084 [NASA ADS] [CrossRef] [Google Scholar]
 Zucker, S., & Mazeh, T. 2000, ApJ, 531, L67 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Zucker, S., & Mazeh, T. 2001, ApJ, 562, 549 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Astrometric solution elements for targets with orbits that have significances above 1σ.
Solution parameters and companion masses for stars with orbits that have significances above 1σ.
Spectroscopic elements of literature targets with orbits of greater than 1σ significance.
All Figures
Fig. 1 Orbital characteristics of the 15 browndwarf candidates observed with CORALIE and HARPS. Filled symbols indicate the companions with stellar mass identified by this study. The error bars are often smaller than the symbol size, though the orbits of HD 4747 and HD 211847 are not covered and therefore less constrained, see the text for details. Top: minimum companion mass and separation. Bottom: eccentricity and orbital period. 

In the text 
Fig. 2 Radial velocities (top) and residuals (bottom) to the best fit, represented by the solid line, of HD 4747 measured with HIRES (orange diamonds), CORALIE C98 (red circles), and CORALIE C07 (blue squares). 

In the text 
Fig. 3 Radial velocities (top) and residuals (bottom) to the best fit, represented by the solid line, of HD 211847 measured with CORALIE C98 (red circles) and CORALIE C07 (blue squares). 

In the text 
Fig. 4 Phasefolded radial velocities of 8 stars with potential browndwarf companions. Red circles and blue squares indicate measurements with CORALIE C98 and C07, respectively. For HD 43848, the orange diamonds show the MIKE measurements. The solid lines correspond to the bestfit solutions. The error bars are smaller than the symbol size. 

In the text 
Fig. 5 Phasefolded radial velocities of 5 stars with potential browndwarf companions. Error bars are smaller than the symbol size. 

In the text 
Fig. 6 Bisector velocity span (top) and FWHM of the crosscorrelation function (bottom) as a function of orbital phase for HD 17289. The corresponding radial velocities are shown in Fig. 4. Both observables clearly show a phasedependent modulation, which is likely to be caused by the companion. 

In the text 
Fig. 7 Eccentricities (left) and velocity semiamplitudes (right) of 10 000 Monte Carlo simulations performed to find the radialvelocity orbit of HD 43848. The bimodal structure is clearly seen in both histograms. The final solution is found by considering only the lower eccentricity realisations, which is done by separating the distributions at m s^{1} and correspondingly at ê ~ 0.78. The probability of obtaining a high or loweccentricity solution is calculated by counting the number of realisations above and below this threshold value and equals 77% and 23%, respectively. 

In the text 
Fig. 8 CaII H emission region for HD 89707 showing slight reemission at the bottom of this absorption line. 

In the text 
Fig. 9 CaII H emission region for HD 211847 showing reemission at the bottom of the absorption line. 

In the text 
Fig. 10 Simulated difference between the derived semimajor axis a_{derived} and its true value a_{true}. Top: a_{derived} as function of a_{true}. Filled circles, open circles, and diamonds designate orbits of 1–2σ, 2–3σ, and >3σ significance, respectively. The solid line has unity slope and the dashed line shows the singlemeasurement precision of σ_{Λ} = 7.2 mas. Bottom: the relative deviation calculated as (a_{derived} − a_{true})/a_{true} as a function of significance. The accuracy and relative precision of the solution improve with significance. At the detection limit of 3σ, the deviation from the true value is smaller than 0.5 mas. 

In the text 
Fig. 11 Simulated difference between the derived inclination i and its true value. The accuracy and the precision of the solution improve with significance. At the detection limit of 3σ, the deviation from the true value is smaller than 0.02°. 

In the text 
Fig. 12 Joint confidence contours on the iΩgrid for the lowsignificance orbit of HD 167665 (left) and the highsignificance orbit of HD 53680 (right). The contour lines correspond to confidences at the 1σ (solid), 2σ (dashed), 3σ (dotted), and 4σ (dashdotted) levels. The crosses indicate the position of the best nonlinear adjustment solution for each of the 100 draws and the star corresponds to the adopted orbit. 

In the text 
Fig. 13 Histogram showing the semimajor axes of the 1000 pseudoorbits for the lowsignificance orbit of HD 167665 (left) and the highsignificance orbit of HD 53680 (right). The vertical solid line indicates the semimajor axis of the nonpermuted, bestfit solution and the dotted line indicates the median Hipparcos singlemeasurement precision σ_{Λ} for the respective star. 

In the text 
Fig. 14 Correlation between the adjusted position and longitude Ω present in the distribution of the combined 100 000 Monte Carlo realisations for HD 154697. The adopted solution is marked with the diamond. 

In the text 
Fig. 15 Constraints on the companion masses of the target sample. Open squares show the minimum companion mass M_{2} sin i derived from radial velocities. Blue filled circles indicate the formal companion masses M_{2} derived from astrometry for orbits with >2σsignificance. Greyshaded areas show the companion mass interval constrained by astrometry. The vertical dashed line indicates the adopted threshold mass of 80 M_{J} between brown dwarfs and stars. Twelve systems could not be constrained by Hipparcos astrometry, either because the derived mass limit was very high (>0.63 M_{⊙}) or because the orbital period was not covered by Hipparcos measurements. 

In the text 
Fig. 16 Comparison of the orbit significance indicators derived from the permutation test and from the Ftest. 

In the text 
Fig. 17 Inclinations and astrometric semimajor axes of the orbital solutions presented in Tables 7 and 8. The inclination is displayed as angular deviation from a faceon orbit (corresponding to i = 0° or i = 180°). Filled symbols indicate the orbits with highsignificance above 3σ, while open circles stand for orbits with moderate significance between 2σ and 3σ. 

In the text 
Fig. 18 Cumulative mass distribution of potential browndwarf companions in the CORALIE survey. The blue dashed line shows the distribution of all 21 candidates. The black solid line shows the distribution for the 11 remaining candidates after removal of the 10 stellar companions. The companion of HD 38529 with M_{2} = 17.6 M_{J} (Benedict et al. 2010) is included. For comparison, the red dashdotted line shows the expected cumulative distribution if all companions had masses of 80 M_{J} based on the assumption that the orbits are randomly oriented in space. 

In the text 
Fig. 19 Histogram of minimum masses of potential browndwarf companions in the CORALIE survey. The lightgrey histogram shows the distribution of all 21 candidates. The darkgrey histogram shows the distribution of the 11 remaining candidates after removal of the 10 stellar companions. 

In the text 
Fig. 20 Visualisation of the highsignificance orbits. Top panels: astrometric stellar orbits projected on the sky. North is up and east is left. The solid red line shows the orbital solution and open circles mark the individual Hipparcos measurements. Dashed lines with orientation along the scan angle ψ and length given by the O–C residual of the orbital solution connect the measurements with the predicted location from our model. The blue solid circles show the normal points for each satellite orbit number. The curl at the lower left corner indicates the orientation of orbital motion. Bottom panels: O–C residuals for the normal points of the orbital solution (filled blue circles) and of the standard 5parameter model without companion (black crosses). The error bars of the normal points correspond to the dispersion of Hipparcos measurements if there are several per satellite orbit and to the individual Hipparcos abscissa error if there is only one measurement. 

In the text 
Fig. 21 Visualisation of the formal solutions for the moderatesignificance orbits. The panels are analogous to Fig. 20. The formal orbit of HD 154697 is not valid and exhibits large residuals (see the text for explanation). It is shown here to illustrate the possible problems that can occur when visualising astrometric orbits. Despite their moderate significance, the four remaining orbits visually appear well constrained. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.