A&A 448, 1023-1030 (2006)
DOI: 10.1051/0004-6361:20054299
J. F. C. Santos Jr.1 - J. J.
Clariá2, - A. V. Ahumada2,
-
E. Bica3 - A. E. Piatti4 - M. C. Parisi2,
1 - Departamento de Física, ICEx, UFMG,
CP 702, 30123-970, Belo Horizonte, MG, Brazil
2 -
Observatorio
Astronómico, Laprida 854, 5000, Córdoba, Argentina
3 -
Depto. de Astronomia, UFRGS, CP 15051, 91500-970, Porto Alegre, Brazil
4 - Instituto de Astronomía y Física del Espacio, CC 67, Suc.
28, 1428, Buenos Aires, Argentina
Received 4 October 2005 / Accepted 15 November 2005
Abstract
Aims. Integrated spectroscopy of a sample of 17 blue concentrated Large Magellanic Cloud (LMC) clusters is presented and its spectral evolution studied. The spectra span the range (3600-6800) Å with a resolution of
14 Å FWHM, being used to determine cluster ages and, in connection with their spatial distribution, to explore the LMC structure and cluster formation history.
Methods. Cluster reddening values were estimated by interpolation, using the available extinction maps. We used two methods to derive cluster ages: (i) template matching, in which line strengths and continuum distribution of the cluster spectra were compared and matched to those of template clusters with known astrophysical properties, and (ii) equivalent width (EW) method, in which new age/metallicity calibrations were used together with diagnostic diagrams involving the sum of EWs of selected spectral lines (K Ca II, G band (CH), Mg I, H,
H
and H
).
Results. The derived cluster ages range from 40 Myr (NGC 2130 and SL 237) to 300 Myr (NGC 1932 and SL 709), a good agreement between the results of the two methods being obtained. Combining the present sample with additional ones indicates that cluster deprojected distances from the LMC center are related to age in the sense that inner clusters tend to be younger.
Conclusions. Spectral libraries of star clusters are useful datasets for spectral classifications and extraction of parameter information for target star clusters and galaxies. The present cluster sample complements previous ones, in an effort to gather a spectral library with several clusters per age bin.
Key words: galaxies: star clusters - Magellanic Clouds - techniques: spectroscopic
Star clusters in the Magellanic Clouds can help us to have a better insight
about the star formation history of the galaxies as a whole
(Piatti et al. 2002a; Mighell et al. 1998; Piatti et al. 2001). In particular, the large number
of young and intermediate-age star clusters in the Large Magellanic Cloud
(LMC) composes an ensemble which allows to probe the galaxy structure and
evolution during the
last few Gyrs. Indeed, the total estimated number of clusters in the LMC
is 4200 (Hodge 1988), which compares to the catalogued 4089
clusters, associations with characteristics of somewhat looser clusters
and newly-formed ones (Bica et al. 1999).
Episodes of enhancement or reduction in cluster formation
can be traced back in the galaxy's history if statistically meaningful
samples of clusters with well-determined properties are considered.
The analysis of continuum and absorption lines of the integrated spectra of star clusters has provided accurate estimates of their ages and metallicities. The determination of these properties for distant and/or compact clusters, in which individual stars are not observable, needs a calibration based on clusters of well-known ages and metallicities. These age and metallicity calibrations are established by means of equivalent width (EW) measurements, in particular for conspicuous lines of H, Fe and Mg in the visible spectra. Based on these data it is possible to characterize clusters of unknown properties as well as to use the integrated spectra as reference populations (templates) to investigate more complex stellar systems. In order to accomplish these two aims, integrated spectra of star clusters in the Magellanic Clouds (e.g. Piatti et al. 2005; Santos et al. 1995; Bica et al. 1990) and in the Galaxy (e.g. Piatti et al. 2002b; Santos & Bica 1993) have been observed, catalogued (Santos et al. 2002), and analyzed (Santos & Piatti 2004, hereafter SP).
The objective of the present work is to derive age for a sample of 17 LMC clusters, all of type II in the sequence defined by Searle et al. (1980, hereafter SWB), which corresponds to the age range (30-70 Myr) according to posterior photometric calibrations (e.g. Bica et al. 1996). The compact nature and high surface brightness of the present sample clusters make them good targets for integrated spectroscopy. For such young clusters we derived ages but it was not possible to derive metallicities using the proposed methods. Even though, cluster ages and deprojected positions provide useful information on the LMC properties. When our age estimates are combined with a larger sample (Piatti et al. 2003b) with ages in the same scale as in the present work, the interpretation of the data regarding the LMC structure and cluster formation history is improved.
The present cluster sample complements previous ones, in an attempt to provide a spectral library with several clusters per age bin. At the same time, we study the clusters themselves individually, determining their ages and analyzing their spatial distribution, in order to explore the LMC structure and cluster formation history. To estimate the clusters' ages, we employ the new calibrations and diagnostic diagrams recently provided by SP for visible integrated spectra, along with template spectra (e.g. Ahumada et al. 2002; Santos et al. 1995). The first method is based on the EW of Balmer H, Mg and Fe lines and the second one on the continuum distribution of clusters with well-known properties. We confirm the reliability of the procedure proposed by SP to determine clusters' ages, since we included in the sample not only unstudied or poorly studied clusters, but also the cluster NGC 1839 which has been previously observed by means of Washington photometry (Piatti et al. 2003a).
The observations were carried out at the Complejo Astronómico El Leoncito
(CASLEO, San Juan, Argentina) with the
2.15 m telescope during four nights in December 2003. We employed a
CCD camera
containing a Tektroniks chip of
pixels attached to a REOSC
spectrograph, the size of each pixel being
m; one pixel
corresponds to 0.94
on the sky. The slit was set in the East-West
direction and the observations were performed by scanning the slit across
the objects in the North-South direction in order to get a proper sampling
of cluster stars. The long slit corresponding to 4.7
on the sky,
allowed us to sample regions of the background sky.
We used a grating of 300 grooves mm
,
producing an average
dispersion in the observed region of
140 Å/mm (3.46 Å/pixel).
The spectral coverage was
(3600-6800) Å. The slit width was 4.2
,
providing a resolution
(FWHM) of
14 Å, as deduced from the comparison lamp lines.
The reduction of the spectra was carried out with the IRAF package at the Observatorio
Astronómico de Córdoba (Argentina) following the
standard procedures. Summing up, we subtracted the bias and used flat-field
frames previously combined to correct the frames for high and low spatial
frequency variations. We also checked the instrumental signature with the
acquisition of dark frames. Then, we performed the background sky subtraction
using pixel rows from the same frame, after having cleaned the background sky
regions from cosmic rays. We controlled that no significant background sky
residuals were present on the resulting spectra. The cluster spectra,
which were
extracted along the slit according to the cluster size and available flux,
were then wavelength calibrated by fitting observed Cu-Ar-Ne
comparison lamp spectra with template spectra. The rms errors involved in these
calibrations are on average 0.40 Å. Finally, extinction correction
and flux calibrations
derived from the observed standard stars were applied to the cluster spectra.
In addition,
cosmic rays on them were eliminated. Table 1
presents the cluster sample including the designations in different
catalogues, the equatorial coordinates, the average diameters according to
Bica et al. (1999) and the averaged signal-to-noise (S/N) ratios of the
spectra. Typically, three spectra were obtained per cluster.
Table 1: The sample clusters.
A colour-colour diagram of the sample is presented in Fig. 1, where the clusters integrated UBV colours are situated among those from a larger sample (Bica et al. 1996) characterizing the whole age range of LMC clusters.
![]() |
Figure 1: Colour-colour diagrams of the 17 clusters (open circles) and a larger sample (plus signs) representing the LMC cluster system, as observed by means of integrated photometry. The lower panel includes only clusters in the present sample with colours corrected for reddening from Table 3. |
Open with DEXTER |
Before measuring EWs in the observed spectra, these were set to the
rest-frame according to the Doppler shift of H Balmer lines.
Next, the spectra
were normalized at approximately 5870 Å.
The flux normalization at this wavelength is meant to
represent the continuum flux around 5870 Å, avoiding spectral
lines eventually
present. In practice, the spectral region around
5870 Å (20 Å wide) is examined and the normalization
applied to a nearby
wavelength that is representative of the continuum flux.
Spectral fluxes at 3860 Å, 4020 Å, 4150 Å, 4570 Å, 4834 Å,
4914 Å and 6630 Å were used
as guidelines to define the continuum according to Bica & Alloin (1986). The
EWs of H Balmer, K Ca II, G band (CH) and Mg I (
5167 + 5173 + 5184) Å were measured within the spectral windows defined by Bica & Alloin (1986) and using IRAF
task splot. Boundaries for the K Ca II, G band (CH), Mg I, H,
H
and H
spectral windows are, respectively, (3908-3952) Å,
(4284-4318) Å, (5156-5196) Å, (4082-4124) Å, (4318-4364) Å, and (4846-4884) Å. Such a procedure has been consistently applied, making the
EWs from integrated spectra safely comparable with those in the well-known
cluster
database. Table 2 presents these measurements as well
as the sum of EWs of
the three metallic lines (
)
and of the three Balmer lines H
,
H
and H
(
).
Typical errors of
10 % on individual EW
measurements were the result of tracing slightly different continua. If
the sums of EWs
and
are separately used, the EW relative
errors are
improved (
smaller than the individual EW errors).
and
prove to be useful in
the discrimination of old, intermediate-age and young systems
(Dutra et al. 1999; Rabin 1982, SP).
Table 2: Equivalent widths (Å).
As a first approach to get cluster ages, Fig. 2 shows the
cluster sample plotted in the diagnostic
diagrams defined by SP, which
are intended to discriminate cluster ages for systems younger than 10 Gyr,
and metallicities for systems older than 10 Gyr.
In order to show how sensitive the EW sums are
on the size of the spatial profile extraction
and on the stochastic effects produced by few bright stars in the integrated
spectrum, we used the clusters NGC 1902
and SL 709. Two spectra were extracted from the spatial profile of
NGC 1902, one of them sampling a larger cluster extent (resulting a
smaller )
than the other. Besides
the observed integrated spectrum of SL 709, another
spectrum was obtained by subtracting a bright
star superimposed on the cluster bulk profile (resulting a larger
).
Subsequently we only used the NGC 1902
spectrum from its larger extraction, since we
judged that it samples more properly the
cluster population and the SL 709 spectrum with the bright star
subtracted, which
also should better represent the cluster bulk stellar content.
The linked symbols in the diagrams correspond to two different spectral
extractions from the spatial profiles of NGC 1902 and SL 709,
which exemplifies the sensitivity of the integrated spectra on the
size of the extraction (for NGC 1902) and on the presence of a bright star
superimposed on the cluster bulk profile (for SL 709).
The integrated spectra
of NGC 1902 resulted in EW sums which place the cluster
in the same region of the diagrams, whereas for SL 709
a transition occurs between two regions with intermediate ages.
Indeed,
SL 709 seems to be still further affected by an additional
superimposed
bright star of a younger age, which was taken into account in the analysis (see
Sect. 3.2).
![]() |
Figure 2: Diagnostic diagrams built from EW sums of H Balmer and metallic lines for the 17 clusters. Continuous lines divide regions of different age ranges. The corresponding EW sums for both spectra of SL 709 and NGC 1902 are linked. |
Open with DEXTER |
The
clusters were then age-ranked according to the calibrations provided by
SP. Since the EW of each Balmer line is a
bivalued function of age with a maximum around 300 Myr, we
used
to get a first age estimate using:
![]() |
(1) |
We then used
to get a second age estimate guided by the
previous
estimate, since from
two solutions are possible:
![]() |
(2) |
The average of these two estimates is listed in Col. 5 of Table 3.
Table 3: Cluster parameters.
The template matching method consists in achieving the best possible match between the analyzed cluster spectrum and a template spectrum of known age and metallicity. In this process we selected, from the available template spectra, the ones which minimize the flux residuals, calculated as the normalized difference (cluster - template)/cluster. Note that differences between cluster and template spectra are expected to be found due to variations in the stellar composition of the cluster, such as the presence of a relatively bright star with particular spectral features or contamination of a field star close to the direction towards the cluster.
All 17 clusters in our sample are well represented by blue stellar populations, according to their spectral properties. For the present sample, the useful template spectra are: Yd (40 Myr), Ye (45-75 Myr), Yf (100-150 Myr), Yg (200-350 Myr) and Yh (0.5 Gyr), which represent young and intermediate-age populations built from Galactic open clusters (Piatti et al. 2002b).
Since the continuum distribution is also affected by reddening, we firstly adopted a colour excess E(B-V) for each cluster, taking into account the Burstein & Heiles (1982, hereafter BH) extinction maps. Secondly, we corrected the observed spectra accordingly and then we applied the template match method. The results are shown in Figs. 6 to 22.
Three clusters (SL 237, SL 508 and SL 709) have their spectra affected by a bright star, which was taken into account according to the following procedure: we matched the spectrum of SL 237 with a combination of templates Yd (40 Myr) and Ye (60 Myr) and an average spectrum of giant stars of early M type from Silva & Cornell (1992). The stellar spectrum contributes with 35% of the total light at 5870 Å (see Fig. 20). As for SL 237, the spectrum of SL 508 shows bands redwards of 5000 Å, characteristic of late spectral types. The same stellar spectrum used in the spectral matching of SL 237 was also used for SL 508 in combination with template Ye. The star contributes with 15% of the integrated light at 5870 Å (see Fig. 21). The same procedure was applied to SL 709, which presents a flat continuum beyond 5000 Å. In the case of SL 709, the spectrum of the Carbon star TT Tau (Barnbaum et al. 1996), with 20% of contribution to the total flux at 5870 Å, was combined with template Yf (see Fig. 22).
The ages determined by the two methods, together with estimates from the literature (whenever available), were used to get final averaged ages (Table 3). Their respective errors take into consideration the dispersion of the values averaged. The reddening adopted in the template match method is presented in Col. 2 of Table 3. The LMC is optically thin, the average foreground and internal E(B-V)colour excesses being 0.06 and 0.06 mag, respectively (Dutra et al. 2001). This explains why the derived E(B-V) values turn out to be relatively small, being all lower than 0.1 mag.
Ten of the clusters in our sample are also present in OGLEII sample (Pietrzynski et al. 1999). Pietrzynski & Udalski (2000) determined their ages by fitting isochrones of Z=0.008 on cluster colour-magnitude diagrams (CMDs) built with OGLEII data (see Table 3). By analyzing the distribution in the BV photographic CMD of cluster member stars, Alcaino & Liller (1987) derived age for NGC 1839, NGC 1870 and SL 237 considering the position of the main sequence turnoff, the position of the brightest blue star and the fitting of Maeder & Mermilliod (1981) isochrones. Their results are displayed in Table 3. Bono et al. (2005) determined the age of NGC 1943 by analyzing the pulsational period of its 9 Cepheids and found it consistent with the value obtained from isochrone fitting on OGLEII data. The value shown in Table 3 is an average of their results.
Within the expected uncertainties, the ages derived in the present work
agree with those given in the literature.
![]() |
Figure 3: Spatial distribution of LMC clusters with ages in a homogeneous scale. Upper panel: the present sample and the sample analyzed in Piatti et al. (2003b) are plotted in equatorial coordinates relative to the LMC center. The positions of the bar and 30 Dor are indicated. Lower panel: the same data is plotted with different symbols discriminating different age ranges. |
Open with DEXTER |
![]() |
Figure 4: Deprojected distance from the LMC center (in degrees) as a function of age for the same sample of Fig. 3. |
Open with DEXTER |
![]() |
Figure 5: Distribution of cluster age with integrated colours for the 17 clusters plus the sample used by Piatti et al. (2003b). |
Open with DEXTER |
Integrated UBV colours of 624 LMC associations and clusters were
obtained by Bica et al. (1996) from photoelectric photometry observations.
We used these data to check the trend of (U-B) and (B-V) with age
according to the present estimates. For clusters with integrated
photometry, Fig. 5 shows
this comparison for the present sample and for the one
presented in Piatti et al. (2003b), totalling 53 clusters with uniform
age estimates. The colour gap seen in both (U-B) and (B-V) is a real
feature first identified by van den Bergh (1981) and supported by
larger cluster samples (Bica et al. 1996). The gap is probably
a natural consequence of cluster evolution with increasing
metallicities towards the present, and epochs of
reduced cluster formation between 300 Myr and
1 Gyr
(e.g. Girardi et al. 1995). At least in the LMC bar, such a period of
reduced cluster
formation is not observed for the field stars (Smecker-Hane et al. 2002).
![]() |
Figure 6:
From top to bottom:
the observed integrated spectrum of NGC 1804, the spectrum corrected
for the adopted reddening E(B-V), the template spectrum which best
matches it, and the flux residuals according to
![]() |
Open with DEXTER |
We have determined ages of blue LMC star clusters using their integrated spectra. The results obtained from two independent methods are in agreement with the literature estimates for clusters studied in common.
The present sample, added to a set of
young and intermediate-age clusters with uniformly determined ages,
was used with the aim of analyzing their spatial distribution. Clusters
closer to the LMC bar tend to be younger than those at the LMC disk.
The colour gap at
and
noticed
in studies with larger cluster samples appears to be
related to the natural cluster evolution, as suggested by the present
analysis. However, the older groups SWB III, SWB IVa and SWB IVb
must be observed spectroscopically and in detail by means of CMDs,
testing also the possibility of decreased cluster formation
contributing to the gap.
Acknowledgements
We would like to thank the CASLEO staff members and technicians for their kind hospitality and support during the observing runs. We are grateful for the use of the CCD and data acquisition system at CASLEO, supported under US National Science Foundation (NSF) grant AST-90-15827. We also thank the anonymous referee for helpful remarks. This work was partially supported by the Brazilian institutions CNPq and FAPEMIG, and the Argentinian institutions CONICET, Agencia Nacional de Promoción Científica y Tecnológica (ANPCyT) and Agencia Córdoba Ciencia.
![]() |
Figure 7: Same as Fig. 6 for NGC 1839. |
Open with DEXTER |
![]() |
Figure 8: Same as Fig. 6 for NGC 1870. |
Open with DEXTER |
![]() |
Figure 9: Same as Fig. 6 for NGC 1894. |
Open with DEXTER |
![]() |
Figure 10: Same as Fig. 6 for NGC 1902 (larger extraction). |
Open with DEXTER |
![]() |
Figure 11: Same as Fig. 6 for NGC 1913. |
Open with DEXTER |
![]() |
Figure 12: Same as Fig. 6 for NGC 1932. |
Open with DEXTER |
![]() |
Figure 13: Same as Fig. 6 for NGC 1940. |
Open with DEXTER |
![]() |
Figure 14: Same as Fig. 6 for NGC 1943. |
Open with DEXTER |
![]() |
Figure 15: Same as Fig. 6 for NGC 1971. |
Open with DEXTER |
![]() |
Figure 16: Same as Fig. 6 for NGC 2038. |
Open with DEXTER |
![]() |
Figure 17: Same as Fig. 6 for NGC 2118. |
Open with DEXTER |
![]() |
Figure 18: Same as Fig. 6 for NGC 2130. |
Open with DEXTER |
![]() |
Figure 19: Same as Fig. 6 for NGC 2135. |
Open with DEXTER |
![]() |
Figure 20: Same as Fig. 6 for SL 237. The matched spectrum is a combination of templates Ye (40 Myr) and Yd (60 Myr) and an average spectrum of early M giants (Silva & Cornell 1992). The flux fraction at 5870 Å is indicated. |
Open with DEXTER |
![]() |
Figure 21: Same as Fig. 6 for SL 508. The matched spectrum is a combination of template Ye (40 Myr) and an average spectrum of early M giants (Silva & Cornell 1992). The flux fraction at 5870Å is indicated. |
Open with DEXTER |
![]() |
Figure 22: Same as Fig. 6 for SL 709. The matched spectrum is a combination of template Yf (100-150 Myr) and the spectrum of a Carbon star (TT Tau from Barnbaum et al. (1996)). The flux fraction at 5870 Å is indicated. |
Open with DEXTER |