Issue 
A&A
Volume 574, February 2015



Article Number  A6  
Number of page(s)  32  
Section  Galactic structure, stellar clusters and populations  
DOI  https://doi.org/10.1051/00046361/201323056  
Published online  15 January 2015 
Are the orbital poles of binary stars in the solar neighbourhood anisotropically distributed?^{⋆}
^{1}
13 rue Beyle Stendhal,
38000
Grenoble
France
^{2}
Laboratoire Lagrange, UMR 7293, Univ. Nice SophiaAntipolis, CNRS,
Observatoire de la Côte d’Azur, 06300
Nice,
France
email: daniel.bonneau@oca.eu
^{3}
Institut d’Astronomie et d’Astrophysique, Université Libre de
Bruxelles, CP 226,
Boulevard du Triomphe,
1050
Brussels,
Belgium
^{4}
CEA/Saclay, DSM, 91191
GifsurYvette Cedex,
France
^{5}
Observatoire de Genève, Chemin des Maillettes, 1290
Sauverny,
Suisse
^{6}
251 vieille rue du Moulin,
1180
Brussels,
Belgium
^{7}
Observatoire Royal de Belgique, Avenue Circulaire 3, 1180
Bruxelles,
Belgique
Received: 15 November 2013
Accepted: 30 October 2014
We test whether or not the orbital poles of the systems in the solar neighbourhood are isotropically distributed on the celestial sphere. The problem is plagued by the ambiguity on the position of the ascending node. Of the 95 systems closer than 18 pc from the Sun with an orbit in the 6th Catalogue of Orbits of Visual Binaries, the pole ambiguity could be resolved for 51 systems using radial velocity collected in the literature and CORAVEL database or acquired with the HERMES/Mercator spectrograph. For several systems, we can correct the erroneous nodes in the 6th Catalogue of Orbits and obtain new combined spectroscopic/astrometric orbits for seven systems [WDS 01083+5455Aa,Ab; 01418+4237AB; 02278+0426AB (SB2); 09006+4147AB (SB2); 16413+3136AB; 17121+4540AB; 18070+3034AB]. We used of spherical statistics to test for possible anisotropy. After ordering the binary systems by increasing distance from the Sun, we computed the falsealarm probability for subsamples of increasing sizes, from N = 1 up to the full sample of 51 systems. RayleighWatson and Beran tests deliver a falsealarm probability of 0.5% for the 20 systems closer than 8.1 pc. To evaluate the robustness of this conclusion, we used a jackknife approach, for which we repeated this procedure after removing one system at a time from the full sample. The falsealarm probability was then found to vary between 1.5% and 0.1%, depending on which system is removed. The reality of the deviation from isotropy can thus not be assessed with certainty at this stage, because only so few systems are available, despite our efforts to increase the sample. However, when considering the full sample of 51 systems, the concentration of poles toward the Galactic position l = 46.0°, b = 37°, as observed in the 8.1 pc sphere, totally vanishes (the RayleighWatson falsealarm probability then rises to 18%).
Key words: binaries: visual / binaries: spectroscopic / techniques: radial velocities / techniques: high angular resolution / methods: statistical / solar neighborhood
Tables 1–3 and Appendices are available in electronic form at http://www.aanda.org† Deceased October 1, 2014.
© ESO, 2015
1. Introduction
In 1838, Mädler announced the probable existence of an anastonishing phenomenon related to the orbital plane of visual double stars: they seem to be aligned with each other. This phenomenon has received quite some attention since then, but the most earliest studies that was based on data secured with old observational techniques and equipments, did not bring significant progress. In the years 1950–1967, Dommanget (with the assistance of Nys) reinvestigated the question after collecting all available data needed to perform the first study based on a sample of orbits for which the orbital poles were determined unambiguously. This research (Dommanget 1968), updated in 1982–1987 (Dommanget 1988), confirmed Mädler’s suspicion that the orbital poles of visual binaries in the solar neighbourhood do not seem to be oriented isotropically, but rather seem to cluster around the positions (l,b) ~ (100°, −15°) and (280°, +15°). This claim was based on a limited number of objects (8 systems within 10 pc of the Sun and 70 systems up to 25 pc), and was not subject to any analysis of its statistical significance. It was later confirmed by Glebocki (2000), who found that the poles of 19 systems within 10 pc of the Sun are concentrated around positions near (l,b) ~ (110°, −10°) and (−115°, −5°), although the global distribution of 252 systems is isotropic. An historical overview of all these studies may be found in Dommanget (2014).
In this paper, we present the results of a new research conducted since 2005 within the framework of a collaboration between professional astronomers and amateur astronomers that are members of the Double Stars Committee of the Société Astronomique de France^{1}.
Following the suggestion of Dommanget & Nys (2006), we reinvestigated the question by concentrating on a very local data set. Moreover, to the best of our knowledge, it seems that the most recent list of orbital poles dates back to Batten (1967). Because the more recent studies did not publish the data sets on which they were based and because of the rapid increase of the number of available orbits, we feel the necessity to collect such data, and provide an updated list of orbital poles for binaries within 18 pc of the Sun.
As already noted by Finsen (1933) and Dommanget (1968, 1988, 2005), a severe pitfall awaits the researcher on this topic: the ambiguity in the position of the ascending node, which, if not correctly assessed (Sect. 2), may ruin the distribution of the poles on the sphere.
To avoid this problem and establish a valid statistical material, we have first a) assembled a list of known visual binary orbits in the vicinity of the Sun for which a set of consistent and accurate radialvelocity measurements is available for solving the ambiguity of the orbital ascending node; and then b) computed the Galactic positions of the orbital poles of all systems for which the ambiguity of the orbital ascending node has been solved.
At first sight, some largescale coherence in the process of binary formation might imply a nonrandom orientation of the poles on the celestial sphere. However, such a largescale coherence is not expected to last after the coeval association has dissolved, since at any given time in the history of the Galaxy, stars in a given volume emanate from very different birth places and times, and their presence in a given volume is a momentary event. Therefore, given these potentially important (and problematic) consequences, any deviation from isotropy in the pole distribution should be firmly established before embarking on any physical discussion on the above questions.
In this new research, emphasis has been place on

increasing the sample size with regard to previous studies;

using dedicated tests to evaluate the statistical significance of any deviation from isotropy in the distribution of the orbital poles;

analysing the evolution, as a function of the volume sampled, of such possible deviations from isotropy;

discussing these results from the standpoint of Galactic kinematics.
This paper is structured as follows: in Sect. 2, we briefly recall the problem of unambiguously determining of the direction of the orbital pole, all mathematical details being provided in Appendix B. The master list of selected systems is presented in Sect. 3. For 51 systems up to 18 pc from the Sun, the distribution of the orbital poles on the sky is presented in Sect. 4. This section describes the available statistical tools used to analyse this distribution, and some indications about its degree of anisotropy are given. Possible selection effects in our sample is discussed in Sect. 5. The problem of possible orbital coplanarity in multiple systems is briefly discussed in Sect. 6. Section 7 discusses the question of the orientation of the orbital poles in the framework of Galactic kinematics. Conclusions and perspectives are given in Sect. 8. In Appendix A, notes are given on individual systems, along with new radialvelocity measurements used in the present study and, in several cases, new combined astrometricspectroscopic orbits that resulted. In Appendix B, we recall the conventions adopted to define the relative orbit of a binary star, and we describe the method used to determine the ascending node and the direction of the orbital pole without any ambiguity.
2. Unambiguous determination of the orbital pole of a binary
The spatial orientation of the orbital planes of double stars is studied based on the knowledge of the orientation of their orbital poles.
It is well known that an essential problem in this type of study is that the oriention of the orbital pole is ambiguous for visual and astrometric binaries because two true orbits are compatible with the apparent orbit calculated from the astrometric observations: although the orbital inclination i is unambiguously determined, it is impossible to distinguish between the two possible values of the position angle Ω of the ascending node. The problem can only be solved if radialvelocity measurements are available for at least one component of the system. Then, two scenarios may occur:

i)
for spectroscopic binaries resolved as visual or astrometric, the ascending node is immediately known unambiguously;

ii)
for visual or astrometry binaries for which radialvelocity measurements are available, the correct ascending node may be selected by comparing the trend of the measured radial velocities with the slope of the relative radial velocity curve computed from the orbit.
The conventions used in defining of the relative orbit of a binary as well as the details on the method used to solve the ambiguity of the ascending node and computating of the orientation of the orbital pole in Galactic coordinates are given in the Appendix B.
3. Data
The selection of the 95 systems in the master sample presented in Tables 1 and 2 was based on the availability of a visual or an astrometric orbit in the 6th Catalog of Orbit of Visual Binary Stars at USNO (6th COVBS^{2}, Hartkopf et al. 2001) and of a Hipparcos parallax (ESA 1997) larger than 50 mas^{3}. Table 1 lists various identifiers for the systems, such as their designations in the WDS (J2000 coordinates and the component designation taken in the 6th COVBS), CCDM, HIP or HIC, HD and SB9 (the 9th Catalogue of Spectroscopic Binary Orbits^{4}), catalogs, and the discoverer name. Table 2 lists various physical data for the systems identified by their WDS designation (Col. 1). The parallax (Col. 2) is for all HIP entries from the revised Hipparcos Catalogue (van Leeuwen 2007), or from the paper listing the orbit otherwise. The spectral type (Col. 3) is generally taken from the Simbad database at the Centre de Données Astronomiques de Strasbourg (CDS; Genova et al. 2000). Columns 4 (orbital period), 5 (orbital semimajor axis), 6 (total mass of the system), 7 (mass ratio), 8 (Grade of the orbit) and 9 (reference of visual orbit) come from the 6th COVBS, or are directly derived from the catalogue data (the total mass of each system – Col. 6 – is computed from the orbital elements and the parallax using the third Kepler law), with the exception of the mass ratio (Col. 7), which is listed when available from SB2 or astrometry. Column 10 (“S.Orb”) lists either the reference of the spectroscopic orbit used to lift the ascendingnode ambiguity, or “New” if the spectroscopic orbit has been computed in the present paper (see Appendix A), or “VD” if only the radialvelocity drift was used to lift the ambiguity, or “SD” if the sign of the velocity difference V_{r}(B) – V_{r}(A) was used to lift the ambiguity. Column 11, labelled “A”, contains a flag indicating whether the system has (or not) been considered in the final analysis: Y stands for yes, NO stands for “no reliable visual orbit”, NV stands for “no radial velocities”, ND stands for “no visible drift in the velocities”, because either the period is too long or the velocities too inaccurate, and IC stands for inconclusive analysis. For systems with a combined visual and spectroscopic solution computed by Pourbaix (2000), the orbital elements were adopted without an additional verification because the method naturally lifts the ambiguity on the ascending node. The last column of Table 2 identifies systems with an explanatory note (“n”) in Appendix A and triple or highermultiplicity systems noted (“t”).
Radialvelocity measurements were searched for in the literature and among unpublished measurements obtained with the CORAVEL spectrovelocimeter (Baranne et al. 1979). A few more radial velocities were obtained especially for this purpose with the HERMES spectrograph installed on the Mercator 1.2 m telescope (Raskin et al. 2011). In the end, 51 systems from the master list comprising 95 binaries (Table 2) have radial velocities available that enable lifting the ambiguity on the ascending node. Hence the number of systems used in our analysis drops from 95 to 51, and the systems analysed all have d< 18 pc.
For each system, the ambiguity on the ascending node was solved using one of the methods described in Appendix B.3. The orbital pole orientation was been computed following the method of Appendix B.4. In several cases, the radial velocities collected from the HERMES spectrograph (Raskin et al. 2011) or from the CORAVEL database (Table A.1) enable computing a new or updated spectroscopic binary orbit, as listed in Tables A.2 and A.3. These solutions were computed by combining the astrometric and spectroscopic data, following the method described by Pourbaix (2000). The orbital solutions are drawn in Figs. A.2–A.8.
The positions of the orbital poles are listed in Table 3. The column “Notes” identifies the entries for which a comment is provided in Appendix A. References to visual and spectroscopic orbits listed in Table 2: A1918b : Aitken (1918); Abt2006: Abt & Willmarth (2006); AST1999: Franz et al. (1999); AST2001: Benedict et al. (2001); B1960c: van den Bos (1960); Bag2005: Balega et al. (2005); Baz1980b: Baize (1980); Chg1972: Chang (1972); CIA2010: Farrington et al. (2010); Cou1960b: Couteau (1960); Doc1985c: Docobo & Costa (1985); Doc2008d: Docobo et al. (2008); Doc2010d: Docobo & Ling (2010); Dom1978: Dommanget (1978); Dru1995: Drummond et al. (1995); Duq1991b: Duquennoy et al. (1991); Egg1956: Eggen (1956); Egn2008: Eggenberger et al. (2008); Fek1983: Fekel & Tomkin (1983); FMR2008b: Rica Romero (2008); Frv1999: Forveille et al. (1999); Gri1975: Griffin & Emerson (1975); Gri1998: Griffin (1998); Gri2004: Griffin (2004); Grr2000: Girard et al. (2000); HaI2002: Han & Gatewood (2002); Hei1974c: Heintz (1974); Hei1984a: Heintz (1984); Hei1986: Heintz (1986); Hei1987b: Heintz (1987); Hei1988d: Heintz (1988); Hei1990c: Heintz (1990a); Hei1990d: Heintz (1990b); Hei1994a: Heintz (1994); Hei1994b: Heintz & Cantor (1994); Hei1996: Heintz (1996); HIP1997d: ESA (1997); Hle1994: Hale (1994); Hnk2011: Hinkley et al. (2011); Hop1958: Hopmann (1958); Hop1973a: Hopmann (1973); Hrt1996a: Hartkopf et al. (1996); Hrt2003: Hartkopf & Mason (2003); Irw1996: Irwin et al. (1996); Jnc2005: Jancart et al. (2005); Kam1989: Kamper et al. (1989); Kiy2001: Kiyaeva et al. (2001); Knc2010: Konacki et al. (2010); Lip1967: Lippincott (1967); Lip1972: Lippincott (1972); Llo2007: Martinache et al. (2007); Maz2001: Mazeh et al. (2001); Mnt2000a: Mante (2000); Msn1995: Mason et al. (1995); Msn1999a: Mason et al. (1999); Msn2011c: Mason & Hartkopf (2011); Mut2010b: Muterspaugh et al. (2010a); Mut2010c: Muterspaugh et al. (2010b); Nid2002: Nidever et al. (2002); Pbx2000b: Pourbaix (2000); Pbx2002: Pourbaix et al. (2002); PkO2006b: Gorshanov et al. (2006); Pop1996b: Popovic & Pavlovic (1996); Sca2002c: Scardia (2002); Sca2007c: Scardia et al. (2007); Sgr2000: Ségransan et al. (2000); Sgr2012: Barry et al. (2012); Sod1999: Söderhjelm (1999); Str1969a: Strand (1969); Str1977: Strand (1977); USN1988a: Dahn et al. (1988); USN1988b: Geyer et al. (1988); USN2002: Seymour et al. (2002); vAb1957: van Albada (1957); Wie1957: WiethKnudsen (1957); Zir2003: Zirm (2003); Zir2011: Zirm (2011).
Fig. 1 Positions of the orbital poles on the sky, with the grid representing Galactic coordinates. The Galactic centre is at the centre of the grid. The systems are labelled by their numbers of Table 3. The circled labels identify the systems up to N = 18 whose distribution is the most significantly different from isotropic according to Fig. 3. The dashed line displays the celestial equator. For comparison, the ecliptic pole is located at l = 96.4°,b = 30°. 
4. Distribution of the orbital poles on the sky
Figure 1 displays the distribution on the sky in galactic coordinates of the orbital poles listed in Table 3, for systems up to 18 pc. After estimating in Sect. 4.1 the typical errors affecting the orbitalpole positions, we describe in Sect. 4.2 the available statistical tools for studying the possible deviation from isotropy of the distribution of orbital poles on the sky.
4.1. Estimating the error on the orbitalpole positions
To evaluate the typical error affecting the orbitalpole positions, we used a Monte Carlo approach because of the complexity of the algebraic relations that link the orbital elements and the system position to the orbitalpole position (Eqs. (B.4) to (B.7)). We used the system WDS 00057+4549 (HD 38) as a test case because it is one of the systems with the longestperiod for which we lifted the ascendingnode ambiguity. Even though the error on the orbital period is large (P = 509.6 ± 100.0 yr), only the uncertainties on the inclination i and on the longitude of the ascending node Ω affect the uncertainty on the pole position. For WDS 00057+4549, we therefore adopted: α = 1.4209°, δ = 45.81°, i = 54.9° ± 2.4°, and Ω = 13.5° ± 2.3°, according to the orbit computed by Kiyaeva et al. (2001). We drew 10^{4} (i,Ω) pairs from this average and standarddeviation values, assuming that they follow a normal distribution; we then derived the orbital pole position for each of these cases from Eqs. (B.4) to (B.7). The results are displayed in Fig. 2. The Gaussian fits to the resulting longitude and latitude distributions (bottom panel), show – as expected – that the uncertainties on i and Ω directly propagate onto the polar position (2.3° and 2.0° in Galactic longitude and latitude, respectively, for WDS 00057+4549). If the uncertainties on i and Ω are unequal, the resulting distribution of the orbital poles will simply be tilted with respect to the axis frame. In the analysis of the pole distribution presented below, one should therefore keep in mind that these positions are subject to an uncertainty of the same order as that on i and Ω.
Fig. 2 Top panel: distribution of the orbitalpole positions on the sky, drawn from the normal distributions for i = 54.9° ± 2.4° and Ω = 13.5° ± 2.3° corresponding to the benchmark system WDS 00057+4549. The intersection of the dashed lines marks the nominal pole position. As expected, the distribution is centred on this nominal position. Bottom panel: histograms corresponding to the longitude and latitude distributions. The dashed lines correspond to Gaussian distributions centred on the nominal pole position, with standard deviations of 2.3° and 2.0° in Galactic longitude and latitude, respectively. 
4.2. Available statistical tools
The current section is mainly based on Gillett (1988), Upton & Fingleton (1989), and Briggs (1993). In this section, we discuss how to test the null hypothesis that the pole distribution is isotropic. The simplest test is the Rayleigh test (Fisher 1953; Briggs 1993). For the purpose of Rayleigh statistics, the i = 1, ...N poles are represented by unit vectors (x_{i},y_{i},z_{i}) (with ), and one computes (1)where ℛ representing the module of the vectorial sum of the N unit vectors (x_{i},y_{i},z_{i}). This vectorial sum gives the orientation of the dipole of the distribution. If the poles are isotropically distributed, the unit vectors representing their orientation will tend to cancel each other out and the resulting ℛ will be close to zero, while if the poles have a preferential orientation, ℛ will tend to be larger than expected for the isotropic distribution. However, there could be very specific anisotropic situations leading to low ℛ values, for instance one with an even number of poles distributed in a bimodal distribution with the modes aligned in exactly opposite directions, and equal numbers of poles populating both modes. This underlines the fact that a test such as the Rayleigh test should always be carried out against an alternative model, and this model must be chosen on physical grounds (Gillett 1988). The most natural alternative to be used with Rayleigh’s test is the Fisher distribution (Fisher 1953; Watson 1956), which is unimodal and assumes symmetry about the mean. Obviously, this bimodal situation does not follow a Fisher distribution, and Rayleigh’s test has therefore no power to distinguish it from isotropy (i.e., the secondkind risk of accepting the null hypothesis of isotropy while it is false is high). The combined statistic of Giné (1975), which are described below (Eq. (17)), is the most powerful against all alternatives (Gillett 1988).
The Rayleigh statistics may be used as a test of isotropy against a unimodal distribution by calculating the probability α_{ℛ} of obtaining a value of ℛ greater than or equal to the observed value ℛ_{obs} assuming that the data are drawn from the isotropic distribution (Watson 1956, 1983): (2)where f_{ℛ} is the probability density function for ℛ (i.e., the probability of finding ℛ between R and R + dR is f_{ℛ}dR). Its detailed expression may be found in Fisher (1953) and Briggs (1993). If α_{ℛ} is small, then the poles are very unlikely to have been drawn from an isotropic parent population. Under such circumstances, it is said that the observed pole distribution is inconsistent with isotropy at the 1 − α_{ℛ} confidence level. It means that there is a (firstkind) risk α_{ℛ} of rejecting the null hypothesis of isotropy while it is actually true. Stephens (1964) provides a table of critical values for ℛ for α_{ℛ} = 1,2,5 and 10%.
While the Rayleigh statistics is valid for any value of N, there exists a more convenient form of it at large N (≥50), known as the RayleighWatson statistics (Watson 1956, 1983): (3)Upton & Fingleton (1989) used the table of Stephens (1964) to devise an improved approximation, which works well even for low values of N: (4) follows a χ^{2} distribution with three degrees of freedom. For convenience, the values of the correction factor c(α) for the Rayleigh test are reproduced from Upton & Fingleton (1989) in Table 4. A value indicates that there is a significant concentration of poles toward one direction (with a confidence level of 99%).
Values of the correction factor c(α) for the Rayleigh test, from Upton & Fingleton (1989).
Instead of considering ℛ as the length of the resulting dipole vector, it may be useful to consider it as the distance from the origin achieved by a random walk of N unit steps. A wellknown result from the randowwalk theory then states that , from which it follows that . It is therefore not surprising that the statistic is asymptotically distributed as a χ^{2} distribution with three degrees of freedom (Watson 1956, 1983), so that (5)where is the probability density of and is the upper cumulative probability distribution function of . Obviously, the latter function is much easier to evaluate than the corresponding one for the Rayleigh statistics (Eq. (2)). An additional advantage of this statistics is that it is asymptotically independent from N.
Again, the isotropy of the data may be tested by calculating the observed statistics and the corresponding probability , which has a similar interpretation as the one given above for α_{ℛ}. For large N, , and the two tests are equivalent.
We stress here an important result from the hypothesistesting theory: although the firstkind risk (of rejecting the null hypothesis of isotropy while it is true, also known as the falsealarm probability) is controlled by the adopted confidence level 1−α_{ℛ}, the secondkind risk (of accepting the null hypothesis of isotropy while it is false, also known as the power of the test) is not controlled as easily, since this secondkind risk depends on the exact nature of the true underlying distribution. This was well illustrated above with the Rayleigh test which is said to have no power against a bimodal distribution. Because different statistical tests have different powers it is important not to use only one if one does not know the actual underlying distribution.
Since the Rayleigh test is not appropriate when the alternative is a bimodal axial distribution, for testing uniformity when the alternative is such, we must turn to tests based upon the properties of statistics involving the eigenvalues of the quadrupolar matrix M_{N} (Watson 1966, 1983; Briggs 1993): (6)Since M_{N} is real and symmetric, it has three eigenvalues λ_{k}, and since it has unit trace, the sum of its eigenvalues λ_{k} is one^{5}. Moreover, the diagonal elements are sums of squares. In the coordinate system in which M_{N} is diagonal, the diagonal elements are the eigenvalues, and therefore the eigenvalues are not negative. If the poles were drawn from the isotropic distribution, then, excepting statistical fluctuations, the eigenvalues λ_{k} should be equal, and since M_{N} has unit trace, λ_{k} = 1 / 3.
If the alternative to isotropy is a rotationally asymmetric axial distribution^{6}, the Bingham statistics ℬ is appropriate: it measures the deviation of the eigenvalues λ_{k} of the quadrupolar matrix M_{N} from the value 1/3 expected for isotropy (Bingham 1974; Watson 1983): (7)The statistics ℬ is asymptotically distributed as a χ^{2} distribution with five degrees of freedom, so that (8)where is the probability density function of and is the upper cumulative probability distribution function of . Again, if α_{ℬ} is very small, the hypothesis of isotropy is contradicted and rejected with a firstkind risk of α_{ℬ}. Briggs (1993) showed that α_{ℬ} as derived from Eq. (8) overestimates the true firstkind risk when N< 20 for α_{ℬ} = 0.05, and when N< 40 for α_{ℬ} = 0.01. In other words, under these circumstances, 1−α_{ℬ} will underestimate the statistical confidence with which the hypothesis of isotropy can be rejected.
If the alternative to isotropy is a rotationally symmetric axial distribution (either a bipolar distribution along a given direction, or a girdle with rotational symmetry, both known as the DimrothWatson distribution; Upton & Fingleton 1989), the Bingham test is still applicable, but two of the three eigenvalues will then be close to each other (Bingham 1974, and below). Assuming that the pole of the distribution is located at φ_{0} = 0° and θ_{0} = 0° ((φ,θ) being a polar coordinate system on the sphere, with φ the azimutal angle, and θ the polar angle), Mardia & Jupp (2000) showed that , the maximum likelihood estimate of the polar direction, is related to the eigenvectors of the M_{N} matrix; a girdle distribution will have its symmetry axis along e_{1}: (9)whereas an axial distribution will have e_{3} as symmetry axis^{7}: (10)In the general case of a rotationally asymmetric axial distribution, one has λ_{2} − λ_{1} and λ_{3} − λ_{1} large and nonequal (girdle), or λ_{3} − λ_{1} and λ_{3} − λ_{2} large and nonequal (polar).
Anderson & Stephens (1972) showed that λ_{1} and λ_{3} themselves are the most appropriate test statistics for the comparison of the hypothesis of uniformity against that of a DimrothWatson distribution. Upton & Fingleton (1989) have used the tables computed by Anderson & Stephens (1972) to devise critical values of λ_{1} and λ_{3} at confidence level α. A significant departure from isotropy towards the girdle form of the DimrothWatson distribution is judged to have occurred ifλ_{1}<λ_{1}(α), where (11)whereas a significant departure towards the bipolar form of the distribution is judged to have occurred if λ_{3}>λ_{3}(α), where (12)where the coefficients a_{1,3}(α), b_{1,3}(α), and c_{1}(α) have values that depend on the significance level, α. These values may be found in Table 10.15 of Upton & Fingleton (1989). The values λ_{1}(α) and λ_{3}(α) have been plotted as dashed lines in the lower panel of Fig. 3 for α = 1% and 10%.
Fig. 3 Upper panel: evolution of the RayleighWatson (blue dashed line, and righthand scale; the correction to – Eq. (4) – has been made for a confidence level α of 1%) and Bingham ℬ (red dotdashed line, and righthand scale) statistics as a function of the number of systems. The horizontal blue dashed curve corresponds to a firstkind risk of 1% of rejecting the null hypothesis of isotropy while it is true, based on the statistics. The horizontal red dashed line is the 10% firstkind risk for the ℬ statistics. The number of systems as a function of distance is given by the black solid line and lefthand scale. Isotropy corresponds to ℬ → 0 and . Lower panel: three eigenvalues of the orientation matrix (Eq. (6)). Isotropy corresponds to λ_{1,2,3} = 1 / 3, as seen at large distances. The critical values for λ_{1,3} corresponding to α = 1% and 10% as given by Eqs. (11) and (12) are displayed by the dashdotted lines. 
In summary, by combining the RayleighWatson statistics with the Bingham statistic, the isotropy of the distribution may be distinguished among several alternatives, since a value reveals that there is a concentration of poles toward one direction, but values close to 3 indicate that the distribution is either isotropic or bipolar. The bipolarity of the distribution is in turn assessed from λ_{3}>λ_{3}(α). To summarize (see also Table 1 of Briggs 1993; and Table 10.7 of Upton & Figgleton 1989):

unipolar distribution around one end of the eigenvector e_{3}:

bipolar distribution around both ends of the eigenvector e_{3}:

girdle orthogonal to the eigenvector e_{1}:
(15)Finally, one should mention the nonparametric statistical tests that do not presuppose any alternative model. They are thus more general than tests against a specific alternative, but for that reason are also less powerful (i.e., have a higher risk of accepting the null hypothesis of isotropy while it is false; Gillett 1988; Mardia & Jupp 2000). To test against alternatives that are asymmetric with respect to the centre of the sphere, Beran’s B test (Beran 1968) is the most powerful: (16)To test against symmetric or axial alternatives, Giné’s test G (Giné 1975) is the most powerful: (17)The combined Giné statistic F = B + G tests against all alternatives. In Eqs. (16) and (17), N is the number of data points, and Ψ_{ij} = cos^{1}(x_{i}·x_{j}) is the smaller of the angles between poles x_{i} and x_{j}. High values of these statistics suggest departure from isotropy. Critical values for various percentage levels of significance for B,G and F have been calculated by Keilson et al. (1983). In the largesample limiting distribution of F, the 10%, 5%, and 1% confidence levels are 2.355, 2.748, and 3.633, respectively.
Various isotropy indicators, listed separately for the 20 systems up to 8.1 pc, and for all 51 systems up to 18 pc.
4.3. Isotropic or anisotropic?
In this section, we apply the methods outlined in Sect. 4.2 to the sample of orbital poles listed in Table 3. The top panel of Fig. 3 displays the evolution of the RayleighWatson and Bingham ℬ statistics as a function of distance, thus for an increasing number of systems.
The various statistical indicators described in Sect. 4.2 are listed in Table 5 separately for the 20 systems up to 8.1 pc, and for all 51 systems closer than 18 pc. For estimating and α_{ℬ}, the asymptotically equivalent and distributions were used, respectively.
For the 20 systems with d ≤ 8.1 pc, the RayleighWatson test signals a very significant deviation from isotropy^{8}, with a (firstkind) risk of % of rejecting the null hypothesis of isotropy while it is actually true (Table 5). The Beran test (for isotropy against alternatives that are asymmetric with respect to the centre of the sphere) provides an identically significant result (B = 3.32 or α ~ 0.7%; Keilson et al. 1983).
As seen from Table 5, the Bingham and Giné tests are less conclusive, but this is easily understood. The Giné test, being “universal” (it tests against all alternatives), is less efficient than RayleighWatson’s in the sense that it requires a stronger deviation from isotropy to yield similarly low α values. Bingham test being devised against symmetric alternatives, which is not the situation encountered in our data (Fig. 1), it is not efficient either in the situation under study, and should not be given much consideration.
Another way to understand why the RayleighWatson and Bingham statistics yield vastly different significance levels is remembering that the RayleighWatson and Bingham statistics test for the presence of a dipole or of a quadrupole, respectively, in the distribution of poles on the sphere (Hartmann & Epstein 1989; Briggs 1993, and Sect. 4.2). As already discussed in relation with Eq. (1), the vector ℛ is the dipole moment of the distribution, whereas the eigenvectors of the quadrupole matrix M_{N} are related to the orientation of the quadrupolar component of the distribution (see the discussion in relation with Eqs. (13)–(15)). Therefore, the dipolar vector and the eigenvector e_{3} associated with the highest eigenvalue of M_{N} ought be different. For the considered orbitalpole distribution, the largest dipole moment is obtained for the N = 20 closest systems and points in the direction l = 46.0°, b = 36.7°, with a clear deficit of poles in the opposite direction. On the other hand, the eigenvector e_{3} associated with the highest eigenvalue λ_{3} = 0.52 for N = 20 points towards l = 77.4°, b = 8.4°. These directions are similar, albeit not identical, and obviously point in the direction of a higher concentration of poles (Fig. 1).
The RayleighWatson (dipolar) test yields a deviation from isotropy that is more significant than the Bingham (quadrupolar) test. Figure 3 indeed signals that the eigenvalues λ_{1,2,3} and fulfill Eq. (13) that corresponds to a unipolar distribution.
To evaluate the robustness of the low firstkind risk obtained for the RayleighWatson statistics, a jackknife approach was used, that is, we repeated the above procedure of computing as a function of distance (or sample size, with systems ordered by increasing distance from the Sun) after removing one system at a time from the full sample. Since peaks around N = 20, it was sufficient to do so for the 25 systems closest to the Sun, removing one system at a time. This way, 25 curves of as a function of N were generated, one for each of the 25 samples thus generated. These curves are presented in Fig. 4 and are similar to the curve presented in the upper panel of Fig. 3. The thick line in Fig. 4 corresponds to all systems being considered (thus identical to the curve in the upper panel of Fig. 3). The lower (red) line is the one corresponding to the sample with system No. 16 removed (since this system is the one weighting heavily in the original peak; the number refers to Table 3 and Fig. 1). As expected, this curve has the lowest significance, peaking at only 1.5% instead of 0.5% in the full sample. The other thin curves are for all the other samples, as explained above. Figure 4 thus provides a fair evaluation of the robustness of the alarm signal: indeed, after removing system No. 16, the falsealarm probability just misses 1% and is thus not especially significant; but removing closer systems (such as Nos. 1 or 2, i.e., systems far away from the pole cluster seen in Fig. 1) makes the curve rocket to 0.1% significance! In summary, the fluctuations caused by the small sample size causes the highest significance level to vary between 1.5 and 0.1%! The reality of the deviation from isotropy can thus not be assessed with certainty at this stage, given the small number of systems available, and despite our efforts to increase it.
The fact that this deviation from isotropy does not appear for subsamples with cutoff distances shorter than 8 pc moreover decreases its physical significance. One could have envisioned for instance a situation where the index maintains statistically significant values all the way to 8 pc, and then fades to insignificant values at larger distances. But the first part of that statement is not fulfilled: the closest 15 systems (up to 6 pc) have % (Fig. 3). At distances larger than 8 pc, the significance of the deviation from isotropy fades away again, at least for the Bingham statistics: ℬ → 0, λ_{1,2,3} → 1 / 3. The rise of ℬ above 2.5 for N> 47 is not at all significant since it corresponds to α_{ℬ} = 0.79. In contrast, the RayleighWatson statistics does not drop to its isotropic value of 3, but rather stays around 4, corresponding to a firstrisk probability of about 18% of rejecting the null hypothesis of isotropy despite being true. Indeed, the addition of systems farther away than 8 pc does not fully destroy the anisotropy that appeared in the 8 pc sphere, since a deficit of poles remain around l = 270°, b = −30° (Fig. 1).
We thus conclude that the poles of the binary systems within 8 pc of the Sun exhibit a weak tendency to cluster around l = 46°, b = 37° (which, incidentally, is not far from the pole of the ecliptic). The systems numbered 15–18 and 20 (see Table 3) are responsible for that clustering, which fades away if shorter or larger sampling distances are considered.
Fig. 4 Evolution of the RayleighWatson statistics as a function of the number of systems in a jackknife approach. Each curve thus corresponds to the full sample with one system removed. The thick line corresponds to the full sample (same as the blue dashed line in the upper panel of Fig. 3). The lowest thin (red) line corresponds to the sample with system No. 16 removed (WDS 14514+1906). Falsealarm probabilities of 2.5, 1, and 0.1% are depicted by the dashed horizontal lines. 
5. Impact of selection effects
Fig. 5 Location of binary systems from Table 2 in Galactic coordinates, with the galactic centre in the middle of the map. The circles identify the systems with an unambiguously defined pole, as plotted in Fig. 1. The dashed line is the celestial equator. 
As Blaauw appropriately commented after Batten’s talk on the topic of the present paper during the IAU symposium 30 (Batten 1967), the possibility that observational selection effects could actually be causing an anisotropy in the pole distribution should be considered. Figure 5 indeed reveals that the studied systems are concentrated in the northern equatorial hemisphere, since lifting the pole ambiguity requires radial velocities that are easier to obtain for northern hemisphere targets (combined with the fact that the astronomy of visual double stars relies on a longer timebase in the northern hemisphere; hence long periods are more commonly derived for northern targets). Could this anisotropic initial distribution be the cause for the pole anisotropy observed in the 8 pc sphere around the Sun?
For this to be the case, there should be a correlation between the position of the binary on the sky and the position of its orbital pole. A nonuniform celestial distribution of the systems studied (as observed in Fig. 5) would then propagate onto their orbitalpole distribution. A correlation between system position and pole position would manifest itself as a nonuniform distribution of the scalar product between the unit vector defining the system position (denoted 1_{∗} hereafter) and the unit vector defining the pole position (denoted 1_{p} hereafter). From Eq. (B.4), it is easy to verify that this scalar product 1_{∗}·1_{p} equals −cosi. Thus, an observational bias acting on our ability to derive the poles for systems with a given range of inclinations i would introduce a correlation between 1_{∗} and 1_{p}. To check whether this is the case, we plot in Fig. 6 the distribution of cosi values for all systems of Table 2 and for those with the poles analysed (thin line and shaded histogram). There is no significant difference between these distributions, meaning that the requirement of having radial velocities available for lifting the pole ambiguity does not introduce a selection bias on the orbital inclination i. On the other hand, the fullsample and poleknown distributions of cosi are close to being uniform^{9} except for two marginal peaks at 0.2 and 0.7. However, given the small amplitude of these peaks, it seems very unlikely that they could be at the origin of the segregation observed in the pole distribution. Furthermore, there is no obvious observational selection effect that we could think of that could cause peaks at cosi ~ 0.2 and 0.7.
Fig. 6 Distribution of cosi for all binaries from Table 2 (thick lines) and for those with the poles analysed (thin line in the bottom panel and shaded histogram in the upper panel). 
6. Coplanarity in multiple systems
When searching for some anisotropy in the space orientation of orbital poles, which could either be signatures of the process of binary formation or of tidal effects from the Galactic disc, one could be misled by evolutionary processes acting within the stellar system itself that modify the orbital inclination. Within a triple system, weak perturbations from the outer body can have strong longterm effects on the inner binary. The simplest of these is precession of the orbital plane, which occurs if the orbital planes of the inner and outer binary are not aligned (Fabrycky & Tremaine 2007). If the inner and outer binary orbits are circular, this precession is analogous to the precession of two rigid rings with the same mass, radius, and angular momentum as the binary orbits; both the mutual inclination and the scalar angular momenta of the rings remain fixed, while the two angular momentum vectors precess around the direction defined by the total angular momentum vector of the triple system.
An unexpected aspect of this behaviour was discovered by (Kozai 1962; see also Kiseleva et al. 1998; Eggleton & KiselevaEggleton 2001; Fabrycky & Tremaine 2007; Muterspaugh et al. 2008). Suppose the inner binary orbit is initially circular, with the initial mutual inclination between inner and outer binaries equal to i_{initial}. Kozai found that there is a critical angle i_{c} such that if i_{initial} is between i_{c} and 180°−i_{c}, then the orbit of the inner binary cannot remain circular as it precesses; both the eccentricity of the inner binary e_{in} and the mutual inclination i execute periodic oscillations known as Kozai cycles. The amplitude of the eccentricity and inclination oscillations is independent of the strength of the perturbation from the outer body (which depends on the mass of that outer body, and on a_{out} and e_{out}), but the oscillation amplitude does depend on i_{initial}; for initial circular orbits with i_{initial} = i_{c} or 180°−i_{c}, the eccentricity is 0, but if i_{initial} = 90°, the highest eccentricity is unity, that is, the two inner bodies collide.
In the framework of the present study, it is important to stress that Kozai cycles occur independent of distances or component masses, their only requirement is that the mutual inclination is between i_{c} = 39.2° and 180°−39.2° = 140.8° (Fabrycky & Tremaine 2007; Muterspaugh et al. 2008).
Hence, for triple systems, the current orbital pole may not at all be representative of its initial orientation if the Kozai mechanism is operating. Therefore, a discussion of the pole orientation should properly identify the triple systems, and possibly exclude them from the analysis. The triple and highermultiplicity systems in our sample are identified by flag “t” in Tables 2 and 3, the notes in the Appendix provide more details in most cases.
Because the Kozai mechanism does alter the orbital inclinations of noncoplanar triple systems, we now investigate the poleanisotropy problem separately for the sample devoid of triple systems (Figs. 7 and 8). These figures reveal that the RayleighWatson test still predicts some deviation from isotropy in the sample restricted to 8 pc, but the firstkind error of the RayleighWatson statistics this time is not better than 5%. We thus conclude that triple systems in the analysed sample do not degrade a possible anisotropy, quite the contrary, since the firstkind error of rejecting the null hypothesis of isotropy even increases when including the triple systems.
Fig. 8 Same as Fig. 3, but restricted to binary systems without evidence for higher multiplicity. The horizontal blue dashed and red dashdotted curves correspond to a firstkind risk of 5% and 10%, respectively, of rejecting the null hypothesis of isotropy while it is true, based on the statistics. The systems have the same number labels as in Table 3 and Fig. 3. 
Nevertheless, the fact that several triple systems are now known not to be coplanar is certainly not in favour of a concentration of orbital poles on the celestial sphere. The quadruple system WDS 11182+3132 belongs to the sample of systems closer than 8 pc to the Sun; the AB and Aa,Ab poles are known (Nos. 18 and 19 in Fig. 1) and are separated by about 100° on the sky. Although such a complete knowledge is not very frequent^{10}, noncoplanarity is also observed in η Vir (HIP 60129: the orbital planes of the two systems are inclined by about 30° with respect to each other; Hummel et al. 2003), and similarly for the triple system V819 Her (HD 157482; O’Brien et al. 2011). In addition to V819 Her and η Vir, unambiguous mutual inclinations have been calculated for several other systems: κ Peg, ϵ Hya, ζ UMa, and Algol (see Table 5 of O’Brien et al. 2011, and references therein), and μ Ori, ξ UMa, ϵ Hya, and 88 Tau (Muterspaugh et al. 2008). These studies, along with the earlier ones by Fekel (1981) and Sterzik & Tokovinin (2002), strongly argued against coplanarity, despite predictions from older theories of formation of multiple stars (Bodenheimer 1978). Modern theories of stellar formation (Sterzik & Tokovinin 2002; Fabrycky & Tremaine 2007, and references therein) indeed predicted that triple systems will form noncoplanarily. Clearly, this result jeopardizes the possibility of observing coplanarity on a larger scale such as that of the solar neighbourhood. Therefore, the weak evidence for pole alignment observed in a sphere of 8 pc around the Sun does not receive support from our analysis of the properties of triple systems. In the next section, we show that neither do theories of binarystar formation coupled with the dispersive effect caused by Galactic rotation provide supportive arguments for a possible alignment of the orbital poles in the solar neighbourhood.
7. Pole alignment in the context of binarystar formation theories and Galactic kinematics
This section addresses the following two questions: (i) Do current theories of binarystar formation predict alignment of orbital poles for all binaries resulting from the same event of star formation; and (ii) assuming that the answer to the first question be positive, is this coherence in the pole orientation of young binary systems preserved by the Galactic kinematics in a volumelimited sample (the solar neighbourhood) consisting of much older systems?
Theories for wide and close binarystar formation have been reviewed by Bodenheimer & Burkert (2001) and Bonnell (2001), respectively. The favoured mechanism for producing most binary stellar systems is the fragmentation of a molecular cloud core during its gravitational collapse. Fragmentation can be divided into two main classes: direct fragmentation and rotational fragmentation (see references in Bate 2000). Direct fragmentation depends critically on the initial density structure within the molecular cloud core (e.g., nonspherical shape or density perturbations are needed to trigger the collapse), whereas rotational fragmentation is relatively independent of the initial density structure of the cloud because the fragmentation occurs due to nonaxisymmetric instabilities in a massive, rotationally supported disc or ring. Not many studies of binarystar formation however include a discussion about the possible spatial alignment of the orbital poles of binary systems resulting from the same fragmentation event. This lack of information results from the difficulty of conducting the numerical simulations till the end of the binary formation, since a phase of accretion follows the initial fragmentation of the molecular cloud. To obtain the final parameters of a stellar system, the calculation must thus be performed until all of the original cloud material has been accumulated by one of the protostars or by their discs. Following simultaneously several among such events is an even more challenging task thus. Although Bate (2000) specifically investigated the properties of binary systems after this phase of accretion, there is no discussion whatsoever of the question of orbitalpole alignment. This was addressed in a very qualitative way by Boss (1988). A short discussion about the question of pole alignment has been found in the context of a variant of the fragmentation scenario, namely the shockinduced gravitational fragmentation followed by capture (SGF+C). That mechanism was initially proposed by Pringle (1989). Collisions between stable molecular cloud clumps produce dense shocked layers, which cool radiatively and fragment gravitationally. The resulting fragments then condense to form protostellar discs, which at the same time collapse and, as a result of tidal and viscous interactions, capture one another to form binary systems. When the initial clumps are sufficiently massive, a large number (>10) of protostellar discs is produced. A prediction of the SGF+C mechanism (Turner et al. 1995) is that the resulting binary systems are likely to have their orbital and spin angular momenta aligned – the preferred direction of alignment being that of the global angular momentum of the original clump/clump pair. However, with more realistic, that is less highly organized, initial conditions, such alignments should still arise locally but are likely to be less well coordinated across large distances.
Even though there is currently no strong evidence from theories of binary formation for pole alignment in binary systems resulting from the same molecularcloud collapsing event, we assumed for now such a coherence and followed what happened as a result of Galactic kinematics dispersing these coeval binary systems across the Galaxy. After several Galactic rotations, this coeval association will (fully or partly) evaporate and form a tube called supercluster. Woolley (1961) has presented an interesting argument that makes it possible to tag stars that formed at the same place and time, and that are still found together today in the solar neighbourhood. In general, stars found today in the solar neighbourhood emanate from various birth locations in the Galaxy because the diversity of stellar orbits in the Galaxy. Woolley’s argument goes as follows. Stars in the supercluster mentioned above still share common V velocities (with the usual notation (U,V,W) for the Galactic velocity components directed towards the Galactic centre, along the direction of Galactic rotation, and towards the north Galactic pole, respectively) when located in the same region of the tube (for example in the solar neighbourhood) for the following reason, as put forward by Woolley (1961): if the present Galactocentric radius of a star on a quasicircular epicyclic orbit equals that of the Sun (denoted r_{⊙}), and if such a star is observed with a peculiar velocity v = V + V_{⊙}, then its guidingcentre radius r_{g} writes (18)where x_{g} is the position of its guidingcentre in the Cartesian reference frame associated with the local standard of rest (LSR) (in the solar neighbourhood approximation, the impact of y_{g} is negligible), and B is the second Oort constant. Woolley (1961) pointed out that disc stars (most of which move on quasicircular epicyclic orbits) that formed at the same place and time, and that because they stayed together in the Galaxy after a few galactic rotations (since they are all currently observed in the solar neighbourhood) must necessarily have the same period of revolution around the Galactic centre, and thus the same guidingcentre r_{g}, and thus the same velocity V = v − V_{⊙} according to Eq. (18).
Therefore, according to this argument, binary systems that formed in the same cluster or association and that stayed together after a few Galactic revolutions are observed today in the solar neighbourhood as sharing the same V velocity. In other words, if there are any orbitalpole concentrations on the celestial sphere that can be ascribed to angular momentum conservation for binary systems originating from the same cluster or association, they should be identified from similar V velocities.
We first computed the (U,V,W) velocity components from the Hipparcos proper motion and from the systemic velocity, according to standard prescriptions for converting equatorial to Galactic coordinates. The resulting distribution in the (U,V) plane is shown in Fig. 9.
Fig. 9 Distribution of the systems in the (U,V) plane, with the same labels as in Fig. 1. System 17 (01083+5455) falls outside the boundaries, at U = −43 km s^{1} and V = −157 km s^{1}. To fix the ideas, two velocity ellipsoïds are shown, both centred on the reflex solar motion in the LSR: −U_{⊙} = −10.2 km s^{1}, −V_{⊙} = −20.5 km s^{1}, and with σ_{U} = 34.5 km s^{1}, σ_{V} = 22.5 km s^{1}, as derived by Famaey et al. (2005) from their most precise sample. Horizontal bands and the corresponding labels are used in Fig. 10. 
As expected, the sample is distributed according to the general velocity ellipsoïd, with overdensities observed at the location of the Sirius (⟨ U ⟩ ~ 6.5 km s^{1}, ⟨ V ⟩ ~ 3.9 km s^{1}) and Hyades (⟨ U ⟩ ~ −40 km s^{1}, ⟨ V ⟩ ~ −20 km s^{1}) streams (e.g., Dehnen 1998; Famaey et al. 2005). Because of Woolley’s argument expressed above, the (U,V) diagram was divided into horizontal bands of width 10 km s^{1}, as indicated in Fig. 9, with a letter assigned to each of these Vbands. These letters were used to label the position of each pole on a Galacticcoordinate map (Fig. 10).
Fig. 10 Same as Fig. 1 but with a letter coding the value of the galactic V velocity (see Fig. 9). A necessary condition for binary systems to originate from the same cluster/association at birth is to have similar V velocities (see text), hence they should be labelled by the same letter. 
The absence of clustering of similar letters in restricted regions of the celestial sphere provides definite proof that any possible pole concentration on the sphere (see Sect. 4.3) cannot be traced back to a similar birth location for those systems. It must therefore either be ascribed to statistical fluctuations in our smallsized sample, or to some as yet unidentified physical process that removes orbital poles from the region 90° ≤ l ≤ 300°, b ≤ 0°, on a time scale of the order of that characterizing the convergence of systems with different Galactic orbits in the solar neighbourhood. Our Galaxy being in general symmetric with respect to its equatorial plane, it is however difficult to imagine a physical process that would break such a symmetry (the fact that the Sun lies somewhat above the Galactic plane? see below), as it seems to be required to account for the observed asymmetry.
Conversely, the poles of the binary systems belonging to the Sirius or Hyades streams do not concentrate either on the sphere, which is a further confirmation that these streams are not evaporated clusters, but rather have a dynamical origin (Famaey et al. 2005, 2007, 2008; Pompéia et al. 2011).
Finally, one should mention, for the sake of completeness, that the location of the Sun slightly above the Galactic plane (~20 pc; Humphreys & Larsen 1995) could be responsible for the (Galactic) North–South asymmetry observed in the orbitalpole distribution, if some effect related to Galactic tides were the cause of the pole anisotropy. That seems unlikely, however, given the low intensity of Galactic tides with respect to the twobody interaction. Heisler & Tremaine (1986) have expressed Galactic tides in terms of the Oort constants. Based on the size of the last closed Hill surface^{11}, they conclude that in the solar neighbourhood, the Galactic tides have a nonnegligible impact only on systems with semimajor axes of several 10^{4} AU (such as Oortcloud comets around the Sun). The binary systems defining the pole statistics are very much tighter than this, so that there seems to be no way for the Galactic tides to imprint their signature on the orbitalpole distribution.
8. Conclusion
Among the 95 systems with an orbit in the 6th Catalogue of Orbits of Visual Binaries that are closer than 18 pc from the Sun, we were able to lift the pole ambiguity for 51 systems, thanks to radialvelocity data either collected from the literature, obtained from the CORAVEL database, or acquired with the HERMES/Mercator spectrograph. Of these 51 systems, several had an erroneous node choice listed in the 6th USNO orbit catalogue and had thus to be corrected. An interesting side product of the present study is a number of new spectroscopic orbits obtained from CORAVEL and HERMES radial velocities, as well as new combined spectroscopic/astrometric orbits for seven systems (WDS 01083+5455aa,Ab; 01418+4237AB; 02278+0426AB (SB2); 09006+4147AB (SB2); 16413+3136AB; 17121+4540AB; 18070+3034AB).
This new sample of 51 orbital poles was subjected to several methods of spherical statistics (RayleighWatson, Bingham, Beran and Giné) to test possible deviations from isotropy. After ordering the binary systems by increasing distance from the Sun, the falsealarm probability (of rejecting the null hypothesis of isotropy while it is true) was computed for subsamples of increasing sizes, from N = 1 up to the full sample of 51 systems. Two tests (RayleighWatson’s and Beran’s) delivered a falsealarm probability of 0.5% when considering the subsample of 21 systems closer than 8.1 pc to the Sun. The poles tend to cluster toward Galactic position l = 46°, b = 37°. This direction is not so distant from that of the ecliptic pole (l = 96°, b = 30°). To evaluate the robustness of this clustering, a jackknife approach was used by repeating the above procedure after removing one system at a time from the full sample. The above falsealarm probability was then found to vary between 1.5% and 0.1%, depending on which system is removed from the sample. The reality of the deviation from isotropy can thus not be assessed with certainty at this stage, given the small number of systems available, despite our efforts to increase it. Our sofar uncertain conclusion should be seen as an incentive to foster further studies on this problem, especially in the Gaia era.
For the full sample extending up to 18 pc, the clustering totally vanishes (the RayleighWatson firstkind risk then rises to 18%). If present, coplanarity of orbits must thus be restricted to small spheres with a radius of the order of 10 pc, but several dynamical arguments were presented to show that the physical origin of any such coplanarity, should it exist, will be difficult to identify. In this context, the recent result of Rees & Zijlstra (2013) on the nonrandom orientation of bipolar planetary nebulae in the Galactic Bulge, using a statistical method much like the ones used in our present study, came as a surprise. Rees & Zijlstra found evidence for an anisotropy in the orientation of the bipolar planetary nebulae, whose specific geometry is usually attributed to the binarity of the central star. Assuming that the symmetry axis of the bipolar nebula coincides with the polar axis of the orbit, the result of Rees & Zijlstra is therefore indirect evidence of a deviation from anisotropy for the poles of a subsample of binaries in the Galactic Bulge. The physical explanation proposed for the Galactic Bulge (related to the impact of the Galactic magnetic field on the binary formation process) is, however, unlikely to hold also for the solar neighbourhood, where the field is of lower intensity.
Online material
Crossidentifications for the master list of 95 systems closer than 20 pc from the Sun and the are known to have an orbit in the 6th COVBS at USNO.
Data and references for the 95 systems of the master list.
Position of the 51 orbital poles in Galactic coordinates for systems closer than 18 pc.
Appendix A: Notes on individual systems and new orbits
WDS 00057+4549AB. The Sixth Catalog of Orbits of Visual Binary Stars lists two orbits for the AB pair of that system (separated by about 6″): one with a period of 1550 yr (Popovic & Pavlovic 1996) and another with a period of 509 yr (Kiyaeva et al. 2001), indicating that the solution is not yet well constrained. Radial velocities obtained by CORAVEL and by Tokovinin & Smekhov (2002), as listed in Table A.4 and displayed in Fig. A.1, correspond to a radial velocity difference V_{B} − V_{A} slightly decreasing from −2.7 km s^{1} in 1980 to −3.8 km s^{1} in 1996. The value −3.0 km s^{1} has been used by Kiyaeva et al. (2001) to determine unambiguously the ascending node. This result is confirmed thanks to the precise radialvelocity measurements of Marcy (2013, priv. comm.) that show a linear variation for V_{B} − V_{A} amounting to −6.4 m s^{1} yr^{1} from 2000 to 2010, while a value of −11.5 m s^{1} yr^{1} is calculated from the visual orbit. Thus we have computed the position of the orbital pole but it is not expected to be precise.
Kiyaeva et al. (2001) also suggested that ADS 48 ABF is a hierarchical triple system whose orbits (AB) and (ABF) are not coplanar. This hypothesis was contradicted by Cvetković et al. (2012), who concluded that the F component has a common proper motion with the AB pair, but that it is not bound by gravity.
WDS 00184+4401AB. Indeterminate visual orbit.
WDS 00321+6715Aa,Ab and AB. The preliminary visual orbits of this triple hierarchical system were computed by Docobo et al. (2008). No radial velocity measurements are available to resolve the ambiguity of the ascending node.
WDS 01083+5455Aa,Ab. Drummond et al. (1995) obtained the first relative astrometric orbit for that system from their direct detection of the faint B component. Combining that orbit with earlier photocentric positions (Lippincott 1981; Russell & Gatewood 1984), masses could be derived: M_{A} = 0.74 ± 0.06 M_{⊙} and M_{B} = 0.17 ± 0.01 M_{⊙}.
Duquennoy & Mayor (1991) derived a preliminary spectroscopic orbit, and Abt & Willmarth (2006) obtained thereafter more radial velocities; however, Abt & Willmarth (2006) did not attempt to compute a combined orbit, and merged the orbital elements ofDuquennoy & Mayor (1991) and Drummond et al. (1995). Because this is not fully satisfactory and is prone to confusion (as further discussed below), especially regarding the ambiguity on ω, we decided to compute a new combined orbit (Fig. A.2), using data of Duquennoy & Mayor (1991) as well as more recent unpublished CORAVEL data (Table A.5), plus velocities reported by Abt & Willmarth (2006). The astrometric positions of B relative to A taken from the WDS database at USNO were used as well. As shown in Table A.2, the value of ω for the astrometric orbit of B around A is 329.4° ± 3.8° (thus a change of 180° from the relative orbit of Drummond et al. 1995; their Table 9), and we stress that the value (332.7° ± 3.1°) listed by Abt & Willmarth (2006) for the spectroscopic orbit (of A around the centre of mass of the system) is incorrect, as it should differ by 180° from the astrometric value, confirming the suspected confusion. It is in fact just the value reported by Drummond et al. (1995), which should not have been copied without change! Instead, the value (147.9° ± 4.9°) of Duquennoy & Mayor (1991) for the spectroscopic orbit is correct.
Fig. A.1 Top panel: radial velocities for the two components of the pair HIP 473 AB = WDS 00057+4549AB. Dots refer to CORAVEL velocities and triangles to velocities from Tokovinin & Smekhov (2002) or Chubak et al. (2012) for the last point, with the time spanned by the measurements represented by the horizontal bar. Open symbols refer to component A and filled symbols to component B. The parallel dashed lines are just a guide to the eye. Bottom panel: radialvelocity difference BA, with their errors taken as the rootmeansquare of the errors on A and B. 
Fig. A.2 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 01083+5455Aa,Ab = HIP 5336 = HD 6582 = μ Cas. The solid line in the lower panel corresponds to the predicted velocity curve for A. Triangles are CORAVEL measurements, squares are measurements from Abt & Willmarth (2006). Individual measurements for the visual orbit were taken from the WDS database. 
WDS 013985612AB. Indeterminate visual orbit. In the absence of recent radialvelocity measurements it is not possible to lift the ambiguity on the position of the ascending node.
WDS 01418+4237AB. Spectroscopic orbits were published by Duquennoy & Mayor (1991) and Abt & Willmarth (2006). The present orbit (Table A.2) improves upon those, as it includes several unpublished CORAVEL velocities given in Table A.6. Lippincott et al. (1983) and Söderhjelm (1999) published astrometric orbits. Here, we obtained a combined astrometric/spectroscopic solution by considering the astrometric measurements from the Fourth Catalogue of Interferometric Measurements of Binary Stars^{12}. This solution agree well with the existing solutions quoted above.
WDS 01425+2016AB. This system was announced as an astrometric binary with a period of 0.567 yr by the Double and Multiple Star Annex (DMSA) of the Hipparcos Catalogue (ESA 1997). However, reprocessing the Hipparcos data along the method outlined by Pourbaix (2000) for a combined solution including the CORAVEL velocities does not allow us to confirm the Hipparcos DMSA/O solution. This system was therefore dropped from our final list.
WDS 02278+0426AB. The most recent and thorough analysis of that system was presented by Andrade (2007), including a revised visual orbit. CORAVEL data (Table A.7) make it possible now to compute a combined SB2/astrometric orbit, which thus gives access for the first time to the masses of the system components. The combined orbital elements and other physical parameters of the system (such as orbital parallax and component masses) are listed in Table A.3, which identifies the components as 1 and 2 rather than A and B. It must be stressed that there is no easy way to identify for sure the spectroscopic component 1 with either visual component A or B, since both components fall in the spectrograph slit and are of almost equal brightness: Hipparcos gives m_{Hp,A} = 9.45 and m_{Hp,B} = 9.63. An indirect way to identify components A and B is by noting that in Table A.3, the velocity semiamplitude K_{2} is higher than K_{1}, hence M_{2} is smaller than M_{1}. Since both components lie on the main sequence, component 2 should thus be less luminous than component 1, hence component 2 should be identified with B. It is possible to check quantitatively this assumption by using the massluminosity relationship provided by Kroupa et al. (1993), namely M_{V,1} = 5.4 for M_{1} = 0.95 M_{⊙} and M_{V,2} = 6.2 for M_{2} = 0.85 M_{⊙}, or M_{V,2} − M_{V,1} = 0.8, somewhat larger than the observed m_{Hp,B} − m_{Hp,A} = 0.18, but certainly consistent with the error bars on M_{1,2}. Although Andrade (2007) classified HD 15285 A as K7V based on several spectral features, among which the presence of weak TiO bands, its inferred mass of 0.95 M_{⊙} is more typical of a G5V type, while the mass of component B would flag it as K0V (Cox 2000).
These masses were derived from the orbital parallax of 52 ± 2 mas, to be compared with the Hipparcos value of 60.2 ± 1.7 mas (ESA 1997) or 58.3 ± 1.1 mas (van Leeuwen 2007), which yield total masses of 1.16 and 1.28 M_{⊙}, respectively, to be compared with 1.8 M_{⊙} for the orbital parallax. From the fractional mass M_{2}/ (M_{1} + M_{2}) = 0.47 (Table A.3), one obtains, (M_{1},M_{2}) = (0.62,0.54) and (0.68, 0.60) M_{⊙}. These values agree much better with a K7V spectral type for component A.
HERMES and CORAVEL radial velocities for components of some visual binaries lacking velocities so far.
WDS 02361+0653A. The 22 CORAVEL measurements available for component A (spanning 6837 d) yield a standard deviation of 0.29 km s^{1}, identical to the average error on a single measurement. The constancy of the radial velocity was confirmed by Chubak et al. (2012), who obtained 62 measurements spanning 1299 d, with a standard deviation of 0.142 km s^{1}. Hence, radial velocities cannot be used to lift the ambiguity on the orbital orientation.
WDS 02442+4914AB. Indeterminate visual orbit and roughly constant radial velocity between 1900 and 1990.
WDS 035750110AB. Indeterminate visual orbit and almost constant radial velocity between 1984 and 1993.
WDS 041530739BC. In the triple system STF 518 ABC, B is the white dwarf omi Eri B and C the flare star DY Eri.
WDS 04312+5858Aa,Ab. Triple system, with the A component being an astrometric binary in the visual binary STI 2051 AB.
WDS 05074+1839AB. Visual orbit highly conjectural and radial velocity from Duquennoy & Mayor (1991) is constant.
WDS 060033102AC. A quadruple system (Tokovinin et al. 2005): HU 1399 AB (P = 67.7 yr, a = 0.912′′), HJ 3823 AC (P = 390.6 yr, a = 3.95′′ by Baize 1980, not confirmed by Tokovinin et al. 2005), TOK 9 C,CE (P = 23.7 yr, a = 0.12′′). A, B and C, CE are coplanar (Tokovinin et al. 2005).
WDS 06262+1845AB. Triple system, with A known as SB2 with P = 6.99 d (Griffin & Emerson 1975).
WDS 062930248AB. Solution for the SB2 orbit includes masses and orbital parallax (Ségransan et al. 2000). Although not clearly stated, the orbital elements listed by Ségransan et al. (2000) correspond to the spectroscopic orbit (hence ω and Ω listed in our Table 3 differ by 180° from theirs).
WDS 065794417AB. Indeterminate visual orbit.
WDS 07100+3832. The spectrovisual orbit of this system was only recently computed by Barry et al. (2012).
WDS 071754659AB. The visual orbit is still very uncertain. The total mass of the system is abnormally low given the spectral type.
WDS 07346+3153AB. A difficult system since both A and B are spectroscopic binaries (Vinter Hansen et al. 1940), with P = 9.21 d and P = 2.93 d, respectively. Batten (1967) has lifted the ambiguity on the ascending node for that system, and we used the orbital elements of the 6th Catalogue of Visual Binary Orbits to compute the position of the pole, agreeing with that of Batten (1967), allowing for the fact that our poles follow the righthand rule, opposite to that of Batten.
WDS07393+0514AB. The pole ambiguity was lifted from the radialvelocity drift present in the Duquennoy et al. (1991) data.
WDS 08592+4803A,BC. Indeterminate visual orbit. It might be a quadruple system, since Aa seems to be a SB1 system with a period P ~ 11 yr (Abt 1965) although the spectroscopic orbit is not of an excellent quality.
WDS 09006+4147AB. The period of 21.05 yr for that visual binary is well established (Muterspaugh et al. 2010a). Based on an older (albeit similar) value of the visual period from Heintz (1967), Abt & Levy (1976) proposed a spectroscopic orbit based on Heintz’ visual elements plus approximate values for the systemic velocity and the semiamplitude, which could account for Abt & Levy’s own velocity measurements and older ones from Underhill (1963). Although it is rather uncertain, this orbit is still the one quoted by the SB9 catalogue (Pourbaix et al. 2004). The spectroscopic orbit is difficult to derive because (i) the components are not separated widely enough on the sky to record their spectra separately; (ii) the velocity semiamplitude is small; and (iii) component A appears to rotate relatively fast. Therefore, the available unpublished CORAVEL measurements appear useless, since the two components appear hopelessly blended in the crosscorrelation function. More recent HERMES/Mercator spectra (Raskin et al. 2011) have a resolution high enough to allow a comfortable doubleGaussian fit of the crosscorrelation function (CCF; Fig. A.5). Although such a fit has not necessarily a unique solution, we found it reassuring that the CCFs for all seven available dates (Table A.1) are well fitted with Gaussians having the same properties: height ~0.08 and σ = 14 km s^{1} for the F4V component A (V = 4.2), and height ~0.04 and σ = 4 km s^{1} for the K0V component B (V = 6.5). It thus appears that component A is a relatively fast rotator. The assignment of either of the two CCF peaks to component A thus relies on the assumption that component A, being of spectral type F4V, should correspond to the deepest peak when computed with a F0V template.
Fig. A.3 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 02278+0426AB. In the lower panel, filled and open triangles correspond to CORAVEL velocities for components 1 and 2, respectively. Individual measurements for the visual orbit were taken from the WDS database at USNO. 
Fig. A.4 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 09006+4147AB. Filled triangles in the lower panel refer to component A. Individual measurements for the visual orbit were taken from the WDS database. 
The validity of our doubleGaussian fits, and of the resulting orbital elements listed in Table A.3, may be assessed a posteriori by the comparison of the orbital parallax (53.8 ± 1.0 mas) with the Hipparcos parallax (60.9 ± 1.3 mas; ESA 1997). The discrepancy probably comes from the fact that the limited phase coverage for the radial velocities does not yet impose strong enough constraints on the velocity semiamplitudes (Fig. A.4). The astrometrically derived masses (M_{A} = 1.37 M_{⊙}, M_{B} = 1.04 M_{⊙}, in relatively good agreement with the spectral types; Martin & Mignard 1998; Söderhjelm 1999) are different from the types derived from our combined orbit (M_{A} = 1.73 M_{⊙}, M_{B} = 1.69 M_{⊙}), which also hints at inaccurate velocity semiamplitudes.
WDS 09144+5241AB. For this longperiod binary, the determination of the orbital pole is based on the agreement between the variation of the radial velocities V_{A} and V_{B} measured between 1910 and 1997 and the value predicted using the elements of the visual orbit. This result is confirmed by the value of the measurement of (V_{B} − V_{A}) performed with the HERMES spectrograph in 2012 (see Table A.1).
WDS 093131329AB. Massratio derived from the astrometry (Söderhjelm 1999). Two HERMES radialvelocity measurements obtained one year apart neither showed any significant drift (Vr = 7.93 km s^{1} on JD = 2 455 657.5 and Vr = 7.89 km s^{1} on JD = 2 456 038.4), nor any line doubling, as would be expected for a system consisting of two components only 0.6 arcsec apart.
WDS 09357+3549AB. Indeterminate visual orbit.
WDS 10454+3831AB. Triple system, Aa being the speckle binary CHR 191. A value of (V_{B} − V_{A}) obtained with the HERMES spectrograph in 2012 is listed in Table A.1.
WDS 11182+3132AB. Quadruple system, A and B are both SB1 with P = 670.24 d and P = 3.98 d respectively (Griffin 1998).
WDS 112476139AB. Abnormally high total mass, inconsistent with the spectral types.
WDS 13100+1732AB. Spectral types and masses of components from ten Brummelaar et al. (2000).
WDS 13198+4747AB. Triple system (Aa = CHR 193; Beuzit et al. 2004). CORAVEL sees this star as SB3, but no satisfactory attribution of the various peaks to the corresponding components could be made, making it very difficult to isolate the trend of Vr(A).
WDS 13473+1727AB. Indeterminate orbit. The component τ BooAb is a star with planet.
WDS 13491+2659AB. The pole ambiguity was lifted from unpublished CORAVEL velocities yielding the sign of Vr(B)− Vr(A) (=–3.3 km s^{1}) during the observation span 1983–1988 (Table A.1).
WDS 13547+1824. Jancart et al. (2005) obtained a combined spectroscopic/astrometric orbit based on Hipparcos data and the spectroscopic orbit of Bertiau (1957). The value of ω from Jancart et al. (2005) has been adopted in Table 3.
WDS 14514+1906AB. The visual orbital period is 151 yr (Söderhjelm 1999), so that no strong trend must be expected from the radial velocities over three decades. CORAVEL observations (Duquennoy et al. 1991) exhibit a weak upward trend between 1977 (1 km s^{1}) and 1991 (1.4 km s^{1}), which is however not significant based on the average measurement uncertainty of 0.34 km s^{1}. Later measurements by Abt & Willmarth (2006) seem to confirm that trend, but the trend is made fragile by a possible zeropoint offset between the two data sets. Nevertheless, the pole quoted in Table 3 is based on the assumption that the radialvelocity trend is real, which then implies that the node given in the 6th USNO Catalogue is the correct one.
WDS 145752125AB. Indeterminate visual orbit.
WDS 15038+4739AB. Triple star: Bb = eclipsing SB2 with P = 0.27 d (Lu et al. 2001).
WDS 15527+4227. The Hipparcos Double and Multiple Star Annex (DMSA; ESA 1997) claims HIP 77760 to be an astrometric binary with an orbital period of 0.14 yr and an inclination of 131.68°. These values seem, however, incompatible with the absence of (CORAVEL) radialvelocity variations (33 measurements with a standard deviation of 0.413 km s^{1}, to be compared with the average instrumental error of 0.380 km s^{1}).
New combined astrometricspectroscopic orbits computed from CORAVEL and/or HERMES velocities.
Fig. A.5 Crosscorrelation function of the spectrum of WDS 09006+4147AB (=HD 76943) obtained with the HERMES/Mercator spectrograph (Raskin et al. 2011) on HJD 2 456 244.717 (November 13, 2012), correlated with a F0V template. Both components (separated by 0.5′′ on the sky) enter the fiber (2.5′′ on the sky). Component A corresponds to the wide peak centred at 31.5 km s^{1}. 
Fig. A.6 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 16413+3136AB. The solid line in the lower panel corresponds to the predicted velocity curve for A. Individual measurements for the visual orbit were taken from the WDS database at USNO. 
Fig. A.7 Combined astrometric (upper panel) and spectroscopic (lower panel) orbit for WDS 17121+4540AB. The solid line in the lower panel corresponds to the predicted velocity curve for A. Individual measurements for the visual orbit were taken from the WDS database at USNO. 
WDS 16413+3136AB. There are many visual orbits available for that system, the most recent ones by Baize (1976) and Söderhjelm (1999). Scarfe et al. (1983) already computed an orbit combining visual data and a collection of Coudé radial velocities. Unpublished CORAVEL velocities given in Table A.8 nicely complement Scarfe’s velocities, yielding a somewhat more accurate orbit as listed in Table A.2. A third component has been detected by IR speckle interferometry (McCarthy 1983).
New combined astrometricspectroscopic orbits (2 observable spectra) computed from CORAVEL or HERMES velocities.
Fig. A.8 Combined astrometric (upper panel) and spectroscopic (lower panel) orbit for WDS 18070+3034AB. The solid line in the lower panel corresponds to the predicted velocity curve for A. Individual measurements for the visual orbit were taken from the WDS database at USNO. In the lower panel, triangles correspond to CORAVEL data, squares to velocities from Abt & Willmarth (2006), diamonds to Tokovinin & Smekhov (2002), and double triangles to the old radial velocities collected by Kamper & Beardsley (1986). 
CORAVEL velocities (labelled C in the last column), Tokovinin & Smekhov (2002) velocities (labelled T), and Chubak et al. (2012) velocities (labelled M) for the AB components of WDS 00057+4549AB = HIP 473.
CORAVEL velocities for the components of WDS 01083+5455Aa,Ab = μ Cas.
New CORAVEL velocities for the components of WDS 01418+4237AB in complement of velocities published by Duquennoy & Mayor (1991).
CORAVEL velocities for the components of WDS 02278+0426AB.
New CORAVEL velocities for the components of WDS 16413+3136AB in complement of velocities published by Duquennoy & Mayor (1991).
WDS 165550820AB. A triple system composed of M dwarfs, in which all three stars are visible in the spectra. Spectroscopic orbits are available for the inner (Bab) and outer pairs (ABab), both as SB2 (Ségransan et al. 2000). The two orbits are probably coplanar (Mazeh et al. 2001). The whole system is probably septuple.
WDS 17121+4540AB. A visual orbit for that system was obtained by Hartkopf et al. (1996). The available CORAVEL velocities shown in Table A.9 make it possible to derive a combined orbit for the first time. Visual observations are from the CHARA database. The new combined orbit is presented in Fig. A.7.
New CORAVEL velocities for the components of WDS 17121+4540AB.
WDS 171532636AB. Irwin et al. (1996) derived precise BA radial velocities (except for an acceleration component that they were unable to explain). Velocity measurements of Vr(A) and Vr(B) were obtained in 2012 by the HERMES spectrograph (Table A.1).
WDS 171914638AB. Indeterminate visual orbit.
WDS 173040104AB. A combined astrometric/spectroscopic solution has been computed by Pourbaix (2000) and has been used to compute the position of the orbital pole.
WDS 17349+1234. Spectral types from Gatewood (2005). Kamper et al. (1989) performed a combined spectroscopic/astrometric analysis of this system, and their orbit has been recently updated by Hinkley et al. (2011), from whom we adopted the orbital elements without requiring any change to comply with our conventions.
WDS 17364+6820Aa,Ab. Velocity measurements of Vr(A) and Vr(B) obtained in 2012 by the HERMES spectrograph are listed in Table A.1.
WDS 17465+2743BC. HD 161797 (μ Her) is a quadruple system (Raghavan et al. 2010). The components of the wide pair STF 2220 AB (35″) have common proper motion. The component A is an astrometric binary (period of 65 years, Heintz 1994). The components B and C form the visual binary AC 7 for which the orbital pole ambiguity has been lifted from unpublished CORAVEL velocities, yielding a drift of Vr(B) over 1975–1995.
WDS 18070+3034AB. Based on available CORAVEL velocities (unpublished measurements are given in Table A.10), a combined astrometric/spectroscopic solution (Table A.2 and Fig. A.8) improves the spectroscopic orbital elements provided by Abt & Willmarth (2006) for that system. AC 15 AB is a triple system, with the separation AaAb amounting to 0.228″ and that of AB to 0.851″ in 2005, according to Scardia et al. (2008).
WDS 18428+5938AB. Velocity measurements (4 for Vr(A) and one for Vr(B)) obtained by the HERMES spectrograph in 2012 are listed in Table A.1.
WDS 18570+3254AB. A combined astrometric/spectroscopic solution confirms the spectroscopic orbital elements provided by Abt & Willmarth (2006) for that system. BU 648 AB is a binary with an exoplanet (Muterspaugh et al. 2010c) in coplanar orbits.
WDS 19121+0254AB. Masses are from Benedict et al. (2001).
Fig. A.9 Crosscorrelation function of the HERMES spectrum of WDS 19255+0307Aa,Ab (obtained at JD 2 455 657.742) with a F0V template. Component Aa shows a rotationally broadened profile with V_{r}sini = 84 km s^{1}, and component Ab is visible as a weak Gaussian (at V_{r}(Ab) = –37.70 km s^{1}) superimposed on the rotationally broadened profile (centred on V_{r}(Aa) = –26.68 km s^{1}). Similar profiles have been observed over several consecutive nights. 
WDS 19255+0307Aa,Ab. Kamper et al. (1989) performed a combined spectroscopic/astrometric analysis of the AaAb pair. We adopted the orbital elements from their Table V, except that 180° was added to their value of ω, to conform with our convention that ω corresponds to the visual orbit of Ab around Aa, rather than to the spectroscopic orbit of Aa around the centre of mass. One HERMES measurement was obtained (Table A.1) that reveals for the first time the companion in the crosscorrelation function as a narrow peak on top of the broad peak (V_{r}sini = 84 km s^{1}) due to the FV primary (Fig. A.9). The companion does not rotate fast. From the orbit described in Kamper et al. (1989), we estimate that the HERMES data point (V_{r}(Aa) = –26.68 km s^{1}) has been taken at orbital phase 0.54 ± 0.14, when the velocity is km s^{1} higher than the centre of mass (CoM) velocity. The CoM velocity is thus estimated to be km s^{1}. Combined with V_{r}(Ab) = –37.70 km s^{1}, this yields a surprisingly high mass ratio M_{Aa}/M_{Ab} of 5.6, or a mass of 0.29 M_{⊙} for the companion if the primary is a F0V star of mass 1.6 M_{⊙}. The corresponding spectral type would then be M4V for the Ab component, 8.6 mag fainter than the primary component. The close pair WDS 19255+0307Aa,Ab, seems to form a triple system with the wide pair (WDS 19255+0307AB = BUP 190, separation 96 to 133″), but the physical nature of this system is yet to be confirmed.
WDS 21000+4004AB. Fekel et al. (1978) found that the visual component A of KUI 103 is itself a spectroscopic binary with a period of 3.27 d, but no information is available about coplanarity (Tokovinin 2008). Pourbaix (2000) has made a combined analysis of the visual and spectroscopic data. However, the authors note that this orbital solution is somewhat uncertain, given the large error bars and inconsistencies in the derived values for stellar masses and orbital parallax. Thus, we decided to base our determination of the orbital pole of this system on the analysis of the published radial velocity measurements (Pourbaix 2000) and the orbital elements of the new visual orbit proposed by Docobo & Ling (2010).
WDS 21069+3845AB. The mass ratio is taken from Gorshanov et al. (2006). The node assignment is based on the difference V_{r}(B) – V_{r}(A) ~ + 1 km s^{1}> 0, as provided by Table 1 of Gorshanov et al. (2006) and confirmed by unpublished CORAVEL measurements.
WDS 22234+3228AB. The total mass seems too high for the spectral type.
WDS 22280+5742AB. Masses are from Delfosse et al. (2000). Among the many radial velocities available in the literature for the components of this system, we give the highest weight to the measurement [V_{r}(A) = –32.7 km s^{1} and V_{r}(B) = –33.2 km s^{1}] by Gizis et al. (2002) because (i) that measurement is precisely dated (JD 2 450 005.6) and on the same night for A and B; and (ii) the 0.5 km s^{1} velocity difference is consistent with the prediction from the visual orbit. The measurement (V_{r}(A) = –34.0 km s^{1} and V_{r}(B) = –35.0 km s^{1}) by Delfosse et al. (1998) is equally useful (its epoch is not as precisely known as that report by Gizis, but it was obtained between September 1995 and March 1997, a short time span with respect to the 44.7 yr orbital period); the 1 km s^{1} velocity difference agrees better with the prediction from the visual orbit than Gizis’ 0.5 km s^{1}. The determination of the sign of V_{r}(B) – V_{r}(A) = –1 km s^{1} is important and it is opposite to the prediction from the visual orbital elements from the 6th COVBS. It is thus necessary to add 180° to the values of Ω and ω listed in the 6th COVBS. It is noteworthy that the velocity values listed in the Wilson – Evans – Batten catalogue (Duflot et al. 1995, V_{r}(A) = –24.0 km s^{1} and V_{r}(B) = –28.0 km s^{1}) or in Reid et al. (1995) (V_{r}(A) = –16.0 km s^{1} and V_{r}(B) = –29.0 km s^{1}), although the differences are much larger than anticipated, always have V_{r}(B) – V_{r}(A) < 0.
New CORAVEL velocities for the components of WDS 18070+3034AB in complement of velocities published by Duquennoy & Mayor (1991).
WDS 23317+1956AB. Indeterminate visual orbit.
WDS 23524+7533AB. Indeterminate visual orbit. A triple system with A being SB2 with a period of 7.75 d (Christie 1934).
Appendix B: Determining the orbital pole of a double star
Appendix B.1: Orbital pole of a visual double star
The true relative orbit (B/A) is related to the classical visual orbit of B relative to A. The pole of this orbit is defined by the unit vector perpendicular to the orbital plane and directed in such a way that an observer placed at its end sees that the B component has a direct motion (counterclockwise).
The spatial motion of the binary components is described using the reference frame (A, x, y, z) centred on the A component with two axes in the plane tangent to the celestial sphere (which is denoted plane of the sky): Ax points north (position angle = 0°), Ay points east (position angle = 90°), and the third axis Az is along the line of sight, pointing in the direction of increasing radial velocities (i.e. positive radial velocity). This reference frame is thus retrograde (i.e., viewed from the positive side of the Azaxis, a rotation of the Axaxis onto the Ayaxis is carried out clockwise).
Fig. B.1 True and the projected relative orbit (B/A) of a visual binary and its geometrical elements. The ascending node of the orbit is located at point N. 
The relative (B/A) orbit is described by means of seven orbital elements:
(i) Four socalled dynamic elements specifying the properties of the Keplerian motion in the true orbit (T, the epoch of passage through periastron; P, the revolution period; a, the semimajor axis and e, the eccentricity).
(ii) Three socalled geometric elements that define the orientation of the true and apparent orbits, the latter being the projection of the true orbit on the plane of the sky (see Fig. B.1):

Ω is the position angle of the line of intersection between the true orbital plane and the plane of the sky. There are two nodes whose position angles differ by 180°, but by convention, Ω is the position angle of the ascending node where the orbital motion is directed away from the Sun. If radialvelocity measurements are not available to identify the ascending node, a temporary value ranging between 0° and 180° is adopted (Finsen 1934).

i, the inclination, is the angle between the planes of the projected orbit and of the true orbit, taken at the ascending node. The motion is direct if 0°<i< 90° and retrograde if 90°<i< 180°.

ω is the argument of the periastron in the true orbit plane, measured in the direction of the orbital motion from the ascending node with a value ranging from 0° to 360°.
Two other orbits may be considered: the (A/AB) orbit of A around the centre of mass AB and the (B/AB) orbit of B around the centre of mass.
For resolved binaries, astrometric observations give then access to the (A/AB) and (B/AB) orbits projected on the plane of the sky. For an unresolved binary, the astrometric measurements capture the position of the photocentre F located between A and AB. The astrometric orbit (F/AB) is similar to that of (A/AB), but with a smaller semimajor axis a_{F}.
For a given projected orbit as derived from astrometric observations, two possible true orbits correspond that are symmetric with respect to the plane of the sky.
Therefore, even though the dynamic elements and the inclination i are determined unambiguously, it is impossible to select between the two possible true orbits (and thus between the two possible ascending nodes Ω) without a radialvelocity measurement. The argument of periastron ω being measured from the ascending node, the ambiguity on Ω thus propagates onto ω. One special case is i = 0° or 180°, the true orbit being then in the plane of the sky, and there is no ambiguity on the ascending node (since there is no such thing as an ascending node in that case!). In all other cases, radialvelocity measurements for at least one component are needed to lift the ambiguity on the value of Ω, before the orientation of the orbital pole can be fixed , as was done by Dommanget (2005) and Dommanget & Nys (2006).
Appendix B.2: Radial velocities
If V_{r}(A), V_{r}(B) and V_{r,sys} are the heliocentric radial velocities of components A and B and of the centre of mass of the binary AB, and denoting V_{A / AB} and V_{B / AB} the radial components of the orbital velocity of A and B around the centre of mass, then (B.1)The orbital radial velocity of each component relates to the orbital elements by (B.2)where v is the true anomaly, a_{A,B} are the semimajor axes, P the orbital period, and ω_{(A,B)} the argument of the periastron of the (A/AB) and (B/AB) orbits.
For doublelined spectroscopic binaries (SB2), the measurements of V_{A / AB} and V_{B / AB} are possible; for singlelined spectroscopic binaries (SB1), only V_{A / AB} measurements are possible.
The relative radial velocity between the components, V_{r} = V_{r}(B) − V_{r}(A), can be calculated from the orbital elements of the (B/A) orbit using Eq. (B.2).
For an astrometricbinary orbit, the radial velocity of the photocenter V_{r}(F) can be computed with the same formula (Eq. (B.2)) where the semimajor axis a_{F} and the argument of the periastron ω_{F} of the photocentric orbit have been inserted instead.
It is easy to show that (B.3)The variation of V_{r} and V_{r}(B) are the same and are opposite to that of V_{r}(A). By definition, the maximum of the V_{r}(A) and V_{r}(B) curves occurs at the passage of the ascending node in the true orbits (A/AB) and (B/AB), respectively, whereas the maximum of the V_{r} curve corresponds to the passage at the ascending node in the relative orbit (B/A).
Appendix B.3: Choosing the correct ascending node
Generally, the ascending node Ω of the relative orbit (B/A) is determined by comparing the variations of the measured radial velocities V_{r} with the shape of the relative radialvelocity curve (or (F)) computed from the visual orbital elements. However, finding the ascending node of the relative orbit (B/A) in practice depends on the type of binary system and on the type of radialvelocity measurements available, as described in Table B.1.
The different ways to determine the ascending node of the relative orbit (B/A).
Fig. B.2 Examples of lifting the ambiguity on the ascending node for a visual binary. Top (cases 2, 5 of Table B.1): the slope of the computed curve is opposite to that of the V_{r}(B) − V_{r}(A) measurements: it is then necessary to add 180° to the values of Ω and ω given by the 6th Catalog of Orbits of Visual Binary Stars (6th COVBS; Hartkopf et al. 2001). Bottom (cases 3, 6 of Table B.1): the slope of the computed curve is opposite to that of V_{r}(A) measurements; there is no need to change the values of Ω and ω listed in the 6th COVBS. 
For spectroscopic systems with a combined spectrovisual orbit available (case 1 of Table B.1), the ascending node of the relative orbit is unambiguously determined (e.g., Pourbaix 2000). The values of the orbital elements i, Ω and ω can then be directly used to compute the orbital pole orientation.
For visual binaries with available radialvelocity measurements, if the slope of the computed curve is the same as that of the V_{r}(B) – V_{r}(A) measurements (cases 2 or 5 of Table B.1), or opposite to that of the V_{r}(A) curve (cases 3 or 6 of Table B.1), there is no need to change the values of the position angle of the ascending node Ω and of the argument of the periastron ω, from the values adopted for computing in Eq. (B.2). In the opposite situation, it is necessary to add 180° to the values of Ω and ω. Examples are given in Fig. B.2.
For astrometric binaries (cases 4 or 7 of Table B.1), if the drift of the computed (F) curve is opposite that of the V_{r}(A) measurements, there will be no need to change the values of the position angle of the ascending node Ω_{F} or of the argument of the periastron ω_{F}. In the opposite case, it is necessary to add 180° to these values. Thus the ascending node of the (B/A) orbit will be Ω = Ω_{F} + 180° and ω = ω_{F} + 180°.
It is important to note that the ambiguityfree values listed in Table 3 are correspond to the relative orbit (B/A). Hence, in Tables A.2 and A.3, which provide newly determined orbital elements of spectrovisual orbits, we list the values of ω_{B}, contrary to usage with spectroscopic orbital elements, to ensure consistency with Table 3.
Appendix B.4: Determining the direction of the pole of the true relative orbit
This section describes the determination of the direction of the pole of the true relative orbit (B/A), including the case of astrometric binaries for which the B component is invisible.
In the reference frame (A, x′, y′, z′) attached to the (B/A) orbit, the direction of the orbital pole is given by the unit vector (x′, y′, z′) = (0, 0, 1).
The ambiguity on the position angle of the ascending node was solved along the method described in Sect. B.3, the geometric orbital elements i, Ω and ω are known, so that the equatorial coordinates of the pole (x_{eq}, y_{eq}, z_{eq}) can be computed from the equatorial coordinates of the binary (Right Ascension α, Declination δ) by a succession of referenceframe changes as shown in Fig. B.3:

From (A, x′, y′, z′) to (A, x′, u, z) by a rotation of angle i around the Ax′ axis;

to (A, x, y, z) by a rotation of angle Ω around the Az axis;

to (O, v, y, z_{eq}) by a rotation of angle δ around the Ay axis;

and finally to (O, x_{eq}, y_{eq}, z_{eq}) by a rotation of angle α around the Az_{eq} axis.
It is then easy to show that in the equatorial reference frame (O, x_{eq}, y_{eq}, z_{eq}), the direction of the orbital pole is given by the unit vector with the coordinates (B.4)We stress that our convention, which leads to a polar direction opposite to that of the system when i = 0°, is opposite to that of Batten (1967). In the terminology of Batten (1967), our poles are thus those that would be derived by applying the righthand rule to the orbital motion.
Fig. B.3 Orientation of the true and apparent relative orbits (B/A) with respect to the celestial sphere and the corresponding reference frames. (α,δ) are the equatorial coordinates of the binary system. 
The Galactic coordinate system is the most appropriate for studying the distribution of the orbital poles in a Galactic framework. Therefore, the coordinates (x_{eq}, y_{eq}, z_{eq}) of the unit vector giving the direction of the pole are converted into the coordinates (x_{g}, y_{g}, z_{g}) in the Galactic frame using the relation (for details see Johnson & Soderblom 1987) (B.5)The terms of the transformation matrix [ T ] are computed from the J2000 equatorial coordinates of the Galactic centre and pole given by the Institut de Mécanique Céleste et de Calcul des Ephémérides (IMCCE):

North Galactic pole: α_{0} = 12^{h}51^{m}26.28^{s} and δ_{0} = 27°7′41.7′′

Galactic center: α_{G} = 17^{h}45^{m}37.20^{s} and δ_{G} = −28°56′10.2′′
This leads to the transformation matrix: (B.6)After the coordinates (x_{g}, y_{g}, z_{g}) are computed, the Galactic longitude l and latitude b are deduced using the relations (B.7)
See also Kirkpatrick et al. (2012), who recently compiled an exhaustive list of stars nearer than 8 pc from the Sun, and Raghavan et al. (2010), who collected orbits for visual binaries with Hipparcos parallaxes larger than 40 mas and with primary spectral types F6K3 (their Table 11).
As apparent from Fig. 3, restricting the sample to 18 systems would yield an even more significant result, but this would imply considering only one orbital pair from the triple system WDS 11182+3132 (Table 3), which is not physically sound. A discussion of the problem of triple systems is presented in Sect. 6.
Sterzik & Tokovinin (2002) revealed that there are only 22 triple systems with both orbits determined visually. Eight of these have questionable elements and only three have had their ascending nodes identified using radial velocity measurements.
Available from http://www.usno.navy.mil/USNO/astrometry/opticalIRprod/wds/int4
Acknowledgments
We dedicate this paper to the memory of J. Dommanget, who initiated this study and could unfortunately not see it accomplished. He would certainly have preferred a more significant result, but he had faith in statistics.
We thank the anonymous referee for his constructive remarks and attitude. We thank A. Lobel, P. Neyskens and C. Siopis, who acted as observers on the HERMES spectrograph and secured the necessary spectra to solve the pole ambiguity for several binary systems. Based on observations made with the Mercator Telescope, operated on the island of La Palma by the Flemish Community, at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. Based on observations obtained with the HERMES spectrograph, which is supported by the Fund for Scientific Research of Flanders (FWO), Belgium, the Research Council of K.U. Leuven, Belgium, the Fonds National de la Recherche Scientifique (FNRS), Belgium, the Royal Observatory of Belgium, the Observatoire de Genève, Switzerland and the Thüringer Landessternwarte Tautenburg, Germany. We gratefully thank D. Pourbaix (Institut d’Astronomie et d’Astrophysique de l’Université Libre de Bruxelles) for his help in the calculation of the new orbits, as well as G. Marcy for sharing his velocity measurements of WDS 00057+4549AB. This research has made use of the Simbad database, operated at the Centre de Données Astronomiques de Strasbourg (CDS), France. Available at http://cdsweb.ustrasbg.fr/. This research has made use of the Washington Double Star Catalog maintained at the US Naval Observatory. Available at http://ad.usno.navy.mil/wds/.
References
 Abt, H. A. 1965, ApJS, 11, 429 [NASA ADS] [CrossRef] [Google Scholar]
 Abt, H. A., & Levy, S. G. 1976, ApJS, 30, 273 [NASA ADS] [CrossRef] [Google Scholar]
 Abt, H. A., & Willmarth, D. 2006, ApJS, 162, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Aitken, R. G. 1918, Lick Observatory Bull., 9, 184 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, T. W., & Stephens, M. A. 1972, Biometrika, 59, 613 [CrossRef] [Google Scholar]
 Andrade, M. 2007, Rev. Mex. Astron. Astrofis., 43, 237 [NASA ADS] [Google Scholar]
 Baize, P. 1976, A&AS, 26, 177 [NASA ADS] [Google Scholar]
 Baize, P. 1980, Inf. Circ. IAU Comm., 26, 80 [Google Scholar]
 Balega, I. I., Balega, Y. Y., Hofmann, K.H., et al. 2005, A&A, 433, 591 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baranne, A., Mayor, M., & Poncet, J. L. 1979, Vistas in Astron., 23, 279 [NASA ADS] [CrossRef] [Google Scholar]
 Barry, R. K., Demory, B.O., Ségransan, D., et al. 2012, ApJ, 760, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Bate, M. R. 2000, MNRAS, 314, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Batten, A. H. 1967, in Determination of Radial Velocities and their Applications, eds. A. H. Batten, & J. F. Heard, IAU Symp., 30, 199 [Google Scholar]
 Benedict, G. F., McArthur, B. E., Franz, O. G., et al. 2001, AJ, 121, 1607 [NASA ADS] [CrossRef] [Google Scholar]
 Beran, R. J. 1968, J. Appl. Probab., 5, 177 [CrossRef] [Google Scholar]
 Bertiau, F. C. 1957, ApJ, 125, 696 [NASA ADS] [CrossRef] [Google Scholar]
 Beuzit, J.L., Ségransan, D., Forveille, T., et al. 2004, A&A, 425, 997 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bingham, C. 1974, Ann. Stat., 2, 1201 [CrossRef] [Google Scholar]
 Bodenheimer, P. 1978, ApJ, 224, 488 [NASA ADS] [CrossRef] [Google Scholar]
 Bodenheimer, P., & Burkert, A. 2001, in The Formation of Binary Stars, eds. H. Zinnecker, & R. Mathieu, IAU Symp., 200, 13 [Google Scholar]
 Bonnell, I. A. 2001, in The Formation of Binary Stars, eds. H. Zinnecker, & R. Mathieu, IAU Symp., 200, 23 [Google Scholar]
 Boss, A. P. 1988, Comm. Astrophys., 12, 169 [Google Scholar]
 Briggs, M. S. 1993, ApJ, 407, 126 [NASA ADS] [CrossRef] [Google Scholar]
 Chang, K. 1972, AJ, 77, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Christie, W. H. 1934, ApJ, 80, 181 [NASA ADS] [CrossRef] [Google Scholar]
 Chubak, C., Marcy, G., Fischer, D. A., et al. 2012 [arXiv:1207.6212] [Google Scholar]
 Couteau, P. 1960, Bull. Astron. Paris, 23, 127 [NASA ADS] [Google Scholar]
 Cox, A. N. 2000, Allen’s Astrophysical Quantities (New York: Springer) [Google Scholar]
 Cvetković, Z., Pavlović, R., Ninković, S., & Stojanović, M. 2012, AJ, 144, 80 [NASA ADS] [CrossRef] [Google Scholar]
 Dahn, C. C., Harrington, R. S., Kallarakal, V. V., et al. 1988, AJ, 95, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Dehnen, W. 1998, AJ, 115, 2384 [NASA ADS] [CrossRef] [Google Scholar]
 Delfosse, X., Forveille, T., Perrier, C., & Mayor, M. 1998, A&A, 331, 581 [NASA ADS] [Google Scholar]
 Delfosse, X., Forveille, T., Ségransan, D., et al. 2000, A&A, 364, 217 [NASA ADS] [Google Scholar]
 Docobo, J. A., & Costa, J. M. 1985, IAU Commission on Double Stars, 96, 1 [Google Scholar]
 Docobo, J. A., & Ling, J. F. 2010, IAU Commission on Double Stars, 171, 1 [Google Scholar]
 Docobo, J. A., Tamazian, V. S., Balega, Y. Y., et al. 2008, A&A, 478, 187 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dommanget, J. 1968, Bull. Astron. Obs. Royal Belgique, 6, 246 [NASA ADS] [Google Scholar]
 Dommanget, J. 1978, Bull. Astron. Obs. Royal Belgique, 9, 47 [Google Scholar]
 Dommanget, J. 1988, Astrophys. Space Sci., 142, 171 [NASA ADS] [CrossRef] [Google Scholar]
 Dommanget, J. 2005, Ciel et Terre, 121, 2 [NASA ADS] [Google Scholar]
 Dommanget, J. 2014, Observations et Travaux, submitted [Google Scholar]
 Dommanget, J., & Nys, O. 2006, Observations et Travaux, 62, 5 [NASA ADS] [Google Scholar]
 Drummond, J. D., Christou, J. C., & Fugate, R. Q. 1995, ApJ, 450, 380 [NASA ADS] [CrossRef] [Google Scholar]
 Duflot, M., Figon, P., & Meyssonnier, N. 1995, A&AS, 114, 269 [NASA ADS] [Google Scholar]
 Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 [NASA ADS] [Google Scholar]
 Duquennoy, A., Mayor, M., & Halbwachs, J.L. 1991, A&AS, 88, 281 [NASA ADS] [CrossRef] [Google Scholar]
 Eggen, O. J. 1956, AJ, 61, 405 [NASA ADS] [CrossRef] [Google Scholar]
 Eggenberger, P., Miglio, A., Carrier, F., Fernandes, J., & Santos, N. C. 2008, A&A, 482, 631 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggleton, P. P., & KiselevaEggleton, L. 2001, ApJ, 562, 1012 [NASA ADS] [CrossRef] [Google Scholar]
 ESA 1997, The Hipparcos and Tycho Catalogues, ESA SP (Nooodwijk: ESA publications Division), 1200 [Google Scholar]
 ESA 1997, VizieR Online Data Catalog I/239 [Google Scholar]
 Fabrycky, D., & Tremaine, S. 2007, ApJ, 669, 1298 [NASA ADS] [CrossRef] [Google Scholar]
 Famaey, B., Jorissen, A., Luri, X., et al. 2005, A&A, 430, 165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Famaey, B., Pont, F., Luri, X., et al. 2007, A&A, 461, 957 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Famaey, B., Siebert, A., & Jorissen, A. 2008, A&A, 483, 453 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Farrington, C. D., ten Brummelaar, T. A., Mason, B. D., et al. 2010, AJ, 139, 2308 [NASA ADS] [CrossRef] [Google Scholar]
 Fekel, Jr., F. C. 1981, ApJ, 246, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Fekel, F. C., & Tomkin, J. 1983, PASP, 95, 1000 [NASA ADS] [CrossRef] [Google Scholar]
 Fekel, Jr., F.,Lacy, C. H., & Bopp, B. W. 1978, AJ, 83, 1445 [NASA ADS] [CrossRef] [Google Scholar]
 Finsen, W. S. 1933, Circular of the Union Observatory Johannesburg, 90, 397 [NASA ADS] [Google Scholar]
 Finsen, W. S. 1934, Circular of the Union Observatory Johannesburg, 91, 23 [NASA ADS] [Google Scholar]
 Fisher, R. 1953, Proc. Roy. Soc. A, 217, 295 [Google Scholar]
 Forveille, T., Beuzit, J.L., Delfosse, X., et al. 1999, A&A, 351, 619 [NASA ADS] [Google Scholar]
 Franz, O. G., Wasserman, L. H., Benedict, G. F., et al. 1999, in Working on the Fringe: Optical and IR Interferometry from Ground and Space, eds. S. Unwin, & R. Stachnik, ASP Conf. Ser., 194, 28 [Google Scholar]
 Gatewood, G. 2005, AJ, 130, 809 [NASA ADS] [CrossRef] [Google Scholar]
 Genova, F., Egret, D., Bienaymé, O., et al. 2000, A&AS, 143, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Geyer, D. W., Harrington, R. S., & Worley, C. E. 1988, AJ, 95, 1841 [NASA ADS] [CrossRef] [Google Scholar]
 Gillett, S. L. 1988, AJ, 96, 1967 [NASA ADS] [CrossRef] [Google Scholar]
 Giné, M. E. 1975, Ann. Stat., 3, 1243 [CrossRef] [Google Scholar]
 Girard, T. M., Wu, H., Lee, J. T., et al. 2000, AJ, 119, 2428 [NASA ADS] [CrossRef] [Google Scholar]
 Gizis, J. E., Reid, I. N., & Hawley, S. L. 2002, AJ, 123, 3356 [NASA ADS] [CrossRef] [Google Scholar]
 Glebocki, R. 2000, Acta Astron., 50, 211 [Google Scholar]
 Gorshanov, D. L., Shakht, N. A., & Kisselev, A. A. 2006, Astrophys., 49, 386 [Google Scholar]
 Griffin, R. F. 1998, The Observatory, 118, 273 [NASA ADS] [Google Scholar]
 Griffin, R. F. 2004, The Observatory, 124, 258 [NASA ADS] [Google Scholar]
 Griffin, R. F., & Emerson, B. 1975, The Observatory, 95, 23 [NASA ADS] [Google Scholar]
 Hale, A. 1994, AJ, 107, 306 [NASA ADS] [CrossRef] [Google Scholar]
 Han, I., & Gatewood, G. 2002, PASP, 114, 224 [NASA ADS] [CrossRef] [Google Scholar]
 Hartkopf, W. I., & Mason, B. D. 2003, IAU Commission on Double Stars, 151, 2 [NASA ADS] [Google Scholar]
 Hartkopf, W. I., Mason, B. D., & McAlister, H. A. 1996, AJ, 111, 370 [NASA ADS] [CrossRef] [Google Scholar]
 Hartkopf, W. I., Mason, B. D., & Worley, C. E. 2001, AJ, 122, 3472 [NASA ADS] [CrossRef] [Google Scholar]
 Hartmann, D., & Epstein, R. I. 1989, ApJ, 346, 960 [NASA ADS] [CrossRef] [Google Scholar]
 Heintz, W. D. 1967, Veröff. Sternw. München, 7, 31 [Google Scholar]
 Heintz, W. D. 1974, AJ, 79, 819 [NASA ADS] [CrossRef] [Google Scholar]
 Heintz, W. D. 1984, AJ, 89, 1063 [NASA ADS] [CrossRef] [Google Scholar]
 Heintz, W. D. 1986, A&AS, 65, 411 [NASA ADS] [Google Scholar]
 Heintz, W. D. 1987, PASP, 99, 1084 [NASA ADS] [CrossRef] [Google Scholar]
 Heintz, W. D. 1988, A&AS, 72, 543 [NASA ADS] [Google Scholar]
 Heintz, W. D. 1990a, A&AS, 82, 65 [NASA ADS] [Google Scholar]
 Heintz, W. D. 1990b, AJ, 99, 420 [NASA ADS] [CrossRef] [Google Scholar]
 Heintz, W. D. 1994, AJ, 108, 2338 [NASA ADS] [CrossRef] [Google Scholar]
 Heintz, W. D. 1996, AJ, 111, 408 [NASA ADS] [CrossRef] [Google Scholar]
 Heintz, W. D., & Cantor, B. A. 1994, PASP, 106, 363 [NASA ADS] [CrossRef] [Google Scholar]
 Heisler, J., & Tremaine, S. 1986, Icarus, 65, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Hinkley, S., Monnier, J. D., Oppenheimer, B. R., et al. 2011, ApJ, 726, 104 [NASA ADS] [CrossRef] [Google Scholar]
 Hopmann, J. 1958, Mitteilungen der UniversitaetsSternwarte Wien, 9, 177 [Google Scholar]
 Hopmann, J. 1973, Mitteilungen der UniversitaetsSternwarte Wien, 14, 18 [Google Scholar]
 Hummel, C. A., Benson, J. A., Hutter, D. J., et al. 2003, AJ, 125, 2630 [NASA ADS] [CrossRef] [Google Scholar]
 Humphreys, R. M., & Larsen, J. A. 1995, AJ, 110, 2183 [NASA ADS] [CrossRef] [Google Scholar]
 Irwin, A. W., Yang, S. L. S., & Walker, G. A. H. 1996, PASP, 108, 580 [NASA ADS] [CrossRef] [Google Scholar]
 Jancart, S., Jorissen, A., Babusiaux, C., & Pourbaix, D. 2005, A&A, 442, 365 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [NASA ADS] [CrossRef] [Google Scholar]
 Kamper, K. W., & Beardsley, W. R. 1986, AJ, 91, 419 [NASA ADS] [CrossRef] [Google Scholar]
 Kamper, K. W., Legget, D., & McCarthy, Jr., D. W. 1989, AJ, 98, 686 [NASA ADS] [CrossRef] [Google Scholar]
 Keilson, J., Petrondas, D., Sumita, U., & Wellner, J. 1983, J. Statistical Computation and Simulation, 17, 195 [CrossRef] [Google Scholar]
 Kirkpatrick, J. D., Gelino, C. R., Cushing, M. C., et al. 2012, ApJ, 753, 156 [NASA ADS] [CrossRef] [Google Scholar]
 Kiseleva, L. G., Eggleton, P. P., & Mikkola, S. 1998, MNRAS, 300, 292 [NASA ADS] [CrossRef] [Google Scholar]
 Kiyaeva, O. V., Kisselev, A. A., Polyakov, E. V., & Rafal’Skii, V. B. 2001, Sov. Astron. Lett., 27, 391 [Google Scholar]
 Konacki, M., Muterspaugh, M. W., Kulkarni, S. R., & Hełminiak, K. G. 2010, ApJ, 719, 1293 [NASA ADS] [CrossRef] [Google Scholar]
 Kozai, Y. 1962, AJ, 67, 591 [Google Scholar]
 Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 [NASA ADS] [CrossRef] [Google Scholar]
 Lippincott, S. L. 1967, AJ, 72, 1349 [NASA ADS] [CrossRef] [Google Scholar]
 Lippincott, S. L. 1972, AJ, 77, 165 [NASA ADS] [CrossRef] [Google Scholar]
 Lippincott, S. L. 1981, ApJ, 248, 1053 [NASA ADS] [CrossRef] [Google Scholar]
 Lippincott, S. L., Braun, D., & McCarthy, Jr., D. W. 1983, PASP, 95, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Lu, W., Rucinski, S. M., & Ogłoza, W. 2001, AJ, 122, 402 [NASA ADS] [CrossRef] [Google Scholar]
 Madler, J. H. 1838, C. R. Acad. Sci., 6, 290 [Google Scholar]
 Mante, R. 2000, IAU Commission on Double Stars, 140, 1 [NASA ADS] [Google Scholar]
 Mardia, K., & Jupp, P. 2000, Directional statistics (New York: Wiley) [Google Scholar]
 Martin, C., & Mignard, F. 1998, A&A, 330, 585 [NASA ADS] [Google Scholar]
 Martinache, F., Lloyd, J. P., Ireland, M. J., Yamada, R. S., & Tuthill, P. G. 2007, ApJ, 661, 496 [NASA ADS] [CrossRef] [Google Scholar]
 Mason, B. D., & Hartkopf, W. I. 2011, IAUDS, 174, 1 [Google Scholar]
 Mason, B. D., McAlister, H. A., Hartkopf, W. I., & Shara, M. M. 1995, AJ, 109, 332 [NASA ADS] [CrossRef] [Google Scholar]
 Mason, B. D., Douglass, G. G., & Hartkopf, W. I. 1999, AJ, 117, 1023 [NASA ADS] [CrossRef] [Google Scholar]
 Mazeh, T., Latham, D. W., Goldberg, E., et al. 2001, MNRAS, 325, 343 [NASA ADS] [CrossRef] [Google Scholar]
 McCarthy, Jr., D. W. 1983, in IAU Colloq. Nearby Stars and the Stellar Luminosity Function, eds. A. G. D. Philip, & A. R. Upgren, 76, 107 [Google Scholar]
 Muterspaugh, M. W., Lane, B. F., Fekel, F. C., et al. 2008, AJ, 135, 766 [NASA ADS] [CrossRef] [Google Scholar]
 Muterspaugh, M. W., Hartkopf, W. I., Lane, B. F., et al. 2010a, AJ, 140, 1623 [NASA ADS] [CrossRef] [Google Scholar]
 Muterspaugh, M. W., Lane, B. F., Kulkarni, S. R., et al. 2010b, AJ, 140, 1631 [NASA ADS] [CrossRef] [Google Scholar]
 Muterspaugh, M. W., Lane, B. F., Kulkarni, S. R., et al. 2010c, AJ, 140, 1657 [NASA ADS] [CrossRef] [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]
 O’Brien, D. P., McAlister, H. A., Raghavan, D., et al. 2011, ApJ, 728, 111 [NASA ADS] [CrossRef] [Google Scholar]
 Pompéia, L., Masseron, T., Famaey, B., et al. 2011, MNRAS, 415, 1138 [NASA ADS] [CrossRef] [Google Scholar]
 Popovic, G. M., & Pavlovic, R. 1996, Bull. Obs. Astron. Beograd, 153, 57 [Google Scholar]
 Pourbaix, D. 2000, A&AS, 145, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pourbaix, D., Nidever, D., McCarthy, C., et al. 2002, A&A, 386, 280 [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]
 Pringle, J. E. 1989, MNRAS, 239, 361 [NASA ADS] [Google Scholar]
 Raghavan, D., McAlister, H. A., Henry, T. J., et al. 2010, ApJS, 190, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Raskin, G., van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rees, B., & Zijlstra, A. A. 2013, MNRAS, 435, 975 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, I. N., Hawley, S. L., & Gizis, J. E. 1995, AJ, 110, 1838 [NASA ADS] [CrossRef] [Google Scholar]
 Rica Romero, F. M. 2008, IAU Commission on Double Stars, 165, 1 [NASA ADS] [Google Scholar]
 Russell, J. L., & Gatewood, G. D. 1984, PASP, 96, 429 [NASA ADS] [CrossRef] [Google Scholar]
 Scardia, M. 2002, IAU Commission on Double Stars, 148, 1 [NASA ADS] [Google Scholar]
 Scardia, M., Argyle, R. W., Prieur, J.L., et al. 2007, Astron. Nachr., 328, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Scardia, M., Prieur, J.L., Pansecchi, L., et al. 2008, Astron. Nachr., 329, 54 [NASA ADS] [CrossRef] [Google Scholar]
 Scarfe, C. D., Funakawa, H., Delaney, P. A., & Barlow, D. J. 1983, J. Roy. Astron. Soc. Can., 77, 126 [NASA ADS] [Google Scholar]
 Ségransan, D., Delfosse, X., Forveille, T., et al. 2000, A&A, 364, 665 [Google Scholar]
 Seymour, D. M., Mason, B. D., Hartkopf, W. I., & Wycoff, G. L. 2002, AJ, 123, 1023 [NASA ADS] [CrossRef] [Google Scholar]
 Söderhjelm, S. 1999, A&A, 341, 121 [Google Scholar]
 Stephens, M. 1964, J. Amer. Statist. Assn., 59, 160 [CrossRef] [Google Scholar]
 Sterzik, M. F., & Tokovinin, A. A. 2002, A&A, 384, 1030 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Strand, K. A. 1969, AJ, 74, 760 [NASA ADS] [CrossRef] [Google Scholar]
 Strand, K. A. 1977, AJ, 82, 745 [NASA ADS] [CrossRef] [Google Scholar]
 ten Brummelaar, T., Mason, B. D., McAlister, H. A., et al. 2000, AJ, 119, 2403 [NASA ADS] [CrossRef] [Google Scholar]
 Tokovinin, A. 2008, MNRAS, 389, 925 [NASA ADS] [CrossRef] [Google Scholar]
 Tokovinin, A. A., & Smekhov, M. G. 2002, A&A, 382, 118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tokovinin, A., Kiyaeva, O., Sterzik, M., et al. 2005, A&A, 441, 695 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Turner, J. A., Chapman, S. J., Bhattal, A. S., et al. 1995, MNRAS, 277, 705 [NASA ADS] [CrossRef] [Google Scholar]
 Underhill, A. B. 1963, Publ. DAO, 12, 159 [Google Scholar]
 Upton, G., & Fingleton, B. 1989, Spatial Data Analysis by Example (New York: Wiley), 2 [Google Scholar]
 van Albada, G. B. 1957, Contributions from the Bosscha Observervatory, 5, 1 [NASA ADS] [Google Scholar]
 van den Bos, W. H. 1960, J. Observateurs, 43, 145 [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vinter Hansen, J. M., Neubauer, F. J., & RoosenRaad, D. 1940, Lick Observatory Bulletin, 19, 89 [NASA ADS] [Google Scholar]
 Watson, G. 1956, Geophys. J. Int., 7, 160 [NASA ADS] [CrossRef] [Google Scholar]
 Watson, G. 1966, J. Geology, 74, 786 [NASA ADS] [CrossRef] [Google Scholar]
 Watson, G. 1983, Statistics on Spheres (New York: Wiley) [Google Scholar]
 WiethKnudsen, N. 1957, IAU Commission on Double Stars, 13, 1 [Google Scholar]
 Woolley, R. 1961, The Observatory, 81, 203 [NASA ADS] [Google Scholar]
 Zirm, H. 2003, IAU Commission on Double Stars, 151, 1 [NASA ADS] [Google Scholar]
 Zirm, H. 2011, J. Double Star Observations, 7, 24 [NASA ADS] [Google Scholar]
All Tables
Values of the correction factor c(α) for the Rayleigh test, from Upton & Fingleton (1989).
Various isotropy indicators, listed separately for the 20 systems up to 8.1 pc, and for all 51 systems up to 18 pc.
Crossidentifications for the master list of 95 systems closer than 20 pc from the Sun and the are known to have an orbit in the 6th COVBS at USNO.
Position of the 51 orbital poles in Galactic coordinates for systems closer than 18 pc.
HERMES and CORAVEL radial velocities for components of some visual binaries lacking velocities so far.
New combined astrometricspectroscopic orbits computed from CORAVEL and/or HERMES velocities.
New combined astrometricspectroscopic orbits (2 observable spectra) computed from CORAVEL or HERMES velocities.
CORAVEL velocities (labelled C in the last column), Tokovinin & Smekhov (2002) velocities (labelled T), and Chubak et al. (2012) velocities (labelled M) for the AB components of WDS 00057+4549AB = HIP 473.
New CORAVEL velocities for the components of WDS 01418+4237AB in complement of velocities published by Duquennoy & Mayor (1991).
New CORAVEL velocities for the components of WDS 16413+3136AB in complement of velocities published by Duquennoy & Mayor (1991).
New CORAVEL velocities for the components of WDS 18070+3034AB in complement of velocities published by Duquennoy & Mayor (1991).
The different ways to determine the ascending node of the relative orbit (B/A).
All Figures
Fig. 1 Positions of the orbital poles on the sky, with the grid representing Galactic coordinates. The Galactic centre is at the centre of the grid. The systems are labelled by their numbers of Table 3. The circled labels identify the systems up to N = 18 whose distribution is the most significantly different from isotropic according to Fig. 3. The dashed line displays the celestial equator. For comparison, the ecliptic pole is located at l = 96.4°,b = 30°. 

In the text 
Fig. 2 Top panel: distribution of the orbitalpole positions on the sky, drawn from the normal distributions for i = 54.9° ± 2.4° and Ω = 13.5° ± 2.3° corresponding to the benchmark system WDS 00057+4549. The intersection of the dashed lines marks the nominal pole position. As expected, the distribution is centred on this nominal position. Bottom panel: histograms corresponding to the longitude and latitude distributions. The dashed lines correspond to Gaussian distributions centred on the nominal pole position, with standard deviations of 2.3° and 2.0° in Galactic longitude and latitude, respectively. 

In the text 
Fig. 3 Upper panel: evolution of the RayleighWatson (blue dashed line, and righthand scale; the correction to – Eq. (4) – has been made for a confidence level α of 1%) and Bingham ℬ (red dotdashed line, and righthand scale) statistics as a function of the number of systems. The horizontal blue dashed curve corresponds to a firstkind risk of 1% of rejecting the null hypothesis of isotropy while it is true, based on the statistics. The horizontal red dashed line is the 10% firstkind risk for the ℬ statistics. The number of systems as a function of distance is given by the black solid line and lefthand scale. Isotropy corresponds to ℬ → 0 and . Lower panel: three eigenvalues of the orientation matrix (Eq. (6)). Isotropy corresponds to λ_{1,2,3} = 1 / 3, as seen at large distances. The critical values for λ_{1,3} corresponding to α = 1% and 10% as given by Eqs. (11) and (12) are displayed by the dashdotted lines. 

In the text 
Fig. 4 Evolution of the RayleighWatson statistics as a function of the number of systems in a jackknife approach. Each curve thus corresponds to the full sample with one system removed. The thick line corresponds to the full sample (same as the blue dashed line in the upper panel of Fig. 3). The lowest thin (red) line corresponds to the sample with system No. 16 removed (WDS 14514+1906). Falsealarm probabilities of 2.5, 1, and 0.1% are depicted by the dashed horizontal lines. 

In the text 
Fig. 5 Location of binary systems from Table 2 in Galactic coordinates, with the galactic centre in the middle of the map. The circles identify the systems with an unambiguously defined pole, as plotted in Fig. 1. The dashed line is the celestial equator. 

In the text 
Fig. 6 Distribution of cosi for all binaries from Table 2 (thick lines) and for those with the poles analysed (thin line in the bottom panel and shaded histogram in the upper panel). 

In the text 
Fig. 7 Same as Fig. 1, but for the binary stars without indications for higher multiplicity. 

In the text 
Fig. 8 Same as Fig. 3, but restricted to binary systems without evidence for higher multiplicity. The horizontal blue dashed and red dashdotted curves correspond to a firstkind risk of 5% and 10%, respectively, of rejecting the null hypothesis of isotropy while it is true, based on the statistics. The systems have the same number labels as in Table 3 and Fig. 3. 

In the text 
Fig. 9 Distribution of the systems in the (U,V) plane, with the same labels as in Fig. 1. System 17 (01083+5455) falls outside the boundaries, at U = −43 km s^{1} and V = −157 km s^{1}. To fix the ideas, two velocity ellipsoïds are shown, both centred on the reflex solar motion in the LSR: −U_{⊙} = −10.2 km s^{1}, −V_{⊙} = −20.5 km s^{1}, and with σ_{U} = 34.5 km s^{1}, σ_{V} = 22.5 km s^{1}, as derived by Famaey et al. (2005) from their most precise sample. Horizontal bands and the corresponding labels are used in Fig. 10. 

In the text 
Fig. 10 Same as Fig. 1 but with a letter coding the value of the galactic V velocity (see Fig. 9). A necessary condition for binary systems to originate from the same cluster/association at birth is to have similar V velocities (see text), hence they should be labelled by the same letter. 

In the text 
Fig. A.1 Top panel: radial velocities for the two components of the pair HIP 473 AB = WDS 00057+4549AB. Dots refer to CORAVEL velocities and triangles to velocities from Tokovinin & Smekhov (2002) or Chubak et al. (2012) for the last point, with the time spanned by the measurements represented by the horizontal bar. Open symbols refer to component A and filled symbols to component B. The parallel dashed lines are just a guide to the eye. Bottom panel: radialvelocity difference BA, with their errors taken as the rootmeansquare of the errors on A and B. 

In the text 
Fig. A.2 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 01083+5455Aa,Ab = HIP 5336 = HD 6582 = μ Cas. The solid line in the lower panel corresponds to the predicted velocity curve for A. Triangles are CORAVEL measurements, squares are measurements from Abt & Willmarth (2006). Individual measurements for the visual orbit were taken from the WDS database. 

In the text 
Fig. A.3 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 02278+0426AB. In the lower panel, filled and open triangles correspond to CORAVEL velocities for components 1 and 2, respectively. Individual measurements for the visual orbit were taken from the WDS database at USNO. 

In the text 
Fig. A.4 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 09006+4147AB. Filled triangles in the lower panel refer to component A. Individual measurements for the visual orbit were taken from the WDS database. 

In the text 
Fig. A.5 Crosscorrelation function of the spectrum of WDS 09006+4147AB (=HD 76943) obtained with the HERMES/Mercator spectrograph (Raskin et al. 2011) on HJD 2 456 244.717 (November 13, 2012), correlated with a F0V template. Both components (separated by 0.5′′ on the sky) enter the fiber (2.5′′ on the sky). Component A corresponds to the wide peak centred at 31.5 km s^{1}. 

In the text 
Fig. A.6 Combined astrometric (upper panel) and spectroscopic (lower panel) orbits for WDS 16413+3136AB. The solid line in the lower panel corresponds to the predicted velocity curve for A. Individual measurements for the visual orbit were taken from the WDS database at USNO. 

In the text 
Fig. A.7 Combined astrometric (upper panel) and spectroscopic (lower panel) orbit for WDS 17121+4540AB. The solid line in the lower panel corresponds to the predicted velocity curve for A. Individual measurements for the visual orbit were taken from the WDS database at USNO. 

In the text 
Fig. A.8 Combined astrometric (upper panel) and spectroscopic (lower panel) orbit for WDS 18070+3034AB. The solid line in the lower panel corresponds to the predicted velocity curve for A. Individual measurements for the visual orbit were taken from the WDS database at USNO. In the lower panel, triangles correspond to CORAVEL data, squares to velocities from Abt & Willmarth (2006), diamonds to Tokovinin & Smekhov (2002), and double triangles to the old radial velocities collected by Kamper & Beardsley (1986). 

In the text 
Fig. A.9 Crosscorrelation function of the HERMES spectrum of WDS 19255+0307Aa,Ab (obtained at JD 2 455 657.742) with a F0V template. Component Aa shows a rotationally broadened profile with V_{r}sini = 84 km s^{1}, and component Ab is visible as a weak Gaussian (at V_{r}(Ab) = –37.70 km s^{1}) superimposed on the rotationally broadened profile (centred on V_{r}(Aa) = –26.68 km s^{1}). Similar profiles have been observed over several consecutive nights. 

In the text 
Fig. B.1 True and the projected relative orbit (B/A) of a visual binary and its geometrical elements. The ascending node of the orbit is located at point N. 

In the text 
Fig. B.2 Examples of lifting the ambiguity on the ascending node for a visual binary. Top (cases 2, 5 of Table B.1): the slope of the computed curve is opposite to that of the V_{r}(B) − V_{r}(A) measurements: it is then necessary to add 180° to the values of Ω and ω given by the 6th Catalog of Orbits of Visual Binary Stars (6th COVBS; Hartkopf et al. 2001). Bottom (cases 3, 6 of Table B.1): the slope of the computed curve is opposite to that of V_{r}(A) measurements; there is no need to change the values of Ω and ω listed in the 6th COVBS. 

In the text 
Fig. B.3 Orientation of the true and apparent relative orbits (B/A) with respect to the celestial sphere and the corresponding reference frames. (α,δ) are the equatorial coordinates of the binary system. 

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.