A&A 433, 659-669 (2005)
DOI: 10.1051/0004-6361:20040396

Metallicities of the $\beta $ Cephei stars from low-resolution ultraviolet spectra[*]

E. Niemczura - J. Daszynska-Daszkiewicz

Astronomical Institute, Wroc\law University, ul. Kopernika 11, 51-622 Wroc\law, Poland

Received 5 March 2004 / Accepted 22 November 2004

Abstract
We derive basic stellar parameters (angular diameters, effective temperatures, metallicities) and interstellar reddening for all $\beta $ Cephei stars observed during the IUE satellite mission, including those belonging to three open clusters. The parameters are derived by means of an algorithmic procedure of fitting theoretical flux distributions to the low-resolution IUE spectra and ground-based spectrophotometric observations. Since the metallicity has a special importance for pulsating B-type stars, we focus our attention in particular on this parameter.

Key words: stars: early-type - stars: abundances - stars: variables: general

1 Introduction

$\beta $ Cephei variables are a well-known group of early B-type pulsating stars. Their pulsations are driven by the classical $\kappa$-mechanism, operating in the layer of the metal opacity bump at $T\approx 2\times 10^5~{\rm K}$ caused by the large number of absorption lines of the iron-group elements. The explanation of oscillations in these stars was possible only after the publication of the new opacity tables by Iglesias et al. (1992), and was given by Cox et al. (1992), Moskalik & Dziembowski (1992), Kiriakidis et al. (1992) and Dziembowski & Pamyatnykh (1993). As was shown by Pamyatnykh (1999), the extent of the theoretical instability domain in the H-R diagram is very sensitive to the heavy element abundance, Z. The instability domain decreases with decreasing value of Z, and for $Z\approx 0.01$ there are no unstable modes in the observed $\beta $ Cephei domain. Therefore, the determination of metallicities is the first step in testing the theory of pulsational instability in massive Main Sequence stars. It would be important to fix the minimum observed value of [m/H] for which the oscillations are still present in these stars. Most $\beta $ Cep pulsators belong to our Galaxy, and they are supposed to have a metal abundance close to the solar value. Recently, Pigulski & Ko\laczkowski (2002) discovered the first extragalactic $\beta $ Cep stars in the LMC, where the average metallicity is much below the solar value. All LMC $\beta $ Cep stars are in or close to young associations and they probably have higher metal abundance than average. Moreover, the question "which pulsational parameters are correlated with the metal abundance?'' is still open.

In this paper we analyze the IUE (International Ultraviolet Explorer) data combined with ground-based spectrophotometric observations of $\beta $ Cephei stars. The ultraviolet (UV) part of the spectra of main-sequence B-type stars is very rich in lines of the iron-group elements. Because the greatest amount of energy for these objects is emitted in the spectral region between 1216 and 10 000 Å, even the low-resolution spectra can give credible values of stellar parameters such as the effective temperature, $T_{\rm eff}$, surface gravity, $\log g$, microturbulent velocity, vt, metallicity, [m/H], angular diameter, $\theta$, and the color excess, E(B-V). The usefulness of the low-resolution IUE spectra for this purpose has already been demonstrated by many authors. Fitzpatrick & Massa (1999) obtained the mean stellar parameters for fourteen main sequence stars in the spectral range O9.5-A1. Morossi et al. (2002) used IUE spectra, combined with the visual spectrophotometric observations, to derive [m/H], $T_{\rm eff}$ and $\theta$ for 53 objects in the solar neighborhood. In Niemczura et al. (2002) we presented the results for some objects from the paper of Code et al. (1976) and for the Magellanic Clouds B-type stars, obtained by the least-squares optimization algorithm. In a previous paper (Daszynska et al. 2002), using the above method, we determined the values of the metal abundance parameter in the atmospheres of 31 field $\beta $ Cephei stars and 16 stars belonging to three open clusters NGC 3293, NGC 4755 and NGC 6231. Here we use the same data as in Daszynska et al. (2002) but we derive the stellar parameters with an improved procedure which allows, in particular, for a more objective error estimation. This method was applied recently by Niemczura (2003) to determine the metallicities, effective temperatures, stellar diameters and interstellar reddenings for Slowly Pulsating B-type stars.

In Sect. 2, all stars we analyze are presented and spectrophotometric observations together with their calibration are shown. The procedure of fitting the theoretical spectra to the observations and the method of bootstrap resampling adopted for determining the errors of parameters are briefly described in Sect. 3. The results and interparametric correlations are presented in Sect. 4. In Sect. 5, we discuss the values of metallicities we obtained. Conclusions are given in Sect. 6.

2 The observational material

The observational material consists of IUE observations obtained with the large or small aperture. We connected the data from both the long-wavelength (LWP and LWR, 1950-3200 Å) and the short-wavelength (SWP, 1100-1950 Å) cameras. We used observations with high (0.1-0.3 Å) and low (6-7 Å) spectral resolutions. The data were processed by means of two reduction packages: IUE/NEWSIPS (New Spectral Image Processing System, Nichols & Linsky 1996; Gahart et al. 1997) and IUE/INES (IUE Newly Extracted Spectra, Rodríguez Pascual et al. 1999; Gonzáles-Riestra et al. 2001, and the references therein). A detailed comparison of the IUE low-dispersion spectra extracted by the INES and NEWSIPS procedures made by Schartel & Skillen (1998) showed an excellent agreement in most cases. The main differences were found for underexposed spectra and strong spectral lines. Niemczura (2003) carried out a comparison of the IUE spectra processed by NEWSIPS and INES packages for SPB and reference stars. A similar comparison for the hotter stars made in this work gives the same results. The systematic differences appear in strong lines and at the wavelengths where the LWP/LWR and SWP spectra overlap. Large differences also occur in the wavelength range 2000-2300 Å for highly reddened stars. Other discrepancies are probably not systematic. As found by Niemczura (2003), better agreement exists for high-resolution data.

All data points with a flag indicating some data-quality problems were excluded from further analysis. Because of the presence of interstellar Ly$\alpha$ absorption, the spectral region around Ly$\alpha$was omitted. The IUE observations for $\lambda > 3000$ Å were also excluded from the analysis because of the very low signal-to-noise ratio at these wavelengths. Thus, we used the IUE fluxes in the spectral range from 1300 to 3000 Å. Whenever more than one spectrum was available for a star, the observations were co-added using a cross-correlation technique. The spectra were re-binned to theoretical wavelength points.

The ultraviolet observations expressed in absolute units were supplemented by ground-based spectrophotometric measurements taken from the Breger Catalogue (Breger 1976), Pulkovo Spectrophotometric Catalogue (Alekseeva et al. 1997), and from catalogues of Glushneva et al. (1998), Glushneva et al. (1992) and Kharitonov et al. (1988). If no spectrophotometric data were available (see Tables 1 and 2), we used Johnson UBV and Strömgren uvby magnitudes. The magnitudes were converted into fluxes by means of the formula $F_\lambda = 10^{0.4(C_\lambda - m_\lambda)}$, with scaling factors $C_\lambda$ taken from Gray (1998) and Bessell (1979) for the uvby and ${\it UBV}$ filters, respectively. All photometric data used in this paper were taken from the GCPD database (General Catalogue of Photometric Data, Mermilliod et al. 1997).

All stars selected for the analysis are presented in the Tables 1 and 2, where their names, HD numbers (ID numbers for cluster stars), numbers of IUE spectra, parallaxes and references for the ground-based observations of field objects are shown. Many $\beta $ Cephei stars we analyze belong to binary or multiple systems. Some of them are wide visual binaries with separation greater than 20 arcsec and can be resolved by the IUE cameras. In this group we have eight stars: $\nu$ Eri, 12 Lac, $\beta $ CMa, $\xi^1$ CMa, $\epsilon$ Cen, $\alpha$ Lup, BU Cir and $\gamma$ Peg.

In some cases the parameters of the companion indicate that it does not influence the spectrum of the much hotter $\beta $ Cephei star. BW Vul is a binary system ( $P_{\rm orb} = 33.3\pm0.3$ y, Pigulski 1992) with the secondary component 6-10 mag fainter than the primary. EN Lac is an eclipsing binary with orbital period 12.0968 d (Pigulski & Jerzykiewicz 1988; Lehmann et al. 2001, and references therein). The secondary component is fainter than the primary by about 5 mag. $\nu$ Cen is a spectroscopic binary (SB) system with an orbital period equal to 2.622$\;\pm\;$0.018 d (Schrijvers & Telting 2002) and a separation of $5.588\pm0.061$ arcsec (Shatsky & Tokovinin 2002). Merezhin (1994) determined the masses of the components to be equal to 10.20$\;\pm\;$0.10 and 1.57$\;\pm\;$0.25 $M_\odot$. $\beta $ Cep is a triple system with a distant visual component (separation equal to 13.4 arcsec). Gezari et al. (1972) found a close component, separated from the primary by about 0.25 arcsec. By 1992 this value decreased to 0.03 arcsec. The orbital elements of this system were determined by Pigulski & Boratyn (1992) who derived a period of 91.6$\;\pm\;$3.7 y from the light-time effect.

In some cases we have not enough information about the secondary component. 27 CMa is a binary system resolved by speckle interferometry (see for example Hartkopf et al. 2000). The separation of the system is growing, the last determined value being equal to $0.133\pm0.008$ arcsec. Abt & Cardona (1984) found the orbital period of the system, $P_{\rm orb} \approx 40$ y. $\theta$ Oph has a close companion with a separation of $0.300\pm0.025$ arcsec (Shatsky & Tokovinin 2002). $\tau^1$ Lup is a visual binary with a separation of 158.2 arcsec and a difference in magnitude of 4.7, therefore the visual component does not influence the spectrum of the $\beta $ Cephei star. But Proust et al. (1981) found a third component with a separation from the primary of the order of 0.3 arcsec. We have no information about this object.

In a few cases the spectra in the multiple systems are strongly influenced by the other components. $\sigma$ Sco is a spectroscopic binary (SB2) with a period of 33 d, and a component of a quadruple system (Pigulski 1992; Chapellier et al. 1992). The separation of the third component is equal to 0.4 arcsec and is growing (Pigulski 1992, and references therein). The brightness of this star is about 2.2 mag lower than that of the SB2 system and its orbital period is longer than 100 y (Evans et al. 1986). The fourth object in this system is separated by about 20 arcsec from the primary. $\beta $ Cru has a distant visual component (separation of about 371.6 arcsec). In addition, it is a primary component of a binary system. The orbital parameters of this system were determined by Aerts et al. (1998), who derived a period of 1828.0$\;\pm\;$2.5 d. The orbital parameters allowed to estimate the masses of the components, to be equal to 16 and 10 $M_\odot$. The spectral type of the secondary is B2 V, and is consistent with the temperature of 22 000-23 000 K and $\log g = 4.0$ dex. $\alpha$ Vir has an orbital period of 4.0145 d and the primary component is the $\beta $ Cephei variable. The magnitude of the secondary is lower by about 2.0 mag. The SB1 system $\beta $ Cen has an orbital period of 1 y (Shobbrook & Robertson 1968), and in addition a close visual component with a separation of 1.3 arcsec and $\Delta m = 3.2$ mag. The separation in the SB1 system is 15.6$\;\pm\;$2 mas (Robertson et al. 1999). Ausseloos et al. (2002), on the basis of observations of lines in the visual part of the spectrum, estimated that both components are $\beta $ Cephei-type stars. $\kappa$ Sco is an SB2 system with orbital period of 195 d (Molenberghs et al. 1999). The orbital parameters of the system and parameters of the components were determined by Uytterhoeven et al. (2001). Effective temperatures amount to 22 000 and 18 000 K, and the masses are equal to 17 and 12 $M_\odot$, respectively. $\lambda$ Sco is also a spectroscopic binary. The separation of the system is equal to 5.595 arcsec (De Mey et al. 1997). Heynderickx et al. (1994) estimated effective temperatures (21 200 and 23 100 K), surface gravities (3.77 and 3.76 dex) and masses (10.81 and 10.50 $M_\odot$) of the components. From the analysis of Berghöfer et al. (2000) of IUE observations it follows that the second component is a white dwarf. This SB system also has a visual component with a separation of 1.3 arcsec and $\Delta m = 3.2$ mag.

3 The analysis

The method of analysis was presented in detail by Niemczura (2003). Here we give only brief description of the method. We analyzed low-resolution spectra by means of an algorithmic procedure of fitting the theoretical flux distributions, $f_\lambda ^{\rm t}$, to the observations, $f_\lambda ^{\rm o}$. The theoretical flux depends on several atmospheric parameters like $T_{\rm eff}$, [m/H], $\log g$, vt. Also the angular diameter of the star, $\theta$ and interstellar extinction described by E(B-V)influence the theoretical spectrum. We used theoretical fluxes calculated by Kurucz (1996) with the standard value of the microturbulent velocity, vt = 2 km s-1. The spectral resolution of the observational and theoretical data were adjusted before the analysis.

We adopted the mean interstellar reddening curve of Fitzpatrick (1999) for the majority of the analyzed stars. For stars with high reddening ( E(B-V)>0.10 mag) the mean curve is not a good approximation because of spatial variability of the extinction law. In these cases we computed five additional parameters specifying the shape of the UV extinction curve (Fitzpatrick 1999). Three parameters describe the Lorentzian-like bump at 2200 Å (its width, $\gamma$, position, x0 and strength, c3), one is the far-UV curvature term (c4) and the last one is the linear term (c2).

Two additional parameters, S1 and S2, were used to adjust the short- and long-wavelength parts of the IUE spectra to the visual flux level. We put S1 = 1 outside the IUE/SWP region, and S2 = 1 outside the wavelength region covered by the long-wavelength IUE observations. These parameters are included in the analysis for several reasons. First, for some stars we used the IUE observations made with a small aperture, and these data are not calibrated to the absolute level with high precision. Second, the flux level changes during the pulsational cycle. Third, the SWP and LWP (or LWR) observations we used were sometimes made some years apart. In such cases, errors in absolute levels caused by instrumental effects can be significant. Finally, different absolute calibrations of the UV and visual parts of the spectrum may result in a difference of levels between these two spectral ranges.

Thus, the vector of all parameters ${\vec p}$ determined simultaneously can be expressed as

\begin{displaymath}\vec p = [T_{\rm eff}, {\rm [m/H]}, \theta, E(B-V), S_1, S_2, [\gamma, x_0, c_2, c_3, c_4]].
\end{displaymath}

To derive the vector of parameters ${\vec p}$, we used the least-squares optimization algorithm (Bevington 1969; Press et al. 1992). For more details see Niemczura (2003).

The problem of simultaneous determination of all parameters during the best-fit procedure can be solved if the parameters produce detectable and different spectral signatures. Fitzpatrick & Massa (1999) illustrated how atmospheric parameters and E(B-V) affect the spectrum of a hot star and concluded that the spectral signatures produced by these parameters are different. One can therefore expect that they could be obtained simultaneously by a best-fit procedure. The same is not true for $\log g$, which cannot be obtained in many cases or is determined with a very large error.

In order to determine surface gravities in the same way for all stars, we made use of photometric calibrations and stellar evolutionary tracks. The value of $\log g$ was derived as the mean from the results from five methods, but gravities lower than 3.00 dex and greater than 4.50 dex were excluded from the mean. We used three methods that employ Strömgren photometry (Napiwotzki et al. 1993; Balona 1994; Dziembowski & Jerzykiewicz 1999) and a method which uses the Geneva photometry (Künzli et al. 1997). In the last method, we estimated the gravities using the formulae obtained from stellar evolutionary tracks. The models were computed for OPAL opacities with X = 0.7 and Z = 0.02without taking into account the effects of rotation and convective overshooting. To this aim we considered Main Sequence evolutionary tracks for masses from 7.0 to 16.0 $M_\odot$ with a step of 1 $M_\odot$. We derived the following relation:

\begin{displaymath}\log g = -12.5894 + 4.4810 \log T_{\rm {eff}} - 0.7870 \log L/L_\odot,
\end{displaymath} (1)

with a standard deviation equal to 0.01 dex. We found that the formulae for Z = 0.01 and Z = 0.02 give a differences in $\log g$ of the order of 0.01 dex. The luminosity, $\log L/L_\odot$, depends on effective temperature, stellar angular diameter and parallax. Here we used Hipparcos parallaxes (ESA 1997). During the best-fit procedure the current luminosity was corrected for the Lutz-Kelker bias (Lutz & Kelker 1973). The correction for the luminosity can be calculated if $\sigma_\pi/\pi < 0.175$. Consequently, we can determine $\log g$by means of this method only for stars with $\sigma_\pi/\pi < 0.175$ (see Table 1). For the $\beta $ Cephei stars in clusters, the distances d to these clusters were taken from the literature: for NGC 3293 $d = 2750\pm250$ pc (Baume et al. 2003), for NGC 4755 $d = 2100\pm200$ pc (Sanner et al. 2001) and for NGC 6231 $d = 1990\pm200$ pc (Baume et al. 1999). This method of determination of $\log g$ was included into the iteration process, and the values of $T_{\rm eff}$ and $\theta$ were calculated during the best-fit procedure. Systematic errors can be caused by such effects as the chemical composition (X,Z), convective core overshooting, rotation etc. For lower masses there may also be a problem with the treatment of stellar convection.

Table 3: The best-fit parameters for the $\beta $ Cephei stars obtained from the IUE/NEWSIPS low-resolution data.


  \begin{figure}
\par\includegraphics[width=16.7cm,clip]{0396fig1.ps} \end{figure} Figure 1: Dereddened stellar energy distribution of $\omega ^1$ Sco (dotted line) in comparison with the best-fit model (solid line). The original flux of $\omega ^1$ Sco is also shown (crosses). The IUE data were calibrated by means of the NEWSIPS reduction package. In the inset, there is a closer view of the short-wavelength part of the spectrum.
Open with DEXTER

Formal errors of best-fit parameters resulting from the least-squares method are definitely underestimated. The technique of bootstrap resampling is probably the most useful method for estimating of confidence levels for complex least-squares solutions (Press et al. 1992, Chap. 15.6; see also Maceroni & Rucinski 1997). This technique uses the input data set, D0, consisting of N points, to generate simulated data sets, Di, with the same number of points. The symbol i denotes the successive simulation. The number of all bootstrap simulations is equal to $N(\log N)^2$ (Babu & Singh 1983). In the generated data, the random fraction of original points (37% in our case) is replaced by the remaining original points by means of random resampling with replacement. The simulated data sets are analyzed in the same way as the original data. This procedure not only makes it possible to determine the parameters and to estimate their reliable uncertainties, but it can be also used to determine correlations between them.

4 Results

4.1 The field $\beta $ Cephei stars


  \begin{figure}
\par\includegraphics[width=16.4cm,clip]{0396fig2.ps} \end{figure} Figure 2: The positions of the field $\beta $ Cep stars ( left panel) and those from the three open clusters ( right panel) on the $\log g$ vs. $\log T_{\rm eff}$ diagram. The evolutionary tracks were calculated without (solid line) and with (dotted line) taking into account the effect of convective core overshooting. We assumed $d_{\rm ov} = 0.2$.
Open with DEXTER

The parameters obtained for the field $\beta $ Cephei stars from the IUE/NEWSIPS spectra analyzed by the bootstrap resampling method are presented in Table 3. For a few stars (19 Mon and $\delta$ Lup) the differences between the parameters obtained on the basis of the IUE data calibrated by the NEWSIPS and INES procedures exceed the error bars. In the case of $\delta$ Lup the main reason for the large discrepancy can be the low quality of the IUE data because there is only a single SWP spectrum and a single LWP/LWR spectrum. There is also only one LWP/LWR spectrum for 19 Mon.

In Fig. 1 we show a comparison of the dereddened spectrum of $\omega ^1$ Sco and the theoretical spectrum resulting from the bootstrap resampling method. The theoretical spectrum was obtained for the parameters listed in Table 3. The data are plotted as a dotted line, the theoretical spectrum as a solid line. As one can see, there is a good overall agreement between these two spectra. In the inset, a closer view of the SWP part of the IUE/NEWSIPS spectrum is shown. All features of the observed spectrum are present in the theoretical fit.

Most stars we analyzed have small interstellar reddening. For stars with E(B-V) > 0.10 mag, the five additional parameters describing the extinction curve were determined (see Sect. 3) and are given in Table 4. For a comparison, parameters of the standard extinction curve of Fitzpatrick (1999) are also shown.

In Fig. 2 we show the $\beta $ Cephei stars we analyzed in the $\log g$ - $\log T_{\rm eff}$ diagram. Evolutionary tracks for M = 7, 9, 12, 16 and 20 $M_\odot$ are also plotted. The stellar models were computed using a standard evolutionary code written by Paczynski (1969, 1970) and improved by Sienkiewicz et al. (private communication). We used the most recent version of OPAL opacity data (Iglesias & Rogers 1996). The effect of rotation was not taken into account. An initial hydrogen abundance X = 0.70 and a metal abundance Z = 0.02 were assumed. Two types of evolutionary tracks are shown: without the effect of convective core overshooting (solid lines), and with the overshooting effect, assuming $d_{\rm ov} = 0.20$ $H_{\rm P}$ (dotted lines). As one can see, all $\beta $ Cephei stars are located inside the Main Sequence area.

4.2 $\beta $ Cephei stars in open clusters

The parameters obtained for $\beta $ Cephei stars belonging to the galactic open clusters NGC 3293, NGC 4755 and NGC 6231 are shown in Table 3. The analysis was performed for six members of NGC 3293, five members of NGC 4755, and five of NGC 6231. Only for two stars (NGC 3293-14, NGC 6231-238) the values of parameters obtained from the two sets of data, IUE/NEWSIPS and IUE/INES, do not agree within their errors. This is mainly caused by the poor quality of the IUE data.

The reddening is in all cases greater than 0.10 mag, so that there is a need to use a specified extinction curve. Because of the poor quality of the IUE spectra for the clusters stars, the extinction curves for NGC 3293 and NGC 6231 were taken from Massa & Fitzpatrick (1986), while for NGC 4755 the standard extinction curve of Fitzpatrick (1999) was used.

The study of $\beta $ Cephei stars in open clusters is important because their location in the colour-magnitude diagram indicates their evolutionary state. There are eleven open clusters in which $\beta $ Cephei variables were discovered (a recent example: three $\beta $ Cep stars in NGC 6910, Ko\laczkowski et al. 2004). A still unexplained problem is why some young clusters contain those variables and others, like NGC 2362 (Balona & Laney 1995), do not. The other problem is the difference in location of the instability strips for various clusters of nearly the same age. According to Balona & Koen (1994), the location of the $\beta $ Cephei variables belonging to NGC 4755 on the HR diagram is shifted to lower effective temperatures compared to NGC 3293. However, we obtained opposite results: the location of $\beta $ Cephei stars in NGC 4755 is shifted to the higher effective temperatures as compared to NGC 3293. In the case of NGC 4755 we found lower values of [m/H] than in the case of NGC 3293. For the lower metallicity the ZAMS line is shifted to the left in the H-R diagram, thus our results are fully consistent.

4.3 Correlations between the parameters

The mean correlations, $\rho _{\rm mean}$, their standard deviations, $\sigma _\rho $, median values, $\rho _{\rm m}$, and ranges of correlation coefficients for the $\beta $ Cephei stars are presented in Table 5. Because of the subject of this work, the most interesting are the correlations between metallicity and other parameters obtained simultaneously. The values of $\rho _{\rm mean}$ and $\rho _{\rm m}$ for metallicity are small and amount to about 0.20. Also the standard deviations are small, which means that even for individual objects correlations are not large. Small values of $\rho _{\rm mean}$ indicate that the metallicity, [m/H], is not correlated with the other parameters and can be determined from the best-fit procedure with good precision. The correlations between the other parameters are larger. The effective temperature is strongly correlated with the interstellar reddening, $\rho_{\rm mean} = 0.73$. A lower value of $\rho _{\rm mean}$was obtained for effective temperature and stellar diameter (0.66). Somewhat lower values of $\rho _{\rm mean}$ and $\rho_m$ are found in the other cases, but the correlations can be strong for individual objects.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{0396fig3.ps} \end{figure} Figure 3: Differences in parameters, $\Delta p_i = p_i({\rm NEWSIPS}) - p_i({\rm INES})$, as a function of the effective temperature obtained from the IUE/NEWSIPS data.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{0396fig4.ps} \end{figure} Figure 4: The same as in Fig. 3, but as a function of the number of IUE spectra, n (IUE).
Open with DEXTER

4.4 Comparison between NEWSIPS and INES

In Fig. 3 we plot the differences between the parameters obtained from the IUE data calibrated with the NEWSIPS and INES procedures, $\Delta p_i = p_i({\rm NEWSIPS}) - p_i({\rm INES})$, as a function of $\log T_{\rm eff}$. In the four panels we can see these differences for $\log T_{\rm eff}$, [m/H], $\theta$, and E(B-V). For the small range of effective temperature shown in the Fig. 3 there is no correlation between $\Delta p_i$ and $\log T_{\rm eff}$.

Small mean differences, $\langle\vert\Delta {\vec p}\vert\rangle$, were found for all parameters, that is $0.014\pm0.001$ for $\log T_{\rm eff}$ ( $\sigma = 0.023$ dex), $0.003\pm0.0001$ ( $\sigma = 0.003$ mas) for $\theta$, $0.046\pm0.001$ ( $\sigma = 0.050$ dex) for [m/H] and $0.010\pm0.001$ ( $\sigma = 0.012$ mag) for E(B-V). For all this mean values, the standard deviations are quite large. This means that the differences for individual objects can be significant.

In Fig. 4 we show the same differences, $\Delta p_i$, as a function of the number of available IUE spectra. As one could expect, the largest differences are obtained for stars with at most two spectra available from each IUE camera. When the number of the IUE spectra grows, the scatter of $\Delta p_i$ decreases.

5 Discussion of the parameter [m/H]

5.1 General properties


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{0396fig5.ps} \end{figure} Figure 5: The distribution of the $\beta $ Cephei stars with the [m/H] parameter. The hatched histogram shows the results obtained from the IUE/NEWSIPS data, while the gray one presents results obtained from the IUE/INES data.
Open with DEXTER

The metallicities of the $\beta $ Cephei stars determined from the IUE/NEWSIPS spectra range from $ -0.47 \pm 0.18 $ dex for 27 CMa to $+0.21\pm0.09$ dex for HN Aqr. The extreme values determined from the IUE/INES data were obtained for $\tau^1$ Lup and, again, HN Aqr; they are equal to $-0.42\pm0.17$ dex and $+0.20~\pm~0.10$ dex, respectively. The distribution of the [m/H] parameter resulting from NEWSIPS and INES spectra is shown in Fig. 5. These two samples follow normal distributions within the error limits. Applying the one-way Anova statistical test we obtained that the two means are not significantly different up to the significance level of $\alpha = 0.83$. Thus we conclude that there is no difference between the values of [m/H] found from the IUE data processed by both reduction packages. The mean values of the metal abundance parameter for the $\beta $ Cephei stars are equal to $-0.15\pm0.03$ dex and $-0.13\pm0.03$ dex for IUE/NEWSIPS and IUE/INES, respectively. The mean metallicities for the field $\beta $ Cephei stars are equal to $-0.14\pm0.03$ dex (IUE/NEWSIPS) and $-0.13\pm0.03$ dex (IUE/INES). The mean values of the metal abundance parameter for the cluster stars, obtained from the analysis of IUE/NEWSIPS spectra are equal to $+0.05\pm0.06$ dex for NGC 3293, $-0.43\pm0.05$ dex for NGC 4755 and $-0.12\pm0.09$ dex for NGC 6231. From the IUE/INES spectra we obtained nearly the same values: $+0.08\pm0.07$ dex for NGC 3293, $-0.40\pm0.05$ dex for NGC 4755 and $ -0.05 \pm 0.10 $ dex for NGC 6231. These determinations are consistent with our previous results (Daszynska et al. 2002), obtained with a less accurate method.


  \begin{figure}
\par\includegraphics[width=7.9cm,clip]{0396fig6.ps} \end{figure} Figure 6: The distribution of the metallicity parameter of the $\beta $ Cephei stars from the IUE/NEWSIPS data as a function of the galactocentric distance (in kpc). The dashed line indicates the mean value of [m/H] obtained for all $\beta $ Cephei stars. Symbols are the same as in Fig. 2.
Open with DEXTER

Metal abundances of hot stars in the solar vicinity lower by about 0.20 dex than the solar value were reported by Gies & Lambert (1992, hereafter GL92), Cunha & Lambert (1994), Kilian (1992), Kilian et al. (1994), Daflon et al. (1999), Gummersbach et al. (1998), Niemczura (2003), and others. These results were obtained from both high- and low-resolution data. All quoted authors used Kurucz's ATLAS models of stellar atmospheres. Metal abundances obtained in the present paper from the low-resolution IUE observations are in agreement with these results.

In Fig. 6 we present the dependence of [m/H] on the galactocentric distance. The distance for the field stars were calculated by using the Hipparcos parallaxes (ESA 1997). For the stars from open clusters we assumed the same values of the distance from the Sun as in Sect. 3. The galactocentric distance of the Sun, $R_{\rm G,\odot} = 7.9$ kpc was adopted from McNamara et al. (2000). As can be seen from the figure, the field stars are spread over a small range of about 0.5 kpc in galactocentric distance, so that no clear indication of a metallicity gradient is apparent. However, we would like to point out that some intrinsic scatter of the metallicity is certainly present.

5.2 Comparison with the literature

For thirteen $\beta $ Cephei stars in our sample, the metallicity has been determined by other authors.

One of the well-studied $\beta $ Cephei stars is $\gamma$ Peg. The majority of determinations was obtained from the visual part of the spectrum. Snijders (1969) determined [Fe/H] to be equal to +0.48 dex, while Peters (1976) found a much lower iron abundance, equal to +0.04 dex. GL92 obtained [FeII/H] = +0.29 dex and Pintado & Adelman (1993) found [Fe/H] = +0.04 dex. From the high-resolution ultraviolet IUE spectra, Proffitt & Quigley (2001, hereafter PQ2001) determined [Fe/H] = -0.15 dex. The values of [m/H] obtained by us, $ -0.04 \pm 0.08 $ dex (IUE/NEWSIPS) and $+0.02\pm0.10$ dex (IUE/INES), are then consistent with the determinations of Peters (1976) and Pintado & Adelman (1993).

The differences between the values of metal abundances for 16 Lac determined by different authors are significant. GL92 obtained the iron abundance equal to -0.23 dex from the visual part of the spectrum, while from the UV spectrum PQ2001 obtained [Fe/H] = -0.09 dex, and Venn et al. (2002) found [m/H] = $+0.16\pm0.18$ dex. Recently, Thoul et al. (2003) determined the abundances of the elements from the visual part of the spectrum. According to these authors, the iron abundance amounts to $-0.08\pm0.09$ dex. Our values are consistent with the last results to within the errors of determination: [m/H] = $ -0.13 \pm 0.13 $ dex and $-0.15\pm0.14$ dex.

For $\delta$ Cet GL92 determined the iron abundance to be equal to -0.24 dex from the analysis of the visual part of the spectrum. Different values were obtained from the ultraviolet spectra. From the IUE spectra PQ2001 determined [Fe/H] to be equal to -0.15 dex. Venn et al. (2002), from HST/STIS data, found [m/H]  $ = -0.15\pm0.15$ dex. All these determinations are in agreement with the present work. For 15 CMa and $\beta $ Cep the metal abundances were determined by the same authors. For 15 CMa GL92 found [Fe/H] = -0.47 dex from the visual part, while from the analysis of the UV data, PQ2001 and Venn et al. (2002) obtained -0.13 dex and $+0.10\pm0.10$ dex, respectively. The values we found are close to the latter determinations. The metal abundances derived for $\beta $ Cep range from -0.31 dex (PQ2001) to $-0.16\pm0.23$ dex (Venn et al. 2002). Both these values were determined from the UV spectra. From the lines in the visual part of the spectrum, GL92 obtained [Fe/H] = -0.23 dex. The values of [m/H] obtained by us are equal to $ -0.07 \pm 0.10 $ dex and $-0.01\pm0.09$ dex, respectively.

For the rest of the $\beta $ Cephei stars there are only two or one determination of the metal abundances in the literature. Recently, Stankov et al. (2003) announced a determination of the abundances of elements for BW Vul, derived from high-resolution spectra in the blue. According to these authors, the iron abundance, [Fe/H], amounts to +0.02 dex. This is close to the value obtained by us, [m/H]  $ = +0.07\pm0.12$ dex.

For HN Aqr, the first abundance determination was done by Ryans et al. (1996). From visual part of the spectrum, the authors obtained [Fe/H] = -0.84 dex. Ramspeck et al. (2001), from the same spectral range, determined a much higher value of the iron abundance, [Fe/H] = -0.12 dex. In this work, [m/H]  $ = +0.21\pm0.09$ dex (IUE/NEWSIPS) and $+0.20\pm0.10$ dex (IUE/INES) were obtained. The metallicities obtained in the present paper for $\beta $ Cru are consistent with the values given by Kilian (1994), [Fe/H] = -0.08 dex (from the visual part of the spectrum) and PQ2001, [Fe/H] = -0.25 dex. For BU Cir, the values determined in this work are lower than the iron abundance, [Fe/H] = -0.12 dex, given by Kilian (1994).

Our metallicities of $\nu$ Eri are close to the [Fe/H] = -0.08 dex given by GL92, but the iron abundance determined by PQ2001 is lower, namely -0.27 dex. For $\xi^1$ CMa we obtained metallicities close to the iron abundance of -0.30 dex determined by PQ2001. The values obtained in this paper are also consistent within the error limits with [Fe/H] = -0.18 dex obtained by GL92. For $\beta $ CMa we obtained [m/H]  $ = +0.04\pm0.10$ dex and $+0.08\pm0.11$ dex from the IUE/NEWSIPS and IUE/INES data, respectively. However, GL92 determined [Fe/H] = -0.41 dex and PQ2001 found that [Fe/H] = -0.22 dex. In the case of 12 Lac, the values determined in this paper: $ -0.20 \pm 0.10 $ dex and $-0.22\pm0.10$ dex are higher than the iron abundance obtained by GL92, -0.41 dex. The iron abundance found by PQ2001 is equal to -0.08 dex.

5.3 Metallicity and pulsations

In Fig. 7 we show the instability domains of $\beta $ Cephei variables for three values of the metal abundance ( Z = 0.01, 0.015, 0.02) and the stars we analyzed. As before, the field $\beta $ Cep stars and those belonging to the three open clusters are plotted in separated panels. The majority of field $\beta $ Cephei stars is located within the instability region defined by Z = 0.015. Only for three stars are the values of surface gravity higher than 4.0 dex. Two stars, 27 CMa and $\omega ^1$ Sco, are located outside the instability regions. The metallicity parameters, [m/H], for both stars are lower than -0.40 dex. This is not a surprise as 27 CMa is a Be star. In the case of $\omega ^1$ Sco there is also the possibility that it is a $\zeta$ Oph-type variable.


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{0396fig7.ps} \end{figure} Figure 7: The location of the $\beta $ Cephei stars on the $\log g$ vs. $\log T_{\rm eff}$ diagrams. The instability regions for Z = 0.02, 0.015 and 0.01 were taken from Pamyatnykh (1999). In the upper panel the position of the field $\beta $ Cephei stars is shown, while in the lower one, the position of variables from the open clusters is presented.
Open with DEXTER

To check which pulsational parameter is sensitive to the metallicity in $\beta $ Cep stars, in Fig. 8 we plot the period of the dominant mode, P, radial velocity amplitude, 2K, light amplitude in the V filter, AV, and a ratio of photometric amplitudes, AU-B/AV, as a function of the [m/H] parameter. We note that the last quantity is independent of the inclination angle, i, and the intrinsic pulsational amplitude, $\varepsilon$. We have checked all passbands and colors available in the literature, and here we give representative examples. Panel a of Fig. 8 is for the pulsational periods, the values of which are marked as filled circles for the field $\beta $ Cep stars, and as diamonds, triangles and squares for $\beta $ Cephei stars in the clusters NGC 3292, NGC 4755 and NGC 6231, respectively. As one can see, there is no correlation between metallicity and mentioned pulsational parameters, as was already pointed out by Daszynska et al. (2002).

Also the projected rotational velocity, $v_{\rm e}\sin i$, is not correlated with the metallicity parameter.


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{0396fig8.ps} \end{figure} Figure 8: Dependence of pulsational parameters: P, 2K, AV and AU-B/AV on the metallicity parameter, [m/H].
Open with DEXTER

In our previous papers (Daszynska 2001; Daszynska et al. 2003) we stated that the mean metallicity parameter, [m/H], is higher in the case of multiperiodic $\beta $ Cep stars than in the case of monoperiodic variables. Here we confirm this result. The mean value of the [m/H] parameter for multiperiodic $\beta $ Cep stars amounts to $-0.07\pm0.03$ dex, while the mean metallicity for monoperiodic variables is equal to $-0.25\pm0.05$ dex. This result is in agreement with the nonadiabatic pulsational theory which predicts a larger number of unstable modes for higher values of the metal abundance.

The dichotomy of $\beta $ Cep stars with regard to metallicity is discussed in detail in a separate paper (Daszynska-Daszkiewicz & Niemczura 2005).

6 Conclusions and prospects

We derived basic stellar parameters such as the effective temperature, gravity, angular diameter, metallicity and interstellar reddening for all $\beta $ Cep stars monitored by the IUE satellite. Among them there were field $\beta $ Cep stars and variables belonging to three young open clusters: NGC 3293, NGC 4755 and NGC 6231. The mean values of the [m/H] parameter for the $\beta $ Cephei stars are equal to $-0.15\pm0.03$ dex and $-0.13\pm0.03$ dex for IUE/NEWSIPS and IUE/INES data, respectively. The mean metallicities for the field $\beta $ Cephei variables are equal to $-0.14\pm0.03$ dex (IUE/NEWSIPS) and $-0.13\pm0.03$ dex (IUE/INES), while the mean values of [m/H] for the cluster stars obtained from the analysis of IUE/NEWSIPS spectra are equal to $+0.05\pm0.06$ dex for NGC 3293, $-0.43\pm0.05$ dex for NGC 4755 and $-0.12\pm0.09$ dex for NGC 6231. From the IUE/INES spectra, we obtained nearly the same values: $+0.08\pm0.07$ dex for NGC 3293, $-0.40\pm0.05$ dex for NGC 4755 and $ -0.05 \pm 0.10 $ dex for NGC 6231. These values are consistent with our previous results (Daszynska et al. 2002), obtained with the less accurate method.

The values of the metallicity parameter we obtained are not correlated to any pulsational parameter or the projected rotational velocity. However, we have noticed that, on the whole, multiperiodic variables have higher values of metallicity than the monoperiodic ones. This is consistent with pulsation theory.

The abundance of metals, especially that of iron, in $\beta $ Cephei variables is one of the fundamental parameters. In the case of main-sequence B stars, metallicities [m/H] obtained from the UV spectra give mainly an estimation of [Fe/H]. The value of [Fe/H] in the Kurucz's models is larger than in the widely used in the pulsation theory OPAL and OP data. One can therefore expect that [Fe/H] determined using the ATLAS9 models should be lower in comparison with [Fe/H] obtained on the basis of the OPAL or OP data (see also Fitzpatrick & Massa 1999). To compare observational metallicities of the stars with the values resulting from the theory, we have to include the correction for [m/H] of about 0.12 dex. The metallicities of the $\beta $ Cep stars predicted by the theory of pulsations are not in contradiction with most values determined by us. Pamyatnykh (1999) showed that the $\beta $ Cephei instability strip vanishes for $Z\approx 0.01$, corresponding to [m/H]  $\approx -0.30$ dex, and to [m/H]  $\approx -0.42$ dex from Kurucz's fluxes. There are several stars with a metallicity parameter lower than this limit. However, our values give information about photospheric metal abundances. Some other effects that we did not take into account such as diffusion or element mixing may also be important. Also, we used a standard mixture of elements. Nevertheless, our results provide important information about the metallicity range for these pulsating stars. In particular, because there is a problem with the estimation of the [m/H] parameter for early type stars from photometry, making use of low-resolution ultraviolet spectra seems to be the best and fastest way to determine the metallicity of B-type pulsating variables.

However, for asteroseismological purposes, there is a need for a detailed analysis of chemical composition, because oscillation frequencies are very sensitive to the adopted mixture, and this is the aim of future work.

Acknowledgements
We thank Profs. W. Dziembowski and A. Pamyatnykh for permission to use the evolution code and non-adiabatic pulsation code.

References

 

  
Online Material

Table 1: The list of $\beta $ Cephei stars. The names, HD numbers, the Hipparcos parallaxes, $\pi $, numbers of the IUE images used, sources of the ground-based data (last column) are given.

Table 2: The list of $\beta $ Cephei stars belonging to open clusters. ID numbers, names and numbers of the used IUE images are given.

Table 4: The best-fit parameters specifying the shape of the UV extinction curve. In the first line parameters of the standard extinction curve (Fitzpatrick (1999) are showed.

Table 5: The mean correlations, $\rho _{\rm mean}$, their standard deviations, $\sigma _\rho $, median values, $\rho _{\rm m}$, and ranges of correlation coefficients between the determined parameters.



Copyright ESO 2005