A&A 416, 661-668 (2004)
DOI: 10.1051/0004-6361:20034445
1 - INAF - Catania Astrophysical Observatory, via S. Sofia 78, 95123, Catania, Italy
2 -
Istituto di Radioastronomia del CNR, Stazione VLBI di Noto,
CP 161 Noto, Italy
Received 6 October 2003 / Accepted 27 November 2003
Abstract
From a new set of spectroscopic observations we determined orbital parameters of six SB2 systems with one or both components being HgMn stars.
We slightly refined the orbital periods for HD 32964, HD 173524, HD 174933
and HD 216494. Our results for HD 358 are in agreement with the previous literature
studies. Regarding HD 33647, our orbital period is shorter than previous
determinations. HD 173524 is a triple system. From the variations of the
-velocity deduced from our and literature data, we refined the orbital period and
we estimated the eccentricity of the third companion, for which we get
e
0.13 and
= 36
3 years. HD 191110 and HD 216494 seem
to be the only synchronous SB2 within our sample. Following the hypothesis of rotational
axes perpendicular to the orbital plane, we estimated the angle of the orbital plane to the
line of sight. We identified in our spectra the signature of the H
of the third
component of HD 216494 and we highlight a slight variation of the
-velocity due to the presence of this component.
Key words: stars: binaries: spectroscopic - stars: chemically peculiar - stars: general
Among the Chemically Peculiar stars of the main sequence, the HgMn sub-group show anomalies in the intensity of mercury and manganese lines. In particular, spectral lines due to those elements are always stronger than the corresponding lines observed in normal main sequence stars.
For binarity among HgMn stars, there is clear evidence that the incidence of SB systems is similar to that observed in normal stars. Abt & Snowden (1973) carried out the first systematic search for binaries among CP stars and found 43% of their sample of HgMn stars belong to SB system. This percentage is very similar to the incidence of binary among normal stars found by Jaschek & Gomez (1970). Aikman (1976) increased the small sample of Abt & Snowden (1973), raising it from 15 to 80 stars, confirming the conclusion. Floquet (1983) analyzed the catalogue of HgMn stars published by Schneider (1981) and concluded that 54% of these objects have variable radial velocity.
On the other hand, it has been found that double-lined binaries (SB2) are more common among HgMn than among normal stars (Wolff & Preston 1978; Gerbaldi et al. 1985). Wolff & Preston (1978) suggested that this excess of SB2, and of close binaries in general, could be explained in term of a strong trend toward synchronism of rotation and revolution, and therefore to low rotational velocities. Slow rotational velocity is one of the most important constraints needed to allow the formation of chemical peculiarity in the stellar atmosphere.
The knowledge of the orbital parameters is of great importance to study the global properties of this group of stars. To deal with homogeneous data sets, we selected all the SB2 systems reported in the catalogue of HgMn stars compiled by Schneider (1981) (23 out of 127 stars). In this paper we present and discuss the results for the first seven orbits, including the SB2 system HD 191110 recently published in Catanzaro et al. (2003).
Radial velocities for the program stars listed in Table A.1 have been derived from
spectra obtained with two different instruments:
To calculate radial velocities for each stars, we first identified all unblended spectral lines using the Kurucz (1993) lines list. For each of these lines we measured the central wavelength with a Gaussian fit of the profile and we computed the radial velocity using the classical Doppler shift formula. Finally we averaged the values obtained and computed the rms. In one spectrum of HD 173524 taken at CFHT (HJD = 2449968.7765) the stars were close to conjunction, so the previous procedure could not be applied. To determine the radial velocity we modeled the observed spectrum as described in Catanzaro et al. (2003). Photospheric parameters have been adopted from Adelman et al. (1998).
Systematic errors have been estimated in Catanzaro et al. (2003),
2.50
1.37 km s-1. At the beginning of 2003, because of
the change in the resolving power of the spectrograph, we needed to re-evaluate the errors.
We observed stars with constant and well known radial velocity: HD 107328, HD 136202,
HD 161096 and HD 204867 extracted from the list of CORAVEL standard stars
(Udry et al. 1999). We then computed the differences between radial velocity values
from the literature and our measurements. Since those differences are smaller than 1 km s-1,
we neglected them and the measured radial velocities of program stars observed in 2003
were not corrected.
The program stars with their measured radial velocities,
and heliocentric Julian Date
are listed in Table A.1.
The radial velocities for a spectroscopic binary system are:
![]() |
(2) |
![]() |
Figure 1: Radial velocity curves of HD 358 A and B. Solid line represents the orbital solution for the primary component, while the dashed curve represents the solution of the secondary orbit. Symbols used to represent the data for HD 358 A are: stars (Abt & Snowden 1973), filled squares (Aikman 1976), filled hexagons (Tomkin et al. 1995), filled triangles (Ryabchikova et al. 1999), asterisks (Adelman et al. 2002) and filled circles (our new measurements). To represent the data for HD 358 B we used the following open symbols: triangles (Ryabchikova et al. 1999) and hexagons (Tomkin et al. 1995). |
| Open with DEXTER | |
For each star, Table 3 reports the determined orbital elements,
the mass and the radii ratios between components. The radii ratio for each SB2
system was evaluated using the formula given by Giménez & Zamorano (1985)
HD 358 is a bright well-known binary with a HgMn primary, whose binarity was
discovered by Ludendorff (1907). As stated by Tomkin et al.
(1995), the brightness of the secondary is about two magnitudes fainter
than the primary. These authors and, later on, Ryabchikova et al. (1999)
using high S/N spectra (300 or more) and spectral resolution of R
40 000,
successfully measured the secondary's radial velocities and calculated the spectroscopic
orbital elements of this SB2 system. Recently Adelman et al. (2002),
performed an extensive study of the spectral variability of the Hg II
3984 Å line, finding the first identified spectrum variation in any HgMn star not
due to the orbital motion of the companion. Due to the low S/N of our data
we did not attempt to identify the secondary's spectral lines and we proceeded in our
analysis using literature data for HD 358 B.
The least-squares fits of Abt & Snowden (1973), Aikman (1976),
Tomkin et al. (1995), Ryabchikova et al. (1999), Adelman et al.
(2002) and our radial velocities using Eq. (1) give the
orbital period of 96.708
0.003 days (Fig. 1) that confirms
the one found by Ryabchikova et al. (1999).
The visual orbit of
And has been determined by Pan et al. (1992).
These authors using the Mark III Stellar Interferometre, Mount Wilson,
found an inclination angle of i = 105.66
0.22 degrees. Combining this value with
the
derived from our orbital solution, we got for the mass of
each components MA = 3.5
1.0
and MB = 1.8
0.5
.
These masses are in agreement, within the errors, with those inferred by
Ryabchikova et al. (1999).
![]() |
Figure 2: Radial velocity curves of HD 32694 A and B. Solid line represents the orbital solution for the primary component, while the dashed curve represent the solution for the secondary orbit. Symbols used to represent the data are: squares for Young (1976), triangles for Yushchenko et al. (1999, 2001) and circles for our data. Filled symbols refer to primary component and open symbols to the secondary. |
| Open with DEXTER | |
![]() |
Figure 3: Radial velocity curves of HD 33647 A and B. Solid line represents the orbital solution for the primary component, while the dashed curve represent the solution for the secondary orbit. Filled symbols are used to represent the data for HD 33647 A: squares (Hube 1970), triangles (Wolff & Preston 1978) and circles (our data). Open symbols are used to represent the data for HD 33647 B: triangles (Wolff & Preston 1978) and circles (our data). |
| Open with DEXTER | |
HD 32964 is a well-known SB2 system since the first orbit computed by Frost & Struve
(1924) in the early 1920 s. Young (1976) used Frost & Struve (1924)
radial velocities in combination with his own data to improve the determination of the
period. He found a new orbit (
= 5.522731 days) and pointed out that the
chemical composition of the components is different. Seven new data sets were collected
by Yushchenko et al. (1999, 2001), and confirmed the orbital solution
computed by Young (1976).
From our spectra, taken at OAC, we collected 10 new radial velocities for both components. Our solution is in agreement with previous works. The radial velocity curves are shown in Fig. 2 and our orbital parameters are reported in Table 3.
The first radial velocity measurements of this spectroscopic binary are by Hube (1970). Wolff & Preston (1978) measured the velocities of this SB2 system and they found a possible period of 25.01 days. A preliminary solution gave an eccentric orbit with e = 0.50.
Combining our new data with the previous, we obtained the orbital period of 24.3174 days, which is shorter than the one calculated by Wolff & Preston (1978). Our new orbital elements are reported in Table 3 and velocity curves for this system are shown in Fig. 3.
![]() |
Figure 4: Bottom panel: radial velocity curves of HD 173524 A and B. The primary component is represented by solid line and filled symbols, the secondary's by dashed line and open symbols. In particular, literature data are represented by: crosses and diagonals for respectively A and B velocities measured by Conti (1970), pentagons (Abt & Snowden 1973), triangles (Aikman 1976), squares (Stickland & Weatherby 1984) and rhombus (Adelman 1998). Our measurements are represented by circles (OAC spectra) and hexagons (CFHT spectra). Upper panel: velocity curve for the motion of the center of mass of HD 173524 A and B due to the presence of the third companion. Open squares are data taken from various literature studies and reported in Table 1, filled circles represent our observations. |
| Open with DEXTER | |
The first orbit for this double-lined system are by Petrie (1935). His data
together with those of Conti (1970) and
Abt & Snowden (1973) were used by Aikman (1976) to refine the
period. He found
9.81 days and a mass ratio of 1.18. Afterwards, other
observation were carried out by Stickland & Weatherby (1984) and
by Adelman et al. (1998).
We used all existing data plus the new velocities measured from our spectra (CFHT and OAC) to refine the period and the orbital solution. The bottom panel of Fig. 4 shows radial velocity curves for HD 173524 A and B; the orbital parameters are reported in Table 3.
The dispersion of velocities observed around points at the same phase has been
interpreted by Adelman et al. (1998) as evidence for the presence of a third
low-mass companion which revolves around the common center of mass with the spectroscopic
binary. These authors gave an estimation of the possible orbital period of the third
star (
= 33.6
2.0 years) and the semi-amplitude of the motion of
the center of mass of the SB2 system (K = 2.6
0.4 km s-1).
Following Adelman et al. (1998), we divided all the data we have at our disposal
in thirteen groups by corresponding epochs, then we computed the
-velocity values
fixing the other orbital elements equal to those we determined and reported in Table 3.
Using these points (see Table 1), we found the velocity curve reported
in the upper panel of Fig. 4 as a possible orbital solution for the center
of mass of the SB2. These velocities have been phased with a period of P = 36
3 years,
an eccentricity
0.13 and
146o. The standard deviation of
the fit for the third component is 0.6 km s-1. The semi-amplitude derived from our solution is
in agreement with the value given by Adelman et al. (1998).
HD 174933 is another well-studied SB2 system discovered by Meyer (1926a). Meyer found a period of 6.3624 days (Meyer 1926b), but he did not observe the secondary's lines in his spectra. A few radial-velocity measurements for the secondary have been reported by Conti (1970) and Seligman (1970). The first orbital solution for both components of this system has been computed by Ryabchikova et al. (1996), combining their own data with those published by the previous authors.
Table 1:
Heliocentric Julian Date and measured
-velocity for
HD 173524, n represents the number of points used. Last two rows represent our data.
Our spectra taken at the OAC do not have a high enough signal-to-noise ratio to show the secondary's lines. Thus, as we did in the case of HD 358, we measured radial velocities only for the primary component and we used literature data for HD 174933 B. For the primary components we considered also data published by Abt & Snowden (1973), Aikman (1976) and Stickland & Weatherby (1984). Derived orbital parameters are reported in Table 3, while the curves are displayed in Fig. 5.
![]() |
Figure 5: Radial velocity curves of HD 174933 A and B. Solid line represents the orbital solution for the primary component, while the dashed curve represent the solution for the secondary orbit. Symbols used to represent the data for HD 174933 A are: filled squares for Seligman (1970) data, filled hexagons for Conti (1970) data, open stars for Abt & Snowden (1973) data, asterisks for Aikman (1976) values, diagonals from Stickland & Weatherby (1984), filled triangles from Ryabchikova et al. (1996) observations and filled circle for our measurements. Open symbols have been used to represent data for the secondary's, namely: squares for Seligman (1970), hexagons for Conti (1970) and triangles for Ryabchikova et al. (1996). |
| Open with DEXTER | |
HD 191110 is a SB2 system with both components being CP stars, precisely the primary HgMn and the secondary Hg. This system has been recently studied by Catanzaro et al. (2003). The authors determined the orbital parameters (reported here again in Table 3) and they also performed a detailed abundance analysis. The rotational period is not reported in the literature but, since this system has null eccentricity, it could be considered synchronous.
![]() |
Figure 6:
Composite H |
| Open with DEXTER | |
In synchronous binary systems the components may be generally considered to rotate around axes perpendicular to their orbital planes and therefore the orbital inclination should correspond to the angle between the line of sight and the rotation axis of both component (Nariai 1971; Kitamura 1980).
According to these hypotheses we attempted to estimate the inclination of the orbital plane
with respect to the line of sight. If the SB2 system is
synchronous, the orbital period is equal to the rotational period. We then could write the
following equation for the equatorial rotational velocity:
HD 216494 is one of a number of HgMn stars discovered in a survey of late B-type stars by Wolff & Wolff (1974). Orbital parameters have been determined by Wolff (1974) who found a period of 3.42981 days and a mass ratio of 1.21. Other data have been collected by Aikman (1976), who did not attempt to refine the previous solution.
In an attempt to refine the orbital solution for the SB2, we observed this system for six nights at OAC and we also used four spectra downloaded from CFHT archive. Combining the measured radial velocities with those from the literature we performed a least-squares fit to Eq. (1). Results are reported in Table 3 and velocity curves are shown in Fig. 7 (bottom panel). We obtained a slightly shorter period than the result published by Wolff (1974) and a mass ratio equal, within the experimental error, to her determination.
Table 2:
Heliocentric Julian Date and measured
-velocity for
HD 216494, n represents the number of points used.
Table 3:
Orbital parameters and their errors (quantities in parentheses). Errors are
defined as the variation in the parameter which increases the
of a unit.
In the last column we report the standard deviation of the fit for both components.
![]() |
Figure 7:
Bottom panel: Radial velocity curves of HD 216494 A and B. The primary component is
represented with solid line and filled symbols, the secondary's with dashed line and
open symbols. In particular, triangles are data from Wolff (1974),
squares are data from Aikman (1976), hexagons and circle are our
measurements (CFHT and OAC, respectively).
Upper panel: |
| Open with DEXTER | |
HD 216494 is a triple system composed of this SB2 system plus a third companion orbiting around the common center of mass in 19.25 years (Tokovinin 1993). Hubrig & Mathys (1994) observed this system at high resolving power (R = 100 000). They suggested that the spectral type of the third component could be similar or slightly earlier than that of the HgMn primary. No more detail could be provided from their spectra.
We identified the H
of HD 216494 C in the spectra we obtained at the
Catania Astrophysical Observatory. Figure 6 shows the composite H
observed at the JD = 2 452 567.4115, when the A and B components were at the orbital phase 0.42.
Measured heliocentric radial velocity for the C component was 3.25 km s-1. We divided all the
data in four groups and, following the same procedure used for HD 173524, we
computed the
-velocities. These data are been reported in Table 2 and
displayed in the upper panel of Fig. 7, where they have been phased with the orbital
period of 19.25 years given by Tokovinin (1993). The small number of points did not
allow us to search for an orbital solution, but a variation in the radial velocity of the center
of mass of the SB2 system is clearly visible.
Manfroid & Mathys (1985) carried out photometric observations of this star. According to their study three rotational periods are possible (0.724/1.27/3.40 days), but given the low projected equatorial velocity observed, they concluded that the most probable period is 3.40 days. Thus, since this value is very close to the orbital period reported in Table 3, we conclude that HD 216494 could be synchronous.
Applying the same procedure as in the case of HD 191110 we attempted to estimate
the orbital inclination. Abt et al. (2002) estimated a rotational velocity for
the primary component
= 20 km s-1. This value is not compatible with
line profiles observed in the high resolution spectra (R = 110 000) downloaded from the
CFHT archive. From these spectra we estimated an upper limit for the velocities for each
components equal to 3 km s-1.
Using this value, from Eq. (4) it follows directly that
.
This inclination angle is not compatible with the high radial
velocity measured for each components. Moreover, this system has been studied by
Guthrie (1986), who found an inclination of the orbit equal to 54
,
a value more compatible with the measured velocities.
Thus, the initial hypothesis of a rotational axis perpendicular
to the orbital plane could be wrong. We can conclude that we are looking at the stars
almost pole-on, while the true orbit is tilted of 54
with respect to the line of sight.
We have carried out spectroscopic observations and determined the orbital elements of six binary systems whose brightest component is a chemically peculiar star of the HgMn sub-group. The results are summarized in Table 3.
In this study we first attempted to compute a possible orbital solution for the motion of the center of mass of HD 173524 A and B due to the presence of the third companion discovered by Adelman et al. (1998). The orbit for HD 33647 has been revised while for the other program stars orbital parameters have been slightly refined.
For HD 358, we computed the orbital solution and the masses of each component
as well. The latter values have been inferred combining
derived by us
with the inclination angle found by Pan et al. (1992). Our results agree
with the ones by Ryabchikova et al. (1999).
![]() |
Figure 8: The eccentricity as a function of the orbital period for our sample of SB2 system. Open circle refer to the data published by Guthrie (1986). |
| Open with DEXTER | |
Two SB2 systems of our program stars seem to be synchronous, HD 191110 and
HD 216494. Under the hypothesis the rotational axis is perpendicular to the orbital
plane and the assumption of a stellar radius equal to 3
,
we attempted to estimate
the inclination of their orbits. For HD 191110, we concluded the inclination angle
varies in the range 21-41 degrees, in agreement with the result of Guthrie (1986).
In the case of HD 216494, previous initial hypothesis fail and we concluded that the
components are almost pole-on while the orbital plane is tilted by 54
(Guthrie
1986) with respect to the line of sight.
HD 216494 is a triple system composed of the SB2 discussed above and a third
companion revolving around the common center of gravity in 19.25 years (Tokovinin
1993). We show that the variability of the center of mass of the SB2 is compatible
with the period found by Tokovinin (1993) for the third companion. Moreover, our
spectra, centered on the H
region, clearly show the presence of the third component.
An other important aspect to address is the membership of HgMn stars in multiple
systems. In this study we studied two of such systems (HD 173524 and
HD 216494). Another two stars are multiple systems: HD 32964 and HD 33647
(Tokovinin 1997). Thus, 4 out of 7 stars are in multiple systems, and if this
percentage is representative of the sample, we could conclude that
57% of
HgMn in SB2 systems have more than one companion.
In a review concerning the binarity among Bp-Ap stars, Gerbaldi et al. (1985)
pointed out that the distribution of eccentricity for the HgMn sample can be considered as
being the same as in any other sample of normal stars. The only difference seems to be a
slight bimodal distribution for 0.4 < e < 0.5, which is less in the normal
star sample. Our results reported in Fig. 8, together with those by Guthrie
(1986), confirm this behavior. This diagram shows that circular orbits (e = 0)
do exist for a wide distribution of orbital periods, up to P
10 days.
Due to the small number of systems observed to date, we do not apply other statistical considerations.
Acknowledgements
This research has made use of the SIMBAD database, operated at CDS, Strasburg, France. We wish to thank the referee, Dr. T. Ryabchikova, for her useful comments that helped us to improve the quality of our manuscript.
Table A.1: Heliocentric Julian Date and measured radial velocities for program stars. n represents the number of spectral lines used for the velocities determination.