Free Access
Issue
A&A
Volume 571, November 2014
Article Number A43
Number of page(s) 16
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/201423802
Published online 06 November 2014

© ESO, 2014

1. Introduction

Stellar clusters are the best approximation in nature of a simple stellar population, they are the clocks of the Universe, and the building blocks of galaxies. More than 90% of massive stars are found in young and massive stellar clusters (de Wit et al. 2005).

In the inner regions of the Galactic disk, the identification of stellar clusters is a difficult task because of their irregular shapes and density profiles, and because of high and patchy interstellar extinction. About 3000 new candidate stellar clusters have been identified as significant peaks of stellar counts in the 2MASS and GLIMPSE catalogs, or with visual inspection of images (e.g., Ivanov et al. 2010; Mercer et al. 2005; Dutra et al. 2003; Froebrich et al. 2007; Borissova et al. 2011, 2014). However, only a few of these candidates have been confirmed with follow-up spectroscopic and photometric studies, mostly on the near-side of the Galactic plane (e.g., Messineo et al. 2009; Davies et al. 2012; Chené et al. 2013). A multiwavelength photometric screening/selection of candidate massive clusters is mandatory for the best use of the observing facilities, since spectroscopic follow-up can only be accessible for a limited number of candidate clusters. Typically, 50% of the detected overdensities are found to be spurious (Froebrich et al. 2007). For example, overdensities may result from regions of low interstellar extinction (e.g., Dutra et al. 2002). An overdensity can be defined as a cluster of stars if it is made of stars at the same distance and approximately the same age. Depending on the dynamical status, a cluster can be a bound cluster or an association. Typically, cluster members have similar interstellar extinction, and may be recognized on color–magnitude diagrams (CMDs) as specific sequences.

A new class of stellar clusters has recently been discovered, the red supergiant clusters (RSGCs), which are characterized by a large number of red supergiants (RSGs; Figer et al. 2006; Davies et al. 2007; Clark et al. 2009; Negueruela et al. 2010, 2011; González-Fernández & Negueruela 2012). RSGC1, RSGC2, RSGC3, RSGC4, and RSGC5 contain 14, 26, ~35, >13, and 7 RSGs, respectively, or collectively ~15% of all known RSGs in the Galaxy (Messineo et al. 2012). These newly discovered RSGCs are all concentrated between longitude l = 25° and l = 30°, i.e., close to the near end of the Galactic bar, where the bar appears to meet the Scutum-Crux spiral arm. However, their census is incomplete (e.g., Messineo et al. 2012); further searches for RSGCs are needed to investigate Galactic structure. RSGCs are dominated by RSGs, which are intrinsically bright at infrared wavelengths. Their near-infrared CMDs are characterized by gaps of several magnitudes between the RSGs and the blue supergiant members. RSGCs are detected as overdensities of solely infrared bright stars, so bright (typically Ks≲ 7.0 mag at 6 kpc) that they can be detected throughout the Galactic plane; infrared dimmer main sequence members are hard to identify in the glare of bright RSGs.

In this paper, we analyze five candidate clusters rich in infrared bright stars that show different types of CMDs by means of a quantitative analysis of near-infrared spectra of their brightest stars. In Sect. 2, we describe the observed candidate stellar clusters; in Sect. 3, we report on the available infrared spectroscopic and photometric data. Spectral and photometric classifications are given in Sect. 4. The cluster properties are discussed in Sect. 6. Finally, our findings are summarized in Sect. 7.

Table 1

List of candidate cluster positions, apparent radii, AKs, and modulus distances (DM).

2. Targeted candidate clusters

thumbnail Fig. 1

GLIMPSE 3.6 μm images of the candidate clusters listed in Table 1: cl1.5, cl9.5, cl16.5, cl49.3, and cl59.8. Spectroscopic targets are labeled as in Tables 3 and 4. Diamonds indicate red giants, squares Mira-like AGB stars, crosses RSG stars, and triangles early- and yellow-type stars. The coordinates of the image centers are in degrees.

A random sample of candidate stellar clusters were selected among overdensities of infrared-bright stars in both GLIMPSE and 2MASS images (see, e.g., Ivanov et al. 2010, 2002). The candidate clusters were selected to have different types of CMDs, with gaps and without gaps (see Sect. 6). The observed candidate clusters are listed in Table 1. The list reports their coordinates, which are the flux weighted centroids of 5× 5Ks band images from 2MASS. The overdensities of infrared-bright stars are shown in Fig. 1, and in the star counts of Table 2.

Table 2

Counts of stars with Ks> 10.0, 8.5, and 7.0 mag in the fields listed in Table 1, and, as a comparison, in control fields of equal area.

3. Data

3.1. Near-IR spectroscopy

3.1.1. UIST spectra

A set of spectroscopic data was taken with the UKIRT 1–5 micron Imager Spectrometer (UIST; Ramsay Howat et al. 1998) on Mauna Kea under program ID H243NS (PI: Kudritzki) on 2008 July 24. We used the long-K grism, which covers from 2.204 μm to 2.513 μm at a resolving power (R) of 1900. Integration times varied from 10 s to 45 s per exposure, and the number of exposures varied from 8 to 20.

Pairs of adjacent frames at different nod-positions were subtracted from each other and flat-fielded, wavelength-calibrated with Ar arc lines, and corrected for atmospheric absorption and instrumental response. Curves of atmospheric absorption and instrumental response were generated by dividing the observed spectra of standard stars (with spectral types from B2 to B9) by blackbody curves. Linear interpolation was used to remove Brγ lines and possible He i lines from the spectra of the standards. A total of 22 stars were observed with UKIRT, and are listed in Tables 3 and 4; their spectra are shown in Fig. 2.

thumbnail Fig. 2

Long-K spectra taken with UIST. Identification numbers refer to Table 4. All are late-type stars, with the exception of star 13.

3.1.2. SOFI spectra

For four of the targeted candidate clusters, additional low-resolution and medium-resolution spectra were obtained with the SofI spectrograph mounted on the NTT telescope (Moorwood et al. 1998). Data were taken under program 60.A-9700(E) at Paranal-La Silla Observatory on 2010 August 3, and under program 089.D-0876 on 2012 June 1. The pixel scale is 0.̋288 pix-1. The low-resolution Red grism in combination with the \hbox{$0\farcs6 \times 290\arcsec$} slit yielded R ≈ 980 over the range ~1.50μm to ~2.49μm. A few spectra were taken with the HR grism, the Ks filter, and a \hbox{$1\farcs0$}-wide slit (R ≈ 1900).

Typically, two slit positions per cluster were observed. The rotator angle was set to optimize the number of stars falling onto the slit. For each slit, frames were taken in a nodding sequence ABBA, which included a small jittering between positions, with detector integration times (DITs) from 1.18 s to 120 s. B-type stars were selected as standards, and observed with the same settings as the targets.

Data reduction was carried out with IDL programs, and the Image Reduction and Analysis Facility (IRAF), using the NOAO/onedspec package1. Pairs of subsequent images at different nod-positions were subtracted and flat-fielded. Wavelength calibration was performed with observations of arc lamps. Only spectra with a signal-to-noise ratio above 30 were used. For each star, typically, four spectra were combined and then corrected for atmospheric transmission and instrumental response. Hydrogen and helium lines were removed from the spectra of the standard stars with a linear interpolation; the intrinsic continuum slopes of the standard spectra were eliminated by dividing them by blackbody curves (calculated at the Teff of the standard stars). A total of 28 low-resolution spectra and 13 medium-resolution spectra were obtained (see Table 4, Figs. 3, and 4).

thumbnail Fig. 3

Low-resolution H and K spectra taken with SofI. Identification numbers are taken from Table 4. Late-type stars are shown in panel a). Early-type stars are shown in panel b) together with a few low-resolution template spectra from the library of Lancon & Rocca-Volmerange (1992), and medium-resolution template spectra from the IRTF library (Rayner et al. 2009).

3.2. Photometric data

For the target stars, simultaneous photometric measurements in bands J, H, and Ks were available in the Two Micron all Sky Survey (2MASS) catalog (Skrutskie et al. 2006), with 1.̋0 resolution; for two fields, additional near-infrared measurements (simultaneous I, J, and Ks) were available from the Deep Near-Infrared Survey (DENIS) catalog (Epchtein et al. 1999). Deep J, H, and K photometry was obtained from the UKIDSS Galactic Plane Survey (Lucas et al. 2008; Hodgkin et al. 2009).

Mid-infrared data were available from the Galactic Legacy Infrared Mid–Plane Survey Extraordinaire (GLIMPSE) with the Spitzer Space Telescope, from the Midcourse Space Experiment (MSX) Spatial Infrared Imaging Telescope (SPIRIT III), and from the Wide-field Infrared Survey Explorer (WISE) satellite (Egan et al. 2003; Benjamin et al. 2003; Fazio et al. 2004; Wright et al. 2010). The Spitzer/IRAC camera acquired simultaneous images in the four channels at 3.6, 4.5, 5.8, and 8.0 μm, with a spatial resolution of 1.′′2 and a sensitivity of 2 mJy. MSX covered from 8 μm to 21 μm, with a pixel scale of 18.̋3, and a sensitivity of 0.1 Jy at the short-wavelength. WISE bands are centered at 3.4, 4.6, 12, and 22 μm; detectability limits are 0.08, 0.11, 1, and 6 mJy; spatial resolutions are 6.′′1, 6.′′4, 6.′′5, and 12.′′0, respectively.

DENIS and 2MASS data were cross-matched using a search radius of 2.′′0; 2MASS astrometry was retained. Mid-infrared counterparts from the MSX catalog were searched within a radius of 5′′; matches from the WISE and GLIMPSE catalogs within a radius of 2′′. In addition, we checked for possible visual counterparts in the NOMAD catalog of Zacharias et al. (2004), using a search radius of 2′′.

Photometric measurements of the spectroscopic targets are listed in Table 3.

3.2.1. UKIDSS photometry

We used JHK photometry from the UKIDSS Galactic Plane Survey (GPS). For every position, JHK images were obtained from the UKIDSS archive. Raw frames were processed by the standard UKIDSS pipeline (Lucas et al. 2008); corrections for linearity, dark, flat-fielding, decurtaining (the removal of a pseudo-periodic ripple), defringing (the removal of interference fringes due to atmospheric emission lines), sky-subtraction, and cross-talk were applied.

The original WFCAM pixel scale is of 0.̋4 pix-1; however, GPS observations were performed with a 2 × 2 micro-stepping technique, which improved the sampling; a microstep of 11.5 pixels was used. Each observation comprises eight individual exposures of duration 10, 10, and 5 s to make up integration times on sources of 80, 80, and 40 s in the J, H, and K filters, respectively. We resampled the exposures of each observation in a finer grid (3 × 3) in order to reduce the image noise, and to improve the detectability of faint sources; we used a bilinear interpolation. Exposures were combined with a three σ clipping to eliminate possible cosmic rays. The astrometric distortion of the final mosaic was corrected by using a set of unsaturated stars with 2MASS photometry, and performing a polynomial spatial de-warping of third order. Source extraction was performed using the DAOPHOT PSF-fitting algorithm by Stetson (1987). For each mosaic, a set of bright and isolated unsaturated stars were selected, and a invariable PSF was modeled. The detection threshold was set to 4σ (the standard deviation). Finally, individual J, H, and K catalogs were photometrically calibrated with overlapping 2MASS sources, positionally cross-correlated, and combined.

For saturated stars that could be uniquely associated with a 2MASS point source, 2MASS magnitudes were retained.

3.3. Information available on SIMBAD

The SIMBAD database reports information only for one target. We classified star 13 as a B0-F0 star; the star coincides with HD166307, which has been classified as an A0III star at optical wavelength (Houk & Smith-Moore 1988).

thumbnail Fig. 4

Medium-resolution K spectra taken with SofI. Identification numbers are taken from Table 4. All stars are late-type stars, with the exception of star 39.

4. Stellar classification

We detected candidate RSGs in two fields. In order to confirm RSGs in the direction of the inner disk of the Milky Way, a complex procedure is required. Spectroscopic observations alone cannot firmly distinguish between RSGs and red giant stars. When analyzing RSGs, it is crucial to combine spectroscopic information (see Sect. 4.1) with photometric quantities (see Sects. 4.1.3 and 4.1.4), and to estimate their fundamental parameters (effective temperatures, Teff , and luminosities). Distances of late-type stars are usually estimated in two ways, either approximated with kinematic distances, or estimated using interstellar extinction as an indicator of distance (e.g., Habing et al. 2006; Messineo et al. 2005; Drimmel et al. 2003). In the inner Galaxy kinematic distances alone cannot be trusted as they rely on the assumption of circular orbits that is not valid in the inner 3.5 kpc (e.g., de Vaucouleurs 1964); our estimates of distances are based on interstellar extinction (Sect. 4.2). Luminosities are described in Sect. 4.3.

4.1. Spectroscopic classification

4.1.1. Late-type stars

We defined as late-type stars those stars with effective temperatures lower than 4500 K (red giant, asymptotic giant branch stars (AGBs), and RSGs). We define a star as a RSG when it has a Teff lower than 4500 K and a luminosity of L/L> 104, which corresponds to an initial mass of 9 M (Ekström et al. 2012). Late-type stars can easily be identified with K band spectroscopy, because of their strong CO band-head at 2.29 μm (e.g., Kleinmann & Hall 1986; Wallace & Hinkle 1996; Rayner et al. 2009; Ivanov et al. 2004). Near-infrared spectra of late-type stars (mostly Mira AGBs) may show continuum absorption by water. Water absorption affects both ends of the H band spectrum, from ~1.4 μm to ~1.55 μm, and from ~1.75 μm to ~1.80 μm, and the blue side of the K band from 2.0 μm to 2.1 μm (e.g., Comerón et al. 2004; Blum et al. 2003; Frogel & Whitford 1987; Alvarez et al. 2000; Rayner et al. 2009). Water absorption in H band is most easily detected, as a curved (pseudo)-continuum. Highly variable water absorption is found in Mira AGBs because of their pulsation (Matsuura et al. 2002). In Fig. 3, the stellar continua of stars 1 and 24 clearly show absorption by water.

For red giants and RSGs, the equivalent width of the CO band-head, EW(CO), can be used as a temperature indicator because it increases linearly with decreasing Teff. Red giants and RSGs follow two different relations (Blum et al. 2003; Frogel & Whitford 1987; Figer et al. 2006); for a given temperature, RSGs have stronger CO bands than red giant stars. EW(CO)s of Miras are variables, do not correlate with Teff values, and their EW(CO)s may be as large as those of late RSGs (Blum et al. 2003). Water absorption decreases with increasing stellar luminosity (Blum et al. 2003; Frogel & Whitford 1987). All these considerations mean that combined information on water absorption and EW(CO) is useful for estimating luminosity classes.

We measured the EW(CO)s using the feature and continuum regions that are specified in Table 5; we obtained spectral types for the targets by comparing their EW(CO)s with those of reference spectra from the atlas of Kleinmann & Hall (1986); the reference spectra were smoothed to match the spectral resolution of the targets (e.g., Figer et al. 2006). Twelve targets were observed with both UIST and SOFI detectors; by comparing their resulting EW(CO)s, the typical accuracy of spectral types is found within two subclasses; star 24 is a Mira-like star that went from >M7 (UIST run) to M0 type (SOFI run). For spectra taken with the medium-resolution mode of SofI, we selected a narrower bandpass for the CO feature (see Table 5). A scaling factor of 1.4 was measured between the EW(CO)s from the medium-resolution SofI and those from UIST. For star 27, which was observed with both the low and medium modes of SofI, we obtained a RSG-type of K4 and K2, respectively. The degeneracy between RSGs and red giants disappears for EW(CO) values larger than ~ 43 ± 4 Å. Miras stars, with their erratic behaviors, may also have EW(CO)s larger than 43 Å. Stars 1, 2, 15, 24, and 28 in Table 4 have EW(CO) values larger than 48; these stars are candidate AGBs or RSGs.

Table 5

Spectroscopic indexes. Spectral regions used as features and continuum are specified.

Water indexes were obtained from the de-reddened low-resolution SofI spectra, which cover H and K bands. We used a water index with continuum and absorption regions defined as in Blum et al. (2003). We measured the water index (100 ∗ (1− ⟨ FH2O/Fcontinuum ⟩)) with a quadratic fit to two continuum regions (see Table 5); FH2O/Fcontinuum is the average ratio of observed flux densities in the water region and estimated flux densities (with a fit). For comparison, we estimated the EW(CO)s and water indexes of a sample of known AGBs and RSGs with spectra available from the IRTF library (Rayner et al. 2009). In Fig. 5, we plot water indexes versus EW(CO)s of the targets, as well as of reference spectra. A combination of water indexes and EW(CO)s allows us to distinguish between RSGs and Mira AGBs; we find that known Mira AGBs have water indexes greater than 15 Å, as in Blum et al. (2003). This confirms that targets 1, 2, and 24 are Mira-like stars. Water indexes for known RSGs and semi-regular AGBs (SR, e.g., Alard et al. 2001) are typically negligible; SR AGBs have EW(CO)s as large as ~44 Å; RSGs up to ~58 Å (see Fig. 5 and Table 4). SR AGBs contaminate the sample of spectroscopic candidate RSGs (EW(CO) ≳ 43 Å); luminosities are crucial for classifying RSGs.

4.1.2. Early-type stars

Low-resolution infrared spectra of O-, B-, A-, and F-type stars are characterized by hydrogen (H) lines (e.g., Hanson et al. 1996, 1998; Lancon & Rocca-Volmerange 1992; Rayner et al. 2009). The low-resolution SofI spectra of stars 9, 11, 13, 19, 20, 21, and 26 show H i lines in absorption at 1.555 μm, 1.570 μm, 1.588 μm, 1.611 μm, 1.641 μm, 1.681 μm, 1.737 μm, and 2.166 μm (see Fig. 3). The strengths of the hydrogen lines indicate types from B0 to F0 for stars 11, 13, 20, and 26; the presence of weak CO bands imply types from F5 to G5 for stars 9, 19, and 21.

A spectrum of star 39 was obtained with the medium-resolution mode SofI in K band. A Brγ in absorption and a He i at 2.1127 μm are detected. The presence of the He i line implies a B0-8 I, or a O9.5-B3V.

The coordinates and spectral types of early-type stars are listed in Tables 3 and 4.

thumbnail Fig. 5

H2O indexes versus the EW(CO)s of targeted late-type stars and comparison samples of AGBs (Mira and semi-regular variables) and RSGs from the IRTF library (Rayner et al. 2009).

4.1.3. Photometric variability index

Mira AGBs are characterized by periodic photometric variations that can reach up to ~1 mag in Ks (Messineo et al. 2004). In contrast, only a small fraction of RSGs (about 20%) are known to vary, and typically RSGs have amplitudes of a few tenths of a magnitude in K band (Yang & Jiang 2011).

Since the two DENIS and 2MASS J filters are similar, we used both measurements for identifying candidate variable stars, i.e., stars with J(DENIS) > 7.5 mag, and | J (DENIS) J(2MASS) | > 3 × 0.15 mag, or with Ks (DENIS) > 6.0 mag, and | Ks (DENIS) Ks (2MASS) | > 3 × 0.15 mag. For nonvariable stars, | Ks (DENIS) Ks (2MASS)| and | J (DENIS) J(2MASS)| are within 0.15 mag (Schultheis et al. 2000; Messineo et al. 2004). We also flagged as variables those targets with indication of variability from the shortest WISE band. The late-type stars 14, 17, 24, and 34 are candidate variables, as well as the early-types 13 and 20, and the F-G stars 19 and 21. Star 24 is a Mira-like AGB, as inferred from the water index; the other variable late-type stars could be SR AGBs (Alard et al. 2001). None of the RSGs found has the variability flag on.

thumbnail Fig. 6

Top panels: interstellar extinction AK versus distance moduli of red clump stars in the fields of cl16.7 (left) and cl49.3 (right). As a comparison, the model of Drimmel et al. (2003) is shown with a dashed line. Bottom panels: UKIDSS JK versus K diagram of data-points within an area of 10× 10 from the estimated centers of cl16.7 (left) and cl49.3 (right).

4.1.4. Target intrinsic colors and interstellar extinction

For red giants and RSGs, we assumed the intrinsic VK, JK, and HK colors per spectral-type provided by Lejeune & Schaerer (2001) and Koornneef (1983). For each spectral type from K0 to M5, the VK (or JK) colors of red giants and RSGs agree within 0.28 mag (or 0.09 mag), respectively (Koornneef 1983); this quantity is comparable to the color change between two neighboring subspectral types of the same luminosity class. We obtained intrinsic IJ colors by interpolating a colored-isochrone of 6 Gyr and solar metallicity (Pietrinferni et al. 2004). The obtained IJ colors of late-type giants agree with those of Lejeune & Schaerer (2001) for RSGs (430 Myr old) within 0.05 mag. The Koornneef photometric system agrees with the 2MASS system within 0.09 mag (Carpenter 2001).

For early-type stars, infrared-colors were taken from Bibby et al. (2008), Humphreys & McElroy (1984), Wegner (1994), Koornneef (1983), and Wainscoat et al. (1992).

For every target, we estimated AKs with the assumed intrinsic colors, and a power-law extinction curve with an index of −1.9 (Messineo et al. 2005). This law yields an excellent agreement between estimates of AKs from the observed (IJ), (JKs), and (HKs) colors. For eight stars with DENIS and 2MASS measurements, we obtained AKs(from JKs) − AKs(from IJ) = − 0.04 mag with σ = 0.14 mag; AKs(from JKs) − AKs(from HKs) = −0.03 mag with σ = 0.09 mag.

Typically, we adopted the AKs values from the (JKs) colors. For every field, an average extinction (AKs) was measured, as the average of the AKs values of the observed late-type stars (with the exclusion of AGB stars). The field results are listed in Table 1. Stars with a value of AKs that differs from the average extinction by more than three standard deviations were classified as foreground (background) sources, as indicated in Table 6. We used the shape of the continuum (quantified by the water index) to distinguish between objects with circumstellar envelopes (e.g., Mira AGB) and obscured, but naked, luminous objects.

4.2. Interstellar extinction as a distance indicator

Distances were derived by matching the target AKs value with a known curve of AKs values versus distances (along the target line of sight). The targets are distributed in five fields along the Galactic plane, at galactic longitudes of 15, 95, 167, 49.°3, and 59.°8, respectively, and approximately zero latitude.

A relation between AKs and distance along a given line of sight can be derived with primary distance calibrators. Thanks to the availability of deep near-infrared surveys, such as VVV (Soto et al. 2013) and UKIDSS (Lucas et al. 2008), red clump stars (analogous to the Galactic horizontal branch stars in metal rich globular clusters) have been detected throughout the Galaxy, and may serve for this purpose (e.g., Gonzalez et al. 2011; Drimmel et al. 2003). We created UKIDSS K versus JK diagrams of datapoints within 10× 10 from the cluster centers, and we visually selected by eyes the locus of a well visible red clump sequence. Per bin of JK color, we analyzed the distribution of K magnitudes (around the clump region), and estimated the peak of the red clump with a Gaussian fit. We assumed an intrinsic (JK)o of 0.68 mag (Babusiaux & Gilmore 2005; Gonzalez et al. 2011) and an absolute magnitude in K band, MKs , of −1.61 mag (Alves 2000). Quoted values in the literature range from MKs = −1.54, −1.55 mag (Groenewegen 2008; Pietrinferni et al. 2004) to MKs = −1.72 mag (Babusiaux & Gilmore 2005). For all fields but one, distances were derived with clump stars in the AKs region of interest (average AKs of the targeted stars); for field cl59.8 it was not possible because there were not enough clump stars. For two fields, cl16.7 and cl49.3, it was possible to derive a curve of variation of AKs with distance (e.g., Drimmel et al. 2003), as shown in Fig. 6.

Alternatively, a curve of AKs versus distance along a given line of sight was derived with a model of Galactic dust distribution from COBE/DIRBE (Drimmel et al. 2003; Drimmel & Spergel 2001). A scaling factor of 0.74 was applied to transform the AKs of Drimmel et al. (2003) and Rieke & Lebofsky (1985) into the AKs derived with a power law of index = −1.9 (Messineo et al. 2005).

Estimated distances are listed in Table 1. In the following, we assume the distances derived with red clump stars. The Drimmel model tends to give shorter distances than those from the red clump fitting (Fig. 6); however, any detected RSGs would still be a RSG at the smaller distance provided by the Drimmel model.

In field cl49.3 it was possible to derive a spectrophotometric distance for star 39, an early B star with Ks = 11.19 mag. Typically, early-B stars have an absolute K of −6.45 ± 0.42 mag, −3.52 ± 1.00 mag, or −2.38 ± 0.77 mag, when they are supergiants (Ia), giants, or dwarfs, respectively (Humphreys & McElroy 1984; Wegner 1994; Bibby et al. 2008; Koornneef 1983). A supergiant class and a distance of 26 kpc can be excluded. By assuming a dwarf or a giant class, we derived a distance modulus (DM) of 13.02 ± 0.77 mag (4.0 kpc) and 14.16 ± 1.00 mag (6.8 kpc), respectively.

4.3. Estimated luminosities

Bolometric magnitudes, Mbol , and luminosities were calculated with Ks magnitudes, AKs, estimated DMs, and bolometric corrections, BCK (see Table 6). Values of BCK per spectral type were taken from Levesque et al. (2005) and Messineo et al. (2011); the solar Mbol was assumed to be +4.74 (Bessell et al. 1998).

For late-type stars, the spectral-energy distribution (Teff from 4500 to 2000 K) peaks in the near-infrared (0.61.5 μm), or longward of this range when circumstellar dust is present (e.g., Ortiz et al. 2002; Figer et al. 2006). For late-types, we additionally integrated available flux densities over frequency ν using a linear interpolation. At short wavelengths, flux densities were extrapolated to zero with a blackbody function, at long wavelengths with a linear interpolation through the last two points. Direct integration yields more robust luminosities for dusty AGB stars, for which BCK of naked stars are not applicable (Messineo 2004; Ortiz et al. 2002). The solar constant, Cbol, was assumed to be −18.90. Our integration method provided Mbol values in agreement within errors with those derived with the BCK values.

Known RSGs have typical Mbol from −5.5 mag to −9.0 mag (Figer et al. 2006; Davies et al. 2007; Clark et al. 2009). Galactic mass-losing AGB stars rarely exceed Mbol≈ −7.1 mag (Vassiliadis & Wood 1993); this theoretical limit may be exceeded with the onset of the hot bottom burning (HBB) process that breaks down the classic relation between core mass and luminosity (see, e.g., Groenewegen et al. 2009; Paczyński 1971). However, the efficiency of the HBB mechanism anti-correlates with metallicity, and it is usually associated with a stellar regime of high mass-loss and superwind (e.g., Wood et al. 1983; Marigo 2001; Vassiliadis & Wood 1993); typically, AGBs have strong water absorption in the bright end (Mbol from −5.5 to −7 mag, Habing & Olofsson 2003; Alard et al. 2001). In the Milky Way, Mira AGBs have mostly late-types (M7-M9, as reported, for example, by Rayner et al. 2009); the distribution of spectral types of known RSGs peaks at M2-M4 (Davies et al. 2007).

For selecting clusters of RSGs, we calculated the photometric properties of the targets at the distance provided by their average AKs per field (as explained in Sect. 4.2); after having selected candidate RSGs, the procedure was reiterated with the AKs of the candidate RSGs. The results are summarized in Table 6. In field cl16.7 (DM = 12.9 mag), stars 22 (K4.5-5I) and 23 (K5.5I) are two RSGs with Mbol = −6.18 mag and −5.71 mag. In field cl49.3 (DM = 14.2 mag), star 27 (K4I) with a Mbol = −7.82 mag is a RSG star, along with the nearby stars 28 (M1-1.5I), 29 (K5.5-M0I), and 30 (K5-K5.5I), which have Mbol = −6.57 mag, −6.07 mag, and −5.76 mag, respectively. In none of these stars did we detect signs of variability or water absorption. There is only one Mira-like star exceeding Mbol = −5.5 mag in field cl1.5, which is star 1. For all other late-type stars, we estimated Mbol typical of giant stars. There are a number of giants (after excluding the detected Mira-like AGBs) with estimated types later than M7III (3, 5, 12, 14, 15, 42, and 45); their Mbol are below the values expected for faint RSGs.

5. Lines of sight towards cl16.7 and cl49.3

In the CMD of field cl16.7, at 1675 of longitude, and −0.°63 of latitude, the red clump trace is clearly visible from JK ≈ 0.8 mag and K ≈ 10 mag to JK ≈ 3.2 mag and K ≈ 13.5 mag (see Fig. 6). Clump stars appear evenly distributed along the line of sight, and their number increases with distance. For the RSGs in cl16.7, we derived an average AKs of 0.49 mag; from this AKs, we obtained a DM of 12.93 mag, which places them on the Scutum-Crux arm (see Fig. 7, as well as the distribution of Galactic massive clusters in Messineo et al. 2009, 2014).

Field cl49.3 is located at a longitude of 4934 and a latitude of +072. The DM inferred from the RSGs places the cluster cl49.3 onto the Sagittarius-Carina arm (see Fig. 7). In the CMD of field cl49.3, red clump stars are traceable from JK ≈ 1 mag and K ≈ 11.5 mag to JK ≈ 2.1 mag and K ≈ 14.5 mag (see Fig. 6). There is a sudden increase in the number of red clump stars at JK ≈ 1.4 mag, and a heliocentric distance larger than ~5 kpc. This overdensity of red clump stars occurs at the shortest Galacto-centric distances along the line of sight, where the stellar density is higher; they are also encountered when tangentially approaching the Sagittarius-Carina arm (heliocentric distance 59 kpc, see Fig. 7). Similar structures appear in the CMDs of neighboring bins of longitudes. Some models predict the Galactic outer Linbland resonance to occur inside the solar circle, at a Galacto-centric distance of about 5 kpc (e.g., Habing et al. 2006); such resonances may generate/enforce extended structures (e.g., rings and pseudo-rings, Buta 1995). The Galactic long bar only extends to about 30° longitude (Cabrera-Lavers et al. 2008).

thumbnail Fig. 7

XY view of the Milky Way. The Galactic center is at (0, 0) and the Sun is at (8, 0). The dashed lines indicate the line of sights to cl16.7 and cl49.3, which are marked with squares. A sketch of the spiral structure is also drawn and the arms are labeled (Cordes & Lazio 2002). The X-axis is oriented along the Sun-Galactic center line of sight.

6. Clusterings

Because of their young age most RSGs are mostly found in stellar clusters. Massive young clusters rich in RSGs have been identified via their near-infrared CMDs, e.g., the RSGC1 (Figer et al. 2006), RSGC2 (Davies et al. 2007), and RSGC3 (Clark et al. 2009). In their (JKs) versus Ks diagrams there are gaps of several magnitudes in Ks between RSGs and blue supergiants (Figer et al. 2006); the gap inversely correlates with the cluster age (Figer et al. 2006; Messineo et al. 2009). For comparison, Mira AGBs are typically isolated on a scale of 1 and are among the brightest stars, but are typically the reddest (Messineo et al. 2005, 2004); in the CMDs of stars surrounding Mira stars gaps may only be present stochastically.

thumbnail Fig. 8

2MASS/UKIDSS JKS versus KS CMDs: cl1.5 (top left), cl9.5 (top right), cl16.3 (middle left), cl49.3 (middle right), and cl59.8 (bottom left). Photometry was obtained with a PSF-fitting technique (see text). For each candidate cluster, a diagram of the candidate cluster is shown in the left panel, and that of a control field of equal area in the right panel. Spectroscopically detected early and yellow stars are marked with filled triangles, Miras with squares, RSGs with crosses, and red giants with diamonds; identification numbers are taken from Table 3. In the CMDs of fields cl16.3 and cl49.3, which are rich in RSGs, a continuous curve displays an isochrone of 20 Myr from Lejeune & Schaerer (2001), reddened with AKs = 0.50 mag and 0.48 mag, and shifted to DM = 12.93 mag and 14.22 mag, respectively; a dashed line shows the sequence of red clump stars (see Sect. 4.2). The CMD of cl59.8 shows three isochrones from Marigo et al. (2008) with solar metallicity and ages of 100 (dashed line), 200 (continuous line), and 300 Myr (dotted line), which were reddened (AKs = 0.55 mag) and shifted to DM = 12.9 mag.

The CMDs of the analyzed candidate clusters are shown in Fig. 8. The CMDs of fields cl16.7 and cl49.3 are characterized by gaps; in these two fields we detected several stars with luminosities typical of RSGs. The infrared CMDs nicely complement and strengthen the detection of six new RSGs, as illustrated in the following.

6.1. RSGs in field cl16.7

The 2MASS/UKIDSS (JKs) versus Ks diagram of field cl16.7 shows two bright RSGs (stars 22 and 23) at (JKs) ≈ 1.9 mag and Ks≈ 4.8 mag (Fig. 8). The two RSGs (stars 22 and 23) have similar extinction values (AKs = 0.46 and 0.52 mag), and spectral types, and so they are drawn from the same population. The CMD shows a group of bright stars that extends over a main sequence to about 9.5 mag. There is a gap of about 4.5 mag between these stars and the two RSGs. The two other bright stars with Ks = 5.91 and 6.64 mag are field giant stars.

Star 26 (early or yellow-type) has AKs~ 0.58 mag, which is consistent with the extinction of the two RSGs. On the CMD, the star is located at (J Ks) = 1.0 mag and Ks = 12.1 mag, over a main sequence (from Ks = 18 mag to 14 mag). It is probably not related to the young RSGs; for DM = 12.9 mag, we obtain MKs = −1.4 mag (typical of old late-B giants, see, e.g., Clampitt & Burstein 1997).

Star 24 is 5′′ from the IRAS 182071441 source, within its positional error ellipse. Its spectrum has strong water absorption. A few rare RSGs have shown water absorption in H band, e.g., My Cep (M7I) in NGC 7419 (Rayner et al. 2009). My Cep, the brightest of the five RSGs in NGC7419, is a pulsating star, and shows OH and water masers (Benson et al. 1990). Star 24 can be a late giant (M8III) or a faint M3I with an integrated Mbol of −4.97 mag. Since the two RSGs (K5I, K5.5I) are bluer and brighter (Mbol = −6.2 and −5.7 mag) than star 24, we concluded that this was most likely an unrelated AGB star. Accurate radial velocities are needed to confirm the proposed scenario.

A colored theoretical isochrone from Lejeune & Schaerer (2001), with increased mass loss, solar metallicity, and an age of 20 Myr, is overplotted on the CMD of field cl16.7 in Fig. 8. The isochrone is shifted to the assumed cluster distance, and reddened. It encompasses the RSGs, as well as the early-type giant star 26.

For the two RSGs, we derived luminosities of 2.3 × 104 and 1.5 × 104L (see Table 6). In Fig. 10, luminosities are plotted versus effective temperatures, and compared to theoretical stellar tracks by Ekström et al. (2012); we used tracks at solar metallicity that include rotation; the RSGs had initial masses from 10 M to 13 M.

6.2. RSGs in field cl49.3

In field cl49.3, we detected four RSGs (27, 28, 29, and 30). The four RSGs have similar AKs(AKs = 0.47 ± 0.11 mag), and spectral types (from K4 to M1.5); their Ks range from 4.1 mag to 6.4 mag. The JKs versus Ks diagram presents a gap of about 2 mag between the RSGs and other fainter stars (Fig. 8).

The O9.5-B3 star 39 is 4.8 mag fainter than the group of RSGs, and has AKs = 0.56 mag. The spectro-photometric properties of star 39 indicate a giant (DM = 14.16 mag) or a dwarf (DM = 13.02 mag). The similarity of reddening suggests that the giant 39 and the four RSGs are at the same distance, and could belong to the same coeval population; an isochrone of 20 Myr (Lejeune & Schaerer 2001) encompass both the RSGs and the colors and magnitudes of star 39.

The bulk of stars at K = 1112.5 mag and JK ≈ 1 mag is made of red clump stars (see Sect. 5 and Fig. 6). Among them we found six blue stars and star 39 (see Fig. 9). These blue stars appear radially concentrated towards the RSGs, and are unrelated to the population of clump stars; late-types and early-types have different intrinsic colors (JK = 0.7−0.8 mag for late, and 0.0 mag for early-type stars), and the red clump of stars at JK ≈ 1 mag is in the foreground of star 39.

For the assumed distance (DM = 14.2 ± 0.1 mag from clump stars), the luminosities of RSGs range from 1.6 × 104L to 1.1 × 105L. From a comparison with the tracks by Ekström et al. (2012), we estimate initial masses from 12 M to 15 M (see Fig. 10).

thumbnail Fig. 9

Panel a): (JK) versus (HK) diagram of stars with 12.5 <K< 11 mag and 0.8 <JK< 1.3 mag in the cl49.3 field. The dashed lines show the locus of a M1 and an O9 star for increasing AKs. Panel b): map of the blue stars (diamonds), along with the detected RSGs (squares) and field stars (dots).

6.3. Other fields

Field cl59.8 contains seven bright stars with Ks from 6.5 mag to 9 mag. For the estimated DM = 12.9 ± 1.4 mag, these stars have typical bolometric magnitudes of giant stars. Nevertheless, they form a statistically significant overdensity of bright stars. The CMD of field cl59.8 is shown in Fig. 8. Stars 44 and 47 are bluer (AKs = 0.06 and 0.11 mag), and are likely in the foreground. Star 45 has AKs = 0.79 mag, and is a possible background object. For the remaining four stars, 42, 43, 46, and 48, AKs range from 0.36 mag to 0.56 mag, and luminosities increase with spectral types; they are consistent with massive giants drawn from a population older than 100 Myr (Marigo et al. 2008). The infrared CMDs of simple stellar populations from 50 Myr to 200 Myr show gaps and a large number of bright giants, e.g., the GLIMPSE13 cluster (Messineo et al. 2009).

The CMDs of fields cl1.5 and cl9.5 do not present peculiar gaps or overdensities of stars brighter than Ks = 7 mag (see Table 2). In conclusion, fields cl1.5 and cl9.5 are populated by older stars (>600 Myr; Ferraro et al. 1995).

Field cl1.5 appears as a window of low extinction at the distance of the Galactic center with an inferred DM = 14.54 ± 0.16 mag, and an average AKs = 0.45 ± 0.08 mag; this result is in agreement with the extinction map by Gonzalez et al. (2012). Spectroscopically detected giant stars within <2 from the center of cl1.5 have Ks from 7.52 mag to 8.55 mag; since the tip of the red giant branch (RGB) is expected at K ≈ 8.2 mag (Messineo et al. 2005; Glass et al. 1995), they are most likely AGBs. The Mira-like star 2 is above the RGB tip (Ks = 6.86 mag). The Mira-like AGB 1 is 3.′3 from the center, and has Ks = 4.86 mag and Mbol = −6.98 mag for the assumed DM, which is compatible with the maximum luminosities predicted for Mira stars (see Fig. 20 in Vassiliadis & Wood 1993).

thumbnail Fig. 10

Luminosities versus effective temperatures of the newly detected RSGs. RSGs in field cl16.7 are marked with squares, RSGs in field cl49.3 with diamonds. Stellar tracks for stars of 9, 12, and 15 M, from the new rotating Geneva models with a solar metallicity (Ekström et al. 2012), are shown with dotted lines.

6.4. The Q1 and Q2 parameters

Messineo et al. (2012) analyzed Galactic stars of known types with the following parameters:

Q1 = (JH) − 1.8 × (HKs),

Q2 = (JKs)− 2.69 × (Ks−[ 8.0 ] ).

These two parameters are interstellar extinction free, and measure the distance of a datapoint from the Galactic reddening vector (passing trough the origin) in the JH versus HKs (Q1), and JKs versus Ks−[ 8.0 ] planes (Q2). Typically, early-type stars have Q1 < −0.20 × (KS − [ 8.0 ] ) + 0.34 mag (Messineo et al. 2012). A fraction of 72% of the spectroscopic targets are correctly classified into late- and early-type stars by applying the above relation. The Q1 and Q2 parameters only allow the selection of stars with colors typical of RSGs. Additional information on distances and luminosities is needed to classify RSGs. Messineo et al. (2012) calculated that 72% of known RSGs fall into the region 0.1 mag <Q1 < + 0.5 mag and −1.1 mag <Q2 < + 1.5 mag; five of the six new RSGs fall in this region, their average Q1 = 0.3 mag with σ = 0.12 mag; in contrast, detected Mira-like stars have an average Q1 = −0.11 mag with σ = 0.14 mag.

7. Summary and discussion

The detection of clusters rich in RSGs is strongly affected by biases. It is difficult to classify and detect RSGs, because of our poor knowledge on distance, and because of the high and patchy extinction in direction of the Galactic plane. Nevertheless, their census is of importance to understanding the history of Galactic chemical enrichment, modality of star formation, and Galactic structure.

We conducted infrared spectroscopic observations of five overdensities with bright stars, which were detected with GLIMPSE and 2MASS data along the Galactic plane. We obtained 48 stellar spectra, and assigned spectral types to 47 new stars.

Late-type stars were classified using photometric data and spectroscopic indexes (water indexes and EW(CO)s). High values of the water index showed that stars 1, 2, and 24 were Mira-like AGBs. We detected six new RSGs (22, 23, 27, 28, 29, and 30) via their large EW(CO)s, low water indexes, and luminosities (> 104L). All but one of the RSGs have EW(CO) > 42 Å.

Clustering of candidate RSGs and gaps in the CMDs are powerful tools for detecting RSGs. The detected RSGs are in two fields, cl16.7 and cl49.3. Each group of RSGs presents narrow ranges of AKs and Ks magnitudes. Stars 22 and 23 fall in field cl16.7, and stars 27, 28, 29, and 30 in field cl49.3. The CMDs of both fields show large gaps between the RSGs and the group of blue giant/supergiants extending over a main sequence; this gap is a typical feature of clusters rich in RSGs (e.g., Figer et al. 2006).

The two RSGs in field cl16.7 are located at a heliocentric distance of ~3.9 kpc, at a Galactic longitude of 16.7° on the the Scutum-Crux arm. The two RSGs had stellar initial masses from 11 M to 13 M.

The cluster of four RSGs in field cl49.3 is at a distance of ~7.0 kpc, at 493 of galactic longitude on the Sagittarius-Carina arm. The RSGs had initial masses of 12 M to 15 M.

By using multiwavelength analysis, and a measure of dust obscuration with distance, it is possible to isolate far distant and luminous objects, and groups of candidate RSGs. However, spectroscopic confirmation is required. Follow-up kinematic and variability studies will provide further constraints on the properties of the newly detected groups/clusters of RSGs. The V band magnitudes of the newly discovered RSGs range from 13.7 mag to 17.2 mag, within the photometric and spectroscopic range covered by the upcoming mission Gaia.

Online material

Table 3

Infrared counterparts of the spectroscopically observed targets.

Table 4

List of stars spectroscopically observed with UKIRT/UIST and NTT/SofI.

Table 6

Photometric properties of the targeted stars.


1

IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.

Acknowledgments

This work was partially funded by the ERC Advanced Investigator Grant GLOSTAR (247078). The material in this work was partly supported by NASA under award NNG 05-GC37G, through the Long-Term Space Astrophysics program. This research was partly performed in the Rochester Imaging Detector Laboratory with support from a NYSTAR Faculty Development Program grant. Qingfeng Zhu thanks for the support by the Fundamental Research Funds for the Central Universities from Educational Department of China, #WK2030220001. M.M. and V.I. thank ESO for supporting a working visit of Dr. Ivanov at RIT. M.M. and Z.Q. acknowledge an awarded support from ESA/ESTEC to host Dr. Zhu. 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 National Science Foundation. This work is based [in part] on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. DENIS is the result of a joint effort involving human and financial contributions of several Institutes mostly located in Europe. It has been supported financially mainly by the French Institut National des Sciences de l’Univers, CNRS, and French Education Ministry, the European Southern Observatory, the State of Baden-Wuerttemberg, and the European Commission under a network of the Human Capital and Mobility program. This research made use of data products from the Midcourse Space Experiment, the processing of which was funded by the Ballistic Missile Defence Organization with additional support from the NASA office of Space Science. This research has made use of the SIMBAD data base, operated at CDS, Strasbourg, France. This publication makes use of data products from WISE, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. We are thankful to the UKIRT team and SofI team for a great support during the observing runs, to the UKIDSS team in Cambridge for providing detailed information on the pipeline products, and to Dr. Stetson for providing the Daophot code, and helping with the installation. M.M. thanks Harm Habing and Ed Churchwell for discussions on GLIMPSE data, and J. Borissova, R. Kurtev, Saurabh Sharma, and G. R. Ivanov for discussions on stellar overdensities. M.M. is grateful to Jos de Bruine and Timo Prusti for useful discussions and support during her ESA fellowship. We thank the referee Ignacio Negueruela for his careful reading of our manuscript.

References

  1. Alard, C., Blommaert, J. A. D. L., Cesarsky, C., et al. 2001, ApJ, 552, 289 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alvarez, R., Lançon, A., Plez, B., & Wood, P. R. 2000, A&A, 353, 322 [NASA ADS] [Google Scholar]
  3. Alves, D. R. 2000, ApJ, 539, 732 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  4. Babusiaux, C., & Gilmore, G. 2005, MNRAS, 358, 1309 [NASA ADS] [CrossRef] [Google Scholar]
  5. Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953 [NASA ADS] [CrossRef] [Google Scholar]
  6. Benson, P. J., Little-Marenin, I. R., Woods, T. C., et al. 1990, ApJS, 74, 911 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 [NASA ADS] [Google Scholar]
  8. Bibby, J. L., Crowther, P. A., Furness, J. P., & Clark, J. S. 2008, MNRAS, 386, L23 [NASA ADS] [CrossRef] [Google Scholar]
  9. Blum, R. D., Ramírez, S. V., Sellgren, K., & Olsen, K. 2003, ApJ, 597, 323 [NASA ADS] [CrossRef] [Google Scholar]
  10. Borissova, J., Bonatto, C., Kurtev, R., et al. 2011, A&A, 532, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Borissova, J., Kumar, M. S. N., Amigo, P., et al. 2014, AJ, 147, 18 [NASA ADS] [CrossRef] [Google Scholar]
  12. Buta, R. 1995, ApJS, 96, 39 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cabrera-Lavers, A., González-Fernández, C., Garzón, F., Hammersley, P. L., & López-Corredoira, M. 2008, A&A, 491, 781 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Carpenter, J. M. 2001, AJ, 121, 2851 [NASA ADS] [CrossRef] [Google Scholar]
  15. Chené, A.-N., Borissova, J., Bonatto, C., et al. 2013, A&A, 549, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Clampitt, L., & Burstein, D. 1997, AJ, 114, 699 [NASA ADS] [CrossRef] [Google Scholar]
  17. Clark, J. S., Negueruela, I., Davies, B., et al. 2009, A&A, 498, 109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Comerón, F., Torra, J., Chiappini, C., et al. 2004, A&A, 425, 489 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Cordes, J. M., & Lazio, T. J. W. 2002 [arXiv:astro-ph/0207156] [Google Scholar]
  20. Davies, B., Figer, D. F., Kudritzki, R.-P., et al. 2007, ApJ, 671, 781 [NASA ADS] [CrossRef] [Google Scholar]
  21. Davies, B., de LaFuente, D., Najarro, F., et al. 2012, MNRAS, 419, 1860 [NASA ADS] [CrossRef] [Google Scholar]
  22. de Vaucouleurs, G. 1964, in The Galaxy and the Magellanic Clouds, ed. F. J. Kerr, IAU Symp., 20, 88 [Google Scholar]
  23. de Wit, W. J., Testi, L., Palla, F., & Zinnecker, H. 2005, A&A, 437, 247 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Drimmel, R., & Spergel, D. N. 2001, ApJ, 556, 181 [NASA ADS] [CrossRef] [Google Scholar]
  25. Drimmel, R., Cabrera-Lavers, A., & López-Corredoira, M. 2003, A&A, 409, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Dutra, C. M., Santiago, B. X., & Bica, E. 2002, A&A, 381, 219 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Dutra, C. M., Bica, E., Soares, J., & Barbuy, B. 2003, A&A, 400, 533 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Egan, M. P., Price, S. D., Kraemer, K. E., et al. 2003, VizieR Online Data Catalog: V/114 [Google Scholar]
  29. Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Epchtein, N., Deul, E., Derriere, S., et al. 1999, VizieR Online Data Catalog: II/2 [Google Scholar]
  31. Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 [NASA ADS] [CrossRef] [Google Scholar]
  32. Ferraro, F. R., FusiPecci, F., Testa, V., et al. 1995, MNRAS, 272, 391 [NASA ADS] [CrossRef] [Google Scholar]
  33. Figer, D. F., MacKenty, J. W., Robberto, M., et al. 2006, ApJ, 643, 1166 [NASA ADS] [CrossRef] [Google Scholar]
  34. Froebrich, D., Scholz, A., & Raftery, C. L. 2007, MNRAS, 374, 399 [Google Scholar]
  35. Frogel, J. A., & Whitford, A. E. 1987, ApJ, 320, 199 [NASA ADS] [CrossRef] [Google Scholar]
  36. Glass, I. S., Whitelock, P. A., Catchpole, R. M., & Feast, M. W. 1995, MNRAS, 273, 383 [NASA ADS] [CrossRef] [Google Scholar]
  37. Gonzalez, O. A., Rejkuba, M., Zoccali, M., Valenti, E., & Minniti, D. 2011, A&A, 534, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Gonzalez, O. A., Rejkuba, M., Zoccali, M., et al. 2012, A&A, 543, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. González-Fernández, C., & Negueruela, I. 2012, A&A, 539, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Groenewegen, M. A. T. 2008, A&A, 488, 935 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Groenewegen, M. A. T., Sloan, G. C., Soszyński, I., & Petersen, E. A. 2009, A&A, 506, 1277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Habing, H. J., & Olofsson, H. 2003, Asymptotic giant branch stars, eds. J. C. Lattanzio, & P. R. Wood, Astronomy and astrophysics library (New York, Berlin: Springer), Chap. 2 [Google Scholar]
  43. Habing, H. J., Sevenster, M. N., Messineo, M., van de Ven, G., & Kuijken, K. 2006, A&A, 458, 151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Hanson, M. M., Conti, P. S., & Rieke, M. J. 1996, ApJS, 107, 281 [Google Scholar]
  45. Hanson, M. M., Rieke, G. H., & Luhman, K. L. 1998, AJ, 116, 1915 [NASA ADS] [CrossRef] [Google Scholar]
  46. Hodgkin, S. T., Irwin, M. J., Hewett, P. C., & Warren, S. J. 2009, MNRAS, 394, 675 [NASA ADS] [CrossRef] [Google Scholar]
  47. Houk, N., & Smith-Moore, M. 1988, Michigan Catalogue of Two-dimensional Spectral Types for the HD Stars, Vol. 4, Declinations −26.̊0 to −12.̊0. [Google Scholar]
  48. Humphreys, R. M., & McElroy, D. B. 1984, ApJ, 284, 565 [NASA ADS] [CrossRef] [Google Scholar]
  49. Ivanov, V. D., Borissova, J., Pessev, P., Ivanov, G. R., & Kurtev, R. 2002, A&A, 394, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Ivanov, V. D., Rieke, M. J., Engelbracht, C. W., et al. 2004, ApJS, 151, 387 [NASA ADS] [CrossRef] [Google Scholar]
  51. Ivanov, V. D., Messineo, M., Zhu, Q., et al. 2010, in IAU Symp. 266, eds. R. de Grijs, & J. R. D. Lépine, 203 [Google Scholar]
  52. Kleinmann, S. G., & Hall, D. N. B. 1986, ApJS, 62, 501 [NASA ADS] [CrossRef] [Google Scholar]
  53. Koornneef, J. 1983, A&A, 128, 84 [NASA ADS] [Google Scholar]
  54. Lancon, A., & Rocca-Volmerange, B. 1992, A&AS, 96, 593 [NASA ADS] [Google Scholar]
  55. Lejeune, T., & Schaerer, D. 2001, A&A, 366, 538 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Levesque, E. M., Massey, P., Olsen, K. A. G., et al. 2005, ApJ, 628, 973 [NASA ADS] [CrossRef] [Google Scholar]
  57. Lucas, P. W., Hoare, M. G., Longmore, A., et al. 2008, MNRAS, 391, 136 [NASA ADS] [CrossRef] [Google Scholar]
  58. Marigo, P. 2001, A&A, 370, 194 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Marigo, P., Girardi, L., Bressan, A., et al. 2008, A&A, 482, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Matsuura, M., Yamamura, I., Cami, J., Onaka, T., & Murakami, H. 2002, A&A, 383, 972 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Mercer, E. P., Clemens, D. P., Meade, M. R., et al. 2005, ApJ, 635, 560 [NASA ADS] [CrossRef] [Google Scholar]
  62. Messineo, M. 2004, Ph.D. Thesis, Leiden Observatory, Leiden University, The Netherlands [Google Scholar]
  63. Messineo, M., Habing, H. J., Menten, K. M., Omont, A., & Sjouwerman, L. O. 2004, A&A, 418, 103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Messineo, M., Habing, H. J., Menten, K. M., et al. 2005, A&A, 435, 575 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Messineo, M., Davies, B., Ivanov, V. D., et al. 2009, ApJ, 697, 701 [NASA ADS] [CrossRef] [Google Scholar]
  66. Messineo, M., Davies, B., Figer, D. F., et al. 2011, ApJ, 733, 41 [NASA ADS] [CrossRef] [Google Scholar]
  67. Messineo, M., Menten, K. M., Churchwell, E., & Habing, H. 2012, A&A, 537, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Messineo, M., Menten, K. M., Figer, D., et al. 2014 [arXiv:1404.7297] [Google Scholar]
  69. Moorwood, A., Cuby, J.-G., & Lidman, C. 1998, The Messenger, 91, 9 [NASA ADS] [Google Scholar]
  70. Negueruela, I., González-Fernández, C., Marco, A., Clark, J. S., & Martínez-Núñez, S. 2010, A&A, 513, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Negueruela, I., González-Fernández, C., Marco, A., & Clark, J. S. 2011, A&A, 528, A59 [Google Scholar]
  72. Ortiz, R., Blommaert, J. A. D. L., Copet, E., et al. 2002, A&A, 388, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Paczyński, B. 1971, Acta Astron., 21, 417 [NASA ADS] [Google Scholar]
  74. Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 [NASA ADS] [CrossRef] [Google Scholar]
  75. Ramsay Howat, S. K., Ettedgui-Atad, E., Bennett, R. J., et al. 1998, in SPIE Conf. Ser. 3354, ed. A. M. Fowler, 456 [Google Scholar]
  76. Rayner, J. T., Cushing, M. C., & Vacca, W. D. 2009, ApJS, 185, 289 [NASA ADS] [CrossRef] [Google Scholar]
  77. Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618 [NASA ADS] [CrossRef] [Google Scholar]
  78. Schultheis, M., Ganesh, S., Glass, I. S., et al. 2000, A&A, 362, 215 [NASA ADS] [Google Scholar]
  79. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  80. Soto, M., Barbá, R., Gunthardt, G., et al. 2013, A&A, 552, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Stetson, P. B. 1987, PASP, 99, 191 [NASA ADS] [CrossRef] [Google Scholar]
  82. Vassiliadis, E., & Wood, P. R. 1993, ApJ, 413, 641 [NASA ADS] [CrossRef] [Google Scholar]
  83. Wainscoat, R. J., Cohen, M., Volk, K., Walker, H. J., & Schwartz, D. E. 1992, ApJS, 83, 111 [NASA ADS] [CrossRef] [Google Scholar]
  84. Wallace, L., & Hinkle, K. 1996, ApJS, 107, 312 [NASA ADS] [CrossRef] [Google Scholar]
  85. Wegner, W. 1994, MNRAS, 270, 229 [NASA ADS] [CrossRef] [Google Scholar]
  86. Wood, P. R., Bessell, M. S., & Fox, M. W. 1983, ApJ, 272, 99 [NASA ADS] [CrossRef] [Google Scholar]
  87. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
  88. Yang, M., & Jiang, B. W. 2011, ApJ, 727, 53 [NASA ADS] [CrossRef] [Google Scholar]
  89. Zacharias, N., Monet, D. G., Levine, S. E., et al. 2004, BAAS, 36, 1418 [Google Scholar]

All Tables

Table 1

List of candidate cluster positions, apparent radii, AKs, and modulus distances (DM).

Table 2

Counts of stars with Ks> 10.0, 8.5, and 7.0 mag in the fields listed in Table 1, and, as a comparison, in control fields of equal area.

Table 5

Spectroscopic indexes. Spectral regions used as features and continuum are specified.

Table 3

Infrared counterparts of the spectroscopically observed targets.

Table 4

List of stars spectroscopically observed with UKIRT/UIST and NTT/SofI.

Table 6

Photometric properties of the targeted stars.

All Figures

thumbnail Fig. 1

GLIMPSE 3.6 μm images of the candidate clusters listed in Table 1: cl1.5, cl9.5, cl16.5, cl49.3, and cl59.8. Spectroscopic targets are labeled as in Tables 3 and 4. Diamonds indicate red giants, squares Mira-like AGB stars, crosses RSG stars, and triangles early- and yellow-type stars. The coordinates of the image centers are in degrees.

In the text
thumbnail Fig. 2

Long-K spectra taken with UIST. Identification numbers refer to Table 4. All are late-type stars, with the exception of star 13.

In the text
thumbnail Fig. 3

Low-resolution H and K spectra taken with SofI. Identification numbers are taken from Table 4. Late-type stars are shown in panel a). Early-type stars are shown in panel b) together with a few low-resolution template spectra from the library of Lancon & Rocca-Volmerange (1992), and medium-resolution template spectra from the IRTF library (Rayner et al. 2009).

In the text
thumbnail Fig. 4

Medium-resolution K spectra taken with SofI. Identification numbers are taken from Table 4. All stars are late-type stars, with the exception of star 39.

In the text
thumbnail Fig. 5

H2O indexes versus the EW(CO)s of targeted late-type stars and comparison samples of AGBs (Mira and semi-regular variables) and RSGs from the IRTF library (Rayner et al. 2009).

In the text
thumbnail Fig. 6

Top panels: interstellar extinction AK versus distance moduli of red clump stars in the fields of cl16.7 (left) and cl49.3 (right). As a comparison, the model of Drimmel et al. (2003) is shown with a dashed line. Bottom panels: UKIDSS JK versus K diagram of data-points within an area of 10× 10 from the estimated centers of cl16.7 (left) and cl49.3 (right).

In the text
thumbnail Fig. 7

XY view of the Milky Way. The Galactic center is at (0, 0) and the Sun is at (8, 0). The dashed lines indicate the line of sights to cl16.7 and cl49.3, which are marked with squares. A sketch of the spiral structure is also drawn and the arms are labeled (Cordes & Lazio 2002). The X-axis is oriented along the Sun-Galactic center line of sight.

In the text
thumbnail Fig. 8

2MASS/UKIDSS JKS versus KS CMDs: cl1.5 (top left), cl9.5 (top right), cl16.3 (middle left), cl49.3 (middle right), and cl59.8 (bottom left). Photometry was obtained with a PSF-fitting technique (see text). For each candidate cluster, a diagram of the candidate cluster is shown in the left panel, and that of a control field of equal area in the right panel. Spectroscopically detected early and yellow stars are marked with filled triangles, Miras with squares, RSGs with crosses, and red giants with diamonds; identification numbers are taken from Table 3. In the CMDs of fields cl16.3 and cl49.3, which are rich in RSGs, a continuous curve displays an isochrone of 20 Myr from Lejeune & Schaerer (2001), reddened with AKs = 0.50 mag and 0.48 mag, and shifted to DM = 12.93 mag and 14.22 mag, respectively; a dashed line shows the sequence of red clump stars (see Sect. 4.2). The CMD of cl59.8 shows three isochrones from Marigo et al. (2008) with solar metallicity and ages of 100 (dashed line), 200 (continuous line), and 300 Myr (dotted line), which were reddened (AKs = 0.55 mag) and shifted to DM = 12.9 mag.

In the text
thumbnail Fig. 9

Panel a): (JK) versus (HK) diagram of stars with 12.5 <K< 11 mag and 0.8 <JK< 1.3 mag in the cl49.3 field. The dashed lines show the locus of a M1 and an O9 star for increasing AKs. Panel b): map of the blue stars (diamonds), along with the detected RSGs (squares) and field stars (dots).

In the text
thumbnail Fig. 10

Luminosities versus effective temperatures of the newly detected RSGs. RSGs in field cl16.7 are marked with squares, RSGs in field cl49.3 with diamonds. Stellar tracks for stars of 9, 12, and 15 M, from the new rotating Geneva models with a solar metallicity (Ekström et al. 2012), are shown with dotted lines.

In the text

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.