Skip to main content

Full text of "Composition of the galactic center star cluster: Population analysis from adaptive optics narrow band spectral energy distributions"

See other formats


Astronomy & Astrophysics manuscript no. ibnew 


© ESO 2009 


March 16, 2009 





o 
o 

(N 



< 
a 
^' 

Oh; 
O 

(N 
> 

(N 

rn 
O 

o 

> 



Composition of tlie galactic center star cluster 

Population analysis from adaptive optics narrow band spectral energy 

distributions 

R. M. Buchholz\ R. Schodel^ i, and A. Eckart^-' 

' I. Physikalisches Institut, Universitat zu Koln, Zulpicher Str. 77, 50937 Koln, Germany 

e-mail: buchholz , eckartOphl . uni-koeln. de 
- Instituto de Astroffsica de Andalucia (lAA)-CSIC, Camino Bajo de Huetor 50, E-18008 Granada, Spain 

e-mail: rainer@iaa.es 
^ Max-Planck-Institut fiir Radioastronomie, Auf dem Hugel 69, 53121 Bonn, Germany 



Received xx.xx.2009, accepted xx.xx.2009 



ABSTRACT 



Context. The GC is the closest galactic nucleus, offering the unique possibility to study the population of a dense stellar cluster 
surrounding a SMBH. 

Aims. The goals of this work are to develop a new method to separate early and late type stellar components of a dense stellar cluster 
based on narrow band filters, to apply it to the central parsec of the GC, and to conduct a population analysis of this area. 
Methods. We use AO assisted observations obtained at the ESO VLT in the NIR H-band and 7 intermediate bands covering the NIR 
K-band. A comparison of the resulting SEDs with a blackbody of variable extinction then allows us to determine the presence and 
strength of a CO absorption feature to distinguish between early and late type stars. 

Results. The new method is suitable to classify K giants (and later) as well as B2 main sequence (and earlier) stars which are brighter 
than 15.5 mag in the K band in the central parsec. Compared to previous spectroscopic investigations that are limited to 13-14 mag, 
this represents a major improvement in the depth of the observations as well as reducing the needed observation time. Extremely red 
objects and foreground sources can also be removed from the sample reliably. Comparison to sources of known classification indicates 
that the method has an accuracy of better than ~87%. We classify 312 stars as early type candidates out of a sample of 5914 sources. 
Several results such as the shape of the KLF and the spatial distribution of both early and late type stars confirm and extend previous 
works. The distribution of the early type stars can be fitted with a steep power law (J3y> = -1.49 ± 0.12, alternatively with a broken 
power law, y6i_io" = -1.08 ± 0.12, y8io-20" = -3.46±0.58, since we find a drop of the early type density at ~ 10")- We also detect early 
type candidates outside of 0.5 pc in significant numbers for the first time. The late type density function shows an inversion in the 
inner 6", with a power law slope of y6i;<6" = 0. 17 ± 0.09. The late type KLF has a power law slope of 0.30±0.01, closely resembling 
the KLF obtained for the bulge of the Milky Way. The early type KLF has a much flatter slope of (0. 14 + 0.02). Our results agree best 
with an in-situ star formation scenario. 

Key words. Galaxy: center - stars: early-type - stars: late-type - infrared: stars 



1. Introduction 

The Galactic center (GC) contains the densest star cluster in 
the Galaxy with a ~4 .0 x 10^ super massive black hole at 
its dynamical center (Eckart et al., 2002; Schodel et al.!, l2002i 
12003; Ghezetal., 2003, 2008; Gillessen et al., 2009). Around 
the black hole, the projected distribution of the stars can best 
be described by a broken power law (break radius Rhreak - 
6".0 + 1".0, all values given here are projected radii), with a 
power-law slope of F = 0.19 + 0.05 within the break rad ius and 
r = 0. 75 +0.10 towards the outside of the cluster (Scho del et al.L 
|2007|) . The Galactic center is l ocated at a distance of ~ 8.0 kpc 



seconds and/or the presence of a large number of bright early 
type stars. Adaptive optics assisted observations with high spa- 
tial resolution have shown that there is indeed an increased num- 
ber of early-type stars in t his region,wh ile the relative number of 
late type stars decrea s es dOenzel et al .'. 2003; Eisenhauer et al.L 
l2005t IPaumard et all 120061; iLu et aLl .2009) . Several authors 



jReidll 1993 1; lEisenhauer et al.L 120051; iGroenewegen et all l2008t 
lGhezetaU l2008; GillessenletaLl |2009|) from the sun. We will 
adopt this value throughout this work. 

The stellar composition of the cluster depends on the dis- 
tance to the center. This has been first observed as a drop in 
CO absorption strength towards Sgr A* in seeing limi ted ob- 
serva tions (Allen et al., 1990; Sellgren et al., 1990; Halle r et al l 
Il996l) . Two explanations have been discussed for this: a signif- 
icantly lower density of late type stars in the central few arc- 



have tried to explain this finding by collisions between stars (or 
between stellar mass black holes and stars), which may lead to 
the destructin of the envelopes of giant stars in the central region 
(Davies e tlll Il99ll; iBailev & Davie si il999l: jRasio & S hapiro , 
1 1990: .Davies et al.L Il998; Alexander et all [l999; .Dale et al . 

Several stellar populations have been detected in the central par- 
sec: the oldest observable objects that make up the bulk of the 
visible sources outside of the innermost few arcseconds are old, 
metal -rich M, K and G type giants with ages of 1 - 10 x 10^ 
years. The helium burning red clump sources are also present, 
although t they have not been discussed in detail until recently 
(Maness et al.L |2007|) . because older works on the stellar popu- 
lation did not reach the necessary lower magnitude limit. 



2 



R. M. Buchholz et al.: Composition of tlie galactic center star cluster 



A number of interm ediate-bright (ma^K ~ 10 - 12) stars 
that a re now on the AGB fcabb e et ailll995l:lBlu"m et all 1 19961 
l2003h have been produced by a star formation event > 100 x 10^ 
years ago. These stars can be distinguished from late type giants 
by the H2O absorption bands in their spectra that are detectable 
even at low spectral resolution. Very few supergiants like IRS 7 
are also present in the central parsec. 

Several objects with featureles s , but v ery red spectra have 
also been detected dBecklin et all 119781: iKrabbe et all 119951: 
iGenzel et all 1 19961 ). namely IRSIW, 3, 9, lOW and 21. With 
high-resolution imaging, most of these sources have been re- 
solved. They are mostly associated with the mini-spiral, and can 
be interpreted as young and bright stars with rapid mass loss in- 
teracting with t he interstella r medium in t he GC, so called bow - 
shock sources ("Tanner et al.! (l2002i |2003|) : iGeballe et all (l2004h . 
see also Perger et al. (2008) ). 

In the central ~0.5 pc there exists yet another distinct stel- 
lar population; massive , young stars create d in a starburst 3 - 
7 X 10^ years ago (e.g.. lKrabbe et al.Lll995h . These stars can be 
found, e.g., i n the IRS 16 and IRS 13 associa tions dEckart et al.L 
[2004; Maillardetali 120041; iLu et al.L I2005L e.g.,).The bright- 
est of those young stars have been described as stars in a 
transitional phase between O supergiants and Wolf-Rayet stars 
(WN9/Ofpe according to e.g. Allen et al. (1990)), with high 
mass-loss during this phase (Naiarro et al., 1994; Krabbe et al., 
1995; Morris et al., 1996; Najarro et al., 1997; Paumard et al., 
200 it [Moultaka et al.l 120051) . These stars account for a large 
part of the luminosity of the central cluster and also con- 
tribute half of the excitation / ionizing luminosity i n this region 
("Rieke, Rieke & Paul ', 1989; Naiarro et al., 1997; ' Eckart et all 
1999; Paum ard et al., 2006; Martins et al., 2007). Recently, 



Muzic et al. I (l2008h have identified a co-moving group of highly 



reddened stars north of IRS 13 that may be even younger objects. 
Besides the most massive early type stars, a large num- 
ber of OB stars with masses of ~ 10-6 Mp-) have been ex - 
amined bv .Levin & Beloborodovl (|2003|) ; iGenzel et all ( |2003|) ; 
iPaumard et all (l2006h : iLu et all (l2009l) . At least 50% of the 
early-type stars in the central 0.5 pc appear to be located within a 
clockwise (i n projection on the s ky) rotating disk, w hich was first 
detected by iL evin & Beloborodov (2003). Later, Genzel et al.l 
( l2003h ; iPaumard et al. (2006) claimed the existence of a sec- 
on d, counter-clock wise rotating disk. A very detailed analysis 
by lLu et al ] (I2009h . based on the fitting of individual stellar or- 
bits, shows only one disk and a more randomly distributed off"- 
disk population (with the numb er of stars in the disk similar to 
that on random orbits). Bartko et al.l (l2009i) . on the other hand, 
claim they at least observe a counter-clockwise structure that 
could be a strongly warped, possibly dissolving second disk. 



Table 1. Details of the observations used for this work. N is the num- 
ber of exposures that were taken with a given detector integration time 
(DIT). NDIT denotes the number of integrations that were averaged on- 
line by the read-out electronics during the observation. The scale of all 
observed images is 0.027" per pixel. 



date 


/l«„rra/[j"m] AA [jum 


N 


NDIT 


DIT[sec] 


29 April 2006 


1.66 


0.33 


31 


28 


2 


09 July 2004 


2.00 


0.06 


8 


4 


36 


12 June 2004 


2.06 


0.06 


96 


1 


30 


12 June 2004 


2.24 


0.06 


99 


1 


30 


09 July 2004 


2.27 


0.06 


8 


4 


36 


09 July 2004 


2.30 


0.06 


8 


4 


36 


12 June 2004 


2.33 


0.06 


99 


1 


30 


09 July 2004 


2.36 


0.06 


8 


4 


36 



distance to IRS 7 ["] 



Fig. 1. The strehl ratio, a value that measures the deviation of the PSFs 
of the sources from an ideal PSF, decreases with the distance to the 
guide star. 



In the immediate vicinity of SgrA*, there is yet another dis- 
tinct group of stars, which for m a small cluster of what ap- 
pear t o be early B-tyt) e stars (lEckart et al.L Il999t iGhez et all 
2OO3I lEise nhauer et al.l (20051) • These so-called "S-stars" stars 
are on closed orbits around Sgr A*, with velocities of up to a 
few thousand km/s and at distances as close as a few lightdays 
(Schodel et al., 2003; Ghez et al., 2003, 2005; Eisenhauer et al.L 
2005; Gh ez et al.LI2008l:IGillessen et al.LI2009l) . Their orbits have 
been used to determine the mass of the black hole and the dis- 
tance to the GC. 

How exactly star formation can take place in the central par- 
sec under the observed conditions is still a debated issue. 
Classical star formation from gas of the observed density is 
severely impeded by the tidal shear exert ed by the black hol e 
and the surrounding dense star cluster (iMorris et al.L Il993h . 
Two scenarios are being discussed to explain the presence of 
the early ty pe stars: [Genzel et al. (2003); Goodman (2003j); 
levin & Belo borodoy[l2003b ; IMilosavlievic & Loeb| (120041) ; 
Navakshin & Cuadral (|2005|) ; [Paumard et al.. (.2006) suggest a 
model of in-situ star formation, where the infall and cooling of a 
large interstellar cloud could lead to the formation of a gravita- 
tionally unstable disk and the stars would be formed directly out 
of the fr agmenting disk. An alternative scenario has been pro - 
posed bv lGerhardl (120011); iMcMil lan & Portegies Z warll (l2003h; 
Portegie s Zwart et al.1 (l2003h; iKim & Morris (2003); iKim et al.1 
(l2004i) ; IGuerkan & Rasiol (l2005b with the infalling cluster sce- 
nario, where the actual star formation takes place outside of 
the hostile environment of the central parsec. Bound, mas- 



sive clusters of young stars can then be transported towards 
the center within a few Myr (dynamical friction in a massive 
enough cl uster lets it sink in much more rapidly than individual 
stars, see iGerhardI (l2001h '). Recent data seem to favor contin- 
uous, in-situ star-form ation (e.g.lNavakshin & SunvaevI (l2005h ; 
IPaumard et al.1 ( l2006l) ; lLu"et al.l ( l2009l) 

The existence of the S stars so close to Sgr A * is yet an- 
other matter, known as the "paradox of youth" (iGhez et al.L 
|2003|) . Two explanations are discussed for the presence of these 
stars, though neither is satisfactory: formation out of collid- 
ing or inter acting giants (Eckart et al., 1993; Genzel et al., 2002t 
iGhez et aU 2003, 2005) or scattering from the disk of young 
stars ([Alexander & Livio i . i2004i) . 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



3 



• 



©58 

Q219 ^ ©12fe254 - 
©195 . . 



®6V 



Detailed spectroscopic studies of the stellar population in the 
central parsec have so far only been conducted in the innermost 
few arcseconds and on small areas in the outer region s of th e 
cluster (see |G hez et al. (2003, 2008); Eis enhauer et alj (l2005h : 
iPaumard et al.i (.2006.) ; ,Maness et al.. (.20070 ). Here the main lim- 
itation is that the high surface density of sources in the GC forces 
the observers to use high spatial resolution observations in or- 
der to be able to examine all but the brightest stars. However, 
the field-of-view of integral field spectrometers is quite small 
at the required angular resolutions (e.g., between 3" x 3" and 
0.8" X 0.8" for the ESQ SINFONI instrument). The aim of this 
work is to provide a reliable method that allows the classification 
of several thousands of stars as late or early type in the central 
parsec down to a magnitude limit of 15.5. These data will then 
be used to obtain constraints on the stellar population and exam- 
ine the distribution of early and late type stars. The distribution 
of very red objects will also be addressed. For this, information 
about the spectra of the stars is necessary. In ^2.11 we present 
our method of intermediate band imaging as well as the pho- 
tometry and calibration. The disadvantage of this method is the 
very low spectral resolution (only 7 datapoints over the K band 
and one H band point) that makes it impossible to detect any but 
the broadest spectral features. Also, line-of-sight velocities can- 
not be measured this way. The big advantage is the possibility 
to cover a large field of view of 40"x40" in a very time efficient 
way. The data analysis in ^2. 2112. 71 will be followed by a popula- 
tion analysis in ^ We summarize and discuss the implications 
of our results in ^ 

2. Observation and data reduction 

2.1. Observation 

The observations used here were carried out with the NAOS- 
CONICA (NACO) instrument at the ESQ VLT unit telescope 4 
on Paranal in June/July 2004 and April 2006 (programs 073. B- 
0084(A), 073.B-0745(A), 077.B-0014(A)). We used an H band 
broadband filter and seven intermediate band filters (see TabUli. 
The seeing varied for the different observations, within a range 
of 0.5 to 1.3". The bright supergiant IRS 7 located about 6" 
north of Sgr A* was used as the guide star for the adaptive op- 
tics (AO) correction, using the infrared wavefront sensor. The 
sky background was sampled by taking several dithered expo- 
sures of a dark cloud near the GC, 713" west and 400" north 
of Sgr A*, a region largely devoid of stars. A rectangular dither 
pattern was used for most observations, while some were ran- 



Table 2. Early type calibration stars, used for the primary calibra- 
tion (n ames, types and K band magnitudes according to Paum ard et alj 
( l2006h . approximate values for T^ff from Allen's Astrophysical 
Quantities) 



name 


ID 


mag K 


type 




E69 


61 


11.32 


early 




E55 


195 


12.45 


BO-II 


~ 20000 


E47 


219 


12.50 


B0-3I 


~ 18000 


IRS 16SSE1 


148 


11.90 


08.5-9.51 


~ 32000 


IRS 16SSE2 


199 


12.10 


B0-0.5I 


~ 20000 


IRS 16SSW 


71 


11.45 


08-9.51 


~ 32000 


IRS 33N 


58 


11.22 


B0.5-1I 


~ 20000 


E22/W10 


220 


12.73 


08-9.5I/III 


~ 32000 


E25/W14 


215 


12.58 


08.5-9.51 


~ 32000 


E43 


126 


12.10 


08.5-9.51 


~ 32000 


E53 


254 


12.31 


BO- 11 


~ 20000 



o 



d 

Q 

<a 
m 

'o 

in 

I 



5 0-5 
offset R.A.["] 

Fig. 2. Sources used for the primary calibration. Numbers correspond 
to the numbers in the common list, as also shown in Tab|2] 



domly dithered. All images were flatfielded, sky subtracted and 
coiTected for dead/hot pixels. 

In order to be able to separate early and late type sources using 
the method described in ^2.71 the photometry has to be accurate 
enough to clearly identify the feature used for the classification 
(see 92.41 1. This means that the typical photometric error should 
be much lower than the typical depth of the classification feature 
(see 92.41 for an estimate of the required accuracy). When ob- 
serving a large field of view like in this case, a good AO coiTec- 
tion can only be achieved within the isoplanatic patch, a region 
of ~10-15" for our dataset. This effect leads to a decrease of the 
Strehl ratio towards the outer regions of the field (see Fig[T]i. The 
values shown in this figure were computed with the strehl algo- 
rithm in the ESO eclipse software packag43from the PSFs deter- 
mined in 12 X 12 subimages of the IB227 image (with the same 
PSF stars that were used in the photometry). The parameters of 
the telescope (like the aperture) were also taken into account. 
The Strehl values exhibit a clear trend towards lower values at 
larger projected distances from the guide star Sources outside 
of the isoplanatic patch are elongated towards IRS 7. This is 
a problem when using PSF fitting photometry, while aperture 
photometry (which is less dependent on the shape of the PSF) is 
faced with the problem of crowding. 

2.2. Photometry 

In order to counter the aforementioned problems and to achieve 
reliable relative photometry over the entire FOV which is con- 
siderably larger than the isoplanatic patch, a two-step deconvo- 
lution process was used. Deconvolution is a good way to reduce 
source confusion in a crowded field, but it can only be applied 
satisfactorily if the PSF is very well known and un iform over the 
whole field. Th is process is described in detail in lSchodel et al.l 
din preparationh . who are using the same photometry. 
The individual lists of detected stars in the filters were merged 
to a common list of stars detected in all 8 filters. The number of 
common sources was limited by the image with the lowest qual- 
ity to 5914. This number includes almost all brighter sources 
(mag;i:.5 < 16), but several extended sources (like the bowshock 



' see N. Devillard, "The eclipse software". The mes- 
senger No 87 - March 1997, publicly available at 
fittp://www. eso. org/projects/aot/eclipse/distrib/index. htmi\ 



4 



R. M. Buchholz et al.: Composition of tlie galactic center star cluster 



+ 



2.1 2.2 2.3 

wavelength in mum 



2.4 




2.1 2.2 2.3 

wavelength in mum 



2.4 




2.1 2.2 
wavelength in mum 



2.4 




1.8 2 
wavelength in mum 



2.4 



Fig. 3. Conversion of continuous spectra into template SEDs. Upper left: K4.5 giant spectrum, the crosses indicate the low resolution SED, 
Ax = 3.3 mag. Upper right: BO main sequence spectrum, = 3.3 mag. Lower left: transmission curves of IB filters. Lower right: K4.5III and 
BOV template SEDs with added H band datapoint calculated from effective temperature and extinction. 



sources in the northern arm of the minispiral) were not detected 
in all filters with sufficiently low photometric and position un- 
certainties. 



2.3. Primary Calibration 

The primary calibration done here served a dual purpose: on the 
one hand, the measured counts for each source in each band were 
converted into a magnitude (absolute calibration). On the other 
hand, the bands had to be calibrated relative to each other to en- 
sure smooth spectral energy distributions (SEDs) as they can be 
expected at this spectral resolution (see Fig|3]l. Only very broad 
spectral features like CO bandheads and H2O absorption bands 
are expected to be observable here. Even the feature around 2.06 
fim that appears in the BOV spectrum only causes a decrease of 
0.014 mag in that filter compared to the K4.5III spectrum that 
lacks this feature. This can be neglected compared to other pho- 
tometric uncertainties. 

Atmospheric features can also influence the shape of the mea- 
sured SEDs and need to be eliminated. But the most important 
parameter that controls the quality of the data turned out to be 
the AO performance, which is why we did not use two of the 
available intermediate band datasets (2.12 and 2.18 //m). 
These goals were achieved by calibrating the common list of 
sources with 11 known OB stars (see Tab|2]and Fig|2]i close to 



Sgr A*. We adopted the classifications of lPaumard et al.l ( l2006l) . 
who provided a list of 90 early type stars in the central parsec. 
Stars of this type have some emission lines in the K band (see 
Fig|3]upper right), e.g. the Hel line at 2.058/im, but these lines 
are narrow enough to be negligible compared to the continuum at 
the spectral resolution of our data. Thus, the spectra of these stars 
can be assumed to be featureless, so they can be described by a 
blackbody spectrum with an effective temperature of ~30000 K. 
The final calibration has only a minimal dependance on the as- 
sumed effective temperature because the Rayleigh- Jeans law is 
a very good approximation for the SED of hot stars in the near- 
infrared (using 20000 K instead leads to a difference of only 0.01 
mag in the reference magnitudes). Therefore, the same T,,// can 
be used for all calibration sources. 

Extinction towards the GC is significa nt even in th e near- 
infrared. An extinction of ~3.0 mag ( Scoville et al 1 120031: 
ISchodel et al.L lin preparationh has been measured towards the 
central parsec in the K band. In order to minimize calibra- 
tion uncertainties, we used individual extinction values for each 
cahbration source, taken from the extinction map provided by 
ISchodel et al.l (lin preparation!) . Since the extinction also depends 
on the wavelength, we used thejOraine ( 1989) extinction law to 
calculate an extinction value for each band from the basic value 
for the K band, i.e. we assumed: 



A) oc A 



(1) 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



-1 — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 
M6III, A-K = 3.4, mag_K = 10.78, CO offset = 0.33 



-I- ftf*^ 



CM Li 1 L. 



2.1 2.2 
wavelength in mum 



2.3 



2.4 



T ' ' ' ■ 1 ' ' ■ ' 1 ' ■ ' ' 1 ■ ' •- 

K2.5III, A-K = 3.4, magj< = 15.37, CO offset = 0.1 B 



2.1 2.2 
wavelength in mum 



2.3 



2.4 




2.1 2.2 
wavelength in mum 



2.4 



in 



-1 — 1 — 1 — 1 — I — I — 1 — 1 — I — 1 — I — 1 — I — 1 — 1 — I — I — 1 — 1- 
G4III, A-K = 3.4, magj< = 16.77, CO offset = -0.03 



_i_ 



_i_ 



2.1 2.2 
wavelength in mum 



2.3 



2.4 



Fig. 4. K band spectra of late type stars expected in the GC. The crosses mark the corresponding SEDs and the green solid line a fitted extincted 
blackbody. As can be seen, the CO band heads of G giants are not deep and wide enough to have much influence on the low-resolution SEDs. 



Using an exti nction law with a different exponent (like /}"^" as 
proposed by iNishivama et al.l ([2008i) ) would lead to a general 
offset of about 0.5 mag in extinction, but neither the relative dis- 
tribution of the extinction values nor the results of our classifica- 
tion would vary, since the same extinction law is used again to fit 
for the individual extinction of each source in the classification 
process. 

The continuous extincted blackbody spectra were converted to 
8 point SEDs by multiplying them with the transmission curve 
of each filter (see Figl3] lower left). By comparing these theo- 
retical SEDs with the measured counts of the calibration stars, 
a calibration factor was calculated for each filter to convert the 
observed counts of each star into a magnitude. This eliminates 
the influence of any atmospheric features, since they should oc- 
cur in all sources and thus also in the calibration sources. The 
magnitude at 2.24/zm was adopted as the Ks band magnitude 
of the source, since this band is the closest available one to the 
center of the Ks band and not affected by absorption features. 
For the classification algorithm and the analysis, we used extinc- 
tion corrected magnitudes, i.e. we modified the measured mag- 
nitudes with the difference between the individual extinction of 
each source and the average extinction. This step eliminates the 
effect that the variable extinction in the observed region has on 
the brightness of the sources, which is necessary since the cut- 
off used for the classification depends on the brightness of the 
source in question. This step also eliminates the effects of spa- 



tially variable reddening on the derived luminosity functions. We 
assumed an average extinction in the central parsec of 3.3 mag 
dSchodel et aLllin preparation!) . 



mage 



magKs - Aks + A 



avg 



(2) 



Here, A^s denotes the measured individual extinction of each 
source in magnitudes (in the Ks band), Agvg = 3.3 mag is the av- 
erage Ks band extinction in the central parsec, magKs represents 
the measured Ks band magnitude at 2.24 jum. 
After the primary calibration, the sources within a few 
arcseconds of Sgr A* show SEDs that agree very well 
wit h SEDs calculated fro m template spectra of typical stars 
( Wa llace & Hinkl j, Il997h . The optimum of the calibration is 
centered around the region where the primary calibration sources 
are located (close to Sgr A*), and not around the guide star. With 
increasing distance from the optimum, the SEDs appear to show 
systematic deviations correlated over areas of a few arcseconds 
size. There are unfortunately not nearly enough known early 
type stars over the whole field to extend our primary calibration 
to the whole central parsec, so an additional local calibration had 
to be introduced that makes use of another class of stars that are 
easy to identify and abundant over the whole field: horizontal 
branch/red clump (HB/RC) stars. We assume that almost all late 
type stars with K band magnitudes between 14.5 and 15.5 are 
part of this population. This leaves the problem of telling apart 



6 



R. M. Buchholz et al.: Composition of tlie galactic center star cluster 



Table 3. Stellar types expected to be observable in the central par- 
sec, considering the distance and the extinction towards the GC. Values 
taken from Allen 's Astrophysical Quantities, the types shown here are 
the ones presented in that work that fall into the observational limits for 
the central parsec. 



early types 


mag K 




late types 


mag K 




B3Ia 


11.80 


-16000 


M6III 


-3200 


10.78 


06I 


12.12 


-36000 


M4III 


-3400 


12.32 


07V 


13.92 


-36000 


M2III 


3540 


12.92 


BOV 


14.65 


30000 


MOIII 


3690 


13.57 


B2V 


16.03 


20900 


K4.5III 


-4100 


14.29 


B3V 


16.33 


-19000 


K2.5III 


-4300 


15.37 


B7V 


17.51 


-13000 


G8.5III 


-4750 


16.42 








G7III 


-4900 


16.50 








G5III 


5050 


16.62 








G4III 


-5100 


16.77 








GOIII 


-5200 


17.27 



late and early type stars. The criterion for this is the same that is 
later used in the final classification; the CO band depth (CBD). 

2.4. CO band depth as a classification feature 

Since only very broad spectral features are visible at the low 
spectral resolution of our data, it is not possible to determine 
the exact spectral type of every single source. The feature used 
to distinguish between early and late types is the region beyond 
2.24 /urn, where late type stars show characteristic CO band head 
absorption. At the spectral resolution available here, no individ- 
ual lines or band heads are visible, but a broad feature like the 
CO absorption alters the shape of the SED in the corresponding 
region significantly (see Fig|4]i. The presence of CO band heads 
causes a significant dip in the SED for wavelengths greater than 
2.24 /urn. In general, the presence of this feature leads to the clas- 
sification of the source as a late type star with our method, while 
its absence makes the star an early type candidate. We did not 
consider young stellar objects with CO band head emission in 
our analysis. If any such objects exist within our data, they were 
treated as early-type stars. 

This feature also sets the limit for the required photometric ac- 
curacy: if a late type source with a Ks band magnitude of 15.5 
is expected to show a CO bandhead feature with a depth of 0. 1 
mag, this corresponds to a difference of in flux at 2.36 

fim and less at shorter wavelengths. Thus, the photometric error 
should be significantly lower than that value, which is achieved 
by our method for most sources. If the eiTor of a single data- 
point exceeded 15%, the source was excluded completely. 
We determined the presence of this feature in the following way: 
an extincted blackbody (Tgff = 5000 K, but see g2.7| for the ef- 
fective temperatures used in the final iteration of the fitting pro- 
cess) was fitted to the first 5 data-points (A < 2.27/vm). The ex- 
tinction was varied in a range of to 8 magnitudes in 0. 1 mag- 
nitude steps. The best fit was chosen based on the reduced of 
the fit. A third order polynomial was then fitted to the complete 
SED, with the first 5 data-points replaced by the fitted extincted 
blackbody to ensure a smooth fit. Although it might appear that 
using only the highest quality filters that are not influenced by 
possible additional spectral features (like the H band and the 
2.24 jum filter) would produce the best results, this would in fact 
be less reliable than the method we used here. Since the GC is a 
very crowded stellar field, variable AO performance and related 
variations of the Strehl ratio can cause variations of the mea- 
sured fluxes of stars. This can easily lead to outliers in the data. 



Therefore, it is safer to use a larger number of measurements in- 
stead of just two or three filters (see g2.7| for the special case of 
AGB stars with intrinsic H2O absorption features). 
For future applications of this method, we would prefer to also 
use the intermediate band filter centered at 2.12 yum. The data 
available at this wavelength had a low Strehl ratio and thus low 
photometric accuracy. Using the 2.18 jum filter is not advisable, 
since the data would be i nfluenced by the stro ng Br y emis- 
sion of the minispiral (e.g. lPaumard et af] (|2004|) ). It would also 
be possible to make use of the available narrow band filters to 
probe regions of special interest in the spectra, e.g. Br y emis- 
sion or absorption features in early type stars or the 2.20 jum 
Nal absorption feature in late type stars. But for that, a very 
careful background subtraction would be nece ssary to eliminate 
the str o ng influence of the minispiral (see e.g. feisenhauer et aU 
(I2005h : iPaumard et ail (2006), the same problem occurs with 
spectroscopy). Most important would be achieving improved 
photometric stability and deeper integration. 
The shape of both the fitted extincted blackbody and the third or- 
der polynomial depend on depend on the local extinction, but the 
difference between the extrapolated values of the two fits at 2.36 
jum is only extremely weekly dependent on extinction. We tested 
this by artificially reddening a typical late type SED (modeling 
extinction values of 1 .6-6.5 mag) and applying the fitting algo- 
rithm to it. The resulting difference of the polynomial and the 
blackbody varied by less than 0.15%, with slightly higher val- 
ues for lower extinction. Thus, we can assume that this value is 
independent of the local extinction within an acceptable margin 
of eiTor, which means it can be used as a good measure for the 
presence and depth of CO band heads and thus as a classification 
criterion. But in order to separate early and late type stars based 
on this value, which we will term CO band depth (CBD) in the 
following, a reliable cutoff is needed. 

2.5. Cutoff determination 

Although one could expect from template SEDs as shown in 
FigSa CBD close to zero for early type sources and greater than 
zero for late type sources, which could be separated by a simple 
general cutoff, the limited signal-to-noise ratio of the data, com- 
bined with additional sources of uncertai nty, such as photomet- 
ric errors due to source crowding (see e.g. lGhez et al.l ( l2008l) ) led 
to considerable scatter in the CBD feature. Therefore, the CBD 
value of early type stars scatter to positive and negative values. 



Table 4. Stars close to Sgr A* used to create the HB/RC template. All 
these sources were classified as late type by our algorithm, which is 
another self-consistency test. 



ID 


mag K 


3421 


15.6 


3369 


15.1 


2608 


15.6 


2888 


15.4 


3896 


16.1 


3565 


15.1 


4215 


15.3 


3253 


15.7 


2260 


15.3 


2248 


15.3 


1440 


14.6 


2206 


14.6 


4448 


14.8 


2158 


15.0 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



7 




-1 -0.5 0.5 1 



CO band depth 

Fig. 5. CO band dep th (CBD) at 2.36 iim for known early type sources 
found in this sample (Paumard et al.l i2006) . green) and known late type 
sources iManess et al. ( 2007), red), with the cutoff lines plotted in black 
(derived from early type sources for inner 12") and blue (based on late 
type sources, used outside of 12"). The cutoffs separate early and late 
type stars except a few outliers. The number of wrongly classified ref- 
erence sources can be used as an estimate for the relative uncertainty of 
the numbers of identified early and late type stars (see 93. 2t . 



But still, this value can be used as a useful tool to reliably sepa- 
rate early- and late-type stars statistically. 
In order to compensate for these difficu lties, we compared th e 
CBD values of the sources classified by iPaumard et alJ (l2006h : 
iManess et alj (lOOf) that are also present in our dataset. As can 
be seen in Fig|5j the early- and late-type sources mostly fall into 
separate regions when we plot the stellar magnitude vs. the CBD 
value. However, there are some (less than 5%) late-type sources 
present in the early-type region and vice versa. This makes it 
impossible to use a single smooth cutoff line, let alone a single 
cutoff value for the whole magnitude range. 
Almost all spectroscopically identified early-type sources are 
concentrated in the innermost 12.9". Therefore, two different 
smooth cutoff lines were calculated: one that reliably separates 
all known early and late type sources in the inner 0.5 pc, and 
one for the outer regions that encompassed the known late type 
sources there (see Fig|5]). This ensures a more reliable detec- 
tion of early type stars in the outer region (where the early type 
density is at best very low. The cutoff line based on the early- 
type stars will underestimate the number of early type stars, just 
as using the less strict cutoff based on the late-type stars might 
overestimate the number of early type stars. Use of both cutoff 
lines allows us to estimate the uncertainty of this method. 
We were able to achieve a clear separation of the reference 
sources with this two-cutoff-method, except a few outliers and 
noisy sources (see ^3.2b . 

The theoretical lower magnitude limit of this method is deter- 
mined by the presence of deep enough CO band heads so that the 
CBD significantly exceeds the photometric uncertainties. This 
was determined from the lWallace & Hinkld (Il997i) spectra to be 
the case for M, K and brighter G giants, corresponding to ~ 15.5- 
16 mag at the distance and extinction towards the GC. Fainter G 
giants and early type stars fainter than B2 have almost identical 
CBD values (see Fig|4]i. The comparison with published spec- 
troscopic identifications and the shape of the CBD vs. magnitude 
plot (Fig|6]) justified the adoption of a magnitude limit of 15.5. 



2.6. Local calibration 

With a classification criterion and a reliable cutoff in place, we 
were able to apply a local calibration that significantly improved 
the results. 

Previously publis hed K band luminosi ty functions of t he cen - 
tral parsec ( e. g. IGenzel et all (l2003h : IManess et all (l2007h : 
ISchodel et al.l (l2007h ') have shown that the magnitude range 
of 14.5 to 15.5 is dominated by HB/RC stars outside of the 
innermost arcsecond. In addition, the density of early type 
stars decreases stee ply outside of the innermost arcsecond 
(iGenzel et al.L l2003t IPaumard et al.L l2006t iLu et al.L l2009h . So 
even without knowing the exact type of every source in this mag- 
nitude range, they can be used as calibration sources assuming a 
typical late type SED. This template was determined as the av- 
erage SED of 14 manually selected late type reference stars (see 
Tab|4|i close to the optimum of the calibration. These stars were 
selected based on the similarity of their SEDs to the expected 
shape. 

After this, the HB/RC sources in the whole field except the in- 
nermost 2" were selected based on two criteria: an extinction- 
coiTected magnitude between 14.5 and 15.5 and a CBD value 
above the cutoff. The SEDs of these sources were compared to 
the template HB/RC SED corrected for the determined individ- 
ual extinction value. This yielded a calibration factor for each 
source in each band. Due to the uniform distribution of HB/RC 
sources over the field, these values map the local deviations and 
can be used to calibrate the SEDs of all sources. The calibration 
factors were checked for outliers first. If a calibration factor ex- 
ceeded the median over the closest 20 sources by more than 3 
cr, it was replaced by that median. On average, 20 sources are 
contained within an area of ~2 x 2", depending on the position 
in the field. This already leads to a spatial resolution of ~2". The 
factors were then processed into 8 calibration maps by smooth- 
ing the raw maps with a Gaussian with a 2" FWHM (see Fig[8]). 
The final cutoff to be used in the actual classification procedure 
was calculated in the same way as described in ^2.51 but with the 
locally calibrated values for the reference sources. This turned 
out to be a minor adjustment, so the local calibration process did 
not have to be reiterated. 

This local calibration eliminates most systematic local devia- 
tions and allows a more reliable classification of the sources to- 
wards the outer edges of the field (see|2]l. 




-0.5 0.5 1 

2.36 mum diff 



Fig. 6. CO band depth at 2.36 yum for the sources in the central parsec 
of the GC. Plotted here is the difference between the fitted blackbody 
and the third order polynomial at 2.36 yum. This serves as a criterion 
for the identification of early type (green) and late type (red) sources, 
separated by the inner (black) and outer cutoff (blue solid line). 



8 



R. M. Buchholz et al.: Composition of tlie galactic center star cluster 
source 133 source 133 



o 
E 
ro 
ri 

I 



1.6 



chr : 0.17 

difl: 0.168 / cutoff: 0.107 
extinction : 2.7 
type : late q1 

dmod : 0.43 
errmod : 0.01 




_i_ 



_i_ 



E 



Chi* : 1.1 

difl: -0.065 / cutoff: 0.112 
extinction : 3.2 
type : early q1 



dmod : 0.39 
errmod : 0.008 




_i_ 



1.8 2 
wavelength in mum 



2.2 



2.4 



1.6 



_1_ 



1.8 2 
wavelength in mum 



2.2 



2.4 



source 1148 



source 1148 



E 
II 



I 



Chi* : 0.27 , 

difl: 0.047 / cutoff: 0.102 
extinction : 4.4 
type : early q1 

dmod : 0.53 
errmod : 0.016 



E 
to 
ri 
II 



Chi* : 8.65 , 

difl: -0.159 / cutoff: 0.097 

extinction : 4.3 

type : early q1 

dmod : 0.51 
errmod : 0.014 

I L 



1.6 



1.8 2 
wavelength in mum 



2.2 



2.4 



1.6 



1.8 2 
wavelength in mum 



2.2 



2.4 



Fig. 7. Two examples for the effects of the local calibration. Upper right: source fitted as early type, but the SED is noisy which is reflected in 
the;^'^. Upper left: same source after local calibration, now fitted as late type, with a much better fit. Lower right: source fitted as early type, noisy 
SED, several data-points do not agree with the fit, large Lower left: same source, still fitted as early type, but with much better fit. 



2.7. Source classification 

The fitting process described in 32.41 was repeated on the cali- 
brated data. The extinction, CBD, the;^^^ of the fitted polynomial 
and the uncertainty of the CBD were calculated for each source 
in several iterations (see below). The reduced value for the 
last five data-points of the polynomial fit was used as a crite- 
rion for the quality of the fit. Using a reduced for all eight 
data-points makes less sense than the method used here, since, 
on the one hand, the first five data-points are replaced by the fit- 
ted blackbody for the polynomial fitting, and on the other hand, 
the first data-points do not have a large influence on the classi- 
fication anyway. In order to allow a realistic comparison to the 
cutoff for the final classification, the uncertainty of the CBD was 
calculated as the average root-mean-square deviation of the last 
five data-points to the polynomial fit. This value deliberately ig- 
nores the individual photometric uncertainties of the data-points, 
because sources with noisy SEDs with large uncertainties for the 
last data-points that are crucial for the classification would oth- 
erwise have a similar CBD uncertainty as sources with smooth 
SEDs with the last data-points close to the fitted blackbody. This 
could lead to erroneous classifications and is not desirable. 
We separated our sample of 5914 stars into the following classes: 

L Noisy sources: if the reduced of a source was higher than 
L5, we rated this source as too noisy for classification (334 



sources). The cutoff of L5 was chosen because it excludes 
the 5% most noisy sources. 

2. Foreground sources: these stars are not part of the popula- 
tion of the central parsec and have to be removed prior to 
any analysis. They are easily recognizable by their low fitted 
extinction value. Every source with a fitted extinction of less 
than 2 magnitudes was rated as a foreground source here and 
excluded. 58 sources were classified as foreground sources. 

3. Very red objects: several very strongly reddened objects (like 
e.g. IRS 3) can be observed in the field. Their SEDs are in- 
fluenced by other effects in addition to their intrinsic stel- 
lar features, e.g. dust shells and bow-shocks. These sources 
were also excluded from the further analysis since they can- 
not be compared easily to normal early or late type stars. 
We rated every source with a fitted extinction of more than 
5 magnitudes as such an object (24 sources). The high fitted 
extinction here is just a selection criterion that results from 
our algorithm. The real extinction toward the sources cannot 
been determined because of the unknown intrinsic SED of 
the very red objects. 

4. Early type sources (quality 1-3): the CBD of every source 
brighter than 16 magnitudes (extinction corrected) was com- 
pared to the applicable cutoff: sources closer than 11.9" were 
compared to the inner cutoff, while any source further away 
than 11.9" was compared to the outer cutoff. Sources with 



R. M. Buchholz et al.: Composition of the galactic center star cluster 




■5 




10 t> -10 

Qftlaet RA["] 

Fig. 8. Examples of the calibration maps used for local calibration (see 92.61 1 in the 8 bands. Left: 2.00 fim. Right: 2.33 fim. A central region 
with a radius of 2" was excluded, the calibration factors were set to there. These maps were generated by comparing a median spectrum of 
known HB/RC sources close to the optimum of the primary calibration with presumed HB/RC sources all over the field (see 92. 6t . Residuals of 
the rectangular dither pattern that was applied in the observations can be seen in the calibration maps. 



a CBD below the cutoff were rated as preliminary early 
type candidates and fitted again with Tj,// - 30000 K which 
comes much closer to the actual effective temperatures of 
early type stars than the 5000 K that were used in the initial 
fit. A 30000 K blackbody has a flatter slope in this wave- 
length regime than a 5000 K blackbody, so the fitted extinc- 
tion and the CBD that results from the fit also differ The 
updated CBD value was compared to the cutoff again. Only 
sources brighter than 15.5 magnitudes (extinction corrected) 
were considered. The extinction fitted in this iteration ex- 
ceeded the previous value by up to ~0.5 due to the flatter 
slope of the blackbody. This explains the higher magnitude 
cutoff used in the first step. A classification of fainter sources 
is not possible with this method, because the dep th of the CO 
band head feature becomes too shallow ( Wallace & Hinkl j, 
Il997h . If the CBD was more than 3 cr lower than the cutoff 
and the SED was not too noisy (reduced;if^ < 1.5), the source 
was rated as an early type with quality 1 (highest quality, 277 
sources). If CBD was between 2 and 3 cr lower than the cut- 
off, the source was rated as an early type with quality 2 (25 
sources). Sources with a CBD 1-2 cr below the cutoff were 
rated as early type quality 3 (10 sources). 
Late type sources (quality 1-2): all other sources that did not 
meet the cutoff criterion to qualify as an early type candidate 
were fitted again with T^jf = 4000 K. This value is typical for 
the expected late type giants. This again resulted in different 
extinction and CBD values. The uncertainties of the CBD 
and the redncedx^ of the fit were calculated in the same way 
as for the early type stars. If the CBD exceeded the cutoff by 
more than 1 cr and the magnitude of the source was brighter 
than 15.5, the source was rated as a late type star, quality 1 
(2955 sources), the others were rated as late type, quality 2 
(2231 sources). The late type quality 2 sources are neglected 
in the analysis of the late type population, since our method 
does not provide a clear identification any more. 
AGB stars: AGB stars are among the brightest sources in 
the field. Due to their prominent H2O absorption feature, 
they are usually fitted with a too flat blackbody/polynomial: 
the absorption feature leads to a higher magnitude at 2.00 
and 2.06 fim that drags the fitted blackbody and the poly- 
nomial "down". This can lead to a CBD below the cutoff 



and a classification as an early type star. To counter this 
effect, all sources were checked whether or not these two 
data-points were more than 1 cr below the fitted blackbody. 
If that was the case, the fitting was repeated without them. 
When the resulting classification changed, we adopted the 
updated class as the better result. We did not deem it neces- 
sary to put these objects in a separate sub-category, but we 
see this modification as a way to further improve the clas- 
sification of early type stars and remove false candidates. 
Source 1049 is a good example for this effect, especially 
since it ha s been identified sp ectroscopically as a late type 
source by iManess et al.l (1200 7). This addition to the algo- 
rithm resulted in 8 less candidates for early type stars in total, 
which does not have a big impact on our results as a whole. 



percentage of noisy sources 




r["] 

Fig. 9. Azimuthally averaged density of stars with noisy SEDs. The 
density profile is practically flat over the inner 15", with a slight increase 
further out. This indicates a good photometry and a reliable mechanism 
to identify noisy sources, since no region seems to have an over-density 
of noisy sources. This is also a good indication that the local calibration 
improved the results. 



10 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



o - 



o 

T 



T 1 1 1 1 1 1 1 r 

^ ■ ■. • ■ . V • -■. . J .. .• 



1 1 1 r 




20 



10 







-10 



Fig. 10. Map of the stars in the GC. Red circles represent late type stars, green circles stand for early type candidates. The blue circles are 
sources not classified, light blue indicates foreground sources, while magenta stands for extremely red objects. The radii of the circles are linearly 
dependent on the extinction corrected magnitudes of the sources. 



although it has a noticeable effect on the density distribution 
in the outer region of the observed area. 

Using only two different values for T^ff to fit early and late type 
candidates turned out to be sufficient, since the variations in the 
fitted extinction stayed within +0.1 mag when we used the min- 
imum resp. maximum values for Te/f as shown in Tab|3]for the 
respective classes of objects. We only determined the extinction 
in steps of 0. 1 mag in the first place, so we considered this slight 
uncertainty acceptable. 



3. Results and discussion 

3. 1 . Stellar classification 

In total, 3349 of 5914 sources have been classified as either early 
type, late type, foreground or very red sources. 223 1 of the re- 
maining sources are too faint to allow a reliable classification 
with our method, although they can be assumed to be mostly 
late type sources. The SEDs of 334 sources were too noisy and 
have been excluded by the classification algorithm. 
Unless otherwise indicated, all results shown here are based on 
the sources with an extinction corrected magnitude brighter than 
15.5 (see Eq|2]i. Foreground sources have also been excluded 
(58 in total), since they do not belong to the population of the 
central parsec. Very red objects (24 sources) have also not been 



included in the number of early or late type stars, but have been 
treated as an extra class of objects. Objects in this class can be 
of different type: highly extincted background objects (and thus 
not very relevant for an analysis of the ceii t ral p arsec), sources 
with dusty envelopes, e.g. IRS 3 (iPott et aI.L l2008l) . young stars, 
whose strong winds i nteract with the interstella r medium in the 
form of bow-shocks dTanner et aI.L 20021 2003h. or even candi- 
dates for young stellar objects ( Muzic et al. , 200 8*) (examples for 
the SEDs of these sources are shown in FigfTSTl. Tab|5] lists the 
stars in each (sub)category, while Fig[TO] shows the spatial dis- 
tribution of the classified sources. 



3.2. Comparison with spectroscopic results and uncertainty 
estimation 

We compared the results of our classification to the lists of 
spectroscopically classified sources provided by 'Pauma rd et al.l 
( 2006) an d Maness et al. (2007). The list published by 
I Paumard et al.l (120061) contains 90 early type stars, 78 of which 
were contained in our list of common sources. This discrep- 
ancy can be explained by the nature of the observations used 
here (different observation dates for the individual filters, 2 years 
between H band and IB observations, different quality of the 
datasets), some sources (especially fast moving objects like S2) 
were not detected at the same position in every image, not de- 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



11 



Table 5. Stars classified in the GC using the method described in 
SectlZTl 



class 


number 


early type quality 1 


277 


early type quality 2 


25 


early type quality 3 


10 


late type quality 1 


2955 


late type quality 2 


2231 


foreground 


58 


very red stars 


24 


noisy sources 


334 



tected through all filters or excluded due to too large photometric 
uncertainties and thus did not make it on the common list. The 
sa me effects are probab ly relevant for the 329 sources classified 
by iManess et al ] (l20()7h . while in addition, parts of this sample 
are outside of the area covered by our data. We only used the 266 
sources that are also present in our dataset for comparison. 
67 of the 78 known early type sources were also classified as 
early type here, 7 had too noisy SEDs and 4 were classified 
as late type. Of these 4 sources, 3 are borderline cases where 
a clear identification is very difficult with our method (sources 
224, 612 and 1772 resp. E36, £89 and E7 in the Paumard Hst). 
Source 3778 (E37 in the Paum ard list) shows a clear CO ab- 
sorption feature. P aumard et al.l ( i200 6) classify this source as a 
potential 08-9 supergiant, but at the same time give it a K band 
magnitude of 14.8 and an absolute magnitude of -3.3. This is in- 
consistent with our expectations, since such a source should be 
at least two magnitudes brighter (see Tab|3]l. We therefore ignore 
this source for the uncertainty estimation. This leads us to 3 out 
of 77 sources classified erroneously and 7 out of 77 sources not 
classified, which corresponds to 3.9% respectively 9.1%. A few 
well known sources like IRS 16SW and IRS 15NE have noisy 
SEDs, which in the case of IRS 16SW is probably due to the 
intrinsic variability of that source. But in general, noisy SEDs 
mostly stem from problems with the photometry: here, too faint 
or saturated sources are the biggest issues. 
Since the known early type sources are concentrated in the in- 
ner 0.5 pc, these values can be adopted as the uncertainties of 
the number of early type stars in the innermost few arcseconds 
identified in this work. 

258 of the 266 known late type sources have been classified as 
late types by our method. The SEDs of 7 sources were too noisy 
and one was classified as early type (source 363, 96 in Maness 
list). Source 363 does not show a clear CO feature despite being 
bright enough, and a comparison of the Maness and Paumard 
lists shows that there is an early type source 0.17" from its po- 
sition (assuming that the positions given in these works use the 
same reference frame). In the imaging data that we have used for 
this work, there is only a single source present at the location of 
source E87 (Paumard) resp. 96 (Maness). It appears sufficiently 
isolated to rule out confusi on with ano t her so urce. This leads us 
to the assumption that Paumard et alj (l2006h and iManess et alj 
(l2007h are looking at the same source there, but classify it dif- 
ferently. 

In order to derive an upper limit for the uncertainty and thus the 
confidence in our method, we assume one erroneous classifica- 
tion in the area covered by the SINFONI observations. In this 
region (north of Sgr A*, mostly outside of 0.5 pc), we find a to- 
tal number of 1 1 sources classified as early type, including the 
one star of disputed type. If we assume one of these classifica- 
tions to be erroneous, this leads us to an uncertainty of ~9% for 
our number of early type stars outside of 0.5 pc. We consider this 



an acceptable level of confidence, considering the low density of 
early type stars we measure this far out. 

To be on the conservative side, we adopt this value, ±9%, as the 
uncertainty of the total number of all early type stars, i.e. in the 
entire field of view. It has to be considered, however, that the 
low total number of sources available for the determination of 
this uncertainty level limits the confidence in it. 

3.3. Structure of the cluster 

FigHn shows the projected stellar density for the total popu- 
lation, the early and late type stars. For comparison, we also 
sho w the projected d ensity of the early type stars provided by 
Pau mard et al.l (l2006l) . Only stars brighter than 15.5 mag have 
been considered here, in order to allow a clear separation of early 
and late ty pe stars and to make a completeness correction unnec- 
essary (see lSchodel et all ( l2007h ). 

The projected density profile of the late type stars is practically 
flat within a radius of -10". Within the innermost 5", it can even 
be fitted with a power law with a positive slope, i.e. the projected 
density increases with the distance to the center. This flattening 
or even inversion of the projected surface density profile of the 
late-type stars combined with the steeply increasing density of 
early-type stars towards Sgr A* explains the dip in CO band head 
absorption strength fou nd in early spectroscopic observat i ons of 
low spatial resolu tion d Allen et all 1 19901; ISellgren et aLl 1 19901; 
iHalleretal .1119961) . We discuss this in detail in S3A\ 
A dip in the density can be observe d at a ra d ius of ~5" that 
has already been observ ed by Schodel et al.| ( |2007|) in deep 
{magK < 17.5) star counts. lZhu et al.l (12008) also find a dip at 0.2 
pc, wh ich corresponds to the 5" given here and inl Schodel et al.l 
(120071) . This dip is a significant feature in the density profiles. It 
is both present in the late- and early-type population, although 
with a low significance in the latter, due to the small number of 
early-type stars. The cause for this feature is probably extinc- 




Fig. 11. Azimuthally averaged stellar surface density plotted against 
the distance to Sgr A* for Ks magnitudes mag^, < 15.5. No complete- 
ness correction has been applied her e, but the data can be assumed to 
be complete down to ma.gxs ~ 15.5 jSchodel et all [2007). The green 
points describe the distribution of early type stars, while red stands for 
late type quality 1 stars and dark blue for all detected stars. This also in- 
cludes stars rated as noisy and bright enough late type quality 2 sources. 
The early type stars given in Paumard et alj ( [2b06) are shown in light 
blue for comparison. The solid lines indicate the power laws fitted to 
the data. 



12 



R. M. Buchholz et al.: Composition of the galactic center star cluster 




20 10 -10 



offset R.A.["] 

Fig. 12. Stellar surface density of all stars brighter than 15.5 mag (lin- 
ear scale, contours trace density in steps of 20, 30 ... 90 percent of max- 
imum density). The density is given in units of sources per arcsec^. 

tion, since there is a ring-like area of high extinction visible in 
our extinction map at this distance to the center (see FigfTTIi. 
Since this dip feature seems to appear in two stellar populations 
that are so different in their age and their dynamical state and 
since it would be very difficult to reproduce such a feature in a 
three-dimensional distribution, extinction seems to be the most 
likely explanation. 

A steep increase of the projected density of early type stars can 
be observed towards the center (except in the innermost arcsec- 
ond). This also confirms previous observations. Paumard et al. 
( 1200 6) proposed that the projected density of the early type stars 
follows a R^^ power law outside of a sharp inner edge at 1" 
(within the disk(s), so this value is not directly comparab le to our 
findings), with R as the projected distance to Sgr A*. iLu et alj 
(2009) also confirmed the R^^ power law within the clockwise 
disk. To allow a comparison with our values, we fitted the pro- 
jected density of the known early type sources that are contained 
in our sample with a power las as well. This yielded a value of 
/3ref = -1.80 + 0.17. In the following, we adopt this value as a 
reference. 

We fitted the projected early type density with a power law 
Pproj. °^ R^'" , excluding the inner 1". But while pro viding a value 
that a yees well with previo usly publi shed results dGenzel et al.L 
12003"; Pauma rd et all 120061: iLu et all 12009*), this single power 
law does not provide a very good fit to our data. For compari- 
son, we also fitted the projected early type density with a broken 
power law with a break at 10". This minimizes the deviations 
from the data, but it introduces a break at 10" that we cannot 
explain. 

/3y' = -1.49 + 0.12 
^1-10" = -1.08 + 0.12 
ySi()-20" = -3.46 ± 0.58 

It has to be kept in mind that the absolute number of early type 
stars in the outer regions is much lower than further towards the 
center, so the uncertainties are considerably larger The value for 
R > 1" agre e s reas o nably well with the values determined by 
iGenzel et all (l2003l) : iPaumard et al.1 (l2006l) : II u et al. I (l2009h . At 
R < 1", the early type density is lower than expected from ex- 



tending the power law inwards. T his is the sign of the inn er edge 
of the clockwise disk reported bv iPaumard et al.l (l2006l) and Lu 
et al. (2008). In this densest part of the cluster, the stellar surface 
density is probably also underestimated because source confu- 
sion will lead to incompleteness of the data (we estimate that 
completeness of sources at magx < 15.5 drops to ~ 80% within 
0.5" of Sgr A* ). Outside of 1", we detect a number of early type 
stars (312 over the whole field compared to 90) that is by a fac- 
tor of ~3 higher than in previous works. Their density profile is 
similar to that of the previously reported early-type sources. The 
power law in the inner few arcseconds becomes flatter with the 
additional early-type stars included, but agrees within the uncer- 
tainties with the previously reported values. Towards the edges 
of th e cluster (~10 - 20"), we observe a steeper density profile 
than IPaumard et al.l (l2006l) . It may be possible that a different 
density law applies outside of 10", but the cause for such a phe- 
nomenon is unknown. It has to be considered, however, that the 
statistics for such a small number of stars at larger distances are 
not very reliable any more. In addition, two different cutoffs were 
used for sources inside and outside of 12". This can also lead to a 
bias here in the way that the early type density is underestimated 
outside of 12" and overestimated within that distance to Sgr A*. 
But this effect should not influence the density in the order of 
magnitude observed here which leads to the different slopes of 
the fitted power laws. 

The relative stellar density of the sources rated as too noisy 
(see Fig|9] appears to be flat over an inner region with a radius 
of ~15". Outside of 15", the ratio increases, but that can be ex- 
pected since the quality of the photometry decreases somewhat 
toward the edge of the FOV. This is probably due to the rectangu- 
lar dither pattern used that leads to shallower integration toward 
the edges of the FOV. 

This gives further evidence that our criteria for excluding noisy 
sources as well as the local calibration we apply are justified, 
since an even distribution of noisy sources can be expected from 
a well calibrated dataset. This is the case here, and it means that 
the exclusion of noisy sources does not lead to a bias in our sur- 
face density profiles. 

FiglE] shows the two-dimensional density distribution of all 




20 10 -10 



offset R.A.["] 

Fig. 13. Stellar surface density of late type stars brighter than 15.5 mag 
(linear scale, contours trace density in steps of 20, 30 ... 90 percent of 
maximum density). The density is given in units of sources per arcsec^. 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



13 




Table 6. Power law indices for late type stars and all classified stars, 
separate fit to sources inside and outside of 6.0" 



10 
offset R.A.["] 

Fig. 14. Stellar surface density of early type stars brighter than 15.5 
mag (logarithmic scale, contours trace density in steps of 2.5, 5, 10, 
20 ... 90 percent of maximum density). The density is given in units of 
sources per arcsec" . 



stars brighter than magK - 15.5 and Fig[T3]that of the late type 
stars of the same magnitude range, while Fig[T4] shows the den- 
sity of the early type stars in the same way. These maps have 
been smoothed with a ~4" Gaussian. While the late type stars 
show a similar distribution as the stars in the cluster viewed as 
a whole (with the exception of the central few arcseconds), the 
early type stars are concentrated in the center This result is not 
surprising since the same can be seen from the azimuthally aver- 
aged density. The area in the immediate vicinity of IRS 7 shows 
a significantly lower stellar density in all our maps, because the 
presence of this extremely bright source impedes the detection 
of other stars close to it. 

The peak of the early type density is, as expected, located close 
to the position of Sgr A*. The distribution of the early type 
stars appears close to circularly symmetric, but indicates exten- 
sions along N-S and E-W. An apparent concentration of early 
-type stars a l ong th ese directions can also be seen in Fig. 1 of 
iBartkoet al.l(l2()09l) . 

The density distribution of the late type stars clearly shows the 
relative lack of late type stars in the very center. The map also 
shows a cor relation with the extin ction map (see extinction map 
provided in ' Schodel et al.l (|2007|) and FigfT?]). Areas of higher 
extinction show a lower density of stars. This effect is not visi- 
ble in the early type density. But in the areas with higher extinc- 
tion (e.g. 5" NW of Sgr A*), the early type density is too low 
for this effect to be relevant, although a trace of it can be seen 
in the slight dip in the early type density in FiglTT] The extinc- 
tion has two separate effects on the projected density: one effect 
is that sources behind a lot of extinction appear fainter, so they 
might be excluded by a simple magnitude cutoff. This effect has 
been corrected here by using extinction corrected magnitudes. 
The other problem is that extinction also impedes the detection 
of sources, especially in filters with lower image quality. This 
effect of the extinction on the stellar density has not been calcu- 
lated here and was not compensated. 







pinner 




/^Oilier 


all stars 


-0.22 


±0.11 


-0.86 


±0.08 


late type stars 


0.17 


± 0.09 


-0.70 


±0.09 



3.4. Evidence for giant depletion in tine center 

iFiger et al.l (120 03') examined the radial velocities of 85 cool stars 
in the GC (mostly M and K giants) and found dynamical evi- 
dence for a flattened dis t ributio n of late type s t ars wi thin 0.4 pc 
(~ 10"). As i Figer et all (l2003h and IZhu et all (|200^ point out, 
the flat projected surface density profile of the late-type stars im- 
plies in fact a hole/dip in the 3-dimensional distribution of the 
late-type stars . 

iSchodel et al.l ( l2007l) described the total population of the cen- 
tral parsec with a broken power law (break radius 'S.break - 
6" .0 + 1 ".0, y8 = -0. 1 9 + 0.05 within Rhreak and/3=-0.75 + 0.10 
outside of the break radius) for a magnitude limit of 17.75 and 
using completeness coiTected data. This does not allow a com- 
parison of absolute densities to our findings, since that dataset 
goes much deeper than ours, but the trend can be confirmed here. 
It also has to be considered that these results were obtained on 
the entire population and not separated into early and late type 
stars. 

We fitted broken power laws to the projected densities of the 
late type stars and all classified stars (see Tab|6]l. We find the 
same break radius of 6.0+1.0" as Schodel et al. (2007), and the 
power law indices for the total population also agree with the 
values given in that work (-0.22±0.11 for the inner region, - 
0.86+0.06 for R>6" in our data). What is new here, however, 
is the possibility to obtain separate values for the late type pop- 
ulation alone. These values give an even stronger support to the 
proposed hole/dip in the center: we find a power law index of 
-0.70+0.09 for the outer region, while the inner region even 
shows a decline towards the center (0.17+0.09). 

This is a very interesting result, since it shows that the pre- 
viously obs erved flattening of the density profil e of the total 
population dGenzel et al.l,ll996l:lF iger et al., 2003t lGenzel et all 
120031: ISchodel et all IIoOTI: IZhu et al.. ,2008 ) is an even stronger 
feature in the late type population. We can therefore assume that 
the stellar population in the innermost ~0.2 pc is indeed de- 
pleted not only of bright giants, but also of fainter giants down 
to our magnitu de limit of 15.5 ma g. Several causes for this have 
been proposed: iDale et al] (l2008h offered an explanation for the 
under-density of late type stars within 1", claiming that colli- 
sions with stellar mass black holes and main sequence stars pre- 
vent 1 -2 M0 giants to evolve so that they are not visible in the K 
band. Their simulations cannot explain the lack of brighter and 
fainter giants. 

Freitag et al. ( 2008h derived collision probabilities for bright gi- 
ants in the GC (see esp. their Fig. 1). They determined that nearly 
all massive stars within 0. 1 pc almost certainly suffer from colli- 
sions during their time on the giant branch. 



Table 7. KLF power law indices of the different classified stellar types, 
calculated for the entire cluster resp. the inner 7" 



type 


complete cluster 


r<7" 


all 


0.26 ± 0.01 


0.21 ± 0.02 


late 


0.31 ± 0.01 


0.27 ± 0.03 


early 


0.14 ± 0.02 


0.13 ± 0.02 



14 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



1000 11 


all sources: slope - 0.261 +/- 0.01 1 




g 






Z 






? 




' 1 , 


O 








I.I.I. 










Fig. 15. Luminosity function of stars in the central parsec. Each data-point represents the center of a magnitude bin (0.25 mag wide bins, resp. 
0.5 mag in lower right frame). Upper left: all stars, the HB/RC bump at ~ 15.0 mag is obvious. Central left: only late type stars, the HB/RC bump 
is visible. Lower left: only early type stars, much flatter luminosity function, no HB/RC bump. Upper right: all stars within 7". Central right: late 
type stars within 7". Lower right: early type stars within 7". The luminosity functions in the central few arcseconds of the cluster appear to have a 
flatter slope than that of the complete cluster, with the exception of the early type stars: here the slopes are almost the same. 



iMerritt & Szelll (l2006l) offered yet another explanation: the in- 
fall of a second SMBH would destroy the stellar density cusp 
present around Sgr A*, which would then be built up again in a 
time-frame of several Gyrs. This process can lead to a practically 
flat density profile, similar to the one observed here. 

3.5. K band luminosity function 

The general K band lumin osity function (FigfTSb agre es well 
with the one presented by e.g lGenzel et al.l (2Q03): Schod elet al.l 
( 12007 ). The red clump can clearly be made out at the expected 
magnitude of ~15.0-15.25. The 15.5 mag limit of the observa- 
tions does not have a significant effect on this feature, since this 
limit only applies to the separation of early and late type stars, 
while the photometric completeness limit lies at ~ 16 mag. 
The luminosity functions can be described to the first order by a 
power law (fitted to the area between 9.0 and 15.5 mag for the 
early type stars, resp. 14.5 to exclude the red clump in the other 
plots): 



The p ower law indices that resulted from this fit are shown in 
TablTl [Alexander & Sternberg(ll999l) : rTiede. Frogel & Terndrupl 



(119951) : IZoccali et all (l2003b measured the power law slope of 
the bulge population of the milky way sev eral degrees from the 
center as /3 ~ 0.3, while Fi ger et aTl (12004) give the same value 
for the KLF on 30 pc scales around the GC, claiming that this 
value agrees very well with the theoretical KLF of an old stellar 
population, reflecting the rate of ev olution of star s along the red 
giant and asymptotic giant branch. iGenzel et al.l (|2003|) give a 
value of p - 0.21 + 0.02 for the central 9". It is expected that for 
greater distances to the center, the slope value approaches that of 
the bulge. 

The power law fitted here for the total population of the central 
parsec is flatter than the one attributed to the bulge population, 
while ou r value fitted for the inner 9" matches the value of 0.21 
given by Genzel et all (|2003[) . 

The individual KLFs for the late and early types give the rea- 
son for this deviation from the bulge power law: since the cen- 
tral parsec (and even more the central few arcseconds) contain a 
significant number of early type stars and their KLF has a much 
flatter slope {Pearh - 0.14 + 0.02), the resulting power law is also 
flatter than that of the late type population alone. The power law 
fitted to the KLF of only the late type stars ifiune = 0.31 + 0.01) 
agrees very wel l with the one observed in th e bulge. 
Agreeing with iLebofskv & Riekd (Il987l) : iBlum et all (Il996l) : 



R. M. Buchholz et al.: Composition of the galactic center star cluster 
source 355 source 363 



15 




_i 1 I 1 I 1 I 1 d t I 1 I 1 I 1 I 1 : 

1.6 1.8 2 2.2 2.4 1.6 1.8 2 2.2 2.4 

wavelength in mum wavelength In mum 



source 613 source 672 



<N F r 




1.6 1.8 2 2.2 2.4 1.6 1.8 2 2.2 2.4 

wavelength In mum wavelength in mum 

Fig. 16. SEDs of stars fitted as early type outside of 0.5 pc. We consider these four examples to be very likely candidates for actual ea rly type 
stars. The source in the upper right frame does not show a CO absorption feature, but still it has been labeled as a late type source in Man ess et alj 
( l2007h . |Paumard et alj j2006h seem to list it as an early type however. The other sources listed in Tab[8]with quality A have similarly smooth SEDs, 
while quality B and C sources are considerably more noisy. 



iDavidge et af] (ll997h : lGenzel et al.1 (l2003h . our data point to the 
population in the central parsec being an old stellar population 
with an admixture of a young, bright component. 
The fitted slope value for the inner 9" of 0.2 1 differs from the one 
for the total cluster, but otherwise, the shape of the KLF there is 
very similar. This difference is due to the flatter slope of the late 
type KLF in this region, while that of the early type KLF stays 
the s ame. 

iPaumard et alj (|2006|) also presented a KLF for the early type 
stars in the disks. They also find a flat KLF, similar to the one 
presented here. The sig nificant difference is the greater sensitiv- 
ity we achieved here. Paum ard et al.l ( 1200^ ) note that their KLF 
has a spectroscopic completeness limit of ~13.5-14 mag, com- 
pared to 15.5 mag here. The fact that the flat slope of the early 
type KLF can be observed down to that magnitude is an impor- 
tant new result and it strengthens assumptions of a top-heavy 
mass function since it improves the statistic relevance of the ob- 
served flatness. 



icant problem for the reliable measurement of apparent mag- 
nitudes, colors and intrinsic reddening of stars. 'Schod el et al.l 
din preparation!) are presenting an extinction map based on H- 
K colors obtained from the same data that we are using here, but 
the individual extinction values obtained here can also be used 
to produc e an extinction map. FigllT] shows this map. It agrees 
well with lSchodel et al.l (1200 7'). except that we find higher over- 
all extinction values. The extinction map produced in this work 
is expected to be more reliable since we use many wavelength 
bands, apply a local calibration and distinguish between hot and 
cool sources. A histogram of the measured extinction values is 
shown in FigllT] The distribution is similar to a Gaussian with 
a mean value and standard deviation of 3.1 ± 0.4. This ag rees 



well with previously published results ( Scoville et al.L 120031) and 
also with ISchodel et al.l (lin preparation !). The asvmmetrv of the 
histogram is due to the exclusion of foreground stars. 



3. 7. Early type stars outside of 0.5 pc 



3.6. Extinction 



The average extinction of 3.3 mag towards the GC is well known 
(e.g. lScovilleetalJ(l2003h'). but this value varies on small scales 
by up to 1.5-2 mag (ISchodel et all |2007|) and poses a signif- 



In addition to not only identifying the well known early type 
stars in the center of the cluster, but indeed more than tripling 
the number of early type candidates in the central 0.5 pc, we can 
also report the identification of 35 early type candidates more 
than 0.5 pc (~12.9") from Sgr A*. If we apply our previously 



16 



R. M. Buchholz et al.: Composition of the galactic center star cluster 




I 




I \ r f J U 

10 

Fig. 17. Left: Extinction map of the central parsec computed from individually fitted extinction values. Features like the mini-cavity and the mini- 
spiral are visible. Right: Distribution of extinction values in the central parsec of the GC. The mean value and the variation agree with previously 
published results. Note that the Draine(1989) extinction law was used here and that a different law might apply to the GC. This may lead to a 
systematic offset as large as Ak = 0.5 mag, which should be taken into account as a possible systematic uncertainty when interpreting these values. 



derived uncertainty of ~9% (see 93.21 ). this leads us to a number 
of 35 + 3 early type candidates outside of the inner region. This 
number is most likely underestimated, due to the stricter cutoff 
criterion used in this region (see ^2.5\ . Within the observed area, 
the distribution of these sources appears to be fairly isotropic, al- 
though one has to be cautious here due to the the asymmetry of 
the observed region and the small number of candidates. The 
nature of these candidates needs to be confirmed with spectro- 
scopic observations, since our method can only provide a first 
estimate for the type. 

We obtain an average density of early type stars of (4.6 ± 0.4) x 
10^2 sources/arc sec^ atR > 12.9". 

The detection of ea rly type stars this far out in significant num- 
bers is a new result: [Paumard et alj ( 120 06') reported no early type 
stars outside of the central 0.5 pc, citing an 1 cr upper limit 
of ~10"^ OB stars per arcsec^ outside of 13" deduced from 
SINFONI data, but referring to an unpublished source. The diffi- 
culty in ideii t ifyin g these stars in existing data is pointed out by 
iTrippe et al.l |2008), who mention very limited coverage of the 
central cluster and data gathered with several instruments and in 
different epochs that is very difficult to compare due to differ- 
ent pixel scales, Strehl ratios and completeness. These obstacles 
have probably impeded the detection of early type stars in the 
outer region of the central cluster until now. Our new method can 
provide targets for spectroscopic confirmation observations, al- 
lowing for a broader search for early type candidates over a large 
area without the need to cover the whole area with integrated 
field spectroscopy. Finally, we would like to point out that we do 
know that at least one early type star exists outsi de of the cen- 
tral 0.5 pc, as has been confirmed by the results of lGeballe et alj 
(120061) on IRS 8. 

As FigHT] shows, there are two possible fitting solutions for the 
early type density. While the single power law fitted to the pro- 
jected density of early type stars agrees with the previously pub- 
lished power law, the slopes of the broken power law show a 
significant difference in the outer and inner region. If a single 
power law distribution can be assumed, this would indicate that 
these stars merely represent the continuation of the disk and the 
off-disk population. But if the much steeper decline outside of 
10" is indeed a significant and real feature, this might point to a 
change in the population respectively so far unknown effects on 



the density distribution. 



4. Conclusions 

Our newly developed method has confirmed several previously 
obtained results and has proven to be able to classify sources of 
known type reliably. It cannot compete with spectroscopic iden- 
tifications of individual sources, but due to the ability to classify 
a large number of sources with a relatively small need for ob- 
servation time, it can provide important statistical information 
about a stellar population which can later be refined by observ- 
ing the early type candidates with spectroscopic methods. Our 
analysis extends the sensitivity limit of stellar classification by 
about 1 .5 magnitudes compared to previous work, providing thus 
a statistically stronger basis for conclusions on the stellar popu- 
lation. 

The following results could be obtained: 

1. A larger number of early type candidates has been de- 
tected than in any previous study, 312 sources compared to 
the known 90. We were able to identify known early type 
sources with 87% accuracy, while 96% of the known late 
type sources were classified correctly. The different percent- 
ages stem from the different selection criteria of the refer- 
ence sources in the cited publications: Maness et al. (2007|) 
already selected their late type sources by the presence of 
sufficiently deep CO band heads, which is the same feature 
that our method makes use of . The early type sources pub- 
lished by Paumard et alj (12006 ) were selected based on nar- 
row emission lines that are not visible at our spectral resolu- 
tion. 

These high rates of correct classifications gives us high con- 
fidence in our results, but the new detections need to be 
confirmed spectroscopically, especially the early type can- 
didates. 

2. The spatial distribution of the early type stars follows a 
power law, with /3i" - 1.49 + 0.12, which is flatter than 
the values of the v alue of P - -1 . 8 we c o mputed f r om th e 
data published by Paumard et all (l2006h : iLu et alj ( |2009|) . 
but still compatible at the 3cr limit. It is also possible to 



R. M. Buchholz et al.: Composition of the galactic center star cluster 
aourea 136 source 5908 



17 




CM ^ I ^ I ^ I ^ I ^ I I ^ I ^ I ^ I 

1.6 1.8 2 2.2 2.4 1.6 1.8 2 2.2 2.4 

wavelength in mum wavelength in mum 

Fig. 18. SEDs of Extremely Red Objects. Upper left: source NW of Sgr A*, no CO feature, but strongly reddened. Upper right: source NW of Sgr 
A*, here a CO feature is visible. This source is located near a local maximum of extinction, so it may just be a normal late type star. Lower left: 
source located in the northern arm of the mini-spiral. Lower right: source located in the far SE of Sgr A*, also in a region of very high extinction. 
The CO feature suggest a late type star, maybe towards the back of the cluster. 



fit the early type density with a broken power law with 
/3i_iO" = -1.08 + 0.12 and ySio-2()" = -3.46 + 0.58. This 
broken power law fits the density distribution better than the 
single power law, but it remains unknown what causes the 
density drop we observe at ~10". The significance of this 
feature should be examined by further observations that ex- 
tend to larger distances from Sgr A*. 
3. We confirm the previously reported flat projected surface 
density profile of the late type stars in the innermost arcsec- 
onds with much higher source numbers and therefore signifi- 
cantly improved statistics. The flat (or even inversed) surface 
density of late type stars must imply a central dip or hole in 
their 3D distribution. This together with the steep early type 
density profile explains the observed drop in CO absorption. 
We would like to emphasize that this result signifies a change 
in our perception of the nuclear star cluster. Separating the 
early- and late-type population shows clearly that the GC 
cluster does not have a cM.;p (see discussions in lGenzel et al.l 
(l2003h : ISchodel et al.l(l2007l) ). To the contrary, the late-type, 
old stellar population, which appears to make up the vast ma- 
jority of stars in the GC cluster, shows a flat or even slightly 
inverted power-law in projection within about 0.2 pc from 
Sgr A*. This means that there is some kind of hole in the 
late-type population near the center, as has been pointed out 



bv lFiger et al.l (l2003h . The exact cause is still not understood, 
but various explanations have been suggested. One often dis- 
cussed possibility is the destruction of the envelopes of giant 
stars by collisions with main-sequence stars, post-main se- 
quence stars, and stellar remnants in the dense environment 



near Sgr A* dRasio & Shapirol 1 19901: iDavies et all 11991 



taenzel et al.', 1996' 
1998; Bailev & Davies 



Alexander et al 



3~r 



1999 

1999; Daleetal., 2008i Anoflier 
possible explanation for the non-existent cusp of late-type 
stars is that the cusp m ay have been destroye d by the infall 
of a second black hole (iMerritt & Szellll2006h and that there 
may not have been sufficient time yet to re-grow the cusp. 
The late type KLF has a power law slope of 0.30+0.01. 
This resembles closely the KLF that has been measured for 
the bulge of the Milky Way. This is surprising considering 
that the nuclear star cluster is pr obably a dyn amically sep- 
arate entity from the bulge (see Boker et akl (l2008l) for an 
overview of the properties of nuclear star clusters in galax- 
ies). It implies a similar star formation history for the NSC 
and the bulge. The early type KLF has a much flatter slope 
of (0.14+0.02). The fact that the flatter early type KLF could 
be confirmed down to our magnitude limit of 15.5 mag is an 
important extension of pr evious works that had a complete- 
ness limit of ~13-14mag. |Paumard et al. l( l2006l) claimed that 



Davies et al. 



18 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



their flat early type KLF agreed best with stellar evolution 
models using a top-heavy initial mass function (IMF), but 
did not give a value for the slope of the KLF itself (Fig. 13 in 
that work). Our early type KLF seems to show a very similar 
shape down to our magnitude limit and thus also appears to 
support the proposed top-heavy IMF, but an accurate com- 
parison is not possible due to the lack of a value for the slope 
in Paumard et al. (2006). 

5. Early type stars have been detected outside of 0.5 pc, in a 
density that still agrees with the power law density profile 
in the inner region, so these sources could well be part of 
the known disk/off'-disk population, if the single power law 
is valid. If the edge we observe at ~10" is a real feature, 
this may point to a more complicated situation. This result 
is of course pending spectroscopic confirmation. We find a 
larger number of early type sources in the whole cluster with 
the density dis tribution following the same power law as in 
iPaumard et al.l (2006 ). This means that our power law still 
agrees with our early type density outside of 0.5 pc that is 
higher than the upper limit provided in that work. 

6. Both foreground stars and strongly reddened objects could 
be detected or excluded easily due to their fitted extinc- 
tion. Unfortunately not all known bow-shock and mini-spiral 



Table 8. Early type stars detected outside of 0.5 pc. x and y denote 
the position of the star in arcsec relative to Sgr A*, magK^g„ the extinc- 
tion corrected K band magnitude and Ak the extinction in the K band. 
Quality indicates the confidence in the identification after an additional 
visual inspection, with A as the highest confidence and C as the lowest. 



x(") y(") magjf.„f Ak CBD quality 

C 
B 
B 
C 
A 
B 
B 
A 
B 
B 
A 
B 
A 
B 
B 
B 
B 
A 
C 
B 
A 
B 
A 
B 
A 
B 
A 
A 
B 
A 
C 
B 
B 
B 
C 



sources could be observed, due to position and photometric 
uncertainties. 

It is still not decided which one of the two main scenarios serves 
best to explain the presence of the early type stars in the central 
half parsec, in situ formation or infall and dissolution of a cluster 
formed at several parsecs distance from the GC. Recent publica- 
tions lea n towards favor i ng the scenario of in-situ star forma- 
tion (e.g. 'Paumard et"aLl (l2006h : iNavakshi n & SunvaevI (l2005h : 
Bartko et al. (200^. Our results also seem to agree best with 
this model: the steep power law decline of the projected early 
type density that agrees much better with the R"^ power law ex- 
pected for that scenario than with the R "-^^ of the in-spiraling 
cluster. It is important to note that the cluster infall scenario pre- 
dicts that l ess massive stars are strip ped from the cluster at larger 
distances dGuerkan & Rasioll2005l) . Our analysis probes, for the 
first time, the density of late-O/early B-type stars in the entire 
central parsec. These stars are less massive tha n the early-type 
stars reported from spectroscopic observations ( Paumar d et al.L 
2006; Bartko et al., 2009). Although we find a somewhat flatter 
power-law for the density of the early-type stars, the improved 
statistics still imply a steep decrease of the early-type stellar den- 
sity with distance from Sgr A*. This supports the in situ scenario, 
but it is not possible to clearly rule out or confirm one scenario 
based on the current data. 

Further observations that cover a larger area should be under- 
taken in the future, as well as spectro scopic confirmations of our 
new early type candidates. iLu et al] (|2009) suggest covering at 
least the inner 5 pc of the GC, and while observations of such 
an area require an extreme amount of time with an instrument 
like SINFONI, our method may be the key to explore large areas 
like this for features like tidal tails of an in-falling cluster or a 
continuation of the early type population of the central parsec. 

Acknowledgements. We are grateful to all members of the NAOS/CONICA and 
the ESO PARANAL team. R. Schodel acknowledges support by the Ramon y 
Cajal programme by the Ministerio de Ciencia e Innovacion of the government 
of Spain. We would also like to thank the referee for his helpful comments. 



References 

Alexander, T. 1999, ApJ, 527, 835-850 
Alexander, T., & Livio, M. 2004, ApJ, 606, L21 
Alexander, T., & Sternberg, A. 1999, ApJ, 520, 137 
Allen, D. A., Hyland, A. R., & Hillier, D. J. 1990, MNRAS, 244, 706 
Bailey, V. C, & Davies, M., B. 1999, MNRAS, 308, 257-270 
Bartko, H., Martins, F., Fritz, T. K., Genzel, R., et al. 2009, ApJ, accepted, 
arXiv:081 1.3903 

Becklin, E. E., Mathews, K., Neugebauer, G., & Willner, S. R 1978, ApJ, 219, 
121 

Blum, R. D., Sellgren, K., & DePoy, D. L. 1996, AJ, 112, 1988 

Blum, R. D., Ramirez, S. V., Sellgren, K., & Olsen, K. 2003, ApJ, 597, 323 

Boker, T., et al. 2008, Journal of Physics: Conference Series, Volume 131, 

Proceedings of "The Universe Under the Microscope - Astrophysics at High 

Angular Resolution", held 21-25 April 2008, in Bad Honnef, Germany. 

Editors: Rainer Schodel, Andreas Eckart, Susanne Pfalzner and Eduardo Ros, 

pp. 012043 

Dale, J. E., Davies, M. B., Church, R. P, Freitag, M. 2008, MNRAS, accepted, 
arXiv:0811.3111 

Davidge, T. J., Simmons, D. A., Rigaut, R, Doyon, R., & Crampton, D. 1997., 
AJ, 1 14, 2586 

Davies, M. B., Benz, W., Hills, J. G. 1991, ApJ, 381, 449-461 
Davies, M. B., Blackwell, R., Bailey, V. C, Sigurdsson, S. 1998, MNRAS, 301, 
745-753 

Diolaiti, E., Bedinelli, O., Bonaccini, D. Close, L., Currie, D., and Parmeggiani, 
G. 2000, A& A Supplement Series, vol. 147 (no. 2), Editions de Physique, 
p335 

Draine, B. T., Interstellar extinction in the infrared (Infrared Spectroscopy in 
Astronomy, Proceedings of the 22nd Eslab Symposium held in Salamanca, 
Spain, 7-9 December, 1988. Edited by B.H. Kaldeich. ESA SP-290, European 
Space Agency, 1989, p.93), 93-+ 



-8.15 


13.50 


10.7 


3.3 


0.07+0.04 


-9.59 


-15.33 


12.5 


3.1 


-0.05+0.01 


15.74 


12.87 


13.1 


2.8 


0.03+0.03 


-8.37 


-11.27 


13.9 


2.7 


-0.095+0.019 


-8.64 


-10.74 


13.5 


3.0 


-0.003±0.003 


-10.57 


-15.90 


13.0 


3.3 


0.05+0.01 


18.92 


0.28 


13.3 


3.2 


0.02+0.01 


-10.62 


-9.76 


13.2 


3.3 


0.008±0.017 


-10.81 


-10.51 


14.1 


3.0 


0.01±0.01 


13.05 


14.29 


13.3 


3.5 


0.04+0.02 


7.39 


12.29 


13.3 


3.5 


0.03+0.01 


13.13 


5.08 


14.6 


2.8 


-0.07±0.04 


19.58 


-16.54 


12.7 


4.0 


0.032±0.005 


4.43 


16.11 


13.9 


3.4 


0.028±0.007 


9.43 


9.29 


13.6 


3.5 


-0.008+0.018 


18.70 


4.44 


14.0 


3.3 


-0.08±0.03 


-10.86 


7.20 


13.3 


3.8 


0.05±0.02 


11.30 


8.47 


14.0 


3.4 


0.018±0.004 


17.69 


3.01 


14.0 


3.4 


0.03±0.01 


10.91 


9.44 


11.9 


4.7 


-0.07+0.03 


-9.05 


-11.77 


15.2 


2.9 


-0.089±0.005 


-12.75 


9.77 


12.8 


4.4 


0.05±0.02 


17.48 


11.36 


13.7 


3.9 


0.036±0.008 


-4.20 


14.06 


14.9 


3.2 


-0.08±0.02 


17.51 


-8.95 


14.8 


3.4 


-0.07±0.01 


9.89 


-14.14 


13.2 


4.4 


0.04±0.03 


14.25 


-5.64 


14.5 


3.8 


-0.033+0.005 


12.55 


-4.19 


13.2 


4.7 


0.03+0.02 


-15.25 


6.68 


14.1 


4.1 


-0.09+0.06 


-13.13 


-12.44 


15.1 


3.6 


-0.079±0.004 


-5.14 


-15.75 


15.4 


3.4 


-0.13±0.02 


4.49 


-15.53 


15.2 


3.6 


-0.14±0.03 


-1.16 


16.91 


13.8 


4.5 


-0.06±0.02 


0.31 


-15.30 


15.1 


3.7 


-0.11±0.03 


8.86 


9.55 


15.4 


3.5 


-0.12±0.02 



R. M. Buchholz et al.: Composition of the galactic center star cluster 



19 



Eckart, A., Genzel, R., Hofmann, R., Sams, B. J., & Tacconi-Garman, L. E. 1993, 

ApJ, 407, 77 

Eckart, A., Ott, T., Genzel, R., & Lutz, D. 1999, Proceedings lAU Symposium 

No. 193 

Eckart, A., Genzel, R., Ott, T., & Schodel, R. 2002, MNRAS, 331, 917 
Eckart, A., Moultaka, J., Viehmann, T., Straubmeier, C, & Mouawad, N. 2004, 
ApJ, 602, 760 

Eisenhauer, P., Genzel, R., Alexander, T., et al. 2005, ApJ, 628, 246^259 
Figer, D., Gilmore, D., Sungsoo, S., Morris, M., Becklin, E. E., McLean, I., 

Gilbert, A., Graham, J., Larkin, J., Levenson, N., Teplitz, H. 2003, ApJ, 599, 

1139-1156 

Figer, D. R, Rich, R. M., Sungsoo, S. K., Morris, M., & Serabyn, E. 2004, ApJ, 
601, 319 

Freitag, M., Dale, J. E., Church, R. R, Davies, M. B. 2008, Proceedings lAU 

Symposium No. 245, 2007 
Geballe, T. R., Rigaut, R, Roy J.-R., Draine, B. T. 2004, ApJ, 602, 770 
Geballe, T. R., Najarro, R, Rigaut, R, Roy J.-R. 2006, ApJ, 652, 370-375 
Genzel, R., Thatte, N., Krabbe, A., Kroker, H., & Tacconi-Gharman, L. E., 1996, 

ApJ, 472, 153 

Genzel, R., Schodel, R., Ott, T., Eisenhauer, F, et al. 2003, ApJ, 594, 812 
Gerhard, O., 2001, ApJ, 546, L39 

Ghez, A. M., Duchene, G., Matthews, K., Hornstein, S. D., Tanner, A., et al. 

2003, ApJ, 586, L127 
Ghez, A. M., Salim, S., Hornstein, S. D., Tanner, A., Morris, M., Becklin, E. E., 

& Duchene, G. 2005, ApJ, 620, 744 
Ghez, A. M., Salim, S., Weinberg, N. N., Lu, J. R., Do, T., Dunn, J. K., Matthews, 

K., Morris, M., Yelda, S., Becklin, E. E., Kremenek, T., Milosavljevic, M., & 

Naiman, J. 2008, ApJ, 689, 1044 
Gillessen, S., Eisenhauer, F, Trippe, S., Alexander, T., Genzel, R., Martins, P., 

Ott, T. 2009, ApJ, 692, 1075 
Goodman, J. 2003, MNRAS, 339, 937 

Groenewegen, M. A. T., Udalski, A., Bono, G. 2008, A&A, 481, pp.441-448 
Guerkan, M. A. & Rasio, F A., 2005, ApJ, 628, 236 

Haller, J. W., Rieke, M. J., Rieke, G. H., Tamblyn, P, Close, L, & Meha, R 1996, 

ApJ, 456, 194 
Holland, S., et al. 1998, AJ, 115, 1916-1920 
Kim, S. S. & Morris, M., 2003, ApJ, 597, 312 
Kim, S. S., Figer, D. F & Morris, M., 2004, ApJ, 607, L123 
Krabbe, A., et al. 1995, ApJ, 447, L95 

Lebofsky, M. J., & Rieke, G. H. 1987, AlP Conf. 155, The Galactic Center, ed. 

D. C. Backer (New York:AlP), 79 
Levin, Y. & Beloborodov, A. M., 2003, ApJ, 590, L33 

Lu, J. R., Ghez, A. M., Hornstein, S. D., Morris, M. R., Becklin, E. E., 2005, 
ApJ, 625, L51-L54 

Lu, J. R., Ghez, A. M., Hornstein, S. D., Morris, M. R., Becklin, E. E., Matthews, 

K., 2009, ApJ, 690, 1463 
Maillard, J. R, Paumard, T, Stolovy, S. R., Rigaut, R 2004, A&A, 423, 155-167 
Maness, H., Martins, F, Trippe, S., Genzel, R., et al. 2007, ApJ, 669, 1024-1041 
Martins, R, Genzel, R., HilUer, D. H., Eisenhauer, P., Paumard, T., et al. 2007, 

A&A, 468, 233-254 
McMillan, S. L. W. & Portegies Zwart, S. R, 2003, ApJ, 596, 314 
Merritt, D., Szell, A. 2006, ApJ, 648, 890-899 
Milosavljevic, M. & Loeb, A., 2004, ApJ, 604, L45 
Morris, M. 1993, ApJ, 408, 496 

Morris, R W., Eenens, R. R J., Hanson,. M. M., Conti, P S., Blum, R. D. 1996, 

ApJ, 470, 597 

Moultaka, J., Eckart, A., Schodel, R., Viehmann, T., Najarro, F. 2005, A&A, 

443, 163 

Muzic, K., Schodel, R., Eckart, A., Meyer, L., Zensus, A. 2008, A&A, 482, 173- 
178 

Najarro, R, Hillier, D. J., Kudritzki, R. P, Krabbe, A., Genzel, R., Lutz, D., 

Drapatz, S., Geballe, T. R. 1994, A&A, 285, 573 
Najarro, F, Krabbe, A., Genzel, R., Lutz, D., Kudritzki, R. R, HiUier, D. J. 1997, 

A&A, 325, 700 
Nayakshin, S. & Cuadra, J., 2005, A&A, 437, 437 
Nayakshin, S., & Sunyaev, R. 2005, MNRAS, 364, L23-L27 
Nishiyama, S., Nagata, T, Tamura, M., Kandori, R., Hatano, Sato, S., Sugitani, 

K. 2008, ApL 680, 1174 
Paumard, T, Maillard, J. P, Morris, M., Rigaut, R 2001, A&A, 366, 466-480 
Paumard, T, Maillard, J.-R, Morris, M. 2004, A&A, 426, 81 
Paumard, T, Genzel, R., Martins, R, et al. 2006, ApJ, 643, 1011-1035 
Perger, M., Moultaka, J., Eckart, A., Viehmann, T, Schodel, R., Muzic, K. 2008, 

A&A, 478, 127 

Portegies Zwart, S. F, McMillan, S. L. W., Gerhard, O., 2003, ApJ, 593, 352 
Pott, J.-U., Eckart, A., Glindemann, A., Schodel, R., Viehmann, T., Robberto, 

M. 2008, A&A, 480, 115 
Rasio, F A., & Shapiro, S. L. 1990, ApJ, 354, 201-210 
Reid, M. J. 1993, ARA&A, 31, 345-372 



Rieke, G. H., Rieke, M. J., & Paul, A. E. 1989, ApJ, 336, 752 
Schodel, R., et al. 2002, Nature, 419, 694 

Schodel, R., Ott, T, Genzel, R., Eckart, A., Mouawad, N., Alexander, T. 2003, 
ApJ, 596, 1015 

Schodel, R., Eckart, A., Alexander, T, et al. 2007, A&A, 469, pp.125-146 

Schodel, R., Eckart, A., Buchholz, R., et al., in preparation 

Scoville, N. Z., Stolovy, S. R., Rieke, M., Christopher, M. & Yusef-Zadeh, F. 

2003, ApJ, 594, 294 
Sellgren, K, McGinn, M. T, Becklin, E. E., & Hall, D. N. 1990, ApJ, 359, 112 
Tanner, A., Ghez, A. M., Morris, M., et al. 2002, ApJ, 575, 860 
Tanner, A., Ghez, A. M., Morris, M., Becklin, E. E. 2003, ASN, 597-603 
Tiede, G. P, Frogel, J. A., & Temdrup, D. M. 1995, AJ, 110, 2788 
Trippe, S., Gillessen, S., Gerhard, O. E., Bartko, H., Fritz, T. K., Maness, H. L., 

Eisenhauer, F, Martins, F, Ott, T, Dodds-Eden, K., Genzel, R. 2008, A&A, 

492, 419 

Wallace, L., Hinkle, K. 1997, ApJS v.lll, p445 
Zoccali, M., et al. 2003, A&A, 399, 931 

Zhu, Q., Kudritzki, R., Figer, D., Najarro, R, Merritt, D. 2008, ApJ, 681, 1254- 
1278