A&A 457, 237-248 (2006)
DOI: 10.1051/0004-6361:20065449

Pulsating pre-main sequence stars in IC 4996 and NGC 6530

K. Zwintz - W. W. Weiss

Institute of Astronomy, University of Vienna, Türkenschanzstrasse 17, 1180 Vienna, Austria

Received 18 April 2006 / Accepted 7 June 2006

Abstract
Context. Asteroseismology of pulsating pre-main sequence (PMS) stars has the potential of testing the validity of current models of PMS structure and evolution. As a first step, a sufficiently large sample of pulsating PMS stars has to be established, which allows to select candidates optimally suited for a detailed asteroseismological analysis based on photometry from space or ground based network data.
Aims. A search for pulsating PMS members in the young open clusters IC 4996 and NGC 6530 has been performed to improve the sample of known PMS pulsators. As both clusters are younger than 10 million years, all members with spectral types later than A0 have not reached the zero-age main sequence yet. Hence, IC 4996 and NGC 6530 are most suitable to search for PMS pulsation among their A- and F-type cluster stars.
Methods. CCD time series photometry in Johnson B and V filters has been obtained for IC 4996 and NGC 6530. The resulting light curves for 113 stars in IC 4996 and 194 stars in NGC 6530 have been subject to detailed frequency analyses.
Results.$\delta$ Scuti-like PMS stars have been discovered in IC 4996 and 6 in NGC 6530. For another PMS star in each cluster, pulsation can only be suspected. According to the computed pulsation constants, the newly detected PMS stars seem to prefer to pulsate in a similar fashion to the classical $\delta$ Scuti stars, and with higher overtone modes.

Key words: techniques: photometric - stars: pre-main sequence - $\delta$ Scuti - open clusters and associations: individual: IC 4996 - open clusters and associations: individual: NGC 6530

1 Introduction

The study of the first stages in the formation of stars is currently one of the most active research fields in stellar astronomy. Pre-main sequence (PMS) stars are contracting from the birthline to the zero-age main sequence (ZAMS) in the Hertzsprung-Russell (HR) diagram. They interact with the circumstellar environment in which they are still embedded; hence they are characterized by a large degree of activity, strong near- or far-IR excesses, and very often by emission lines. Two major groups can be distinguished: T Tauri and Herbig Ae/Be objects. Members of both groups show photometric and spectroscopic variabilities on timescales from minutes to years, indicating that stellar activity begins in the earliest phases of stellar evolution, prior to the arrival on the main sequence.

T Tauri stars are newly formed low-mass stars that have recently become visible in the optical range. As they have primarily spectral types ranging from late F to M, they are not the prime targets to search for PMS pulsation. The Herbig Ae/Be (HAEBE) stars are intermediate-mass (1.5-10 $M_{\odot}$) objects (e.g., Herbig 1960; Finkenzeller & Mundt 1984), which are either in their PMS phase ( $M < 6~M_{\odot}$) or have already arrived on the ZAMS ( $M \geq 6~M_{\odot}$). They have spectral types B and A, are characterized by emission lines, and mostly show excess radiation in the infrared. PMS stars with masses between 1.5 and 4 $M_{\odot}$ have the right combination of temperature, luminosity, and mass to become vibrationally unstable, and hence are the main targets in which to search for pulsation.

During their evolution to the main sequence, young stars move across the instability region in the HR diagram, which suggests that part of their activity is due to stellar pulsations. PMS stars differ from their counterparts with the same effective temperature, mass, and luminosity, but that have already evolved off the main sequence, mostly in the inner regions, although their atmospheres are quite similar (Marconi & Palla 1998). The discovery of pulsating PMS stars is important because it allows to constrain the internal structure of young stars and to test evolutionary models.

As the evolutionary tracks for pre- and post-main sequence stars intersect each other several times (e.g., Breger & Pamyatnykh 1998), the determination of the evolutionary stage of a field star may be ambiguous. Additional information, like the age or distance of the star, is needed to resolve this ambiguity. Young clusters are most suitable to search for PMS pulsators because all members have the same age and distance, hence confusion with more evolved objects can be reduced.

The first two pre-main sequence pulsators, V 588 Mon (HD 261331, NGC 2264 2) and V 589 Mon (HD 261446, NGC 2264 20), were discovered by Breger (1972) in the young open cluster NGC 2264. It took more than 20 years before the discovery of the next PMS pulsating star. Kurtz & Marang (1995) observed pulsations in the well-studied Herbig Ae field star HR 5999 (HD 144668, V856 Sco). Since then, several detections of pulsations in Herbig Ae field stars (e.g., Donati et al. 1997; Marconi et al. 2000; etc.), as well as in members of young clusters (e.g., Zwintz et al. 2005) have been made.

In Sect. 2, the general characteristics of the two open clusters IC 4996 and NGC 6530 are summarized; Sects. 3 and 4 describe the observations, data reduction, and frequency analysis. The pulsating PMS stars in IC 4996 and NGC 6530 are shown in Sects. 5 and 6, respectively, and the calculation of the according pulsation constants is given in Sect. 7. As all stars in the fields of the clusters were analyzed, the other detected variable stars are listed in Sect. 8. A summary of the cluster properties derived from our analysis can be found in Sect. 9, while the conclusions of this paper are specified in Sect. 10.

2 The young open clusters IC 4996 and NGC 6530

IC 4996 ( $\alpha_{2000} = 20^{\rm h} 16\hbox{$.\!\!^{\rm m}$ }5, \delta_{2000} = +37\hbox{$^\circ$ }38\hbox{$.\mkern-4mu^\prime$ }0$) is located in the direction of Cygnus, 40 pc above the galactic plane, and is part of a large region with active star formation that contains other young open clusters and Wolf-Rayet stars. An IRAS map of the region (Lozinskaya & Repin 1990) shows the presence of a dusty shell that surrounds the cluster.

The age and distance estimates for IC 4996 from different authors (e.g., Alfaro et al. 1985; Delgado et al. 1998) agree: the cluster is slightly younger than 107 years and is located $\sim $1.7 kpc from the Sun. In particular, the inferred age indicates the likely existence of pre-main sequence members in the range of spectral types A and F. Delgado et al. (1999) performed spectroscopic observations with the aim of estimating radial velocities and spectral types for 16 proposed PMS stars to confirm or reject their cluster membership. Additionally, the authors found spectral features indicative of the PMS nature of the stars with spectral types later than A0.

NGC 6530 ( $\alpha_{2000} = 18^{\rm h} 04\hbox{$.\!\!^{\rm m}$ }5, \delta_{2000} = -24\hbox{$^\circ$ }22\hbox{$.\mkern-4mu^\prime$ }0$) is located in the central part of the HII region M8, the Lagoon nebula. Since the first study performed by Trumpler in 1930, several investigations have been devoted to the study of this cluster. The authors agree to an age of 1.0-3.0 Myr (Kilambi 1977; Sung et al. 2000), a distance of 1.86 $\pm$ 0.07 kpc (Mc Call 1990), and a color excess E(B-V) = 0.39 $\pm$ 0.09 mag (Chini & Neckel 1981).

Kilambi (1977) obtained $\it UBV$ photographic photometry of NGC 6530 and found that all stars fainter than V = 12.0 mag have not reached the ZAMS yet. Van den Ancker et al. (1997) obtained Walraven WULBV, Johnson/Cousins UBV(RI), and near-IR JHK photometric data, and performed spectroscopy of NGC 6530.

3 Observations and data reduction

3.1 IC 4996

IC 4996 was observed with the 1.5 m telescope (f/8) at the Sierra Nevada Observatory, Spain, between Sept. 2 and Sept. 15, 2002, in Johnson B and V filters using a 1k $\times $ 1k CCD chip with a scale of 0.33 $\hbox {$^{\prime \prime }$ }$/px providing a field of view of 6 $\hbox {$^\prime $ }$ $\times $ $6\hbox {$^\prime $ }$. In total, 69.76 h of time-series photometry have been obtained in 12 out of 14 granted nights. Only the first ten nights (corresponding to 62.38 h of observations and 1990 science frames) were used for the analysis because during the last two nights the weather conditions were not good enough, resulting in an extremely high scatter in the light curves.

The images were already bias- and dark-corrected with the standard procedure used at the observatory. The flat field correction was performed within the reduction software PODEX written by Kallinger (2005). PODEX allows to extract the photometric signal of CCD time series photometry, using a combination of aperture photometry and point-spread function fitting. After the light curve is computed for each selected star, the mean value of the comparison light curve is subtracted, where the stars used for the comparison light curve can be selected arbitrarily. Optionally flat field and (color-dependent) extinction corrections can be applied. The output of the program not only contains the photometric signal for each integration, but also the values of the airmass and the comparison light curve. The advantage of PODEX is that variable stars can be identified easily and immediately discarded for the comparison light curve. This is specifically important if - as in our case - no a priori information about the variability of the analyzed stars is available.

The 113 stars selected for the further analysis (see Fig. 1) lie in the range between V = 11-18 mag. Nightly means were subtracted to correct for zero-point changes and long-term irregular light variations, which most likely are due to variable extinction by circumstellar dust. Color-dependent extinction was taken into account, as described in Zwintz et al. (2005).

Our own star numbers are used, but cross references with numbers used in the WEBDA database and given by Delgado et al. (1998, 1999) are listed, too.


  \begin{figure}
\par\includegraphics[width=6.55cm,clip]{5449fig1.eps}\end{figure} Figure 1: Schematic map of the observed field of IC 4996 (fov 6 $\hbox {$^\prime $ }$ $\times $ $6\hbox {$^\prime $ }$) with all stars measured in Johnson  B and  V, where 1 pixel (px) corresponds to 0.33 $\hbox {$^{\prime \prime }$ }$. Identifiers refer to the objects discussed in the text.
Open with DEXTER

3.2 NGC 6530

For NGC 6530, CCD photometric time series in Johnson B and V filters were obtained with the 0.9 m telescope (f/13.5) at the Cerro Tololo Interamerican Observatory (CTIO), Chile. Between Aug. 1-7 and Aug. 9-15, 2002, NGC 6530 was observed using the dedicated 2048 $\times $ 2046 SITe CCD chip, which provides a field of view of $\sim $13.5 $\hbox {$^\prime $ }$ $\times $ 13.5 $\hbox {$^\prime $ }$ with a scale of 0.396 $\hbox {$^{\prime \prime }$ }$/px. In total, 80.16 h of time-series photometry was acquired in 12 out of 14 nights. As the cluster was slightly too large to be observed on a single frame, two overlapping regions were chosen for which the observing time was shared. Out of 3437 scientific frames, 1601 were observed for field 1 and 1336 for field 2.

Bias subtraction and the creation of superflat images were performed using the IRAF ared.quad package. Flat field correction was performed with the reduction software PODEX, which was also used to extract the photometric signal, as it was for the data of IC 4996.

In total, 194 stars have been identified in both fields where 43 stars lie in the overlapping region, 79 only in field 1 and 72 only in field 2 (see Figs. 2 and 3). Nightly means were subtracted from the light curves to correct for zero-point changes and long-term irregular light variations caused by variable extinction due to circumstellar material. Again, the effect of color-dependent extinction was encountered and corrected using the same method as for IC 4996 and NGC 6383 (Zwintz et al. 2005).


  \begin{figure}
\par\includegraphics[width=6.4cm,clip]{5449fig2.eps}\end{figure} Figure 2: Schematic map of the observed field 1 of NGC 6530 (fov $\sim $ 13.5$^\prime $ $\times $ 13.5$^\prime $), South is at the top and East is to the left, with all stars measured in Johnson  B and  V, where 1 pixel (px) corresponds to 0.369 $^{\prime \prime }$. The area marked with the solid lines corresponds to the overlapping region. Identifiers refer to the objects discussed in the text.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.4cm,clip]{5449fig3.eps}\end{figure} Figure 3: Schematic map of the observed field 2 of NGC 6530, similar to Fig. 2.
Open with DEXTER

4 Frequency analysis

For all selected stars, detailed frequency analyses were performed using primarily P ERIOD98 (Sperl 1998), which is based on the Discrete Fourier Transformation (DFT, Deeming 1975) and provides a multi-sine fit option. Here, a signal was considered to be significant if it exceeds four times the noise level in the amplitude spectrum (Breger et al. 1993; Kuschnig et al. 1997).

Additionally, S IGS PEC (Reegen 2006), which computes significance levels for amplitude spectra of time-series with arbitrarily given sampling and white noise, has been used, providing a second independent method for detailed frequency analyses. S IGS PEC's advantage over classical techniques is the employment of analytically clean statistics instead of coarse signal-to-noise ratio estimations. A S IGS PEC-significance of $\sim $5.46 corresponds to a signal-to-noise ratio of 4 (Breger et al. 1993) mentioned before. Hence, a signal was considered to be significant within the S IGS PEC analysis, if its significance was higher than 5.46.

Only if frequencies appeared significant in both filters and using both methods are they believed to be intrinsic and not artefacts of the reduction.

The errors of frequencies and amplitudes (Table 2) were calculated using the relations given by Montgomery (1999):

 \begin{displaymath}%
\sigma(a) = \sqrt{\frac{2}{N}}~\cdot~\sigma(m)
\end{displaymath} (1)


 \begin{displaymath}%
\sigma(f) = \sqrt{\frac{6}{N}}~\cdot~\frac{1}{\pi~T}~\cdot~\frac{\sigma(m)}{a},
\end{displaymath} (2)

where $\sigma(m)$ is the rms magnitude of the data set, a the corresponding amplitude, N is the number of data points and T is the time base of the observations. It has to be noted that the errors calculated with the above formulae can only provide a lower limit of the true errors. If the errors in magnitude are correlated in time and the correlation time is of the order of the period of the signal, the error estimates are too low by a factor $\sqrt{A}$, where the function $A \approx 0.24 P / \Delta t$ with the separation of data points of $\Delta t$ (Montgomery 1999).

5 Pulsating PMS stars in IC 4996

40 PMS cluster members have been selected as prime candidates for pulsation in IC 4996 due to their position in the HR diagram, where the major criterion has been the location of the stars in the region of the classical instability strip (e.g., Pamyatnykh 2000). Two bona fide, IC 4996 37 and 40, and one suspected, IC 4996 46, pulsating pre-main sequence stars have been discovered in this cluster. The previously known parameters of these stars are given in Table 1 and their corresponding frequencies and amplitudes are listed in Table 2.

Table 1: Parameters of the bona fide and suspected PMS pulsators in IC 4996.

5.1 IC 4996 37

For IC 4996 37, Delgado et al. (1998) measured V = 15.30 mag, (B-V) = 0.8 mag, and (U-B) = 0.44 mag using CCD photometry. This coincides well with the observations performed by Vansevicius et al. (1996), who give V = 15.207 mag and (B-V) = 0.784 mag, but also (V-R) = 0.428 mag and (V-I) = 0.921 mag. Delgado et al. (1999) obtained long-slit spectra for 16 stars of IC 4996, among those also star 37, which was found to have a spectral type of A5. Hence, it is located in the region of the HR diagram, where pulsation can be expected, and this made it an ideal target to search for pulsation.

The frequency analysis with P ERIOD98 yielded a single intrinsic significant frequency in both filters at 31.875 c/d corresponding to a period of 45 min. This frequency was also detected using S IGS PEC with significances of 14.7 in V and 5.9 in B. In the amplitude spectra of both filters, the peak at 31.875 c/d can clearly be noticed; as an example, Fig. 4 shows the amplitude spectrum of the star in the V filter. Several frequencies related to the one-day alias appear to be significant in the amplitude spectra of both filters. However they were omitted in the analysis, as they result from the daily gaps in the observations due to the day-night-cycle (see the appropriate spectral window).

Table 2: Frequencies and amplitudes determined for the two bona fide PMS pulsators in IC 4996, where the errors in the last digits of the corresponding quantities are given in parentheses.


  \begin{figure}
\par\includegraphics[width=6.5cm,clip]{5449fig4.eps}\end{figure} Figure 4: Amplitude spectrum of IC 4996 37 in the V filter; the identified frequency is marked with an arrow.
Open with DEXTER

5.2 IC 4996 40

For IC 4996 40, Delgado et al. (1998) give V = 15.03 mag, (B-V) = 0.75 mag, and (U-B) = 0.43 mag. A spectral type of A4 was determined by Delgado et al. (1999), which places the star inside the instability region in the HR diagram.

The frequency analysis with P ERIOD98 resulted in the detection of a single significant frequency in both filters at 33.569 c/d corresponding to a period of 43 min (Fig. 5). The same frequency was also detected using S IGS PEC with significances of 33.1 in V and 22.2 in B. Again, a number of frequencies related to the one-day alias appear in the amplitude spectra and have been omitted in the analysis.

5.3 IC 4996 46

With (B-V) = 0.72 mag and V = 15.3 mag (Delgado et al. 1998), IC 4996 46 is also located in the region of the classical instability strip in the HR diagram. The frequency analysis resulted in the detection of different periods in V and B filters: two periods of approximately 6.0 and 6.5 h with amplitudes of 5.4 and 3.8 mmag in V and two periods of 4.8 and 3.7 h with amplitudes of 4.2 mmag each in B. Hence, the star is considered as a suspected PMS pulsator.

6 Pulsating PMS stars in NGC 6530

30 PMS cluster members have been selected as prime candidates for pulsation in NGC 6530 according to their position in the HR diagram. Again, the major criterion has been the location of the stars in the region of the classical instability strip. Six PMS pulsators were discovered with classical $\delta$ Scuti type frequencies (NGC 6530 5, 82, 85, 263, 278, and 281), and for a seventh star, pulsation can only be suspected (NGC 6383 288).

Table 3 shows an overview of all new pulsating PMS stars in NGC 6530, where V, (B-V), and (U-B) are taken from the WEBDA database, and the spectral types were derived by van den Ancker et al. (1997).


  \begin{figure}
\par\includegraphics[width=6.5cm,clip]{5449fig5.eps}\end{figure} Figure 5: Amplitude spectrum of IC 4996 40 in the V filter; the identified frequency is marked with an arrow.
Open with DEXTER

Table 3: Parameters of the bona fide and suspected PMS pulsators in NGC 6530.

6.1 NGC 6530 5

NGC 6530 5 is situated in the region where the two observed fields overlapped. For star 5 (WEBDA #159), no information about its spectral type or cluster membership is available, but $\it UBV$ CCD photometry gives V = 13.59 mag, (B-V) = 0.43 mag, and (U-B) = 0.26 mag (Sung et al. 2000). Hence, from the location in the cluster HR diagram, it is most likely a pre-main sequence member lying in the region of the classical instability strip.

Two significant frequencies at 46.596 c/d with amplitudes of 1.4 mmag in V and 1.8 mmag in B and at 53.417 c/d with amplitudes of around 1 mmag in both filters have been found and are clearly visible in the amplitude spectra (for the V filter the amplitude spectra are shown in Fig. 6).

6.2 NGC 6530 82

NGC 6530 82 (WEBDA #78) is located in field 1 of our observations, and again no spectral type or membership information is available in the literature. $\it UBV$ CCD photometry (Sung et al. 2000) yields V = 13.97 mag, (B-V) = 0.61 mag, and (U-B) = 0.30 mag, which makes it a likely cluster member and places the star inside the classical instability strip in the HR diagram.

In the analysis, three frequencies between 24.829 c/d and 38.531 c/d with amplitudes between 1.7 mmag and 2.8 mmag have been found to be significant (see Table 4). A simultaneous multi-sine fit with these three frequencies represents the shape of the light curve well. The frequencies are clearly detectable in both filters (the amplitude spectra of the V filter are given in Fig. 7).


  \begin{figure}
\par\includegraphics[width=6.8cm,clip]{5449fig6.eps}\end{figure} Figure 6: Amplitude spectrum of NGC 6530 5 in the V filter; "f1'' and "f2'' mark the identified frequencies.
Open with DEXTER

Table 4: Frequencies and amplitudes determined for all PMS pulsators in NGC 6530, where the errors in the last digits of the corresponding quantities are given in parentheses.

6.3 NGC 6530 85

NGC 6530 85 (WEBDA #53) was also situated in field 1 of our observations and was studied by different authors. Kilambi (1977) obtained V = 12.96 mag, (B-V) = 0.62 mag, (U-B) = 0.30 mag, and confirmed its cluster membership. These values match the $\it UBV$ CCD measurements of Sung et al. (2000), who found V = 13.07 mag, (B-V) = 0.65 mag, and (U-B) = 0.35 mag quite well.

The proper motion study by van Altena & Jones (1972) yielded a membership probability of 78% for NGC 6530 85 (i.e., their star #173). They also reported an E(B-V) = 0.35 mag and absolute magnitude MV = +0.81 mag. Chini & Neckel (1981) determined a spectral type of A0 for NGC 6530 85.

The frequency analysis using P ERIOD98 and S IGS PEC allowed us to identify simultaneous pulsation with five frequencies between 10.585 c/d and 31.148 c/d and amplitudes in the range of 39.1 mmag to 1.8 mmag (Table 4, Fig. 9). The shape of the light curve demonstrates the multi-periodic nature of the star well (Fig. 8).


  \begin{figure}
\par\includegraphics[width=7cm,clip]{5449fig7.eps}\end{figure} Figure 7: Amplitude spectrum of NGC 6530 82 in the V filter; the identified and sequentially prewhitened frequencies are marked with arrows.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.65cm,clip]{5449fig8.eps}\end{figure} Figure 8: Differential light curve of NGC 6530 85 with a five frequency multi-sine fit (solid line); top: V filter, bottom: B filter (shifted for better visibility).
Open with DEXTER

6.4 NGC 6530 263

NGC 6530 263 is one of the stars located in field 2 of our observations. $\it UBV$ CCD photometry (Sung et al. 2000) yielded V = 13.67 mag, (B-V) = 0.63 mag, and (U-B) = 0.42 mag for star 263 (WEBDA #57), which agree well with earlier measurements (e.g., from Kilambi 1977). No information on the spectral type of the star is available in the literature.

NGC 6530 263 was also contained in the proper motion study by van Altena & Jones (1972), but it belonged to the 135 stars that have not been included in their "primary'' absolute parameter solution. They indicate that the star (their star #178) might only be considered a probable cluster member.

The frequency analysis of both filters showed a single significant frequency at 19.223 c/d, corresponding to a period of 1.25 h, with an amplitude of 8.3 mmag in B and 7.1 mmag in V. This $\delta$ Scuti-like pulsation frequency clearly shows up in the amplitude spectra (e.g., for the V filter see Fig. 10).


  \begin{figure}
\par\includegraphics[width=7cm,clip]{5449fig9.eps}\end{figure} Figure 9: Amplitude spectra of NGC 6530 85 in both filters; the identified and sequentially prewhitened frequencies are marked with arrows.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.8cm,clip]{5449fi10.eps}\end{figure} Figure 10: Amplitude spectrum of NGC 6530 263 in the V filter; the identified frequency is marked with an arrow.
Open with DEXTER

6.5 NGC 6530 278

NGC 6530 278 (WEBDA #38) has also been one of the stars observed in field 2. Kilambi (1977) obtained V = 12.16 mag, (B-V) = 0.44 mag, and (U-B) = 0.39 mag for this star, reported its cluster membership, and found indication for variability. The more recent values for V = 12.17 mag, (B-V) = 0.53 mag and (U-B) = 0.36 mag from CCD photometry (Sung et al. 2000), confirm these values.

The proper motion study by van Altena & Jones (1972) gave a membership probability of 68% for NGC 6530 278 (their star #157). They also report a value of E(B-V) = 0.35 mag and absolute magnitude MV = -0.12 mag. Chini & Neckel (1981) used their $\it UBV$ and $H_{\beta}$ observations to derive a spectral type of A0 for this star.

The light curve of NGC 6530 278 (Fig. 11) clearly illustrates its multi-periodicity. The formal solution of the detailed frequency analysis yielded nine significant frequencies in both filters between 4.0 c/d and 15.6 c/d with amplitudes ranging from 13.8 mmag down to 2.6 mmag (see Table 4, Fig. 12). The identification and prewhitening of intrinsic frequencies was difficult due to the numerous peaks caused by aliases in the frequency range between 0 and 15 c/d. At least three more frequencies between 20 c/d and 30 c/d might exist, but the variability of NGC 6530 278 is proven.


  \begin{figure}
\par\includegraphics[width=7cm,clip]{5449fi11.eps}\end{figure} Figure 11: Differential light curve of NGC 6530 278 with a nine frequency multi-sine fit (solid line); top: V filter, bottom: B filter (shifted for better visibility).
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.9cm,clip]{5449fi12.eps}\end{figure} Figure 12: Amplitude spectra of NGC 6530 278 in the V filter; the identified and sequentially prewhitened frequencies are marked with arrows.
Open with DEXTER

6.6 NGC 6530 281

NGC 6530 281 (WEBDA #13) lies in field 2 of our observations. $\it UBV$ CCD photometry (Sung et al. 2000) gives V = 13.35 mag, (B-V) = 0.45 mag, and (U-B) = 0.27 mag for this star, which agree well with values derived earlier by different authors. No spectral classification is available, but Kilambi (1977) reported that the star is most probably a cluster member. Its location in the cluster CMD and in the HR diagram confirm this suggestion.

Pulsation with seven periods was found as a formal solution. The frequencies lie between 30.6 c/d and 43.4 c/d, corresponding to periods of 47 to 33 min, where the amplitude spectrum in the B filter is shown in Fig. 13. The amplitudes range from 5.1 mmag to 1.4 mmag (see Table 4). Similar to NGC 6530 278, the identification and prewhitening of intrinsic frequencies was difficult due to the numerous peaks caused by aliases in the frequency range between 35 c/d and 45 c/d.


  \begin{figure}
\par\includegraphics[width=7cm,clip]{5449fi13.eps}\end{figure} Figure 13: Amplitude spectrum of NGC 6530 281 in the B filter; the identified and sequentially prewhitened frequencies are marked with arrows.
Open with DEXTER

6.7 NGC 6530 288

NGC 6530 288 is also located in field 2 of our observations. With V = 13.23 mag and (B-V) = 0.41 mag, NGC 6530 288 (WEBDA #28) falls into the region of the classical instability strip in the HR diagram. Hence, it was also one of the prime candidates to search for PMS pulsation. In the V filter data, a frequency of 17.996 c/d corresponding to a period of 1.33 h with an amplitude of 0.8 mmag is significant. The same frequency can be found in the B filter data, but it is not above 4$\cdot$S/N and has a S IGS PEC significance of only 5.04. Longer time-series are needed to confirm the variability of this star.

7 Pulsation constants

The main objective of this work was to discover new members of the group of pulsating pre-main sequence stars, not to provide a complete asteroseismic investigation. Additional photometric time-series from multi-site campaigns or from space are needed to characterize the stars asteroseismologically in more detail.

Table 5: Pulsation constant Q for the newly discovered PMS pulsators (see text for detailed explanation).

Nevertheless, for classical $\delta$ Scuti stars, the pulsation constant, Q, can be used to distinguish radial fundamental or higher overtone modes. As shown by Breger (1979), Q can be computed using the following relation:

\begin{displaymath}%
{\rm log}~Q = -6.454 + {\rm log}~P + 0.5~{\rm log}~g + 0.1~M_{{\rm bol}} + {\rm log}~T_{\rm eff}.
\end{displaymath} (3)

Using the empirical calibration by Reed (1998), which is described below, it is possible to derive $M_{\rm bol}$, $\log$ $T_{\rm eff}$, and log L/$L_{\odot}$ from (B-V)0 and Mv for all the newly discovered pulsating pre-main sequence stars. For stars with $\log$ $T_{\rm eff}$ $\le$ 3.961 (i.e., $T_{\rm eff}$ $\le$ 9140 K) the relation

\begin{displaymath}%
(B-V)_0 = -3.648~(\log T_{\rm eff}) + 14.551
\end{displaymath} (4)

can be used to determine $\log T_{\rm eff}$. The bolometric magnitude, $M_{\rm bol}$, is further defined as:

\begin{displaymath}%
M_{{\rm bol}} = M_v + BC(T),
\end{displaymath} (5)

where the bolometric correction, BC(T), can be determined from
                                   BC(T) = $\displaystyle -8.499~(\log T_{\rm eff}-4)^4 + 13.421~(\log T_{\rm eff}-4)^3$  
    $\displaystyle - 8.131~(\log T_{\rm eff}-4)^2 - 3.901~(\log T_{\rm eff}-4) - 0.438.$ (6)

To compute the pulsation constant, Q, the mass ratios and $\log g$ values are needed. Voigt (1991), e.g., gives:

\begin{displaymath}%
{\rm log} (M/M_{\odot}) = 0.59 - 0.13~M_{\rm bol}
\end{displaymath} (7)


\begin{displaymath}%
g/g_{\odot} = (M/M_{\odot})\cdot(T_{\rm eff}^4 / T_{\rm eff,\odot}^4) / (L/L_{\odot}),
\end{displaymath} (8)

with $T_{\rm eff,\odot}$ = 5780 K and $g_{\odot} = 2.7 \times 10^4~{\rm cm~s}^{\rm {-1}}$.

The radial Q values are 0.033 d for the fundamental mode, 0.025 d for the first, 0.020 d for the second, and 0.017 d for the third overtones; values smaller than 0.017 d are identified as higher overtones. Regarding the resulting Q values given in Table 5, it becomes evident that most of the detected periods correspond to higher overtone modes. Only the first period of NGC 6530 278 (P = 0.139 d) and possibly the fourth period of NGC 6530 85 (P = 0.094 d) seem to be radial fundamental modes.

8 Other variable stars

Several other cluster members show variability on very different timescales. Some other objects are variable, but may not be cluster members, or their variability remains inconclusive in the data. For completeness, the most interesting other variable stars, which are not $\delta$ Scuti like PMS pulsators, are discussed below in detail. All definitive and suspected variables are listed in the corresponding tables (Tables 6 and 7).

8.1 Other variable stars in IC 4996

8.1.1 IC 4996 67

IC 4996 67 has V = 12.82 mag and (B-V) = 0.52 mag (Delgado et al. 1998), but no spectral classification exists for this star. The light curves in both filters look variable with relatively high amplitudes (Fig. 14). The frequency analyses using P ERIOD98 and S IGS PEC yielded two significant frequencies, which seem to be intrinsic. The first frequency is at 1.928 c/d with amplitudes of 10.3 mmag in V and 10.5 mmag in B, and the second at 2.088 c/d with amplitudes of 9.7 mmag in V and 10.2 mmag in B. Permitting the star to be less distant than the cluster, it could be a foreground $\gamma$ Doradus type star. Further studies are needed to determine the nature of its variability, and a reliable spectral classification would be most important.

Table 6: Other variables and suspected variables in the field of IC 4996: "Star'' denotes our star number and "Ref'' the cross-reference with the literature (D ... Delgado et al. 1998, W ... number given in the WEBDA database). Values marked with asterisks were derived from our calibration, when no values were available in the literature.

Table 7: Other variables and suspected variables in the field of NGC 6530: "Star'' denotes our star number and "Ref'' the cross reference with numbers given in the WEBDA database, V and (B-V) values were taken from the literature.

8.1.2 IC 4996 80

IC 4996 80 is located at the blue edge of the instability region in the HR diagram. V = 14.04 mag, (B-V) = 0.58 mag, and (U-B) = 0.15 mag were derived by Delgado et al. (1998), but again no spectral classification is available. In the data sets of both filters, one significant frequency at 3.578 c/d, i.e., a period of $\sim $6.7 h, with an amplitude of 2.0 mmag in V and 2.9 mmag in B was detected (Fig. 15).

If star 80 is slightly hotter than indicated from its color, it could fall into the region of the slowly pulsating B star (SPB-) instability domain in the HR diagram.


  \begin{figure}
\par\includegraphics[width=7cm,clip]{5449fi14.eps}\end{figure} Figure 14: Differential light curve of IC 4996 67 in V ( top) and B filters ( bottom, shifted for better visibility). The solid line corresponds to the multi-sine fit with the two significant frequencies.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=6.3cm,clip]{5449fi15.eps}\end{figure} Figure 15: Phase plots for IC 4996 80 in V ( top) and B ( bottom) filters with a frequency of 3.578 c/d.
Open with DEXTER

8.1.3 IC 4996 86

IC 4996 86 is again bluer than the blue edge of the classical instability region and has V = 13.87 mag, (B-V) = 0.47 mag, and (U-B) = 0.08 mag (Delgado et al. 1998). Two frequencies at 3.294 c/d (i.e., a period of 7.3 h) and 3.380 c/d (i.e., a period of 7.1 h) indicate a possible $\beta$ Cephei- or SPB-type variability.

8.1.4 IC 4996 95

For IC 4996 95, no magnitudes and colors are available from the literature. Our calibration yields V = 13.72 mag and (B-V) = 0.61 mag, which puts the star outside the classical instability strip. The two significant frequencies at 2.807 c/d (with 5.4 mmag amplitude in V and B) and 2.942 c/d (with 6.9 mmag amplitude in V and 5.1 mmag in B) lie very close to each other and produce a beat in the light curve (Fig. 16). The corresponding periods are 8.6 and 8.2 h, respectively. It also seems that the star is not a cluster member, and the origin of its variation can only be suspected. It may be caused by $\beta$ Cephei- or SPB-type pulsations.


  \begin{figure}
\par\includegraphics[width=7.25cm,clip]{5449fi16.eps}\end{figure} Figure 16: Differential light curve of IC 4996 95 in V with a multi-sine fit (solid line) of the two significant frequencies.
Open with DEXTER

8.1.5 IC 4996 104

With V = 14.03 mag, (B-V) = 0.64 mag, and (U-B) = 0.02 mag the star is also near the blue edge of the classical instability region. Two frequencies at 2.209 c/d (i.e., a period of 10.9 h) with an amplitude of 7.5 mmag in V and 10.4 mmag in B, and at 1.786 c/d (i.e., a period of 13.4 h) with an amplitude of 6.6 mmag in V and 7.7 mmag in B were found to be significant, and a simultaneous multi-sine fit to the data reproduces the shape of the light curve well (Fig. 17). There is no spectral class available, but it could be a SPB star.


  \begin{figure}
\par\includegraphics[width=7cm,clip]{5449fi17.eps}\end{figure} Figure 17: Differential light curves for IC 4996 104 in V ( top) and B ( bottom, shifted for better visibility) filters; the solid line corresponds to a multi-sine fit with the two significant frequencies.
Open with DEXTER

8.2 Other variable stars in NGC 6530

8.2.1 NGC 6530 13

NGC 6530 13 (V = 12.65 mag) was located in the overlapping region of the CCDs. It was classified as a B8 star by van den Ancker et al. (1997). Van Altena & Jones (1972) estimate its cluster membership probability to be only 25%. The frequency analysis yielded a single significant frequency at 3.432 c/d, corresponding to a period of 7.0 h, with amplitudes of 3.6 mmag in V (Fig. 18) and 4.7 mmag in B. The star could be a SPB-type variable.


  \begin{figure}
\par\includegraphics[width=6.9cm,clip]{5449fi18.eps}\end{figure} Figure 18: Amplitude spectrum of NGC 6530 13 in the V filter, where the significant frequency is indicated with an arrow.
Open with DEXTER

8.2.2 NGC 6530 21

NGC 6530 21 was observed in the overlapping field of the CCDs. With V = 14.17 mag and (B-V) = 0.97 mag, it is located redwards of the instability region in the HR diagram. No spectral type or membership information is available for this star, but from its (B-V)0 = 0.62 mag (dereddened with E(B-V) = 0.35 mag, adopted from Sung et al. 2000) it should be an early G spectral type with log  $T_{\rm eff}$ = 3.77. Fourier analysis with P ERIOD98 yields two significant frequencies in the data sets of both filters, which are confirmed by S IGS PEC. A simultaneous multi-sine fit to the data with the two significant frequencies detected with both methods, at 0.277 c/d (i.e., a period of 3.613 days) with amplitudes of 8.7 mmag in V and 9.7 mmag in B and at 0.166  c/d (i.e., a period of 6.037 days) with amplitudes of 2.8 mmag in V and 2.3 mmag in B reproduces the shape of the light curve well (Fig. 19).


  \begin{figure}
\par\includegraphics[width=6.55cm,clip]{5449fi19.eps}\end{figure} Figure 19: Part of the light curve of NGC 6530 21 in the V ( top) and B filters ( bottom).
Open with DEXTER

8.2.3 NGC 6530 55

NGC 6530 55 is situated bluewards of the instability region and was classified as spectral type B6 by van den Ancker et al. (1997). The shape of the light curve clearly indicates variability (see Fig. 20). The frequency analyses using P ERIOD98 and S IGS PEC detect a single significant frequency of 2.285 c/d corresponding to a period of 10.5 h. Such a variability could be explained by SPB type pulsation.


  \begin{figure}
\par\includegraphics[width=6.85cm,clip]{5449fi20.eps}\end{figure} Figure 20: Part of the differential light curve of NGC 6530 55 in the V ( top) and B filters ( bottom, shifted for better visibility).
Open with DEXTER

8.2.4 NGC 6530 73

NGC 6530 73 (V = 11.94 mag) was classified as spectral type B5e by van den Ancker et al. (1997). From the shape of its light curve, it is most likely a binary star: as the eclipse is deeper in the V than in the B filter, the redder star seems to occult its bluer companion (see Fig. 21).


  \begin{figure}
\par\includegraphics[width=6.55cm,clip]{5449fi21.eps}\end{figure} Figure 21: Light curve of the suspected eclipsing binary NGC 6530 73 in the V (open circles) and B (solid circles) filters.
Open with DEXTER

9 Cluster properties

9.1 IC 4996

40 stars lie within the region of the classical instability strip in IC 4996, and hence have been prime candidates to search for PMS pulsation. Two stars, IC 4996 37 and IC 4996 40, had large enough amplitudes to allow for an unambiguous detection of pulsation, while for IC 4996 46 pulsation can only be suspected. From the location of the PMS pulsators in the HR diagram, they are most likely members of the cluster (Fig. 22). A definite decision on their cluster membership can only be drawn after a detailed proper motion or radial velocity study. However, the fraction of pulsating pre-main sequence stars corresponds to $\sim $10%. As the observations have been performed in service mode and the integration time per frame was not long enough for stars in the magnitude range between 14 < V < 16 mag, additional PMS pulsators may be discovered in a follow-up observing run for data with a lower noise level.

Of the total 113 stars analyzed in the field of the cluster, 16 stars have been found to be variable, including the PMS pulsators, which means only $\sim $15%. The reason for such a low percentage lies in the higher noise level of the measurements: Especially towards the fainter, and hence redder stars, the noise level increased dramatically and made it impossible to detect variability. Figure 22 shows the location of all detected variables and suspected variable stars in IC 4996, where the lack of variable, red stars is clearly illustrated.


  \begin{figure}
\par\includegraphics[width=6.55cm,clip]{5449fi22.eps}\end{figure} Figure 22: Observational HR diagram of all stars in the field of IC 4996: bona fide PMS pulsators (filled circles), candidate PMS pulsator (open circle), other variables (filled triangles), and other suspected variables (open triangles); the ZAMS values are taken from Schmidt-Kaler (solid line), and the borders of the classical instability strip (dashed lines) have been transformed into the MV-(B-V)0 plane.
Open with DEXTER

Figure 23 shows HR diagrams for IC 4996 (top: V and bottom: B filters), where the stars have different symbols and sizes according to their amplitude noise levels ( $amp_{\rm noise}$) based on their time-domain point-to-point scatter ( $\sigma_{\rm pt2pt}$). The amplitude noise level in the Fourier domain is given as:

\begin{displaymath}%
amp_{{\rm noise}} = 2 \cdot \sigma_{{\rm pt2pt}} / \sqrt{N}
\end{displaymath} (9)

where N denotes the number of data points per star. The brightest star, IC 4996 34, has a V = 10.45 mag and was partly saturating the CCD chip, which is the reason why it has one of the highest values for the amplitude noise.


  \begin{figure}
\par\includegraphics[width=6.65cm,clip]{5449fi23.eps}\vspace*{3mm}
\includegraphics[width=6.65cm,clip]{5449fi24.eps}\end{figure} Figure 23: HR diagram of all measured stars in the field of IC 4996, where the stars have different symbols and sizes according to their expected amplitude noise levels; top: V filter, bottom: B filter; the values for the ZAMS (solid line) are taken from Schmidt-Kaler (1965).
Open with DEXTER

9.2 NGC 6530

In NGC 6530, six bona fide $\delta$ Scuti-like pulsating pre-main sequence stars could be discovered. Of those pulsators, one star is located in the overlapping region, two in field 1, and three in field 2, which justified the decision to split the observing time for two different fields. For NGC 6530 288, which is situated in field 2 of the observations, the data quality is not sufficient to unambiguously decide on its pulsation. Hence, it is only a likely candidate for a pulsating PMS star. Taking into account that 25 stars have been primary candidates for pulsation because they are located in the region of the classical instability strip, the six pulsating stars correspond to 24%. This percentage is higher than for the other clusters we analyzed. It is also closer to the fractions of 1/3 to 1/2 pulsating stars in the lower instability strip reported by Breger (2000) for the classical, post-, and main-sequence $\delta$ Scuti stars. Whether this discrepancy is of instrumental origin or has an astrophysical background connected with the evolution of such young stars, cannot yet be decided.

Among the 194 observed stars in the fields of NGC 6530, at least 12 other stars were found to be variable. Figure 24 shows the location of all detected variable and suspected variable stars in NGC 6530 in the HR diagram, and it becomes evident that some stars in the field of the cluster might not be members, but are fore- or background field stars.


  \begin{figure}
\par\includegraphics[width=6.6cm,clip]{5449fi25.eps}\end{figure} Figure 24: Observational HR diagram of all stars in the field of NGC 6530: bona fide PMS pulsators (filled circles), candidate PMS pulsator (open circle), other variables (filled triangles), and other suspected variables (open triangles); the ZAMS values are taken from Schmidt-Kaler (solid line) and the borders of the classical instability strip (dashed lines) have been transformed into the MV-(B-V)0 plane.
Open with DEXTER

Analogous to IC 4996, the amplitude noise level based on the point-to-point scatter in the time domain has been calculated. Figure 25 shows the HR diagrams (top: V and bottom: B filters), where the stars have different symbols and sizes according to their amplitude noise level bins. Please note that again some blue stars are too bright and were partly saturated on the CCD chip. Hence, their number of data points is somewhat lower and their accuracy is worse than for their slightly fainter counterparts.


  \begin{figure}
\par\includegraphics[width=6.65cm,clip]{5449fi26.eps}\vspace*{3mm}
\includegraphics[width=6.65cm,clip]{5449fi27.eps}\end{figure} Figure 25: HR diagram of all measured stars in the field of NGC 6530, where the stars have different symbols and sizes according to their expected amplitude noise level bins; top: V filter, bottom: B filter; the values for the ZAMS (solid line) are taken from Schmidt-Kaler (1965).
Open with DEXTER

10 Conclusions

CCD photometric time-series of 113 stars in IC 4996 and 194 stars in NGC 6530 have been analyzed with the aim of searching for pulsating pre-main sequence members of the clusters using classical Fourier techniques as well as the new concept S IGS PEC (Reegen 2006).

40 stars are situated in the instability region of the HR diagram in IC 4996, and hence have been primary candidates in which to search for pulsations. Only two cluster members, IC 4996 37 and IC 4996 40, have been identified as bona fide PMS pulsators, whereas pulsation can only be suspected for IC 4996 46. According to the mid-A spectral types of the two bona fide pulsating PMS stars and their location in the HR diagram, the probability of their membership in the cluster is very high. Those three stars correspond to only $\la$8% pulsators within the region of the classical instability strip in IC 4996 down to an amplitude noise level of 0.2 mmag in V and 0.6 mmag in B in the Fourier domain. The percentage of detected pulsating PMS stars in the instability region is somewhat lower than is observed for their post-ZAMS counterparts, even if the not optimal data quality is taken into account. However, this finding is still consistent with earlier investigations of pulsating PMS stars in young clusters (e.g., Zwintz et al. 2005).

In NGC 6530 the total number of potential candidates for PMS pulsation of 30 is less than for IC 4996. Still, significantly more pulsating pre-main sequence stars have been discovered in NGC 6530. Six stars (NGC 6530 5, NGC 6530 82, NGC 6530 85, NGC 6530 263, NGC 6530 278, and NGC 6530 281) definitely show $\delta$ Scuti-like pulsation periods, while for one star (NGC 6530 288) pulsation can only be suspected. This corresponds to $\sim $24% pulsating stars in the region of the classical instability strip down to an amplitude noise level of 0.2 mmag in both filters in the Fourier domain. Three of the PMS pulsators in NGC 6530 have been included in the proper motion study by van Altena & Jones (1972), who report a high probability of cluster membership. For the other stars, the membership was estimated from their position in the HR diagram and from their location on the sky with respect to the center of the cluster.

Often not much additional information is available for a star, e.g., for most of the objects spectral classification would urgently be needed. For the couple of pulsating stars with known spectral types, mostly emission lines have been reported in the literature, and very often an excess in the infrared was also observed. From stellar evolution theory it is suspected that such young stars are fast rotators, but for nearly all known PMS pulsators no $v\cdot\sin{i}$ determination has been performed yet, making it difficult to confirm this concept. Additionally, the detected periodicities due to pulsation always lie on top of longer, sometimes irregular light variations caused by circumstellar material. All these features are characteristics of the Herbig Ae/Be stars indicating their pre-main sequence nature.

From the pulsational analysis of the 8 pulsating pre-main sequence stars discovered in the two clusters, it is evident that the higher overtone modes seem to have higher amplitudes in pre-main sequence stars and hence are easier to detect. This is confirmed by the asteroseismic investigations of other detected PMS pulsators (e.g. Ripepi et al. 2006). At the moment, the reason for this behavior can only be speculated about: the amplitudes of fundamental and lower overtone modes seem to be damped by convection, which favor higher overtones.

A total of 25 PMS $\delta$ Scuti-like stars and 10 PMS $\delta$ Scuti-like candidates are known in the literature, including the 10 bona fide and 2 candidate PMS pulsators found in this study. Such a number will allow to probe the instability strip for pre-main sequence stars empirically (Zwintz 2006) for the first time, but also supports the ongoing development of non-radial pulsation models for pre-main sequence stars.

Acknowledgements
This project was supported by the Austrian Fonds zur Förderung der wissenschaftlichen Forschung (P14984). Frequency analysis was performed using P ERIOD98 written by M. Sperl (1998) and S IGS PEC developed by P. Reegen (2006). Use was made of the WEBDA database, developed by J.-C.Mermilliod (Laboratory of Astrophysics of the EPFL, Switzerland) and operated by E. Paunzen at the University of Vienna, Austria. Special thanks goes to Rafa Garrido and the team of the Sierra Nevada Observatory for carrying out the observations of IC 4996. Finally, it is a pleasure to acknowledge T. Kallinger, A. Pamyatnykh, and P. Reegen for fruitful discussions on data reduction and analysis.

References

 

Copyright ESO 2006