A&A 457, 949-961 (2006)
R. E. Mennickent1 - L. Cidale2 - G. Pietrzynski1,3 - W. Gieren1 - B. Sabogal1
1 - Universidad de Concepción, Departamento de Física, Casilla 160-C, Concepción, Chile
2 - Facultad de Ciencias Astronómicas y Geofísicas, Univ. Nacional de La Plata, Paseo del Bosque S/N, 1900 CGA, La Plata, Argentina
3 - Warsaw University Observatory, Al. Ujazdowskie 4,00-478, Warsaw, Poland
Received 4 April 2006 / Accepted 27 May 2006
Aims. We investigate the nature of a sample of 17 long-term periodic variables in the direction of the Small Magellanic Cloud.
Methods. Based on new spectroscopic data, we determined spectral types, radial velocities, absolute magnitudes, and colors for these stars. We present a refined discussion of their OGLE light curves along with an analysis of their 2MASS photometry.
Results. Most stars turned out to be B-A giants members of the Small Magellanic Cloud. We find a new interacting eclipsing binary with a period of 184 days and two new early-type ellipsoidal variables. One of our objects is the ROSAT source RX J0058.2-7231. We analyzed 11-years of data for this Be X-ray binary finding that their photometric period varies by 4% quasiperiodically, on a time scale of 1200 days. We find evidence of multiple photometric periods in 2 Ae-type and 1 late-Be type stars. The case of OGLE00445466-7328029 is especially interesting, as this late-type Be star shows a beating phenomenon primarily caused by two closely-spaced frequencies, 0.05733 c/d (17.44 d) and 0.06347 c/d (15.76 d). Four other emission-line objects show stable long-term periodicities and probably correspond to Be-star binaries. Transient photometric periods are only found in 4 non-emission early-type stars.
Key words: stars: binaries: general - stars: binaries: eclipsing - stars: early type - stars: emission-line, Be - stars: circumstellar matter - stars: mass-loss
From an inspection of the OGLE-II photometric database (Zebrun et al. 2001), while searching for Be stars in the Small Magellanic Cloud (SMC), Mennickent et al. (2002, hereafter M02) found 78 variable stars characterized by periodic photometric variability. As noted by these authors, the nature of these variables is still unclear. They are redder than typical Be stars, usually show variable amplitudes, and subsequent research has shown that their periods clearly show a bimodal distribution (Mennickent et al. 2003) with a main broad peak between 20 and 130 days and a narrower secondary peak between 140 and 210 days. The same kind of type-3 stars were found in the Large Magellanic Cloud, and they also show a bimodal period distribution (Sabogal et al. 2005). As noted by M02, the photometric variability is very different to that observed in any known pulsating variable in the same color range, like SPBs or Cyg stars.
From these investigations several interesting questions arise, for instance: What kind of phenomenon is triggering the photometric activity in these stars? Why is a bimodal period distribution for type-3 variables observed in both galaxies? In order to contribute to answering these very fundamental questions we present here the first spectroscopic and detailed photometric study of a sub-sample of bright type-3 stars. We selected objects with V magnitudes in the range 13.4-15.5 mag. In Sect. 2, we give an overview of the photometric and spectroscopic observations, along with a brief description of the methods employed in the analysis. In Sect. 3 we present the results of our study. We discuss our results in Sect. 4, and give conclusions in Sect. 5.
|OGLE00445466-7328029||SMC3-2||14.82||137 26||7||B7-9 IIIe||M|
|OGLE00455414-7314043||SMC3-3||15.48||-2.66||174 23||7||F0 Ib-II||S|
|OGLE00463376-7312043||SMC3-4||14.06||126 15||1||F5 Ie + G0 I||E1|
|OGLE00475014-7313164||SMC3-5||15.48||-2.05||139 5||7||B7-8 III||T|
|OGLE00492141-7258449||SMC3-6||14.61||-0.84||191 23||13||A5 II||T|
|OGLE00502564-7258071||SMC3-7||14.01||127 15||10||F4 IV||T|
|OGLE00510018-7253039||SMC3-9||13.61||172 32||12||B1 II-IIIe||S|
|OGLE00535922-7235089||SMC3-12||14.40||-0.01||54 17||12||A3 III||E2|
|OGLE00552027-7237101||SMC3-14||14.88||141 16||2||A0 Ibe||M|
|OGLE00574525-7235321||SMC3-15||13.82||-0.23||-4.33||195 37||22||B1 II-III||T|
|OGLE00581258-7230485||SMC3-16||14.90||164 24||4||B0 IIIe||Be X-ray|
|OGLE01000078-7255229||SMC3-18||14.44||222 32||14||B8-9 IIIe||S|
|OGLE01045121-7246469||SMC3-20||13.42||-0.21||-4.32||179 32||16||B1 Ib-IIe||S|
|OGLE01045299-7159188||SMC3-21||14.22||-0.18||-4.36||143 24||15||B2 IIIe||S|
Spectroscopic observations were conducted at CTIO during October 8-12, 2003 (UT). We used the 1.5 m telescope with the Cassegrain Spectrograph and the Loral 1K detector. The first two nights, the grating #26 tilted at 16.2 degree and a slit width of 2 arcsecond yielded a spectral range of 3700-5470 Å, a dispersion of 1.3 arcsec/pix, and an effective resolution of 2 Å. The Loral CCD had a few dead pixels around wavelength 4051 Å; the artifacts produced by these pixels in the spectra were removed by interpolating the flux of nearby pixels. During the third and fourth nights, we used grating #36, with the same slit width, which yielded a dispersion of 0.77 Å/pix, a 1.9 Å effective spectral resolution, and a wavelength range of 5765-6694 Å. For this setup, the CCD dead zone was between 5948 and 5951 Å. As for the blue spectra, bad pixels were removed by interpolating the flux between nearby pixels. The He-Ar comparison spectra were obtained between several science exposures typically every half hour. These comparison spectra provided wavelength calibration functions with typical rms of 0.1 Å for grating #26, and 0.01 Å for grating #36 (6 km s-1 and 0.5 km s-1 at H and H, respectively). Observations of the standard star LTT 2415 were used to calibrate the spectra in absolute flux units. The rms of the sensitivity function was 0.05 mag for grating #26 and 0.02 mag for grating #36. Corrections for atmospheric extinction were done for every spectrum using average extinction coefficients available for CTIO. Spectroscopic reductions were carried out in the usual manner, using several IRAF routines like ccdproc and doslit. A summary of the spectroscopic observations is given in Table 1.
In this section we present the methods used for spectral classification, give details on our search for photometric periodicities, provide a brief discussion of every object and discuss the available infrared photometric data.
We measured radial velocities and equivalent widths of several lines in the spectra. The spectral classification of each star was derived by comparing it with the MK standard classification system, using the digital spectra available on-line (http://stellar.phys.appstate.edu/Standards/stdCTIO.html). In the case of early-type stars, we also determined the spectral classification by using the measurements of the equivalent widths of the strongest lines (Didelon 1982). We mainly used lines of He I, Si II, Si III, Fe II, and Mg II. In a few cases we used the absorption lines of H and H to confirm the classification, but in general we avoided them since they were usually filled by emissions. For many stars we were able to estimate absolute magnitudes (MV) and intrinsic colors, (B-V)0. For B-type stars we adopted the equivalent widths of some selected pairs of He I lines (Stock et al. 2002) and for A, F, and G-types we chose the metallic lines and the H line (Stock & Stock 1999). We could not determine fundamental parameters for some stars since the lower members of H lines displayed emission features. Our results are listed in Table 2. There N is the number of lines averaged for the radial velocity determination and a note about photometric variability is provided: M stands for multiperiodic, T for transient periodicity, E1 for eclipsing binary, E2 for ellipsoidal variable and S for single period variable.
We observe that the determinations of MV and (B-V)0 agree with the spectral types that were derived independently. The major uncertainties may be present in late-B, A, and early-F types due to the presence of weak metallic lines and emissions in H and H. We note that the five early B-type stars show similar continuum slopes (see Sect. 3.3); a ratio between 2 and 2.7, except the enigmatic star SMC3-9, that shows a remarkable down turn in flux downwards of about 4500 Å. We observe this fact as a confirmation of the similar spectral type and the presence of reddening. As a reference, black-bodies with temperatures of 20 000, 13 000, and 10 000 K have the aforementioned ratios equal to 3.8, 2.7, and 2, respectively.
Our sample consists of stars whose photometric periods were first determined by M02. As that was a seminal paper, focused on the discovery of new types of variables in the Small Magellanic Cloud, no further detailed information was given about the possible multiperiodicity of the stars. In this section we describe our method of searching for multiple frequencies in the OGLE time series and list the results. The photometric data discussed here is part of the OGLE (Zebrun et al. 2001) database, and it can be accessed at http://www.astrouw.edu.pl/~ogle. The OGLE dataset consists of I-band time series observed between HJD 2 450 466 and HJD 2 451 757 (3.5 years), typically with 271-306 observations per star. To test the stability of our solutions, we also checked if the frequencies we found were also present in the OGLE-III database (Pietrzynski, private communication). This database now consists of four additional years of observations, which have almost been reduced at the moment of the writing of this paper, but lacking reliable zero points.
During our investigation of frequencies in the time domain, we used a computer program especially dedicated to the statistical analysis of large astronomical time series containing time gaps, the Period04 program first developed by Martin Sperl (Lenz & Bregel 2005, http://www.astro.univie.ac.at/~dsn/dsn/Period04/).
First, we got the stronger frequency in the signal
by applying a fast Fourier transform to the data. The analysis was done between 0 c/d and
the Nyquist threshold. The main peak of the periodogram, i.e. those with the largest amplitude,
was selected as our first frequency, labeled F1.
A least-square sinus fit was then performed on the data with the function:
We calculated uncertainties in the fit parameters using Monte Carlo simulations. We constructed 100 trials of synthetic datasets, determined solutions for every one, and calculated the standard deviation of frequency, amplitude, and phase. The results of our analysis are summarized in Table 3, which gives the frequencies, amplitudes (A, in mmag), and phases, along with their uncertainties. The table also gives the noise level for every star and the standard deviation of the fit for every new frequency added, and a flag indicating whether the frequency passed the stability test.
In some cases, we obtained unrealistically high errors associated to the amplitudes of the sinusoidal solutions. We interpreted these results as unstable solutions, which could be due to formal, rather than physical, reasons. For these cases, we calculated uncertainties from the error matrix of the least-square calculation.
We cross-correlated our targets with stellar catalogues using the SIMBAD on-line query form. A number of matches were found that correspond to 7 previously catalogued emission-line objects (from H imaging surveys), one Wolf Rayet star (OGLE00504344-7327053 = AzV 81), one Be X-ray binary (OGLE00581258-7230485 = RX J0058.2-7231), and the non-emission object SMC41243 (Massey 2002). A brief description of the results for individual objects and a summary of the literature search is given below.
|Figure 1: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and the spectrum of SMC3-1 ( bottom). b) Fourier spectrograms illustrating the prewhitening process for SMC3-1. The dotted line indicates the significance level determined as explained in the text. This figure and subsequent ones should be read from up to down and from left to right.|
|Figure 2: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-2. b) Fourier spectrograms illustrating the prewhitening process for SMC3-2. The dotted line indicates the significance level determined as explained in the text.|
This member of the SMC open cluster NGC 242 was previously known as the H emission object LHA 115-S 3 (Henize 1956). It was classified "type-3'' by M02 due to the periodic pattern in the OGLE light curve with a periodicity of 118 7 days (Fig. 1a). The spectrum shows double-peaked H emission with equivalent width -99 Å and the red peak higher than the violet one (V < R). We also observed H emission with equivalent width -8.3 Å, higher-order Balmer lines filled by emission, and many metallic absorption lines. All this evidence supports the existence of an optically thin gaseous envelope around the star. From the blue spectral region, we derived a spectral type between A7 and F0. The red spectrum is similar to those of SMC3-4 (classified F5 I, see below). We conclude that the spectral type is between A7 and F5 but that the absolute magnitude and luminosity are uncertain. A Gaussian fit to the Na D lines observed in the red spectrum gives equivalent widths of Na D 1 = 1.25 Å and Na D 2 = 0.86 Å, i.e., a ratio marginally closer to 1.1. This could suggest that the star is embedded in gas of high optical depth (e.g. Munari & Zwitter 1997). At our resolution, the Na D lines appear as single components, so we could apply the calibration between their equivalent width and the reddening found by Munari & Zwitter (1997). However, the extremely large equivalent width should imply E(B-V) > 1, yielding an very unrealistic blue color. We conclude that the Na D line likely include multiple, not resolved, lines of probably interstellar origin.
Our method of searching periodicities yielded 7 significant frequencies (Fig. 1b); the most significant are 0.00844 0.00044 c/d, 0.06498 0.00002 c/d, and 0.01261 0.00005 c/d. The other four frequencies, 0.02115 0.00007 c/d, 0.00416 0.00204, 0.02932 0.00975, and 0.05439 0.00565 c/d, improved the fit only marginally. We first suspected that only the existence of the first three frequencies was well supported; F2 improved by 26% and F3 by 14%. However, all frequencies, except F7, are also present in the OGLE-III dataset, confirming the presence of 6 periods in this variable star.
The H emission with equivalent width of around -100 Å could suggest a Herbig Ae/Be object. The pre-main sequence classification is in principle consistent with the membership of this star in the open cluster NGC 242 and the relatively young cluster age, viz. 20 10 Myr (Ahumada et al. 2002). We collected J,H,K magnitudes from the 2MASS survey and found J-H = 0.20 0.06 and H-K = 0.07 0.06. Although these figures could indicate a small or null infrared excess (a galactic F0 supergiant has J-H = 0.15 and H-K = 0.03, Koornneef 1983), this is not enough evidence to reject the Herbig Ae/Be classification in the low metallicity environment of the SMC (de Wit et al. 2003, 2005). However, we note that the amplitude of variability (0.02 mag) is small for Herbig Ae/Be objects, which typically show irregular variability with amplitude of 0.1 mag.
No previous study was found for this object. The OGLE data indicate a periodic star with a light curve of modulated amplitude, produced primarily by the beating of two closely-spaced frequencies, 0.05733 c/d (17.44 d) and 0.06347 c/d (15.76 d) (Fig. 2a, M02). No other star in our sample shows such a remarkable phenomenon. The spectrum shows H emission with an equivalent width of -6.5 Å and other Balmer lines filled by emission. We found a spectral type of B7-9 III e, but the spectrum lacked the signatures to determine the color and absolute magnitude. The infrared colors for this star compiled from the 2MASS catalogue are J-H = 0.32 0.09 and H-K = 0.11 0.09, which when compared with typical colors for a (galactic) B9 III star, namely J-H = -0.02 and H-K = 0.02 (Winkler 1997), indicate significant reddening at the J-H color.
Our frequency-searching method yielded 5 frequencies for this star below the Nyquist threshold of 0.92 c/d. These are listed in Table 3. Although all these frequencies satisfy our criterion having amplitudes more than 4 times the noise level, only F1, F2, and F3 produce a rapid and significant improvement in the fit quality (as measured by ). The calculated fit reproduces the beating observed in the light curve very well (Fig. 2a). In spite of the fact that F4 and F5 improve the fit only marginally, all frequencies are confirmed by OGLE-III data, which confirms the multiperiodicity of this star.
We note the similarity between the light curve of SMC3-2 and those of RX UMa, a galactic RV-Tau type semiregular variable (Kiss et al. 2000). The variability of these stars differs only by a scale factor, the variability of RX UMa having larger amplitude and longer time scale than that of SMC3-2. Kiss et al. (2000) explain the behavior of RX UMa as a beating of two closely-separated modes of pulsation. According to these authors, the period ratio (close to the ratio 1.107 found in SMC3-2) suggests either high-order overtone or radial+non-radial oscillation. Contrary to SMC3-2, however, semiregular variables are usually highly unstable in coherence and amplitude of variability and have cooler spectral types.
|Figure 4: a) OGLE light curve ( upper) and spectrum ( bottom) of SMC3-4. b) Fourier spectrogram for SMC3-4. c) Light curve folded with the period of 184 days. d) O-C diagram for local minima outside eclipse; the test period is 23.5658 days.|
We see three possible interpretations for the beating phenomenon: (i) it could be caused by modulated light variations in a three-body system; (ii) it could be caused by two protoplanets surrounded by dust shells orbiting at nearby radii as measured from the central star; and (iii) it could reflect some kind of stellar pulsation. Alternative (i) is unlikely since the volume of the putative B-type binary probably excludes another possible binary with a period on the order of 15 days. Alternative (ii) is attractive and predicts larger amplitudes in redder passbands. This is what we really observe. The analysis of the MACHO instrumental b and r photometry indicates that the b-r color oscillates with the same two main frequencies as the I-band magnitude, but with amplitudes of 0.007 0.001 mag (F1) and 0.010 0.001 mag (F2), and the star is redder when brightest.
The problem with this view is that two protoplanets orbiting at close orbits would be gravitationally unstable, hence the probability of finding such a system is quite low. Finally, alternative (iii) is hard to admit, since there is no pulsation mechanism known in B-type stars acting with such a time scale of variability (Pamyatnykh 1999). If some kind of unknown pulsation is assumed, along with binarity, then the modulation could reflect the misalignment between the orbital axis and the pulsational axis, produced by the tidal force of the faint stellar component. In fact, Reed & Brondel (2005) report modulated light curves in these cases, very similar to that shown in Fig. 2. Although their analysis was performed for close binary systems with short orbital periods, we can scale their results to our system. Their results show that for some inclinations (see their Fig. 1) two closely spaced frequencies are obtained as the orbital aliases of the original pulsational frequency. In our case, the pulsational frequency should be (0.05733 + 0.06347)/2 c/d, i.e. 0.0604 c/d (16.556 days), and the orbital frequency 0.00307 c/d, i.e. 325.73 days. In this scenario, the doubt about the origin of the pulsation still remains, as does the question of which component could tilt the pulsational axis of the B-type star in spite of the large orbital separation. It is also possible that the pulsation axis is tilted with respect to the rotational axis without invoking the assumption of a binary companion, as in rapidly oscillating Ap (roAp) stars (Baldry et al. 1998). However, this case should produce frequency triplets, not doublets.
We note that the red spectrum of SMC3-2 shows peculiar absorption lines at 6125.6, 6197.1, 6300.1, 6319.2, and 6431.5 Å, which we have tentatively identified (considering a stellar velocity of 137 km s-1) as Ca I 6122.2 (it could also be Mn II 6122.4, 6122.8 and 6123.2), Zr III 6193.7 (or Sc I 6193.7), Fe I 6297.8 (or Ti I 6296.6, 6298.0), Ni I 6316.6 (or Fe I 6315.8), and Fe I 6428.8. These lines could be formed in the gaseous envelope, but they could also indicate enhanced chemical abundances as in magnetic Bp stars. Further high-resolution spectroscopy is needed to test the possible magnetic nature of this variable.
The OGLE light curve shows a noisy periodic pattern (Fig. 3a), and M02 report a low amplitude periodicity of 99 4 days. No emission is detected in our spectrum, which indicates a spectral type of F0 Ib-II. In agreement with M02, we found only one significant frequency, F1 = 0.00983 0.00005, which is also present in the OGLE-III database. Fourier spectrograms are given in Fig. 3b.
This variable was previously classified as the emission line object Lin-81 (Lindsay 1961). The OGLE light curve reveals an unusual eclipsing binary with a period of 184 18 days (Figs. 4a,b, M02). The depth of the minima and the general appearance of the light curve are quite variable. The presence of additional oscillations outside the principal eclipse makes it difficult to identify the secondary eclipse. However, an inspection of the phase diagram, constructed with the ephemeris P = 184 days and = 50 523.7769 (MJD), suggests that the secondary minimum, if present, should be >0.18 mag shallower than the main minimum. This indicates that the fainter component should have at most 83% of the temperature of the brighter component. The phase curve also shows that the oscillations outside the main minimum are not modulated with the binary period (Fig. 4c). Figure 4d shows the O-C diagram (observed minus calculated, see Sterken 2005) for these local minima, with a test period of 23.5658 days; some outliers in the graph suggest quasi-periodicity for this short-term variability.
The H emission is double-peaked, with V R, and an equivalent width of -31 Å. Their wings extend to -760 km s-1 and 700 km s-1 from the line centre. The double H emission is more consistent with emission from a rotationally supported, disc-like envelope, than from a stellar wind. As we do not have blue spectra for this star, we could not use the method described in Sect. 3.1 to find the spectral type. As an alternative method, we compared the red spectrum with the library of stellar spectra of Jacoby et al. (1984), excluding the regions around the H and the Na D lines. This yielded a spectral type of F5 I. A double Gaussian fit to the sodium D lines gives equivalent widths of Na D 1 = 1.15 Å and Na D 2 = 0.93 Å, i.e., a ratio close to 1.1. This could suggest that the star is embedded in gas of high optical depth (Munari & Zwitter 1997). We do not know if these lines are single or multiple, so we cannot use the relationship between equivalent width of single Na D lines and reddening given by these authors. If the lines were single, then their equivalent widths would indicate a very large reddening, larger than E(B-V) = 1, but this yields an unrealistically blue color for the star. This suggest that the lines are formed at different regions and are multiple. This is consistent with the fact that the equivalent width of the Na D lines is larger than expected for a normal F5 I star. A contribution could come from the envelope responsible for the H emission. However, the velocity of the Na lines, viz. 79 km s-1, is different from those of H (126 km s-1), suggesting the presence of another medium, maybe interstellar, or the companion star (see below), in the formation of the sodium D doublet.
That the eclipses are of variable depth and shape (Fig. 4c), and the strong variability of the light curve together suggest that the system could be an interacting spectroscopic binary. If we assume that the F5 supergiant is surrounded by an accretion disc fed by mass transferred from a secondary star filling its Roche-lobe, then the orbital period constrains the mean density of the mass donor star. Using standard formulae for close binary systems (e.g. Frank et al. 2002), we obtain 5.8 10-6 g cm-3 for the mean density of the secondary star, which corresponds to a G0 supergiant.
M02 report a low-amplitude periodicity of 30.0 0.3 days (Fig. 5a). No emission was detected in our spectrum, which indicates a spectral type of B7-8 III. We found two frequencies, 0.03321 0.00003 and 0.00135 0.00322. F2 accounts for the long-term variability, while their large error and low amplitude, close to the noise level (Fig. 5b), suggest that F1 is the only real frequency in this star. However, neither F1 nor F2 are present in OGLE-III data. This suggests that we have detected a transient periodicity in this star.
The OGLE light curve shows long-term periodic variability (Fig. 6a). M02 report a period of 197 20 days. No emission is detected in our spectrum, which indicates a spectral type of A5 II. In agreement with M02, we found only one frequency, 0.00509 0.00155 c/d (Fig. 6b). However, this frequency is not reproduced by OGLE-III data, so we conclude that this star probably shows a transient periodicity.
The OGLE light curve shows a noisy light curve with a rapid low-amplitude periodic oscillation (Fig. 7a). M02 gives a periodicity of 27.32 0.40 days, which agrees with our finding of only one frequency, 0.03661 0.00105 c/d (Fig. 7b), but this is not reproduced by OGLE-III data, so it could be a transient periodicity. No emission is observed in our spectrum, which indicates a spectral type of F4 IV, locating this star in the Milky Way halo.
This object was previously classified as the Wolf Rayet star AzV 81 (Massey 2002), and it corresponds to the star WR 4 in the list of Massey & Duffey (2001), who classify it as WN6. From RV studies, the star seems to be single, rather than a binary (Foellmi et al. 2003; Moffat 1988). Its wind properties were studied by Crowther (2000) from optical and infrared spectra. The OGLE light curve shows a periodic variability of variable amplitude. Our spectrum shows typical emission lines of WR stars, especially the strong He II 4686 line (Fig. 8a). A list of emission line strength for this star is given in Table 4. For this object we found 6 significant frequencies, which are listed in Table 3. M02 report a periodicity of 62 2 days, but a re-analysis of their data with the IRAF pdm routine gives 51.8 days, in agreement with our first frequency. The difference probably corresponds to a typo in the M02 table. Even with these frequencies, the fit could not account for the variability in several epochs (Fig. 8b); this suggests the presence of aperiodic components in the signal.
|Figure 8: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-8. b) Fourier spectrograms illustrating the prewhitening process for SMC3-8. The dotted line indicates the significance level determined as explained in the text. The high frequency at 0.41 c/d is probably an artifact.|
Foellmi et al. (2003) reported a periodicity of 6.55 days in the blue MACHO data for this star, and no periodicity in the OGLE light curve. We have no explanation for why they missed the 52-day periodicity and the other associated frequencies (we do not find the 6.55-day period in our analysis). It is then possible that the wrong light-curve was analyzed by these authors. Unfortunately, we had no additional data from OGLE-III to track the long-term behaviour of this star. We investigated the possibility that the 52-day periodicity could be the binary period. In this case it would be the longest period among WR stars in the SMC (Foellmi et al. 2003). Consequently, the radial-velocity amplitude could be low and eventually undetectable. However, using masses of 10 and 30 for the WR and O-type component, and an eccentricity of zero, the radial velocity half amplitude should still be detectable, between 60-150 km s-1 for inclination angles of 30-90 degres.
This object was classified as the H emission object LHA 115-N 41 (Henize 1956) and is also named AzV 85 (Azzopardi & Vigneau 1982). The OGLE light curve shows irregular variability that consists of modulated oscillations prior to HJD 2 451 600 and a sudden drop in brightness after this epoch (Fig. 9a). M02 report a period of 63 3 days. We found Hydrogen and Helium lines in emission and determined a spectral type of B1 II-IIIe. The H equivalent width is -2.2 Å (Fig. 9a), which is typical of Be stars. Fading events as the observed in SMC3-9 have also been reported in some galactic Be stars. For instance, the Be star 28 CMa showed, during March/April 1988, a sudden fading of a few weeks duration that was preceded by a quasi-periodic oscillation in time scales of 25 days (Mennickent et al. 1994). Our analysis of the data prior to HJD 2 451 600 indicates the presence of three frequencies, F1 = 0.00199 0.000272 c/d, which removes the long-term tendency, F2 = 0.01594 0.00288 c/d, and F3 = 0.01743 0.00006 c/d (Fig. 9b). Only F2 is recurrent in the OGLE-III time series, and it is also the frequency reported by M02.
The OGLE light curve reveals the periodic variability, which M02 reported as of low amplitude with a period of 185 17 days. Our spectra cover only the red range, H is seen as a weak absorption, and the continuum slope is rather blue (Fig. 10a). We detected only the frequency reported by M02, viz. 0.00541 0.00018 c/d (Fig. 10b), which is corroborated by OGLE-III data. However, the star could be an ellipsoidal variable if the true periodicity were twice F1. Based on the appearance of the phase diagrams, we favor this solution; we find an ephemeris for the times of a main minimum of P = 369.69 days and = 50 303.3773 (MJD) (Fig. 10a).
The OGLE light curve shows oscillations with variable amplitudes (Fig. 11a). M02 derive a period of 73 3 days. No emission is observed in our spectrum. We estimate the spectral type A3 III (Fig. 11a). The low radial velocities could indicate that the object is a possible member of the galactic halo, but they could also be due to orbital motion in a binary star. Our analysis gives two frequencies satisfying the significance criterion, F1 = 0.01378 0.00002 and F2 = 0.00688 0.00006 (i.e. 0.5 F1, Table 3). Both frequencies are confirmed by OGLE-III data. The solution with minima of alternating depth is compatible with a binary showing ellipsoidal variations with orbital frequency 0.00688 c/d (145.349 days) (Fig. 11a). We find an ephemeris for times of a main minimum of = 50 330.26981 (MJD).
|H + He II||-4.5|
|He II 4200||-1.0|
|H + He II||-1.9|
|He II 4517||-0.3|
|He II 4542||-1.3|
|N V 4610||-0.4|
|N III 4642||-2.0|
|He II 4686||-43.4|
|H + He II||-8.5|
|He II 5411||-6.9|
|C IV 5802+5812||-1.9|
|He I 5875||-1.5|
|H + He II||-33.3|
This object was previously classified as the H emission line star MA93 836 (Meyssonnier & Azzopardi 1993). The OGLE light curve shows complex periodic oscillations with a total amplitude of about 0.1 mag (Fig. 12a). M02 reported a period of 107 3 days. The spectral type is compatible with A0 Ibe. The H line shows an incipient emission at the red wing of the absorption line, i.e. a weak P-Cygni type profile, with maximum emission at velocity compatible with the H emission velocity (see Table 2). The Balmer absorption lines are blue-shifted by 190 km s-1 as regards the emission lines, which could be explained by unseen red emission in their absorption wings. The H equivalent width is -18.5 Å. We observe strong absorption lines at 5833.1, 5862.8, 5911.1, 5949.9, 5999.5, and 6081.1 Å (Fig. 12a). Using the velocity system of the H emission, we identified these lines as Fe I (except 5833.1 Å, which could not be identified). We also tried the velocity system of the blue lines, but it yielded inconsistent and unreliable identifications.
We found six frequencies, which fit the overall photometric behavior of the star relatively well, but do not reproduce some excursions in some epochs, revealing the presence of aperiodic or non-deterministic components in the light curve (Fig. 12a). We note that F2 = 2F1 for this star. F5 and F6 are not present in the OGLE-III data, hence they probably correspond to transient or spurious frequencies. The prewhitening process is illustrated in Fig. 12b.
photometry was presented by Massey (2002, its object SMC 41243). The OGLE light curve shows a slow simple oscillation (Fig. 13a). Our spectrum does not show emission lines, and it is compatible with a spectral type of B1 II-III. We found only one frequency, 0.00412 0.00132 c/d, which prewhites the data well (Fig. 13b). This frequency is not found in OGLE-III data. Their very low amplitude (3 mmag) suggests that is an spurious detection. The period given by M02 is erroneous.
This is the Be X-ray binary RX J0058.2-7231 (e.g. Schmidtke et al. 2003, hereafter S03). The OGLE light curve shows a series of low-amplitude ( mag) outbursts every 60-days and a sustained increase in the mean brightness (Fig. 14a). S03 suggest that the brightening may be related to changes in the Be star disc at periastron, assuming that the neutron star is in an eccentric orbit. They noted the narrow eclipse at maximum light, which explained as the eclipse of an excited region in the disc, by the B-type star, and provided the following ephemeris for minimum light (subtracting the long-term tendency): = 59.72 days and T0 = MJD 50 030. They also report evidence that the period may be decreasing at a rate of 0.01 days per cycle. Our analysis of OGLE data confirms the narrow eclipse at maximum and shows that the period is changing (Fig. 14a); furthermore, the joint analysis of MACHO and OGLE data, covering 74 orbital cycles, indicates a quasicyclic period change with a time scale of variability of 1200 days (Fig. 14a). Our spectra were taken during the bright phase, at phases 0.30 (blue) and 0.33 (red). At these phases we observe H and H emission, with equivalent widths of -15.5 Å and -1.7 Å, respectively. The continuum drops dramatically in a lapse of two days (Fig. 14a), more than the 0.1 mag expected from the light curve. We have doubts about the reality of this finding, since it could indicate a drop in the atmospheric transparency during our observation obtained in the red wavelength range, probably by the passage of a thin cloud. However, our records indicate that the weather was good and the sky clear that night, so we have no argument to support this explanation.
The full width at zero intensity of the H line is 976 km s-1. S03 predicted, according to their proposed eccentric binary model, the appearance of He II 4686 during the bright phase, and also the strengthening of the H emission. However, we observed neither He II in our spectra nor stronger H emission. Comparing with the emission line data given by S03 (based on two spectra taken at phases 0.04 and 0.70), we conclude that the Hequivalent width seems to be constant during the 60-day cycle, while their width does not. The three available measurements suggest that the H emission width is maximum at minimum light, just before starting the outburst, and then decreases during the slow outburst decline. This could be explained as the formation of a reservoir of high-velocity gas around the compact object prior to outburst. The optical outburst could be due to the accretion of this gas reservoir onto the compact object and the subsequent reprocessing of high-energy photons in the Be-star envelope. More observations are needed to confirm this scenario. We determined a spectral type of B0 III e from the appearance of the blue absorption features in the spectra.
We removed the long-term tendency of the light curve by fitting a straight line to the data, and worked with the residuals. The four frequencies we found (Table 3) reproduce the shape of the light curve well, but only F1 and F2 are present in OGLE-III data, and they probably represent the fitting to a non-sinusoidal light curve rather than true physical frequencies (Fig. 14a). The prewhitening process is illustrated in the periodograms of Fig. 14b.
|Figure 14: a) Detail of the light curve and fit ( upper), the frequency as a function of time ( middle) and the spectrum ( bottom) of SMC3-16. In the middle panel filled squares correspond to our data and open circles to data from Schmidtke et al. (2003). b) Fourier spectrograms illustrating the prewhitening process for SMC3-16. The dotted line indicates the significance level determined as explained in the text.|
The period change in this binary could be compared with the X-ray variability predicted by recent models of Be X-ray binaries by Hayashaki & Okazaki (2005, 2004). X-ray variability could modulate the optical light curve through re-processing of high energy photons by the Be star envelope. The aforementioned authors studied the long-term evolution of accretion discs around the neutron star in Be X-ray binaries using a three-dimensional, smoothed particle hydrodynamics (SPH) code. They showed that a time-dependent accretion disc is formed around the neutron star and that the disc has a one-armed spiral structure that is induced by a phase-dependent mass transfer from the Be disc. Hayashaki & Okazaki (2005) find that, after the disc is fully developed, the peak mass-accretion rate is distributed at phases around apoastron, which they interpret as evidence that the mass accretion rate is basically caused by an inward propagation of the one-armed spiral wave. According to their model, the time scale of variability for the X-ray luminosity peak is around 24 orbital cycles. This number is remarkably close to our number (1200 days corresponds to 20 orbital cycles). However, the distribution of the luminosity peak in the models of Hayashaki & Okazaki amounts to a few tenths of the orbital cycle, much more than our detected variations of a few percents of the orbital cycle. The Hayashaki & Okazaki simulations also show the interesting result that a double-peaked X-ray maximum is observed when the disc is growing to maximum size. In contrast to the persistent double maxima observed in the optical light curve of RX J0058.2-7231, the structure seen in the simulations disappears when the disc is fully developed, after a few tenths of the orbital cycles.
This is the emission line object AzV 241 or Lin 363 (Azzopardi & Vigneau 1982; Lindsay 1961). The OGLE light curve shows complex, presumably periodic, variability (Fig. 15a). M02 report a period of 35 3 days. Our spectrum shows H emission with an equivalent width of -2.0 Å, with the other Balmer lines probably partly filled by emission. Weak helium lines are also visible in the spectrum. We determined a spectral type of B8-9IIIe. Our search for frequencies yielded five significant frequencies for this star (Table 3). In general these frequencies fit the overall shape of the light curve well, but several outliers at different epochs indicate aperiodic components in the data. The only frequency of this set that survives during the OGLE-III period is F1. The prewhitening process is shown in the periodograms of Fig. 15b.
This is the emission line object AzV 366 or Lin 437 (Azzopardi & Vigneau 1982; Lindsay 1961). The OGLE light curve shows periodic variability (Fig. 16a). M02 reported a periodicity of 44 3 days. The spectrum shows prominent H and H emission lines and the spectral type is B1 Ib-IIe. The H equivalent width is -46.2 Å. Only one frequency is found in our analysis, 0.02249 0.00006 c/d, which is also observed in the OGLE-III dataset, confirming the M02 report. Aperiodic activity is also visible in some epochs (see outliers in Fig. 16a). The Fourier periodograms are shown in Fig. 16b.
This is the emission line star AzV 364, also named MA93 1492 (Azzopardi & Vigneau 1982; Meyssonnier & Azzopardi 1993). The OGLE light curve shows a periodic variability of 44 3 days (Fig. 17a, M02). Our spectrum shows weak H emission and H partly filled by emission. The H line has an equivalent width of -12 Å. We derived a spectral type of B2 IIIe. We found four significant frequencies (Table 3), but only F1, those found by M02, survives the stability test. The Fourier periodograms are shown in Fig. 17b.
In this section we present and discuss the infrared photometry extracted from the 2MASS catalogue using the VizieR on-line query form (Ochsenbein et al. 2000). We estimated the total J-H color excess using the expected colors for the spectral types and luminosity classes given in Sect. 3.1. The results are shown in Table 5. We did not attempt to determine absolute magnitudes from 2MASS photometry, since the large amount of infrared free-free and free-bound emission for stars with gaseous envelopes is hard to estimate with our data.
In Fig. 18 we show the color-color diagram for all objects, along with the position for classical Be stars and the main sequence and supergiant galactic tracks. We observe that the stars showing multiple periods and those with transient periodicities lie close to the main sequence and supergiant tracks, whereas those with single stable periods and the eclipsing and ellipsoidal binaries show much more reddening in general. The star SMC3-5 (those with the largest error bar in the H-K color in Fig. 18) has a flag in the 2MASS catalogue indicating photometric confusion probably by a nearby companion, therefore its infrared photometry should be taken with caution.
|Figure 18: The color-color diagram for the stars studied in this paper. The objects showing transient periods are indicated by open circles, those with single stable periods by filled circles, those with multiple periods by open squares, and the ellipsoidal and eclipsing binaries are indicated by asterisks. The colors for galactic main sequence and supergiant O9-G0 stars are also shown by the upper and bottom lines, respectively (Koornneef 1983). The dashed oval indicates the region occupied by the classical Be stars of Dougherty et al. (1994). The arrow shows the reddening vector for RV = 3.3 and a visual extinction of 2.5 mag.|
In the preceding section we have given a detailed discussion of the observed properties of each star in our sample. In this section, we focus on the general picture we can extract from this study. For this reason, in the following we exclude the objects SMC3-8 and SMC3-16, as well as the eclipsing and ellipsoidal variables, which were extensively discussed in the previous section, and will focus our attention on the remaining objects of our sample.
A general pattern seems to arise from the analysis in Sect. 3. We identify three general groups of stars according to their photometric and spectroscopic properties. The first group corresponds to those stars showing transient periods; they are 4 non-emission stars showing very small photometric amplitudes, except SMC3-6 displaying an amplitude of 0.05 mag. The second group consists of stars showing only one stable frequency; it corresponds to 4 giant Be stars and one F0Ib-II star. Finally, the third group displays multiperiodicities and consists of 2 A-type and 1 B7-9e stars.
It is possible that the variability for the first two groups arises from emitting gaseous envelopes and not from the star itself. Transient objects could be explained as those having low density, highly variable, low-emission envelopes. This should explain the low amplitude of the variability and their transient nature. It is well known that episodes of mass ejection in Be stars produce brightness variability. If these ejections consist of gas blobs taking angular momentum away from the star, they should rotate producing a quasi-periodic pattern of brightness before dissolving in the interstellar medium. This view requires gas ejections from the stellar surface, a process normally associated to rapid rotation, which should work in the B-A stars of our sample (given the required conditions), but probably not for SMC3-7 having a spectral type of F4 IV. This explanation is supported by the finding that a much larger fraction of massive stars at lower Z environments (like that found in the ) reaches break-up velocities. The reason is that at low Z, mass loss is smaller and the removal of angular momentum during evolution is much weaker (Maeder & Meynet 2001).
On the other hand, well-developed envelopes, which are a signature of classical Be-type stars, could be responsible for the stable, long-term variability observed in most Be stars of Table 2. This variability could be explained as changing aspects of an elliptical envelope in a binary system (by precession, for instance), or it eventually could be caused by ellipsoidal variability masked by additional phenomena related to the formation and dissipation of the Be-star envelope. The amplitudes of the variations are in principle compatible with ellipsoidal variability. Therefore, these objects are serious candidates for binarity. While Be star binaries amount to 20% of the total Be-star population in our galaxy (e.g. Porter & Rivinius 2003), a corresponding figure for the has not yet been established. Therefore, these objects are ideal targets for RV studies aimed at confirming their binary nature and at testing binary formation models in different metallicity environments (e.g. de Loore & Vanbeveren 1995).
Finally, the last group of variables in this study corresponds to stars showing multiple periods. As noted before, the time scale of variability does not fit any known stellar pulsation mechanism. However, all of them are emission line stars, so the variability could be associated to the envelope. Binarity could be the right scenario for these objects, and future time-resolved radial velocity studies are needed to check this view.
In this paper we have presented a photometric and spectroscopic study of 17 bright blue and yellow variable stars in the direction of the . These stars were previously classified type-3 by Mennickent et al. (2002); i.e., they correspond to photometrically periodic variables, Be star candidates, whose nature was not completely clear. During our analysis, we detected peculiar behavior in several light curves and interesting new features in the spectra of these stars. We summarize our findings as follows:
Further studies are clearly needed to improve our knowledge about these objects. For instance, monitoring radial velocities over extended periods of time should be desirable to test the binary hypothesis. To fully understand the variable stars discovered by programs like OGLE and MACHO, time-resolved spectroscopic and photometric data must be analyzed together. Part of our future research programs will focus on this goal.
We acknowledge the anonymous referee who helped us to improve the first version of this article. R.M. acknowledges support by Fondecyt grant 1030707. R.M., W.G. and G.P. acknowledge support from the Chilean FONDAP centre for Astrophysics 15010003. B. Sabogal acknowledges support by Programa MECE Educación Superior UCO0209. This research has made use of the SIMBAD database, operated at the CDS, Strasbourg, France. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis centre, funded by the National Aeronautics and Space Administration and the National Science Foundation.
|Figure 3: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-3. b) Fourier spectrograms illustrating the prewhitening process for SMC3-3. The dotted line indicates the significance level determined as explained in the text.|
|Figure 5: a) OGLE light curve and the best fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-5. b) Fourier spectrograms illustrating the prewhitening process for SMC3-5. The dotted line indicates the significance level determined as explained in the text.|
|Figure 6: a) OGLE light curve and the best-fitting function ( upper) and spectrum ( bottom) of SMC3-6. b) Fourier spectrograms illustrating the prewhitening process for SMC3-6. The dotted line indicates the significance level determined as explained in the text.|
|Figure 7: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-7. b) Fourier spectrograms illustrating the prewhitening process for SMC3-7.|
|Figure 9: a) OGLE light curve and the best-fitting function ( upper) and spectrum ( bottom) of SMC3-9. b) Fourier spectrograms illustrating the prewhitening process for SMC3-9. The dotted line indicates the significance level determined as explained in the text.|
|Figure 10: a) OGLE light curve and the best-fitting function ( upper), the light curve folded with the period 184.84 days and 369.69 days ( middle panels) and spectrum ( bottom) of SMC3-10. We suggest that the second periodicity is the right one. b) Fourier spectrograms illustrating the prewhitening process for SMC3-10. The dotted line indicates the significance level determined as explained in the text.|
|Figure 11: a) OGLE light curve and the best-fitting function ( upper), light curve folded with the orbital period ( middle, the x-axis units are cycles) and spectrum ( bottom) of SMC3-12. The peaks at 4659 Å and 4772 Å are artifacts not removed by the reduction process; they probably correspond to hits of cosmic rays. b) Fourier spectrograms illustrating the prewhitening process for SMC3-12. The dotted line indicates the significance level determined as explained in the text.|
|Figure 12: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-14. b) Fourier spectrograms illustrating the prewhitening process for SMC3-14. The dotted line indicates the significance level determined as explained in the text.|
|Figure 13: a) OGLE light curve and the best-fitting function ( upper) and spectrum ( bottom) of SMC3-15. b) Fourier spectrograms illustrating the prewhitening process for SMC3-15. The dotted line indicates the significance level determined as explained in the text.|
|Figure 15: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-18. b) Fourier spectrograms illustrating the prewhitening process for SMC3-18. The dotted line indicates the significance level determined as explained in the text.|
|Figure 16: a) OGLE light curve and the best-fitting function ( upper) and spectrum ( bottom) of SMC3-20. b) Fourier spectrograms illustrating the prewhitening process for SMC3-20. The dotted line indicates the significance level determined as explained in the text.|
|Figure 17: a) OGLE light curve and the best-fitting function ( upper), detail of the light curve and fit ( middle) and spectrum ( bottom) of SMC3-21. b) Fourier spectrograms illustrating the prewhitening process for SMC3-21. The dotted line indicates the significance level determined as explained in the text.|
|Object||ID||HJD (blue-spec)||exptime (s)||HJD (red-spec)||exptime (s)|
|OGLE00445466-7328029||SMC3-2||56.49586||600||58.54036, 58.54399, 58.54763||300|
|OGLE00455414-7314043||SMC3-3||56.50746||900||59.64039, 59.64402, 59.64766||300|
|OGLE00463376-7312043||SMC3-4||59.50340, 59.51515, 59.52640||900|
|OGLE00504344-7327053||SMC3-8||56.54504||900||58.52697, 58.53061, 58.53425||300|