A&A 483, 171-182 (2008)
Y. Nazé1, - G. Rauw1, - J. Manfroid1,
Institut d'Astrophysique et de Géophysique, Université de Liège, 17 allée du 6 août, Bât. B5C, 4000-Liège, Belgium
Received 15 October 2007 / Accepted 23 December 2007
Context. Westerlund 2 is a young and massive, obscured stellar cluster of our Galaxy. It harbors the most massive star with well-determined parameters, WR20a (82 + 83 ), a dozen very early O-type stars (O3-7), and a wealth of PMS stars. Although of clear significance, this star cluster has not been well-studied.
Aims. The high-energy properties of this cluster, especially those of its early-type stars are examined in detail. The variability of the X-ray sources is investigated.
Methods. A monitoring of the field was performed using three Chandra observations. This dataset probes daily as well as monthly to yearly timescales and provides the deepest X-ray view of the cluster to date.
Results. The two Wolf-Rayet stars WR20a (WN6ha+WN6ha) and WR20b (WN6ha) were analyzed in detail. They are both luminous and display hard spectra, but WR20b does not appear to vary. In contrast, WR20a, a known eclipsing, colliding-wind binary, brightens in the X-ray domain during eclipses, i.e. when the collision is seen face-on. This can be explained by the properties of the wind-wind collision zone, whose high density leads to a large absorbing column ( cm-2).
All twelve O-type stars previously classified spectroscopically, two eclipsing binaries previously identified and nine newly-identified, O-type star candidates are detected in the high-energy domain; ten of those could be analyzed spectroscopically. Four are overluminous, but the others have typical ratios, suggesting that several O-type objects are actually binaries. Variability at the level is detected for a majority of the sources, of unknown origin for single objects.
Faint, soft, diffuse emission pervades the entire field-of-view but no clear structure can be identified, even at the position of a blister proposed to correspond to the TeV source HESS J1023-575.
The X-ray properties of PMS objects, in particular the brightest flaring ones, are also investigated. They provided an additional argument in favor of a large distance (8 kpc) for the cluster.
Key words: Galaxy: open clusters and associations: individual: Westerlund 2 - X-rays: stars - stars: individual: WR20a - stars: individual: WR20b
The young (1-2 Myr) stellar cluster Westerlund 2 lies in the H II region RCW49. It is an intense star forming region with numerous Pre Main Sequence (PMS) objects, as enlighted by recent Spitzer data (Churchwell et al. 2005, and references therein). The cluster also contains many early-type stars: at least two Wolf-Rayet stars, WR20a and WR20b (van der Hucht 2001), and twelve very hot O-stars with spectral types ranging from O3 to O6.5 (Rauw et al. 2007). In the IR and radio domains, these stars appear surrounded by two shells, one enclosing WR20a and the core of Westerlund 2, the other surrounding WR20b.
There are still many uncertainties about the properties and content of this cluster. For example, its distance is still a matter of intense debate. To cite only the most recent papers: Ascenso et al. (2007) and Tsujimoto et al. (2007) proposed a ``low'' value of 2-5 kpc, based on the X-ray and IR luminosities of PMS stars; while Rauw et al. (2007) argue in favor of kpc, based on the spectrophotometric properties of the O-type star population and on the eclipsing characteristics of WR20a. Dame (2007) recently discussed the distance of the Westerlund 2 cluster on the basis of the CO emission of the giant molecular cloud MC8 (likely associated with the cluster) and the H I 21 cm absorptions along the line of sight towards Westerlund 2. This author found that the cluster must be located on the far side of the Carina spiral arm and inferred a kinematic distance of kpc in coarse agreement with the value of Rauw et al. (2007). An indirect argument can be added in favor of a ``high'' value, using the properties of the diffuse emission detected by Townsley et al. (2004). The latter authors found an X-ray luminosity of about 8- erg s-1 (for an adopted distance of 2.3 kpc). Considering that the diffuse emission from massive-star forming regions typically amounts to 2- erg s-1 per early-type star (Townsley et al. 2004,2003), we would expect values of erg s-1 for a cluster containing at least 14 objects earlier than O7 (including the two Wolf-Rayet systems). This is compatible with the observed emission only if the distance is rather large (4-7 kpc rather than 2.3 kpc). In this paper, we will thus adopt 8 kpc as the distance to Westerlund 2.
Studying the X-ray emission of young clusters provides crucial information, especially if the cluster is heavily obscured as in the case of Westerlund 2. PMS objects, mainly T Tauri stars, are more easily detected in the X-ray and IR domains, and their high-energy flares can be analyzed to determine their decay times, which are related to the physical size of the flaring region. In the high-mass range, the interacting colliding-wind binaries and young magnetic objects are more easily identified and their properties constrained in the X-ray domain. Furthermore, the large-scale interaction with the surrounding medium is better studied by the analysis of the diffuse X-ray emission.
The X-ray properties of Westerlund 2 have been studied in increasing detail through the years. The first detection of the cluster was made by the Einstein observatory: Goldwurm et al. (1987) reported the presence of one X-ray emitter at the position of RCW49, interpreted as unresolved sources or truly diffuse emission (from a supernova remnant or a wind-blown bubble). Using ROSAT, Belloni & Mereghetti (1994) detected possible diffuse emission and 24 point sources (7 of them being in our ACIS-I field-of-view), among which HD 90273, MSP18 and WR20a. The high X-ray luminosity of WR20a was suggested as possible evidence of an ongoing colliding-wind phenomenon. More recently, Townsley et al. (2004) and Tsujimoto et al. (2007) reported the results of a first Chandra observation. Thanks to its high spatial resolution, this facility is the first to resolve the individual sources in the core of Westerlund 2. More than 400 X-ray sources were identified by Tsujimoto et al. (2007), mostly PMS stars and a few early-type objects. Faint diffuse emission is also present; it extends to the southeast of the cluster and presents a soft spectrum plus a very faint hard tail (Townsley et al. 2004).
To investigate further the high-energy properties of Westerlund 2, especially those of its massive star population, we have obtained additional X-ray observations, which combine to provide the deepest X-ray view of the cluster. It enables us to study the variability of the X-ray sources on several timescales, short-term variations within one exposure and monthly to yearly changes between the datasets.
The paper is organized as follows. Section 2 presents the datasets and their data reduction. Sections 3 and 4 focus on the O-stars and WR-stars content, respectively. The search for a counterpart to the recently discovered TeV source and the presence of diffuse emission are discussed in Sect. 5, while Sect. 6 concentrates on the flaring sources. Finally, we conclude in Sect. 7.
Westerlund 2 was observed three times using the ACIS-I camera onboard the Chandra X-ray Observatory, first on 2003/08/23 for 38 ks (PI Garmire), then on 2006/09/05 and 2006/09/28 both for 50 ks (PI Rauw); for details see Table 1. The first dataset was taken in very faint mode, whereas the last two were taken in faint mode; in all cases, the frame time was the 3.2 s standard. For all three observations, the ACIS-I chips covered a similar field-of-view but the position angle was slightly different from one exposure to another.
Table 1: Journal of observations. Phases refer to the ephemeris of WR20a (Rauw et al. 2007).
|Figure 1: Three-color images of Westerlund 2 in X-rays ( a), left: full field; b), right: central 8 in the middle and c), below: core of the cluster). The red, green and blue colors correspond to soft (0.3-1 keV), medium (1.-2. keV) and hard (2.-10. keV) photon energies, respectively. The images were corrected by the exposure maps and adaptively smoothed prior to RGB combination.|
|Open with DEXTER|
The datasets were downloaded from the CXC website and we have followed their recommendations for the reduction and analysis of the data. As they had been processed by the latest version of the CIAO pipeline, we did not completely reprocess them.
Background flares were only detected for the first observation, and times with background count rates larger than 0.02 cts s-1 were therefore discarded. The astrometry was improved by comparing the position of the bright X-ray point sources WR20a and MSP18 with their counterparts in the NOMAD catalog. Shifts in pixels of (-0.19, 0.20), (-0.10, 0.04), and (0.14, 0.20) were applied using the CIAO task WCS_UPDATE to the first, second and third observation respectively. The event lists were then reprojected without randomization, first to apply this astrometric correction, then to match the X, Y coordinates from the last observation. A merged event list was created using the CIAO task DMMERGE (with [subspace -expno] to avoid an incorrect merging of the good time intervals). Source detection was performed on this merged event list using the CIAO wavelet algorithm. It was applied separately to 1024 1024 images of the whole field (pixel of 1.4 ), the central 8 (pixel of 0.5 ) and the core of the cluster (pixel of 0.25 ). The derived source positions were used to define source/background regions for use in estimating count rates or in extracting spectra.
Exposure maps were then calculated for each observation and each energy band, assuming a monochromatic source of energy 1 keV (total band, 0.3-10. keV), 0.7 keV (soft band S, 0.3-1. keV), 1.5 keV (medium band M, 1.-2.keV) and 3 keV (hard band H, 2.-10. keV), respectively. Merged exposure maps were obtained using the CIAO task REPROJECT_IMAGE. The resulting fluxed images (in units photons cm-2 s-1 px-2) were used to derive background and exposure-corrected photon fluxes. The merged data were also used to make the three-color images shown in Fig. 1.
Spectrum extraction was performed for each source detected with more than 100 counts in the merged event list. It was performed separately on each observation. As recommended by the CXC, we extracted the spectral response files with the tasks MKACISRMF and MKARF (the latter including a mask file when the considered source was close to a gap and affected by the spacecraft's dithering). When there were less than 100 counts per observation, the individual spectra were merged using the FTOOLS tasks MATHPHA, ADDRMF, ADDARF. Lightcurves were extracted for the brightest sources using the CIAO task DMEXTRACT.
Throughout our analysis, we have used CIAO 3.4, CALDB 22.214.171.124 and Xspec 11.2.
From the photometry of Rauw et al. (2007), we could also select O-star candidates: they have colors typical of early-type stars belonging to the cluster (B-V between 1.3 and 1.6 mag, see Fig. 7 of that paper) and absolute magnitudes smaller than -3.9, a value typical of O9.5V stars (Martins et al. 2005). Some of these candidates are also definite X-ray emitters. The optical properties of this subgroup, composed of 9 objects, are presented in Table 2. Their color excess was derived from the observed color using an intrinsic color of -0.27, typical of O-type stars (Martins & Plez 2006). Then an estimation of their spectral type was made using their absolute V magnitude at an assumed distance of 8 kpc: they are late O-type stars, O7-O9.5. The relevant bolometric correction from Martins & Plez (2006) was then applied to find an estimate of their bolometric luminosities (Table 3) .
We have added the O7 star HD 90273 to the list of O-type objects to analyze. Although this star is much more luminous and less absorbed than the stars of Westerlund 2, we have calculated its luminosity by assuming that it is an O7 supergiant at 8 kpc.
Spectral types, bolometric luminosities and absorbing columns (derived from the color excess using , Bohlin et al. 1978) for the known O-type stars (from Rauw et al. 2007) and the candidates objects (calculated from the photometry as explained above) are shown in Table 3.
|Figure 1: continued.|
|Open with DEXTER|
|Figure 2: The core of the Westerlund 2 cluster seen in visible ( a), left) and X-rays ( b), right). X-ray contours superimposed on the visible image ( c), below; the X-ray data were adaptively smoothed).|
|Open with DEXTER|
|Figure 2: continued.|
|Open with DEXTER|
Table 2: Properties of candidate late-O stars. The last two lines reproduces the data of two eclipsing binaries (from Rauw et al. 2007).
The brightest objects had sufficient counts to provide a detailed X-ray spectrum and they were thus also analyzed spectroscopically. The derived best-fit parameters are shown in Table 4. The second column provides the origin of the spectrum: Obs 1, 2, 3 or Mgd (for the case where a spectrum had insufficient counts to be analyzed individually). The fitted model has the form wabs( )*wabs()*mekal(kT) or wabs( )*[wabs( )*mekal(kT1)+wabs( )*mekal(kT2)], with wabs( ) fixed to the interstellar column (Rauw et al. 2007; Bohlin et al. 1978). In the table, quoted fluxes and luminosities are in the 0.5-10.0 keV energy range and the unabsorbed luminosities are corrected only for the interstellar absorbing column. We note that lightcurves and spectra were extracted in the same regions as defined above. A plot of X-ray luminosity against bolometric luminosity if provided in Fig. 4.
Table 3: Properties of the early-type stars in the optical and X-ray domains. Background- and exposure-corrected source fluxes (in 10-5 photons cm-2 s-1) are estimated for each observation from the fluxed images in a region of radius r (in px) and for the 0.3-10. keV band. The last columns give the average hardness ratios and dereddened X-ray luminosity in 1032 erg s-1 and for the 0.5-10. keV energy band (see text for details). A colon (:) indicates an uncertain value.
Table 4: Best-fitting models and X-ray fluxes at Earth for the brightest early-type stars in Westerlund 2. For each parameter, the lower and upper limits of the 90% confidence interval (derived from the ERROR command under XSPEC) are noted as indices and exponents, respectively. The normalisation factors are defined as , where D, and are respectively the distance to the source (8 kpc), the electron and proton density of the emitting plasma. Abundances were fixed to solar except for the WR stars (for more details, see text).
|Figure 3: Full band background-corrected lightcurves in units of counts per 5 ks bin for the three Chandra observations ( earliest at the top, latest at the bottom). Data from MSP18 are shown in the left panel, those of MSP188 in the middle one and those of WR20b in the right one. Note that the first and last bins, when affected by uncompleteness, were discarded and that the first observation of MSP18 was affected by the spacecraft's dithering (thereby explaining the lower number of recorded counts).|
|Open with DEXTER|
The other bright objects have more typical , with values between -6.5 and -7.3. Among these, MSP157 and 203 display unusually high temperatures. Some of these objects also exhibit 2-3 variations: HD 90273 is 36% less luminous in the last observation, MSP183 and 203 are 35% more luminous in the last observation, the flux of MSP157 almost doubles in the first observation, whereas that of Src2 increases by 35%. If these stars are really single, as their modest luminosity suggests, it is difficult to determine the cause of their variability.
In conclusion, most O-type stars appear to follow the relation derived for the somewhat later O-type stars of NGC 6231. The larger scatter for Westerlund 2 can largely be explained by the fainter X-ray fluxes and the poorer knowledge of their properties, notably multiplicity, compared e.g. to NGC 6231 (Sana et al. 2006). The situation appears definitely more complicated than the faint+soft vs. bright+hard dichotomy suggested by Tsujimoto et al. (2007, see also Fig. 4).
|Figure 4: Left: X-ray luminosity vs. bolometric luminosity. The solid line corresponds to the canonical relation of Sana et al. (2006); the stars represents the Wolf-Rayet objects, the filled circles the brightest O-type stars (see Table 4) and the open circles the faintest ones (see Table 3). Right: hardness ratios, symbols as in the left panel.|
|Open with DEXTER|
|Figure 5: Left: visible lightcurve of WR20a ( top figure), X-ray lightcurves for the full energy band with bins of 2 ks ( middle figure) and 5 ks ( bottom figure). Right: lightcurve in the medium (M, top) and hard (H, middle) bands, together with evolution of the hardness ratio HR2=(H-M)/(H+M) with phase ( bottom), all shown in 5 ks bins. The X-ray data from the second observation are shown in solid black circles and those from the third one in open red triangles; uncomplete bins were discarded. The uncertainty on the phase is about 0.01, taking into account the uncertainty on the period and the time elapsed since our optical observations (700 d, Rauw et al. 2007).|
|Open with DEXTER|
The X-ray spectrum of WR20a appears overall rather hard. This is partially due to the absorption by the large interstellar column density, but reflects also the intrinsic hardness of the spectrum. Indeed, our two-temperature fits of the ACIS spectra of this source indicate temperatures around 1.3-2.0 keV for the hotter component. We emphasize that plasma at temperatures that high is usually not found in the winds of single early-type stars, thus suggesting that a significant part of the X-ray emission of WR 20a arises in the wind interaction region. Assuming that the emission arises from a radiatively-cooling, shock-heated plasma, we find that pre-shock velocities of about 1000 km s-1 are required to account for this hot plasma. This is twice the value of the on-axis pre-shock velocity of the wind-wind interaction that we estimated from a simple velocity law (Rauw et al. 2005). This probably indicates that such a simple approximation of the pre-shock wind velocity law is not appropriate in an interacting-wind, close binary system such as WR20a. It might also mean that the contribution from the off-axis collision, which occurs at higher speed, is not negligible. In what follows we shall thus assume that the pre-shock velocity amounts to 1000 km s-1.
The lightcurve not only reveals the absence of eclipse, it also indicates a brightening of the emission near conjunction phase, i.e. when the collision region is seen face-on. More precisely, on the fluxed images, WR20a appears 20-30% more luminous in the last observation while the photon fluxes in the first two exposures are compatible within 2-3 (Table 3). During each exposure, increasing trends are detected, but individual bins are generally within one or two sigma of a constant value. The last two observations, during which WR20a was unaffected by dithering, follow each other in phase, and the behaviour of WR20a can thus be more easily quantified when the lightcurves are folded in phase (Fig. 5). At the beginning, WR20a presented about 175 counts per 5 ks-bin, while at the end it reached 250 counts per 5 ks-bin. If we compare the lightcurve in the hard and medium bands (see Fig. 5), the increase appears slightly larger in the latter, thereby leading to a slightly smaller hardness ratio in the last observation. The binary is thus significantly X-ray brighter and somewhat softer during the photometric eclipse.
Antokhin et al. (2004) provided a semi-analytical model for a radiative wind interaction zone. These authors propose to calculate the cooling length l0 as where r is the distance between the centre of the star and the stagnation point in , v is the local wind velocity in 103 km s-1, is the angle between the tangent to the interaction front and the local wind direction and is the mass-loss rate in 10-6 yr-1. In the part of the interaction region around the binary axis, that is of most relevance for X-ray emission, the deflection due to the orbital motion is relatively small and hence . For WR 20a, we further adopt r = a/2 = 27 , v = 1000 km s-1 and yr-1 (see Rauw et al. 2005, and Sect. 4.1.1. for v). This then leads to . Due to the steep temperature dependence of the radiative cooling, the width of the postshock cooling layer is of order 2 l0/7 (Antokhin et al. 2004), which then yields about 0.25 . Antokhin et al. (2004) further showed that the surface density of the shock region on the binary axis should be equal to . For WR 20a, this number amounts to about 0.11 g cm-2, implying a mean density in the postshock cooling layer of about g cm-3. If we assume that the wind interaction region within a radius of 12 around the binary axis is delimited by two parallel planar shocks, integrating the density along the direction perpendicular to the binary axis (using equation  of Antokhin et al. 2004) and assuming the same chemical composition as Rauw et al. (2005) yields a (highly ionized) hydrogen column density of order cm-2. This extremely large value over (probably) only a limited fraction of the interaction region shows that the absorption of the X-ray emission by the material of the interaction region must indeed be quite substantial around quadrature phases. Of course, this large value of the column density only applies to a single point (the stagnation point on the binary axis). Since the X-ray emission arises over an extended part of the wind interaction region, the actual (observable) absorption will be lower. In addition, in a highly radiative wind interaction, various kinds of instabilities are expected to affect the shape of the interaction region (Stevens et al. 1992), making it deviate significantly from a simple planar region. The latter effect will reduce the peak value of the column density towards the stagnation point, but will also broaden the episode of high absorption around quadrature phases.
The photon rates and hardness ratio of WR20b are presented in Table 3 while its spectral properties for each observation can be found in Table 4. Overall, the system is 37-38% less luminous than WR20a both in the visible and in X-rays. The spectral fits reveal a large extinction, typical of Wolf-Rayet objects, but also a high temperature. The hardness ratio is also large, even larger than for WR20a (Fig. 4). Such properties are only expected for interacting binaries. However, the X-ray flux appears rather constant on both short and long timescales (Fig. 3 and Table 3), suggesting that the star is single. This contradictory result is important in the context of the current debate about the level of intrinsic X-ray emission from WN stars: while some apparently single WN stars, such as WR40 (Gosset et al. 2005), remain undetected in X-rays, other objects, e.g. WR6 and 110 (Skinner et al. 2002), clearly emit X-rays. WR20b seems to belong to the latter category. An optical monitoring is under way to better constrain the nature of WR20b.
Aharonian et al. (2007) recently reported the discovery of TeV high-energy -ray emission in the vicinity of Westerlund 2 with the HESS Cherenkov telescope array. This emission appears extended with a Gaussian fit yielding , equivalent to a diameter of 28 pc for a distance of 8 kpc. The source center lies 40s (or 5.5) west of WR 20a. Therefore, it is unlikely to be associated with WR 20a alone, especially since the photon-photon absorptions in the wind interaction zone of the colliding-wind binary are likely to be extremely severe, preventing TeV emission from escaping from the inner regions of the wind interaction zone of WR 20a. The size of the source is more typical of a wind-blown bubble. Aharonian et al. (2007) accordingly propose the TeV source to be associated with a blister-like structure detected in the radio range west of the cluster (Whiteoak & Uchida 1997).
To explain the TeV emission, Bednarek (2007) proposed a hybrid model where both leptons and hadrons contribute to the -ray emission from young, open clusters. In this model, the leptons are responsible for X-ray and GeV -rays through synchrotron and inverse Compton radiation, whilst the hadrons produce the TeV -ray emission through pion decay. The electrons and the hadrons would be accelerated at the shocks resulting from the interaction of the stellar winds of the massive stars with the surrounding matter. Whilst the relativistic electrons lose their energy inside the lower density cavity created by the winds of the massive stars, the relativistic protons travel over larger distances and will only lose their energy when interacting with the dense ambient giant molecular cloud. Hence the TeV emission produced by these particles should be more extended than the X-ray and GeV emission.
An alternative explanation has been proposed by Anchordoqui et al. (2007) who suggest that the TeV emission seen towards Westerlund 2 could be the result of Lorentz-boosted MeV -rays emitted in the de-excitation of daughter nuclei produced by the photo-disintegration of PeV cosmic-ray nuclei in the radiation field of the early-type stars inside the cluster core.
In some of the scenarios discussed above, the TeV emission is expected to have an extended counterpart in the X-ray domain. However, in the Chandra observation, the diffuse emission pervades most of the field and appears rather featureless: our X-ray data provide no clear indication of a blister-like structure (Fig. 1). The brightest emission is to be found to the southeast of the cluster, not towards the radio blister, and no bright X-ray point source is detected in the close vicinity of HESS J1023-575. Moreover, the X-ray spectrum of the diffuse emission appears mainly thermal although a faint, hard component might (or might not) be non-thermal (Townsley et al. 2004).
A photometric analysis was carried out on the 0.3-10.keV fluxed images of each observation using the standard procedures for aperture photometry available using IRAF. The positions of the sources have been derived from the merged data in order to increase the signal-to-noise ratio. Apertures with radii between 1.5 to 5 were used. The relative errors were estimated from the photon statistics derived from fluxes in the uncorrected images.
|Figure 6: Light curves of the eight sources that display a strong flare during one of the Chandra observations of Westerlund 2. The sources are identified as in Table 5 and the number of the observation where the flare occured is also indicated. For each panel the x-axis shows the time (in ks) elapsed since the beginning of the observation. The lightcurves were binned with a time bin of 5 ks but note that the first and last bins might be affected by incompleteness.|
|Open with DEXTER|
Table 5: Spectral analysis of the flaring PMS candidates, ordered by RA. Columns 2 and 3 yield the coordinates (J2000) of the sources, whilst the fourth column indicates the number of the observation during which the flare occured. Columns 5-9 list the best fit parameters of the spectral fits with a wabs mekal model. For each parameter, the lower and upper limits of the 90% confidence interval are noted as indices and exponents, respectively. Finally, Col. 10 yields the 1/e decay time of the flare; is the maximum plasma temperature during the flare (derived from the fitted ) and l is the loop half-length.
To identify flaring sources, we have compared the photon rates of each source in the three observations. Many objects appear variable (about 6% at the 99% significance level), but we decided to focus on those that display a pronounced variability between observations ( ) and are bright enough to allow for a detailed analysis of their lightcurve and spectrum. In this way, we have identified eight objects that display a strong flare during one of the observations (see Fig. 6). It should be emphasized that we are not aiming to establish a complete census of variable sources, but to characterize the properties of the brightest flaring candidates. Note that in some cases we find that the observation does not cover the entire flare but only the decay phase.
Correlating the positions of the flaring sources with the 2MASS point-source catalogue (Cutri et al. 2003), we find that sources A, B, D, and G have a near-IR counterpart within less than one arcsecond. However, the only source with good quality 2MASS photometry is Src A (J = 15.47, H = 14.80 and ), the other three have highly uncertain measurements.
Assuming that the flaring plasma is confined in a closed, coronal loop as is the case for solar flares, we can use the simple analytical approach proposed by Serio et al. (1991, for an application of the method see e.g. Favata et al. 2005# to express the loop half-length l as a function of the maximum temperature at the top of the loop
and the thermodynamic decay time
Our X-ray data might however provide hints for constraining the distance of Westerlund 2. A total of 202 X-ray sources have an optical counterpart within 1 of the positions derived using the photometric data of Rauw et al. (2007) (see Fig. 7). The color-magnitude diagram of these sources indicates that most of the X-ray sources with an optical counterpart appear along the main sequence of the cluster or are foreground field stars. The latter objects provide a large fraction of the counterparts, a fact apparently underestimated by Tsujimoto et al. (2007). In addition, the color-magnitude diagram shows no sign of the ZAMS turn-on. None of our flaring sources appear to be associated with an optical source in our photometry. The optical counterparts of these sources must be fainter than V = 21. For very young clusters such as e.g. NGC 6231 (Sana et al. 2007), the optically-brightest PMS stars that display X-ray emission are found to have absolute magnitudes of -1.5, in good agreement with the predictions of evolutionary models (Siess et al. 2000). For a distance of 8.0 kpc ( DM = 14.52) and interstellar extinction AV = 4.7-5.8 (Rauw et al. 2007), the optically-brightest PMS stars in Westerlund 2 would be expected to occur somewhere in the range V = 19.7-21.9. If the cluster is at 2 kpc, the turn-on would be brighter by 3 mag and would be easily spotted in our data, and counterparts to our sources could be detected. Their absence therefore suggests a large distance for the cluster.
|Figure 7: (B-V, V) color-magnitude diagram of the X-ray sources that have an optical counterpart within a correlation radius of less than 1 . The dashed lines correspond to the main-sequence between spectral types O3 and B2 for a distance modulus of 14.52 and reddened with AV = 4.68 (leftmost curve) or 5.84 (rightmost curve) and assuming RV = 3.1.|
|Open with DEXTER|
In this paper, we have analyzed three Chandra observations aimed at monitoring the massive objects of Westerlund 2. In this heavily-obscured region, 24 O-type stars previously identified or candidate were detected and analyzed. Most of these stars display typical values of the ratio, although the scatter is rather large (about 60%) due to the low signal-to-noise ratio of the X-ray spectrum, poorer knowledge of the bolometric luminosities and multiplicities. The most X-ray luminous objects (MSP18, 165, 167 and 188) are probably binaries, but further spectroscopic follow-up observations are required for this to be confirmed. About half of the O-type stars also exhibit variability at the 2 level or more, which is puzzling given their apparent single status.
WR20a, the most massive binary, displays a hard and heavily absorbed X-ray spectrum. It brightens by about 40% in the X-ray domain during the photometric eclipses. The X-ray emission is believed to originate mainly from the wind-wind collision region, which would easily explain the presence of a high temperature in the best-fit models. If the colliding-wind region is sufficiently large and opaque, a modulation will appear when the system rotates and it should then be brighter when seen ``face-on'', as is the case during eclipses. The other Wolf-Rayet star of the field, WR20b, is thought to be single and its flux is accordingly constant. However,the star also exhibits a high temperature and a high X-ray luminosity, and its hardness ratio is even larger than for WR20a. These latter characteristics are more typical of colliding-wind binaries.
The photometry of the other X-ray sources was examined for each observation. Many objects appear variable, and we focused on the eight brightest, flaring ones. Their lightcurve decay time and spectral properties enabled us to infer rather compact flaring regions. The comparison between the optical and X-ray properties of the X-ray sources again suggest a large distance for Westerlund 2.
We acknowledge support from the Fonds National de la Recherche Scientifique (Belgium) and the PRODEX XMM and Integral contracts. We also wish to thank Drs. M. De Becker, H. Sana and I. Stevens for their careful reading of the manuscript and their useful comments.