A&A 488, 935-941 (2008)
DOI: 10.1051/0004-6361:200810201
M. A. T. Groenewegen
Koninklijke Sterrenwacht van België, Ringlaan 3, 1180 Brussels, Belgium
Received 15 May 2008 / Accepted 26 June 2008
Abstract
Context. Over the past decade the use of the red clump (RC) as distance indicator has increased in importance as this evolutionary phase is well populated and a good local calibration exists.
Aims. The absolute calibration of the RC in the I and K band is investigated again based on the recently published revised Hipparcos parallaxes.
Methods. A numerical model is developed that takes the various selection criteria and the properties of the Hipparcos catalogue into account. The biases involved in applying certain selections are estimated with this model.
Results. The absolute magnitudes that are found are
and
(on the 2MASS system). The I-band value is in good agreement with previous determinations, the K-band value is fainter than previously quoted, and this seems to be related to a selection bias whereby accurate K-magnitudes are only available for relatively few bright stars. Applying population corrections to the absolute K magnitude of RC stars in clusters supports the fainter magnitude scale.
Key words: stars: distances - cosmology: distance scale
Paczynski & Stanek (1998) were the first to advocate the use of red clump (RC) stars in the I-band as distance indicators. Two advantages are that this phase of stellar evolution is usually well populated, and that many RC stars exist with good parallaxes in the Hipparcos catalogue.
They consider stars with an error in parallax smaller than 10%, located
in the box 0.8 < (V-I) < 1.25 and
-1.4 < MI < 1.2 and with
Hipparcos flag H42 equal to one of ``A, C, E, F, G'', which means that the
star has one or more I-band measurements in the Cousins, Johnson, or
Kron-Eggen system. They find no significant trend of MIon colour and derived
for the entire sample.
They then considered the 228 stars with d<70 pc and average distance
of 50 pc to find
,
and
for the rest of the sample at an average distance
of 106 pc. The difference is ascribed to reddening, and when extrapolating
to 0 pc, they arrive at their final result of find
M I = -0.279and assign an error bar of 0.088. In deriving these number, they also
applied a correction for ``a distance bias''.
Stanek & Garnavich (1998) considered the same sample, and their final
value of
is based on the sub sample of
stars within 70 pc and they argue that no correction for reddening is
needed, and they do not apply a correction for selecting the sample on parallax.
Udalski (2000) considered the effect of metallicity. He took the same
cut in parallax error and Hipparcos H42 flags, but also included stars
with flag ``H'' to increase the sample and considered stars with a
metallicity determination from high-resolution spectroscopy by
McWillam (1990). Based on this sample of 284 stars, he finds a weak
dependence on metallicity:
[Fe/H]
.
Alves (2000) was the first to consider the K-band for RC stars as
a standard candle. He considered the 284 stars from Udalski (2000) and
finds that 238 of them have a K-magnitude in the Two Micron Sky
Survey (TMSS, Neugebauer & Leighton 1968; stars are designated by
``IRC'' numbers). Assuming no reddening and considering stars with
-2.5 < MK < -0.8 and all (V-K) colours, he finds
and a weak dependence on metallicity
[Fe/H]
.
The RC has gained importance as a distance indicator since the Paczynski & Stanek (1998) paper as shown for example by distance estimates based on the RC to Fornax (Rizzi et al. 2007), M 33 (Kim et al. 2002), or IC 1613 (Dolphin et al. 2001).
With the recent publication of the revised Hipparcos parallax data, it is timely to investigate the absolute magnitude of the RC again. In Sect. 2 the revised Hipparcos data is introduced with auxiliary data used in the analysis. Section 3 describes the numerical model used for estimating selection biases. Section 4 describes the results, and Sect. 5 presents a brief summary and discussion.
The Hipparcos data used are the parallax and error in the parallax from the re-reduction of the raw data by van Leeuwen (2007). For other data (V, I, Hipparcos flags), the originally published catalogue is used (ESA 1997).
Using the script submission tool within VizieR (Ochsenbein et al. 2000) the list of Hipparcos objects was correlated with selected measurements available, in particular, spectroscopic [Fe/H] values (available for 2445 stars, and adding 11 stars from Zhao et al. 2001; 64 from Liu et al. 2001; and 178 from Mishenina et al. 2006), the TMSS catalogue (Neugebauer & Leighton 1969, available for 3167 stars), and Strömgren photometry (available for 36 337 stars). In case of multiple [Fe/H] determinations the mean was taken.
Independently, the Hipparcos objects are correlated on position with
the DENIS (Epchtein et al. 1999, for 42 696 stars) and 2MASS (Skrutskie
et al. 2006, for 117 358 stars) NIR catalogues.
The TMSS was the main source of K-magnitudes used by Alves (2000).
This is the only source of uniform K-magnitude for bright objects
that often saturate in the DENIS and 2MASS surveys. In the present
study TMSS, DENIS and 2MASS are being used but this requires
transformation to the 2MASS system, which will be the system of choice.
Using 749 objects with a 2MASS photometric quality flag, phqual-flag,
of ``A'' in the K-band and that not saturate in DENIS I (implying
I >10) it is found that DENIS K - 2MASS
with no colour term on
.
There are no stars in common between the TMSS survey and non-saturated
DENIS or 2MASS objects with phqual-flags of ``A'' or ``B'' in K. The
transformation to the 2MASS system is done indirectly via the SAAO
Carter (1990) system. The advantage over a comparison to the Koornneef
(1983) system, which was used by Alves (2000), is that the colour term
and offset with respect to the 2MASS system are smaller (Carpenter 2001, and
updated on the IPAC website
).
Using 62 TMSS stars that are listed in Carter (1990) and that have both
an IRC K and I magnitude available, the following relations are derived:
with a dispersion of 0.04 mag, and eliminating 5 outliers,
with no significant colour term, and with a dispersion of 0.04 mag. The transformation from the Carter to the 2MASS system follows the relation on the quoted IPAC website:
.
The adopted K-band magnitude is, in order of preference,
(1) the 2MASS value, but only when phqual is ``A'';
(2) The DENIS value (for non-saturated stars with
)
corrected by +0.018 mag;
(3) the TMASS/IRC value, transformed to the 2MASS system as described above.
If no
is available a mean colour for RC stars of
is used in the transformation.
Of all Hipparcos stars that have a K-magnitude assigned, 3116 are based on TMSS, 17 on Denis and 103827 on 2MASS.
In this section a model is described to construct synthetic samples of stars and apply various selection criteria in order to generate samples of stars that closely resemble in nature the observed data. The model is follows largely Groenewegen & Oudmaijer (2000) but is based on the properties of the revised Hipparcos catalogue.
The coordinate system used is cylindrical coordinates centred on the
Galactic centre. The galactic distribution is assumed to be a double
exponential disk with a scale height H in the z-direction (the
coordinate perpendicular to the galactic plane), and a scale length
in the galacto-centric direction. The distance to the
Galactic centre is taken to be 7800 pc (Zucker et al. 2006).
The luminosity function can be arbitrary but, in the case of simulating
RC stars, is assumed to be a Gaussian. In I with a mean
MI = -0.27
and
of 0.20. The (V-I) colours is assumed to be a Gaussian
with mean 0.98 and
.
In K the Gaussian has a mean
MK = -1.60 and
of 0.22,
while the (V-K) colours is taken as a Gaussian with mean 2.32 and
of 0.21.
Visual extinction in the simulation and to de-redden the observations
is based on several 3-dimensional reddening models available in the literature.
Marshall et al. (2006) presents reddening in the K-band with
15
sampling along 64 000 lines-of-sight in the direction
and
based on 2MASS data.
Drimmel et al. (2003) presents a reddening model based on the dust
distribution model of Drimmel & Spergel (2001), which is based on COBE/DIRBE data.
Arenou et al. (1992) presents a reddening model based on a comparison of
observed B,V photometry with predicted photometry from a star's spectral type.
Finally, a simple law following Parenago (1940) was considered:
![]() |
(1) |
However, this model is only available over a limited range in galactic coordinates, therefore a Monte Carlo simulation was made to randomly generate galactic coordinates and distances assuming a constant number density of stars within 3 kpc distance of the Sun. If the galactic coordinates were within the range of applicability of the Marshall et al. model, the visual extinction from the other models was determined. In the end, using several 1000 positions, the average and dispersion in the various ratios between extinction models were calculated, and the results are in Table 1, in the sense model listed in the row divided by model listed in the column.
Although the dispersion in all the ratios is quite large, they mostly agree in the mean. It only seems that the Arenou et al. model gives higher extinction than the Drimmel et al. and Parenago model. The finally adopted visual extinction was the average of the models by Parenago, Drimmel et al., and Arenou et al. multiplied by 0.84.
Table 1: Comparison of different reddening models.
The adopted completeness function in Hp magnitude is
![]() |
(2) |
The completeness functions quoted below are derived by inspecting the
ratio of stars that fulfil a particular selection to all stars, as a
function of Hp magnitude. The completeness function of stars with I-band flags ``A+C+E+F+G'' is approximated as
![]() |
(3) |
![]() |
(4) |
![]() |
(5) |
As outlined above, spectroscopic [Fe/H] values have been retrieved
using Vizier and recent literature. The completeness function is
approximated as
![]() |
(6) |
![]() |
(7) |
![]() |
(8) |
![]() |
(9) |
The ``measured'' parallax is based on the true parallax (from the true distance) and the Gaussian distributed parallax error.
Table 2: Results based on different sample selections.
To begin with, a number of input parameters need to be set:
The consecutive steps in the model are the following:
From the properties of the observed sample (e.g. from the observed parallax, or a photometrically determined parallax), it is already possible to eliminate stars beyond a certain limiting distance.
Table 2 lists the observed absolute magnitudes depending
on various selections. These selections largely follow the selections
made by previous works as quoted in the Introduction.
In all cases, only stars where Hipparcos flag H30 (the goodness-of-fit) is <5,
,
and single stars (isoln = 5 as listed in the revised
Hipparcos catalogue) are considered.
Following earlier work, stars are selected inside the box (all colours
and magnitudes are understood to be de-reddened according to the model
described in Sect. 3.3) 0.75 < (V-I) < 1.25,
-1.6 < MI < 1.2.
The data is binned in absolute magnitude and then fitted with the
formula introduced by Paczynski & Stanek (1998):
![]() |
Figure 1: Colour-magnitude diagrams for model 3. The outer box is used to select the data to fit Eq. (10) to. The inner box is used to fit the dependence of the absolute magnitude against colour or metallicity. |
Open with DEXTER |
Several of the entries in Table 2 are put between parentheses, to indicate a lower quality. In the I-band this refers to selections that include lower-quality (V-I) data. In Kthis refers to model 5 that includes too few stars for an accurate fit of the Gaussian.
As the selection box and Eq. (10) are designed to include a
``background'' of stars, the fit of the absolute magnitude against
colour or metallicity is made for stars that are located inside a
smaller box. Based on the fit with Eq. (10) and the coefficients
,
the number of red clump stars can be estimated.
The location of the smaller box is chosen such that it a pproximately
includes this number of stars. In the I-band this is
0.8 < (V-I) < 1.2
and
-0.6 < MI < 0.1. In the K-band this is
1.9 < (V-K) < 2.7
and
-1.9 < MK < -1.35.
First the dependence on colour and metallicity is investigated,
fitting a linear relation. When
is fitted against (V-I)colour, the slope is
,
,
and
,
respectively, when all, ``ACEFGH'' and ``ACEFG'' Hipparcos
I-band flags are considered.
For these three selections, the slope is
,
,
and
when MI is fitted against [Fe/H].
The three determinations are consistent, as one would expect since any dependence
on metallicity should not depend on the origin of the I-band data.
No significant dependence on colour is found, in accordance with
Paczynski & Stanek (1998). Also the dependence on metallicity is marginal.
![]() |
Figure 2: The fit of Eq. (10) to the data in the I-band ( upper panel), and K-band for model 3. |
Open with DEXTER |
Table 3: Results of the simulations.
When MK is fitted against (V-K) colour, the slope is
,
,
and
,
respectively,
when all, ``ACEFGH'' and ``ACEFG'' Hipparcos I-band flags are considered.
Fitted against [Fe/H], the slopes become, respectively,
,
,
and
.
In K there is therefore no evidence of any metallicity dependence
and only a marginal one on (V-K) colour.
For the selections in Table 2, the results of the
simulations are listed in Table 3. As mentioned in
Sect. 3.2, the input values are
MI = -0.27 with
,
and
MK = -1.60 with
= 0.22.
One can observe that the bias can be up to 0.07 mag in I depending
on the selection bias. The fifth column in Table 3 lists
the ``true'' value, which is the observed value from Table 2,
corrected by the bias, estimated from the difference between the input
value and the output value from the simulation.
The average of the 4 available models or the one with smallest correction
(model 5) gives the identical result:
.
In K the biases are larger, up to 0.1 mag, and the average of 7
estimates is
.
![]() |
Figure 3:
Fit of ![]() |
Open with DEXTER |
Using revised Hipparcos parallaxes and extensive numerical simulations that take into account de properties of the Hipparcos catalogue in terms of completeness and parallax error, and also taking into account the selection functions that are involved in selecting stars with certain properties (K-band, high-quality I-band, spectroscopic metallicities), the absolute magnitude of the red clump in I and K is investigated.
![]() |
Figure 4: Fit of MI versus metallicity and (V-I) colour for model 8, with best fits indicated by the dashed line. The slope of the fit against colour is not significant. |
Open with DEXTER |
Table 4: The RC in old clusters.
In the I-band, the value of
essentially
agrees with the value quoted by Stanek & Garnavich (1998).
There is no dependence of MI on (V-I) colour, and
the dependence on metallicity is marginal:
([Fe/H] + 0.15).
In the K-band, the influence of selection biases appears stronger and
the value derived
differs significantly
from the value originally derived by Alves (2000), which is almost 0.1 mag brighter,
when his value of -1.61 on the Bessell & Brett (1988) system is
converted to the 2MASS system (Carpenter 2001). There is no
dependence of
on metallicity, and the dependence on colour
is weak:
(
(V-K)0 - 2.32).
The RC stars in clusters have also been used to derive MK,
originally by Grocholski & Sarajedini (2002). Absolute magnitudes are
derived by combining the observed mean magnitude in a box in colour
and magnitude with a reddening estimate and main-sequence fitting
distances. Grocholski & Sarajedini (2002) used the second
incremental data release of 2MASS and 14 clusters to find an absolute
K-magnitude that is in agreement with the value in Alves (2000).
Using the All-Sky data release and increasing the sample to 24
clusters, Van Helshoecht & Groenewegen (2007; hereafter vHG) find a
value fainter by 0.05 mag, or
on the
Bessell & Brett system. Both Grocholski & Sarajedini (2002) and vHG
used averages over the cluster sample to arrive at the quoted means.
To compare these values to the absolute magnitude of the local Hipparcos sample, a ``population correction'' has to be made, to account for the difference in metallicity and age of the RC population. The calculation of this correction is outlined and tabulated in Girardi & Salaris (2001) and Salaris & Girardi (2002). In the K-band, Salaris & Girardi (2002) show that the correction is a strong function of age for ages below 3 Gyr. The top panel in their Fig. 3 shows that the correction is well behaved for ages above about 4 Gyr, and for each metallicity a linear function was fitted.
From vHG the five clusters older that 4 Gyr were taken and the
population correction was determined by linear interpolation in age
and [Fe/H] in the results in Salaris & Girardi. The error in
is based on a 1 Gyr error in age and 0.1 dex in
metallicity. The results are listed in Table 4,
which lists the age and metallicity (see vHG for the references for
age and metallicity determination), the MK value and its
error that takes into account the error in reddening and the assumed
distance to the cluster (from vHG), the population correction with
error, and the corrected value.
The weighted mean of the corrected M K values is
with a dispersion of 0.2 mag. Correcting to the 2MASS system
(Carpenter 2001) this becomes
.
Although the
dispersion is large, the mean value is also fainter than the straight
mean for the cluster sample and agrees within 2-sigma with the
determination from the Hipparcos sample.
To settle the issue on the importance of the bias and the absolute
K-magnitude of RC stars would require accurate NIR magnitudes of a
100 to a few hundred (cf. Table 2) bright (
)
RC stars.
Given the brightness, this represents a challenge to modern
instrumentation because of saturation.
The absolute magnitudes derived here are in better agreement with theory than previously.
Using a plausible star formation history for the solar neighbourhood, Salaris & Girardi (2002) derived absolute magnitudes of
MI= -0.17
and
MK= -1.584 (on the 2MASS system), which they compared to the Hipparcos-based results quoted in Alves et al. (2002) of
and
(on the 2MASS system), which indicates differences at
the 2-3 sigma level. The new results of
and
agree with theory at the 1-2 sigma level.
Finally, some implications for existing distance determinations using the RC are discussed.
The derived value of the absolute magnitude in the I-band of
is not very different from previously adopted values in the
literature of -0.23 or -0.26, so the impact on derived distances is not so great.
The distance to the LMC quoted by Pietrzynski & Gieren (2002) of
18.501
0.049 (random+systematic error and assuming no population
correction) based on JK observations of two fields in the bar would
become 18.35
0.05 for an assumed population correction of
-0.03 (Salaris & Girardi 2002).
Babusiaux & Gilmore (2005) obtained infrared data on some fields in
the direction of the Galactic bulge. Their distance scale is
tied to one field for which they obtain
MK = -1.72 for an
assumed distance of 8 kpc. Taking
MK = -1.54would result in a distance of 7.4 kpc (0.2 kpc judging from the
information in their Table 2). Adding the error in MK results in a total error of 0.3 kpc.
Using similar type of observations Nishiyama et al. (2006) obtained
7.52
0.1 kpc for
,
which translates
to 7.33
0.14 kpc with the new DM.
This can be compared to a distance of 7.94
0.37 (random)
(systematic) kpc
from Groenewegen et al. (2007) based on K-band monitoring observations of Population- II
cepheids and RR Lyrae stars (distance tied to an assumed LMC distance modulus (DM) of 18.50) and
the direct determination of
kpc based on the orbit of
the star called S2 around the central black hole when post-Newtonian
physics is taken into account (Eisenhauer et al. 2005; Zucker et al. 2006), although neglect of the space motion of Sgr A* may
result in systematic errors of about 0.1-0.45 kpc (Nikiforov 2008).
All these independent observation are consistent and seem to indicate a
DM to the GC of slightly less than 8 kpc.
Acknowledgements
This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France.