A&A 453, 535-545 (2006)
DOI: 10.1051/0004-6361:20052891

Recent star formation in the inner Galactic Bulge seen by ISOGAL[*],[*],[*],[*]

I. Classification of bright mid-IR sources in a test field

F. Schuller1,2 - A. Omont1 - I. S. Glass3 - M. Schultheis1,4 - M. P. Egan5 - S. D. Price6


1 - Institut d'Astrophysique de Paris, CNRS and Paris-VI University, 98bis Bd. Arago, 75014 Paris, France
2 - Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
3 - South African Astronomical Observatory, PO Box 9, Observatory 7935, South Africa
4 - Observatoire de Besançon, CNRS, 41bis Av. de l'Observatoire, BP 1615, 25010 Besançon Cedex, France
5 - Air Force Research Laboratory, ODASD, 1851 s. Bell St., CM3, Suite 7000, Arlington, VA 22202, USA
6 - Air Force Research Laboratory, Space Vehicles Directorate, 29 Randolph Road, Hanscom AFB, MA 01731, USA

Received 16 February 2005 / Accepted 23 February 2006

Abstract
Context. The stellar populations in the central region of the Galaxy are poorly known because of the high visual extinction and very great source density in this direction.
Aims. To use recent infrared surveys for studying the dusty stellar objects in this region.
Methods. We analyse the content of a ${\sim}20\times20$ arcmin2 field centred at $(l,b)=(-0.27^{\circ},-0.06^{\circ})$ observed at 7 and 15 $\mu $m as part of the ISOGAL survey. These ISO observations are more than an order of magnitude better in sensitivity and spatial resolution than the IRAS observations. The sources are cross-associated with other catalogues to identify various types of objects. We then derive criteria to distinguish young objects from post-main sequence stars.
Results. We find that a sample of about 50 young stellar objects and ultra-compact H II regions emerges, out of a population of evolved AGB stars. We demonstrate that the sources colours and spatial extents, as they appear in the ISOGAL catalogue, possibly complemented with MSX photometry at 21 $\mu $m, can be used to determine whether the ISOGAL sources brighter than 300 mJy at 15 $\mu $m (or [15] $\leq$ 4.5 mag) are young objects or late-type evolved stars.

Key words: Galaxy: centre - stars: formation - stars: variables: general - infrared: stars

   
1 Introduction

The ISOGAL survey is a large set of observations of the inner Galaxy at 7 and 15 $\mu $m (Omont et al. 2003; Schuller et al. 2003), obtained with the ISOCAM instrument on board the ISO satellite (Cesarsky et al. 1996). The ISOGAL fields are distributed along the Galactic plane, mostly in the $\vert b \vert \la 1^{\circ}$ range. The achieved sensitivity is about 10 mJy in most fields in the Galactic Disk, and of order 20-30 mJy in the fields observed with narrow filters and 3'' pixel, including the Galactic Centre region, and thus the observations discussed in the present paper. The catalogues of more than 105extracted point sources have been systematically cross-identified with the near infrared DENIS data (Epchtein et al. 1997), resulting in the ISOGAL-DENIS Point Source Catalogue (PSC), which provides measurements at 7 and/or 15 $\mu $m, together with up to three near infrared bands (I, J and $K_{\rm s}$). The astrometric accuracy is better than 0.5'' for sources with a DENIS counterpart, and of order 2'' for sources without.


  \begin{figure}
\par {\includegraphics[width=16.8cm,clip]{2891f1.eps} }\end{figure} Figure 1: ISOCAM raster images used for the FC-027 field: LW5 (7 $\mu $m, left) and LW9 (15 $\mu $m, right). The rectangular frames show the limits of the "regular'' field, as defined in Schuller et al. (2003) - see also Eq. (1). The dotted lines show the limits of the area included in the Glass et al. (2001) survey (see Sect. 3.1): only the area to the left of these lines was observed by Glass et al. (2001). North and east directions are indicated.
Open with DEXTER

The analysis of various individual ISOGAL fields in the Galactic disk and bulge have shown that the bulk of the detected stellar population is composed of AGB and late RGB stars (see references in Omont et al. 2003). They show a rather well defined sequence in the [15] vs. [7]-[15] colour magnitude diagram, in good agreement with the colours expected for such objects (Fig. 2). On the other hand, a few thousand ISOGAL sources can be interpreted as young stellar objects or young massive stars (hereafter young objects) deeply embedded in dusty envelopes or thick disks. Theoretical models, as well as a systematic analysis of ISOGAL images and catalogues, enabled us (Felli et al. 2000, 2002; Schuller 2002) to define criteria to roughly distinguish candidate young objects from the bulk of post-main sequence stars in the [15] vs. [7]-[15] diagram.

The ISOGAL FC-00027-00006 field (hereafter called FC-027), centred at $(l,b)=(-0.27^{\circ},-0.06^{\circ})$, is located very close to the Galactic Centre, between Sgr A* and Sgr C, in a region of strong diffuse emission at infrared and radio wavelengths. The visual extinction in this field ranges from 20 to more than 35 mag (Schultheis et al. 1999), and the photometry must be corrected carefully in order to interpret the observed magnitudes and colours properly, especially in the near infrared. However, the [7]-[15] colours cannot be affected by more than $\sim$1 mag of interstellar extinction (see Fig. 2).

We will show in the present paper that, in spite of the extinction, the nature of the bright sources (with magnitudes at 15 $\mu $m [15] $\leq$ 4.5, equivalent to flux densities above 280 mJy) can be inferred from the raw data, as they appear in the ISOGAL PSC. The bulk of the brightest sources can be identified with various classes of objects (long period variable stars, OH/IR stars, ultra-compact H II regions). The mid-infrared properties (especially the [7]-[15] colour and spatial extension at 15 $\mu $m) of young objects (YSOs and H II regions) are quite different from those observed for late-type stars. This leads to robust selection criteria for extracting the young population from the ISOGAL PSC (see also Lumsden et al. 2002, for a similar analysis based on the MSX data). These criteria will be used to make a census of such objects in all ISOGAL fields covering the inner Galactic Bulge in a forthcoming paper, where we will estimate their luminosities, taking into account extinction corrections.

The ISOGAL data are briefly presented in Sect. 2. We then show in Sect. 3 that most ISOGAL bright sources can be associated with known objects, in particular long period variable stars, radio-sources, or infrared sources from the IRAS and MSX point source catalogues. Our main results, consisting of criteria to select young sources from the ISOGAL catalogue, are summarised and used to extract luminous young candidate objects without known counterparts from the whole FC-027 field in Sect. 4.

  
2 ISOGAL observations

A complete description of the ISOGAL observations, their processing and the creation of the Point Source Catalogue (PSC) can be found in Schuller et al. (2003). In the present paper, we analyse the source content of the FC-027 field. We confine our analysis to the formal limits of this field, or the "regular'' catalogue (as opposed to the "edge'' catalogue, see Schuller et al. 2003), which covers $18.5\times17$ arcmin2 (0.09 deg2) centred at $(l,b)=(-0.27^{\circ},-0.06^{\circ})$. Namely, the present paper deals with the area defined by:

 \begin{displaymath}
\left\{ \begin{array}{l}
-0.424 \; \leq \; l \; \leq \; -0.115 \\
-0.194 \; \leq \; b \; \leq \; 0.089.
\end{array} \right.
\end{displaymath} (1)

This field was observed with the narrow filters LW5 (6.5-7.1 $\mu $m) and LW9 (13.9-15.9 $\mu $m), using pixels of 3''field of view (see the raster images in Fig. 1). The 7 $\mu $m image was observed on 1996, Sep. 24, while the 15 $\mu $m observation was done on 1998, Feb. 19, so that variability can affect the observed colours. Point sources were extracted from the 7 $\mu $m and 15 $\mu $m observations down to magnitudes 8.4 at 7 $\mu $m (35 mJy), and 7.0 at 15 $\mu $m (28 mJy), corresponding to roughly 50% completeness (Schuller et al. 2003). About 90% of the sources also have DENIS (Epchtein et al. 1997) near-infrared counterparts, with $K_{\rm s}$ magnitudes brighter than 11.3. Fainter DENIS sources were not considered in order to reduce the fraction of chance associations, so that a few faint near-infrared counterparts may have been missed. Altogether, we estimated that no more than $\sim$1% of ISO-DENIS associations should be wrong.
  \begin{figure}
\par {\includegraphics[width=7.6cm,clip]{2891f2.eps} }\end{figure} Figure 2: [15] vs. [7]-[15] colour magnitude diagram for the FC-027 field. Sources located inside the area also observed by Glass et al. (2001) are shown with black symbols, while those outside the Glass et al. (2001) area are indicated with gray symbols. The sources associated with LPVs are indicated by plus signs, OH/IR stars are shown with crosses, and sources where an 86 GHz SiO maser was detected are shown with square symbols. The thick lines show the loci of giant stars (continuous line), supergiants (dotted line), M-type AGB stars and L-type stars with different mass-losses (dashed line) in this diagram, as computed using synthetic spectra (Cohen, private communication). The extinction vector in the bottom right corner corresponds to $A_{\rm V}=30$ mag; it was computed using the mean values of $A_\lambda /A_{\rm V}$ for the LW5 and LW9 filters derived by Jiang et al. (2006).
Open with DEXTER

Table 1: Associations between uncatalogued radio sources visible in the 1616 MHz map of Liszt & Spiker (1995) and bright ISOGAL sources. The 1950 coordinates (Cols. 1 and 2) are extracted from Fig. 1 in Liszt and Spiker (1995), and were precessed to 2000 epoch in Cols. 3 and 4.

  
3 Identification of known sources

In this section, we will show how ISOGAL point sources can be associated with sources from various catalogues, taken from the abundant literature dealing with the centre of the Galaxy. In particular, we will focus on: long period variable stars (LPVs) and OH/IR stars (Sect. 3.1), radio sources detected in the cm continuum (Sect. 3.2), the IRAS Point Source Catalogue (Sect. 3.3) and the MSX Point Source Catalogue (Sect. 3.4). Additional results derived from specific near infrared spectroscopic observations are discussed in Sect. 3.5. The main results outlined from this compilation of data are summarised in Sect. 3.6.

  
3.1 Late evolved stars

As in the other ISOGAL fields, the huge majority of the detected point sources in the FC-027 field are luminous late-type giants toward the end of the red giant phase, mostly in the AGB phase (Omont et al. 1999; Glass et al. 1999; Ojha et al. 2003; van Loon et al. 2003). Their photospheres are bright emitters in the near-infrared, and many of them are surrounded by dusty circumstellar envelopes due to mass loss. The dust grains absorb the stellar radiation and re-emit in the mid-infrared, so that these stars are bright in mid-infrared images. Most of these sources are variable (Alard et al. 2001), and we have been able to find ISOGAL counterparts to variable stars of different classes (Miras, OH/IR, other long period variables).

Glass et al. (2001) have monitored an area of $24\times24$ arcmin2centred at the Galactic Centre over four years, in the K band, in search of variable stars. Their surveyed area overlaps the FC-027 field by $\sim$0.042 deg2, or nearly one half the ISOGAL field (see Fig. 1). In this overlap area, they found 110 long period variable stars (hereafter LPVs).

Among these 110 LPVs, 97 are associated with ISOGAL sources within a 5''search radius, 89 of which are detected at both 7 and 15 $\mu $m. These have [7]-[15] colours mostly around 1 mag, along a sequence spanning the range 3 to 6 mag at 15 $\mu $m (Fig. 2). However, a few 0.1 mag in the observed colours may be due to variability, because of the 1.5 year delay between the 7 $\mu $mand 15 $\mu $m observations. Blending effects may also introduce errors in the [7]-[15] colours. Most of the 13 LPVs without association appear as faint objects at 7 $\mu $m. Confusion with the bright diffuse background, added to the high source density in this field, explain these peculiar cases. In addition, 10 out of these 13 sources have mean K-band magnitudes above 9.5 (Glass et al. 2001), while only one third of those with an ISOGAL counterpart have such faint near-IR magnitudes.

Three catalogues of OH/IR stars (Lindqvist et al. 1992; Sevenster et al. 1997; Sjouwerman et al. 1998) were compiled and associated with ISOGAL data in the inner Galactic Bulge by Ortiz et al. (2002). They found 37 OH/IR stars in the FC-027 field, 34 of which are associated with ISOGAL sources within 6''. As shown with crosses in Fig. 2, they are found at the brightest end of the main sequence in this diagram, with magnitudes [15] $\leq$ 4.5 and colours 0.8 $\leq$ [7]-[15] $\leq$ 2.3. These are typical colours and magnitudes for OH/IR stars with bolometric magnitudes in the range $-5.0 \pm 1.0$ mag and mass-loss rates between ${\sim} 10^{-6}$ and a few 10-5 $M_{\odot}$/yr (Ortiz et al. 2002; Groenewegen 1993).

Another seven sources (indicated with square symbols in Fig. 2) were detected in the 86 GHz SiO maser line by Messineo et al. (2002), who observed 12 sources selected in the FC-027 field. This detection rate of 58% is consistent with the 48% detection rate found by Deguchi et al. (2004), who surveyed the complete sample of LPV stars published by Glass et al. (2001) in SiO 43 GHz lines.

  
3.2 Radio continuum sources

Large fractions of the Galactic Disk, including the central region, have been surveyed at radio wavelengths by different teams, using the VLA at 1.4 GHz (Zoonematkermani et al. 1990), 5 GHz (Becker et al. 1994), and 1.6 GHz (Liszt 1985; Liszt & Spiker 1995). In addition, Downes et al. (1979) conducted pointed observations in the Galactic Centre region with several telescopes.

A number of compact or small extended radio continuum sources from these catalogues are found in the FC-027 field. We consider as their infrared counterparts the nearest ISOGAL neighbour brighter than 5 mag at 15 $\mu $m. Infrared counterparts are found within 10'' for all sources, except for source number 3 in Downes et al. (1979) for which the separation is 18'', and source number 7 in Liszt & Spiker (1995) which coincides with a compact group of ISOGAL sources discussed in Sect. 4.3.

In addition to these catalogued sources, a number of faint radio sources can be seen in the 1616 MHz map of Liszt & Spiker (1995, see their Fig. 1). They were not included in their table, but some of them can be associated with bright ISOGAL sources, as indicated in Table 1. The positions of all these radio sources in the [15] vs. [7]-[15] colour magnitude diagram are shown in Fig. 3.


  \begin{figure}
\par {\includegraphics[width=7.6cm,clip]{2891f3.eps} }\end{figure} Figure 3: [15] vs. [7]-[15] colour magnitude diagram for all sources with 7 $\mu $m and 15 $\mu $m detections and reliable 7-15 $\mu $m associations in the FC-027 field. The sources associated with radio sources or with H II regions detected by IRAS are indicated with various symbols, as indicated in the top left corner.
Open with DEXTER

The sensitivity of the Becker et al. (1994) catalogue in the central region of the Galaxy is of order 40 mJy, and that of the Downes et al. (1979) survey is close to 10 mJy. Liszt & Spiker (1995) report a noise level of 0.7 mJy/beam in their data, but their observations cover only a fraction of the FC-027 field. In addition, all these radio data show extended emission at various positions in the field, so that confusion may limit the actual sensitivity of these surveys. It is clear that any H II region powered by a single B-type star at the distance of the Galactic Centre would have been missed by these radio observations.

  
3.3 IRAS sources

The IRAS observations were highly limited by source confusion in the very dense regions of the inner Galactic Bulge, so that no clear and unique association can be found between the IRAS and ISOGAL catalogues in most cases. We therefore looked for ISOGAL counterparts at least within the 3$\sigma $ uncertainty ellipse of each IRAS point source. In most cases, both the positions and the photometry of the IRAS sources probably result from a combination of several infrared sources, on top of diffuse emission (Fig. 4). We will thus consider the brightest source at 15 $\mu $m as the most likely main counterpart to the source responsible for the emission detected by IRAS.

Eleven IRAS PSC sources are located in this field. Four of them (17418-2914, 17419-2907, 17420-2902 and 17422-2911) have F25/F12 and F60/F25 flux ratios all higher than or close to 10, making these sources very likely to be associated with compact H II regions (Wood and Churchwell 1989). Two other sources (17416-2916 and 17417-2904), with F25/F12 > 3 and no detection at 60 $\mu $m, could well be compact H II regions. The ISOGAL counterparts to these ultra-compact H II regions show peculiar colours: they are very bright at 15 $\mu $m, as was expected from their detection by IRAS, and very red in [7]-[15] (see triangle symbols in Fig. 3).

  
3.4 MSX sources

The Midcourse Space Experiment (Price et al. 2001) performed observations in the mid-infrared of the complete Galactic plane, in the $\vert ~ b ~ \vert ~ \leq 5^{\circ}$range, in six bands centred on wavelengths between 4 and 21 $\mu $m, at a spatial resolution of ${\sim} 18''$. Here we will use version 2.3 of the MSX PSC (Egan et al. 2003). It was built by extracting point sources from co-added images, resulting in a 400 000 source catalogue with a sensitivity of about 40 mJy in the A band (8.3 $\mu $m). In particular, three raster-scanned experiments were combined at the Galactic Centre, greatly improving the sensitivity. However, as many as 10% of the point sources extracted from the densest region of the Galactic Centre may actually be diffuse emission on the scale of the MSX pixel rather than true point sources, as we will show below.

  \begin{figure}
\par {\includegraphics[width=8cm,clip]{2891f4.eps} }\end{figure} Figure 4: Positions of IRAS point sources found in the FC-027 field. The gray scale background shows the 15 $\mu $m ISO image. Every IRAS source is shown with its name, and the 3$\sigma $ uncertainty ellipse. The IRAS field of view at 12 and 25 $\mu $m ( $0.76\times 4.55$ arcmin2) is also shown with dotted lines.
Open with DEXTER

We find 333 MSX sources in the FC-027 field. After adding small offsets $(\Delta \alpha , \Delta \delta)
= (+1.44'',$ +0.32'')to the MSX coordinates, 257 sources (77%) can be associated with ISOGAL sources within an 8'' search radius. Two among the remaining 76 are located near the edge of the field and can be associated with ISO sources outside the limits. A few other ones, including the two brightest ones at 15 $\mu $m among the non-associated MSX sources, are found at the barycentre of two or more bright sources that appear resolved with ISO. This is clearly an effect of the higher spatial resolution of ISO as compared to the 18'' pixels of MSX.

However, at least one ISOGAL point source can always be found within 30'' of any MSX source. We have performed a $\chi^2$-analysis, similar to that used to cross-corelate the MSX and 2MASS catalogues in the Large Magellanic Cloud (Egan et al. 2001), to estimate the probability that an MSX source and its closest ISOGAL neighbour can be real matches. It turned out that, among the 76 sources with no counterpart at less than 8'', 55 associations have $\chi^2$ values below 18.6, where 99.99% of the true matches must lie. All of them have ISO-MSX separations less than 18'', i.e. are within one MSX pixel, so that these associations may be real, except for about ten of them that show ISO-MSX magnitude differences greater than 2 mag. Altogether, about one half of the latter associations with separations greater than 8'' may still be real. Most of the 76 MSX sources without ISOGAL associations within 8'' are located in regions of diffuse emission, and the other half may well correspond to patches of diffuse emission rather than genuine point sources.

To estimate the fraction of chance associations between MSX and ISOGAL, we also performed cross-identifications between the ISOGAL PSC and the MSX catalogue, adding an arbitrary offset to the coordinates of the latter. The results show that nearly one third of the MSX sources can be randomly associated with ISO sources within an 8''search radius, as expected from the average ISO source density. However, when taking into account the fact that a fraction of these sources are well associated with ISOGAL, we estimate that no more than 40 ISOGAL-MSX associations (out of 257) should be false. Also note that the number of MSX sources found within the FC-027 field corresponds to 0.1 source/pixel in average, so that the extraction was certainly limited by confusion. Therefore, the reality of the weakest sources and their associations with ISOGAL should be regarded with extreme caution.

Nevertheless, we find a fair general agreement between the ISO and MSX magnitudes, both at 7 and at 15 $\mu $m, as shown in Fig. 5, although with quite large dispersion partly due to the variability of some sources, and to a few wrong associations for the faintest sources. The distributions of magnitude differences show clear tails toward positive differences, that can be explained by MSX including more flux from diffuse extended emission, or from several point sources resolved by ISO, while only the magnitudes of the nearest ISO counterparts have been taken into account. In addition, the false ISO-MSX associations, as well as the Malmquist bias on the faintest MSX sources, certainly contribute to the tail in that direction (MSX magnitudes seem brighter than ISOGAL), due to the lower sensitivity of MSX compared to ISOGAL.


  \begin{figure}
\par {\includegraphics[width=7.3cm,clip]{2891f5.eps} }\end{figure} Figure 5: Distributions of differences in magnitudes between ISOGAL and MSX: ISO 7 $\mu $m-MSX A-band ( top), ISO 15 $\mu $m-MSX D-band ( middle), and ISO 7-15 colour-MSX A-D colour ( bottom).
Open with DEXTER

The MSX survey provides photometry up to 21 $\mu $m (E band, 18.2-25.1 $\mu $m), which can be used to extract candidate young massive stars with a colour criterion similar to the Wood and Churchwell (1989) criteria for IRAS colours. Indeed, the ISOGAL sources identified as young objects from their IRAS or radio counterparts all have MSX flux ratios FE/FD above 2 (or colours [15]-[21] > 1.55 mag, see Fig. 6). Therefore we will use the criterion $F_E/F_D\geq2$ as an indication of the young nature of ISOGAL-MSX sources in Sect. 4. The interpretation as young objects from the MSX D-E colour alone is somewhat less secure than the Wood and Churchwell (1989) criteria, since objects of this type emit the bulk of their energy at wavelengths longer than 25 $\mu $m. Also, not all sources with $F_E/F_D\geq2$ show other signs of high mass star formation, and the nature of those with moderate [8]-[15] colours ($\la$2 mag) remains unclear (see also Sect. 4.2 and Fig. 10). However, the limit $F_E/F_D\geq2$ is more conservative than the criteria FE>FD>FA and $F_E>2\times F_A$used by Lumsden et al. (2002) to extract massive young objects from the MSX PSC. We will therefore use the criterion $F_E/F_D\geq2$ to confirm the young nature of sources selected on the basis of their red ISOGAL colours.


  \begin{figure}
\par {\includegraphics[width=8cm,clip]{2891f6.eps} }\end{figure} Figure 6: MSX colour-magnitude diagrams of the FC-027 field: [15] vs. [8]-[15] ( top) and [15] vs. [15]-[21] ( bottom). Sources identified with known H II regions, as derived from IRAS or radio continuum measurements, are shown with square symbols, and those associated with LPV or OH/IR stars are shown with $\times $ symbols or diamonds, respectively. The grey crosses show the positions of MSX sources without ISOGAL counterparts within 8''. The dotted line in the bottom panel corresponds to a flux ratio FE/FD=2 (see text).
Open with DEXTER

  
3.5 Near-infrared spectroscopic observations

We have conducted near-infrared (H and K bands) spectroscopic observations with the SOFI instrument on the NTT in July 2000; the observations and their analysis have been extensively described by Schultheis et al. (2003). The purpose of these observations was to derive the nature of sources with various positions in the [15] vs. [7]-[15] colour magnitude diagram from their spectra. Indeed, near infrared spectra contain enough features to distinguish late type evolved stars from young objects: the former show strong molecular bands (mainly CO and H2O), while the latter generally show only atomic lines (hydrogen Br$\gamma$, He I, Fe II, Na I and Ca I, among others), or red featureless continuum, possibly with signatures of shocks (e.g. H2 at 2.12 $\mu $m), characteristic of massive YSOs (Hoare et al. 2004). It is therefore possible to infer the nature of the sources from the equivalent widths of the observed molecular bands, as shown in Schultheis et al. (2003, see in particular their Fig. 5).

The observed sample included 65 sources in the FC-027 field. Only two of them show near-IR spectra consistent with a young stellar nature, or with hot stars in a more evolved stage. These are 174516.2-290315 (labelled A33 in Schultheis et al. 2003), also associated with a radio and an IRAS sources; and 174504.9-291146 (labelled A40), which also has a radio counterpart. The near-IR spectrum of the latter can possibly be interpreted as tracing a Wolf-Rayet star. As can be seen in Fig. 7, these observations tend to confirm that most sources with [7]-[15] colours below or close to 2 mag are evolved stars in the AGB phase.

Five sources with colours [7]-[15] $\ga$ 2.5 mag show near-IR spectra interpreted as AGB (Fig. 7). In addition, the weak ([15] = 4.87, [7]-[15] = 1.74) source 174521.9-291344 has an AGB-type near-IR spectrum, but also has a radio continuum counterpart (Fig. 3). A visual inspection of the images shows that, for these six sources, either a second source or an extended patch of diffuse emission shows up only at 15 $\mu $m. Thus, these are certainly blends between AGB stars, mainly detected in K-band and at 7 $\mu $m, and unrelated diffuse emission that does not appear at wavelengths shorter than 15 $\mu $m, and could possibly trace recent star formation.

  
3.6 Summary

The results of all identifications discussed in the previous sections are shown in the summary [15] vs. [7]-[15] colour magnitude diagram (Fig. 8). It can be seen that sources interpreted as young objects (i.e. with IRAS H II regions or radio continuum counterparts, and the two sources with NIR-spectra of hot stars), occupy a well-defined region in the right-upper part of this diagram, while sources identified with AGB stars are mostly found to the left of [7]-[15] = 2. Only seven OH/IR stars are located in the [7]-[15] > 2, [15] < 4.5 region. Six of them also have MSX counterparts, all with flux ratios FE/FD less than 2, while all identified young objects have MSX flux ratios FE/FD >2. Therefore, FE/FD >2 seems to be a good criterion to confirm the young nature of YSO candidates.

  \begin{figure}
\par {\includegraphics[width=8.5cm,clip]{2891f7.eps} }\end{figure} Figure 7: Colour magnitude diagram for all sources in the FC-027 field, showing the results of near-infrared spectroscopic classification. Sources with spectra of YSO or hot stars are shown with square symbols, while those with AGB or RGB spectra are shown with $\times $ or large plus symbols, respectively. Sources with AGB spectra also associated with LPV or OH/IR stars are indicated with diamond or triangle symbols, respectively.
Open with DEXTER


  \begin{figure}
\par {\includegraphics[width=16.8cm,clip]{2891f8.eps} }\end{figure} Figure 8: Summary colour magnitude diagram for all sources detected at 7 and 15 $\mu $m in the FC-027 field. The symbols show the positions in this diagram of sources identified with known objects (LPVs, OH/IR, IRAS and radio H II regions), and those with near-infrared spectral classification. Sources with MSX associations with flux ratios FE/FD>2 are indicated with open circles, also when detected in E but not in D. The dotted lines show the limits [7]-[15] = 2 and [15] = 4.5 (see text).
Open with DEXTER

Also note that the majority of sources brighter than [15] = 4.5 have known counterparts in either of the young object or AGB classes. Indeed, in the area covered by Glass et al. (2001) observations, the nature of 69 out of 88 sources (almost 80%) with [15] < 4.5 and a 7-15 $\mu $m association can be derived from existing data, or from their red MSX D-E colours for five of them. The fraction of unidentified sources becomes much larger beyond this magnitude limit, where the distinction between young and evolved objects also becomes less clear.

  
4 Extraction of the young stellar population

  
4.1 Selection criteria

According to the summary colour magnitude diagram (Fig. 8 and Sect. 3.6), the simple colour criterion [7]-[15] > 2 seems appropriate to select young objects, especially for the brightest objects ([15] < 4.5 mag). However, a moderate number of evolved stars (e.g. deeply embedded OH/IR stars) also show [7]-[15] colours between 2 and 3 mag. Carbon stars could also have [7]-[15] colours in excess of 2 mag, but it is known that such objects are very rare in the direction of the inner Galactic Bulge  (e.g. Blanco et al. 1984; Glass 1986). Consistently, we found no carbon stars among the $\sim$100 sources in the inner Galactic Bulge for which we obtained near-IR spectra (Schultheis et al. 2003). Planetary nebulae may also show red mid-IR colours, but Jacoby & Van de Steene (2004) estimated their total number in the central 16 deg2 of the Galaxy to be of order 250, so that we expect a very small number of such objects in this $\la$0.1 deg2 field.

Another important parameter of the ISOGAL sources, related to their spatial extension, appears to give a powerful criterion to distinguish massive star forming regions from stars in their late evolution. The uncertainty in the measured magnitude, especially at 15 $\mu $m, noted $\sigma _{15}$ in the ISOGAL PSC (Schuller et al. 2003), is derived from the residuals between the source profile and the analytical PSF used for the extraction. While the values of $\sigma _{15}$ show a larger dispersion for faint sources, and can reach typically 0.1-0.15 mag, they are generally low ($\leq$0.05 mag) for point sources brighter than [15] = 4.5 mag.

The distribution of the $\sigma _{15}$ values is shown versus the [7]-[15] colours for all sources with [15] < 4.5 in Fig. 9. Only one LPV star is found with [7]-[15] > 2 and $\sigma _{15}$ > 0.05, but its mid-infrared magnitudes are affected by a blend with a bright nearby extended source. On the other hand, all sources interpreted as young objects have high values of both [7]-[15] and $\sigma _{15}$. In other terms, stars in the AGB phase appear point-like in the ISO observations, while more extended emission is detected around young sources: these may be small groups or clusters of YSOs, or single objects embedded in diffuse emission. Indeed, ground-based mid-IR images of similar objects selected from the MSX survey show, in most cases, a point-like source embedded in extended diffuse emission, probably arising from PAH (Hoare et al. 2004).


  \begin{figure}
\par {\includegraphics[width=8.5cm,clip]{2891f9.eps} }\end{figure} Figure 9: $\sigma _{15}$ vs. [7]-[15] diagram for all sources detected at both wavelengths and brighter than [15] = 4.5 in FC-027. The sources associated with AGB stars are shown with cross symbols, and those associated with young objects (IRAS H II regions or radio sources) are indicated with square symbols. Sources with MSX associations with FE/FD>2 are indicated with open circles. The dotted line shows the criteria that we have used to extract candidate young objects from the ISOGAL catalogue (Sect. 4.1).
Open with DEXTER

Thus, we can combine a colour and an extension criteria in order to extract the population of ISOGAL sources associated with young objects, with very little contamination by evolved stars. As shown in Fig. 9, the criteria that we have defined are:

 \begin{displaymath}
\left\{ \begin{array}{l}
2 ~ \leq ~ [7]-[15] ~ < ~ 2.5 ~ \mb...
...mbox{no constraint on~} ~ \sigma_{15} .\\
\end{array} \right.
\end{displaymath} (2)

Most sources with an MSX association with a flux ratio FE/FD>2 are found in the same region of this diagram, as shown with open circles in Fig. 9. This gives an additional argument showing that the combination of the proposed colour and spatial extension criteria are appropriate to select candidate young objects.

As compared with the selection criteria that we defined in Felli et al. (2002), the addition of a constraint on the spatial extension at 15 $\mu $m for the sources with moderately red colours (2 $\leq$ [7]-[15] < 3) greatly reduces the contamination of our sample by anomalously red AGB stars (e.g. OH/IR stars with large mass loss, Fig. 9). Also, we consider only sources redder than [7]-[15] = 2 as candidate young objects, while Felli et al. (2002) used [7]-[15] $\geq$ 1.8 to select all candidates from the ISOGAL database. However, both the interstellar extinction and confusion issues due to a high source density are more extreme in the FC-027 field, located near the Galactic Centre, than in the general Galactic Disk environment of most ISOGAL fields. In addition, OH/IR and similar evolved stars are relatively rare outside of the Galactic Bulge: in the $\vert b \vert
< 2^{\circ}$ range, the number density of IRAS sources with a probability of variability greater than 90% (i.e. mostly AGB stars) already drops by a factor $\sim$2 from $l \sim 0^{\circ}$ to $\vert l \vert \sim 20^{\circ}$, while the number density of sources detected at 25, 60 and 100 $\mu $m (mostly young stars) is roughly constant over $\vert l \vert < 80^{\circ}$(see Fig. VII.Ap.7 and VII.Ap.16 in the IRAS Explanatory Supplement). Therefore, we do not expect many evolved sources in the Felli et al. (2002) sample extracted with somewhat looser criteria. On the other hand, we may have missed a few young objects with [7]-[15] colours below 2 mag in the present study.

We also want to stress that our criterion on $\sigma _{15}$ selects sources marginally resolved with ISO, i.e. with angular sizes of order 30''. This roughly corresponds to a 1 pc linear size at the distance of the Galactic Centre. Thus, our combination of colour, magnitude and spatial extension criteria is well adapted to select young objects in the direction of the Galactic Bulge, where most sources are expected to be at similar distances (see also Sect. 4.4 below), but it would be much less efficient in sight lines where distances to sources are spread over a wide range.

  
4.2 Extraction of all YSO candidates

Using the proposed selection criteria (see Eq. (2)), 33 candidate young objects with [15] < 4.5 mag (or F15 > 280 mJy) were selected, including 13 with IRAS and/or radio continuum counterparts. In addition, 15 sources are found with [15] < 4.5 mag, but with no association at 7 $\mu $m in the ISOGAL PSC (which has a 50%-completeness level around 8.0 mag at 7 $\mu $m, see Sect. 2). One of them (174529.1-290413, with [15] = 4.41 and $\sigma _{15}$ = 0.13) is associated with an LPV star; the lack of magnitude at 7 $\mu $m and the high $\sigma _{15}$value can be explained by confusion issues. The other 14 sources with [15] < 4.5 all have $\sigma_{15} \geq 0.07$, and can also be considered as candidate young objects.

Among these 46 candidates, 30 (or 65%) are also associated with MSX sources. All of them are detected in the A, D and E bands, and all but two have flux ratios FE/FD>2, as can be seen in Fig. 10. On the other hand, only four LPV stars and one OH/IR star have FE/FDabove 2. This colour-colour diagram therefore seems to be an efficient means for distinguishing young objects from evolved stars, whereas the mid-IR selection criteria used by Lumsden et al. (2002) do not reject evolved stars as efficiently, as can be seen in Fig. 10 - however, they used additional criteria based on 2MASS near-IR magnitudes to eliminate many of the evolved sources.


  \begin{figure}
\par {\includegraphics[width=7.5cm,clip]{2891f10.eps} }\end{figure} Figure 10: MSX colour colour diagram for all sources associated with ISOGAL in FC-027. The YSO candidates that we selected from this field are indicated with square and triangle symbols. Sources associated with LPV or OH/IR stars are shown with crosses and diamonds. The dotted horizontal line corresponds to FE/FD=2, and the dashed lines show the criteria used by Lumsden et al. (2002).
Open with DEXTER

The most relevant infrared properties of our 46 YSO candidates are given in Tables A.2 and A.3 in Appendix A. Tables with all available magnitudes (ISO, DENIS, 2MASS, MSX and IRAS when relevant) can be retrieved electronically through CDS. These tables also contain aperture magnitudes at 15 $\mu $m. These were computed by integrating a number of pixels in an inner radius, and subtracting the background estimated as the median in an outer annulus; the radii that we used were adapted for every individual source, in order to avoid contaminations by nearby sources.

  
4.3 Comments on the cluster around LS95 7

A compact group of nine ISOGAL point sources (Fig. 11) is located in the vicinity of the radio source No. 7 detected at 1.6 GHz by Liszt & Spiker (1995), who reported an extended size of $31''\times48''$ (Sect. 3.2). Sources extracted in a 30'' radius around 17$\rm ^h$45$\rm ^m$05$\rm ^s$, $-29^{\circ }16'50''$ are reported in Table 2. At least two of them, sources C and H, are probably young objects, since they have [7]-[15] colours above 3.5 mag, and signs of spatial extension as indicated by value of $\sigma _{15}$ above 0.1 mag. These two sources are in the samples extracted with our criteria discussed in Sect. 4.2.


  \begin{figure}
\par {\includegraphics[width=8.2cm,clip]{2891f11.eps} }\end{figure} Figure 11: Contour plots of the ISOCAM LW5 (7 $\mu $m, black contours) and LW9 (15 $\mu $m, white contours) around ($\alpha $, $\delta $) = ( $\rm 17^h45^m05^s$, $-29^{\circ }16'50''$). The LW5 contours correspond to flux densities of 40, 60, 80, 110, 150, 200 and 300 mJy/pixel, and the LW9 contours correspond to 60, 100, 140, 180, 220 and 260 mJy/pixel. The crosses and letters indicate the positions of extracted ISOGAL sources, and details are given in Table 2. The ellipse drawn with thick dashed line shows the position and size (FWHM) of the radio source LS95 7. The very bright source at the limit of the field to the north is associated with LS95 8.
Open with DEXTER

Table 2: Coordinates, 7 and 15 $\mu $m magnitudes, extension at 15 $\mu $m $\sigma _{15}$ and 7-15 $\mu $m association quality flag for the sources located within the cluster. The Id. column gives the identification letters, as they appear in Fig. 11.

Diffuse emission is also visible around these sources at 15 $\mu $m, suggesting that this could be a group of young objects still in their parent cocoon. We can also note that if this group is at the distance of the Galactic Centre, the projected distance between sources C and H would be of order 1 pc, which could correspond to a very young OB-association. Complementary observations at higher spatial resolution, with a sensitivity at least comparable to ISOGAL, would greatly help in understanding the nature of this object.

  
4.4 Distance to the sources

Most radio sources that we considered in the present paper were only observed in the continuum, so that no radial velocity data are available. Only for the source LS95 7 that we discussed in the previous section, Liszt & Spiker (1995) report a radial velocity $V_{\rm LSR}=-44.1$ km s-1, as derived from observations in the H70$\alpha $ radio recombination line. Thus, this source is probably associated with the Sgr C complex, in the vicinity of the Galactic Centre.

Using the available data, we cannot exclude that some of the candidates that we have selected are foreground. However, given the high over-density of such sources in the inner Galactic Bulge as compared to the rest of the Galactic disk, we expect a very small fraction of foreground sources: the Besançon model (Robin et al. 2003) predicts $\sim$8% contamination by foreground objects. It is still possible that a few of our candidates are much closer than the Galactic Centre, especially among the brightest sources. Additional observing programmes are ongoing to further constrain the distance to these sources.

4.5 Summary


  \begin{figure}
\par {\includegraphics[width=15.4cm,clip]{2891f12.eps} }\end{figure} Figure 12: Distribution in galactic coordinates of our 46 candidates. All 2043 sources from the ISOGAL catalogue are represented with small dots. The 32 candidates with 7 and 15 $\mu $m detections are shown with plus signs, and the 14 candidates only seen at 15 $\mu $m with crosses. Sources with IRAS counterparts are indicated by triangles and those with radio continuum by squares (see also the IRAS uncertainty ellipses in Fig. 4). The sources associated with MSX point sources are shown with small diamonds. The contours show the radio continuum emission at 1.6 GHz (data were kindly provided by Dr. H. S. Liszt), with levels at 0., 0.01, 0.02, 0.04, 0.06, 0.08, 0.1, 0.15, 0.2 and 0.25 Jy/beam.
Open with DEXTER

We have been able to extract a sample of 46 candidate massive young objects brighter than 4.5 mag at 15 $\mu $m (or $F_{\nu} > 280$ mJy) from a 0.09 deg2 area at about 40 pc in projection from the Galactic Centre, using mid-infrared photometry from the ISOGAL survey. We have adopted a conservative approach, since this survey goes down to magnitude 7.0 at 15 $\mu $m with good completeness (Sect. 2); the distinction between young and evolved stars becomes more subtle when going to fainter levels than this. A source located at 8.5 kpc with magnitude [15] = 4.5 would have a luminosity around 104 $L_{\odot}$ according to Felli et al. (2002Sect. 3.2), consistent with a single B0 star.

Only 13 of these 46 sources show other signatures of young massive stars, such as IRAS or radio continuum counterparts (including one not seen by MSX), and 16 do not appear in the MSX PSC, while most of the others display MSX colours typical of young stars. The distribution of these YSO candidates is shown in galactic coordinates in Fig. 12. It appears that most of them are located in the $-0.1^{\circ}\leq b
\leq 0^{\circ}$ range, and that their distribution is tilted with respect to the Galactic Plane, going from $b \sim
0^{\circ}$ at $l \sim -0.1^{\circ}$ to $b \sim -0.1^{\circ}$at $l \sim -0.4^{\circ}$. In the mid-infrared, diffuse emission probably arising from PAHs and small dust grains is seen in the same latitude range (Fig. 1). Moreover, strong radio continuum emission at 1.6 GHz is observed at the positions of most candidates. Some of them are coincident with peaks of radio emission, and others are located in regions of extended radio emission, including eight candidates without known counterparts located in the extended component around Sgr A. This adds to the confidence that we are interpreting these sources correctly as young massive objects.

5 Conclusion

We have shown in this paper that the bright end of the stellar populations detected by ISOGAL can be well characterised, using already existing data and complementary mid-infrared observations. A sample containing 182 sources brighter than 280 mJy at 15 $\mu $m, complete over a 0.09 deg2 area close to the Galactic Centre, has been associated with several catalogues of evolved stars on one hand, and with high-mass young objects previously detected by IRAS or in the radio continuum on the other. With the addition of infrared photometry at 15 and 21 $\mu $m from MSX, and spectroscopic near-infrared data for some sources, a majority (61%) could be identified with evolved stars or with young objects.

The mid-infrared magnitudes and colours and the spatial extension of the 15 $\mu $m emission of the young objects are sufficiently peculiar to provide robust criteria for selecting good candidate young massive objects from the ISOGAL data. Using a colour criterion [7]-[15] $\geq$ 2 mag, combined with a minimum spatial extension indicated by $\sigma_{15}>0.1$ mag for the least red sources ([7]-[15] between 2 and 2.5 mag), we could extract 46 candidates from the analysed field. Only 13 of them were previously reported as ultra-compact H II regions. About 65% are also detected by MSX, and they also have red colours in the longer 15 to 21 $\mu $m wavelengths. Their positions in the field are almost always coincident with diffuse extended infrared emission, with peaks or more extended emission in the radio continuum in most cases. They are mainly distributed in the $-0.1^{\circ}~ \leq ~ b ~ \leq ~ 0^{\circ}$ range, or very close to the Galactic Plane, where one expects to find the largest incidence of star formation. All this tends to confirm their interpretation as young objects, embedded in vast structures of dust and gas heated by recently formed stars.

The exact nature of our objects remains unclear. Following Felli et al. (2002), their bolometric luminosities can be roughly derived from their flux densities at 15 $\mu $m, and one finds values in the range 104-106 $L_{\odot}$, assuming that they are at the distance of the Galactic Centre. Such luminosities could come from single O or early B stars. However, most of them show spatial extension of order 1 pc (as resolved by ISO) and only very few are associated with known ultra-compact H II regions, suggesting that there is no very early type star present in most cases. Thus, these objects could be small groups or clusters of young stars of intermediate masses. To be detected with such red colours in the mid-infrared, they would have to be in a protostellar phase, possibly in the compact molecular core stage, prior to the formation of a compact H II region (e.g. Brand et al. 2001), or an evolutionary stage somehow similar to class I for low-mass YSOs. Some candidates may also be compact H II regions falling below the sensitivity of current radio continuum surveys, thus corresponding to the next stage in the early evolution of massive stars.

Several follow-up projects, including spectro-imagery with the Spitzer Space Telescope, have been initiated to better understand the nature of these objects. Moreover, the ISOGAL data cover most parts of the inner Galactic bulge, in the ${\pm} 1.5^{\circ}$longitude range. The criteria that we have presented here can be used to extract the young population from the other ISOGAL fields, to get a more complete picture of the recent star formation in the central $\sim$400 pc of the Galaxy; this will be the subject of a second paper.

Acknowledgements
We are very grateful to Dr. H. S. Liszt for providing us with his VLA 1.6 GHz data of the Galactic Centre region. This research made use of data products from the Midcourse Space Experiment. Processing of the data was funded by the Ballistic Missile Defense Organization with additional support from NASA Office of Space Science. We gratefully thank the MSX team and in particular Dr. D. Mizuno for making us parts of the MSX catalogue version 2.3 available prior to its publication. MS is supported by the APART programme of the Austrian Academy of Science.

   
Appendix A: Complete list of YSO candidates

Table A.1 shows two examples of the full data for our 46 candidate young objects. This table contains photometry and quality flags from the 2MASS, DENIS, ISOGAL, MSX and IRAS catalogues. The separations between the ISOGAL objects and the other sources are also given, as well as the separations between the 7 $\mu $m and 15 $\mu $m sources. The ISOGAL 15 $\mu $m flux densities have been determined by PSF-fitting (Col. F15) and by aperture photometry (Col. ap15). The complete version of this table is available through VizieR at the CDS.

Table A.1: 2MASS and DENIS magnitudes, ISOGAL, MSX and IRAS fluxes for two among our 46 candidate young objects. The full table is available through Vizier at the CDS1.

Tables A.2 and A.3 show only the ISOGAL and DENIS data for the candidates detected at 7 and 15 $\mu $m, and the ones only detected at 15 $\mu $m, respectively. These tables are available in the online version of the present paper.

References

 

  
Online Material

Table A.2: Main characteristics of the 32 YSO candidates detected at 7 and 15 $\mu $m: J and K DENIS magnitudes, ISO-DENIS association quality flag $q_{\rm ID}$and separation; magnitudes, uncertainties and quality flags at 7 and 15 $\mu $m, 7-15 $\mu $m association quality flag $q_{\rm II}$ and separation; aperture magnitude at 15 $\mu $m; MSX: this column contains a "+'' if the source is associated (within 8'') with an MSX source with FE/FD>2, and a "X'' if there is an association, but with FE/FD<2.

Table A.3: Main characteristics of the 14 YSO candidates detected only at 15 $\mu $m. Columns are the same as in Table A.2, except for the 7 $\mu $m and 7-15 $\mu $m association data.



Copyright ESO 2006