{'fraction_non_alphanumeric': 0.04164931278633902, 'fraction_numerical': 0.03611590408758256, 'mean_word_length': 4.216635630043451, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
HD 5388 b is a 69 M_Jup companion instead of a planet

We examined six exoplanet host stars with non-standard Hipparcos astrometric solution, which may be indicative of unrecognised orbital motion. Using Hipparcos intermediate astrometric data, we detected the astrometric orbit of HD 5388 at a significance level of 99.4 % (2.7 sigma). HD 5388 is a metal-deficient star and hosts a planet candidate with a minimum mass of 1.96 M_J discovered in 2010. We determined its orbit inclination to be i = 178.3 +0.4/-0.7 deg and the corresponding mass of its companion HD 5388 b to be M_2 = 69 +/- 20 M_J. The orbit is seen almost face-on and the companion mass lies at the upper end of the brown-dwarf mass range. A mass lower than 13 M_J was excluded at the 3-sigma level. The astrometric motions of the five other stars had been investigated by other authors revealing two planetary companions, one stellar companion, and two statistically insignificant orbits. We conclude that HD 5388 b is not a planet but most likely a brown-dwarf companion. In addition, we find that the inclinations of the stellar rotation axis and the companion's orbital axis differ significantly.


Introduction
For most extrasolar planets detected with the radial-velocity (RV) method, only their minimum masses are known because of their unknown orbit inclinations. Although the large number of planets allows us to draw statistically sound conclusions about their mass distribution (Udry & Santos 2007;Howard et al. 2010), it is highly desirable to derive this distribution without the inclination incertitude. Astrometric measurements can resolve this ambiguity. For instance, Sahlmann et al. (2011) used Hipparcos data to analyse the astrometric motion of stars with potential brown-dwarf companions and found that about half of the candidates are low-mass stars. Here, we extend this work to stars with planet candidates and non-standard Hipparcos solution, and discover the astrometric orbit of HD 5388, which hosts a recently announced planet candidate.

Target selection and astrometric analysis
In September 2010, the list of RV-planets at exoplanet.eu contained 461 entries around 389 stars. Only 286 host stars are included in the new Hipparcos reduction (F. van Leeuwen 2007). Many host stars of transiting exoplanets discovered in transit surveys (e.g., Kepler, CoRoT, HAT, WASP) and confirmed by RV are fainter than the completeness limit of Hipparcos. We selected the stars for which the new Hipparcos reduction found a nonstandard solution and that were not flagged as member of a multiple system, thus having solution types '1', '7', or '9'. Type '1' solutions are termed stochastic and adopted when the standard five-parameter solution (type '5') is not satisfactory and neither orbital nor acceleration models improve the solution in terms of χ 2 . The types '7' and '9' are given, when the model has to include proper-motion derivatives of first and second order to obtain a reasonable fit. Six stars satisfied these criteria and are listed in Table 4.   Santos et al. (2010). Red circles indicate the HARPS measurements and the solid line corresponds to the best-fit solution.
The astrometric analysis was performed as described in Sahlmann et al. (2011), where a detailed description of the method can be found. We briefly recall the main elements of the analysis. Using the orbital parameters known from RV measurements, the intermediate astrometric data of the new Hipparcos reduction was fitted with a seven-parameter model depending on the inclination i, the longitude of the ascending node Ω, the parallax , and offsets to the coordinates (∆α , ∆δ) and proper motions (∆µ α , ∆µ δ ) given in the published catalogue of F. van Leeuwen (2007). A two-dimensional search grid in i and Ω defined starting values for a standard nonlinear χ 2 -minimisation procedure identifying the global minimum. The uncertainties in the RV parameters are propagated to the astrometric solution by means of Monte Carlo simulations. The statistical significance of each astrometric orbit was determined by the distribution-free  permutation test employing 1000 pseudo orbits. Uncertainties in the solution parameters were derived by Monte Carlo simulations. This approach has proven to be reliable in detecting orbital signatures in the Hipparcos astrometric data and efficiently distinguishing a significant orbit in the present low signal-to-noise ratio regime (Sahlmann et al. 2011).

The orbit of HD 5388
HD 5388 (HIP 4311) is listed with a stochastic solution in the catalogue of the new Hipparcos reduction. Santos et al. (2010) used the HARPS spectrograph to discover a massive planet candidate around this F6V-star. Figure 1 shows the radial velocities from the discovery paper. The companion has a minimum mass of M 2 sin i = 1.96 M J and orbits its host in an eccentric orbit (e = 0.4) with an orbital period of P = 777 days. The time of periastron passage T 0 , the longitude of periastron ω, the radialvelocity semi-amplitude K 1 , and the minimum semimajor axis of the star's astrometric motion a sin i are given in Table 1. The residuals of the RV solution (3.3 ms −1 ) are larger than the average measurement error (2.8 ms −1 ). By considering the star's spectral type and projected rotational velocity and analysing the bisector inverse slope, Santos et al. (2010) conclude that the excess RV noise is insignificant. We detected the astrometric orbit with a significance of 99.4 % (2.7 σ), which means that Hipparcos measured the stellar orbit induced by a companion with the orbital characteristics given in Table 1. The derived astrometric orbit with a = 1.7 ± 0.5 milli-arcsec (mas) is small compared to the median single-measurement precision σ Λ = 3.8 mas given by the new Hipparcos reduction. However, the large number of measurements N mes = 191 covering more than one orbital revolution (N orb = 1.5) resulted in an effective signal-to-noise ratio of S/N = a · (σ Λ / √ N mes ) −1 = 6.2 allowing the detection to be possible. With a reduced chi-square value of χ 2 red = 1.02 the final  Table 2. Reffert & Quirrenbach (2011) found an upper mass limit of 124 M J for the companion of HD 5388, which is compatible with our result.
Assuming that the orbit orientations are randomly oriented in space, the probability of measuring an inclination smaller than 2 • is only 0.06 %. However, the number of planets detected in RV-surveys reached several hundreds making our detection reasonably probable. HD 5388 is a metal-deficient star with [Fe/H] = −0.27 ± 0.02. The hydrogen-burning mass limit, commonly used to distinguish between stars and brown dwarfs, increases with decreasing metallicity and reaches 87 M J (0.083 M ) for [Fe/H] < −1 (Chabrier & Baraffe 1997). Because of the moderate metal-deficiency of HD 5388, we used a limit of 80 M J and found that the probability that the companion of HD 5388 has  Fig. 3. Probability density functions of orbit inclination (i orbit , dark-grey histogram) and stellar spin axis inclination (i rot , lightgrey histogram). The orbit inclination PDF is very sharply peaked with its maximum at 0.18 and is truncated for display clarity. The PDF of ψ = i orbit − i rot is shown as dashed line. The bin size is 0.2 • . a mass below this limit is 71 %. The companion of HD 5388 is thus most likely a brown dwarf. This adds HD 5388 to the very few known Sun-like stars that have a brown-dwarf companion and a well-determined astrometric orbit.

Spin-orbit alignment
After determining the orbital inclination from astrometry, we compared it to the orientation of the stellar spin axis. The inclination i rot of the stellar spin axis is defined with respect to the line of sight (i rot = 0 • or 180 • for a pole-on view) and can be derived from the spectroscopic estimate of υ sin i rot = 4.2 km s −1 (Santos et al. 2010). The authors did not give an error bar for this measurement and we assumed a conservative uncertainty of 1 km s −1 . On the basis of the activity indicator log R H,K determined from the HARPS spectra, we used the calibration of Mamajek & Hillenbrand (2008) to derive the stellar rotation period P rot . Assuming an effective temperature of T eff = 6297 ± 32 (Santos et al. 2010), the apparent visual magnitude m V = 6.839 ± 0.001 (F. van Leeuwen 2007), and the parallax given in Table 2, we derived the star's absolute magnitude, luminosity (L = 4.5 ± 0.3 L ), and radius (R = 1.8 ± 0.1 R ), using standard formulae and Monte Carlo resampling. Bolometric corrections were computed using the Flower (1996) parameters given by Torres (2010). Using the stellar radius and rotation period, we derived the equatorial rotation velocity to be υ = 10 ± 3 km s −1 , and then the inclination of the spin axis by calculating i rot = arcsin(υ sin i rot /υ).
Because we cannot determine the star's sense of rotation, the angle i rot has a 180 • ambiguity and we first considered the value that falls into the same quadrant as i orbit , i.e. a prograde configuration. The astrometric analysis yielded the distribution of i orbit , obtained from 100 000 Monte Carlo simulations. To obtain the i rot -distribution, we performed 10 6 Monte Carlo simulations. We finally compared these two distributions by drawing 2 · 10 7 pairs of values [i orbit , i rot ] and obtained the distribution of the orbit obliquity or spin-orbit angle ψ by computing ψ = i orbit − i rot . Figure 3 shows the probability density functions (PDF) of i orbit , i rot , and ψ. The i orbit -distribution is very narrow and peaks around Table 3. Spin inclination, orbital inclination, and obliquity.
The median values and confidence intervals of i rot , i orbit , and the obliquity ψ are given in Table 3. We found that ψ is larger than 4.4 • and 19.3 • with a probability of 99.7% and 68.3 %, respectively, and its nominal value with 1 σ confidence intervals is ψ = 23 +10 • −8 in prograde configuration, which indicates that the orientations of stellar spin and orbital axis differ substantially. In retrograde configuration, the obliquity would be ψ = 153 +8 • −10 . We note that these values of ψ and ψ are lower limits, because the ascending node Ω rot of the spin axis is not constrained and we assumed that Ω = Ω rot in their derivation. As illustrated in Fig. 4, ψ increases if the two ascending nodes Ω rot and Ω do not coincide, creating an additional uncertainty in ψ of 2 · (180 • − i orbit ) = 3.4 • in the prograde configuration. For a retrograde orbit, this uncertainty is larger at 51 • . In summary, the obliquity is ψ = 23 +10 for prograde and retrograde orbits, respectively, where the additional uncertainties are denoted in square brackets.
We performed these calculations using the P rot -calibration by Noyes et al. (1984), which yielded a value of 8.8 ± 1.6 days and an obliquity of ψ = 22 +9 • −7 , in agreement with the result obtained with the calibration of Mamajek & Hillenbrand (2008). The isochrones of Marigo et al. (2008) 1 indicated that HD 5388 is at the end of the core-H burning phase just past the turn-off point. Therefore, HD 5388 may be about to leave the main sequence or even be slightly evolved. This would introduce additional uncertainties and possibly biases into our estimates of log R H,K , P rot , the stellar radius, and eventually i rot , and may therefore alter the above findings.

The orbits of five additional selected stars
In addition to HD 5388, we examined five other stars that satisfied our selection criteria in Sect. 2. They are listed in Table 4 and although their astrometric motions were previously studied in the literature, we describe our findings below: -55 Cnc has five planetary companions announced by Fischer et al. (2008). The planet 55 Cnc d has a minimum mass of M 2 sin i = 3.84 M J . Using HST astrometry, McArthur et al. (2004) derived an orbit inclination of 53±7 • for the d-planet, an astrometric perturbation size of a = 1.9 ± 0.4 mas, and a mass of M 2 = 4.9 ± 1.1 M J . Compared to the measurement precision of σ Λ = 1.4 mas, the stellar orbit is large enough to have been measured by Hipparcos, which explains the nonstandard solution type. Because the Hipparcos data only covers 20 % of the 5200 day orbital period, we did not attempt to fit an astrometric orbit. -HD 81040 hosts a massive planet with M 2 sin i = 6.86 M J (Sozzetti et al. 2006). We did not detect orbital signature because the obtained orbit has a very low significance of 68 %, in agreement with the result of Sozzetti et al. (2006). As in that work, we used the Hipparcos astrometric data to set an upper limit to the companion mass of 47.5 M J , thus excluding that the companion has a stellar nature. -HD 179949 hosts a planet candidate with M 2 sin i = 0.84 M J (Tinney et al. 2001). We did not detect an orbital signature because the obtained orbit has a very low significance of 71 %, which confirms the result of Zucker & Mazeh (2001). -HD 195019 hosts a planet candidate with M 2 sin i = 3.70 M J (Wright et al. 2007). We detected the astrometric orbit with a high significance of 98 % and a very low inclination of 0.3 ± 0.1 • , which renders a stellar companion and confirms the result of Zucker & Mazeh (2001). -γ Cep (HD 222404) hosts a planet candidate with M 2 sin i = 1.7 M J (Hatzes et al. 2003). We did not treat this system because an exhaustive analysis, including the use of Hipparcos astrometry, was performed by Torres (2007), who set an upper mass limit of 16.9 M J to the planetary companion.

Discussion
Among the potential planets around the six stars selected in Sect. 2, one turned out to be a star (HD 195019 b), two are likely planets (55 Cnc d and γ Cep b), and one is a brown-dwarf companion (HD 5388). The astrometric orbits of two stars were not detected with Hipparcos and the companions remain planetary candidates. The significant fraction of one third of non-planetary companions found in this small sample is a selection effect and does not reflect the properties of the planetary population. That a large population of non-planetary companions is observed with low orbital inclinations, hence mistaken for planets, is statistically unlikely and would have been detected with Hipparcos (Zucker & Mazeh 2001). At the individual level, however, a planet detected by RV remains a planet candidate until its mass range is confined by complementary constraints derived for instance from a transit lightcurve, dynamical interactions, by direct imaging, or from astrometric data, as in the case of HD 5388. In the future, as the data of the GAIA astrometry satellite becomes available, we will be able to measure the astrometric orbits of many planetary systems, thus obtaining an accurate census of the planetary mass distribution. A few more planet candidates may then turn out to be brown-dwarf companions.
We have detected the astrometric orbit of HD 5388, which harboured a planet candidate of minimum mass M 2 sin i = 1.96 M J . We found an almost face-on orbit with an inclination of i orbit = 178.3 +0.4 • −0.7 . Consequently, we determined the mass of the companion HD 5388 b to be M 2 = 69±20 M J . A mass lower than 13 M J was excluded at the 3 σ confidence level. Thus, HD 5388 b can no longer be considered a planet. It is instead a close browndwarf companion, which is remarkable because at most 0.6 % of Sun-like stars have close brown-dwarf companions (Sahlmann et al. 2011). In particular, it is one of the very few brown-dwarf companions for which the astrometric orbit could be measured and whose orbital motion is therefore fully characterised, such as HD 38529 (Benedict et al. 2010).
The measurement of the spin-orbit alignment in planetary and binary systems, e.g., Triaud et al. (2010) and Albrecht et al. (2011), gives insight into the formation of these systems. The most frequently applied technique uses spectroscopic measurements during the transit of eclipsing systems and only determines the projected spin-orbit angle. Because it is restricted to eclipsing systems, it favours small-separation, short-period systems (another method relies on interferometric imaging of the stellar surface, see Hinkley et al. 2011). We determined the stellar spin axis inclination on the basis of the spectroscopic estimate of υ sin i rot and the rotation period calibration using log R H,K . We found evidence of non-zero spin-orbit angles with numerical values of ψ = 23 +10 [+3.4] • −8 and ψ = 153 +8 [+51] • −10 for prograde and retrograde orbits, respectively. Thus, for both orientations of the unknown stellar spin, we found substantial non-coplanarity between the stellar and orbital axes. Spin-orbit angle measurements obtained from the determination of the astrometric orbit and the spectroscopic υ sin i offer the opportunity to access this information for non-eclipsing systems.