A&A 433, 117-125 (2005)
DOI: 10.1051/0004-6361:20041748
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
70
)
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 0.006),
i.e.
,
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 0.17
for emission-line stars and -0.05 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
The Galactic Center stellar cluster shows some intriguing characteristics: it is extremely dense, with an unusual observed stellar population consisting mainly (80% of the 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 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 70 ) 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 13 (Clénet et al. 2001; Davidge et al. 1997), appear to have neglected the M-band (4.66 m), which shows some interesting properties, which we discuss in this paper.
Figure 1: ISAAC M-band final image of the Galactic Center. The field of view is 72.6 72.6 . | |
Open with DEXTER |
Figure 2: Part of the ISAAC L-band final image of the Galactic Center showing most of the named sources. | |
Open with DEXTER |
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).
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).
Notes:
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 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 , 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 mag would result in even bluer K-L colors, requiring an even larger offset to the L-band calibrations mentioned above.
For de-reddening purposes, the extinction law of Rieke & Lebofsky (1985) was adopted, i.e. , , and (initially) . 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 mag was assumed, instead of the traditional average value of 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.
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 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 70 . 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).
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:
= | |||
= |
Assuming a Gaussian distribution of the L-M values for
simplicity, we find a corresponding standard deviation of
,
which allows the calculation of the contribution
of screening and source intrinsic effects to the L-M distribution.
The uncertainties of the distribution are given by
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 |
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 |
Figure 5: Two-color diagram showing versus . The triangles represent sources from the inner 18 18 , the squares represent brighter sources (mL<9.5) and the diamonds a selection of fainter sources from outside this inner region. | |
Open with DEXTER |
Figure 6: Two-color diagram showing versus . The triangles represent sources from the inner 18 18 , 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) and colors of the He-stars. These appear to be slightly but systematically redder (0.2 mag) than expected for hot stars, whose and 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.
In our analysis, the M-band photocenter of IRS 3 is shifted to the northwest
by 0.16
with respect to the H-, K- and L-band positions.
This angular separation is well above the relative positional uncertainty of 0.035
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 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).
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 |
We have presented multi-band infrared photometry of the Galactic Center, covering the large field of view (70 70 ) 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 18 region with sources located further out, up to the 70 70 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 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 , i.e. , 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.