A&A 424, 125-132 (2004)
DOI: 10.1051/0004-6361:20040552
S. Demers 1 - P. Battinelli 2 - B. Letarte3
1 - Département de Physique, Université de Montréal,
CP 6128, Succursale Centre-Ville, Montréal,
Québec H3C 3J7, Canada
2 -
INAF - Osservatorio Astronomico di Roma,
Viale del Parco Mellini 84, 00136 Roma, Italia
3 -
Kapteyn Astronomical Institute, Posbus 800, Groningen, 9700 AV,
The Netherlands
Received 29 March 2004 / Accepted 26 May 2004
Abstract
We present wide field CFH12K observations of the highly
reddened Local Group galaxy IC 10. Using R, I, CN and TiO filters we
identify 676 carbon stars in the field of IC 10. After mapping the
reddening in the field, from the colours of the G dwarfs seen along the
line of sight, we determine the mean apparent magnitude of the
C star population to be
,
leading to a
true modulus of
(
kpc). The old
red giant stars define an asymmetric halo.
With a halo diameter of at least 30', IC 10 is among the largest
dwarfs of the Local Group.
The surface density of C stars
follows a radial power law with a scale length of 2.36', a value nearly
identical to the scale length defined by the old giants thus IC 10 has
a stellar halo where the old and intermediate-age populations are
well mixed.
Key words: galaxies: individual: IC 10 - stars: carbon - galaxies: local group - galaxies: structure
According to van den Bergh's (2000) interesting Local Group review,
IC 10 was recognized as an extragalactic object by Mayall (1935). One
year later Hubble (1936) suggested that it might be a Local Group
member. However, his other two suggested members IC 342 and NGC 6946, never
joined the Local Group club. IC 10 is located at
thus
is it not only
dimmed considerably by interstellar extinction but its outlines
are confused by the high foreground stellar density. This unfortunate
fact handicapped so much the optical investigations of IC 10 that during
the first decades, following its recognition, only one photometric result
appears.
De Vaucouleurs & Ables (1965) obtained its integrated
magnitude (
)
and determined its reddening (
E(B -V) = 0.87)
from its integrated (U-B) and (B-V) colours. During that period, IC 10 was
targeted by radio astronomers,
Huchtmeier (1979) found that the HI envelope of IC 10
reaches
,
a size much larger than its optical image whose
dimension are
(Massey & Armandroff 1995).
In this
respect, IC 10 and NGC 6822 are two Local Group irregulars with a
huge HI envelope, like NGC 3109 which belongs to the nearby
Antlia - Sextans group (van den Bergh 1999).
Table 1 presents a 40-year compilation of the published distance and reddening estimates. This list shows that recent distances range from 500 kpc to 1 Mpc. The reddening E(B -V) inside of IC 10, determined from its Population I objects, seems to be larger than 1.0, while the foreground reddening is closer to 0.80.
Our Local Group carbon star survey has shown (Demers et al. 2003) that the mean I magnitude of late C stars can be used as a standard candle. Thus we can adopt this approach to obtain the distance of IC 10, providing that the C star apparent magnitudes can be corrected for the interstellar extinction. To achieve this goal, the reddening toward IC 10 is evaluated from the whole CMD of the observed area.
Table 1: Historic summary of IC 10 distance and reddening determinations.
The data distributed by the CFHT have been detrended. This means that the images have already been corrected with the master darks, biases, and flats. Fringes have been removed on I exposures under 60 s and large scale structures such as the "Skyring'' effect have been removed when relevant. This pre-analysis produces 12 CCD images, of a given mosaic, with the same zero point and magnitude scale.
We carefully combine images taken with the same filter, making sure that the final FWHM was similar to the average FWHM. The photometric reductions were done by fitting model point-spread functions (PSFs) using DAOPHOT-II/ALLSTAR series of programs (Stetson 1987, 1994). Details of the photometric calibration have been described in recent papers of this series (see, for example Demers et al. 2003).
![]() |
Figure 1: Colour-magnitude diagram of the whole CFH12K field. Stars with colour errors <0.10 are plotted. A substantial differential reddening is responsible of the fuzziness of the features. |
Open with DEXTER |
![]() |
Figure 2: Colour-magnitude diagram of a representative cell showing better defined features. The line corresponds to the mid-ridge, see text for details. |
Open with DEXTER |
In order to investigate the differential reddening across the IC 10 field,
we divide the whole field into 24 cells of 2k 2k pixels (smaller
cells were often found to have too few stars in each magnitude interval
for robust statistics). The widening of the sloping feature visible in the CMD
of Fig. 1 is more
clearly seen in Fig. 2 which corresponds to a representative cell, outside
of the core of IC 10. This widening is undoubtedly produced by the combined
effects of differential reddening
and photometric errors which become larger toward the limiting magnitudes.
The broken line
corresponds to the mid-ridge, defined as the peak of the colour distribution
at various magnitude intervals.
This ridge, which can be easily followed down to
,
represents Galactic G dwarfs located on their main sequence turnoff
seen at increasing distance and reddening along the line
of sight. These stars
have intrinsic colour
(R-I)0 = 0.35 (Cox 2000).
For I = 21, with a reddening of
0.60 a G dwarf of
MI = + 3.6 would
be located at
14 kpc. At this distance, the line of sight in the
direction of IC 10, would be some 800 pc above the Galactic plane.
Is it reasonable to interpret the
slope of the ridge as a continuous increase in reddening, up to this distance
above the plane? The Galactic longitude of IC 10 is 118
,
we are thus looking at some 25 kpc from the Galactic center.
IRAS observations (of warm dust) yield a scale height of 120 pc for the exponential dust layer, see Deul & Burton (1992) for review. Subsequent observations at longer wavelengths by COBE/DIRBE (Sadroski et al. 1994), ISO (Alton et al. 1998) and submillimetric observations (Alton et al. 2001) have shown that spiral galaxies have extended distribution of cold dust which reaches further out than the stellar disk. Optical and near infrared modeling of edge-on spirals by Xilouris et al. (1999) reveals that the dust scale heights are about half the stellar scale heights. Furthermore, the modeling of NGC 891 (Xilouris et al. 1998), a spiral quite similar to the Milky Way, requires two stellar and two dust disks to fit the surface brightness of the galaxy. It is then reasonable to assume that the Galaxy possesses two layers of dust, one corresponding to the thin stellar disk and a more tenuous one corresponding to the thick disk. The stellar scale heights of the thin and thick stellar disks have been evaluated by Méndez & Guzmán (1998) from star counts in the Hubble deep field. They quote 250 pc for the thin disk and two solutions for the thick disk: 1300 pc with a 2% contribution or 750 pc with a 6% contribution.
We determined that the observed reddening increase
with apparent magnitudes is compatible, within 0.05 mag, with a
single exponential dust layer with a scale height of 250 pc. We stress
that better solutions are certainly possible, however we did not investigate
multi-layer solutions because the number of free parameters involved makes
this exercise unproductive.
Furthermore, the likelihood that the scale height
varies with galactocentric distances introduces a serious complication to the
modeling.
In order to reduce, as much as possible, the reddening discontinuity
between the cells we also evaluate the reddening in 15
k cells
centered on the intersections of the grid. We deredden each star by
calculating the weighted mean distance of the star to its nearest three
cells,
we then calculate its extinction from this reddening.
We summarize the reddening determinations by presenting a reddening
map in Fig. 3. The reddening varies, in a irregular way,
by as much as 0.30 mag across the field. The mean reddening in the
field is
E(R-I) = 0.63 which translate into
E(B-V) = 0.79 and
AI = 1.51, using the reddening law by Cardelli et al. (1989) adjusted
for the CFHT Mould's filters.
This average colour excess represents the reddening of the periphery
of IC 10 since its estimate is based mostly on stars away from the
central core of IC 10.
It agrees well with some recent determinations, listed in
Table 1. The application of these reddenings and extinctions yields
the reddening corrected colour-magnitude diagram displayed in Fig. 4;
it is evident how the sloping ridge is now narrower than that shown
in Fig. 1. The slope of the ridge is still present because by applying
the maximum reddening we over-correct the nearby brighter stars.
![]() |
Figure 3:
The mean E(R-I) in each of the
![]() |
Open with DEXTER |
![]() |
Figure 4: The reddening corrected colour-magnitude diagram of the whole region. |
Open with DEXTER |
Carbon stars are selected from their location in the
(CN-TiO) vs. (R-I)0 plane. The colour-colour diagram is displayed in Fig. 5.
Throughout this series of papers we have selected as an adoption criterion,
stars with
,
where
.
Over 66 000 stars in the field satisfy this criterion.
The zero point of the (CN-TiO) index is set according to the procedure
outlined by Brewer et al. (1995). We set the mean of
(CN-TiO) = 0.0 for all blue stars since hot stars are expected to have
a featureless spectra in the CN and TiO regions. We thus define a blue star,
as in Letarte et al. (2002), i.e. as a star
in the color range
0.0 < (R - I)0 < 0.45. The CN-TiO approach
requires a previous estimate of the local colour excess.
![]() |
Figure 5: Colour-colour diagram to identify C stars from M stars. |
Open with DEXTER |
The upper branch of the colour-colour diagram shows that IC 10 contains several hundreds C stars. Since the CN-TiO technique confuses M and C stars for spectral types earlier than M0, we restrict our definition of C stars as those with (R -I)0 > 0.90, a colour corresponding to M0: 676 stars are in the C star box. The equatorial coordinates (J2000), the dereddened magnitudes and colours of these stars are listed in Table 2. It is, however, evident from Fig. 5 that in addition numerous bluer C stars are also present in IC 10.
Among the C stars retained there is a bright star with I0 = 16.33 (I=17.86). This star is well outside the central core of IC 10. It is however very close to an even brighter star, saturated in our images, whose profile may be overlapping with the C star. The photometric errors of the C star are abnormally high for its magnitude, they are barely below the threshold of acceptability. A search through the 2MASS database yields a star, 3.4'' from it, with K = 12.125 and J-K = 1.07. It is not clear if this corresponds to the bright star or to our C star. Nevertheless, the 2MASS source is too blue to correspond to a C star, as we define them (Demers et al. 2002).
Table 2: Full table is available at the CDS. A portion is shown here for guidance regarding its form and content. Units of right ascensions are hours, minutes and seconds, and units of declination are degrees, arcminutes and arcseconds.
The luminosity function of the 676 C stars, located within the C star box
of the colour-colour diagram, is presented in Fig. 6.
A Gaussian, with a
mag. is eye fitted to the distribution.
The mean magnitude of these C stars is
while
their mean colour,
.
Demers et al. (2003)
have shown that the
of C stars in galaxies having
more than one hundred of them, shows very little dispersion. Since that paper
we have observed Wolf-Lundmark-Melotte (Battinelli & Demers 2004a), NGC 185
(Battinelli & Demers 2004b) and NGC 147 (Battinelli & Demers 2004c). Taking
into account these new data we obtain, for galaxies with more than 100 C stars:
.
The uncertainty attached to
depends
on the uncertainty of the zero point of the photometry and on the accuracy of the
adopted mean E(R-I). Assuming that both of these are less than 0.05 mag leads
to a true modulus of
,
corresponding to 741 kpc.
This distance places IC 10 at about 240 kpc from M 31. The absolute
magnitude of IC 10 is then
Mv = -15.6, adopting de Vaucouleurs &
Ables (1964) integrated apparent magnitude. Its luminosity is similar to NGC 3109 but brighter than NGC 6822 or IC 1613.
![]() |
Figure 6: Luminosity function of C stars. A Gaussian is eye fitted to the histogram. |
Open with DEXTER |
IC 10 has been often called the nearest starburst galaxy. It is even
considered a "blue compact dwarf'' by Richer et al. (2001). The recent
planetary nebulae survey by Magrini et al. (2003) has however revealed that IC 10 must possess a much larger halo than the starburst core.
Our C star survey confirms this finding and extends even further the
distribution of the stellar halo. Figure 7 presents a map of the
C stars and of the 16 known planetary nebulae. We see C stars over most
of the CFH12K
field.
The obvious central hole in the distribution corresponds to the densest
core of IC 10 which can be seen in the I-band image published by
Sakai et al. (1999). Note, however that their East-West orientation is
reversed to ours.
The horizontal gap and a less evident vertical one
are due to
8'' gaps between CCDs.
![]() |
Figure 7: The CFH12K field of IC 10. The 676 C stars are shown by solid dots while the known planetary nebulae are represented by circles. North is on top, East is on the left. |
Open with DEXTER |
In order to determine the size of IC 10 we can use two tracers to study
the surface density behaviour: 1) the bulk of the IC 10 stars, seen
in the CMD of Fig. 4 or 2) C stars identified in Sect. 3.2. While the former sample
consists in a much larger number of stars (20 000 stars)
the latter has the advantage of being free from foreground contamination.
One of the major problem is to determine the center of the IC 10 because the galaxy is highly patchy near the center. For the following analysis we assume that: a) the center of the IC 10 corresponds to the coordinates given in NED; b) both the bulk of red stars and C stars have a spherical distribution around this center.
Having determined the reddening and the distance of IC 10 we are
now in position to establish the evolutionary status of the the
bulk of its giants seen in our CMD. Since the oxygen abundance of IC 10,
determined from its HII regions by Garnett (1990), is log(O/H
) + 12 = 8.26we select Girardi et al. (2002) isochrones for 10 Gyr and Z = 0.004.
The isochrone, as can be seen on a close-up of the CMD, displayed in
Fig. 8, does not match the RGB location. A lower abundance, z = 0.001,
fits better the bulk of the observed RGB. Hunter (2001) has determined the
mean colour of the TRGB to be
(V -I)0 = 1.74, this corresponds to
(
R -I)0 = 0.77 using Battinelli et al. (2003) relation. Therefore,
z = 0.001 is obviously much more representative of the old population abundance.
The parallelogram outlines the 15 000 stars we select as
representative of the bulk of the old stellar population seen in our images.
![]() |
Figure 8: A close-up of the CMD presented in Fig. 4 with 10 Gyr isochrones with different abundances. The parallelogram defines the sample of red giants. |
Open with DEXTER |
The spatial distribution of these RGB stars is shown by isodensity contours displayed in Fig. 9. They clearly reveal how complex and irregular the structure of IC 10 is. In particular, all the northern half is highly asymmetric, in the inner as well as outer regions. Such irregularities cannot be ascribed to an extra dust absorption not properly accounted by our dereddening procedure (described in Sect. 3.1) since the average (R - I)0 for the C stars in the Northern and Southern halves of the galaxy are identical. Certainly, at least in the NW quarter of IC 10, the presence of a cospicuous population of young bright stars may prevent the detection of the fainter RGB component thus resulting into an artificial lower surface stellar density. On the other hand, the southern half appears much more regular and the contours, at least in the outer regions, are well matched by circular arcs centered on the adopted IC 10 center. These arcs demonstrate that the adopted center and the shape are well justified. We have to stress however that it is not possible to judge if a slight ellipticity is present.
![]() |
Figure 9:
Isodensity contours of red stars in IC 10 (see text)
over the
![]() |
Open with DEXTER |
The RGB stars in the box of Fig. 8 are used to derive the
surface density in circular annulii, for the southern half of the galaxy.
This profile is shown in Fig. 10.
IC 10 seems to extend
up to 16' where a clear plateau, corresponding to the foreground,
is reached. The average density of the last 8 points leads to a foreground
density of
arcmin-2.
A radius of 16 arcmin, at the distance
of IC 10 corresponds to a diameter of 7 kpc, making IC 10 one of the
largest dwarf irregular of the Local Group.
![]() |
Figure 10: Density profile (stars per arcmin2) from RGB star counts in the southern half of IC 10. The dotted line indicates the foreground density plateau. |
Open with DEXTER |
From a least square fit of the
foreground corrected logarithmic counts (see Fig. 11) we estimate a radial scale length
of
arcmin (560 pc).
We remind again that such scale length refers only to the southern half of IC 10 and
that the galaxy is much more irregular and possibly is less extended on the northern
side
(see Fig. 9).
The mean X and Y coordinates of the 676 C stars correspond to a point 35'' off the presumed center, again justifying our adopted center.
Since the foreground contamination is essentially nil for C stars we can derive
the radial scale length directly from the observed counts. Similarly to the
previous subsection, we limit our analysis to the southern half of the galaxy.
The radial
surface density profile is displayed in Fig. 11.
The surface density, from 4' to 15'follows closely a power law. A least square solution through these points
yields a radial scale length of
.
This scale length, which
is in fair agreement to what obtained in the previous subsection, implies
that the old and the intermediate-age populations are - at least up to a radius
of 15 arcmin - spatially distributed
in the same way.
![]() |
Figure 11: Radial profile of the surface densities of old giants (filled circles) and C stars (squared dots), For an easy comparison 2 has been added to the logarithmic density of C stars. Surface density units are star arcmin-2. Lines are least square fits to the points in the radial range 4' to 15'. |
Open with DEXTER |
To select AGB M stars among the stars in the M-box of the colour-colour diagram
in Fig. 5, we follow the procedure adopted by Battinelli
et al. (2003) which consists in rejecting stars with
fainter
than the TRGB luminosity as well as those brighter
than the brightest members of the galaxy. The bright cutoff for the
selection of AGB M stars is
I0 = 18.5which corresponds
to the value used by Battinelli et al. (2003) for M 31 scaled for the
the different distance and absorption of IC 10. The whole CFH12K
area contains 11 204 such AGB stars. We evaluate the foreground contribution
by assuming that two 5' wide strips at both extremities of the field
are essentially free of IC 10 stars. From these counts we estimate the total
foreground contribution to be
stars resulting in
M AGB stars belonging to IC 10. The C/M ratio is then
.
A comprehensive discussion of the
C/M ratios as a metallicity indicator will be presented in the next paper of
this series.
Acknowledgements
This research is funded in parts (S.D.) by the Natural Science and Engineering Council of Canada.