Articles

A FIRST LOOK AT THE X-RAY POPULATION OF THE YOUNG MASSIVE CLUSTER VVV CL077

, , , , and

Published 2015 March 2 © 2015. The American Astronomical Society. All rights reserved.
, , Citation Arash Bodaghee et al 2015 ApJ 801 49 DOI 10.1088/0004-637X/801/1/49

0004-637X/801/1/49

ABSTRACT

Multi-wavelength analysis of the young massive cluster VVV CL077 is presented for the first time. Our Chandra survey of this region enabled the detection of three X-ray emitting stellar members of the cluster, as well as a possible diffuse X-ray component that extends a few arcseconds from the cluster core with an intrinsic flux of (9  ±  3) × 10−14 erg cm−2 s−1 in the 0.5–10 keV band. Infrared spectra we obtained for two of these X-ray point sources show absorption lines typical of the atmospheres of massive O stars. The X-ray spectrum from the visible extent of VVV CL077  i.e., a 15'' radius around the cluster, can be modeled with an absorbed power law with NH = $(6_{-3}^{+4})\times 10^{22}$ cm−2 and Γ = 2  ±  1. In addition, the X-ray core of VVV CL077 coincides with diffuse emission seen in the infrared band and with a local maximum in the radio continuum map. A possible association with a neighboring H ii region would place VVV CL077 at a distance of around 11 kpc, on the far side of the Norma Arm. At this distance, the cluster is 0.8 pc wide with a mass density of (1–4) × 103M pc−3.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Young massive clusters (YMCs) are loose collections of tens to thousands of stars more massive than the Sun (see Portegies Zwart et al. 2010, for a recent review). While they are a common feature of starburst galaxies, the Milky Way contains relatively few YMCs, only a small fraction of which are known X-ray emitters or have been (or can be) resolved into individual objects over several wavebands (e.g., Townsley et al. 2011a, and references therein). The most massive Galactic YMC is Westerlund 1 which is located 4 kpc away in the direction of the Norma Arm (e.g., Clark et al. 2005). Members of its stellar population have been resolved in several wavebands, including in the X-rays (by Chandra) where evidence was found for X-ray emission from massive, evolved (i.e., post-main-sequence) stars such as OB supergiants and Wolf–Rayet stars, from colliding wind binaries, and from an anomalous X-ray pulsar (Skinner et al. 2006; Muno et al. 2006a; Clark et al. 2008).

X-ray emitting YMCs with resolvable populations represent rare opportunities for observational tests of stellar evolution models including those involving the initial mass function and stellar and cluster dynamics, as well as for understanding the energetic feedback processes between massive stars and their environment. As the more massive stars tend to sink to the center of the cluster, they can pair up into massive binaries. Dynamical interactions between this core population can lead to the tightening of the binary (some of which can evolve into an X-ray binary) or even to its ejection from the cluster (e.g., Portegies Zwart et al. 2007, 2010, and references therein). The distribution of the X-ray emitting point sources compared with the normal stellar content can be used to glean clues to the evolutionary history of the cluster.

We used Chandra to construct a high-resolution X-ray map of a section of the Norma Arm. Among the new objects discovered in the survey were multiple X-ray sources consistent with a YMC named VVV CL077 (Borissova et al. 2011) from the Vista Variables in the Via Lactae Survey (VVV: Minniti et al. 2010). In this paper, we describe the multi-wavelength properties of VVV CL077 by complementing the soft X-ray data with dedicated and archival observations gathered in the radio and infrared (IR) bands. The data and analytical methods are described in Section 2, results are presented in Section 3, and we discuss their implications in Section 4. A summary of our key findings is given in Section 5.

2. OBSERVATIONS AND DATA ANALYSIS

2.1. X-Ray Observations

In 2011 June, the field of VVV CL077 was observed with Chandra during a survey of the Norma Arm (PI: Tomsick; observation IDs 12529 and 12530). An additional observation (ObsID 15625) was performed in early 2013 as part of Chandra Director's Discretionary Time in which VVV CL077 was serendipitously in the field. Table 1 contains the journal of these observations.

Table 1. Chandra Observations

Observation ID Start Date End Date Exposure Time Off-axis Angle of VVV CL077
(s)
12529 2011 Jun 16 06:58:05 2011 Jun 16 12:31:35 19014 3farcm6
12530 2011 Jun 16 12:31:35 2011 Jun 16 18:09:14 19260 8farcm5
15625 2013 Mar 23 08:17:49 2013 Mar 23 11:54:43 9839 6farcm3

Download table as:  ASCIITypeset image

Reprocessing and reduction of the Chandra data relied on CIAO v.4.5 and HEASoft 6.13. Detailed descriptions of the imaging analysis including background, exposure, and astrometric corrections, as well as the source-detection procedure are presented in Fornasini et al. (2014). These cleaned event files and their associated lists of detected sources were used as the inputs for this study. Source coordinates are given in the J2000.0 epoch with uncertainties quoted at 90% confidence.

Spectra were extracted from each event file in the 0.5–10 keV band for source regions (of variable sizes, see Section 3) and from a source-free background region (a rectangle with dimensions: 200'' × 100'' = 400 pixels × 200 pixels) on the same detector chip. Parameters from absorbed power laws fit to the individual spectra were statistically consistent with each other (i.e., there was no significant X-ray variability between 2011 and 2013). Therefore, the X-ray spectra described in the following sections are based on summing these individual spectra with combine_spectra.

2.2. Infrared Observations

The IR data consist of observations we gathered (PI: Rahoui) at the Cerro Tololo Inter-American Observatory. Photometry was performed on 2011 July 19 with the NOAO/NEWFIRM instrument on the 4 m Blanco Telescope. Simultaneous J, H, and Ks-band spectra (R  ∼  1200) were collected on 2012 June 1–4 for five stars near the core of the VVV CL077 using the Ohio State Infrared Imager/Spectrometer (OSIRIS) mounted on the 4.1 m Southern Astrophysical Research Telescope.

A set of standard stars was observed in similar conditions in order to remove telluric absorption lines. The ABBA nodding technique was employed for background subtraction. Data reduction relied on the IRAF echelle package which included bad-pixel correction, dark subtraction, linearity correction, flat-fielding, sky subtraction, and spectral extraction along the dispersion axis. The spectra were then wavelength-calibrated by comparing them to an argon lamp. Additional details regarding the observation and analysis procedures will be described in Rahoui et al. (2014).

Unfortunately, weather and atmospheric conditions were poor with thin clouds and a seeing higher than 1farcs3. Useful spectra were obtained for only two of the five targets, and only in the K band. These data are supplemented by archival imaging observations from the VIRCAM instrument of the VVV (Minniti et al. 2010).

3. RESULTS

3.1. X-Ray Results

In Figure 1, we present images of the field of VVV CL077 in the X-rays (panels (a) and (b): 0.5–10 keV from Chandra) and in the near-IR (Ks-band image from the VVV Survey (Minniti et al. 2010)). Detected X-ray sources are indicated by circles with radii equal to the position uncertainty at 90% confidence (see Table 2). The large circle (radius =15'') denotes the visual extent of the YMC in the IR according to Borissova et al. (2011); it is centered at R.A. =16h34m48s and decl. =−47°32'49''.

Figure 1.

Figure 1. Images of the young massive cluster (YMC) VVV CL077 in X-rays (0.5–10 keV) from Chandra-ACIS-I (panels (a) and (b)) and in the infrared (Ks) band from the VVV Survey (panel (c); Minniti et al. 2010). Each image covers the same region of the sky presented in equatorial (J2000.0) coordinates where North is up and east is left. The X-ray mosaic images combine ObsIDs 12529 and 12530, they are astrometrically corrected, background subtracted, rebinned with 2'' wide pixels, and logarithmically scaled. The X-ray image in panel (b) has been smoothed with a Gaussian (σ = 3 pixels) while the infrared image is scaled with a histogram equalization. Detected X-ray sources are shown with their CXOU designations and circle sizes representing the position uncertainty (at 90% confidence, see Table 2). The large circle (radius =15'') represents the visual extent of the YMC according to Borissova et al. (2011). Contours overlaid on the smoothed X-ray image denote the outline of diffuse infrared emission from the cluster (and bright point sources), while the contours on the Ks-band image signify radio intensities at 843 MHz from the MGPS-2 Survey (Murphy et al. 2007).

Standard image High-resolution image

Table 2. X-Ray Sources Detected in the Vicinity of VVV CL077

Name R.A. Decl. 90% Conf. Radius Det. Sig.
(deg) (deg) (σ)
CXOU J163447.8−473243 248.69950 −47.54550 1farcs1 5.8
CXOU J163448.0−473246 248.70008 −47.54619 1farcs6 6.9
CXOU J163448.5−473229 248.70221 −47.54161 0farcs9 6.2

Download table as:  ASCIITypeset image

This figure illustrates the complex morphology of the multi-wavelength emission from VVV CL077. In the IR, the cluster contains 21 visually identified members (Borissova et al. 2011), plus an extended component from unresolved stars and diffuse emission. In the X-rays, there are three detected sources not including extended emission. The size and shape of the extended emission in X-rays and IR are similar and both have positions that are consistent with a local maximum in the radio continuum map from the Molonglo Galactic Plane Survey 2nd Epoch (MGPS-2; Murphy et al. 2007).

Two candidate X-ray sources are detected near the core of VVV CL077 (Table 2). The first source (in R.A.) is CXOU J163447.8−473243, is detected at a significance of 5.8σ at R.A. $=16^{\mathrm{h}}34^{\mathrm{m}}47\buildrel{\mathrm{s}}\over{.}88$ and decl. =−47°32'43farcs8 with an error radius of 1farcs1. The other, CXOU J163448.0−473246, is detected at 6.9σ at R.A.$=16^{\mathrm{h}}34^{\mathrm{m}}48\buildrel{\mathrm{s}}\over{.}02$ and decl. =−47°32'46farcs3 with an error radius of 1farcs6.

An angular distance of 2farcs9 separates the two X-ray objects, which is somewhat larger than the 90% encircled energy radius (EER; i.e., 2farcs4 at 2.3 keV) of the Chandra point-spread function (PSF). Their X-ray positions are statistically incompatible at 90% confidence, and they are consistent with two bright nodes (stars belonging to VVV CL077) at opposite ends of a clump of unresolved emission in the IR images. Given the positional coincidence with stars seen in the IR, we conclude that these are likely to be two distinct X-ray sources rather than a single extended source.

To confirm whether any of the X-ray emission emanating from the VVV CL077 is extended or simply arises from two point sources, we used ray-tracing simulation data from ChaRT (part of the CIAO package) to generate reprojected PSFs for the observation in which the source had the smallest off-axis angle (ObsID 12529). For this simulation, we adopted an exposure time of 19 ks, and an absorbed power-law spectrum with Γ = 2.2 and NH fixed at 6 × 1022 cm−2 (see below, and in Table 3). The CIAO tool srcextent returns an intrinsic source size of (7farcs0  ±  2farcs5) for this emission, assuming the object is placed on-axis. Using an observation with a larger off-axis angle (ObsID 12530) leads to a larger intrinsic source size (∼15''), as expected.

Table 3. Parameters from Absorbed Power-law and Apec Models Fit to the Unbinned X-Ray Spectra of Selected Regions

NHa Γ b kTc Norm.d C-stat./Goodnesse Sf Hg HRh Obs. Fluxi Unabs. Flux j
VVV CL077: 15'' extraction radius centered on IR position of Borissova et al. (2011)
$6.4_{-2.9}^{+3.7}$ $2.2_{-1.0}^{+1.2}$  ⋅⋅⋅  5.4 73% 48 ± 6 59 ± 10 0.10 ± 0.11 $6.8_{-1.8}^{+2.4}$ $21.6_{-5.9}^{+10.9}$
$5.6_{-2.1}^{+3.2}$  ⋅⋅⋅  $4_{-2}^{+8}$ 9.9 78%  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅  $6.6_{-1.7}^{+2.3}$ $15.0_{-4.5}^{+16.1}$
VVV CL077 diffuse emission: annulus with 6''–12farcs5 radii centered between CXOU J163447.8−473243 and CXOU J163448.0−473246
$5.4_{-3.5}^{+5.0}$ $1.7_{-1.3}^{+1.6}$  ⋅⋅⋅  1.4 73% 16 ± 5 28 ± 7 0.27 ± 0.20 $4.6_{-1.7}^{+2.8}$ $8.9_{-2.5}^{+3.1}$
$5.9_{-2.4}^{+3.7}$  ⋅⋅⋅  $5_{-3}^{+26}$ 5.4 82%  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅  $4.1_{-1.3}^{+1.7}$ $8.9_{-3.2}^{+8.8}$
VVV CL077 inner core: circle with 6'' radius centered between CXOU J163447.8−473243 and CXOU J163448.0−473246
$8.0_{-4.6}^{+8.9}$ $3.1_{-1.6}^{+2.4}$  ⋅⋅⋅  9.2 61% 18 ± 5 22 ± 5 0.10 ± 0.18 $2.8_{-0.9}^{+1.5}$ $27.7_{-6.9}^{+8.1}$
$9.0_{-6.3}^{+8.5}$  ⋅⋅⋅  $1.6_{-0.4}^{+1.0}$ 16.8 70%  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅  $2.4_{-0.7}^{+1.3}$ $19.6_{-14.8}^{+21.1}$
CXOU J163448.5−473229: circle with 6'' radius
$4.5_{-2.5}^{+4.3}$ $3.4_{-1.6}^{+2.9}$  ⋅⋅⋅  4.5 3% 20 ± 5 8 ± 4 −0.43 ± 0.25 $1.3_{-0.5}^{+0.9}$ $13.1_{-6.9}^{+22.3}$
$4.9_{-2.0}^{+2.2}$  ⋅⋅⋅  $1.2_{-0.5}^{+1.2}$ 8.2 2%  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅  $1.1_{-0.4}^{+0.6}$ $12.0_{-8.6}^{+36.6}$

Notes. aColumn density in units of 1022 cm−2. bPhoton index of the power-law model. cPlasma temperature (keV) for the apec model. dModel normalization (× 10−5). eGoodness of fit defined as the percentage of simulated spectra (of 10,000 trials) which return a statistic smaller than the observed value. fNet source counts in the soft (S) band: 0.5–3 keV. gNet source counts in the hard (H) band: 3–10 keV. hHardness ratio defined as (HS)/(H + S). iObserved flux (i.e., not corrected for absorption) in units of 10−14 erg cm−2 s−1 in the 0.5–10 keV band. jAbsorption-corrected flux in units of 10−14 erg cm−2 s−1 in the 0.5–10 keV band.

Download table as:  ASCIITypeset image

After subtracting off 95% of the photons from the two X-ray point sources in the core via aperture restriction (see below), we find that there is still some residual emission with as many net counts outside the core (44) as from the two point sources combined (40). Along with the large intrinsic size of the X-ray emitting region, this indicates the presence of an extended and diffuse X-ray component, in addition to the emission from the two point sources.

Just outside the YMC is CXOU J163448.5−473229 which coincides spatially with a star that is likely associated with the YMC given the similar IR colors to those of other cluster members. The best-fitting X-ray position (6.2σ detection significance) is R.A. $=16^{\mathrm{h}}34^{\mathrm{m}}48\buildrel{\mathrm{s}}\over{.}53$ and decl. =−47°32'29farcs8 with an error radius of 0farcs9. Its intrinsic size, estimated via ChaRT, is consistent with 0'' in ObsID 12529.

Note that in ObsID 12530, a second candidate X-ray source is detected outside the YMC, CXOU J163448.6−473234, located at R.A. $=16^{\mathrm{h}}34^{\mathrm{m}}48\buildrel{\mathrm{s}}\over{.}62$ and decl. =−47°32'34farcs1 with an error radius of 2farcs5 (significance =5.1σ). The large estimate of the source extent from wavdetect (=10''), the fact that it was only detected in this ObsID when the YMC was 8farcm5 off axis, and the lack of a clear IR counterpart suggest that this is the same source as CXOU J163448.5−473229 (4farcs4 away), and so it will not be discussed further.

Spectral data were fit with power-law, apec, and bremsstrahlung models. However, the observed photon counts within the source extraction regions are low (i.e., ≲ 100 counts), and so we caution that any derived spectral parameters are merely estimates, and may not be indicative of the true values. Given the low source counts and very low background counts, spectral fitting relied on Cash (1979) statistics and Pearson (1900) χ2 test statistics. All fits assumed Wilms et al. (2000) abundances and photoionization cross sections from Balucinska-Church & McCammon (1992). The hardness ratio is defined as (HS)/(H + S) where H is the count rate in 3–10 keV and S is the count rate in 0.5–3 keV.

The X-ray spectrum extracted from a 15'' radius (30 pixels) around the IR position of VVV CL077 is presented in Figure 2. This spectrum includes the two X-ray point sources in the core of the YMC, plus diffuse X-ray emission. There are 107 net counts in this region for a hardness ratio of 0.10 ± 0.11. A power-law fit to the data yields NH $=(6.4_{-2.9}^{+3.7})\times 10^{22}$ cm−2 and $\Gamma = 2.2_{-1.0}^{+1.2}$, while an apec model yields NH $=(5.6_{-2.1}^{+3.2})\times 10^{22}$ cm−2 and $kT = 4_{-2}^{+8}$ keV.

Figure 2.

Figure 2. X-ray spectrum extracted from a 15'' radius (30 pixels) around the infrared position of VVV CL077. The lower panel shows residuals from fitting an absorbed power law (solid curve) to the spectral data. For visual clarity, each spectral bin contains a minimum of 15 counts. Spectral parameters from the unbinned analysis are listed in Table 3.

Standard image High-resolution image

To compare the properties of the X-ray emission in the core of the YMC with that of the diffuse emission, we extracted a spectrum from a circle of 6'' radius (12 pixels) that encompasses the 95% EER of both point sources in the core (CXOU J163447.8−473243 and CXOU J163448.0−473246) giving 40 net source counts, as well as a spectrum from an annulus (centered at the same position) with inner and outer radii of 6'' and 12farcs5, respectively, that represents the diffuse emission (44 net counts). The results of models fit to these spectra are tabulated in Table 3. There are hints that the diffuse component has a slightly higher plasma temperature (with apec models) and a harder photon index (with power-law models) than the stellar sources in the core region, but the low photon counts results in large uncertainties, which leads to some overlap in the range of possible temperatures and photon indices at the 90% confidence level. We attempted to extract spectra for each X-ray point source in the core by using non-overlapping regions but the sources ended up with too few counts (nine net counts each) to permit reliable fitting.

A source spectrum was made for CXOU J163448.5−473229 by using an extraction radius of 6'' (12 pixels) centered on its X-ray position. While the spectral parameters of CXOU J163448.5−473229 are compatible with those of the other source regions (within the statistical uncertainties), it is clearly softer than the others given its hardness ratio of −0.43 ± 0.25. This is illustrated in Figure 3 which compares the 68%, 90%, and 99% confidence contours of spectral parameters from the power-law model fit to the data of CXOU J163448.5−473229 (dashed lines), and for the YMC (solid lines). There is some overlap between the parameter spaces at the 99% confidence level, but the X-ray emission from the cluster tends to be harder and more absorbed than that of a point source situated in the cluster's outskirts.

Figure 3.

Figure 3. Parameter space of NH and Γ resulting from absorbed power-law fits to the X-ray spectra (0.5–10 keV) of selected extraction regions. Two extraction regions are considered: a 15'' radius circle centered on the IR position of VVV CL077 (solid lines), and a 6'' radius circle centered on the X-ray position of CXOU J163448.5−473229 (dashed lines). The best-fitting parameters are represented by crosses and the contours denote confidence levels of 68%, 90%, and 99%.

Standard image High-resolution image

3.2. Infrared Results

Figure 4 shows the Ks-band image of the region gathered with NOAO/NEWFIRM. Infrared spectra for Source #2 and Source #3 are presented in Figure 5. In the Ks-band spectra of Source #2, we detect a Brγ line in absorption at 21660 Å. The equivalent width is 11.5 Åwith a typical error of ∼1–1.5 Å. There is a marginal detection of the He i line in absorption at 20600 Å. Source #3 also presents a Brγ line in absorption with an equivalent width of 6.5 Å. With standard deviations being less than 0.06, this leads to signal-to-noise ratio >10 in both sources. Comparing these spectra to those of a sample of high-mass stars provided by Hanson et al. (1996) indicates that Sources #2 and #3 are probably late-O giants.

Figure 4.

Figure 4. Image of VVV CL077 in the infrared Ks band from NOAO/NEWFIRM. Circles denote the positional uncertainty of the X-ray point sources described in this work (Table 2), and the visual extent of the YMC according to Borissova et al. (2011). Infrared objects targeted for follow-up spectroscopy are numbered 1–5.

Standard image High-resolution image
Figure 5.

Figure 5. Infrared Ks-band spectra of two stars (labeled "2" and "3" in Figure 4) in the YMC VVV CL077 from observations gathered with NOAO/OSIRIS.

Standard image High-resolution image

In addition to these spectra, we searched the ViZier archives for multi-wavelength counterparts consistent with the X-ray error circles of the sources detected by Chandra. Table 4 lists the reported magnitudes of these counterpart candidates. The derived JH and HKs colors are compatible with the colors observed in O stars (e.g., Ramírez Alegría et al. 2012). Comparing CXOU J163448.0−473246 and CXOU J163448.5−473229, we find consistent values (within the uncertainties) for their JH (1.8 ± 0.2 versus $2.0_{-0.2}^{+0.1}$) and JKs ($3.3_{-0.2}^{+0.1}$ versus $3.0_{-0.2}^{+0.1}$) colors, but the HKs color is ∼50% higher for CXOU J163448.0−473246 (1.5 ± 0.2 versus $1.0_{-0.1}^{+0.2}$) which suggests it has a slight near-to-mid-IR excess, although the source is blended in the Two Micron All Sky Survey (2MASS) Ks-band image. We point out the presence of a mid-IR source located 2farcs9 away from CXOU J163448.0−473246 (i.e., 1farcs4 from the IR position of the cluster): WISE J163448.13−473248.9 has reported magnitudes of 5.563 ± 0.110 (at 3.35 μm), 4.756 ± 0.055 (4.6 μm), −0.896 ± 0.121 (11.6 μm), and −3.190 ± 0.020 (22.1 μm).

Table 4. Infrared Counterpart Candidates to X-Ray Sources Near VVV CL077

Name IR Counterpart Offset J H Ks 3.6 μm 4.5 μm 5.8 μm
('')
CXOU J163447.8−473243  ⋅⋅⋅  a  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅ 
CXOU J163448.0−473246 2MASS J16344796−4732452 1.2 14.485 ± 0.087 12.700 ± 0.137 11.230 ± 0.081  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅ 
  DENIS J163447.9−473245 0.8 14.10 ± 0.11  ⋅⋅⋅  11.12 ± 0.07  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅ 
  VVV 515726976188 0.2 14.731 ± 0.005 13.466 ± 0.005  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅ 
CXOU J163448.5−473229 2MASS J16344850−4732297 0.3 13.729 ± 0.075 11.772 ± 0.075 10.730 ± 0.043  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅ 
  VVV 515726978416 0.6 13.647 ± 0.002 11.898 ± 0.001  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅ 
  GLIMPSE G336.8914+00.0524 0.4  ⋅⋅⋅   ⋅⋅⋅   ⋅⋅⋅  10.052 ± 0.048 9.844 ± 0.053 9.818 ± 0.103

Notes. Infrared magnitudes of counterpart candidates in: the J, H, and Ks bands from 2MASS (Skrutskie et al. 2006), DENIS (DENIS Consortium 2005), and VVV (Minniti et al. 2010); and at 3.6 μm, 4.5 μm, and 5.8 μm from Spitzer–GLIMPSE (Spitzer Science 2009). aThe nearest IR source with cataloged magnitudes, 2MASS J16344796−4732452, is 1farcs7 away and is likely associated with CXOU J163448.0−473246.

Download table as:  ASCIITypeset image

4. DISCUSSION

At least 21 stars belonging to VVV CL077 were visually identified in the IR (Borissova et al. 2011), and this YMC is now detected for the first time in the X-rays with Chandra. A complete census of the X-ray emitting sources in VVV CL077 cannot be achieved with the current data. Nevertheless, recent X-ray observations of this region reveal a number of intriguing results.

Imaging data from Chandra suggest that there are three X-ray point sources associated with IR stellar counterparts that belong to the cluster. Two of these sources are located near the cluster core, and the third is located in the cluster's outskirts. In addition to these point sources, there is extended and diffuse X-ray emission: the estimated intrinsic (on-axis) size of the X-ray emitting core of the YMC is larger than expected for the combined 95% energy-containment circles of the two closely spaced point sources in the core; and this component, which contains as many net counts as the core X-ray point sources combined, cannot be directly associated with any resolved sources.

The presence of diffuse emission, with X-ray spectral properties statistically consistent with that of the resolved stellar X-ray sources, would suggest a population of unresolved X-ray sources in the YMC whose emission lies below the sensitivity limit of the survey (∼10−15 erg cm−2 s−1; Fornasini et al. 2014).

Diffuse X-ray emission has been observed in other YMCs including in Wd 1 and in active star-forming regions like the Carina Nebula (Muno et al. 2006b; Townsley et al. 2011b, 2011c). In addition to the combined emission from unresolved stars, the diffuse component could originate from putative collisions between the winds of massive stars, although the temperature and absorption in VVV CL077 are higher than expected for such processes (e.g., Stevens & Hartwell 2003).

Two of the stars VVV CL077 display IR spectral signatures typical of late-type O stars (Hanson et al. 1996). There is good overall agreement between the outlines of the emission as seen in the X-rays and IR, and this multi-wavelength consistency continues into the radio band where the X-ray and IR emission coincide spatially with a local maximum in the 843 MHz continuum map from MGPS-2 (Murphy et al. 2007). The coincidence of the multi-wavelength emission from VVV CL077 could indicate a common physical origin; however, the observed X-ray flux of the diffuse component is larger than expected from an extrapolation of the IR flux density (see below).

The nearest cataloged ultra-compact H ii region is IRAS 16311−4726 (Helou & Walker 1988; Bronfman et al. 1996) which is 6'' from the IR position of VVV CL077. It has a velocity of 121 km s−1 with respect to the local standard of rest, which translates to a source distance of approximately 11 kpc (e.g., Brand & Blitz 1993). There is also a cataloged H ii region 7' away from VVV CL077 with a listed distance of 10.9 ± 0.2 kpc (Russeil 2003). If these regions are related, then we can assign a distance of around 11 kpc to VVV CL077. This would make it one of the most distant Galactic YMCs known, comparable to Mercer 81 (Davies et al. 2012), at a location consistent with the far Norma Arm.

Taking the visual extent of the bulk of the stellar population in the IR as the angular scale (15''; Borissova et al. 2011) gives a cluster diameter of 0.8 pc for an assumed distance of 11 kpc. This value is consistent with the intrinsic sizes of other Galactic YMCs (Figer 2008, and references therein).

Corrected for absorption, the diffuse emission from VVV CL077 has an intrinsic flux of (9 ± 3) × 10−14 erg cm−2 s−1 (0.5–10 keV). This is more than an order of magnitude above the flux expected from the Galactic Ridge diffuse emission (Ebisawa et al. 2001, 2005). If we assume, nevertheless, that the Ridge emission accounts for 20% of the flux from VVV CL077, as it does for Wd 1 (Muno et al. 2006b), then the residual flux of ∼7 × 10−14 erg cm−2 s−1 can be attributed to the diffuse cluster emission. For a distance of 11 kpc (assumed), this translates to an X-ray luminosity of 1034 erg s−1. Flux densities reported for IRAS 16311−4726 are 20.1 ± 0.6 Jy at 12 μm and 153 ± 5 Jy at 25 μm (Helou & Walker 1988). Extrapolating a power law with a photon index of Γ = 2.8 ± 0.1 leads to an X-ray flux density of 1.3 × 10−10 Jy at 1 keV (3.0 × 10−15 erg cm−2 s−1 over 0.5–10 keV). This is an order of magnitude smaller than the lower 90% confidence boundary of the intrinsic flux density (3 × 10−9 Jy at 1 keV). In order to match the X-ray flux, the photon index should be no steeper than 2.5, which is consistent with what we found in Section 3.1.

A simple estimation of the cluster mass density yields ρ ∼(1–4) × 103M pc−3. This assumes that all 21 visually identified members of the cluster are similar to the sources we identified as O type stars with masses in the range of 10–50 M. While the mass density is consistent with those of other YMCs, the cluster mass implied by this first-order approximation (∼103M) is an order of magnitude below the inferred masses of other YMCs (Figer 2008). This indicates a large population of unresolved stars hiding in and around the cluster which is not surprising given the extinction in this direction.

From the observed IR excess E(HKs) = 1.5 ± 0.2, we estimate the extinction in the optical band as AV = 15 ± 2 assuming a near-IR extinction law with α = −2.1 (Stead & Hoare 2009). Applying the relation of Güver & Özel (2009) yields a column density of NH $=(3.3_{-1.0}^{+1.1})\times 10^{22}$ cm−2. When the contribution from molecular hydrogen is included, the total column density of the interstellar medium toward the Norma region ranges from 4.9 × 1022 cm−2 to 9.4 × 1022 cm−2 (Rahoui et al. 2014). This is consistent with our X-ray-derived value of $(6.4_{-2.9}^{+3.7})\times 10^{22}$ cm−2.

Plasma temperatures of around 1–2 keV are found for the X-ray emitting point sources in the cluster which is higher than average for massive stars (e.g., Rauw et al. 2014). The X-ray source located in the outskirts of the cluster (CXOU J163448.5−473229) has a temperature that is statistically compatible (at 90% confidence) with 0.7 keV which is a more typical value. As we move toward the center of the cluster, we find that the X-ray derived temperatures of these stars are higher, possibly due to increasing column density (which preferentially absorbs low-temperature photons). With its proximity to IRAS 16311−4726, this could indicate that the stars in the core of VVV CL077 are still embedded within their dusty natal cocoons; similar to the hot, young stars of W49A (Conti & Blum 2002), a cluster that resides within an ultra-compact H ii region. However, we cannot exclude that the unusual spectral values are the result of poor photon statistics.

Based on the IR spectra, we narrowed the range of spectral types to O5–O9 giants for two of the stars in the cluster. We derived the distance modulus mM = 5log (d) + A − 5 using the observed J-band magnitude of CXOU J163448.0−473246 in Table 4, the absolute J-band magnitude of −4.6 for an O giant (Martins & Plez 2006), and AJ = 0.28AV. This gives a distance to the object of 10 ± 3 kpc, consistent with the distance assumed based on its possible link with IRAS 16311−4726.

In massive stars, the ratio between the X-ray and bolometric luminosities is ∼10−7 (e.g., Berghoefer et al. 1997; Rauw et al. 2014). Such stars have bolometric luminosities of log (L/L) ∼ 5.2–5.7 (Martins et al. 2005). To test whether the objects we classify as massive stars obey this rule, we began with the observed X-ray flux of the "inner core" region and corrected this for the interstellar absorption alone (NH =3 × 1022 cm−2), then halved this value to obtain the absorption-corrected flux of one of the stars of the pair. For an assumed distance of 10 ± 3 kpc, this yields an X-ray luminosity of 3.5 × 1032 erg s−1. A ratio LX/Lbol of 10−7 is obtained for either the high end of the bolometric luminosity range (log (L/L) ∼5.7), or for distances at the lower end of the estimate (∼7 kpc). Fixing the temperature to a more reasonable value of 0.6 keV raises the column density to unlikely values (NH $=(22_{-5}^{+7})\times 10^{22}$ cm−2), which, once corrected for interstellar absorption, results in LX  ∼  2.5 × 1032 erg s−1, and so it remains consistent (within the uncertainties) with the expected LX/Lbol value for massive stars.

Finally, we consider the question of whether there are supernova remnants and accreting compact objects in VVV CL077. There are no obvious shell-like features in the IR and X-ray data that would indicate recent supernova activity, although the low level of soft X-ray emission and lack of hard X-ray or gamma-ray detections appear to rule out such scenarios. Luminous (LX ≳ 1038 erg s−1) X-ray point sources are missing from VVV CL077 which excludes high-mass X-ray binaries (HMXBs) where the accretion proceeds through Roche-lobe overflow. Wind-fed HMXBs could be present but they would have to be weakly accreting or quiescent (LX  ∼  1033 erg s−1), such as is the case for Supergiant Fast X-ray Transients (SFXTs) outside of outbursts. However, the photon indices measured in the core and for the external X-ray sources are too soft to support the weakly accreting SFXT scenario (e.g., Romano et al. 2011). An additional point to consider is that the supernova that leads to the formation of the compact object in such systems can impart a velocity to the nascent X-ray binary that is several times larger than the cluster escape velocity (van den Heuvel et al. 2000; Clark et al. 2005). The apparent lack of supernova remnants or accreting compact objects would tend to favor a relatively young age for the YMC, i.e., a few Myr at most.

5. SUMMARY AND CONCLUSIONS

We presented results from recent Chandra observations of VVV CL077, a YMC located toward the Norma Arm. Three X-ray point sources were detected that correspond to visually identified stellar members of the cluster. Infrared spectra we obtained for two members of the cluster show absorption lines typical of the atmospheres of massive stars. Plus, the IR colors and LX/Lbol relation also indicate massive stars. However, X-ray spectral fits of these stars yield temperatures that are higher than average for massive stars.

In addition to these point sources, there are hints that an extended region of diffuse X-ray emission permeates the cluster core. The X-ray column density and optical extinction, as well as an estimate of the cluster mass, suggests a population of unresolved stars that contribute to the diffuse emission. The X-ray spectrum of the 15'' radius visible extent of the VVV CL077 is best modeled with an absorbed [NH =(3–10) × 1022 cm−2] power law (Γ = 2 ± 1).

The X-ray core of the cluster coincides with diffuse emission seen in the IR band and with a local maximum in the radio continuum map. A tentative link with a neighboring H ii region would place VVV CL077 at a distance of around 11 kpc; on the far side of the Norma Arm. At this distance, the cluster is 0.8 pc wide with a mass density of (1–4) × 103M pc−3, i.e., within the range of values seen for other YMCs.

A complete picture of the X-ray population and diffuse emission from VVV CL077 is not possible at this time, but continued X-ray observations of YMCs are important for the study of massive stars, the formation and evolution of compact objects, and the recent star formation history of the Galaxy.

The authors thank the anonymous referee whose careful, patient, and constructive review led to significant improvements in the manuscript. This work was supported in part through Chandra Award Number GO1-12068A issued by the Chandra X-Ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. F.E.B. acknowledges support from Basal-CATA (PFB-06/2007) and CONICYT-Chile (under grant FONDECYT 1141218 and Anillo ACT1101). The scientific results reported in this article are based on observations made by the Chandra X-Ray Observatory. This research has made use of: software provided by the Chandra X-ray Center (CXC) in the application packages CIAO, ChIPS, and Sherpa; data obtained from the High Energy Astrophysics Science Archive Research Center (HEASARC) provided by NASA's Goddard Space Flight Center; NASA's Astrophysics Data System Bibliographic Services; the SIMBAD database operated at CDS, Strasbourg, France; data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by NASA and the NSF; data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by NASA and the NSF.

Facilities: CXO - Chandra X-ray Observatory satellite, Blanco (NEWFIRM) - Cerro Tololo Inter-American Observatory's 4 meter Blanco Telescope, SOAR (OSIRIS) - The Southern Astrophysical Research Telescope

Please wait… references are loading.
10.1088/0004-637X/801/1/49