A&A 441, 1135-1148 (2005)
DOI: 10.1051/0004-6361:20053298

On metal-deficient barium stars and their link with yellow symbiotic stars[*],[*]

A. Jorissen1,[*] - L. Zacs2 - S. Udry3 - H. Lindgren4 - F. A. Musaev5,6,7


1 - Institut d'Astronomie et d'Astrophysique, Université Libre de Bruxelles, CP 226, Boulevard du Triomphe, 1050 Bruxelles, Belgium
2 - Institute of Atomic Physics and Spectroscopy, University of Latvia, Raina bulvaris 19, R{\={\i}}\kern.15emga, LV 1586 Latvia
3 - Observatoire de Genève, 1290 Sauverny, Suisse
4 - Lund Observatory, Box 43, 221 00 Lund, Sweden
5 - Special Astrophysical Observatory and Isaac Newton Institute of Chile, SAO Branch, Nizhnij Arkhyz 369167, Russia
6 - ICAMER, National Academy of Sciences of Ukraine, 361605 Peak Terskol, Kabardino-Balkaria, Russia
7 - Shamakhy Astrophysical Observatory, National Academy of Sciences of Azerbaijan, Azerbaijan

Received 25 April 2005 / Accepted 20 May 2005

Abstract
This paper addresses the question of why metal-deficient barium stars are not yellow symbiotic stars (YSyS). Samples of (suspected) metal-deficient barium (mdBa) stars and YSyS have been collected from the literature, and their properties reviewed. It appears in particular that the barium nature of the suspected mdBa stars needs to be ascertained by detailed abundance analyses. Abundances are therefore derived for two of them, HD 139409 and HD 148897, which reveal that HD 148897 should not be considered a barium star. HD 139409 is a mild barium star, with overabundances observed only for elements belonging to the first s-process peak (Y and Zr). It is only moderately metal-poor ([Fe/H] = -0.4). The evidence for binarity among mdBa stars is then reviewed, using three different methods: (i) radial-velocity variations (from CORAVEL observations), (ii) Hipparcos astrometric data, and (iii)  a method based on the comparison between the Hipparcos and Tycho-2 proper motions. An orbit is obtained for HIP 55852, whereas evidence for the (so far unknown) binary nature of HIP 34795, HIP 76605, HIP 97874 and HIP 107478 is presented. No conclusion regarding the binary nature of HIP 11595, HIP 25161 could be reached. Two stars with no evidence for binarity whatsoever (HIP 58596 and BD +3$^\circ $2688) are candidates low-metallicity thermally-pulsing asymptotic giant branch stars, as inferred from their large luminosities. The reason why mdBa stars are not YSyS is suggested to lie in their different orbital period distributions: mdBa stars have on average longer orbital periods than YSyS, and hence their companion accretes matter at a lower rate, for a given mass loss rate of the giant star. The definite validation of this explanation should nevertheless await the determination of the orbital periods for the many mdBa stars still lacking periods, in order to make the comparison more significant.

Key words: binaries: symbiotic - stars: abundances - stars: AGB and post-AGB - binaries: spectroscopic

1 The problem

Our understanding of the link between chemically-peculiar red giants like barium stars or CH stars, and yellow symbiotic stars (YSyS) has made substantial progress in the last decade (see the reviews by Jorissen 2003b,a), mainly with the realisation that likely all yellow symbiotics (i.e., involving a giant of spectral type G or K as primary component) involve barium stars (Smith et al. 1996, 1997; Pereira & Porto de Mello 1997; Pereira et al. 1998). The metal-deficient nature of the giant is a key factor, because it implies a rather large luminosity for the giant. Since evolutionary tracks of metal-deficient stars are shifted towards the blue, metal-deficient giants of spectral type K must lie on the upper part of the (asymptotic) giant branch (see Fig. 11 of Smith et al. 1996), where they suffer strong mass loss. If such metal-deficient K giants are in binary systems, their strong wind will interact with the companion and trigger symbiotic activity.

Table 1: The barium syndrome among YSyS. The spectral type of the cool component is taken from Mürset & Schmid (1999), or references therein. In column labelled ``nebula'', ``y'' means that an optical nebula has been detected, and  ``PN'' that, based on its emission line spectrum, the star has traditionally been included in planetary nebulae catalogues, even though no optical nebula may be visible. The column labelled "P'' lists the orbital period, from Mürset & Schmid (1999).

The other facet of this problem, namely whether all metal-deficient barium stars are symbiotic stars, is not yet fully answered. The present paper offers a first step in that direction.

We have collected a list of candidate metal-deficient barium stars and have assembled new observations to check (i) whether these stars are binaries, (ii) whether they are barium stars and (iii) whether they exhibit symbiotic activity.


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{3298f1a.eps}\hspace*{4mm}
\includegraphics[width=7.8cm,clip]{3298f1b.eps}\end{figure} Figure 1: Left panel a): evolutionary tracks of Schaller et al. (1992) compared with the locations of classical barium stars (stars labelled "G'' in Mennessier et al. (1997); open dots) and the yellow SyS AG Dra (filled circle), BD  $-21^\circ 3873$ (filled triangle) and Hen 2-467 (filled square). The bolometric magnitudes were taken from the references listed in Table 1. These bolometric magnitudes were combined with bolometric corrections from Bessell et al. (1998) and B-V indices from Munari et al. (1992) and Munari & Buson (1992) to yield the absolute visual magnitudes. Right panel b): same as a) but for YSyS (filled symbols as in the left panel) and metal-deficient barium stars [open triangles: stars flagged as "H'' by Mennessier et al. (1997); crosses: CH stars also flagged as "H'' by Mennessier et al. (1997); small open squares: additional metal-deficient, s-process-rich stars from Table 2; open circle: HD 104340, large open square: HD 206983 from Junqueira & Pereira (2001)]. The dashed horizontal line represents the luminosity ( $M_{\rm bol} = -3$, corresponding to $M_{\rm V} \sim -2$) at the first thermal pulse in a 1 $M_\odot $ AGB star of metallicity [Fe/H] = -1.8according to Lattanzio (1991).
Open with DEXTER

2 The samples

Before discussing metal-deficient barium stars, it is useful to first summarize the properties of YSyS, to which metal-deficient barium stars may be compared.

2.1 Yellow symbiotic stars

All known YSyS are listed in Table 1, which shows that all the stars studied so far exhibit the barium syndrome. YSyS with a stellar infrared continuum (s-type, as opposed to the dusty d'-type; see below) are clearly metal-deficient objects, as revealed by their low metallicities and high space velocities (CD  $-43^\circ14304$ may be an exception; however, it is of spectral type K7, and should perhaps not be included in the family of YSyS). The presence of the barium syndrome among a family of binary metal-deficient stars fully supports the commonly accepted hypothesis that the s-process is more efficient at low metallicities (Clayton 1988; Jorissen 2003a). s-Type YSyS, with their metallicities lower than classical barium stars, may be expected to be, on average, more luminous than the latter (see Fig. 11 of Smith et al. 1996, comparing the luminosity function of Pop.I and Pop.II K giants). This is a direct consequence of the fact that evolutionary tracks shift towards the blue in the Hertzsprung-Russell (HR) diagram as metallicity decreases, as shown in Fig. 1b. Figure 1a confirms that the YSyS AG Dra and BD  $-21^\circ 3873$ are indeed more luminous than classical barium stars. This difference in the average luminosity - and hence mass-loss rate - of the two populations thus explains why YSyS, despite hosting a K giant, exhibit symbiotic activity whereas barium stars do not. The larger mass-loss rates for the cool components of s-type YSyS - as compared to Ba stars - may be inferred from the comparison of their IRAS [12]-[25] color indices, which reflect the amount of dust present in the system: ([12]-[25]) $_{\rm Ba} < $ 0.1, as compared to 0.45 for AG Dra (Smith et al. 1996). Mürset et al. (1991) and Drake et al. (1987) provide direct measurements (or upper limits) for the mass loss rates of AG Dra and of Ba stars, respectively, which confirm the above conclusion.

YSyS with a dusty infrared continuum (d'-type; Schmid & Nussbaumer 1993; Allen 1982) differ from their s-type counterparts in several respects (Table 1): they host a complex circumstellar environment (including cool dust, bipolar outflows, extended optical nebulae or emission-line spectra closely resembling those of planetary nebulae), the cool components have early spectral types (F to early K), they are often fast rotators (with the possible exception of M 1-2 = V471 Per; Grauer & Bond 1981) and, finally, they belong to the galactic disk unlike s-type YSyS which belong to the halo.

All these arguments suggest that the hot component in d'-type SyS has just evolved from the AGB to the WD stage. The rather cool dust (Schmid & Nussbaumer 1993) is a relic from the mass lost by the AGB star. The optical nebulae observed in d'-type SyS are most likely genuine planetary nebulae rather than the nebulae associated with the ionized wind of the cool component (Corradi et al. 1999). This is especially clear for AS 201 which actually hosts two nebulae (Schwarz 1991): a large fossil planetary nebula detected by direct imaging, and a small nebula formed in the wind of the current cool component. Finally, the rapid rotation of the cool component has likely been caused by spin accretion from the former AGB wind like in WIRRING systems (Jeffries & Stevens 1996; Jorissen 2003b). The fact that the cool star has not yet been slowed down by magnetic braking is another indication that the mass transfer occurred fairly recently (Theuns et al. 1996). Corradi & Schwarz (1997) obtained 4000 y for the age of the nebula around AS 201, and 40 000 y for V417 Cen.

2.2 Metal-deficient barium stars

Metal-deficient barium stars (with metallicities in the range -1.1 to -1.8, comparable to that of YSyS) were identified by Luck & Bond (1991), Mennessier et al. (1997) and Zacs et al. (2000), and occupy the same region of the HR diagram as YSyS (Fig. 1b). The question thus arises why metal-deficient barium stars are not YSyS. Different answers must be sought, depending upon their absolute visual magnitudes $M_{\rm V}$. The most luminous systems, with $M_{\rm V} < -2$, are likely located on the thermally-pulsing AGB, so that their Ba syndrome may be explained by internal nucleosynthesis. They thus should not be binaries, and therefore cannot be YSyS! HD 104340 (open circle in Fig. 1b), a metal-deficient Ba star studied by Junqueira & Pereira (2001), and BD + $03^\circ 2688$ (Table 2) provide good illustrations of this situation, since they both lie above the TP-AGB threshold and CORAVEL radial-velocity measurements spanning several years do not reveal any clear orbital motion (Figs. 9 and 10, as well as Sect. 4).

Table 2: Stars classified as metal-deficient barium stars by Mennessier et al. (1997). The absolute magnitude MV is a maximum likelihood estimate obtained by Mennessier et al. (1997), except for the additional stars where it is derived from a straight inversion of the Hipparcos parallax (BD  $+75^\circ 348$), from a spectroscopic estimate of gravity and an educated guess for the mass (BD + $3^\circ 2688$) or from a fit to the M 92 isochrone (CS stars). The column labelled "Ba'' indicates whether detailed chemical analyses have confirmed the Ba nature of the star.

The less luminous and warmest among metal-deficient Ba stars, clumping around $M_{\rm V} \sim +1$ in the HR diagram, are also sometimes classified as CH stars (crosses in Fig. 1b). They are not losing mass at a large enough rate to trigger any symbiotic activity, as revealed by their small [12]-[25] color indices (<0.3; Smith et al. 1996).

Finally, at intermediate luminosities ( $-2 \le M_{\rm V} \le
+1$) where YSyS are located, metal-deficient Ba stars are not luminous enough to be TP-AGB (hence they should be binaries), but yet their mass loss rates must be large enough to trigger symbiotic activity (provided that the orbital separation is not too large, since it is the mass accretion rate by the compact companion which is in fact the key parameter; see Sect. 5 and Jorissen 2003a). It is thus of great interest to check (i) the Ba nature of those metal-deficient stars with intermediate luminosities, (ii) their binary nature, and (iii) their suspected symbiotic activity. The first two issues are addressed in Sects. 3 and 4, respectively.

As far as a possible symbiotic activity is concerned, there is no indication from their photometric U-B and B-V indices that the metal-deficient stars in Table 2 have a strong blue continuum which could betray their symbiotic nature. It is thus very likely that none among these stars is a full-fledged symbiotic star. No signature of weak symbiotic activity (of the kind exhibited by some binary S stars; see Fig. 14 and Van Eck & Jorissen 2002) is observed in the H$_\alpha $ line profile either (Fig. 2).


  \begin{figure}
\par\includegraphics[width=7.6cm,clip]{3298f2.eps}\end{figure} Figure 2: $H_{\alpha }$ profiles for two candidate metal-deficient Ba stars (HD 139409 and HD 148897; see however the discussion in Sect. 3), compared to that of the K giant Arcturus. The wavelength scale has been corrected for the stellar radial velocity. The profiles show no indication of symbiotic activity.
Open with DEXTER

   
3 Abundances

The classification of the stars in Table 2 as metal-deficient Ba stars is subject to caution, as it does not rely on spectroscopic abundance analyses, but rather on a maximum-likelihood assignment based on kinematic, spatial and luminosity properties (Mennessier et al. 1997) for barium stars from the list of Lü (1991). Nevertheless, when a metallicity determination is available, it confirms the metal-deficient nature of the object (see Table 2). HD 139409 is an exception, though, since detailed spectral analyses reveal that it is neither metal-poor nor strongly enriched in s-process elements (see below). It may nevertheless be hoped that the metal-deficient assignment made by Mennessier et al. (1997) is valid in all the other cases. Regarding the Ba nature of these stars, it is known that the Lü (1991) catalogue of barium stars, from which the sample of barium stars used by Mennessier et al. (1997) was drawn, is contaminated by many non-barium stars (Jorissen et al. 1996; Griffin & Keenan 1992), especially among those stars having a Ba index smaller than 1.

It would therefore in principle be necessary to re-evaluate the Ba nature of all the stars listed in Table 2. So far, spectra could be obtained for two of them, HD 139409 and HD 148897, which are discussed in detail in the present section. Of these, only HD 139409 appears to be a mild barium star, thus confirming the suspicion about the Lü (1991) catalogue expressed above.

3.1 The case of HD 148897 and HD 139409

HD 148897 (=HR 6152) has been tagged as a "likely marginal barium star'' by Boyle & McClure (1975) and as a marginal CH star by Vilhu et al. (1977). It therefore found its way into the barium-star catalogs of Lü et al. (1983) and Lü (1991), as well as the catalog of CH and metal-deficient Ba stars of Sleivyte & Bartkevicius (1990). The star was classified as G8.5III CN-2 Fe-1 CH-1 by Keenan & McNeil (1989), and this classification as CN- and CH-weak contradicts the earlier assignments.

A detailed abundance analysis has thus been performed to clarify the situation, and its results are compared with previous analyses by Kyröläinen et al. (1986) and Luck (1991) in Sect. 3.1.2.

HD 139409 (=HIP 76605) has been classified as a marginal barium star by MacConnell et al. (1972), as G5 III Ba1 by Yamashita & Norimoto (1981) and as K0 III/II Ba 0.5 by Lü (1991).

3.1.1 Observations

A high-resolution spectrum of HD 148897 was obtained using the Coudé Matrix Echelle Spectrometer (MAESTRO; Musaev et al. 1999) delivering a resolving power of 45000 and installed on the 2-m Zeiss telescope of the Terskol Observatory (located in Northern Caucasus at an altitude of 3100 m). The spectrometer is equipped with a Wright Instruments CCD detector with 1242 $\times$ 1152 pixels (22.5 $\times$ 22.5 $\mu$m). A total exposure of 1800 s was taken on February 18, 2003. The spectrum covers the range 365 to 1020 nm spread over 85 spectral orders.

A high-resolution spectrum of HD 139409 was obtained on the HARPS spectrograph (Mayor et al. 2003), delivering a resolution of 115000 and installed on the ESO La Silla 3.6 m telescope. A total exposure of 200 s was obtained by Xavier Bonfils on May 31, 2004. In order to check equivalent widths delivered by HARPS, a spectrum of the standard star Arcturus was obtained as well by Fabien Carrier. The HARPS spectra were reduced by the observers using standard pipeline processing. Equivalent widths for the same set of lines as those studied in HD 139409 have been measured by one of us (L.Z.) in the HARPS spectrum of Arcturus and compared to those from the Arcturus spectral atlas (Hinkle et al. 2000). The agreement between the two sets of equivalent widths is excellent (Fig. 3), thus qualifying HARPS for abundance analyses. Spectra around the $\lambda614.172$ nm Ba II line are shown in Fig. 4 for the two analyzed stars and for Arcturus (Hinkle et al. 2000).


  \begin{figure}
\par\includegraphics[width=7.4cm,clip]{3298f3.eps}\end{figure} Figure 3: Comparison of equivalent widths measured on the HARPS spectrum of Arcturus with those of the Arcturus spectral atlas (Hinkle et al. 2000), for the same set of lines as that measured in HD 139409.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.6cm,clip]{3298f4.eps}\end{figure} Figure 4: The spectra of the target stars in the region of the Ba II line at 614.172 nm. Measured lines are marked. Also shown is the spectrum of the standard star Arcturus (Hinkle et al. 2000).
Open with DEXTER

   
3.1.2 Reduction and analysis

HD 148897

The reduction of the CCD frames (subtraction of bias, dark and scattered light, flat fielding, extraction of echelle orders and wavelength calibration) was performed with the DECH20T software (Galazutdinov 1992). More than 500 weak to medium-strong atomic lines, free of blends, were identified and their equivalent widths were measured with the DECH routines. The radial velocity was measured using a large number of symmetric absorption lines. The observed velocity has been brought to the heliocentric system by adding +22.7 km s-1. The mean heliocentric radial velocity for HD 144897 was found to be +16.7 km s-1. The atmospheric parameters for that star cannot be derived from photometry, since the standard temperature calibrations refer to stars of normal chemical composition. To obtain a colour-independent estimate of the temperature, a spectroscopic temperature has been derived from the excitation equilibrium of Fe I lines. The surface gravity $\log g$ was determined using the FeI/FeII ionization balance, whereas the microturbulent velocity $v_{\rm t}$ was derived by forcing the abundances of individual Fe I lines to be independent of the reduced equivalent width. The resulting atmospheric parameters for HD 148897 are as follows: $T_{\rm eff} = 4350$ K, $\log g = 1.0$ (cgs), and $v_{\rm t} = 2.0$ km s-1. An independent determination of the surface gravity, using the above $T_{\rm eff}$ value, Mennessier et al. (1997) absolute visual magnitude ( MV = -1.8), a bolometric correction of -0.5 and a mass of 1 $M_\odot $, yields $\log g = 1.1$, in agreement with the adopted value (see Luck 1991, for a discussion of the discrepancy usually observed between the photometric and spectroscopic gravities). Comparison with atmospheric parameters from the literature is presented in Table 3.

The abundance analysis has been performed with the standard LTE line analysis program WIDTH9 developed by Kurucz. The model atmospheres were taken from Gustafsson et al. (1975). The synthetic spectra were generated using the spectral synthesis code STARSP (Tsymbal 1996). Oscillator strengths have been taken from the VALD database (Piskunov et al. 1995). The resulting abundances (normalized by the solar-system abundances of Grevesse & Sauval 1998) are listed in Table 3, from which it may be concluded that HD 148897 appears to be a rather typical metal-deficient star (McWilliam 1997), and should certainly not be considered as a (metal-deficient) barium star.

Interestingly enough, there are no indications whatsoever from the Hipparcos and Tycho data (applying the methods described in Sect. 4) that this star is binary, in agreement with the fact that it is not a barium star.

Table 3: Averaged absolute abundances (in the scale where $\log \epsilon$(H) = 12) and normalized abundances [X/Fe] for HD 148897 (relative to the Sun; Grevesse & Sauval 1998). The standard deviations $\sigma $ and the number of lines used in the analysis (n) are also given. Two sets of abundances from Luck (1991) are given: one is based on a "photometric $\log g$'' (derived from the absolute magnitude, the effective temperature and an estimate for the stellar mass) and the other on the "spectroscopic $\log g$'' (Fe I/Fe II ionization balance).

HD 139409

The effective temperature of HD 139409 has been derived from the excitation equilibrium of Fe  I, Ti  I and Cr  I lines (see Fig. 5, for Fe). The surface gravity $\log g$ was determined from the Fe  I/Fe  II ionization balance, and the microturbulent velocity by forcing the abundances of individual Fe  I, Ti  I and Cr  I lines to be independent of equivalent width (see Fig. 6, for Fe). Although spectroscopic gravity and temperature determinations in late-type, metal-deficient stars are probably affected by non-LTE effects, these effects remain small when [Fe/H] $\ge -1.0$ (see, for example, Allende Prieto et al. 1999). The stellar parameters of Arcturus ([Fe/H] = -0.6, $\log g = 1.3$) derived by the spectroscopic method (applied on the HARPS spectrum) are in good agreement indeed with those derived by other (non-spectroscopic) methods.


  \begin{figure}
\par\includegraphics[width=7.9cm,clip]{3298f5.eps}\end{figure} Figure 5: The iron abundances derived from Fe I (open circles) and Fe II (filled circles) lines are displayed as a function of the excitation potential for the lower energy level of the line. The absence of a trend in these data is used to derive the spectroscopic temperature of HD 139409.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.9cm,clip]{3298f6.eps}\end{figure} Figure 6: The abundances derived from the Fe I (open circles) and Fe II (filled circles) lines are displayed as a function of the line equivalent widths. The absence of a trend in these data is used to derive the microturbulent velocity of HD 139409.
Open with DEXTER

Table 4: Abundances for HD 139409, in the scale where $\log \epsilon$(H) = 12, and normalized with respect to the solar abundances (Grevesse & Sauval 1998). The standard deviations $\sigma $ and the number N of lines used in the analysis are also given.


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{3298f7.eps}\end{figure} Figure 7: Abundance pattern in HD 139409 (abundances derived from ionized species are represented by black dots, and from neutral species by crosses).
Open with DEXTER

The resulting atmospheric parameters for HD 139409 are as follows: $T_{\rm eff} = 5000$ K, $\log g = 2.8$ (cgs), and $v_{\rm t}$ = 2.0 km s-1. The spectroscopic gravity, combined with a mass of 1 $M_\odot $, leads to $M_V \simeq +1.5$. Thus our calculations indicate that HD 139409 is less luminous than predicted by Mennessier et al. (1997). The derived iron abundance is [Fe/H] = -0.42 [adopting $\log(\epsilon(\rm Fe)) = 7.50$].

The basic conclusion from the set of abundances listed in Table 4 and displayed in Fig. 7 is that HD 139409 appears to be a mild barium star[*]. The s-process overabundances observed in HD 139409, although quite moderate, are not much smaller than those observed in a yellow symbiotic star like BD $-21^\circ$3873, which exhibit overabundances of s-process elements in the range 0.3-0.8 dex (Smith et al. 1997). With [Fe/H] = -0.4, HD 139409 is, however, not as metal-deficient as the other stars considered in this paper.

   
4 Statistics of binarity among metal-deficient Ba stars

The list of confirmed or suspected metal-deficient barium stars remaining after the screening process based on abundances (as described in Sect. 3) is given in Table 5. The possible binary nature of those stars may be assessed using three different methods:

These methods are now described in turn.

4.1 Radial-velocity variations

Several stars from Table 5 (namely HIP 4347, HIP 29740, HIP 34795, HIP 55852, HIP 58956, HIP 76605 and BD +3 $^\circ2688$) have been monitored for many years using the CORAVEL spectrovelocimeter (Baranne et al. 1979), as part of a larger program aiming at finding the frequency of spectroscopic binaries among s-process-rich late-type giants (see Jorissen et al. 1998; Jorissen & Mayor 1992,1988, for details and other results from this CORAVEL monitoring). Individual radial-velocity measurements for those stars (in the CORAVEL-ELODIE system as defined in Udry et al. 1999) are given in Table 7 or in Udry et al. (1998). For a few other stars (HIP 43042, CS 22942-019 and CS 22948-027), radial velocities were monitored using other instruments, and their results were taken from the literature (Preston & Sneden 2001; Zacs et al. 2005).

Table 5: Binary properties of confirmed or suspected metal-deficient barium stars. The column labeled "Ba'' indicates whether detailed chemical analyses have confirmed the Ba nature of the star. The columns labeled $\chi ^2_{\mu _{\rm HIP} - \mu _{\rm Tyc}}$ and Prob( $\chi ^2 > \chi ^2_{\mu _{\rm HIP} - \mu _{\rm Tyc}})$ provide the $\chi ^2$ (and its associated probability) involved in the comparison of Hipparcos and Tycho-2 proper motions (see text for details). The column "Bin.'' has been set to "y'' if the first kind risk of rejecting the null hypothesis that the Tycho and Hipparcos proper motions are equal is smaller than 10%. The column labeled IAD indicates whether the signature of an orbital motion is present in the Hipparcos Intermediate Astrometric Data, according to the various tests described in the text. The column labeled  $\sigma (Vr)$ provides the radial-velocity standard deviation, $\Delta t$ and N correspond to the time span and number of observations, respectively. The columns labeled P and "Ref'' list the orbital period (when available) and the reference for the radial velocity and/or orbital data. The column "Binary'' gives the final binarity diagnostics.

Orbits were already available for HIP 4347 and HIP 29740 (Udry et al. 1998), HIP 43042 (Zacs et al. 2005), as well as for CS 22942-019 and CS 22948-027 (Preston & Sneden 2001). A new orbital solution has been derived for BD +04$^\circ $2466 (Table 6 and Fig. 8). The binary nature of those stars is therefore beyond doubt.

The radial-velocity standard deviation of HIP 58596 is larger than expected based on the uncertainty on one measurement (Fig. 9), but no satisfactory orbital solution could be found. A 3-d orbit (with an eccentricity of 0.30) is possible, but this short orbital period is not consistent with the giant nature of HIP 58596, which imposes orbital periods of at least 20 d (see Fig. 4 of Pourbaix et al. 2004). The large standard deviation exhibited by HIP 58596 is therefore very likely another example of the large intrinsic jitter often observed for metal-deficient stars, as discussed by McClure (1984) and Carney et al. (2003).

A 14 y radial-velocity monitoring for HIP 34795 with the northern and southern CORAVELs is not very conclusive either, mostly because there are difficulties in finding the zero-point offset between the two instruments for such large radial velocities (Udry et al. 1999). When a -1 km s-1 offset is applied to the northern velocities (with respect to the values listed in Table 7), a long-term trend seems to be present, albeit with some superimposed jitter (Fig. 11). The analysis of the Hipparcos astrometric data presented in Sect. 4.2 suggests that the star might be binary, although the evidence is not very conclusive.

Finally, there is no sign of radial-velocity variations for BD + $3^\circ 2688$(Fig. 10).

   
4.2 Orbital motion from Hipparcos astrometric data

A tailored reprocessing of the Hipparcos Intermediate Astrometric Data (hereafter IAD; van Leeuwen & Evans 1998) makes it possible to look for a possible orbital signature in the astrometric motion, following the method outlined by Pourbaix & Jorissen (2000), Pourbaix & Boffin (2003), Pourbaix (2004) and applied to barium stars by Jorissen et al. (2004b). We give here only a brief summary of the method.

Table 6: Orbital elements for BD +04$^\circ $2466.


  \begin{figure}
\par\includegraphics[width=7.45cm,clip]{3298f8.eps}\end{figure} Figure 8: Radial velocities as a function of heliocentric Julian Day for BD +4$^\circ $2466, superimposed on the orbital solution corresponding to the orbital elements listed in Table 6.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{3298f9.eps}\end{figure} Figure 9: Radial velocities as a function of heliocentric Julian Day for HIP 58596 = HD 104340 (open circles correspond to measurements obtained with the southern CORAVEL, and filled circles with the northern one).
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.45cm,clip]{3298f10.eps}\end{figure} Figure 10: Radial velocities as a function of heliocentric Julian Day for BD +3$^\circ $2688.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{3298f11.eps}\end{figure} Figure 11: Radial velocities as a function of heliocentric Julian Day for HIP 34795 = HD 55496. The symbols are as in Fig. 9. With respect to the data listed in Table 7, an offset of -1 km s-1 has been applied to the northern measurements to make them consistent with the southern ones.
Open with DEXTER

The basic idea is to quantify the likelihood of the fit of the Hipparcos astrometric data with an orbital model. For that purpose, Pourbaix & Arenou (2001) (see also Jancart et al. 2005) introduced several statistical indicators which allow us to decide whether to keep or to discard an orbital solution. Those indicators relevant to our purpose are the following:

With the above notations, the requirements for a star to qualify as a binary is then

 \begin{displaymath}%
\alpha \equiv \left(Pr_2 + Pr_3\right)/\epsilon \le 0.02,
\end{displaymath} (6)

where the threshold value of 0.02 has been chosen to minimize false detections (Jorissen et al. 2004b).

Hipparcos data are, however, seldom precise enough to derive the orbital elements from scratch. Therefore, when a spectroscopic orbit is available beforehand, it is advantageous to import e, P, T from the spectroscopic orbit and to derive the remaining astrometric elements (as done by Pourbaix & Boffin 2003; Pourbaix & Jorissen 2000). If a spectroscopic orbit is not available, trial (e, P, T) triplets scanning a regular grid (with $10 \le P (\rm d) \le 5000$ imposed by the Hipparcos scanning law and the mission duration) may be used. The quality factor $\alpha $ is then computed for each trial (e, P, T) triplet, and if there exist triplets yielding $\alpha < 0.02$, the star is flagged as a binary. To test its success rate, this method has been applied by Jorissen et al. (2004b) on a sample of barium stars. These authors show that, when $\varpi > 5$ mas and $100 < P (\rm d)< 4000$, the (astrometric) binary detection rate is close to 100%, i.e., the astrometric method recovers all known spectroscopic binaries (see also Jancart et al. 2005). When the orbit is not known beforehand, the method makes it even possible to find a good estimate for the orbital period, provided, however, that the true period is not an integer fraction, or a multiple, of one year. Here the method is applied to the sample of metal-deficient barium stars listed in Table 5.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{3298f13a.eps}\hspace*{3mm}
\includegraphics[width=8cm,clip]{3298f13b.eps}\end{figure} Figure 13: Same as Fig. 12 for HIP 34795 ( left panel, flagged as binary) and HIP 58596 ( right panel, non-binary).
Open with DEXTER

The method flags as definite binaries the stars HIP 29740, 34795, 43042, 97874 and 107478 (Figs. 12 and 13). In two cases (HIP 29740 and 43042), the IAD method thus confirms the conclusion from the radial-velocity monitoring, but yields as well three new binaries (HIP 34795, 97874 and 107478). Two spectroscopic binaries (HIP 4347 and HIP 55852) are not detected by the IAD method because of their small parallax or long orbital period. The non-binary nature of HIP 58596, already suspected from the radial-velocity data, is confirmed by the analysis of the IAD (Fig. 13).

   
4.3 Orbital motion from a comparison of Hipparcos and Tycho-2 proper motions

Kaplan & Makarov (2003) suggested that the comparison of Hipparcos and Tycho-2 (Høg et al. 2000b) proper motions offers a way to detect binaries with long periods (typically from 2000 to 4000 d). The Hipparcos proper motion, being based on observations spanning only 3 y, may be altered by the orbital motion, especially for systems with periods in the range of 2000 to 4000 d whose orbital motion was not recognized by Hipparcos. On the other hand, this effect should average out in the Tycho-2 proper motion, which is derived from observations covering a much longer time span. This method, already used by Makarov (2004), Pourbaix (2004) and Jancart et al. (2006), works best when applied to stars with parallaxes in excess of about 5 mas.

The method evaluates the quantity

\begin{displaymath}%
\chi^2_{\rm obs} = \left(\vec{\mu}_{\rm HIP} - \vec{\mu}_{\...
...-1}\; \left(\vec{\mu}_{\rm HIP} - \vec{\mu}_{\rm Tyc}\right),
\end{displaymath} (7)

where $\vec{\mu}_{\rm HIP}$ and $\vec{\mu}_{\rm Tyc}$ are the vectors of $\alpha $ and $\delta$ components of the Hipparcos and Tycho-2 proper motions, respectively, and W is the associated 2 $\times$ 2 variance-covariance matrix. The covariance between  $\mu_{\alpha,{\rm HIP}}$ and  $\mu_{\delta,{\rm HIP}}$, as provided by field H28 of the Hipparcos catalogue (ESA 1997) and the correlation between Tycho-2 and Hipparcos proper motions, as encapsulated in the quantity R of Table 1 of Høg et al. (2000a), have both been considered (see Jancart et al. 2006, for details).

Since the above quantity follows a $\chi ^2$ probability distribution function with 2 degrees of freedom, it is then possible to compute the probability Prob $(\chi^2 > \chi^2_{\rm obs})$, giving the first kind risk of rejecting the null hypothesis $\vec{\mu}_{\rm Tycho} = \vec{\mu}_{\rm HIP}$ while it is actually true. This probability is listed in Table 5, along with  $\chi^2_{\rm obs}$, and the star is flagged as binary if Prob < 0.1.

Only HIP 29740, HIP 76605 and HIP 97874 satisfy the test at the 10% threshold. Note, however, that all the other stars have parallaxes smaller than 5 mas, which make the test less efficient.

   
5 Summary of the binary criteria and discussion

The situation may be summarized as follows (see also last column of Table 5):

The latter two non-binary stars are in fact good candidate thermally-pulsing AGB stars, as revealed by their location in the HR diagram (Table 2 and right panel of Fig. 1). Hence, they must not be binaries.


  \begin{figure}
\par\includegraphics[width=7.3cm,clip]{3298f14.eps}\end{figure} Figure 14: Comparison of the period distributions for samples of binary systems with different kinds of red giant primaries: S stars (Jorissen et al. 1998), red SyS (excluding symbiotic novae and symbiotic Miras; Mürset & Schmid 1999), M giants (Jorissen et al. 2004a) and K giants (Mermilliod 1996). The lower two panels present the orbital-period distribution for barium stars (Jorissen et al. 1998) and yellow SyS (Mürset & Schmid 1999). In the S star panel, the arrow marked H$_\alpha $ indicates the period range where binary S stars exhibit H$_\alpha $ emission as a signature of weak symbiotic activity. The shaded area marks stars with only a lower limit available on their orbital period. The numbers in the upper right corner of each panel correspond to the sample size.
Open with DEXTER

Finally, we come to the central question of this paper: why do the metal-deficient barium stars, despite being binaries and occupying the same location of the HR diagram as YSyS, do not exhibit symbiotic activity? Three possible answers have been suggested in this paper:

(i)
Some among the stars listed in Table 2 and displayed in Fig. 1 are in fact not barium stars (especially HIP 80843 = HD 148897).

(ii)
Among those which are barium stars, some seem to lie on the TP-AGB, and thus need not be binaries (HIP 58596 = HD 104340, BD +3$^\circ $2688). They therefore cannot exhibit symbiotic activity.

(iii)
Finally, there remain a few genuine metal-deficient barium stars in the sample. Why are they not symbiotic stars? It seems that the answer to that question lies in the different period distributions for YSyS and metal-deficient barium stars: YSyS have shorter orbital period than metal-deficient barium stars, as seen in Fig. 14. This argument seems to apply especially to HIP 29740 (=HD 43389), which has been assigned a very bright absolute visual magnitude of -3.5 by the maximum likelihood method of Mennessier et al. (1997). It is therefore expected to have a rather strong mass loss rate, and be a good candidate YSyS. However, with its orbital period of 1689 d, it lies at the long-period edge of the period distribution of YSyS (Fig. 14). The same difference seems to exist between the period distributions of red symbiotics and binary S stars (Van Eck & Jorissen 2002, and Fig. 14).
However, a firm conclusion on this issue should await the determination of the orbital periods for the metal-deficient barium stars flagged as binaries by the IAD method (especially HIP 11595, HIP 25161, HIP 69834, HIP 97874, HIP 107478), so as to make the comparison between the orbital period distributions of metal-deficient barium stars and YSyS more meaningful.

Acknowledgements
This work was performed in the framework of the NATO Collaborative Linkage Grant SA (PST.CLG.979128)6774/FP. We thank D. Pourbaix for the processing of the Hipparcos data of the metal-deficient stars. F. Carrier and X. Bonfils are thanked for obtaining the Arcturus and HD 139409 spectra. L.Z. thanks I. Platais for valuable discussions and support. The Fonds National Suisse de la Recherche Scientifique has funded the operations of the CORAVEL spectrometer and of the Swiss 1-m telescope installed at the Haute-Provence Observatory.

References

 

  
Online Material

Table 7: Individual CORAVEL radial velocities and associated errors for HIP 34795, HIP 58596, HIP 76605, BD + $3^\circ 2688$ and BD +04$^\circ $2466. The last column indicates which one of the two CORAVEL spectrovelocimeters has been used (NO = CORAVEL north at the Observatoire de Haute Provence; SO = CORAVEL south on the Danish 1.54-m telescope at ESO).



Copyright ESO 2005