text:
# Long Period Variables detected by ISO in the SMC

## 1 Introduction

The stellar content of the Small Magellanic Cloud (SMC) can be studied in great detail because of the low galactic extinction and a known distance m$``$M$`=18.99\pm 0.03`$(formal)$`\pm 0.08`$(systematic), Cioni et al. (cionitip (2000)). However, it is well known, but not yet quantified, that the SMC is extended along the line of sight much more than its companion galaxy, the Large Magellanic Cloud (LMC). Together with our own Galaxy, they constitute an interacting system (e.g. Gardiner & Noguchi garden (1996), Moore & Davis moore (1994), Heller & Rohlfs heller (1994)) and a metallicity scale environment in which the Galaxy is the most metal rich and the SMC the most metal poor of the three. Fields covering almost the whole of the two Magellanic Clouds and a considerable number of fields in the Galaxy have been repeatedly observed by the micro–lensing projects: OGLE, MACHO and EROS. This is providing a wealth of data to identify and study variable stars in every evolutionary phase. Partial or all sky survey data in the near–infrared (IR) provide essential counterparts for the specific analysis of the reddest objects: red (RGB) and asymptotic giant branch (AGB) stars. Most AGB stars are pulsating with an amplitude of up to a few magnitudes in the optical wave bands and somewhat less in the near–infrared bands. Long Period Variables (LPVs) have a period of at least $`100`$ or $`150`$ days and longer. Many are also multi–periodic objects. More than one dependence between period and luminosity has been detected in the LMC (Wood et al. woodiau (1999), Cioni et al. cionivar (2001), Noda et al. moa (2002), Lebzelter et al. agape (2002)). Each period–luminosity relation (PLR) is probably associated with a different mode of stellar pulsation. Because most of these stars are probably experiencing the thermal pulsing AGB phase (TP–AGB) they are losing mass at an average rate of $`10^6`$Myr-1 and they have highly extended atmospheres. The pulsation as the driving mechanism of the stellar wind is a fundamental element in the study of the mass–loss rate of AGB stars. They are surrounded by an initially optically thin, and later optically thick, circumstellar envelope and thus sometimes become detectable only in the infrared wave bands. Many obscured AGB stars were discovered by the IRAS satellite. More recently the Midcourse Space Experiment (MSX - Price et al. msx (2001)) observed selected high density regions in the Magellanic Clouds. Improved source position were derived. Egan et al. (egan (2001)) present the cross-identification between MSX and 2MASS data in the LMC. However, the Infrared Space Observatory (ISO) obtained more sensitive data in some fields in the Magellanic Clouds and the Galaxy that increased the number of obscured AGB stars and allowed us to learn in more detail about their dust properties.

In order to study the evolved AGB population it is thus important to utilize observations that cover the spectral range from the optical ($`I`$ band) to the mid–IR. It is then possible to derive fundamental quantities such as the bolometric magnitude and the mass–loss rate of the stars. In addition, from the period, the amplitude and the regularity of the variation, it is possible to characterise the stars as Mira or Semi–Regular (SR) pulsators and to indicate, from the PLR, the mass of the stellar progenitor. In this article the properties of a sample of AGB stars in the SMC extracted from the ISO–Mini-Survey catalogue of the Magellanic Clouds (ISO–MCMS, Loup et al. in preparation), which already contains the cross–identification with the near–IR catalogues DENIS and 2MASS, are discussed. This is done by additionally analysing the light–curves extracted from the MACHO database. A similar work was performed in the LMC bar–west field (Cioni et al. cionivar (2001)) and will be further extended to other fields were ISO measurements are now available. A detailed comparison of more or less obscured AGB variables between the two Clouds and the Baade’s window NGC6522 will be given in another paper (Glass, Schultheis & Cioni, in preparation).

The data and the cross–identification with the light–curves are described in Sect. 2. In Sect. 3 the analysis of the period and amplitude of variability and the determination of the bolometric magnitude is presented. A discussion of the different type and variability class of the sources in the PLRs and in the near– and mid–IR colour–magnitude diagrams is given in Sect. 4. Sect. 5 concludes the article. Appendix A discusses the re–analysis of the LPVs studied by Feast et al. (feastpl (1989)) and B compares the bolometric correction obtained in this work with the relations by Alvarez et al. (alva (2000)) and Montegriffo et al. (monte (1998)).

## 2 Data

The sample of data in the SMC, analysed here, is based on the cross–identification between the ISO–MCMS (Loup et al. in preparation), and the MACHO light–curves publicly available111http://wwwmacho.mcmaster.ca.

The ISO mini–survey in the direction of the SMC covers an area of $`0.28`$ square degrees. Observations were performed with the infrared camera ISOCAM on board the ISO satellite in the LW2 ($`6.61\mu `$m) and LW10 ($`11.18\mu `$m) wavebands. The catalogue lists $`1333`$ sources brighter than about $`12.4`$ and $`11.3`$ mag in LW2 and LW10, respectively. The coordinates of the ISO sources are better than $`4`$″ with a faint tail up to $`6`$″. The catalogue also lists the cross–identification with the near–infrared catalogues 2MASS (2nd Incremental Release PSC222http://www.ipac.caltech.edu/2mass/) and the DENIS (Epchtein et al. nicolas (1997)) catalogue towards the Magellanic Clouds (DCMC – Cioni et al. 2000a ), and a few smaller optical catalogues of carbon stars, HII regions, PNe and emission–line stars. According to the association likelihood criteria defined by Loup et al. (in preparation) $`48\%`$ of the sources have no, or an unlikely, association with DENIS/2MASS; $`13\%`$ have a questionable association; $`25\%`$ have a confident association and the remaining $`14\%`$ are “known” sources with a DENIS/2MASS and optical counterpart. It is perhaps surprising that there are many ISO sources without a near–infrared counterpart. However, most of these sources are among the faintest ISO detections in both ISO wavebands. The histogram of their magnitudes and colour peak at: $`LW2=12`$, $`LW10=11`$ and $`(LW2LW10)=1.5`$. Moreover, the ISO–MCMS includes sources well detected in one band only. Among sources detected only in two bands about $`86\%`$ have a near–IR counterpart. The ISO photometric magnitudes are part of the final version of the catalogue released to the Co–Is and associates on March 2002.

The MACHO (Alcock et al., alcock (2000)) light–curves were obtained from repeated measurements in two, blue and red, broad pass–bands from the 50 inch telescope at the Mount Stromlo Observatory in Australia. Observations of fields covering most of the Magellanic Clouds were performed from $`1992`$ to $`2000`$. The light–curves are easily accessible via the web interface giving in input the right ascension (RA) and the declination (DEC) of each source.

Fig. 1 shows the regions observed by ISO superimposed on the MACHO fields. The area covered by ISO is entirely contained in MACHO fields number $`211`$ and $`212`$; a few sources are also present in field number $`207`$. These MACHO fields were observed for the whole period of eight years while the observations of fields in the external parts of the SMC cover just five years. The OGLE light–curves cover a period of about five years (observations started in 1997). Because of the nature of LPVs, sources with periods from $`100`$ to $`1000`$ days, the MACHO light–curve database was preferred.

### 2.1 Match between ISO–MCMS and MACHO

Each ISO source with a confident DENIS/2MASS association has been searched for a light–curve in the MACHO database using, as input, the 2MASS coordinates. Variables were confidently identified within a radius of $`1`$′. Fig. 2 shows that most of the sources were actually associated within $`1`$″. In a very few cases the MACHO counterpart, a source with a typical LPV light–curve, was not the closest to the 2MASS source – but still within no more than $`4`$″. We consider a light–curve to be that of a LPV when a more or less periodic behaviour is detected by eye. It is expected that most of the sources detected in the near– and mid–IR are associated to AGB stars. This is a consequence of the mass–loss process occurring in AGB stars that builds up a circumstellar envelope observable in the mid–infrared.

Among the $`502`$ sources with a confident DENIS/2MASS counterpart, $`458`$ have been found to have a MACHO light–curve. Table 1, electronically available, lists the MACHO coordinates at the epoch J2000 (Column 1 & 2), the MACHO identifier field, tile and number (Column 3, 4 & 5), the distance in arcsec between the 2MASS and the MACHO source (Column 6) and the ISO–MCMS identifier (Column 7). The first $`10`$ lines are given as an example (Table 1).

The sources listed in Table 2 are those for which we could not find a convincing MACHO counterpart. The ISO–MCMS name is given in Column 1, the distance in arcsec to the first possible MACHO counterpart is given in Column 2 and in some cases the reason of the missing cross–identification is written in Column 3. Too few data points means that there are not enough points to characterise the shape of the light–curve. This would be useful to recognize a MACHO counterpart within the ISO astrometric precision.

## 3 Period analysis

The Schwarzenberg–Czerny method (periods (1996)) has been used to detect periodicities in the light–curves of the selected sources, independently for the measurements in the blue and red bandwidth of MACHO. The algorithm was executed thirty times searching for periods in a range from $`10`$ to $`100`$, $`200`$, …, $`3000`$ days with a step of $`0.005`$ days. The result of the program is the period of the variation and a few other parameters among which the quality of the variability “q”. We identified as the main period the value corresponding to the highest value of “q” among the thirty different determinations. We noticed that at the first passage through the periodic detection algorithm the dominant period did not give, in some cases, a sufficiently defined phase–curve (distribution of the measurements folded on a given period). By restricting or enlarging the range of searching periods it was always possible to optimize the phase curve to a certain period. Therefore we decided to execute the algorithm more times by changing the period searching window. The reliability of the assigned period was then inspected by looking at the phase curve of each source. Sources of regular variability and with a large amplitude of variability have a clean phase plot with a large S/N ratio that allows the variation of the position of the maximum or minimum brightness and/or the presence of bumps to be distinguished (sometimes identified as the secondary period). In other cases the S/N was not good enough to detect these secondary variations but the phase–curve was clear enough to confirm the main period. Then we introduced the following criteria to flag the sources in each band. Flag$`=0`$ for sources that do not have data points in a given MACHO band. Flag$`=5`$ for sources with $`q15`$ or a period that coincides with the extremes of the searching window in all $`30`$ cases (i.e. sources with a period of $`10`$ or $`100`$, $`200`$,… $`2900`$, $`3000`$ days). Their light curve is judged to be of poor quality. Flag$`=1`$ and Flag$`=2`$ for sources with $`q>15`$ and for which the assigned period (P) differs not more than $`10`$ days from the mean of all the measurements (among the thirty possible) with $`q>15`$ in the blue and red band width, respectively. With these criteria we are confident that a single periodicity is the dominant source of variation in the light–curve. Flag$`=6`$ for sources with $`q>15`$ and a difference from the same mean larger than $`10`$ days. These sources are likely to have more than one significant periodicity. In this work we are not aiming at a detail analysis of the secondary periods which will be discussed, together with the comparison of LPVs with an ISO counterpart between the Magellanic Clouds and the Galaxy , in Glass, Schultheis & Cioni (in preparation). Aliases were identified from the diagram $`Log(P)`$ versus $`K_S`$ as those periods creating clear vertical patterns. These correspond to periods equal to $`29`$ days and in the range from $`340`$ to $`390`$ days, included. These sources were assigned Flag$`=9`$. Examples of MACHO light–curves for LPVs can be found in Wood et al. (woodiau (1999)). Table 3 summarises the detected periodicities. It lists the description of the Flag (Column 1), the value of the flag and the number of sources corresponding to a given flag in the red (Column 2 & 3) and blue bandwidths (Column 4 & 5).

Table 4 provides the magnitudes from the different databases, the bolometric magnitude (see next Section) and the parameters that describe the dominant pulsation of all the sources in the sample with an extracted light–curve ($`458`$). This table is also made electronically available and the first $`10`$ lines are given as an example (Table 4). The content of the table is the following: ISO–MCMS name (Column 1), $`IJK_S`$ DCMC magnitudes (Columns 2, 3 & 4), $`JHK_S`$ 2MASS magnitudes (Columns 5, 6, & 7), LW2 and LW10 ISO magnitudes (Column 8 & 9), the apparent de–reddened bolometric magnitude (Column 10), period, amplitude, flag and quality parameter derived for the light–curve in the MACHO–red bandwidth (Column 11, 12, 13 & 14) and in the MACHO–blue bandwidth (Column 15, 16, 17 & 18).

Values of $`99.99`$ for the magnitudes and of $`99`$ for period and amplitude indicate missing data. Delmotte et al. nausicaa (2002) have shown that $`J_{DCMC}=J_{2MASS}(0.11\pm 0.06)`$ and $`K_{sDCMC}=K_{s2MASS}(0.14\pm 0.05)`$ but we did not apply this correction to Table 4. Amplitudes (A) are the difference between the minimum and the maximum value of the MACHO intensity and are expressed in $`mag\times 100`$. In a handful of sources this calculation might be biased by spurious data points with a magnitude excursion compared to the sinusoidal behaviour of the light–curve. Perhaps these points are real and represent a flare. The period is in days.

### 3.1 Bolometric Magnitude

With the availability of magnitudes from the $`I`$ band to the mid–IR band LW10, it is possible to determine quite precisely the bolometric magnitude of the sources in the sample by integrating under the spectral energy distribution (SED). This integration has been done by fitting a spline to the SED profile in the plane $`(\nu ,F_\nu )`$, extrapolating linearly to zero flux at zero frequency, and also extrapolating linearly at the other extreme of the SED through the two near most points. The SED of all the sources detected by DENIS and ISO were inspected one–by–one to exclude those which are either too red or too blue to describe satisfactorily their SED with the available measurements.

The adopted extrapolations represent the SED very well, almost independently from the circumstellar opacity for all M spectral subtypes (Loup et al. in preparation). That article gives a detailed determination of the multi–band bolometric correction for AGB stars by combining data from the major survey instruments: DENIS, 2MASS, ESO, IRAS, ISOCAM and MSX. The bolometric correction is obtained from the different combination of pass–bands and a theoretical model for the SED. This is compared with the correction obtained from the integration under the SED using a spline or a linear fit through the flux points. A fairly well sampled SED, that detects most of the energy emitted by the objects in a set of characteristic bands, is a necessary condition to derive good bolometric luminosity ($`M_{BOL}`$) for stars of different M spectral subclass and C stars.

The best wavelength coverage for the data in the sample studied here is obtained by combining DENIS and ISO magnitudes. Though 2MASS $`K_S`$ magnitudes exist for $`455`$ sources, DCMC–$`JK_S`$ magnitudes exist for $`366`$ sources and DCMC–$`IJ`$ mag exist for about $`450`$ sources. On one hand it would seem natural to combine $`IJ`$ DENIS magnitudes with 2MASS–$`K_S`$ magnitudes and the ISO magnitudes, for the same sources, to derive the bolometric magnitude. On the other hand the DCMC–$`IJK_S`$ magnitudes were obtained by simultaneous observations and this is a unique constraint to the SED profile.

Note that the amplitude of variation in the mid–IR bands can be as high as $`1`$ mag (e.g. van Loon et al. jacco98 (1998), van Loon et al. jacco01 (2001)), this effect may influence the shape of the SED during the variability cycle. The constraint of the $`I`$ band to the shape of the SED at higher frequencies is very important for optically thin AGB stars because the peak of the SED is at about $`K_S`$. In Loup et al. (in preparation) we show that the linear extrapolation to zero flux at the longer frequencies starting from the flux in the $`J`$ band subtends a larger area than the extrapolation starting from the flux in the $`I`$ band. This effect overestimates the bolometric magnitude.

Prior to the calculation of the flux emission in the different bands and the integration under the SED we de–reddened the data adopting the extiction law by Glass (gext (1999)). For the DENIS pass–bands $`[A_V:A_I:A_J:A_{K_S}=1:0.592:0.256:0.089]`$, and we assumed zero absorption in the ISO wave bands. We used $`R_V=3.1`$ and $`E(BV)=0.065`$, the latter is the average of the measurements, in the SMC, discussed in Westerlund (westbook (1997)). The zero–point of the bolometric scale is computed with: $`M_{BOL,\mathrm{}}=4.74`$ and $`f_{tot,\mathrm{}}=1.371\times 10^6`$ erg cm-2 s-1. A comparison of our bolometric magnitudes with the relations obtained by Alvarez et al. (alva (2000)) and Montegriffo et al. (monte (1998)) is given in Appendix B.

## 4 Discussion

### 4.1 Period–luminosity relation(s)

Fig. 3a shows the relation between the logarithm of the period (days) and $`K_S`$. Variable stars are distributed on four different parallel relations. These relations were first discovered by Wood et al. (woodiau (1999)) from the analysis of a sample of MACHO stars in the LMC. They were subsequently found by Cioni et al. (cionivar (2001)) using EROS2 (Lasserre et al. lass (2000)) light–curves, by Noda et al. (moa (2002)) using the MOA database and by Lebzelter et al. (agape (2002)) using the AGAPEROS data. These latter measurements refer also to samples of stars in the LMC. A preliminary comparison between the relations in the LMC and in the SMC has been shown by Cioni et al. (2002a ), concluding that sources in both galaxies occupy the same relations in the diagram ($`K_S`$, $`Log(P)`$).

There are a few differences that depend on the way the two samples were selected and the light–curves extracted: in the SMC there are sources with longer periods because the MACHO data cover a longer time range than the EROS2 data. Thus it is possible to identify longer periodicities. Also there is a relative lack of SMC sources around the tip of the red giant branch (TRGB) because the LMC sources were selected from the DCMC and not among those sources detected at $`7`$ and $`12`$ $`\mu `$m. It may be that AGB stars that are not yet in the TP–AGB phase or are in a very early phase of their TP–AGB evolution have a mid–IR emission below the ISO detection limit. The ISO survey was indeed designed to find mass–loss, not to detect sources without circumstellar extiction.

The analysis of the light–curves of LMC sources detected by ISO has just started. The $`K_S`$ magnitudes are always a one epoch measurement, and coupled with the amplitude of variation, that in the $`K_S`$ band can be as high as $`1.5`$ mag for OH/IR stars, prevent us from detecting a difference (if any) in the width of the PLRs due to metallicity. The same situation holds for stars in the Bulge (Alard et al. alard (2001)), which have the additional disadvantage that the distance to the AGB stars is too uncertain to find all four PLRs. Furthermore, AGB stars in stellar clusters are too few to outline the presence of one and possibly two PLRs. There is evidence that in the elliptical galaxy NGC5128 one and maybe two PLRs exits for LPVs (Rejkuba marina (2002)).

The atmosphere of AGB stars is dominated by molecules of TiO and VO if the star is O–rich (C/O$`<1`$) or by molecules of CN, C2 and other carbonaceous molecules if the star is C–rich (C/O$`>1`$). Stars of both types are easily statistically distinguished in the $`(JK_S,K_S`$) colour–magnitude diagram (Fig. 4a) because of the different strength of the molecular absorption bands at these particular wavelengths (Cioni et al. cionivar (2001), Loup et al. loupbmb (2002)). The dividing line at $`(JK_s)=1.33`$ takes into account the shifts calculated by Delmotte et al. (nausicaa (2002)) between 2MASS and DCMC magnitudes and has been defined in Cioni & Habing (cionicm (2002)) for the whole AGB population of the SMC. It was also formerly checked by Cioni et al. (cionitip (2000)) via the cross–identification with the large sample of spectroscopically identified C–rich stars of Rebeirot et al. (reb (1983)). Filled symbols in Fig. 3 indicate C–rich objects. In Fig. 3b they define relation B and C, and dominate relation D (possibly offset from the oxygen–rich stars).

### 4.2 Mass and Age

According to the Vassiliadis and Wood (vass (1993)) prescription for the evolution of low– and intermediate–mass stars carbon stars are less massive than about $`3M_{\mathrm{}}`$. They form after a certain number of shell flashes when the atoms of carbon overnumber the oxygen atoms. Hot bottom burning may prevent massive stars from becoming C–rich. Each track drawn in Fig. 3b starts at the first shell–flash cycle on the TP–AGB phase. Most of the C–rich stars in the SMC sample have a mass between $`1`$ and $`3M_{\mathrm{}}`$. The brightest stars ($`M_{BOL}<7`$) are probably supergiants (Wood et al. woodsg (1992)). The most massive AGB stars with a mass of about $`5M_{\mathrm{}}`$ are about $`0.1`$ Gyr old. The oldest detected AGB stars are about $`9`$ Gyr old, but most of the sample is from $`0.6`$ to $`2`$ Gyr old.

In the LMC (Cioni et al. cionivar (2001)) C–rich stars occupy only the brightest part of the PLRs, thus in a lower metallicity environment, such as that of the SMC, C–rich stars form at lower masses, unless these stars are in the luminosity dip that follows a thermal pulse or at the minimum of their light–curve. Menzies et al. (menzies (2002)) recently showed that the top magnitude of the AGB stars in $`K_S`$ in the dwarf galaxy Leo I (\[Fe/H\]=-2.0) are mostly carbon stars. They have $`JK_S>1.1`$. The value of the $`JK_S`$ that discriminates between O–rich and C–rich stars is a function of metallicity. The lower the metallicity, the bluer the colour. Obscured AGB star candidates ($`JK_S>2`$) might also be of C type ($`50\%`$ level) and the observation at two different epochs indicates that at least one of these object might be a Mira variable ($`\mathrm{\Delta }K_S=0.87`$).

### 4.3 Large Amplitude Variables

Fig. 3c shows the amplitude of the light–curves as a function of $`M_{BOL}`$ for C–rich and O–rich objects. The sources with $`A100>110`$ and $`f=1`$ have $`q>100`$ and are distributed in the PLR(s) as shown in Fig. 3d. They occupy relation C which is the location of Mira variables in the LMC and the Galaxy. Over plotted is the sequence derived by Feast et al. (feastpl (1989)) and the single LPVs (crosses) as re–analysed in this work (Appendix A). This sets the minimum level of the quality parameter for Mira variables. This sequence was obtained including both O–rich and C–rich AGB stars with a period $`P<420`$ days. A distance modulus of $`18.99`$ has been used (Cioni et al. cionitip (2000)). The single carbon at $`Log(P)=2.1`$ and $`M_{BOL}=6`$ is of large amplitude because of flare–like events super imposed on a regular variation of smaller amplitude, thus this source should not be regarded as a Mira star.

The scatter in the PLR, of the sources in this study, results from by the unavailability of bolometric light–curves. However, it is consistent with the sigma associated with the relation derived by Feast et al. feastpl (1989). Estimating the magnitude of this effect would be valuable for studies that plan to use AGB stars as tracers of the star formation history of systems as far as $`250`$ kpc (Kuĉinskas et al. kucin (2002)). This preliminary investigation on the use of AGB stars concerns mainly early–type AGB stars which vary with amplitudes much below that of Mira–type stars. On the other hand Kuĉinskas et al. combined broad–band observations from the literature. It should not be neglected that the amplitude of variability of AGB stars is different in each photometric band. Only simultaneous observations in different bands, as those that will be available with GAIA, provide a coherent SED from which a single–epoch bolometric magnitude can be derived.

Finally we conclude that a large amplitude and a large value of $`q`$ indicate a clear and regular light–curve; thus a Mira star. There are some sources of large amplitude and with a large value of the quality parameter suspected of multi–periodicity. Their main period does not correspond to relation C. A few Mira variables in sequence D have been found by Alard et al. (alard (2001)). On the contrary, this PLR could be populated by binary stars with an AGB companion (Wood et al. woodiau (1999)). In what follows we mainly discuss the properties of the stars with a light–curve in the red band of MACHO and $`f=1`$ or $`f=6`$.

### 4.4 Mira & Semi–Regular Variables

The distinction between Mira and Semi–Regular (SR) stars has become more and more subtle with the availability of more sensitive long term observations and large samples of stars. The stellar light–curve can be as regular in SR as in Mira stars. Lebzelter et al. (agape (2002)) define a classification criterium for LPVs based only on the regularity and multi–periodicity of the light–curves. Most of the Mira stars fall into the class of regular variables but not all regular variables have large amplitudes of variation. SR belong to a second class and have on average a smaller amplitude than the stars in the first class. Into the third class fall irregular variables, that do not show a clear periodicity. The Mira phase is believed to be just a moment during the AGB evolution and the star can thus become or have been SR. This scenario is supported by observations of very nearby variable stars over a period of about $`90`$ years (i.e. Percy & Bagby perbag (1999), Kiss et al. kiss (2000)). However the evidence, e.g. from globular clusters, suggests that most SRs are in the pre–Mira stage.

Mira variables lie on one PLR (Feast et al. feastpl (1989)) and SR occupy this and other PLRs according to their mode of pulsation (Wood et al. woodiau (1999), Cioni et al. cionivar (2001), Noda et al. moa (2002), Lebzelter et al. agape (2002)). If Mira variables are selected as discussed in Sect. 4.3 then they are clearly concentrated on the C PLR. This selection criterium is based on the amplitude of variation. However, the sources are the most regular as well ($`q>100`$).

Whitelock et al. (wfmo (1997)), by analysing the light–curves of a dozen C-rich AGB stars in the Milky Way Galaxy, showed that the distinction between Mira and SR variables is not necessarily very clear among carbon stars, and may not be useful as a distinction for C–rich stars as it is for O–rich stars. Fig. 4b shows the histogram of the distribution of the amplitude of variation separately for O–rich and C–rich stars. Perhaps there is a distinction between Mira, the pulsators at the tail of the amplitude distribution(s), and SR. Both O–rich and C–rich sources have a histogram with a dominant peak at about $`A\times 100=40`$ and $`A\times 100=80`$, respectively, and a tail to larger amplitudes. The fact that C–rich stars have on average larger amplitudes than O–rich stars is different from what was found in the LMC. There, both histograms had the same distribution (Fig. 4 in Cioni et al. (cionivar (2001))). The same histogram but for the amplitudes derived from the BLUE light–curve shows that O–rich stars peak at about $`A100=50`$ and C–rich stars at about $`A100=100`$. This confirms that the average amplitude is larger in the blue than in the red spectral domain.

### 4.5 Periodicity peaks

The histogram of the distribution of the variable sources versus their period is shown in Fig. 4c. The distribution of O–rich stars has a main peak at roughly $`Log(P)=2.80`$ and two lower peaks at about $`Log(P)=1.95`$ and $`Log(P)=2.30`$. C–rich stars have well defined peaks at about $`Log(P)=2.45`$ and $`Log(P)=3`$. Note that most of the sources with $`2.53<Log(P)<2.59`$ have been excluded because their period determination is suspected of aliasing.

Fig. 4 in Cioni et al. (cionivar (2001)) shows that O–rich AGB stars in the LMC have a strong broad peak between $`Log(P)=1.5`$ and $`2.0`$, and a lower peak at about $`Log(P)=2.45`$, while C–rich AGB stars are homogeneously distributed from $`Log(P)=1.5`$ to $`2.5`$. The first LMC peak of O–rich AGBs is probably a consequence of the selection criterium used to define the two samples. The SMC sample is extracted from sources with an emission in the mid–IR, these stars are brighter and with longer periods than the stars in the LMC sample which were selected from the DCMC database. The difference in strength is probably due to the mean difference in the AGB type between the two Clouds. There are more carbon stars in the SMC which is on average more metal poor than the LMC (Cioni & Habing cionicm (2002)). This effect is also important for the sources populating relation C which are predominantly C–rich in the SMC and associated to the peak at about $`Log(P)=2.45`$. That the majority of accepted Miras in the SMC are carbon stars was first recognized by Lloyd Evans et al. (llglca (1988)) as the result of a search for large–amplitude red variables in the Radcliffe Variable Star Field. The longer time range spanned by the MACHO observations allows to detect periodicities above $`Log(P)=2.8`$ which was about the limit reached in the LMC study using EROS light–curves. The large number of sources detected in the SMC at these periods are likely to be multi–periodic. Their secondary period will be part of relation A, B or C, thus they will enhance the number of O–rich stars between $`Log(P)=1.5`$ and $`2.5`$. In more detail we found that in our sample, among those with a red light–curve and $`q>15`$, $`24`$ sources have $`flag=6`$. These sources have $`2.0<Log(P)<2.3`$ which, taking into account the differences between the various period determinations and the mean, gives a period ratio from $`1.5`$ to $`2`$.

The same sample of LMC stars analysed by Cioni et al. (cionivar (2001)) was afterwards studied by Lebzelter et al. (agape (2002)) using the AGAPEROS datasets. Ninety per cent of their sources were classified into: regular, SR, irregular and other variability classes. Their analysis covers a time window of up to $`900`$ days and is thus incomplete for the longest periods, but it is more complete for sources with an amplitude of variation below $`0.1`$ mag. The authors found that the period distribution of regular variables has two peaks: one at $`Log(P)=2.0`$ and the other at $`Log(P)=2.5`$. The latter is probably formed by Mira stars. SR variables peak at about $`Log(P)=1.8`$. The peak at the shortest periods is probably the same as the one found here. The peak of mostly Mira stars corresponds to the location of the peak of the C–rich SMC stars. In the Baade’s window, Glass & Schultheis (gs (2002)) found that most of the late–type M giants have $`Log(P)<2.0`$; those few with a secondary longer period have about $`Log(P)>2.4`$. The histogram that can be derived from their Fig. 2 is very similar to that of O–rich stars in the SMC (Fig. 4c).

### 4.6 ($`JK_S`$, $`K_S`$) diagram

In a similar fashion to the LMC sample (Fig. 7 in Cioni et al. cionivar (2001)), variable stars of a different type and variability class are plotted in the colour–magnitude diagrams (CMDs) in Fig. 5. The distinction has been made as follows: Miras have $`A\times 100>110`$ independently of their type, and a regular light–curve with a well defined single periodicity. Similar sources but with a smaller amplitude of variation are SRa stars. SRbs are those stars which indicate the presence of a secondary period. Note that in this last category there is no distinction based on the amplitude of variation. As in the LMC, Mira stars are uniformly distributed in the regions where optically thin O–rich and C–rich, and also more optically thick stars, are located (Fig. 5a). SRa and SRb are almost equally distributed within the two branches of O–rich and C–rich optically thin AGB stars. Fig. 5d shows the distribution of the sources that could not be classified in one of the three categories either because the light–curve has a too low S/N or it is extremely irregular. Most of these stars ($`180`$ out of $`194`$) have $`A>10`$, thus they are variable stars and lie above the TRGB, $`K_S=12.62`$ (Cioni et al. cionitip (2000)). It is possible that stars considerably brighter than $`K_S=10`$ are supergiants and that fainter stars at about $`(JK_S)=0.4`$ are foreground objects. Schultheis et al. (sch (2000)) found variable star candidates at $`(JK_S)<1`$. They conclude that it is probably a result of uncertainties in the determination of extinction. However, they find a displacement of the Bulge stars detected as variables in both $`J`$ and $`K_S`$ to bluer colours than the sequence indicated by all detected ISO stars in the diagram ($`LW2LW3`$, $`LW3`$). It would be interesting to know if some of these stars are supergiants and of which type.

Finally we conclude that all the sources in the sample are variable stars because they have amplitudes above or equal to $`0.1`$ mag either in the RED or BLUE wave band. A similar conclusion was reached by Alard et al. (alard (2001)) for a sample of sources in Baade’s window detected at $`7`$ and $`15`$ $`\mu `$m. The emission at the longest wavelength is a signature of mass–loss rate. Using their dependence between $`M_{BOL}`$ and the mass–loss rate and assuming that the mass–loss rate of highly evolved AGB stars does not depend very strongly on metallicity (van Loon jacco (2000)), we derive that the sources in the present sample lose mass at a rate from $`10^7M_{\mathrm{}}yr^1`$ ($`M_{BOL}=4`$) to $`10^4M_{\mathrm{}}yr^1`$ ($`M_{BOL}=7`$), but most of the objects have $`M_{BOL}=5`$ which corresponds to a mass–loss rate of about $`10^6M_{\mathrm{}}yr^1`$. Note that a complete study of the dependence of mass–loss includes temperature, mass and outflow velocity as well.

### 4.7 ($`K_SLW10`$, $`LW10`$) diagram

The $`7\mu `$m LW2 filter does not include a large dust contribution for non–Mira sources and it is mostly a measure of the photospheric emission. On the other hand, for stars with a relatively thick circumstellar dust envelope, this band is strongly affected by the presence of dust (Alard et al. alard (2001), Omont et al. omont (1999)). The $`12\mu `$m LW10 filter, which is very similar to the IRAS $`12\mu `$m band, includes the $`9.7\mu `$m silicate and the $`11.3\mu `$m SiC dust features. Thus the colour ($`KLW10`$) provides a strong indication of the amount of circumstellar dust for the entire sample.

Fig. 6 shows the colour magnitude diagram ($`K_SLW10`$, $`LW10`$) for the sources in the present sample, distinguishing between Mira and SR variables. Except for a few SR all the sources with $`(K_SLW10)>2`$ are Mira variables (Fig. 6a). From the location of some K– and M–type supergiant (SG) stars in the ISO–MCMS one derives that they have $`(K_SLW10)<1`$ and $`LW10=8`$. They follow a sequence that extends almost vertically to fainter magnitudes, but the nature of those sources with $`LW10>10`$ is not confirmed. Bright O–rich LPVs would have $`(K_SLW10)=3`$ and $`LW10=7`$. These stars would have a late M subtype and not many of them are expected in the SMC. A quantitative statement will be possible once the study of the LMC LPVs detected by ISO is completed. A few SRb stars populate the fainter part (Fig. 6b, c). There are many supergiant stars along the whole sequence that were not classified in any of the three categories (Fig. 6d). It is thus not surprising that their light–curve is more difficult to classify because of a low S/N and irregularities. Thus, except for a few stars overlapping the region occupied by SR and Mira stars we are able to distinguish sources with a poor or non–regular light curve. A more careful look at the near–IR colour–magnitude diagram (Fig. 5) confirms that the unclassified stars with $`(JK_S)=1`$ or a bluer colour are experiencing the supergiant phase because AGB stars have overall $`(JK_S)>1`$. This is quite a good criterion to discriminate between supergiant and AGB stars, though the distinction is better achieved in the combined near– and mid– infrared diagram (Fig. 6).

Most carbon stars are located at about $`(K_SLW10)=1`$ and $`8<LW10<11`$ (Fig. 6–filled points); the sequence of dust obscured C–type stars extends to redder colours but the $`LW10`$ mag follows a shallower slope such as $`(K_SLW10)=4`$ at $`LW10=8`$. On average C–rich stars are redder than O–rich stars.

Planetary nebulae and emission line objects, included in the ISO–MCMS catalogue, also have red colours ($`K_SLW10>2`$) but they have fainter magnitudes with respect to Mira stars and are roughly located in a parallel sequence. Three of these objects are present in Fig. 6d. The object at $`LW10<8.5`$ and $`K_SLW10>4`$ is a known carbon star while the object with the same $`LW10`$ magnitude and a bluer colour is a known variable M giant, because of its red colour we erroneously classified it as a carbon star.

### 4.8 ($`K_SLW2`$, $`LW2LW10`$) diagram

Fig. 7 shows the colour–colour diagram of the variables identified in this study. The O–rich sources with $`(LW2LW10)>0.4`$ are AGB stars or red supergiants. These AGB stars have a significant circumstellar extinction. This is also the colour of obscured O–rich and C–rich stars, and those with the highest amplitude of variability have $`(K_SLW2)>1.5`$. O–rich stars without a considerable mass loss have approximately $`(LW2LW10)=0`$ but can have $`0<(K_SLW2)<1`$. There is not a remarkable distinction between SRa and SRb pulsators. At about $`(K_SLW2)>3.5`$ one can find PN and emission line objects.

It seems (Fig. 7d) that we may have missed the classification of a few Mira variables, O–rich and C–rich SR stars. At this stage it is worth going back to the specific light–curves to find out the reason of this missed classification and to perhaps identify some peculiar objects. Among the objects with $`K_SLW2>1.5`$ (Fig. 7d) J005113.6$``$731035, a known carbon star (RAW 658), has a flat light-curve, while the light–curve of J005136.5$``$732016, identified as a PN, is irregular. Another PN is J005157.9$``$731421 which may have a regular light–curve with an amplitude of about $`0.8`$ mag. The periodogram of the other objects, J004849.4$``$731627 (emission line star), J005212.9$``$730852 (variable M giant) J005304.7$``$730409 and J010020.6$``$730648 (BMB-W 29 carbon star) shows a very low amplitude with flare–like events.

### 4.9 Miscellaneous

SRb, as have been defined in this work, do not constitute a separate class of variables. They are clearly multi–periodic objects and from their photospheric and dust colours they do occupy the same regions of SRa and Mira stars. Until a confirmation of the specific theoretical models (such as those developed by Winters et al. (winter (1994)) to explain the multi–periodicity of the light–curve variation by coupling of stellar pulsation with the dust formation) or follow up radial velocity observations (to infer the binary nature of the sources) are available, the most appropriate class to which assign these stars is that of general SR variables. This is also strengthened by the fact that most of the objects in the sample indicate the presence of more than one periodicity. A sophisticated observational and theoretical analysis of the multiple periods is required to further address this point. Because we did not discriminate on the amplitude of SRb stars some of them do have Mira–like amplitude but are likely to have more than one significant period.

Most of the DENIS sources in the present sample were observed in December 1996 while 2MASS observations date August 1998. This difference in time of about two years may allow us to recognize variable stars by comparing their detection in the near–IR bands of DENIS and 2MASS. Fig. 8 shows the distribution of these differences for all the sources in the sample. We checked that the separate histogram of sources with an assigned MACHO periodicity is the same as that of sources with Flag$`=5`$ (a period was not assigned but the source is probably a variable). The histogram does not account for the systematic difference between the 2MASS and DCMC magnitudes as measured by Delmotte et al. (nausicaa (2002)). Selecting only the variable stars with the largest amplitude produces a similar histogram. The displacement of the peak from zero indicates that the sources are variable stars. Note that the difference is approximately the same in the $`J`$ and $`K_S`$ bands.

## 5 Conclusions

The stars selected from the ISO–MCMS with a confident DENIS/2MASS counterpart have been cross–identified with the MACHO database. The extracted light–curves have been analysed to derive the main period and amplitude of pulsation. The bolometric magnitude has been calculated by integrating the spectral energy distribution from the $`I`$ DENIS to the $`LW10`$ ISO wave band. This, in combination with the correlation of Alard et al. (alard (2001)) between $`M_{BOL}`$ and the mass–loss rate, indicates that most of the stars in the present sample lose $`10^6M_{\mathrm{}}yr^1`$ of matter. A distinction between Mira and SR is obtained by applying an amplitude and regularity selection criterion. The histogram of the amplitude of pulsation of O–rich and C–rich stars has a similar distribution, but on average C–rich stars have a larger amplitude. This may indicate that either most of the C–rich stars in the SMC are of Mira type, or that the metallicity affects the amplitude in such a way that in a lower metallicity environment the amplitude of pulsation is larger. This effect cannot yet be checked for in other metal–poor galaxies in the Local Group, because despite the fact that many AGB stars have been discovered there is not enough information on their variability and type.

The most obscured stars have large ($`K_SLW2`$) and ($`LW2LW10`$) colours. This indicates sources with a thicker circumstellar envelope and allows to easily distinguish AGB stars from supergiant stars and other classes of objects such as PN and emission line sources. The difference between the 2MASS and DCMC magnitudes suggests that most of the sources in the sample are variables.

###### Acknowledgements.

The first author is very grateful to J.-P. Beaulieu for providing the algorithm to detect the periodicities. This paper utilises public domain data obtained by the MACHO Project, jointly funded by the US Department of Energy through the University of California, Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48, by the National Science Foundation through the Center for Particle Astrophysics of the University of California under cooperative agreement AST-8809616, and by the Mount Stromlo and Siding Spring Observatory, part of the Australian National University. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the Mational Science Foundation.

## Appendix A The period–luminosity relation of Mira stars

In order to improve our understanding of the identification of Mira–type variables we decided to re-analyse the LPVs studied by Feast et al. feastpl (1989). These are the Mira stars that define the period–luminosity relation in the LMC. We could find the B1950 coordinates of $`42`$ out of the $`55`$ sources with a “good” or “fair” light–curve as classified by Feast et al. (feastpl (1989)). We could not find in the literature the coordinates of the objects from the sample of Glass & Lloyd Evans (glallo (1981)) and for three objects from the list of Glass & Reid (gr1985 (1985)) or Reid et al. (reidetal1988 (1988)). However, there still remain enough sources for our purpose.

We searched in the MACHO database for the light–curve of the closest LPV to these $`42`$ sources. The identification was easy in most of the cases because of the highly regular and large amplitude light–curve typical of a Mira star. Five sources (W158, GR28, GR4, RGC91, RGC89)333The designation of the sources follows that of Feast et al. feastpl (1989): W=Wood et al. wood1985 (1985), GR=Glass & Reid gr1985 (1985) and we introduce RGC=Reid et al. reidetal1988 (1988). are unfortunately outside the fields observed by MACHO. For one source the MACHO web site gives a persistent initialization error (W94) and for four additional sources (GR13, GR21, RGC37, RGC53) we could not find a nearby counterpart (a first possible candidate is further than $`8`$″). Only one source (W30) has two possible MACHO counterparts and due to the low accuracy of the initial coordinates compared to the MACHO coordinates we decided to remove this source from the discussion. Finally we could successfully extract the light–curves of $`31`$ sources.

We determined the main periodicity, as described in Sect. 3, of the sources with a MACHO–red light–curve. We obtained periods that are at most $`20`$ days different from those derived by Feast et al. (feastpl (1989)) except for $`4`$ sources which result to have a period coincident with one of the extremes of the searching window (GR2, GR3, GR18 and GR29, they have been removed from the following analysis) and source GR17 for which we derive a period of $`51`$ days shorter. The period of source RGC20 falls in the region of aliases (Sec. 3), thus the source has been removed from the following discussion. The quality parameter of the $`26`$ sources with a convincing period determination is always above $`100`$ and the difference between the assigned period and the mean of all the determinations with $`q>15`$ is of the order of $`1`$ day for most of the sources but always below $`10`$ days. These numbers may be different in the case of MACHO–blue light–curves.

Table 5 lists the name of the sources (Column 1), the coordinates at the epoch J2000 precessed from the original coordinates (Columns 2 & 3), the MACHO identification code (Columns 4, 5 & 6), the distance in arcsec between the Feast et al. (feastpl (1989)) source and the MACHO counterpart (Column 7), the period, amplitude and the value of the quality parameter for the red light–curve (Columns 8, 9 & 10).

We conclude that by selecting, from our sample, large amplitude sources with a single periodicity and $`q>100`$ we are identifying Mira variables of the same type as those studied by Feast et al. feastpl (1989). They define a tight PLR and the small scatter in Fig. 3d is probably due to the availability of only one epoch bolometric magnitude. Note that with this selection criterion we also find sources at very bright magnitudes ($`M_{BOL}=7`$) which deviate from the extrapolation of the PLR defined only for sources with $`P<420`$ days. We estimate from Fig. 3b that these stars have $`M>5M_{\mathrm{}}`$ and as previously noted by Feast et al. feastpl (1989) the PLR breaks at high masses. These stars are likely to be “over–luminous” as a result of Hot Bottom Burning (Zijlstra et al. albi (1996), Whitelock et al. whiti (2003)). However, because of their very high luminosity they could also be supergiants.

## Appendix B Bolometric correction(s)

### B.1 BCI versus $`IJ`$

Fig. 9 compares the bolometric correction obtained from the integration under the SED (Sect. 3.1) and the relation derived by Alvarez et al. (alva (2000)) for cool O–rich stars. The agreement is very good for the O–rich stars of our sample. This confirms that the relation is valid for the SMC as well as for the LMC and the Bulge. However this relation overestimates the correction for C–rich variables in the SMC. On the other hand we should check that the integration under the SED does not underestimates the bolometric correction (if the stars have a circumstellar envelope and dust extinction the correction for energy re–radiated by dust becomes important). Our measurements have been corrected for interstellar extinction while the relation has not. The agreement for O–rich stars might be a consequence of the systematic difference between the DENIS filters, used in this work, and the broad pass–bands by Bessels & Brett (besbre (1988)) and Bessel (besse (1990)) used to derive the relation. The disagreement for C–rich stars might be a consequence of the SiC dust feature which causes the spline not to follow the continuum of the SED. However the data of Alvarez et al. (alva (2000)) include only $`6`$ C–rich stars and a total of four stars belonging to the LMC or the SMC. In addition the spectral range covered by their measurements is limited to $`2.5\mu `$m. The zero–point of the bolometric scale is the same in both studies.

### B.2 BCK versus $`JK`$

Fig. 10 compares the bolometric correction obtained as in Sect. 3.1 with the relation derived by Montegriffo et al. (monte (1998)) for Population II giants of low (\[Fe/H\]$`<1.0`$) and ”high” metallicity (\[Fe/H\]$`>1.0`$). They used $`M_{BOL,\mathrm{}}=4.75`$. Their relation was defined up to $`(JK)=0.8`$ but being the metal–poor relation very close to the metal–rich relation at redder colours we consider the latter as a possible extrapolation. This results in a fairly good fit for the O–rich variables while C–rich variables lie systematically below the relation. This is a similar effect as that observed for $`BC_I`$ versus $`IJ`$ (Fig. 9). The giants in the sample of Montegriffo et al. (monte (1998)) are members of galactic globular clusters. C–rich stars are rare members of these systems. The spectral range covered by the observations is limited to $`3.5\mu `$m ($`L`$ band) and they derive the bolometric magnitude by integrating the SED which is previously re–constructed by means of a set of Plank functions at a given wavelength. These effects may influence the correction to be applied to the flux measured only in the $`K`$ band. Besides the $`2.2\mu `$m filters used in both studies are not strictly the same. We conclude that we cannot at present compare the correction for C–rich stars while O–rich stars are qualitatively in agreement. It would be useful to compare stars of other galaxies in the Local Group at a different metallicity.