A&A 433, 117-125 (2005)
DOI: 10.1051/0004-6361:20041748

L- and M-band imaging observations of the Galactic Center region[*]

T. Viehmann1 - A. Eckart1 - R. Schödel1 - J. Moultaka1 - C. Straubmeier1 - J.-U. Pott1,2


1 - I. Physikalisches Institut, Universität zu Köln, Zülpicher Str. 77, 50937 Köln, Germany
2 - ESO, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany

Received 28 July 2004 / Accepted 16 November 2004

Abstract
In this paper we present near-infrared H-, K-, L- and M-band photometry of the Galactic Center from images obtained at the ESO VLT in May and August 2002, using the NAOS/CONICA (H and K) and the ISAAC (L and M) instruments[*]. The large field of view (70 $^{\prime \prime }$ $\times $ 70 $^{\prime \prime }$) of the ISAAC instrument and the large number of sources identified (L-M data for a total of 541 sources) allows us to investigate their colors, infrared excesses and the extended dust emission. Our new L-band magnitude calibration reveals an important offset to the traditionally used "standard'' calibrations, which we attribute to the use of the variable star IRS 7 as a flux calibrator. Together with new results on the extinction towards the Galactic Center (Raab 2000; Scoville et al. 2003), our magnitude calibration results in stellar color properties expected from standard stars and removes any necessity to modify the K-band extinction. The large number of sources for which we have obtained L-M colors allows us to measure the M-band extinction to AM = (0.056 $\pm$ 0.006)$A_{\rm V}$, i.e. $A_M\approx A_L$, a considerably higher value than what has so far been assumed. L-M color data has not been investigated previously, due to lack of useful M-band data. We find that this color is a useful diagnostic tool for the preliminary identification of stellar types, since hot and cool stars show a fairly clear L-M color separation, with average L-M colors of 0.46 $\pm$ 0.17 for emission-line stars and -0.05 $\pm$ 0.27 for cool red giants/AGB stars. This is especially important if visual colors are not available, as in the Galactic Center.
For one of the most prominent dust embedded sources, IRS 3, we find extended L- and M-band continuum emission with a characteristic bow-shock shape. An explanation for this unusual appearance is that IRS 3 consists of a massive, hot, young mass-losing star surrounded by an optically thick, extended dust shell, which is pushed northwest by wind from the direction of the IRS 16 cluster and Sgr A*.

Key words: Galaxy center - infrared: stars - techniques: photometric - stars: imaging

   
1 Introduction

The Galactic Center stellar cluster shows some intriguing characteristics: it is extremely dense, with an unusual observed stellar population consisting mainly (80% of the $K\la14$ stars according to Ott et al. 1999) of late-type red giants, many of which are suspected to lie on the asymptotic giant branch (AGB), as well as young massive stars with energetic winds. Spectroscopic measurements (e.g. Najarro et al. 1997; Krabbe et al. 1995) allow identification of these two components, the late-type red giants and supergiants showing strong 2.3 $\mu$m CO bandhead absorption and the massive, hot and windy stars ("He stars'') exhibiting He and H emission lines. The emission line stars appear to dominate energetically in the central few arcseconds, where the bright IRS 16 cluster is located. These stars are generally classified as Ofpe/WN9, although some of them might be luminous blue variables (LBV) and a few show characteristics of Wolf-Rayet stars. A third, less numerous component of the Galactic Center stellar cluster consists of luminous, extended objects with steep, red and featureless (K-band) spectra and a strong infrared excess. These sources were considered as potential YSO candidates (Clénet et al. 2001; Ott et al. 1999; Krabbe et al. 1995). Recently, however, they have been successfully described as bow-shock sources (Rigaut et al. 2003; Eckart 2004; Tanner et al. 2002,2003), with the implication that the stars powering these luminous infrared sources are windy stars of the types mentioned above.

We present multi-band (H, K, L and M) photometry of the Galactic Center, obtained with the ESO VLT, in order to further investigate the properties of the Galactic Center stellar cluster. The large number of sources available for color analysis allows for improved statistics, resulting in a new L-band calibration which eliminates anomalous color effects found in other surveys of the Galactic Center stellar cluster (e.g. Simons & Becklin 1996; Clénet et al. 2001). The large field of view (70 $^{\prime \prime }$ $\times $ 70 $^{\prime \prime }$) of the ISAAC instrument allows us to compare the immediate Galactic Center cluster population with the stellar population further out, with separations from SgrA* of up to about half an arcminute. This comparison results in a confirmation of the new extinction values of Scoville et al. (2003). Previous (multi-band) imaging studies of the Galactic Center, besides being limited in their field of view to e.g. 13 $^{\prime \prime }$ $\times $ 13 $^{\prime \prime }$ (Clénet et al. 2001; Davidge et al. 1997), appear to have neglected the M-band (4.66 $\mu$m), which shows some interesting properties, which we discuss in this paper.

   
2 Observations and data reduction

2.1 ISAAC observations

On May 23rd-30th 2002, L- and M-band (3.78 $\mu$m and 4.66 $\mu$m) imaging observations of the Galactic Center were obtained with the ISAAC infrared imager at the ESO VLT (Unit Telescope 1, Antu), as part of a joint IR/X-ray monitoring program of Sgr A* (Baganoff et al. 2003; Eckart et al. 2003b). The seeing-limited images were taken in ISAAC's higher resolution mode, with a pixel scale of 0.071 $^{\prime \prime }$ and a field of view of 72.6 $^{\prime \prime }$ $\times $ 72.6 $^{\prime \prime }$ (see Cuby et al. 2002). Seeing conditions varied during the campaign, resulting in angular resolutions between 0.4 $^{\prime \prime }$ and 1.2 $^{\prime \prime }$. The observations were made in batches of 25 to 40 pairs of frames (with chopping and nodding), with an integration time of 0.99 s per frame.

2.2 NAOS/CONICA observations

During the science verification phase of the commissioning and observatory preparation of the NAOS/CONICA adaptive optics assisted imager/spectrometer at the ESO VLT (UT4, Yepun), H- and $K_{\rm s}$-band (1.6 $\mu$m and 2.1 $\mu$m) images of the Galactic Center were obtained on 30 August 2002, with a FWHM of 43 and 56 mas, respectively, and a pixel scale of 0.0132 $^{\prime \prime }$/pixel (Schödel et al. 2002,2003). The total integration times were 1500 s and 1200 s, respectively. Since the $K-K_{\rm s}$ color difference is an order of magnitude smaller than the photometric accuracy that we reach, no further distinction between K (2.2 $\mu$m) and $K_{\rm s}$ will be made in this paper. Since the purpose of the NAOS/CONICA images in this context was to gain additional color information for analysing the seeing-limited ISAAC L- and M-band images, the pixel scale of the NAOS/CONICA images was reduced to the corresponding 0.071 $^{\prime \prime }$/pixel via re-binning. Further data reduction was then performed using the same methods as with the ISAAC images.
  \begin{figure}
\par\includegraphics[width=15.6cm,clip]{1748fig1.eps}
\end{figure} Figure 1: ISAAC M-band final image of the Galactic Center. The field of view is 72.6 $^{\prime \prime }$ $\times $ 72.6 $^{\prime \prime }$.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=15.6cm,clip]{1748fig2.eps}
\end{figure} Figure 2: Part of the ISAAC L-band final image of the Galactic Center showing most of the named sources.
Open with DEXTER

2.3 Data reduction

All image processing was performed using routines from the IRAF software package. The data reduction procedures employed were standard: after flat field correction, the paired images were subtracted from each other and subsequently shifted to a common reference frame. All such subtracted images from each batch were then combined by calculating the median, which also effectively removed cosmic rays and bad pixels. The final images were then created by coaddition of the median images with the best effective seeing, i.e. better than 0.6 $^{\prime \prime }$, resulting in total integration times of 229.7 s for the L-band image and 162.4 s for the M-band image. The final M-band image is shown in Fig. 1. Finally, PSF-fitting photometry was performed on all four final images using the IRAF/DAOPHOT software package.

   
2.4 Photometry

The observing conditions for both the ISAAC and the NAOS/CONICA observations were generally photometric and are adequate for relative photometry. For the L- and M-band case, we linked the calibration of the IRAF/DAOPHOT PSF-fitting results to observed flux density reference stars. For the H- and K-band, we relied on data taken from the literature (Depoy 1990; Clénet et al. 2001).

L- and M-band magnitude calibration
For calibration purposes, the standard stars HD 130163 (mL=6.856, $m_L-m_M\approx0$, van der Bliek et al. 1996) and HR 5249 (mM=4.53, van der Bliek et al. 1996) were imaged together with science observations on May 28th 2002, using identical techniques. The airmass difference was small (maximum 0.2), resulting in corrections smaller than 0.04 magnitudes, which are insignificant given the uncertainties in the photometry (see Sect. 3.1). These observations were used to calibrate the final L- and M-band images, using aperture photometry on approximately 30 isolated stars.

K-band magnitude calibration
Reference K-band magnitudes for some of the bright IRS sources were taken from Clénet et al. (2001), who fixed their zero-point with the photometry of Ott et al. (1999). The differences between the DAOPHOT magnitudes and the values of Clénet et al. (2001) were averaged, resulting in a calibration value which was subtracted from the DAOPHOT magnitudes. Saturated sources were excluded from the determination of the calibration magnitudes.

H-band magnitude calibration
The main problem encountered in calibrating the H-band data was saturation. All the bright H-band sources (for which magnitudes are available) are saturated in the NAOS/CONICA image, whereas no magnitudes for the weaker unsaturated sources could be found in the literature.
The solution employed here was secondary calibration using an H-band mosaic made from images taken with the SHARP speckle imager in 1992 (Eckart et al. 1992). This image was flux calibrated using data for the brighter IRS sources (Depoy 1990). Subsequently, the NAOS/CONICA image was calibrated with the flux densities from weaker, isolated (i.e. clearly separated from their neighbors) sources in the SHARP mosaic image.

The total uncertainty of the photometry results is estimated as 0.25 mag in the H- and K-band and 0.15 mag in the L- and M-band. The difference is due to the different calibration procedures and the different shape of the respective PSFs of the ISAAC and NAOS/CONICA images, resulting in less accurate fit results for the latter. The limiting magnitudes of this analysis are 12.73 (M-band), 14.47 (L-band), 17.66 (K-band) and 20.34 (H-band).

2.5 Additional data

The field of view of the NAOS/CONICA images (approximately 17.8 $^{\prime \prime }$ $\times $ 17.8 $^{\prime \prime }$) is much smaller than the ISAAC field of view. In order to check the properties of sources outside this region, additional H- and K-band images were analysed.

ISAAC K-band image
During the 1999 ISAAC spectroscopy campaign (Eckart et al. 1999), a K-band image of the Galactic Center was obtained, using ISAAC's lower resolution mode (pixel scale 0.148 $^{\prime \prime }$, field of view approximately 150 $^{\prime \prime }$ $\times $ 150 $^{\prime \prime }$, see Cuby et al. 2002).

Gemini data
For H-band, we also used publicly released Gemini data from the 2000 Galactic Center Demonstration Science Data Set (using the Quirc camera and the Hokupa'a adaptive optics system, see Cotera & Rigaut 2001).
In order to limit the amount of work required, photometry on these datasets was performed only for the brighter sources in the L-band image (mL<9.5) and an exemplary selection of approximately 15 fainter sources in order to verify that the faint stars do not show systematically different color properties than the bright stars. The images were calibrated relative to our NAOS/CONICA data.

   
3 Results

   
3.1 Photometry

Table 1 summarizes the photometry of those sources within the NAOS/CONICA field of view (17.8 $^{\prime \prime }$ $\times $ 17.8 $^{\prime \prime }$) that are detectable in at least one of the ISAAC images. Table 2 summarizes the photometry of the sources outside the NAOS/CONICA field of view, for which therefore only L- and M-band data are available.

Notes:

   
3.2 New L-band calibration

Since our L-band calibration introduces an offset of approximately 0.5 mag to the frequently used calibration of Blum et al. (1996), some comments are necessary.

Calibrating the L-band image with reference magnitudes taken from Clénet et al. (2001), as for the NAOS/CONICA K-band image (see Sect. 2.4), resulted in $(K-L)_{{\rm obs}}$ colors which were at least 0.3 mag too blue. An examination of the color-color-diagrams showed this effect to be systematic, since the whole population was blueshifted in K-L, whereas the H-K colors lie within the expected range for stellar sources. It therefore seems unlikely that the K-band calibration is wrong, especially considering that it is based on stars shown to be non-variable (Clénet et al. 2001; Ott et al. 1999). Directly calibrating the L-band image with the observed calibration star instead of adjusting the zero-point with already published photometry does not show this effect. We note that Clénet et al. (2001) also found unusually blue K-L values, which they attributed to a lower K-band extinction. If the extinction law remains unchanged, this implicitly suggests a lower extinction $A_{\rm V}$, which our observations can confirm (see Sect. 3.3).

The L-band data and thus the calibration given by Blum et al. (1996) are taken from Depoy & Sharp (1991). Their calibration (which is not discussed in detail in their 1991 publication) was fixed with the star IRS 7, for which Tamura et al. (1996); Blum et al. (1996) and Ott et al. (1999) have shown that it is variable with an amplitude of of approximately 1 mag in K-band. It seems that the offset to the other frequently applied calibration by Simons & Becklin (1996) is also due to their use of the variable source IRS 7 as flux calibrator. Our new L-band calibration, however, provides a clear improvement to this point, especially since it does not rely on any (non-)variability assumptions for the observed sources and does not result in unexpected K-L color effects. We finally note that a higher extinction than the assumed $A_{\rm V}=25$ mag would result in even bluer K-L colors, requiring an even larger offset to the L-band calibrations mentioned above.

   
3.3 Colors and extinction

For de-reddening purposes, the extinction law of Rieke & Lebofsky (1985) was adopted, i.e. $A_H=0.175~A_{\rm V}$, $A_K=0.112~A_{\rm V}$, $A_L=0.058~A_{\rm V}$ and (initially) $A_M=0.023~A_{\rm V}$. Since recent work by Scoville et al. (2003) indicates that the extinction toward the Galactic Center exhibits a broad minimum over the area of the IRS 16 Cluster and the Northern Arm, their value of $A_{\rm V}=25$ mag was assumed, instead of the traditional average value of $A_{\rm V}=27$ or even 30 mag (Rieke 1987). This choice is supported by the colors of the stars observed, especially K-L (see Sect. 3.2), since they would otherwise be too blue.

In the HKL color-color diagram (Fig. 5), the majority of sources form a cluster in the area populated by stars (for intrinsic colors see Koornneef 1983), with a fairly sharp boundary towards bluer colors. This is satisfactory, since the majority of sources exhibit colors that lie in the possible color range for stellar sources. Only very few objects show bluer colors, which can be seen as exceptions. Quite a few objects appear to be strongly reddened, especially within the inner area covered by the NAOS/CONICA images. This is not surprising, since the dust and gas of the "mini-spiral'' passes through this region.

The conclusion to be drawn here is that the newer findings concerning extinction (Scoville et al. 2003), together with our new L-band calibration, provide a consistent picture without requiring modifications of the extinction law at this point.

3.3.1 Comparison of outer and inner field
With the exception of some highly reddened sources (e.g. IRS 1W, IRS 3, IRS 10, IRS 21) found almost exclusively within the inner field (i.e. the 18 $^{\prime \prime }$ $\times $ 18 $^{\prime \prime }$ of the NAOS/CONICA field of view), the stars from the inner and outer regions cluster in the same area of the two-color diagrams, with the clustering much more pronounced for the stars from the outer region. The highly reddened sources from the inner field are apparently obscured/reddened by gas and dust from the mini-spiral (either directly or by a bow-shock type interaction).

The H-K and K-L colors obtained from our data, especially if one excludes unusual stars such as Wolf-Rayet stars etc., confirm the results of Scoville et al. (2003) for the lower extinction (we assumed $A_{\rm V}=25$ mag) towards the IRS 16 cluster and the Northern Arm of the mini-spiral, since the colors - after dereddening - agree with the intrinsic colors of stars given by Koornneef (1983). Any higher extinction would lead to significantly bluer intrinsic colors, which cannot be explained in a straightforward manner. Our data indicates that this relatively low extinction applies to the complete region covered by the (high-resolution) ISAAC images, i.e. the inner 70 $^{\prime \prime }$ $\times $ 70 $^{\prime \prime }$. This can be inferred from the two-color diagrams (Figs. 5 and 6), since the populations of the inner and outer region cluster in the same area of these diagrams. If the extinction were to increase significantly within the ISAAC high resolution field of view, the "outer'' population (or at least a significant part of it) should show systematically redder colors than the main part of the "inner'' population, which is clearly not the case. The higher extinctions of 35 to 50 mag reported by Scoville et al. (2003) thus apply to areas outside the inner 70 arcsec, which is consistent with the location of the circum-nuclear molecular ring, which has a sharp inner edge (Latvakoski et al. 1999).

3.3.2 Extinction in M-band
Our M-band calibration, using two calibration stars of significantly different magnitude (see Sect. 2.4), provides rather surprising results: the stars appear too blue by about 0.7 mag in L-M, if one applies the reddening law of Rieke & Lebofsky (1985), i.e. $A_M=0.023~A_{\rm V}$. However, the ISO/SWS Galactic Center extinction measurements of Lutz et al. (1996) and Lutz (1999), while generally in agreement with the extinction law of Rieke & Lebofsky (1985), permit considerably higher (or lower) extinctions in the 4-5 $\mu$m regime within their uncertainties. This is especially true if one considers that their extinction values increase again longward of 5 $\mu$m (see also Raab 2000). Visual inspection of the distribution of L-M colors in the color-color diagram (Fig. 6) suggests that $A_M\approx A_L$, since the bulk of the L-M colors observed then agree with the expected L-M colors for stellar sources (Koornneef 1983).

Our extensive L-M color data available for a total of 541 sources allows us to quantify the M-band extinction more precisely, assuming that the theoretical intrinsic colors given by Koornneef (1983) are valid. From the mean observed and theoretical intrinsic L-M colors averaged over all types of stars we can calculate the M-band extinction via:

                                $\displaystyle \frac{A_M}{A_{\rm V}}$ = $\displaystyle \frac{1}{A_{\rm V}}(\langle L-M\rangle_{{\rm Koornneef}}-\langle L-M\rangle_{{\rm obs}}+A_L)$  
  = $\displaystyle 0.056\pm0.006.$  

Here we quote a 3$\sigma$ uncertainty. This result is plotted in Fig. 7 and agrees very well with the above estimate of $A_M\approx A_L$, with $A_L=0.058~A_{\rm V}$ following Rieke & Lebofsky (1985). The possible luminosity-induced bias (observed stars are bright) is <0.003 and therefore not significant given the uncertainties.

Assuming a Gaussian distribution of the L-M values for simplicity, we find a corresponding standard deviation of $\sigma_{L-M}=0.51$, which allows the calculation of the contribution $\sigma_{{\rm ext}}$ of screening and source intrinsic effects to the L-M distribution. The uncertainties of the distribution are given by

\begin{displaymath}\sigma_{L-M}^2=\sigma_L^2+\sigma_M^2+\sigma_{{\rm Koornneef}}^2+\sigma_{{\rm ext}}^2\
,\end{displaymath}

where $\sigma_{{\rm Koornneef}}=0.16$ represents the scattering of theoretical intrinsic colors over the whole range of stellar types as listed by Koornneef (1983). Here $\sigma_L$ and $\sigma_M$ are the uncertainties in the L- and M-band magnitudes. We find $\sigma_{{\rm ext}}=0.44$, which implies that the main contribution to the width of the L-M color distribution is due to screening effects and source-intrinsic reddening (i.e. local dust concentrations etc.). The main caveat here is that the L-M color distribution is not Gaussian (there are two principal stellar types with different intrinsic colors, i.e. red giants and young, massive emission-line stars). However, the width of the distribution due to different stellar type components is considerably smaller than the scattering observed. Nevertheless, this analysis is useful as a preliminary estimate, and further M-band studies of the Galactic Center are necessary to clarify the situation.
  \begin{figure}
\par\includegraphics[width=16.4cm,clip]{1748fig3.eps}
\end{figure} Figure 3: H-, K- and M-band view of IRS 3 (encircled in the H-band view). The bright source to the northeast ( upper left) is the supergiant IRS 7. A histogram equalization color scale is used, while the contour lines follow a logarithmic scale in order to show the IRS 3 source structure as clearly as possible.
Open with DEXTER

3.3.3 Colors of He-stars and red giants
A comparison of the L-M colors of the named sources of known type (i.e. He-stars or cool red giants/AGB) shows a fairly clear color separation, with the cool stars significantly bluer (see Table 1 and Fig. 6). This result is in general agreement with the colors given by Koornneef (1983): hot stars have on average an $(L-M)_{{\rm intrinsic}}$ of $\approx$-0.05, while cool giants have $(L-M)_{{\rm intrinsic}}$ values of typically -0.2 to -0.3. In our sample, the average observed $(L-M)_{{\rm obs}}$ values are 0.46 $\pm$ 0.17 for the emission-line stars and -0.05 $\pm$ 0.27 for the cool stars (ten of each, see Table 1), with the standard deviation as uncertainty. Therefore, the L-M color, and more particularly the H-K versus L-M color-color diagram is a reasonably good diagnostic tool for preliminary identification of these stellar types. Applying this L-M criterion to our sample, we find that the emission-line stars show a clear concentration in the inner region, with only about ten emission-line-star candidates in the area outside the 18 $^{\prime \prime }$ $\times $ 18 $^{\prime \prime }$ NAOS/CONICA field of view (containing 166 objects with full photometry data), as opposed to at least forty such objects in the inner region (197 objects in total).


  \begin{figure}
\par\includegraphics[width=8.4cm,clip]{1748fig4.eps}
\end{figure} Figure 4: L-band view of IRS 3 with NAOS/CONICA (Eckart 2004). The circular structures surrounding the bright sources are due to the adaptive optics. Note the asymmetry in the lower contour lines producing a bow-shock-like effect at the southern edge of the extended emission.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=16.4cm,clip]{1748fig5.eps}
\end{figure} Figure 5: Two-color diagram showing $(K-L)_{{\rm obs}}$ versus $(H-K)_{{\rm obs}}$. The triangles represent sources from the inner 18 $^{\prime \prime }$ $\times $ 18 $^{\prime \prime }$, the squares represent brighter sources (mL<9.5) and the diamonds a selection of fainter sources from outside this inner region.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=16.4cm,clip]{1748fig6.eps}
\end{figure} Figure 6: Two-color diagram showing $(L-M)_{{\rm obs}}$ versus $(K-L)_{{\rm obs}}$. The triangles represent sources from the inner 18 $^{\prime \prime }$ $\times $ 18 $^{\prime \prime }$, the squares represent brighter sources (mL<9.5) and the diamonds a selection of fainter sources from outside this inner region. The marked areas fully contain the spectroscopically identified He stars and red giants / AGB stars from Table 1.
Open with DEXTER

K-L colors are much less useful in this context, since the K-L color separation is comparatively low. Naturally, spectroscopic evidence is required for a definite classification.

An additional comment can be made on the dereddened (i.e. intrinsic) $(K-L)_{{\rm intrinsic}}$ and $(L-M)_{{\rm intrinsic}}$ colors of the He-stars. These appear to be slightly but systematically redder ($\approx$0.2 mag) than expected for hot stars, whose $(K-L)_{{\rm intrinsic}}$ and $(L-M)_{{\rm intrinsic}}$ colors are approximately 0 to -0.1. This indicates the presence of an infrared excess. A possible explanation for such an excess could be that the He-stars all show signs of an interaction with the dusty Galactic Center ISM. This excess could be of a similar nature but much less pronounced than the bow shock sources along the mini-spiral (Rigaut et al. 2003; Eckart 2004; Tanner et al. 2002; Clénet et al. 2004; Tanner et al. 2003). Our result is in general agreement with the findings of Davidge et al. (1997), who indicate that the K-L colors of the He-stars are clearly redder than those of known emission line stars in the LMC.

   
3.4 IRS 3

One of the most intriguing mid-infrared sources in the Galactic Center region is IRS 3. In the L- and M-band, this is the most extended source observed here, and also one of the brightest (mM=3.4). It is much fainter in the K-band, hardly visible in the H-band (mK=10.6, mH=15.0), and also hardly extended at these shorter wavelengths (see Figs. 3 and 4). Consequently, IRS 3 appears to be dominated by dust emission. According to Krabbe et al. (1995), IRS 3 has an almost featureless, red continuum K-band spectrum. Also IRS 3 emits strongly at 10 $\mu$m (e.g. Gezari et al. 1985).

In our analysis, the M-band photocenter of IRS 3 is shifted to the northwest by $\sim$0.16 $^{\prime \prime }$ with respect to the H-, K- and L-band positions. This angular separation is well above the relative positional uncertainty of 0.035 $^{\prime \prime }$ we estimated above and thus not a mere artifact. A visual inspection of the maps clearly shows, however, that the two sources detected by the PSF-fitting algorithm are not distinct sources, but a consequence of the extended flux density distribution. This shift of the photocenter of the extended M-band emission relative to that at shorter wavelengths, combined with the clearly visible asymmetry of the lower contour lines in the NAOS/CONICA L-band image (Fig. 4) results in a bow-shock-like appearance, similar to the sources discussed by Tanner et al. (2003), although less distinct than the remarkable source mentioned by Clénet et al. (2004). The asymmetry is the consequence of a compression of the extended emission at its southern edge, resulting in a steeper flux density gradient ("bunched'' contour lines). A possible explanation for this unusual appearance is that IRS 3 consists of a hot mass-losing star, e.g. a young, dust-embedded O star (Krabbe et al. 1995) or a dusty protostar or Wolf-Rayet star (Horrobin et al. 2004; Tanner et al. 2003), surrounded by a very thick, extended dust shell, which is pushed northwest by wind from the direction of the IRS 16 cluster and Sgr A*. Gezari et al. (1985) give a temperature of $\approx$400 K for this strong L-band source. The H- and K-band emission is dominated by the stellar source and therefore not extended. The dust shell must be very thick, however, since IRS 3 is a faint source in K and especially H. The shift in the M-band emission is due to part of the dust shell being extended to the northwest. It should be possible to clarify the nature of IRS 3 through L- and M-band (and possibly longer wavelength) studies at higher resolution, i.e. using adaptive optics or interferometric methods that are becoming available in the infrared (e.g. VLTI or the upcoming LBT, see Eckart et al. 2003a).

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{1748fig7.eps}
\end{figure} Figure 7: Extinction towards Sgr A* according to Lutz et al. (1996). The data points are derived from SWS H recombination line data, while the solid curve is the extinction law of Draine (1989) for standard graphite-silicate mixes. The data point marked by an arrow is our new M-band result.
Open with DEXTER

   
4 Conclusion

We have presented multi-band infrared photometry of the Galactic Center, covering the large field of view (70 $^{\prime \prime }$ $\times $ 70 $^{\prime \prime }$) of the ISAAC instrument at the VLT. Our independent L-band magnitude calibration, relying on the robust method of observing standard calibration stars, reveals an offset to the traditionally used "standard'' calibrations (Depoy & Sharp 1991; Blum et al. 1996; Simons & Becklin 1996). We attribute this offset to the use of the variable star IRS 7 (Tamura et al. 1996; Blum et al. 1996; Ott et al. 1999) as a flux calibrator. This offset is important when investigating colors.

Together with new results on the extinction towards the Galactic Center (Scoville et al. 2003), our magnitude calibration makes the stars redder in K-L, which results in normal stellar properties. There is no longer any necessity to modify the K-band extinction in order to explain strange blue colors.

Color analysis of our large sample of sources is consistent with the results obtained by Scoville et al. (2003) concerning the extinction towards the Galactic Center, since we can compare the color distributions of sources lying within an inner 18 $^{\prime \prime }$ $\times $ 18 $^{\prime \prime }$ region with sources located further out, up to the 70 $^{\prime \prime }$ $\times $ 70 $^{\prime \prime }$ limit of the field of view of the ISAAC camera. This comparison reveals no systematic effects, indicating that the average extinction remains at the lower value of $A_{\rm V}=25$ mag up to the sharp inner edge of the circum-nuclear molecular ring.

Our M-band data indicates that the extinction in the M-band is higher than predicted by the extinction law of Rieke & Lebofsky (1985). The large number of sources for which we have obtained L-M colors allows a measurement of the M-band extinction, assuming that the intrinsic colors given by Koornneef (1983) are valid. This measurement gives $A_M =0.056A_{\rm V}$, i.e. $A_M\approx A_L$, a considerably higher value compared to what has been assumed so far.

Hot and cool stars show a fairly clear L-M color separation (Koornneef 1983), with the hotter stars exhibiting redder colors. Our results show that the Galactic Center sources follow this behavior, for which the H-K versus L-M color-color diagram (if available) is a useful diagnostic tool for preliminary identification of stellar types.

Several sources within the central few arcseconds are dust embedded and show extended L- and M-band emission. For one of the most prominent dust-embedded sources, IRS 3, we find extended L- and M-band continuum emission with a characteristic bow-shock-like appearance, including a shifted M-band photocenter relative to the other wavelengths. A possible explanation for this unusual appearance is that IRS 3 consists of a massive, hot, young and mass-losing star surrounded by an optically thick, extended dust shell, which interacts with the strong wind from the direction of the IRS 16 cluster and Sgr A*. As a result of this the dust shell is pushed towards the northwest by this wind.

It would appear that further studies in the few micron wavelength range could add considerably to the knowledge of the properties of the stellar population in the Galactic Center. This is especially true for the M-band, since it is a fairly long wavelength that is not completely dominated by dust emission. The stars therefore remain accessible for observation.

Acknowledgements

This work was supported in part by the Deutsche Forschungsgemeinschaft (DFG) via grant SFB 494. We are grateful to all members of the NAOS/CONICA team from MPIA/MPE, Meudon/Grenoble Observatories, ONERA, ESO. In particular, we thank N. Ageorges, K. Bickert, W. Brandner, Y. Clénet, E. Gendron, M. Hartung, N. Hubin, C. Lidman, A.-M. Lagrange, A. F. M. Moorwood, C. Röhrle, G. Rousset and J. Spyromilio.

References

 

Copyright ESO 2005