arxiv: v2 [astro-ph.ga] 22 Dec 2009

Transcription

1 Near-Infrared Counterparts to Chandra X-ray Sources Toward the Galactic Center. II. Discovery of Wolf-Rayet Stars and O Supergiants arxiv: v2 [astro-ph.ga] 22 Dec 2009 J. C. Mauerhan 1, M. P. Muno 2, M. R. Morris 3, S. R. Stolovy 4, A. Cotera 5 ABSTRACT We present new identifications of infrared counterparts to the population of hard X-ray sources near the Galactic center detected by the Chandra X-ray Observatory. We have spectroscopically confirmed 16 new massive stellar counterparts to the X-ray population, including nitrogen-type (WN) and carbon-type (WC) Wolf-Rayet stars, and O supergiants. These discoveries increase the total sample of massive stellar X-ray sources in the Galactic center region to 30 (possibly 31). For the majority of these sources, the X-ray photometry is consistent with thermal emission from plasma having temperatures in the range of kt = 1 8 kev or non-thermal emission having power-law indices in the range of 1 Γ 3, and X-ray luminosities in the range of L X erg s 1 ( kev). Several sources have exhibited X-ray variability of several factors between observations. These X-ray properties are not a ubiquitous feature of single massive stars but are typical of massive binaries, in which the high-energy emission is generated by the collision of supersonic winds, or by accretion onto a compact companion. However, without direct evidence for companions, the possibility of intrinsic hard X-ray generation from single stars cannot be completely ruled out. The spectral energy distributions of these sources exhibit significant infrared excess, attributable to free-free emission from ionized stellar winds, supplemented by hot dust emission in the case of the WC stars. With the exception of one object located near the outer regions of the Quintuplet cluster, most of the new stars appear isolated or in loose associations. Seven hydrogen-rich WN and 1 Spitzer Science Center, California Institute of Technology, Pasadena, CA 91125, USA; mauerhan@ipac.caltech.edu 2 Space Radiation Laboratory, California Institute of Technology, Pasadena, CA 91125, USA 3 Department of Physics and Astronomy, University of California, Los Angeles, CA , USA 4 Spitzer Science Center, California Institute of Technology, Pasadena, CA 91125, USA 5 SETI Institute, 515 N. Whisman Rd., Mountain View, CA, USA

2 2 O stars are concentrated near the Sagittarius B H ii region, while other similar stars and more highly evolved hydrogen-poor WN and WC stars lie scattered within 50 pc, in projection, of Sagitarrius A West. We discuss various mechanisms capable of generating the observed X-rays and the implications these stars have for massive star formation in the Galaxy s Central Molecular Zone. 1. Introduction The Galactic center harbors the largest concentration of X-ray sources in the sky. Indeed, a total of 9017 X-ray sources have been detected in deep observations of the central of the Galaxy with the Chandra X-ray Observatory (Muno et al. 2009). The majority of these sources are highly absorbed by interstellar gas and dust, indicative of distances near the Galactic center (D 8 kpc; Reid 1993). A wide variety of stellar phenomena comprise the Galactic center X-ray population, including accreting magnetic white dwarfs in cataclysmic variables (CVs; e. g., see Muno et al. 2003) and symbiotic binaries; accreting neutron stars and black holes in low- and high-mass X-ray binaries (e. g., see Pfahl, Rappaport & Podsiadlowski 2002); Wolf-Rayet and O-supergiant binaries with colliding supersonic winds; and, perhaps, rogue compact objects accreting from the dense Galactic center medium (Agol & Kamionkowski 2002). All such objects are rare among stellar populations, and involve exotic forms of matter and high-energy radiation mechanisms, of which our knowledge is fragmentary. Thus, the characterization of the Galactic center X-ray population is an important step toward improving our understanding of the end stages of stellar evolution, and the zoology of X-ray sources contained within the only galactic nucleus we are currently able to resolve in such detail. X-ray properties alone, however, are insufficient to characterize the Galactic center X- ray population. The detection of counterparts at other wavelengths is necessary, although such a pursuit is limited by the high extinction suffered by starlight traversing 8 kpc from the Galactic center to Earth (A V 30 mag). Fortunately, the brightness distribution of infrared counterparts can be used to constrain the relative contributions of low- and high-mass objects to the sample. For instance, CVs with cool dwarf donors in the Galactic center should have very faint IR counterparts with K 22 mag, while main-sequence stars earlier than B0V, and Wolf-Rayet stars, will have K 15 and 12 mag, respectively. It was demonstrated in Mauerhan et al. (2009; hereafter referred to as Paper I) that only 6%±2% of the 6760 absorbed X-ray sources have real infrared counterparts with K s 15.6 mag. Although this result is consistent with earlier studies that concluded the population to be widely dominated by CVs (Muno et al. 2003, Laycock et al. 2005), it suggests, nonetheless, that

3 3 of the X-ray sources in our sample have real infrared counterparts within a brightness range consistent with that of the population of late-type giants and hot massive stars observable in the Galactic center. Indeed, several X-ray sources with hot supergiant counterparts have already been discovered, including O supergiants, and nitrogen-type (WN) and carbon-type (WC) Wolf-Rayet stars (Muno et al. 2006, Mikles et al. 2006, Mauerhan, Muno & Morris 2007, Hyodo et al. 2008). The successful detection of massive-star counterparts to Chandra sources is not surprising, given the environmental conditions at the Galactic center. The region boasts the highest star formation rate density in the Milky Way (10 7 M yr 1 pc 3 ; Figer et al. 2004), owing to the immense reservoir of molecular gas that occupies the Central Molecular Zone (e.g., see Morris & Serabyn 1996). This material is vigorously forming massive stars near the Sagitarrius B H ii region (e. g., see depree et al. 1998), and within the last several Myr has produced at least three extraordinary stellar clusters that are among the most massive and dense in the Galaxy: the Arches and Quintuplet (Nagata et al. 1995; Cotera et al. 1996, Figer et al. 1999, 2002), and the Central parsec cluster (Krabbe et al. 1995). These clusters are rich in Wolf-Rayet stars and O supergiants. Such stars are typically sources of soft, thermal X-rays with kt 1 kev, while only a small subset of massive stars (those in close binaries) will produce hard X-rays with kt > 1 2 kev. Since the softer X-ray photons are heavily absorbed by interstellar gas and dust, only the harder sources are likely to be detected at the Galactic center. This explains why the Arches and Quintuplet clusters contain only several detectable X-ray sources among their hundred or so massive stars (Yusef-Zadeh et al. 2002; Law & Yusef-Zadeh 2004; Wang et al. 2006). Outside of these clusters, the relatively isolated X-ray-emitting supergiants discovered by Mauerhan, Muno & Morris (2007), which reside within 10 projected pc of the Arches and Quintuplet clusters, may have been dynamically ejected from one of these systems. Alternatively, these and other isolated massive stars throughout the Galactic center may be the products of an isolated mode of massive star formation operating in the region, in tandem with the formation of stellar clusters. Thus, the identification of infrared counterparts to the Chandra X-ray population can provide insight into the evolution of stellar clusters, and highlight previously unknown regions of massive star formation in the Galactic center. In this work, we present the results of our pursuit of counterparts via infrared spectroscopy. In Section 2, we list the various facilities and instrumentation used for this work and describe our observations. In Section 3, the spectra of the newly found massive stars are presented, including new spectra of several previously discovered objects. In Section 4, the stellar parameters of confirmed sources are examined via their photometry and infrared spectral energy distributions (SEDs). In Section 5, the X-ray photometric properties of the massive stars are listed and compared with the greater field population of X-ray sources in

4 4 the Galactic center. In Section 6, we discuss potential hard-x-ray emission mechanisms, consider the spatial distribution, origin, and formation mode of these stars, and discuss the implications for the greater population of massive stars in the Central Molecular Zone. 2. Observations We performed infrared spectroscopy of candidate counterparts to X-ray sources that were selected via the cross-correlation of the Chandra catalog of Muno et al. (2009) and the SIRIUS and Two Micron All Sky Survey (2MASS; Cutri et al. 2003) near-infrared catalogs. The details of the sample and the cross-correlation analysis are presented in Paper I. We targeted sources that (1) have highly absorbed soft X-ray fluxes and reddened near-infrared photometry (the so-called red infrared matches to hard X-ray sources from Paper 1), which is the characteristic of objects lying at the distance to the Galactic center, (2) lie at angular distances >7 from Sgr A, thus avoiding the high stellar confusion and large number of accidental infrared/x-ray matches that lie within this radius (3) have K s < 12 mag, and (4) have excellent X-ray astrometry, with positional uncertainties 1. We have currently obtained spectra of 52 infrared matches to X-ray sources. We observed in the K-band since it contains many diagnostic features of both hot and cool stars (e.g., see Hanson et al. 1996, 2005; Morris et al. 1996; Figer, McLean & Najarro 1997), and suffers the least amount of extinction in the near infrared. Eighteen sources, including those from Mauerhan, Muno & Morris (2007), exhibit characteristics of massive stars (see Section 3), while the remainder exhibit CO band-head absorption features near λ2.3 µm, and sodium absorption features near λ2.21 µm, typical of cool, late-type stars (spectra not presented here). In general, we assumed that a late-type match to an X-ray source was a spatial coincidence, since latetype stars vastly outnumber hot stars in the field. Thus, late-type counterparts will not be presented in this paper. However, we note that symbiotic X-ray binaries consisting of latetype giants or supergiants with accreting compact companions could be indistinguishable from random coincidences if emission-line accretion signatures were too weak to be detected in our spectra. Thus, we cannot completely rule out the presence of symbiotic stars in the late-type sample. We do not claim that the true statistics of infrared counterparts to X-ray sources are indicated by the number of massive-star detections presented below; there are several reasons for this: our spectroscopy campaign began before the completion of the Chandra survey and the final source catalog, which in the end allowed us to refine the X-ray point-source positions and uncertainties, and eliminate spurious detections (the method of refining the X-ray astrometry is discussed in Muno et al. 2009). Thus, several matches for which we had

5 5 obtained prior spectroscopy did not remain in the final candidate counterpart list that was based on the refined X-ray catalog. In fact, only 15 of the 35 sources which we confirmed spectroscopically to be late-type stars remained in the final list. Furthermore, candidate counterparts were also given priority for spectroscopy based upon their location in particular regions of interest. For instance, if a candidate s position on the sky places it near an interesting diffuse structure, such as a bubble-like or shell-like feature in the Spitzer/IRAC λ8 µm images of the region, then it was given priority over other candidate counterparts in the list. The motivation for this additional selection criterion lies in the fact that intense winds from massive stellar X-ray sources can sweep up and heat interstellar material. Thus, we interpreted the spatial association of candidate infrared/x-ray counterparts with diffuse mid-infrared structures as indicative of a higher likelihood that the X-ray source is associated with a massive star. Although this technique probably increased our success rate, it almost certainly imposed a selection effect that hampers our ability to statistically analyze our spectroscopic success rate for finding counterparts to the Chandra population. Table 1 lists the X-ray sources and associated massive-star counterparts that we targeted for spectroscopy, along with the instrument used, and the dates of the observations. The table also includes information on seeing and sky conditions. All sources are presented with their full Chandra designation in Table 1, and are hereafter referred to in the text by an abbreviated version (e.g., CXOGC J will hereafter be referred to as X ) NASA Infrared Telescope Facility (IRTF) and SpeX Spectra of four stars were obtained using the SpeX medium-resolution spectrograph on the 3 m Infrared Telescope Facility (IRTF) telescope (Rayner et al. 2003), located on the summit of Mauna Kea in Hawaii. For X and X , SpeX was used in the high-throughput, low-resolution prism mode, which provides coverage of the entire JHK bandpass, although we will only include the K-band in this work. A slit width of 0. 3 was used, providing a spectral resolution of R 250. For the infrared counterparts to X and X , medium-resolution spectra were obtained using the short cross-dispersed mode (SXD) of SpeX, with slit widths of 0. 5 and 0. 3 for the respective stars, providing spectral resolutions of R 1200 and R 2000 in the K-band. All spectra were acquired in an ABBA nodding sequence in order to subtract the sky background, and to suppress the contribution of bad pixels. The spectra were reduced and extracted using the IDL-based software package Spextool, specially designed for the reduction of data obtained with SpeX on the IRTF (Cushing et al. 2004). For telluric correction, spectra of the A0V standard

6 6 star HD were obtained and applied to all of the low-resolution data, while the A0V star HD was used for the medium-resolution spectra. Telluric correction was executed using the IDL package xtellcor (Vacca et al. 2003), which applies and removes model H i absorption lines from the A0V standard star before application to the science data The United Kingdom Infrared Telescope (UKIRT) The 3 m United Kingdom Infrared Telescope (UKIRT) on Mauna Kea was used to obtain spectra of eight infrared counterparts: X , X , X , X , X , X , X , and X The UKIRT 1 5 µm Imager Spectrometer (UIST; Ramsay et al. 2004) was used in service mode as part of the UKIRT Service Programme. The short-k grism and 4-pixel slit (0. 48) were used, which provided a spectral resolution of R 2133 and a wavelength range of µm. Spectra were acquired in an ABBA nodding sequence for sky subtraction. Basic reductions were executed by the Starlink ORACDR pipeline and the spectra were extracted using the IRAF routine APALL. The telluric spectrum of HD was applied to all UKIRT/UIST data using the xtellcor general program, a generalized version of the same program used for IRTF SpeX data (used for all stars not observed at the IRTF) The Anglo-Australian Telescope and IRIS2 Spectra of X , X , and X were obtained using the 4.1 m Anglo- Australian Telescope (AAT) on Siding Spring Mountain (Mount Woorat) in New South Wales, Australia. The IRIS2 instrument (Tinney et al. 2004) provided a spectral resolution of R 2400 in the K-band, using the 1 slit. The data images were reduced using the Starlink ORACDR pipeline. The spectra were extracted using the IRAF routine APALL. The telluric spectrum was obtained from observations of the F0V star BS The Southern Observatory for Astrophysical Research and OSIRIS Spectra of X , X and X were obtained at the Southern Observatory for Astrophysical Research (SOAR), located on Cerro Pachon in Chile. The Ohio State Infrared Imager Spectrometer (OSIRIS; Depoy et al. 1993) was used in high-resolution longslit mode, which provided R 3000 in the K-band. Stellar spectra were acquired in a slitscan sequence of five positions separated by 5 each. The five spectral images were median

7 7 combined with a suitable bad-pixel rejection algorithm to produce a master sky spectrum, which was subtracted from each individual exposure. The A0V standard HD was used for telluric calibration. 3. Spectral Classification In the following subsections, we estimate the spectral types of the infrared counterparts to Galactic center X-ray sources. All of the spectral criteria we use to classify these stars were adopted from previous spectroscopic studies of massive stars in the near-infrared, conducted by Morris et al. (1996), Hanson et al. (1996, 2005), Figer, McLean & Najarro (1997), and Martins et al. (2008). The central wavelengths of all spectral lines referred to in the text are adopted from Figer, McLean & Najarro (1997, their Table 2) Late O Supergiants Figure 1 shows the K-band spectra of the infrared counterparts to X , X , and the Quintuplet member X (star qf344 of Figer, McLean & Morris 1999). Each of these stars exhibits absorption lines of Brγ at λ µm and He i at λ2.058, λ2.113, and λ µm, which are typical features of late O and early B supergiants. In addition, each star also exhibits a weak emission feature of N iii, which appears at λ2.115 µm, on the red side of the He i absorption line. The main feature distinguishing the stars of Figure 1 is the varying depth of the Brγ line, which tends to increase with later spectral type. The Brγ feature is weakest for X , deepest for X , and is of intermediate strength for X The latter star exhibits a weak absorption line of He ii as well, which indicates that it is somewhat hotter and of slightly earlier spectral type than X and X Basedonacomparison ofthese spectra withtheob starsof Hansonet al. (1996, 2005), we classify X and X as supergiants within the O9I B0I range, and X as an O8 O9I star. We note that the X-ray detection of X is questionable. The source detection was flagged as spurious and not included in the catalog of Muno et al. (2009), but owing to the source s position near an interesting shell-like feature in Spitzer λ8 µm images, it was targeted for spectroscopy nonetheless. This implies that either some legitimate X-ray detections near the threshold of significance have been discarded by our selection algorithm, or the spatial correlation of the OB star with a truly spurious X-ray detection was a very unusual coincidence.

8 Early O Supergiants The infrared counterparts of X , X , X , and X all exhibit very similar K-band spectra, presented in Figure 2. The spectra are dominated by a complex of blended emission lines near λ µm that is dominated by He i and N iii, and might also may contain contributions from C iii and O iii. Blueward of this line complex is weaker emission lines of C iv at λ2.069 and λ2.078 µm. Brγ is seen in emission in X , although it may be nebular feature, while there is marginal detection of Brγ absorption in X Both X and X exhibit He ii absorption at λ2.189 µm, which is either absent or below the noise level in the other two stars. All of these features are consistent with the K-band spectra of early O supergiants, specifically O4 6I stars (e. g., see Martins et al. 2008) and we classify them as such. The Arches cluster of the Galactic center is rich with stars of O4 6I spectral type (Figer et al. 2002, Najarro et al. 2004, Martins et al. 2008). For comparison, accompanying the spectra in Figure 2 is a model spectrum for the O4 6I star F18 of the Arches cluster from Martins et al. (2008). This model was computed using the stellar atmosphere code CMFGEN (Hillier & Miller 1998). The model star has the following values for effective temperature, luminosity, mass-loss rate, and terminal wind velocity of these stars: T eff = K, L = L, Ṁ = M yr 1, and v = 2150 km s 1. It is reasonable to assume that the stars in Figure 2 have similar properties WNh and OIf + Stars The K-band spectra of X , X , and X are presented in Figure 3. The spectra of these stars are dominated by emission lines of Brγ and the λ µm complex of He i, N iii, C iii and O iii. The spectrum of X was first presented in Mauerhan, Muno & Morris (2007; hereafter referred to as MMM07), where it was classified as an O6If + star, but we present a refined version of it here. The Brγ line of this particular star exhibits an asymmetry on its blue side owing to a contribution from He i, although the detectability of this feature may be due to the fact that the spectrum of X has significantly higher spectral resolution than the other stars in Figure 3. An absorption line of He i at λ2.058 µm is present for all three stars, although the detection is marginal for X Weak N iii emission appears at λ2.247 µm in the low-resolution spectra of X and X , while this line appears as a strong doublet feature in the higher resolution spectrum of X Weak N v emission is also observed in all three stars, although its detection in the higher resolution spectrum of X is the most secure. X also exhibits a weak He ii line (λ2.189 µm) in absorption.

9 9 These features are consistent with O4 6If + and core-h burning WNL stars of the hydrogen-rich variety (WNh, Smith et al. 1996). OIf + and WN8 9h stars with relatively weak lines can be very difficult to distinguish from one another using near-infrared spectroscopy alone, as these two spectral types occupy opposite ends of a continuous morphological sequence (Bohannan & Crowther 1999). WN8 9h stars are usually distinguished by He ii (λ2.189 µm) in emission or as a P Cygni profile. However, several OIf + have been shown to exhibit this line in emission as well, while several WN8 9h stars are known to exhibit this line in absorption (Conti et al. 1995; Bohannan & Crowther 1999). An additional classification criterion, more relevant for the K band, is provided by the Brγ emission line, which is usually stronger than the λ µm complex for WN8 9h stars, relatively weak in OIf stars, and of comparable strength in OIf + stars; based upon this criterion, we classify X and X as WN8 9h stars, and maintain the original O6If + classification from MMM07 for X WN8 9h stars are also well represented in the Arches cluster. A model spectrum of the WN8-9h star F7 from Martins et al. (2008) is included in Figure 3 for comparison with the other WN8 9h stars. The basic parameters of this model are T eff = 32,900 K, L = L, Ṁ = M yr 1, and v = 1300 km s 1. Figure 4 presents the K-band spectra of X , X , and X These stars are distinguishable from the WN8 9h stars of Figure 3 by their more prominent He ii emission at λ2.189 µm, indicative of slightly earlier spectral types in the range of WN7 8h (e. g., see Martins et al. 2008). The higher signal-to-noise spectrum of X also exhibits emission features of N v near λ2.10 µm, and C iv emission near λ2.078 µm. In conclusion, we classify X , X , and X as WN7 8h stars Hydrogen-Deficient WN Stars Figure 5 presents K-band spectra of WN stars X , X , X , and X These stars have hydrogen-deficient spectra containing very broad, flat-topped emission lines of He i and He ii, indicative of fast, extended stellar winds. The discovery of X was first reported in Cotera et al. (1999), and it was first identified as an X-ray source in Muno et al. (2006). We present a new, higher resolution K-band spectrum of this star and examine its spectral type. For X we measure EW(λ2.189 µm)/ew(λ µm)=0.52 and EW(λ2.189 µm)/ew(λ2.112 µm)=0.55, which is consistent with stars of subtype WN7, according to Figer, McLean & Najarro (1997); this classification updates the earlier WN6 classification in

10 10 Cotera et al. (1999). By direct comparison, the spectrum of X is very similar to that of X , also exhibiting EW ratios consistent with WN7 classification. Slight differences in the spectrum of X include a stronger He ii emission at λ2.189 and λ2.314 µm, and a deeper P Cygni absorption component from He i at λ2.058 µm. The He i absorption component for X also has larger blueshift than that of X , implying that the wind of X may be faster and more extended than that of X The counterpart to X was first presented in MMM07, where it was classified as a WN6b star. However, a refined spectral analysis of this star demonstrates that the original telluric correction and continuum fitting resulted in an underestimate of the He i (λ µm) emission-line flux. The more reliable result presented here exhibits a stronger emission line at this wavelength, with a blueshifted absorption profile from the same He i transition superimposed on the emission. This absorption line forms in the outer parts of the stellar wind where He ii has recombined, and so provides an estimate of the wind velocity there. We measure the absorption line centroid at λ= µm, which implies a wind velocity of 2400 ± 100 km s 1, typical for WN5 6 stars (Crowther & Smith 1996). Furthermore, the broad emission feature near λ2.11 µm appears to extend further into the blue end of the line, when compared with the same feature in the WN7 stars X and X We attribute this to N v emission, which is a characteristic of WN4 5 stars (Crowther et al. 2006a). Therefore we relax our previous classification of WN6b to WN5 6b for X Finally, we classify X as a WN5 6b star as well, owing to the strong similarity between its spectrum and that of X Weak-lined WC9 Stars Figure 6 presents the K-band spectra of X , X and X The former two stars exhibit broad emission lines from the complex of He i, N iii, C iii and O iii (λ µm), He i and He ii (λ2.058 µm), C iii and C iv (near λ µm), and weak He ii emission near λ2.189 µm. These spectral characteristics are consistent with those of late carbon-type (WC) Wolf-Rayet stars, specifically WC9 stars, whose values of EW(λ2.08 µm)/ew(λ2.112 µm) are below 1 (e. g., see the spectrum of WR 112 in Figer, McLean & Najarro 1997). However, the peaks of the strongest emission lines are only 10% of the continuum flux. WC9 stars typically exhibit much stronger lines. However, WC stars associated with hot thermal dust emission (WC9d stars) typically exhibit relatively weak emission lines in the K-band, much like those of X and X (e. g., see Crowther et al. 2006a, their Figure 11). This is presumed to be the result of a bright thermal dust continuum which competes with the line emission of stellar wind.

11 11 The lower spectrum in Figure 6 is for X , and appears completely featureless within the noise level. Based on this star s infrared SED and hard X-ray emission, it was identified by Hyodo et al. (2008) a dust-enshrouded, late-type WC binary, also known as a DWCL star (Williams et al. 1987), the same spectral type as the Quintuplet proper members (e. g., see Figer, McLean & Morris 1999). Figure 6 provides the first near-infrared spectrum of this star, and it is consistent with the DWCL spectral type based on its featureless character. Although the weakness or absence of emission lines in the near-infrared spectra of late-type WC stars is mainly attributable to dilution by continuum emission from hot dust, line dilution may also result from the stellar continuum emission from a very luminous companion. Later, in Section 4, we will present additional evidence for the presence of hot dust associated with X and X CXOGC J X is one of the brightest emission-line stars known in the Galactic center (K s = 7.84 mag). This source was first identified as a young stellar object by Felli et al. (2002), and it was subsequently characterized as an Ofpe or LBV emission-line star by Muno et al. (2006), who also showed this star to have a 2.3 mjy radio counterpart at ν8.4 GHz. Stars of these spectral types are hydrogen rich, yet represent the beginning of post mainsequence evolution for the most massive stars (Morris et al. 1996). Figure 7 shows our K-band spectrum of X at medium resolution. Strong emission lines of Brγ and He i are accompanied by weak low-ionization emission lines of Mg ii. The λ2.058 µm line is the strongest emission feature andappears as a P-Cygni profile, while thehe i line at λ2.112µm, also exhibiting a P-Cygni profile, is much weaker and may be blended with weak emission from N iii. There also appears to be a weak P-Cygni feature of He ii near λ2.19 µm. Based on the above characteristics, X appears to be a member of Ofpe/WN9 class (e. g., see Morris et al. 1996), also resembling several stars that lie near the central parsec of the Galactic center, namely IRS16C, IRS16NW, and IRS33E, the spectra of which were modeled in Martins et al. (2007). Owing to the strong similarity of our spectra to these stars, we classify X as an Ofpe/WN9 star. From the P-Cygni profile of He i λ2.058 µm, we may infer the velocity at the radius in the wind where He i becomes an absorption line. The trough in the λ µm feature is centered at µm and the emission peak is centered at λ µm, so the implied velocity in the absorption zone is 320 km s 1.

12 12 4. Stellar Properties 4.1. Bolometric Luminosity The near-infrared photometry of each massive star is presented in Table 2, which also includes mid-infrared photometry from the Spitzer Infrared Array Camera (IRAC) pointsource catalog of the Galactic Center (Ramírez et al. 2008). To derive the absolute photometry, we must first determine the effects of interstellar extinction. Nishiyama et al. (2006) derived a near-infrared extinction relation for stars of the Galactic center using photometry from the SIRIUS survey; they report the following relations: A Ks = 1.44±0.01E H Ks and A Ks = ± 0.006E J Ks, where E(H K s ) = (H K s ) obs (H K s ) 0. For the stars presented inthis work, the observed JHK s photometry wasgenerally taken fromthesirius catalog, except for very bright stars that were saturated (K s < 8 mag), in which case we used 2MASS photometry. The intrinsic colors, (H K s ) 0 and (J K s ) 0, were taken from the literature using comparison stars of similar spectral types. For the WN stars, we used the intrinsic colors from Crowther et al. (2006a, their Table A1). For the O supergiants we used the synthetic photometry for O supergiants from Martins & Plez (2006). Using the adopted intrinsic colors, we were able derive two values of A Ks using the observed J K s and H K s colors, which were then averaged to produce a final value for the K s -band extinction. The absolute K-band magnitude of each star was derived assuming a distance of 8 kpc to the Galactic center. To calculate the total luminosity of each star, bolometric corrections of comparison stars having similar spectral types were also taken from the same literature sources referenced above. In the interest of a thorough comparison, we not only performed these calculations for new WN and O X-ray sources, but included previously known sources as well. The results are presented in Table 3. The most luminous stars in the sample are the WN8 9h stars X and A7 (Figer et al. 2002) of the Arches cluster (CXOGC J ), with L bol = L and L, respectively; the least luminous star is the broad-lined WN7 star X , which has a luminosity of L bol = L. In contrast to WNh stars which are still burning hydrogen, the WN7 stars, such as X , are in an evolved state of core-he burning. Such stars have already blown off their hydrogen envelope, and as a result, have a significantly smaller radius and lower bolometric luminosity than the less-evolved WNh sources. The situation is more difficult in the case of the WC stars; we were unable to derive an intrinsic color for these stars, owing to their significant near-infrared excess. Instead, the photometry of X was extinction-corrected by adopting the average extinction value for the Quintuplet cluster (A K = 2.7 mag; Figer et al. 1999), since X appears to be an outer member of the Quintuplet. For X , we adopted the extinction value measured for a nearby B2Ia star, S from MMM07, which has A K = 3.2 mag.

13 13 Although these extinction approximations will be useful in the analysis of the WC star SEDs, we did not compute bolometric luminosities for them, and they are not included in Table Infrared Excess To determine whether the massive stellar X-ray sources are distinguishable from Galactic center field stars in the infrared, we constructed color color diagrams. Figure 8 illustrates that almost every star exhibits a significant infrared excess, apparent as an offset from the main reddening locus that is traced by the roughly linear scatter distribution of field stars. In H K s versus J H color space, the WN stars and most of the O supergiants occupy the same space, while the WC9 and DWCL stars appear more separable as a group, owing to their significantly larger K s -band excess. With the inclusion of λ8 µm photometry, the WN, O, WC9, and DWCL stars are even easier to distinguish as separate groups in K s [8.0] versus J K s color space. The WN stars and one O supergiant appear to cluster together, while the WC stars, again, exhibit a relatively large excess compared with the WN stars, especially the DWCL stars qf211 and qf231 of the Quintuplet cluster. Three O supergiants exhibit an extremely large excess at λ8 µm, which likely originates from unresolved, warm gas and dust surrounding these young stars (evidence in support of this hypothesis provided by mid-infrared field images is presented in Section 6.2). There are several ways in which Wolf-Rayet and O stars can exhibit an infrared excess. Within the ionized winds of hot supergiants, scattering of electrons off of H + and He + ions produces free-free emission, which contributes a significant amount of continuum flux in the infrared, becoming more significant with increasing wavelength as the underlying photospheric continuum becomes fainter with respect to the free-free component (Wright & Barlow 1975). Thermal emission from hot circumstellar dust (T = K) is another source of infrared excess that is often associated with WC stars (Cohen et al. 1975, Williams et al. 1987). This form of excess emission can exhibit significantly higher flux density than free-free emission. Finally, a large infrared excess may also be detected from massive stars surrounded by unresolved clouds of gas and dust local to the star. To further investigate the nature of the infrared excesses observed in our sample, we generated λ1 10µm SEDs. For thephotometry ofthe WN andostars we applied extinction corrections to the photometry in each band according to the extinction ratios in Indebetouw et al. (2005), using the derived values of A Ks in Table 3 and the adopted values for the WC stars. To convert the mid-infrared data to flux density units, we applied zero points of 280.9, 179.7, 115.0, and Jy for IRAC channels 1, 2, 3, and 4, respectively (IRAC

14 14 Data Handbook Table 5.1). In the absence of confusion-induced photometric errors, the near-infrared SIRIUS photometry matches 2MASS very closely; the mean and rms variance of the magnitude differences are and mag in the J band, respectively (Nishiyama et al. 2008). Therefore, we obtained near-infrared fluxes for SIRIUS using the 2MASS zero points of , , and Jy for the J, H, and K s bands, respectively. Below, we examine the SEDs of several WN and WC stars to investigate the nature of the infrared excess emission in more detail Free-Free Emission from WN/O Stars To test the hypothesis that free-free emission is responsible for the infrared excess observed toward the massive stars in our sample, we employed a modified version of the free-free emission model of Cohen et al. (1975). In this model the optically-thin, free-free emission spectrum for an unresolved source may be written as F ν (Jy) = R3 s d 2 N e N i Z 2 Te e C/λTe (1) where Z is the mean electric charge per ion, N e and N i are the respective electron and ion densities in units of cm 3, T e is the electron temperature in Kelvins, R s is the size of the emitting region in cm, d is the distance to the source in cm, λ is the wavelength in µm, and C=14400 µm K. Free-free emission from evolved massive-star winds typically becomes optically thick in the mid-infrared at some wavelength, λ t. If this turnover wavelength is observed, one can also solve for the emission measure: N e N i R s = T e λ 3 tz 2 (1 e C/λtTe ) cm 5. (2) By combining Equations (1) and (2), we constructed model SEDs for several stars: X (Ofpe/WN9), X (WN5 6b), X (WN8 9h), and X (WN7 8h), all of which have an observable turnover wavelength in their mid-infrared SEDs. For each star we determined the flux density attributable to free-free emission by subtracting off a blackbody continuum. The continuum flux was calculated using the luminosity values in Table 3 and the Stefan Boltzmann relation L = 4πR 2 σt 4, where L is the bolometric luminosity, R is the stellar radius, and T is the effective temperature. Unfortunately, there is a degeneracy between temperature and radius on the Rayleigh Jeans portion of the SED, which we are limited to. So, initially we guessed values of T and R that are appropriate for

15 15 the spectral type of the star being analyzed. The values were then adjusted until we obtained a satisfactory match to the extinction-corrected near-infrared photometry. Of course, a slight underestimate of T could be compensated by an overestimate of R. This uncertainty, however, is unimportant for our purposes because the free-free emission spectrum is not affected by it. Once the resulting Rayleigh Jeans component was subtracted from the photometry, the remaining free-free flux, described by Equations (1) and (2), was fit using a least squares method. By performing a least squares fit to the excess spectrum, we were ultimately able to derive T e, N e, λ t, and R s. For the models we assume N e = N i and Z = 1 for the WNh and Ofpe/WN9 stars, and Z = 2 and N e = 2N i for the hydrogen-deficient WN5 6b star. The resulting model parameters are listed in Table 4, and are illustrated graphically in Figure 9. The model curves match the data reasonably well, and the derived values for T e, N e, and R s are similar to those of the WN stars modeled in Cohen et al. (1975). Reliable errors in our parameter fits are difficult to ascertain; uncertainties in the extinction law and the effect of strong emission lines on the photometry are likely to play a role. Nonetheless, a very precise free-free emission model of the stars in our sample is not the focus here. The goal is to unambiguously demonstrate that the infrared excess emission from the WN stars in Figure 8 is well described by a simple free-free emission model having physical parameters that are typical for such stars Hot Dust Emission from WC Stars The infrared excess emission from the WC stars in our sample (X and ) has a significantly higher intensity than that of the WN stars, which is apparent in Figure 8. Figure 10 shows the extinction-corrected SEDs of the new WC9 discoveries X and X , accompanied by the SEDs of the Quintuplet proper member X (qf231 from Figer, McLean & Morris 1999), a putatively single WC9 star HDM13 (Hadfield et al. 2007), and the WC9d star WR 59 for comparison. For the photometry of the comparison stars HDM13 and WR 59, we applied extinction corrections using the A K or A V values from Hadfield et al. (2007) and van der Hucht (2001), respectively. For the WC9 stars X and X , we were unable to use our free-free emission model to simultaneously reproduce the near-infrared and mid-infrared flux densities. No sensible combination of the stellar radii, effective temperature, and wind parameters (T e and R s ) resulted in a reproduction of the rise in flux density between the H and K s bands. Fitting the bright mid-infrared portion of the WC9 SEDs without overestimating the nearinfrared portion would require insensibly low values of T e. So, even though WC stars are

16 16 indeed important sources of infrared free-free emission, a stronger additional emission component is required to explain the structure of their infrared SEDs. Thermal emission from hot dust has been observed from many WC9 stars and binaries, and we suspect this is the case for X and X This hypothesis is strengthened by the following comparison with the other dusty WC stars in Figure 10. qf231 exhibits the largest infrared excess. This star is an X-ray source (X ), and has been classified as a DWCL star (Figer et al. 1999; Moneti et al. 2001). Such stars exhibit nearly featureless near-infrared spectra, as demonstrated in Section 3.5 (Figure 6) with the spectrum of the DWCL star X (the Spitzer/IRAC photometry of X is saturated, so its SED was not included in Figure 10). By comparison, the infrared excess from the single WC9 star HDM13 is the smallest of the stars in Figure 10, exhibiting characteristics very similar to the WN stars in Figure 9. Thus, the excess from HDM13 can perhaps be attributed to free-free emission alone. The magnitude of the infrared excess in the SEDs of X and X appears to lie in-between the extremes defined by the SEDs of qf231 and HDM13, but bear the closest resemblance to the WC9d star WR 59. Compared with qf231, the smaller infrared excess of X and X provides a natural explanation for why the emission lines in their spectra are only partially diluted, instead of being completely diluted as in the case of the DWCL spectra. Still, the mid-infrared SEDs of these sources exhibit an excess component that rises significantly farther above the RJ continuum than the free-free components of the single WC9 star HDM13 and the WN stars in the previous section. Furthermore, the excess emission peaks between λ5 and λ10 µm and is detectable at wavelengths as short as the J-band wavelength. Thus, these characteristics can be attributed to the presence of hot dust with temperatures near 1000 K. Realistically, the collective excess of these objects is probably due to multiple hot dust components lying at different radii from the central star(s), stellar continuum photons, and free-free emission. A more detailed modeling of the complex excess emission for these stars is beyond the scope of this paper. We simply intend to demonstrate that hot dust is responsible for the bulk of the mid-infrared excess detected from X and X , and provides a simple explanation for the continuum diluted emission lines observed in their K-band spectra. WC stars have long been known to be producers of dust if they exist in a binary system. In fact, it has been suggested that all dusty WC stars (WCd) are members of binaries (Williams et al. 2005), in which the dust is produced by the collision of carbonrich and hydrogen-rich winds from the WC star and an OB companion, respectively. Wind collision also naturally explains the detection of these sources at hard X-ray energies (see Section 6.1.1). For several known WCd stars, the binary hypothesis has been strengthened beyond any reasonable doubt by high-resolution imaging experiments, which have revealed rotating Archimedian spirals of dust surrounding them (e. g., seetuthill et al. 2006, 2008).

17 17 Similar structures may be observable, modulo unfortunate projection effects, in association with X , X , and X X-ray Emission Table 5 lists the X-ray photometry for all 30 confirmed massive stellar X-ray sources in the Galactic center region, including new and previously identified objects. The data for each star were compiled from the catalog of Muno et al. (2009). Table 5 includes X-ray source astrometry, total on-source integration time, photon number counts in the hard and soft bands, broadband photon flux, hardness ratios for soft and hard energy bands, the average energy of all detected photons, the derived X-ray luminosity, and the spectral types of their stellar counterparts. Most of the sources have photon fluxes ranging from F X = to F X = cm 2 s 1. An approximate energy flux can be obtained by multiplying F X by the average photon energy E, which yields values ranging from to erg cm 2 s 1. In order to constrain the nature of the X-ray emission, it is useful to examine the hardness ratio, defined as HR=(h s)/(h + s), where h and s are the fluxes in the hard and soft energy bands, respectively. The soft color, defined as HR0, is defined by letting h and s be the fluxes in the respective kev and kev energy bands. In Paper I, we used HR0 to select X-ray sources that are likely to be located near the Galactic center, since the softest X-ray photons are highly absorbed by the ISM. For this reason, many of the sources in Table 5 are not detected in the softest kev energy band, which results in them having an HR0 value of 1. The hard color, HR2, is defined by letting h and s be the fluxes in the kev and kev energy bands, respectively. Compared with HR0, the more energetic photons used to calculate HR2 suffer relatively little absorption from intervening gas and dust, so they are more useful for constraining the hard components of the X-ray SEDs. Figure 11 is a scatter plot of HR2 versus broadband ( kev) photon flux for absorbed X-ray sources that are likely to be located near the Galactic center (i.e., those sources whose soft X-ray colors imply a hydrogen absorption column of N > cm 2 ). The majority of the X-ray sources that comprise the cloud of small data points in Figure 11 are likely to be CVs, a hypothesis that is consistent with both the X-ray properties of the sample (Muno et al. 2003, 2009) and the rarity of near-infrared counterparts with K s 15.6 mag, discussed in Paper I. By comparison to the field population of presumed CVs, the confirmed massive stars appear to be systematically softer and brighter in X-rays. To estimate the plasma temperatures of thermal sources, or the power-law indices of non-

18 18 thermal sources, we constructed models of HR2 versus luminosity for optically-thin thermal plasmas with temperatures of kt = 0.9, 2.2, 3.4, and 6.8 kev, and power-law emission models for non-thermal sources with photon indices of 1 < Γ < 3. The models were constructed for sources lying at a distance of 8 kpc that are absorbed by a hydrogen column of N H = cm 2. This value of N H is generally adopted as a reasonable average value for X-ray sources near the Galactic center (Muno et al. 2003), and is also consistent with the average extinction of the stars listed in Table 3, according to the relation between N H and A V given in Predehl & Schmitt (1995). All models were computed using the XSPEC 1 program. The HR2 values of the massive stars imply the presence of plasma components with temperatures of kt > 1 kev (in most cases kt > 2 kev), if the emission is thermal, and luminosities in the range L X erg s 1 ( kev). The X-ray sources with the highest fluxes and hardest colors are the DWCL star X (Hyodo et al. 2008) and several of the Arches cluster members, all of which appear to have relatively high luminosities near L X erg s 1. The softest sources, undetected in the kev band, are the O9I B0I star X and the Quintuplet members qf231(dwcl star) and qf242(obi). The hardest massive stellar X-ray source in the Galactic center is the WN7 star X , with HR2=0.3, and an average photon energy of 5.3 kev. We note that a constant value of N H is a gross simplification; more so is the assumption of a single temperature component for the X-ray plasma. Below, we address the effect that these simplifications have on the derived properties of the X-ray-emitting plasma. The fiducial conversion from photon flux to unabsorbed kev luminosity is such that erg s 1 equals cm 2 s 1, assuming Γ = 1.5 or kt = 7 kev, which are average values for the entire X-ray population. Assuming kt = 2 kev instead, which is more applicable to our sample of massive stellar X-ray sources, the conversion is such that erg s 1 equals cm 2 s 1. However, the uncertainty in these conversion factors is substantial. Indeed, interstellar absorption has extinguished most of the X-ray flux at energies below 1 kev, which is where massive stars emit the bulk of their X-rays (e.g., see Oskinova 2005, Skinner et al. 2002a), so we could be significantly underestimating the total intrinsic X-ray luminosities in the kev range. As a particularly fitting example of the errors imposed by assuming a single temperature plasma, we refer to the X-ray analysis of the source X by Mikles et al. (2006). Owing to relatively bright soft and hard X-ray fluxes, the authors were able to fit a detailed thermal plasma model to the X-ray spectrum of this source. Approximating with a single temperature plasma, as we do for our sample as a whole, the authors derived kt 3.5 kev and L X = erg s 1 ( kev) for X , which is reasonably close to the plasma temperature implied by the 1

19 19 HR2 value of this source, and to our derived luminosity of erg s 1. However, the authors obtained a more accurate spectral fit using a two-temperature model, having a very soft component kt 0.7 kev and an additional hard component with kt 4.6 kev, each with respective luminosities of L X = erg s 1 ( kev) and erg s 1 ( kev). Thus, the resulting X-ray luminosity difference between the 1-T and 2-T plasma models is a sobering 1.4 dex. Clearly, in assigning a single temperature plasma we may be neglecting a significant amount of soft X-ray luminosity, if it lies beneath our flux threshold due to absorption. Unfortunately, only two of the new detections reported in this work were detected at energies below 2 kev, the brightest of which has only 13 photon counts in that range. Without a significant flux of soft X-rays to subject to a model fit, assuming a single temperature plasma is the best we can currently do. So, we assumed that 2 kev is a reasonable average intrinsic plasma energy and based our luminosity calculations on this assumption. The resulting X-ray luminosities are included in Column 12 of Table 5. As a means of estimating the errors, we compared our values to the results of other authors who were able to derive luminosities for relatively bright, Galactic center X-ray sources based on a modeling of their X-ray spectra, which enable them to solve for kt directly. We did not use the above-mentioned X in this exercise because it is a variable X-ray source (see Section 5.1). The comparison is facilitated by Table 6, which includes the three bright WN8 9h X-ray sources from the Arches cluster (Wang, Dong & Lang 2006), the DWCL star X (Hyodo et al. 2008), the WN7 star X , and the Of star X (Muno et al. 2006). The average and standard deviation of the difference in the X-ray luminosity between our results and those of the other authors is 0.25 (0.24) dex. Thus, we assign 3σ errors of 0.75 dex for our derived luminosities. The results should still be met with caution, however, since the comparison measurements from the other authors listed in Table 6 also assumed a single temperature plasmas. Finally, using our derived X-ray luminosities and the stellar bolometric luminosities in Table 3, we recovered a relation between these values, such that L X /L bol = ±0.6. Figure 12 illustrates this result. In conclusion, the main points we intend to convey with X-ray data are that (1) the X-ray photometry of the massive stars as a group is systematically softer and brighter than the field population presumed to be dominated by CVs, (2) each source contains a hard X-ray component, requiring either a thermal plasma with kt > 1 kev (in most cases kt > 2 kev), or a power-law source with 1 < Γ < 3, (3) the majority of sources have kev luminosities between L X and erg s 1, while several others appear to have higher luminosities closer to erg s 1, and (4) the L X /L bol ratio for the Galactic center sources has an approximately constant value of ±0.6.

20 X-ray Variability Variability of the Galactic center X-ray population was addressed in Muno et al. (2009). Variable sources were grouped according to three types of variability: flux variations between individual observations, separated by day- to year-long time scales (long-term variability), variations that occurred within individual observations (short-term variability), and periodic variations that occurred within an individual observation. Four of the massive stellar X-ray sources near the Galactic center exhibit long-term variability. None of them exhibits shortterm or periodic variability. The variable sources include the WN7-8h stars X and X , the Ofpe/WN9 star X , and the OBI Quintuplet member X We reproduce the variability data for these sources from Muno et al. (2008, their Table 6) in our Table 7. The table includes the source record numbers from the original X-ray catalog, the observation ID numbers in which the largest and smallest fluxes were observed, the values of the largest and smallest fluxes, and the ratios of those fluxes. The most variable source among the massive stars is the WN7-8h star X , which varied by a factor of 3.6 between 2000 March 29 and 2001 July 16. A variety of phenomena may give rise to X-ray variability. Colliding-wind binaries with eccentric orbits will increase in X-ray luminosity near periastron. Eclipses of a wind collision zone may also result in variable fluxes if the orbit is close to edge-on. However, colliding-wind binaries may not exhibit detectable X-ray variations, as in the case of WR 25 (Pollock &Corcoran2005), which remained stableforatimespanoftenyears(raassenet al. 2003). For accreting compact objects in HMXBs, eclipses are also a possible source of X-ray variability, as is an intrinsic variation of accretion rate. Finally, X-rays from hot spots within an equatorial disk surrounding a massive star (see Section 6.1) may be modulated with the stellar rotation period. The applicability of these models in explaining the characteristics of the X-ray emission of our sample is discussed in the following section. 6. Discussion 6.1. Physical Origin of the X-ray Emission The X-ray photometry of the majority of massive stellar X-ray sources near the Galactic center indicates that the emission is from thermal plasmas having temperatures of kt 1 2 kev, or non-thermal sources with power-law indices of 1 < Γ < 3, andx-ray luminosities in the range L X = erg s 1 ( kev). Hard-X-ray emission at these luminosities is not a ubiquitous feature of single O supergiants or Wolf-Rayet stars; these features are typically regarded as an indication of binarity, since it is well established that hard X-rays

21 21 can be generated within the shocked interface of the binary components opposing supersonic winds, or by accretion of stellar wind material onto a compact companion. However, there are extraordinary mechanisms by which some single massive stars can intrinsically generate hard X-rays as well; so, we entertained several interpretations of the hard X-ray emission from our sample Intrinsic X-ray Emission from Massive Stars The intrinsic X-ray emission from single massive stars throughout the Galaxy is typically thermal and soft, having kt 0.6 kev (e.g., see Oskinova 2005 and references therein). The X-ray luminosities of O stars follow an observed trend in which they scale with stellar bolometric luminosity as L X /L bol 10 7 (Long & White 1980; Berghoefer et al. 1997; Sana et al. 2006; Broos et al. 2007, Albacete Colombo et al. 2007). This trend has been interpreted as a manifestation of the physical link between stellar luminosity, which drives supersonic winds, and the X-ray-emitting shocks that result from those winds via line-driven instabilities (Lucy & White 1980; Owocki et al. 1988; Feldmeier, Puls, & Pauldrach 1997). In Section 5, we derived the average and standard deviation of the bolometric to X-ray luminosity ratio for our sample, and obtained L X /L bol = ±0.6, which is consistent with the known trend for single massive stars. However, the fact that the L X /L bol relation is obeyed does not necessarily imply that the X-rays are solely generated via radiative shocks within a single stellar wind; colliding-wind binaries have also been shown to obey the L X /L bol relation as well, although many colliding-wind binaries have L X /L bol closer to 10 6 (e.g., see Oskinova 2005). Still, the emission of hard X-rays is uncommon property of single massive stars. There is, however, an extraordinary mechanism by which massive stars can genrerate hard X-rays intrinsically. The extraordinary O5.5V star θ 1 Orionis C, for example, has an average X-rayluminosity ofl X ergs 1 ( keV), andahardspectral component from thermal plasma with kt 3 kev (Gagné et al. 2005). The hard X-ray emission is modulated with the star s 15 day rotation period; so, the hard X-rays have been interpreted as a result of the star s strong, fossil magnetic field (B 1000 G, Donati et al. 2002), which channels and confines the outflowing stellar wind into an equatorial collision/cooling disk (the magnetic wind-confinement model). The process results in strong shocks that generate plasma temperatures in excess of 10 7 K, and hence, hard X-rays. The peculiar Of?p star HD is another example, having a measured magnetic field of 1500 G (Donati et al. 2006) and an X-ray spectrum containing thermal plasma components with kt kev and an X-ray luminosity of erg s 1 (Nazé et al. 2007). HD exhibits periodic spectroscopic variations (a time-variable equivalent width of the Hα line) with a period of 538 days. Walborn et al. (2004) discussed the possibility that the variability is induced

22 22 by the periastron passage of a binary companion, while Donati et al. (2006) attributed the variations to the star s rotational period and, hence, its strong magnetic field. Since θ 1 Orionis C and HD both exhibit periodic changes in their spectra, we should consider the potential to observe the same phenomena in the spectra of the O supergiants of our sample, if they represent similar stellar phenomena. The magnetic confinement model has been discussed in attempts to explain the detection of hard X-rays from the nitrogen-rich Wolf-Rayet stars as well, in particular, for WR110, EZ CMa, and the oxygen-rich Wolf-Rayet star WR142, all of which have stellar and X-ray properties that are very similar to our sample and exhibit no evidence for a companion (Skinner et al. 2002a, 2002b; Oskinova et al. 2009). However, the magnetic wind-confinement model faces difficulties when applied to Wolf-Rayet stars; the wind momenta of Wolf-Rayet stars are over an order-of-magnitude stronger than that of O stars such as θ 1 Orionis C and HD , and are more difficult to channel effectively. According to Babel & Montmerle (1997), a stellar wind is confineable if the magnetic energy density exceeds the wind kinetic energy density (B 2 /4π > ρv 2 ). Assuming ρ = Ṁ/4πvr2, and v = 0.5v at r = 1R star, this would imply that magnetic fields greater than 5 kg in strength are necessary to confine the winds of the Wolf-Rayet stars and O supergiants that dominate our sample. Currently, there is no evidence for such strong fields in Wolf-Rayet stars. Strong fields have been considered in the interpretation of hard X-ray emission from the WO star WR 142 (Oskinova et al. 2009); however, this source is intrinsically 10 times fainter in X-rays than the faintest source in our sample. Thus, we do not currently regard the magnetic wind-confinement scenario as a highly viable model for Galactic center Wolf-Rayet stars, although it may explain the X-ray emission from some of the O stars in our sample. Still, we can not completely rule out the possibility that Wolf-Rayet stars can generate hard X-rays intrinsically by some other exotic mechanism, such as inverse Compton scattering (e.g., see Chen & White 1991). However, since there are so many Wolf-Rayet and O stars known in the Galactic center that have not been detected in X-rays (e.g., in the Arches and Quintuplet clusters), intrinsic hard X-ray generation from our sample, which represents the minority fraction of known massive stars in the Galactic center, would imply that our particular stars are special, having physical properties not ubiquitous to all Wolf-Rayet stars and O supergiants in the region, which would be puzzling Colliding-Wind Binaries The collision of opposing, supersonic winds in massive binaries creates a hot shock (T 10 7 K) that is a potential source of hard thermal X-rays with kt > 1 kev. It is

23 23 reasonable to adopt this model in interpreting the hard X-ray emission from our sample, since many known colliding-wind binaries exhibit hard X-rays with luminosities that are comparable to our Galactic center sample, such as WR 147 (WN8+OB; Skinner et al. 2006), WR 25 (WN6ha+O4f ; Raassen et al. 2003; Albacete Colombo et al. 2008), γ 2 Velorum (WC8+O7.5; Skinner et al. 2001), V444 Cyg (WN5+O6; Maeda et al. 1999), Cyg OB2 8 (O6If+O5.5III(f); Albacete Colombo et al. 2007), and Plaskett s star (O6I+O7.5I; Linder et al. 2006), all of which have hard spectral components with kt 2 4 kev and L X = erg s 1. With regard to the WC stars of our sample (X , X , and X ), the case for colliding winds is particularly strong, since most single WC stars are not detected in X-rays at all (Skinner et al. 2006), and hot dust is present (see Section 4.3). It has been proposed that the production of dust from carbon-rich WC winds requires collision and mixing with a hydrogen-rich companion star s wind (Crowther 2003). While the explanation for hard X-ray emission in the WN/O stars of our sample is not as clear, the colliding-wind hypothesis may be strengthened by the fact that binarity is common in Wolf-Rayet and O stars (Wallace 2007). If the colliding-wind binary hypothesis is correct, one might expect the signature of a companion star to be evident in the infrared spectra, presented in Section 3. However, the current infrared data are probably of insufficient spectral resolution and signal to noise to distinguish the contributions from individual stellar components of a massive binary. In the case where a Wolf-Rayet star or O supergiant has a fainter OB dwarf companion, the magnitudes may differ by a factor of 3. With the implied flux ratio of 16 for the stellar components, HandHeemissionorabsorptionlinesfromtheOBcompanionwillbedifficultto detect in a moderate-resolution spectrum. In the case where the binary consists of two similar Wolf-Rayet stars of comparable brightness, the broad emission lines of both components may blend together. In this case, the spectral lines of individual binary components might only be distinguishable with multi-epoch observations, where radial velocity shifts could betray their presence. Thus, the colliding-wind binary scenario is consistent with the current infrared spectra. Assuming our sample is mainly represented by binaries, it is worthwhile to speculate whether their global X-ray properties can be used to constrain the physical parameters of the binary system. For colliding-wind binaries, the emergent X-ray luminosity from the shocked wind-wind interface is a function of the wind velocities and mass-loss rates of the stellar components, and the binary separation. The exact functional form depends on the nature of the wind collision. For instance, if the wind momenta of two stars in a binary are comparable, then the shock zone will lie somewhere near the center of the space between them, where the opposing winds achieve momentum balance. The collision interface will form a bow shock structure that curves toward the star with the weaker wind. According

24 24 to Usov (1992; see their Equations (89) and (95)), the total X-ray luminosity in this case can be expressed as the sum of the contributions from two separate layers, one that lies on the inside (concave) surface of the bow shock, composed of material from the wind of the weaker star, and one that lies on the outside (convex) surface of the bow shock, composed primarily of material from the wind of the stronger star. Using mass-loss rates and terminal wind velocities that are typical for the Wolf-Rayet and O stars of our sample (Ṁ M yr 1, v = km s 1 ; see Section 3), this particular model predicts X-ray luminosities of L X ergs s 1 for binary separations of D = AU. Given our uncertainty in X-ray luminosity, the brightest several sources in our sample could be binaries in which the components have comparable wind momenta. However, the majority of our sources have fainter X-ray luminosities in the range L X erg s 1. The emergent X-ray luminosity in this particular colliding-wind model has a dependence of L X D 1, so, this scenario could result in X-ray luminosities consistent with those of our sample only if the binary separations are very large (D AU). This would make the Galactic center sources analogous to the WN8h+OB binary WR 147, which has an estimated binary separation of D 400 AU, and an X-ray luminosity and plasma temperature of L X (0.5 7 kev) erg s 1 and kt = 2 4 kev, respectively (Skinner et al. 2007). With periods of 1000 years, binaries like these would not exhibit detectable radial velocity variations. However, at a distance of 8 kpc, they could have apparent separations of 50 mas, which could be resolvable with adaptive optics on 8 10 m class telescopes. Alternatively, Usov (1992) also considers the case where the ram pressure of the stronger stellar wind suppresses the weaker wind of the companion (see their Equation (81)), in which case the companion star is modeled as a hard sphere and the wind collision occurs very near or on its surface. For this particular model, which has a much stronger dependence on the binary separation (D 4 ), the expected X-ray emissivity of the shock zone yields lower values of X-ray luminosity in the range of L X ergs s 1. Thus, if this particular model is applicable, then perhaps the Galactic center sources are binaries where the components have significantly different wind strengths. This model has been considered to explain the X-ray emission ofthewn5 6star WR 110(Skinner et al. 2002a), which hasl X (0.5 7keV) erg s 1 and a hard spectral component with kt 3 kev. It is interesting that the most X-ray-luminous stars in our sample, with L X erg s 1, are mainly comprised of objects within the Arches cluster. These particular stars could be binaries that have relatively small separations. Frequent stellar encounters within dense stellar clusters will tend to harden close binaries while ionizing wide ones; so, it is reasonable to speculate that the relatively high X-ray luminosity of the Arches X-ray sources is the result of their relatively small orbital separations that were imposed by the dense cluster environment. If so, radial velocity variations will occur on shorter timescales

25 25 than for the rest of our sample and may exhibit a larger amplitude in the radial velocity variations, similar to that of the WN6h+O binary NGC 3603-A1, which lies near the center of the starburst cluster NGC 3603 and exhibits radial velocity amplitudes of 400 km s 1 (Schnurr et al. 2008) Supergiant High-Mass X-Ray Binaries We must also consider the possibility that wind-accreting neutron stars or black holes are responsible for the X-ray emission from some of the massive stellar X-ray sources near the Galactic center, making them supergiant X-ray binaries (SGXBs). Although there is no reason to expect such objects to follow the observed trend of L X /L bol = 10 7, there are several quiescent SGXBs whose L X /L bol ratios that are consistent with this trend. So, the wind-accretion model still deserves consideration. According to Pfahl, Rappaport & Podsiadlowski (2002), the time-averaged 1 10 kev accretion luminosity of a wind-accreting neutron star is given by ( ) L X erg s η Ṁ ( w a ) ( ) 2 4 v w M yr AU 1000 km s 1 (1 e 2 ) 1/2, (3) where η is the conversion efficiency of gravitational energy into radiation, Ṁ w and v w are the mass-loss rate and wind velocity of the donor star, respectively, a is the orbital separation in AU, and e is the orbital eccentricity of the system. This formula is based on the Bondi Hoyle Lyttleton accretion model (Hoyle & Lyttleton 1941; Bondi & Hoyle 1944). For the massive stars of our sample, most of which have Ṁ M yr 1 and v 2000 km s 1, Equation (3) implies that the emergent X-ray luminosity from a wind-accreting SGXB with these wind parameters is L X erg s 1 for binary separations ranging from 0.5 AU to several AU, assuming a radiative efficiency of 1. Currently, there are very few physical examples of SGXBs with X-ray luminosities of L X erg s 1 or lower, presumably because X-ray surveys of the Galactic plane are biased toward bright SGXBs (e. g., Roche-lobe overflow systems with L X erg s 1 ; Liu et al. 2000). However, a new class of SGXB has recently been unveiled in obscured regions of the Galactic plane by hard X-ray surveys, which may provide some insight into X-ray production in low-accretion-rate systems: the superfast X-ray transients (SFXT; Negueruela et al. 2006a; Sguera et al. 2006) are SGXBs containing hot supergiant primary stars. These systems undergo very brief X-ray outbursts, in which the X-ray flux increases by a factor of for durations on the order of hours, in-between much longer states of quiescence (L X erg s 1 ) that last on the order of weeks to a month. Approximately 10 such

26 26 objects are currently known, many of which contain X-ray pulsars. The outbursts are thought to occur when a neutron star in an eccentric orbit nears periastron, and plunges into dense inhomogeneities within the hot supergiant s wind. The long quiescent periods ensue as the neutron star travels further from the primary; it is this particular state of SFXTs which we intend to compare with our sample of Galactic center X-ray sources. For example, during their quiescence states the SFXTs XTE J and IGR J exhibit hard power-law spectra having Γ 1 2 and X-ray luminosities of L X (3 4) erg s 1 within the kev energy range (Sidoli et al. 2008; Romano et al. 2007, 2009). If these particular values of Γ and L X for XTE J and IGR J were plotted in Figure 11, they would be practically indistinguishable from the stars of our sample. Thus, the X-ray properties of the Galactic center sources are consistent with the quiescent states of SFXTs. Furthermore, quiescent X-ray emission from an SFXT is also consistent with Equation (3): IGR J , which contains the B1Ia star HD , exhibits evidence for an orbital period of P 165 days (Romano et al. 2009). B1Ia stars have M M, Ṁ w 10 6 M yr 1, and v w 1000 km s 1 (Crowther et al. 2006a, 2006b). So, Kepler s law implies that this system should have a separation of a 2 AU. Thus, Equation (3) implies an X-ray luminosity of L X erg s 1 for this system, which is a close approximation to this source s lowest observed quiescent luminosity of L X = erg s 1 (Romano et al. 2007). However, although there has been moderate X-ray variability observed for the several sources listed in Table 7, there is no evidence that any of the stars of our sample underwent eruptions or outbursts that are on the level known to occur in SFXTs. Thus, we are not suggesting that the stars in our sample are SFXTs; rather, we are examining the X-ray properties of quiescent SFXTs to provide some insight into the nature of faint X-ray emission from SGXBs with relatively wide binary separations. Indeed, SGXBs with wider and more circular orbits could exhibit the quiescent X-ray properties of SFXTs without the periodic outbursts. In conclusion, the current data are consistent with both the colliding-wind and windaccretion hypotheses. The most effective method of discriminating between these two possibilities is to derive a constraint on the system mass functions. Colliding-wind binaries will exhibit variable doppler shifts in their spectra on the order of km s 1, modulo unfavorable binary inclinations or large binary separations, while the doppler signature of the gravitational tug of a wind-accreting neutron star (M 2M ) in the spectrum of a supergiant donor star would be more challenging to detect. We plan to perform such experiments in the near future.

27 Spatial Distribution, Environments, and Implications for Star Formation Mid-infrared images from Spitzer/IRAC of the fields containing the massive stellar X- ray sources are presented in Figures For Figures 13 and 14, three-color images were assembled using data from IRAC channels 1 (λ3.6 µm), 3 (λ5.8 µm), and 4 (λ8.0 µm), and are displayed with a histogram-equalization contrast stretch to enhance subtle structure in surrounding nebulosity. Each image is 144 on a side, which corresponds to a projected physical distance of 5.6 pc, assuming a distance of 8 kpc to the Galactic center. The images bear the name of the X-ray source that lies at the field center as well as alphabetical labels that show the location of the field in the wide-field λ8.0 µm image presented in Figure 15. Several fields in Figures contain multiple massive stellar X-ray sources or additional, known massive stars: field D of the WC9d star X also contains the Ofpe/WN9 star X onthewestern side of theimageandtheb2 Ia star frommmm07 ontheeastern side; field F of the O4 6I star X (the central source of the H1 nebula) contains the Of star X on the western side of the image and the WN7 8h star X (Mikles et al. 2006) on the eastern side; field G of the O9I B0I star X contains the OBe star H5 (Cotera et al. 1999) to the northeast of the field center and the WC8 9 star WR 101p (Homeier et al. 2003) to the far northeast; finally, field K of the WC9d star X also contains the central portion of the Quintuplet cluster toward the north-west of the image, which is dominated by the infrared emission of very bright DWCL Quintuplet proper members, including X (qf231) and X (qf211). The LBV Pistol Star and its spherical ejection nebula are also apparent in field K. Several massive stellar X-ray sources appear to have diffuse mid-infrared nebulae associated with them. In particular, each of the three O supergiants that exhibit the largest λ8.0 µm excess in Figure 8 (X , X , and X ) appears to be surrounded by such nebulae, which is extended in the case of , and more compact in X , and X The excess point-source emission may be natal gas and dust in the process of being heated and cleared away by the strong radiation fields and winds of the supergiants. Although we did not present the SEDs of these stars in Section 5, we did examine them, and determined that they appear very similar to the bow-shock source IRS8 near the central parsec, which although unresolved by Spitzer, was clearly resolved by Geballe et al. (2006) via adaptive-optics imaging. Hence, we speculate that the bright λ8.0 µm excess from X , X and X may be unresolved dust emission from bow-shocks or heated dust structures that are similar to that surrounding IRS8. The WN7-8h stars X and X also appear to be associated with bright mid-infrared nebulae, although in these cases the nebulae are well resolved from the central stellar point

28 28 source. The O9I B0I star X and the WN7 8h star X also appear to be having strong impact on their surroundings, as evidenced by the bow-shock-like or shell-like morphology of nearby nebulae. The origin of the massive stellar X-ray sources near the Galactic center is currently unclear. With the exception of several rather obvious associations with the Arches and Quintuplet clusters, it is possible that the relatively isolated stars near the clusters (e. g., X and X ) originated within them but were dynamically ejected via gravitational interactions with other massive stars and binaries within the cluster, or were ejected as the result of a companion supernova. Indeed, a noteworthy fraction of O stars ( 10% 30%) in the Galactic disk are known to be runaways, 10% or more of which are binaries (Gies & Bolton 1986; Gies 1987). Runaway Wolf-Rayet+O binaries are also known to exist, such as WR 21 and WR 22 in Carina (Moffat et al. 1998). Thus, it is likely that some massive stars and binaries have been dynamically ejected from the Arches and Quintuplet clusters over their respective lifetimes of 1 2 Myr (Figer et al. 2002) and 4 6 Myr (Figer, McLean & Morris 1999). It is also tempting to speculate on the ejection scenario for the DWCL star X (Figure 15, subfield M), owing to its apparent solitude and relatively high latitude above the Galactic plane, compared with the other massive stellar X-ray sources. Unfortunately, there is no obvious potential point of origin for this evolved star, so its birthplace remains elusive. Proper motion measurements performed with adaptive optics, and high-resolution spectra, would be required to obtain the kinematic information necessary to place meaningful constraints on the origins of the isolated massive stars. Alternatively, some of the stars from our sample may be the products of a mode of isolated massive star formation, operating in tandem with the formation of dense stellar clusters. The relatively loose aggregations of early-o and WNh stars near the H1 H8 H ii regions(figure13,fieldseandf)andsagittariusb,andthemoreevolvedwnandwcstars southwest of Sagittarius A West might be the products of such a formation mode. However, it is also possible that some of these stars are members of stellar clusters less extreme in mass and density than the Arches and Quintuplet, which have gone unnoticed owing to confusion with the dense stellar background. Indeed, clusters near the Galactic center are subject to strong tidal forces, causing them to dissolve on time scales of 10 Myr (Kim et al. 1999), while they may fade beyond detectability as a surface density enhancement within only a few Myr (Portegies-Zwart et al. 2001). If such hidden clusters exist, narrow-band imaging in the vicinity of the newly discovered stars of our sample could betray their presence with the detection of additional, coeval emission-line stars. Although the origin of the massive stellar X-ray sources reported in this work is puzzling, it is clear that significant massive star formation has been occurring in the Central Molecular

29 29 Zone, outside of the well known clusters, for several Myr, including the formation products of the Sagittarius B molecular cloud complex and the newly discovered regions southwest of Sagittarius A West. The Chandra observations of the Arches and Quintuplet clusters indicate that only a small fraction of massive stars in the Galactic center are detectable as hard X-ray sources. Indeed, the Arches and Quintuplet contain over 100 massive stars each (M 20M ), yet only several ( 5%) are detected by Chandra as X-ray sources. If roughly the same fraction of massive stars is detectable in X-rays throughout the entire Central Molecular Zone, we speculate that the existence of 20 massive stellar X-ray sources outside of the stellar clusters implies the presence of 400 additional massive stars that have yet to be identified. It is noteworthy that the total inferred Lyman-α photon production rate that is required to explain the far-ir luminosity emerging from the central half-kiloparsec ( photons s 1 ) is about twice the amount generated by the three clusters collectively (Cox & Laureijs 1989, Figer et al. 2004); the several hundred massive stars implied to exist by the results of this work would account for the current deficit. However, dense clusters may dynamically produce more hard X-ray sources than less dense associations, in which case the number of massive X-ray-emitting stars detected outside of the clusters would imply an even larger population of undetected massive stars. 7. Summary and Concluding Remarks We have confirmed 16 new massive stars in the Galactic center through infrared spectroscopy of near-infrared counterparts to Chandra X-ray sources. The spectral types include early- and late-o supergiants, as well as WN and WC Wolf-Rayet stars. All of the stars exhibit infrared excess, which is attributable to free-free emission from ionized stellar winds, supplemented by thermal emission from hot dust in the case of several WC stars. For several cases we modeled the free-free emission and derived stellar wind parameters, producing results that are in excellent agreement with previous studies of similar stars. The X-ray photometry of the sources exhibit significant flux at energies above E > 2 kev, which we regard as an indication of binarity since this property is not typical of single massive stars. Our favor currently leans toward the colliding-wind binary hypothesis for the majority of our sample, mostly owing to similarities in their X-ray properties with known objects from the literature, and the commonality of such objects among massive-star associations in the Galaxy. However, the current data are also consistent with the presence of accretion-powered SGXBs in our sample. It should be noted, however, that without direct evidence for companions, intrinsic hard X-ray generation from single stars, either via magnetically confied winds or some other mechanism, cannot be completely ruled out. Only

30 30 a constraint on system mass functions via photometric and radial velocity monitoring will allow further discrimination between these alternate hypotheses, and we plan to conduct such experiments in the near future. The newly discovered Wolf-Rayet stars and O supergiants bring the total number of massive stellar X-ray sources to 31 (30 if we regard the X-ray detection of X as spurious; see Section 3.1), including those within the Arches and Quintuplet clusters. The detectionof 20suchsourcesoutsideoftheclusterssuggeststheremaybealargepopulation of several hundred massive stars that are currently unidentified in the Galactic center region, since we expect no more than 5% of massive stars in the Galactic center to be detectable in X-rays within our current sensitivity limit. If such a large population exists, it would account for the missing Lyman-α photon flux required to explain the far-ir luminosity emerging from the central half-kiloparsec. Finally, in Paper I we estimated that there are likely to be absorbed X-ray sources having real infrared counterparts with K s < 15.6 mag. Expanding the spectroscopic search to include fainter infrared matches to X-ray sources will undoubtedly reveal more massive stellar counterparts in the Galactic center, including additional Wolf-Rayet and O stars, and HMXBs. Such discoveries will elucidate the origin of these stars and their overall distribution throughout the Galaxy s Central Molecular Zone. The next generation of near-infrared multi-object spectrographs on 8 10 m telescopes, such as FLAMINGOS2 on Gemini-South (Eikenberry et al. 2008) and MOSFIRE on Keck (McLean et al. 2008), will be well suited for this task. This research was based upon observations with the Chandra X-ray Observatory. Support was provided by the National Science Foundation (grant AST ). J. M. thanks Fabrice Martins for providing model stellar spectra of the Arches stars for comparison, and comments on the spectral classification of the massive stars presented in this work. We also thank Shogo Nishiyama and Tetsuya Nagata for providing near-infrared photometry from the SIRIUS survey.

31 31 REFERENCES Agol, E., & Kamionkowski, M. 2002, MNRAS, 334, 553 Albacete Colombo, J. F., Damiani, F., Micela, G., Sciortino, S., & Harnden, F. R., Jr. 2008, A&A, 490, 1055 Albacete Colombo, J. F., Flaccomio, E., Micela, G., Sciortino, S., & Damiani, F. 2007, A&A, 464, 211 Arendt, R. G., et al. 2008, ApJ, 682, 384 Babel, J., & Montmerle, T. 1997, A&A, 323, 121 Berghoefer, T. W., Schmitt, J. H. M. M., Danner, R., & Cassinelli, J. P. 1997, A&A, 322, 167 Bohannan, B., & Crowther, P. A. 1999, ApJ, 511, 374 Bondi, H., & Hoyle, F. 1944, MNRAS, 104, 273 Broos, P. S., Feigelson, E. D., Townsley, L. K., Getman, K. V., Wang, J., Garmire, G. P., Jiang, Z., & Tsuboi, Y. 2007, ApJS, 169, 353 Chen, W., & White, R. L. 1991, ApJ, 366, 512 Clark, J. S., Crowther, P. A., & Mikles, V. J. 2009, arxiv: Cohen, M., Hammersley, P. L., & Egan, M. P. 2000, AJ, 120, 3362 Cohen, M., Kuhi, L. V., & Barlow, M. J. 1975, A&A, 40, 291 Conti, P. S., Hanson, M. M., Morris, P. W., Willis, A. J., & Fossey, S. J. 1995, ApJ, 445, L35 Cotera, A. S., Erickson, E. F., Colgan, S. W. J., Simpson, J. P., Allen, D. A., & Burton, M. G. 1996, ApJ, 461, 750 Cotera, A. S., Simpson, J. P., Erickson, E. F., Colgan, S. W. J., Burton, M. G., & Allen, D. A. 1999, ApJ, 510, 747 Cox, P., & Laureijs, R. 1989, in IAU Symp. 136, The Center of the Galaxy, (Dordrecht: Kluwer), Ed. Mark R. Morris, p121 Crowther, P. A. 2003, Ap&SS, 285, 677

32 32 Crowther, P.A., Hadfield, L.J., Clark, J.S., Negueruela, I., &Vacca, W. D.2006a, MNRAS, 372, 1407 Crowther, P. A., Lennon, D. J., & Walborn, N. R. 2006b, A&A, 446, 279 Crowther, P. A., & Smith, L. J. 1996, A&A, 305, 541 Cushing, M. C., Vacca, W. D., & Rayner, J. T. 2004, PASP, 116, 362 Cutri, R. M. et al. 2003, Explanatory Supplement to the 2MASS All Sky Data Release (Pasadena, CA: IPAC/Caltech) ibitem[depoy et al.(1993)]1993spie d Depoy, D. L., Atwood, B., Byard, P. L., Frogel, J., & O Brien, T. P. 1993, Proc. SPIE, 1946, 667 de Pree, C. G., Goss, W. M., & Gaume, R. A. 1998, ApJ, 500, 847 Donati, J.-F., Babel, J., Harries, T. J., Howarth, I. D., Petit, P., & Semel, M. 2002, MNRAS, 333, 55 Donati, J.-F., Howarth, I. D., Bouret, J.-C., Petit, P., Catala, C., & Landstreet, J. 2006, MNRAS, 365, L6 Eikenberry, S., et al. 2008, Proc. SPIE, 7014, 28 Feldmeier, A., Puls, J., & Pauldrach, A. W. A. 1997, A&A, 322, 878 Felli, M., Testi, L., Schuller, F., & Omont, A. 2002, A&A, 392, 971 Figer, D. F., Kim, S. S., Morris, M., Serabyn, E., Rich, R. M., & McLean, I. A. 1999, ApJ, 525, 750 Figer, D. F., McLean, I. S., & Najarro, F. 1997, ApJ, 486, 420 Figer, D. F., Rich, R. M., Kim, S. S., Morris, M., & Serabyn, E. 2004, ApJ, 601, 319 Figer, D. F. et al. 2002, ApJ, 581, 258 Gagné, M., Oksala, M. E., Cohen, D. H., Tonnesen, S. K., ud-doula, A., Owocki, S. P., Townsend, R. H. D., & MacFarlane, J. J. 2005, ApJ, 628, 986 Geballe, T. R., Najarro, F., Rigaut, F., & Roy, J.-R. 2006, ApJ, 652, 370 Gies, D. R. 1987, ApJS, 64, 545 Gies, D. R., & Bolton, C. T. 1986, ApJS, 61, 419

33 33 Hadfield, L. J., van Dyk, S. D., Morris, P. W., Smith, J. D., Marston, A. P., & Peterson, D. E. 2007, MNRAS, 376, 248 Hanson, M. M., Conti, P. S., & Rieke, M. J. 1996, ApJS, 107, 281 Hanson, M. M., Kudritzki, R.-P., Kenworthy, M. A., Puls, J., & Tokunaga, A. T. 2005, ApJS, 161, 154 Hillier, D. J., & Miller, D. L. 1998, ApJ, 496, 407 Homeier, N. L., Blum, R. D., Pasquali, A., Conti, P. S., & Damineli, A. 2003, A&A, 408, 153 Hoyle, F., & Lyttleton, R. A. 1941, MNRAS, 101, 227 Hyodo, Y., Tsujimoto, M., Koyama, K., Nishiyama, S., Nagata, T., Sakon, I., Murakami, H., & Matsumoto, H. 2008, PASJ, 60, 173 Indebetouw, R., et al. 2005, ApJ, 619, 931 Kim, S. S., Morris, M., & Lee, H. M. 1999, ApJ, 525, 228 Krabbe, A., et al. 1995, ApJ, 447, L95 Law, C., & Yusef-Zadeh, F. 2004, ApJ, 611, 858 Laycock, S., Grindlay, J., van den Berg, M., Zhao, P., Hong, J., Koenig, X., Schlegel, E. M., & Persson, S. E. 2005, ApJ, 634, L53 Liermann, A., Hamann, W.-R., & Oskinova, L. M. 2009, A&A, 494, 1137 Linder, N., Rauw, G., Pollock, A. M. T., & Stevens, I. R. 2006, MNRAS, 370, 1623 Liu, Q. Z., van Paradijs, J., & van den Heuvel, E. P. J. 2000, A&AS, 147, 25 Long, K. S., & White, R. L. 1980, ApJ, 239, L65 Lucy, L. B., & White, R. L. 1980, ApJ, 241, 300 Maeda, Y., Koyama, K., Yokogawa, J., & Skinner, S. 1999, ApJ, 510, 967 Martins, F., Genzel, R., Hillier, D. J., Eisenhauer, F., Paumard, T., Gillessen, S., Ott, T., & Trippe, S. 2007, A&A, 468, 233 Martins, F., Hillier, D. J., Paumard, T., Eisenhauer, F., Ott, T., & Genzel, R. 2008, A&A, 478, 219

34 34 Martins, F., & Plez, B. 2006, A&A, 457, 637 Mauerhan, J. C., Muno, M. P., & Morris, M. 2007, ApJ, 662, 574 (MMM07) Mauerhan, J. C., Muno, M. P., Morris, M. R., Bauer, F. E., Nishiyama, S., & Nagata, T. 2009, ApJ, 703, 30 McLean, I. S., Steidel, C. C., Matthews, K., Epps, H., & Adkins, S. M. 2008, Proc. SPIE, 7014, 99 Mikles, V. J., Eikenberry, S. S., Muno, M. P., Bandyopadhyay, R. M., & Shannon, P. 2006, ApJ, 651, 408 Moffat, A. F. J., et al. 1998, A&A, 331, 949 Moneti, A., Stolovy, S., Blommaert, J. A. D. L., Figer, D. F., & Najarro, F. 2001, A&A, 366, 106 Morris, P. W., Eenens, P. R. J., Hanson, M. M., Conti, P. S., & Blum, R. D. 1996, ApJ, 470, 597 Morris, M., Serabyn, E. 1996, ARA&A, 34, 645 Muno, M. P., Bower, G. C., Burgasser, A. J., Baganoff, F. K., Morris, M. R., & Brandt, W. N. 2006, ApJ, 638, 183 Muno, M. P., et al. 2003, ApJ, 589, 225 Muno, M. P., et al. 2009, ApJS, 181, 110 Nagashima, C. et al. 1999, in Nakamoto T., ed., Proceedings of Star Formation 1999, Nagoya, Japan, p. 397 Nagata, T., Woodward, C. E., Shure, M., & Kobayashi, N. 1995, AJ, 109, 1676 Nagayama, T. et al. 2003, Proc. SPIE4841, p. 459 Najarro, F., Figer, D. F., Hillier, D. J., & Kudritzki, R. P. 2004, ApJ, 611, L105 Nazé, Y., Rauw, G., Pollock, A. M. T., Walborn, N. R., & Howarth, I. D. 2007, MNRAS, 375, 145 Negueruela, I., Smith, D. M., Reig, P., Chaty, S., & Torrejn, J. M. 2006a in ESA SP, 604, The X-ray Universe 2005, ed. A. Wilson, 165

35 35 Negueruela, I., Smith, D. M., Harrison, T. E., & Torrejón, J. M. 2006b, ApJ, 638, 982 Nishiyama, S. et al. 2006, ApJ, 638, 839 Nishiyama, S., Nagata, T., Tamura, M., Kandori, R., Hatano, H., Sato, S., & Sugitani, K. 2008, ApJ, 680, 1174 Oskinova, L. M. 2005, MNRAS, 361, 679 Oskinova, L. M., Hamann, W.-R., Feldmeier, A., Ignace, R., & Chu, Y.-H. 2009, ApJ, 693, L44 Owocki, S. P., Castor, J. I., & Rybicki, G. B. 1988, ApJ, 335, 914 Pfahl, E., Rappaport, S., & Podsiadlowski, P. 2002, ApJ, 571, L37 Pollock, A. M. T., & Corcoran, M. F. 2006, A&A, 445, 1093 Portegies Zwart, S. F., Makino, J., McMillan, S. L. W., & Hut, P. 2001, ApJ, 546, L101 Raassen, A. J. J., van der Hucht, K. A., Mewe, R., Antokhin, I. I., Rauw, G., Vreux, J.-M., Schmutz, W., G&udel, M. 2003, A&A, 402, 653 Ramírez, S. V., Arendt, R. G., Sellgren, K., Stolovy, S. R., Cotera, A., Smith, H. A., & Yusef-Zadeh, F. 2008, ApJS, 175, 147 Ramsay Howatt S.K. et al., 2004, Proc Spie 5492, 1160 Rayner, J. T., Toomey, D. W., Onaka, P. M., Denault, A. J., Stahlberger, W. E., Vacca, W. D., Cushing, M. C., & Wang, S. 2003, PASP, 115, 362 Reid, M. J. 1993, ARA&A, 31, 345 Romano, P., Sidoli, L., Mangano, V., Mereghetti, S., & Cusumano, G. 2007, A&A, 469, L5 Romano, P., Sidoli, L., Cusumano, G., Vercellone, S., Mangano, V., & Krimm, H. A. 2009, ApJ, 696, 2068 Sana, H., Rauw, G., Nazé, Y., Gosset, E., & Vreux, J.-M. 2006, MNRAS, 372, 661 Schnurr, O., Casoli, J., Chené, A.-N., Moffat, A. F. J., & St-Louis, N. 2008, MNRAS, 389, L38 Sguera, V., et al. 2006, ApJ, 646, 452

36 36 Sidoli, L., et al. 2008, ApJ, 687, 1230 Skinner, S. L., Güdel, M., Schmutz, W., & Stevens, I. R. 2001, ApJ, 558, L113 Skinner, S. L., Zhekov, S. A., Güdel, M., & Schmutz, W. 2002a, ApJ, 572, 477 Skinner, S. L., Zhekov, S. A., Güdel, M., & Schmutz, W. 2002b, ApJ, 579, 764 Skinner, S., Güdel, M., Schmutz, W., & Zhekov, S. 2006, Ap&SS, 304, 97 Skinner, S. L., Zhekov, S. A., Güdel, M., & Schmutz, W. 2007, MNRAS, 378, 1491 Smith, L. F., Shara, M. M., & Moffat, A. F. J. 1996, MNRAS, 281, 163 Stolovy, S., et al. 2006, J. Phys. Conf. Ser., 54, 176 Tinney, C. G., et al. 2004, Proc. SPIE, 5492, 998 Tuthill, P., Monnier, J., Tanner, A., Figer, D., Ghez, A., & Danchi, W. 2006, Science, 313, 935 Tuthill, P. G., Monnier, J. D., Lawrance, N., Danchi, W. C., Owocki, S. P., & Gayley, K. G. 2008, ApJ, 675, 698 Usov, V. V. 1992, ApJ, 389, 635 Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389 van der Hucht, K. A. 2001, New Astronomy Review, 45, 135 Walborn, N. R., et al. 2004, ApJ, 617, L61 Wallace, D. J. 2007, in ASP Conf. Ser. 367, Massive Stars in Interactive Binaries, ed. Nicole St.-Louis, & Anthony F.J. Mo?at (San Francisco, CA: ASP), 37 Wang, Q. D., Dong, H., & Lang, C. 2006, MNRAS, 371, 38 Williams, P.M., vanderhucht, K.A.,&Rauw, G.2005,inProc.JENAM2005MassiveStars and High-Energy Emission in OB Associations, ed. G. Rauw, Y. Nazé, R. Blomme, & E. Gosset, 65 Williams, P. M., van der Hucht, K. A., & The, P. S. 1987, A&A, 182, 91 Wright, A. E., & Barlow, M. J. 1975, MNRAS, 170, 41

37 37 Yusef-Zadeh, F., Law, C., Wardle, M., Wang, Q. D., Fruscione, A., Lang, C. C., & Cotera, A. 2002, ApJ, 570, 665 This preprint was prepared with the AAS L A TEX macros v5.2.

38 38 Table 1. Infrared Spectroscopic Observations of Counterparts to X-Ray Sources Associated X-Ray Source Observation Date (UT) Telescope/Instrument λ/δλ Seeing Sky Cond. CXOGC J May 8 11:22 IRTF/SpeX 2000 clear CXOGC J Jul 9 11:06 AAT/IRIS clear CXOGC J Jul 20 08:22 KeckII/NIRC clear CXOGC J Jul 20 08:48 KeckII/NIRC clear CXOGC J May 8 13:53 IRTF/SpeX thin clouds CXOGC J May 9 11:28 IRTF/SpeX thin clouds CXOGC J Jun 15 08:26 UKIRT/UIST thin clouds CXOGC J Jun 15 09:02 UKIRT/UIST thin clouds CXOGC J Jun 16 09:45 UKIRT/UIST clear CXOGC J Jun 28 06:55 UKIRT/UIST thin clouds CXOGC J Jun 28 07:50 UKIRT/UIST thin clouds CXOGC J Jul 20 06:52 UKIRT/UIST clear CXOGC J Jul 22 06:38 UKIRT/UIST thin clouds CXOGC J Jul 22 07:36 UKIRT/UIST thin clouds CXOGC J May 14 12:51 AAT/IRIS clear CXOGC J May 16 16:07 AAT/IRIS clear CXOGC J Jun 16 04:01 SOAR/OSIRIS % cover CXOGC J Jun 16 07:29 SOAR/OSIRIS % cover CXOGC J Jun 17 05:33 SOAR/OSIRIS % cover CXOGC J Jun 14 01:34 SOAR/OSIRIS thin clouds CXOGC J Aug 05 08:45 IRTF/SpeX clear

39 Table 2. Infrared Photometry of Massive Stellar X-ray Sources X-Ray Source R.A. IR Decl. IR J H K s M 3.6 M 4.5 M 5.8 M 8.0 (CXOGC J) (deg, J2000) (mag) (mag) (mag) (mag) (mag) (mag) (mag) a ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± b ± ± ± ± ± ± ± c ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± saturated saturated saturated 0.63 d ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± Note. All known massive stars in the Galactic center with X-ray counterparts are listed, in addition to those discovered in this work. IR source positions and JHK s photometry are taken from Paper I. The λ3.6 λ8.0 µm photometry is from the Spitzer/IRAC point-source catalog

40 of Ramírez et al. (2008). Uncertainties are included where available. a X-ray source considered spurious and not included in main X-ray catalog of Muno et al. (2009). This may be an X-ray source near our limit of sensitivity (see text). b This source did not make the final candidate counterpart list in Paper I; it has an IR positional error greater than 0.1. We suspect this large error is a result of hot dust emission contribution to the K s-band, since this source also has a very bright λ8 µm counterpart. c JHK photometry for this source from Figer et al. (1999). Note that in this case a K-band measurement is presented, rather than the K s band. d λ8 µm data from the MSX Point Source Catalog, assuming the zero-point flux of 58.4 Jy from Cohen et al. (2000). 40

41 Table 3. Spectral Types and Absolute Photometry for WN and O Supergiant X-ray Sources X-Ray Source Spectral K s (J K s) 0 (H K s) 0 J K s H K s A J Ks K s A H Ks K s A Ks M Ks BC Ks Log L bol (CXOGC J) Type (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (L ) O9I B0I WN Ofpe/WN WN7 8h WN5 6b WN8 9h O4 6I WN8 9h O9 B0I OB WN8 9h WN8 9h WN8 9h WN WN5 6b OBI OBI O8 9I O6If O4 6I WN8 9h O4 6I WN8 9h WN7 8h WN7 8h O4 6I Note. Intrinsic colors and BC K values were adopted from Crowther et al. (2006a) for all WN X-ray sources. For the O4 6I stars, intrinsic colors and BC K values were adopted from Martins & Plez (2006) by averaging the values for spectral type O4 6I, and adopting the values for an O9.5I star to represent that of the O9I B0I stars. We did not include WCd or DWCL stars in this analysis, since we do not have knowledge of their intrinsic stellar colors.

42 42 Table 4. Free-free Model Parameters Star Spectral T eff R F ff (3.6, 4.5, 5.8, 8.0 µm) T e N e R ff λ τ=1 Type (kk) (R ) (Jy) (K) (cm 3 ) (R ) (µm) X Ofpe/WN : 1.74 : 1.87 : X WN8 9h : 0.26 : 0.21 : X WN7 8h : 0.31 : 0.26 : X WN5 6b : 0.16 : 0.14 :

43 Table 5. Chandra X-ray Data for Massive Stellar X-ray Sources in the Galactic Center X-Ray Source R. A. X Decl. X σ X Exposure C soft net C hard net F tot (0.5 8 kev) HR0 HR2 E LogL X Spectral Ref. a Ref. b (CXOGC J) (deg, J2000) (arcsec) (sec) (cm 2 s 1 ) (kev) (erg s 1 ) Type X IR (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) < < < < < < < < < < < < < < < < < < < < WN Ofpe/WN WN7 8h WC9d WN5 6b Of 2 6, O4 6I WN8 9h O9I B0I OB WN8 9h 4 6, 8, WN8 9h 4 6, 8, WN8 9h 4 6, 8, WN WN5 6b OBI DWCL OBI DWCL O8 9I O6If WC9d O4 6I DWCL WN8 9h O4 6I WN8 9h WN7 8h WN7 8h O4 6I Note. The table columns contain the following: (1) The source name, which was derived from the coordinates of the sources based on the IAU format, in which least-significant figures are truncated (as opposed to rounded). The names should not be used as the locations of the sources. (2 3) The right ascension and declination of the X-ray source, in degrees (J2000). (4) The 95%-confidence uncertainty in the position (the error circle). (5) The total exposure time in seconds. (6) The net number of counts in the soft kev energy band, with 90%-confidence uncertainties. (7) The net numbers of counts in the hard kev energy band, with 90%-confidence uncertainties. (8) The total kev flux in units of photons cm 2 s 1 (9) The soft color, defined as HR0 = (h s)/(h + s), where h and s is the flux in the kev and kev bands, respectively. The upper and lower 90%-confidence uncertainty is included. An HR value of 9.00 indicates that the sources was not detected in either of the h or s bands. HR uncertainties of 9.00 are considered null, and are present whenever a source is detected in only one band. (10) The hard color, defined as HR2 = (h s )/(h +s ), where h and s is the flux in the kev and kev bands, respectively. The upper and lower 90%-confidence uncertainty is included. (11) The average energy per photon in units of kev. This value can be multiplied by the photon flux in column 8 for an estimate of energy flux.

44 (12) The estimated X-ray luminosity of the source ( kev), assuming N H = cm 2 and kt = 2 kev. Owing to our insensitivity to soft X-ray photons between 0.5 the uncertainties are significant, so, these values should be met with caution (see text). (13) The spectral type of the stellar counterpart. (14) References to the first published association of a massive star with an X-ray source. (15) References to the spectral classification of a source. a The following authors first recognized a massive-star as an X-ray source: (1) This work or Mauerhan et al. 2007; (2) Muno et al. 2006b; (3) Mikles et al. 2006; (4) Yusef-Zad 2002; (5) Hyodo et al b The following authors determined the spectral type of the associated massive star from infrared observations: in addition to the authors listed in a), (6) Figer et al. 2002; (7) al. 1999; (8) Cotera et al or Cotera et al. 1999; (9) Martins et al. 2008; (10) Clark et al. 2009; (11) Liermann et al

45 45 Table 6. L X Compared with Other Authors Source kt (us) N H (us) LogL X (us) kt (them) N H (them) LogL X (them) L X Ref. (CXOGC J) (kev) (cm 2 ) (erg s 1 ) (kev) (cm 2 ) (erg s 1 ) (dex) a Note. This table compares the X-ray luminosities derived via (1) X-ray photometry, adopting constant average values for plasma temperature (kt) and interstellar absorption (N H ), and (2) X-ray spectral fitting with temperature and absorption as free parameters. Only non-variable X-ray sources were used for this comparison. References: (1) Hyodo et al. 2008; (2) Muno et al. 2006; (3) Wang et al a The luminosity of this source is for the kev band.

46 Table 7. Massive Stars Exhibiting X-Ray Variability Record Source Min. ObsID UT Date F min Max. ObsID UT Date F max F max/f min No. (CXOGC J) (10 6 cm 2 s 1 ) (10 6 cm 2 s 1 ) Feb 27 06:26:24 < Jul 06 22:29:57 9.4± Mar 23 12:25:04 7.1± Jul 24 19:58: ± Mar 29 09:44:36 1.1± Jul 16 11:52:55 4.1± Jul 07 19:05:19 1.7± Jun 09 08:50:32 2.7± Note. The columns give: the record number from the catalog of Muno et al. (2009), the source name, the ObsID and UT date in which the minimum flux was observed, the minimum flux, the ObsID and UT date in which the maximum flux was observed, and the maximum flux. The final column contains the ratios of the maximum and minimum observed fluxes, with upper and lower uncertainties. Uncertainties of 1.0 imply that the flux ratio is a lower limit.

47 47 Fig. 1. K-band spectra of late-o/early-b supergiant X-ray sources in the Galactic center, exhibiting absorption lines of H i and He i, and a weak emission feature fromniii near λ2.11 µm. The spectra of X and X resemble that of O9I B0I stars, while the presence of He ii absorption and a relatively shallow Brγ line in the spectrum of X are indicative of a slightly earlier spectral type in the range of O8 9I.

48 48 Fig. 2. K-band spectra of early O4 6I supergiant X-ray sources, dominated by emission from the complex of He i, N iii, C iii, and O iii near λ µm. Weak features of Brγ emission or absorption and He ii absorption at λ2.19 µm are also present. A model spectrum of an O4 6I star is plotted near the bottom of the figure for comparison (Martins et al. 2008).

49 49 Fig. 3. K-band spectra of WN8 9h X-ray sources X and X , and the O6If + star X from Mauerhan, Muno, & Morris (2007). Strong Brγ emission is accompanied by the emission complex of He i, N iii, C iii and O iii near λ µm. A model spectrum of a WN8 9h star is plotted near the bottom of the figure for comparison (Martins et al. 2008).

50 50 Fig. 4. K-band Spectra of WN7 8h X-ray sources, which exhibit stronger He ii emission at λ2.189 µm, relative to the WN8 9h stars presented in Figure 3.

51 51 Fig. 5. K-band spectra of hydrogen-deficient Wolf-Rayet X-ray sources. Broad emission lines of He i and He ii dominate the spectra, indicating high temperatures and fast extended winds. X and X have spectral types of WN7; X and X have spectral types in the range of WN5 6b.

52 52 Fig. 6. K-band spectra of WC counterparts to Galactic center X-ray sources. The weakness of the emission lines is attributable to dilution by strong excess continuum emission from hot dust, and perhaps from a supergiant companion. The spectrum of X is an extreme example of dust dilution, and is completely featureless within the noise level.

53 53 Fig. 7. K-band spectra of the counterpart to X-ray source CXOGC J Strong emission lines of Brγ and He i P-Cygni at λ2.058 µm, are accompanied by weaker lines of Mg ii near λ2.14 µm, and weak He i and He ii P-Cygni profiles near λ2.11 µm and λ2.19 µm, respectively. The spectrum resembles that of Ofpe/WN9 stars.

54 54 Fig. 8. Infrared color-color diagram of massive stellar X-ray sources near the Galactic center, plotted with 1σ error bars. Almost all of the massive stars (large symbols) exhibit significant long-wavelength excess compared to field stars from the SIRIUS near-infrared survey (small dots).

55 55 Fig. 9. Extinction-corrected λ1 10 µm photometry (filled squares) and model SEDs for selected WN stars exhibiting free-free emission in the infrared. For the model SEDs, the dashed lines represent the Rayleigh Jeans portion of the stellar continuum; the dash-dotted lines represent the free-free emission spectrum; and the solid lines are the model sum.

56 56 Fig. 10. Extinction-corrected λ1 10 µm SEDs for new WC9d stars X and X (filled triangles), compared with the dust-enshrouded, extreme DWCL star qf231 of the Quintuplet cluster (squares), the single WC9 star HDM13 and the WC9d star WR 59 (circles). Rayleigh Jeans curves (dotted lines) appropriate for WC9 stars (T eff = 40 kk) are scaled to the J-band data points for comparison with the SED curves. The SEDs of X and X exhibit an infrared excess similar to that of WR 59, consistent with the WC9d class.

57 57 Fig. 11. Photon flux ( kev) vs. hard X-ray color (HR2) for confirmed massive-star counterparts to Chandra X-ray sources. Sources not detected at energies above 3.3 kev were assigned HR2= 1.1 for the figure, but have HR2= 9.00 in Table 5. The massive stars (large dots) are systematically brighter and softer in X-rays than the field population of X-ray sources (small grey points). Sources associated with the Arches (pentacles within circles) and Quintuplet (pentacles) clusters are marked separately for comparison. The data are accompanied by models of optically thin thermal plasma (solid lines) and non-thermal power-law emission (dotted lines) absorbed by a hydrogen column of NH = cm 2. The implied X-ray luminosities of the massive stars are in the range LX erg s 1 ( kev).

58 58 Fig. 12. Scatter plots illustrating the relationship between the X-ray luminosity ( kev) and the stellar bolometric luminosity of the massive stellar X-ray sources near the Galactic center. A linear fit to the data (dashed line) yields a value of L X /L bol = ±0.6, which is consistent with the canonical value of 10 7 (solid line).

59 59 Fig. 13. Spitzer/IRAC images of fields containing massive stellar X-ray sources. Red, green, and blue correspond to 8.0, 5.8, and 3.6 µm images, respectively. Each image spans 144 on a side, and is centered on the labeled sources. Other known massive stars, if present within these fields, are marked by perpendicular pointers. In all images, north is up and east is to the left. Image contrast was scaled by histogram equalization to help reveal faint structure, so the colors of the images are not necessarily illustrative of source SEDs.

60 60 Fig. 14. Spitzer/IRAC images of fields containing additional massive stellar X-ray sources.

61 M ARCHES CLUSTER SGR B1 SGR B2 G F C S R P O N L K I J H E D B A 61 Q QUINTUPLET CLUSTER SGR A WEST 1 DEG. (~140 PC) Fig. 15. Spitzer/IRAC λ8 µm image of the Galactic center region (Stolovy et al. 2006; Arendt et al. 2008), oriented so the Galactic plane is horizontal. The field of each massive stellar X-ray source studied in this work is marked, matching the individual, zoomed color images in Figures 13 and 14.