A&A 463, 589-600 (2007)
DOI: 10.1051/0004-6361:20066242
P. Székely1,2 - L. L. Kiss3 - R. Jackson3 - A. Derekas3 - B. Csák4 - K. Szatmáry1
1 - Department of Experimental Physics & Astronomical Observatory, University
of Szeged, 6720 Szeged, Dóm tér 9., Hungary
2 -
Hungarian Eötvös Fellowship, School of Physics, University of Sydney, NSW 2006,
Australia
3 - School of Physics, University of Sydney 2006, NSW 2006, Australia
4 - Department of Optics and Quantum Electronics, University of Szeged, 6720 Szeged,
Dóm tér 9., Hungary
Received 13 August 2006 / Accepted 14 November 2006
Abstract
Aims. NGC 362 is a bright southern globular cluster in the foreground of the Small Magellanic Cloud (SMC), for which no extensive variability survey has ever been done. Finding regularly pulsating RR Lyrae stars in the cluster can lead to improved metallicity and distance estimates of the system, while other types of variable objects may be used to confirm the results.
Methods. Time-series CCD photometric observations have been obtained. Light curves have been derived with both profile fitting photometry and image subtraction. We developed a simple method to convert flux phase curves to magnitudes, which allows the use of empirical light curve shape vs. physical parameters calibrations. Periods and light curve parameters of the detected variable stars have been determined with Fourier analysis, phase dispersion minimization and string-length minimization. Using the RR Lyrae metallicity and luminosity calibrations, we have determined the relative iron abundances and absolute magnitudes of the stars. The color-magnitude diagram has been fitted with Yale-Yonsei isochrones to determine reddening and distance independently. For five RR Lyrae stars we obtained radial velocity measurements from optical spectra.
Results. We found 45 RR Lyr stars, of which the majority are new discoveries. While most of them are cluster members, as shown by their radial velocities and positions in the color-magnitude diagram, we also see a few stars in the galactic field and in the outskirts of the SMC. About half of the RR Lyraes exhibit light curve changes (Blazhko effect). The RR Lyrae-based metallicity of the cluster is
,
the mean absolute magnitude of the RR Lyrae stars is
mag implying a distance of
kpc. The mean period of RRab stars is
days. These properties place NGC 362 among the Oosterhoff type I globular clusters. The isochrone fit implies a slightly larger distance of
kpc and an age of
Gyr. We also found 11 eclipsing binaries, 14 pulsating stars of other types, including classical Cepheids in the SMC and 15 variable stars with no firm classification.
Conclusions. NGC 362 hosts a large number of RR Lyrae stars, which makes the cluster a potentially important test object for studying the Blazhko effect in a chemically homogeneous environment.
Key words: stars: variables: general - stars: variables: RR Lyr - globular clusters: individual: NGC 362
Although globular clusters (GCs) play an important role in testing stellar evolutionary models, one can still find suprisingly bright neglected GCs, especially in the southern hemisphere. NGC 362 is such a cluster, having been target of relatively few studies and so far lacking any comprehensive search for variable stars, particularly in the central regions of the cluster. Apart from the pioneering work of Sawyer (1931), who presented data for fourteen variable stars, recent works mainly aimed at obtaining deep color-magnitude diagrams (CMDs, e.g. Alcaino 1976; Harris 1982; Bolte 1986, 1987; Bellazzini et al. 2001). In the catalog of globular cluster variables (Clement et al. 2001, Clement 2002) there are only 16 stars from NGC 362 and many are far from the cluster center and have ambiguous properties.
We have been carrying out a CCD photometric survey of southern GCs since
mid-2003. In this paper we discuss the results for NGC 362, which is located in front
of the outer edge of the Small Magellanic Cloud (SMC). NGC 362 lies at the position
,
(J2000.0),
,
,
located at 8.5 kpc from the Sun and
9.4 kpc
from the galactic center; the metallicity ([Fe/H]) is around -1.16, while
the Horizontal Branch is at
.
The reddening is almost negligible,
E(B-V)=0.05 (Harris 1996).
Models for globular clusters suggest that all the member stars were formed from the same gas cloud at approximately the same time. Physical parameters of member stars should therefore be representative of the cluster itself. Recently, many empirical relations have been calibrated for RR Lyrae stars (Jurcsik & Kovács 1996; Kovács & Jurcsik #Kov&, Jurcsik 1998, Kovács & Walker #Kov&), which can be used to determine metallicities, absolute magnitudes, reddenings and other physical parameters from the light curve shapes of these variables. Here we present an analysis of RR Lyrae stars in NGC 362 and briefly discuss properties of variable stars of other types.
Table 1: Log of time series observations.
Table 2: Log of spectroscopic observations.
Time-series V-band CCD photometric observations were carried out on 17 nights
between July 2003 and October 2004, using the 1 m ANU telescope at Siding Spring
Observatory, Australia. For imaging we used three chips of the Wide Field Imager
giving
field of view. In addition to the photometry, we
measured radial velocities for 5 RR Lyrae stars from optical spectra taken with the
2.3 m ANU telescope and the Double-Beam Spectrograph. We used only the red arm of the
spectrograph with the 1200 mm-1 grating, covering the wavelength range
5800-6800 Å. The nominal resolving power was
.
With 20-min exposures, the spectra had S/N ratios of 20-50, depending on the
weather conditions. The full journals of observations are summarized in Tables 1 and 2.
All data were processed with standard IRAF routines, including bias and flat corrections. Cluster images were then
analysed both with PSF photometry, using the daophot package in IRAF, and image
subtraction photometry, using the ISIS 2.1 package of Alard (2000). For standard
photometric calibrations we observed selected equatorial fields of Landolt (1992): SA
109-954, SA 110-503, SA 111-1965, Mark A and SA 113-241. The transformation
coefficients were determined with the task photcal. Radial velocities were measured
with cross-correlation, using spectra of the IAU radial velocity standard HD 187691 and
the task fxcor, leading to an estimated accuracy of
7 km s-1.
Light curves from PSF-fitting and image subtraction were analysed for finding periods. For this, we applied a combination of Fourier analysis, phase dispersion minimization (Stellingwerf 1978) and string-length minimization (Lafler & Kinman 1965). The latter gives better results when the light curve contains sharp features, such as narrow minima of Algol-type eclipsing binaries. Every light curve was then inspected visually and classified according to its period, amplitude and shape. Variable star designations were assigned in an increasing order of right ascension.
Image subtraction is far more efficient than PSF photometry in detecting variable stars in crowded fields. Globular cluster studies are heavily affected by this, where large incompleteness factors can arise from the highly concentrated distributions of stars (Kaluzny et al. 2001). Even for the best-studied clusters, image subtraction can yield new discoveries (e.g. Clementini et al. 2004). There is, however, a major practical difficulty in using image subtraction results for the empirical light curve shape vs. physical parameters calibrations. Whereas ISIS gives time variations of brightness expressed in differential fluxes (relative to a master CCD frame), all the empirical relations were calibrated with light curves in magnitudes, which have a different shape, especially for the largest-amplitude objects. To convert fluxes to magnitudes, one has to determine a reliable magnitude estimate for every variable star in the master CCD frame, otherwise the conversion (and employing the calibrations) is not possible (see, e.g., Corwin et al. 2006).
Since we detected about twice as many RR Lyrae stars with ISIS 2.1 than with daophot, we developed a simple method for a flux-to-magnitude conversion. The main
assumption is that there is a unique scaling between the mean-subtracted ISIS fluxes and
flux values calculated from the mean-subtracted daophot magnitudes, so that one can
write the following equation:
![]() |
(1) |
![]() |
Figure 1: The conversion between mean-subtracted flux values. |
Open with DEXTER |
![]() |
Figure 2: Mean-subtracted magnitude curves obtained using daophot (light grey/turqoise points) and Eq. (1) applied to ISIS fluxes (black dots). The large scatter in V29 and V27 is caused by crowding. |
Open with DEXTER |
In Fig. 1 we plot data for two high-amplitude RR Lyrae stars which are
located in the outer regions of the field. V71 is the largest-amplitude variable with
no crowding, so its data were used to determine c. The solid dots in Fig. 1 yield a slope of
.
Since these
calibration points cover the whole range of ISIS fluxes for all the RR Lyraes, no
extrapolation was needed for the other stars. For comparison, we also show data for
another uncrowded star, V35. The two sets do not differ by more than 1% anywhere, so
that Eq. (1), calibrated by V71 alone, can predict magnitudes from fluxes within 0.01
mag (without the zeropoint). For a few other RR Lyrae stars we also compared the
magnitude curves that were calculated from the ISIS flux data to the daophot
magnitudes and always found excellent agreement (Fig. 2), thus
confirming the applicability of our approach. We therefore converted ISIS flux
curves to magnitudes for every star with good phase coverage.
All observations presented in this paper are available electronically at the CDS.
![]() |
Figure 3: Radial velocity curves for two RR Lyrae stars. |
Open with DEXTER |
![]() |
Figure 4: Sample phased light curves of RR Lyr stars, in order of increasing pulsational period from V26 (P=0.333 d) to V57 (P=0.679 d). |
Open with DEXTER |
![]() |
Figure 5: Sample light curves of cluster RR Lyr stars. |
Open with DEXTER |
From the periods and light curve shapes we identified 45 RRab/RRc-type variable stars. In
some cases only partial phase coverage was obtained and only the RR Lyrae-class could be
safely concluded. The majority have very similar mean apparent brightnesses
(
mag). Five of these have been observed spectroscopically (see Fig. 3 for two examples) and their mean radial velocities (
220 km s-1)
confirm membership, since the cluster's radial velocity is
km s-1(Fischer et al. 1993). Figures 4 and 5 show sample phased/unphased
data for cluster member RR Lyraes with well-defined mean light curve shapes and with
partial coverage, respectively. Table 3 contains the fundamental photometric
data for all the detected RR Lyrae stars. Since photometric standards were observed only
on one night, V brightnesses and V-I colors are subject to uncertainty due to
photometric variability. To address this, we calculated the pulsation phase of the
standard V observation of the cluster for each RR Lyrae star, so that we know from the
V-band light curve the difference between the magnitude in that phase and the mean. The
measured standard V magnitudes were then corrected for these differences. With no
I-band light curve a similar correction for the standard I magnitudes was not
possible, so that the V-I indices may suffer from larger uncertainties (up to several
tenths of a magnitude). Stars with missing V and V-I in Table 3 values were
not detected as individual objects in the master standard frames, which means they were
merged with other stars. Variables in the core of the cluster were checked for blending
using an archived WFPC2 image of the Hubble Space Telescope (taken for the project "Hot
stellar populations in globular cluster cores'', HST Proposal 6095 by S. Djorgovski, 160 s exposure in the F439W filter). In the last column of Table 3 we give comments
and identifications of previously known variables from the literature.
In Fig. 6 we plotted the positions in the CMD of RR Lyrae stars
with reliable
calibrated photometry. Most of the stars fall on the
Horizontal Branch (HB) clearly establishing membership. There are
also some stars with deviant positions, for example, V80 and V81 at
,
both in the background of the cluster with distances well beyond 20 kpc. Further
RR Lyraes are seen between V=19-20 mag, which are located in the SMC.
![]() |
Figure 6: RR Lyrae stars in the color-magnitude diagram of NGC 362. |
Open with DEXTER |
![]() |
Figure 7: Upper panel: period distribution of cluster RR Lyrae stars. The darker columns refer to the RRc stars. Lower panel: Bailey-diagram for the cluster RR Lyrae stars. Solid lines: Clement & Rowe (2000); dashed line: Cacciari et al. (2005). |
Open with DEXTER |
We could determine periods, amplitudes and epochs for 28 RR Lyraes, of which 4 are
first-overtone (RRc) stars. This fairly low percentage implies that the cluster is
RRc-poor (cf. Fig. 1 of Castellani et al. 2003), which is characteristic for
Oosterhoff type I globular clusters (Oosterhoff 1939). The computed specific frequency of
RR Lyr variable stars, including the newly found ones, is
.
The histogram of the periods (upper panel in
Fig. 7) shows that large number of RRab stars have periods between 0.5 and 0.6 d. The mean period of 24 RRab stars is
or
,
while the period-amplitude relation
closely follows the Oosterhoff I lines given by Clement & Rowe (2000) and Cacciari
et al. (2005). Both properties are consistent with the Oosterhoff I class (see the
lower panel in Fig. 7).
In order to determine metallicities and absolute magnitudes of the RR Lyrae stars, we
fitted the phase diagrams with Fourier polinomials of up to the 10th order. Using Eq. (3)
of Jurcsik and Kovács (1996):
,
we estimated iron
abundances for 13 RR Lyrae stars. These objects have phase curves with low scatter that
covered both maxima and minima and did not show light curve modulations. We then
determined absolute magnitudes using the metallicity-luminosity relation of Kovács &
Jurcsik (1996):
.
Upon comparing this to our phase corrected
apparent magnitudes, we calculated distance moduli for the individual stars. We
determined other physical parameters (e.g. effective temperature, luminosity) using
further equations calibrated by Morgan et al. (2005), Kovács (#Kov&) and Jurcsik (1998).
The results are presented in Tables 4 and 5.
The numbers clearly show that most of the studied RR Lyraes define reasonably small ranges of the parameters. The distance to V58 is erroneous because the star is a blend, but its other, distance-independent, parameters are similar to those of the other variables. The calculated mean metallicity of the cluster is in good agreement with values found in the literature, e.g. Caldwell (1988) and Davidge (2000), while the distance is slightly smaller, but still compatible with, the 8.5 kpc in Harris (1996). In conclusion, RR Lyraes in NGC 362 define a consistent set of physical parameters for the cluster.
Table 3:
RR Lyrae stars in the field of NGC 362. The given coordinates are accurate
within -
.
Remarks on types: 1. strong Blazhko
modulation, the whole light curve is changing; 2. slight modulation in certain phases;
3. blend; 4. unsure. The last column
contains other designations of already known variables: SAW xxx -- Sawyer (1931);
SMC Vxxx -- Samus et al. (2004); MA 1315 --
Meyssonnier & Azzopardi (1993); SSP xx and HV xx --
Sharpee et al. (2002).
Table 4:
Fourier and the physical parameters of selected RRab stars. Notes:
Blazhko effect present. x Blend. * Outliers, not included in the cluster average.
Table 5: Fourier and physical parameters of selected RRc stars.
![]() |
Figure 8: RR Lyr stars with Blazhko effect. From bottom to top: V7, V67, V68. |
Open with DEXTER |
It is known that a large fraction of RR Lyr stars (20-30% of RRabs and 2% of RRcs, Kovács #Kov&) show periodic amplitude and/or phase modulations, the so-called Blazhko effect, which is one of the greatest mysteries in classical pulsating star research. In our sample, three variables, V7, V67 and V68, have strong light curve modulations, while a further 14 stars (V15, V23, V26, V27, V29, V35, V39, V40, V43, V55, V58, V64, V78, V80) exhibited subtler changes during the 15 months of observations. Furthermore, V33, V34 and V50 also showed detectable variations but these objects are in the very crowded core of the cluster and their photometry is of lower quality.
Assuming that these light curve modulations are representative for the whole sample, the occurrence rate is rather high: about half of the cluster RR Lyrae stars show some kind of changes in the phased data over the course of our monitoring. The three stars with the strongest modulations are shown in Fig. 8, where one can see that changes in the light curve shape occurred in almost all phases. Very recently, Jurcsik and co-workers obtained high-quality observations for a number of field RR Lyraes with the Blazhko effect (Jurcsik et al. 2006; Sódor et al. #S&; Jurcsik et al. 2005), revealing many puzzling features of the observed modulations, such as the confinement of the light curve changes to a certain range of phases or the presence of doubly-periodic modulations. Our data are not extensive enough to study these phenomena, but it is obvious that there is a valuable set of RR Lyraes in the cluster. With no firm theoretical explanation of the effect, a globular cluster with many modulated RR Lyraes clearly deserves further investigations. Based on our observations, NGC 362 may be an excellent target for obtaining long-term, homogeneous photometry of variables showing the Blazhko effect.
![]() |
Figure 9: V, V-I CMDs with different star samples within increasing apertures centered on NGC 362. The aperture radii were 3, 5, 10 and 30 arcmin. The outer regions clearly show a secondary faint RGB/AGB caused by stars in the SMC. The isochrones were fitted to data in the upper right (R<5') panel. |
Open with DEXTER |
Figure 9 shows CMDs for stars within an increasing
area around the cluster. The radii were chosen as 3,
5
,
10
and 30
,
leading to 2450, 4060, 5870 and 8567 stars in each plot.
As we increase the radius, the CMD exhibits an increasing field contamination and the
well-defined faint secondary Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB)
of the SMC. We selected the second diagram
(4060 stars) for isochrone fitting.
![]() |
Figure 10:
V, V-I CMD of the cluster with fitted Z=0.001 isochrones for 10, 11 and
12 Gyr (from top). The derived cluster parameters are
![]() |
Open with DEXTER |
To determine the distance modulus, age, metallicity and reddening of the cluster we
fitted Yonsei-Yale (Y2) isochrones (Demarque et al. 2004).
The best fit was given by the 11 Gyr isochrone with Z=0.001,
Y=0.232,
and
.
It yields a distance modulus
of
and reddening of
E(V-I)=0.05. The two next best fits are also shown in
Fig. 10, which gave very similar distance and reddening estimates
but have an age of 10 and 12 Gyr, respectively. These two fits bracket the 11 Gyr
isochrone in Fig. 10.
Between the HB and the RGB we also detect the Red Clump (RC) as a distinct concentration
of stars at
mag and
.
From the intrinsic color of the RC
(V-I=0.92, Olsen & Salyk, 2002) and absolute brightness
(for local Red Clump
,
Alves et al. 2002) we can
independently estimate the reddening and distance to the cluster:
E(V-I)=0.04 and
,
both in good agreement with the results of the isochrone fitting.
We corrected the distance modulus for extinction using
E(V-I)/E(B-V)=1.28
(adopting
AI/E(B-V)=1.962 and
AV/E(B-V)=3.24, Savage & Mathis 1979;
Schlegel et al. 1998), thus
AI=0.08. In this way we arrived at
E(B-V)=0.04 and,
using
AI/AV=0.601, we obtained
AV=0.13, resulting in a true distance
modulus of
.
While the red HB is heavily populated, there is a noticeable lack of stars
in the blue wing of the HB (see a recent discussion of this feature by Bellazzini et al. 2001 and Catelan et al. 2001). From Fig. 10 we measured the
following relative age horizontal and vertical parameters (Rosenberg et al. 1999):
and
.
Both values agree within the uncertainties with those measured by Rosenberg
et al. (1999). The brightness of the horizontal branch is
,
which is consistent with the RR Lyr-based relation
relation (Gratton et al. 2003), which predicts
.
We also calculated the
Horizontal Branch morphology parameter L, as defined by Lee (1990):
L=(B-R)/(B+R+V),
where B, R and V are the numbers of blue, red, and variable HB stars,
respectively. The result is L=-0.855, which can be considered as an improvement to
the currently available value (-0.87, Harris 1996), thanks to our better statistics
of variable HB stars.
![]() |
Figure 11: RA-Dec positions of the variables stars projected on the positions of all the stars for which WCS transformation could be done. Blank sections are due to gaps between CCD segments. Note the increasing stellar density towards the SMC (in southern direction). |
Open with DEXTER |
Table 6: Other variable stars in the field of NGC 362. The last column contains other designations of already known variables: SAW xxx -- Sawyer (1931); SMC Vxxx -- Samus et al. (2004); SSP xx and HV xx -- Sharpee et al. (2002).
We found a reasonably large number of variable stars other than RR Lyraes in the field of NGC 362. In Fig. 11 we plot the celestial positions of 8616 stars for which transformation to World Coordinate System could be achieved, with all variable stars shown separately. Apparently, most of the cluster variables are located in the central regions. The excess of variables in the lower part of the figure is caused by the Small Magellanic Cloud, located just outside of the field of view towards the south.
Selected light curves are shown in Fig. 12, while the basic
parameters are listed in Table 6. Here the V-I colors are again only
indicative because of the single-epoch standard I observation. V brightnesses were
corrected for variability in the same way as for the RR Lyraes. Besides eclipsing
binaries, field
Scuti stars and long-period variables, the most interesting
objects include several Cepheids (three new discoveries) that can be identified as
members of the Small Magellanic Cloud. Their positions in the sky fall closer to the
SMC, which strengthens the identification.
![]() |
Figure 12:
Sample light curves of field variable stars. Top and second rows: Cepheids;
third row: eclipsing binaries, bottom row: ![]() |
Open with DEXTER |
In Fig. 13 we plotted the locations of non-RR Lyr variables in the
CMD. Blue squares show the Cepheids, which indeed occupy the
expected instability strip of the SMC. Moreover, their mean V magnitudes strictly
follow the Cepheid period-luminosity relation
,
based
on OGLE observations of the SMC (Udalski et al. 1999). In addition, we can clearly
identify four cluster member red giant variables on the AGB and RGB and
one in the SMC.
Table 7:
Frequencies for two Scuti/SX Phe stars.
Sct-type variables show a large range in apparent magnitude, showing that they
belong to the galactic field. Based on its colour, magnitude and proximity to
NGC 362, V63 may be member of the cluster thus belonging to the SX Phe class of variable
stars. Two of these stars, V63 and V69, exhibited complex multiperiodic variations,
from which we determined several pulsation frequencies with iterative sine-wave fitting
using Period04 of Lenz & Breger (2005). Although a thorough pulsation analysis of these
stars is beyond the scope of this paper, we list the resulting frequencies in Table 7. We included all frequencies with S/N ratio (Breger et al. 1993) higher than
4. For V63, we suspect that the low-frequency components can be caused by instrumental
effects and/or poor sampling. Otherwise, the frequency ratios of the strongest modes
suggest low-order radial pulsation, possibly flavoured with non-radial modes.
Finally, we found a number of variable stars, for which neither periods nor types could be determined. They are listed in Table 8.
Table 8: Variable stars with no period and secure classification.
![]() |
Figure 13: Variable stars of other types than RR Lyr in the color-magnitude diagram of NGC 362. |
Open with DEXTER |
The main results of this paper can be summarized as follows.
A brief overview of the main physical parameters of NGC 362 is given in Table 9. Of these, the only significant difference is between the distance
estimates using the RR Lyraes and the color-magnitude diagram. The difference is likely to
originate predominantly from the RR Lyraes, whose absolute magnitude-metallicity
relation is not yet settled, with a number of suggestions about nonlinear
calibrations of
[Fe/H] (see Sect. 5 of Sandage & Tamman 2006).
Table 9: Physical parameters of NGC 362.
Acknowledgements
This project has been supported by the Hungarian OTKA Grant #T042509, a Hungarian Eötvös Fellowship to PSz and the Australian Research Council. L.L.K. is supported by a University of Sydney Postdoctoral Research Fellowship. The NASA ADS Abstract Service was used to access references. This work used SIMBAD and Vizier databases operated at CDS-Strasbourg, France.
In this section we give additional information for some of the detected variable stars.