Free Access
Issue
A&A
Volume 508, Number 3, December IV 2009
Page(s) 1279 - 1283
Section Galactic structure, stellar clusters, and populations
DOI https://doi.org/10.1051/0004-6361/200912697
Published online 17 September 2009

A&A 508, 1279-1283 (2009)

Model-independent diagnostics of highly reddened Milky Way star clusters: age calibration[*] (Research Note)

Y. Beletsky1 - G. Carraro1,2 - V. D. Ivanov1

1 - ESO, Alonso de Cordova 3107, Santiago de Chile, Chile
2 - Dipartimento di Astronomia, Università di Padova, Vicolo Osservatorio 2, 35122 Padova, Italy

Received 15 June 2009 / Accepted 11 August 2009

Abstract
Context. The next generation near- and mid-infrared Galactic surveys will yield a large number of new highly obscured star clusters. Detailed characterization of these new objects with spectroscopy is time-consuming.
Aims. Diagnostic tools that will be able to characterize clusters based only on the available photometry will be needed to study large samples of the newly found objects.
Methods. The brightness difference between the red clump and the main-sequence turn-off point have been used as a model-independent age calibrator for clusters with ages from a few 108 to 1010 yr in the optical. Here we apply for the first time the method in the near-infrared.
Results. We calibrated this difference in K-band, which is likely to be available for obscured clusters, and we apply it to a number of test clusters with photometry comparable to the one that will be yielded by the current or near-future surveys.
Conclusions. The new calibration yields reliable ages over the range of ages for which the red clump is present in clusters. The slope of the relation is smoother than that of the corresponding V-band relation, reducing the uncertainty in the age determinations with respect to the optical ones.

Key words: open clusters and associations: general - surveys

1 Introduction

The recent advance in the infrared (IR) instrumentation made it possible to carry out large scale sky surveys. The Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006) was the first one to provide a deep look into the most obscured regions of the Milky Way. A number of new clusters were found in 2MASS (i.e. Dutra & Bica 2000; Ivanov et al. 2002; Bica et al. 2003; Borissova et al. 2003). Currently, the UKIRT Infrared Deep Sky Survey (UKIDSS; Lawrence et al. 2007) is creating a map of selected Northern sky regions, typically $\sim$5 mag deeper than the 2MASS, and in the near future the Visible and Infrared Survey Telescope for Astronomy (VISTA; Emerson et al. 2006; Dalton et al. 2006) will provide similar coverage in the South. These projects include specialized multicolor Galactic plane surveys: UKIDSS Galactic plane surveys (GPS; Lucas et al. 2008) and VISTA variables in the Via Lactea (VVV; Minniti et al. 2006) which are expected to reveal a large number of new star clusters.

Their characterization is a time-consuming process, requiring spectroscopy in the near-IR, because of the strong intervening dust absorption ( $A_K\sim0.1\times A_V$; Rieke & Lebofsky 1985). This prompted us to attempt to develop various methods to determine the cluster parameters from the photometry alone. For example, the metal abundances of globular clusters is easily measured with an error of $\sigma({\rm [Fe/H]})\sim0.2$ dex (i.e. Ferraro et al. 1999; Ivanov & Borissova 2002) from the slope of the RGB on the color-magnitude diagrams (CMD).

In this Note we focus on a new age calibrator based on the K-band. difference between the clump and the turn off point (TO) in the CMD of star clusters older than about 300 Myr, which will be widely used to get a quick and approximate age values for newly discovered IR star clusters.

The layout of the paper is as follows. In Sect. 2 we introduce the concept of age calibrator as derived and used in the optical regime. Section 3 describes the sample of very well know star clusters on which we base our new age calibrator. A few analytical relations are given in Sect. 4. These relations are then applied to a few test-cases in Sect. 5. Section 6, finally, summarizes out findings.

2 Age calibrations for open star clusters

The magnitude difference $\Delta V$ between the main sequence (MS) TO and the clump of He burning stars is a well known age indicator for Galactic open clusters. The very requirement to have a well-defined red clump limits it to clusters older than $\sim$300 Myr when the He burning stars appear. However, the method is distance and reddening independent, as demonstrated early on by Cannon (1970), and it is not surprising that it has been widely used ever since (see, for instance, Carraro & Chiosi 1994; Salaris et al. 2004).

Different, more sophisticated versions of this calibrator have been devised in the past, in the optical regime, by Anthony-Twarog & Twarog (1985) and Phelps et al. (1994). They also include the color difference between the TO and the red giant branch (RGB). However, their application is more limited because usually the RGBs of open clusters are poorly populated, and the RGB colors are ill-defined.

The simple difference $\Delta V$ between the red clump and TO - the two major features in the CMD - has the advantage of providing a quick age estimate even from un-calibrated CMDs (Phelps et al. 1994; Salaris et al. 2004). As described by Carraro & Chiosi (1994) this indicator shows a slight metallicity dependence in the optical: $\Delta V$ at fixed age is larger for higher metallicity clusters than for lower metallicity ones.

The aim of the present work is to extend this calibrator to the near-IR, where the metallicity effect should be weaker. We will use the 2MASS $K_{\rm S}$ filter as the least affected by the extinction, and throughout the Note will mark it as K, for simplicity.

3 Sample selection and indicator definition

We selected 14 open clusters with 2MASS photometry, with reliable age and metallicity estimates, spanning a wide range of the parametric space (Table 1). The data come from the latest edition of the WEBDA[*] database maintained by E. Paunzen at Vienna University. Throughout the paper we adopted a conservative estimates of the uncertainty in the age of 15% and a factor of 2 for the abundance [Fe/H].

Table 1:   Parameters of the clusters' sample.

Some of the clusters suffer from significant fore- and background contamination since they are located low onto the Galactic plane. To minimize it, we built up homogeneous CMDs considering only the stars within the accepted cluster radii estimates from Dias et al. (2002). The CMDs for M 67 and NGC 2506 - two clusters with rather different ages and abundances - are shown in Fig. 1.

\begin{figure}
\par\includegraphics[width=7.8cm,clip]{AA12967fig1.eps}
\end{figure} Figure 1:

Definition of the parameter $\Delta K$. The left panel shows the CMD of M 67 (solar metallicity, 4 Gyr old), the right shows the CMD of NGC 2506 (significantly subsolar metallicity, twice younger than M 67).

Open with DEXTER

Next, we identified the clump and TO magnitudes, respectively at $K^{\rm clump}=8.00\pm0.05$, and $K^{\rm TO}=11.80\pm0.10$ mag, for M 67. This implies $\Delta K=3.8\pm0.1$ mag. For NGC 2506 we found $K^{\rm clump}=10.85\pm0.15$, and $K^{\rm TO}=14.00~\pm~0.20$ mag, yielding $\Delta K=3.15\pm0.25$ mag. These examples illustrate that for an old cluster like M 67 the definition of the clump magnitude is easier despite the smaller number of He burning stars, while the clump at younger ages is wider both in color and magnitude due to evolutionary effects (Girardi et al. 2000).

The measured values of $\Delta K$ are also listed in Table 1. Their uncertainties were obtained tentatively, by eye. These take into account the broadening of features in the CMD due to photometric errors, patchy extinction and unresolved binaries. The experience shows that such estimates can be even more reliable than the automated ones (i.e. see the discussion in Sect. 4 of Phelps & Janes 1994). The age (logarithm) versus $\Delta K$ relation is plotted in Fig. 2.

We emphasize that this IR relation looks smoother than the optical counterpart (Carraro & Chiosi 1994) mainly because of the significantly weaker metallicity dependence (see below). Furthermore, the smooth points distribution in Fig. 2 makes any analytical fit through the data easier and more robust.

\begin{figure}
\par\includegraphics[width=7.8cm,clip]{AA12967fig2.eps}
\end{figure} Figure 2:

Age vs. $\Delta K$ relation for the clusters listed in Table 1. In the same Table we report the uncertainties in age and $\Delta K$.

Open with DEXTER

4 Calibration

An analytical representation of the age versus $\Delta K$relation is needed to facilitate the derivation of cluster ages. The relation appears non-linear, prompting us to fit a higher order polynomial:

\begin{displaymath}\log({\rm Age})\!=\!9.09(0.29)+0.24(0.20)\times\Delta K+0.10(0.03)\times(\Delta K)^2~~
\end{displaymath} (1)

with rms = 0.28. The reduced $\chi^2$ of the fit is 0.08.

For comparison, the linear fit to the data is:

\begin{displaymath}\log({\rm Age})= 8.27(0.1)+0.35(0.03)\times\Delta K
\end{displaymath} (2)

with rms = 0.36. Numbers in brackets are corresponding uncertainties.The reduced $\chi^2$ of this fit is 0.04.

The age in these relations is expressed in Gyr. The two fits are nearly indistinguishable for $2.5\leq \Delta K \leq3.5$ mag - the deviations from the linearity occur in fact at the extremes of the age range. The fits are plotted in the left panel of Fig. 3, while in the right panel we show the optical relation from Carraro & Chiosi (1994), for the sake of comparison.

As anticipated in the previous Section, the data point distribution in the IR plot (left panel) is smoother than the optical one (right panel), and allow a better fit through the points.

\begin{figure}
\par\includegraphics[width=9cm,clip]{AA12967fig3.eps}
\end{figure} Figure 3:

Left panel: age vs. $\Delta K$ relation for the 14 clusters listed in Table 1. The polynomial fit (Eq. (2)) is showed as a dashed line, and the linear one (Eq. (1)) - as a solid line. Right panel: the optical relation from Carraro & Chiosi (1994). The solid line is a linear fit through the data.

Open with DEXTER

Next, we include the metallicity in the fit, obtaining a 3-dimensional relation:

                $\displaystyle \log({\rm Age})$ = $\displaystyle 8.28(0.08)+0.34(0.03)\times\Delta K$  
    $\displaystyle +0.00005(0.1)\times{\rm [Fe/H]}$ (3)

with rms = 0.06 (Fig. 4). The reduced $\chi^2$ of the fit is 0.004.

\begin{figure}
\par\includegraphics[width=9cm,clip]{AA12967fig4.eps}
\end{figure} Figure 4:

Age vs. $\Delta K$ for the 14 clusters listed in Table 1, including the metallicity dependence. The projection on the age vs. $\Delta K$ is also shown.

Open with DEXTER

As expected, the dependence on metallicity is negligible, making $\Delta K$ very useful, since in most cases the newly discovered clusters will not have metallicity estimates available. In comparison with the optical relation (Carraro & Chiosi 1994, Eq. (2)), the linear relation in the IR is less sensitive to metallicity, and also shows a somewhat shallower dependence on $\Delta K$ (0.34 versus 0.45, see Carraro & Chiosi 1994).

As for the uncertainty on the age, it is easy to show - by simple propagation - that in the linear cases (Eqs. (2) and (3)), the uncertainty in the age does not depend on the value of $\Delta K$. This is not the case of course for Eq. (2), which has a quadratic dependence on $\Delta K$.

\begin{figure}
\par\includegraphics[width=14.5cm,clip]{AA12967fig5.eps}
\end{figure} Figure 5:

Application of the new age indicator to NGC 6791, NGC 2141 and Berkeley 17 ( from left to right). $\Delta K$ is marked on each CMD.

Open with DEXTER

5 Validation of the new calibrations

To demonstrate the power of the new calibrations, we applied them to estimate the ages of three clusters with shallow 2MASS photometry, but that have IR photometry from other sources: NGC 6791, NGC 2141, and Berkeley 17. Spectroscopic metallicity estimates are also available for them, allowing us to use all three relations in a metallicity range which encompasses the typical metal content of Galactic open clusters. The results are listed in Table 2. The age uncertainties were derived propagating the measurement errors through the calibrations.

Table 2:   New age estimates for the three test-cases discussed in this paper.

5.1 NGC 6791

IR photometry for this cluster is taken from Carney et al. (2005). It was acquired with the IRIM camera at the 4-m Mayall telescope at Kitt Peak National Observatory. The data were photometrically calibrated with the UKIRT standards (Casali & Hawarden 1992), which were also used to calibrate the 2MASS itself. Not surprisingly, Eq. (37) from Carpenter (2001):

                      $\displaystyle K_{\rm 2MASS}$ = $\displaystyle K_{\rm UKIRT}+0.004(0.006)\times(J-K)_{\rm UKIRT}$  
    +0.002(0.004) (4)

indicates that the color difference $\Delta (J-K) \sim0.3$ mag between the TO and clump leads to only $0.0012\pm0.0018$ mag difference between the K-band magnitudes of the two systems - much smaller than our measurement errors. The constant term is negligible for calculating the $\Delta K$ parameters, and it is statistically indistinguishable from zero, anyway.

The CMD of NGC 6791 is shown in the left panel of Fig. 5. Widely accepted values for age and metallicity are in the range of 7-9 Gyr and $\rm [Fe/H]=+0.35$ to +0.45 dex, respectively (Carney et al. 2005, and references therein). Our Eqs. (1)-(3), yield ages of 9.3, 6.5 and 6.0 Gyr, respectively.

5.2 NGC 2141

IR photometry for this cluster is taken from Carraro et al. (2001). It was acquired with the ARNICA camera at the Tirgo telescope. The data were photometrically calibrated with the ARNICA standards (Hunt et al. 1998). Equation (8) from Carpenter (2001):

                        $\displaystyle K_{\rm 2MASS}$ = $\displaystyle K_{\rm ARNICA}-0.024(0.011)\times(J-K)_{\rm ARNICA}$  
    +0.012(0.006) (5)

indicates that the color difference $\Delta (J-K) \sim0.25$ mag between the TO and clump leads to only ${\sim} 0.006\pm0.003$ mag difference between the K-band magnitudes of the two systems - again, much smaller than our measurement errors.

The CMD of NGC 2141 is shown in the middle panel of Fig. 5. For this cluster, metallicity is $\rm [Fe/H]=-0.26$, and the age is $\sim$2.5 Gyr (Carraro et al. 2001). We obtained ages of 2.2, 2.4 and 2.3 Gyr, respectively from Eqs. (1)-(3). These results are in good agreement with the isochrone fitting estimates (Carraro et al. 2001).

5.3 Berkeley 17

IR photometry for this cluster is taken from Carraro et al. (1999). It was acquired with the ARNICA camera at the Tirgo telescope, and photometrically calibrated with the ARNICA standards (Hunt et al. 1998). As shown in the previous subsection, our age measuring parameter $\Delta K$is not affected significantly by using a different photometric system than the 2MASS.

The CMD of Berkeley 17 is shown in the right panel of Fig. 5. This cluster is known to have marginally subsolar abundance $\rm [Fe/H]=-0.10$ dex and is it 9-10 Gyr old (Friel et al. 2005). Our calibrations yielded ages of 10.0, 7.6 and 7.0 Gyr, respectively for Eqs. (1)-(3).

This exercise shows that the age of a cluster within the linear regime of our calibrator ( $2.5\leq \Delta K \leq3.5$ mag) is closer to the isochrone-derived one when simple linear (Eq. (1)), or the metal dependent (Eq. (3)) fitting are used. The most deviating relation is the polynomial one. In the other two cases of very old clusters, the polynomial relation (Eq. (2)) is the most accurate, and provides age estimates very close to the widely accepted in the literature. Finally, the new calibrations provide, in the linear regime, age estimates that agree quite well with isochrone-based ages, while outside this range the difference might be of the order of 30%.

6 Summary and conclusions

In this paper we described relations of the magnitude difference between the TO and the clump $\Delta K$ in the IR CMD versus age for 800-1000 Myr old open star clusters. The age limits are set by the life time of the red clump stars and available photometric data. Our goal is to develop a simple tool to measure the age of newly detected star clusters from the next generation IR sky surveys. The relations are derived in the widely used 2MASS photometric system. We provide linear and second order polynomial fits of the age versus $\Delta K$ relationship, and a multi-linear fit which also takes into account the metallicity of the clusters, when available. However, it is unlikely that the metallicity will be immediately available for most of the new objects, so Eqs. (1) and (2) will be mostly used used for en-masse analysis of the new clusters.

Beside, we show that the dependence on metallicity is very shallow, which makes these relations of more robust use than their optical counterpart.

We then illustrate with three examples that the simple linear relation works better for intermediate $\Delta K$ values, while the polynomial relation is preferable at the extremes of the validity range.

The method can potentially be useful for a wider age range, and these calibrations should be extended in the future to cover both younger and older clusters.

Acknowledgements
This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, and of the WEBDA database maintained by E. Paunzen at Vienna University.

References

Footnotes

... calibration[*]
Original data are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/508/1279
... WEBDA[*]
http://www.univie.ac.at/webda/navigation.html
Copyright ESO 2009

All Tables

Table 1:   Parameters of the clusters' sample.

Table 2:   New age estimates for the three test-cases discussed in this paper.

All Figures

  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{AA12967fig1.eps}
\end{figure} Figure 1:

Definition of the parameter $\Delta K$. The left panel shows the CMD of M 67 (solar metallicity, 4 Gyr old), the right shows the CMD of NGC 2506 (significantly subsolar metallicity, twice younger than M 67).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{AA12967fig2.eps}
\end{figure} Figure 2:

Age vs. $\Delta K$ relation for the clusters listed in Table 1. In the same Table we report the uncertainties in age and $\Delta K$.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{AA12967fig3.eps}
\end{figure} Figure 3:

Left panel: age vs. $\Delta K$ relation for the 14 clusters listed in Table 1. The polynomial fit (Eq. (2)) is showed as a dashed line, and the linear one (Eq. (1)) - as a solid line. Right panel: the optical relation from Carraro & Chiosi (1994). The solid line is a linear fit through the data.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{AA12967fig4.eps}
\end{figure} Figure 4:

Age vs. $\Delta K$ for the 14 clusters listed in Table 1, including the metallicity dependence. The projection on the age vs. $\Delta K$ is also shown.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=14.5cm,clip]{AA12967fig5.eps}
\end{figure} Figure 5:

Application of the new age indicator to NGC 6791, NGC 2141 and Berkeley 17 ( from left to right). $\Delta K$ is marked on each CMD.

Open with DEXTER
In the text


Copyright ESO 2009

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.