A&A 490, 1055-1070 (2008)
DOI: 10.1051/0004-6361:200809999

An X-ray survey of low-mass stars in Trumpler 16 with Chandra[*]

J. F. Albacete-Colombo1,2 - F. Damiani2 - G. Micela2 - S. Sciortino2 - F. R. Harnden, Jr.3


1 - Centro Universitario Regional Zona Atlántica - Univ. COMAHUE, Monseñor Esandi y Ayacucho (CP 8500), Viedma, Argentina
2 - INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy
3 - AF Smithsonian Astrophysical Observatory, 60 Garden St., Cambridge, MA 02138, USA

Received 17 April 2008 / Accepted 23 July 2008

Abstract
Aims. We identify and characterize low-mass stars in the $\sim $3 Myr old Trumpler 16 region by means of a deep Chandra X-ray observation, and study their optical and near-IR properties. We compare the X-ray activity of Trumpler 16 stars with the known characteristics of Orion and Cygnus OB2 stars.
Methods. We analyzed a 88.4 ks Chandra ACIS-I observation pointed at the center of Trumpler 16. Because of diffuse X-ray emission, source detection was performed using the PWDetect code for two different energy ranges: 0.5-8.0 keV and 0.9-8.0 keV. Results were merged into a single final list. We positionally correlated X-ray sources with optical and 2MASS catalogs. Source events were extracted with the IDL-based routine ACIS-Extract. X-ray variability was characterized using the Kolmogorov-Smirnov test and spectra were fitted by using XSPEC. The X-ray spectra of early-type, massive stars were analyzed individually.
Results. Our list of X-ray sources consists of 1035 entries, 660 of which have near-IR counterparts and are probably associated with Trumpler 16 members. From near-IR, color-color, and color-magnitude diagrams we compute individual masses of stars and their $A_{\rm v}$ values. The cluster median extinction is $A_{\rm v}$ = 3.6 mag, while OB-type stars appear less absorbed, having $A_{\rm v}$ = 2.0 mag. About 15% of the near-IR counterparts show disk-induced excesses. X-ray variability is found in 77 sources, and typical X-ray spectral parameters are $N_{\rm H}\sim$ 5.37 $\times $ 1021 cm-2 and $kT\sim$ 1.95 keV. The OB stars appear, softer with a median $kT\sim$ 0.65 keV. The median X-ray luminosity is 6.3 $\times $ 1030  \ensuremath{{\rm erg~s^{-1}}}, while variable sources show a larger median $L_{\rm x}$ value of 13 $\times $ 1030  \ensuremath{{\rm erg~s^{-1}}}. OB-stars have an even higher median $L_{\rm x}$ of 80 $\times $ 1030  \ensuremath{{\rm erg~s^{-1}}}, about 10 times that of the low-mass stars.
Conclusions. The Trumpler 16 region has a very rich population of low-mass X-ray emitting stars. A large fraction of its circumstellar disks have survived the intense radiation field of its massive stars. Stars with masses 1.5-2.5 $M_\odot $ display X-ray activity similar to the Cyg OB2 stars, but much less intense than observed for Orion nebula cluster members.

Key words: stars: formation - stars: early-type - stars: pre-main sequence - X-rays: stars

   
1 Introduction

The Carina nebula region (NGC 3372) is one of the most massive star formation regions of the Galaxy. It is associated with a giant H II region spanning about 4 deg2 of the sky and is bisected by a prominent V-shaped dark gas and dusty lane. This prominent young structure is not as compact as some of the other young galactic clusters, but seemingly to be related to a spiral feature. In this direction, we are looking almost tangentially toward the now recognized Carina-Sagittarius spiral arm, at the edge of a giant molecular cloud extending over about 130 pc, which has a content in excess of 5 $\times $ 105 solar masses (Grabelsky et al. 1988). The concentration of massive stars (i.e. $M\geq$ 20 $M_\odot $) interacts with the parent giant molecular cloud of the region, leading to triggered star formation events on intermediate to lower masses (e.g. Smith et al. 2004).

This region harbors several open clusters and/or star concentrations (Trumpler 14-16; Collinder 228 and 232; Bochum 10 and 11) containing more than 60 known O-type stars (Feinstein 1995). Large cavities within the giant molecular cloud are supposed to be carved out by the Tr 14 and 16 open clusters, which contain most of massive stars of the region. In particular, Tr 16 includes three rare main-sequence O3 stars, the Wolf-Rayet (WR) star HD 93162 and the famous luminous blue variable (LBV) $\eta $ Carinae. There is a historical controversy about the distance and age of Tr 14 and Tr 16 (Walborn 1995). For instance, from extensive spectroscopy and photometry Massey & Johnson (1993) find 3.2 kpc for both clusters. However, photometric studies are strongly affected by differential extinction in the region and peculiar reddening, so the derived distance are different. An example is the Carraro et al. (2004) work, who compute distances of 4.0 kpc and 2.5 kpc, for different R= $A_{\rm v}$/E(B-V) values (3.48 and 4.16 for Tr 16 and Tr 14 regions). A more reliable distance (2250  $\pm$ 180 pc), was derived from proper motion and Doppler velocities of the expanding $\eta $ Carinae Homunculus using HST-STIS[*] observations (Davidson & Humphreys 1997). Recent work (Tapia et al. 2003) derives a common distance DM = 12.14 (2.7 kpc) and an age between $\sim $1 Myr and 3 Myr, for Tr14 and Tr16, respectively. For this study, we adopt a distance for Trumpler 16 of 2250 pc and an age of 3 Myr. This young age agrees with the Smith et al. (2000) results, which report the existence of several embedded IR sources where star formation might be active. Also, DeGioia-Eastwood et al. (2001) confirm clear evidence of pre-main sequence (PMS) stars in the region, while Brooks et al. (2001) have identified two compact H II regions possibly linked to very young O-type stars. Finally, Hägele et al. (2004) report a compact cluster of infrared PMS-stars in Tr 16.

Of the existing methods of identifying young stellar populations, the use of X-ray emission is perhaps the least biased (Feigelson et al. 2002). While in main-sequence (MS) stars, from late A to M dwarfs, X-rays are believed to originate from the hot coronal gas that is heated by stellar dynamo magnetic fields (Maggio et al. 1987), for late type Pre-MS stars (T Tauri stars (TTSs)) X-ray emission is attributed to solar-like coronal activity but elevated by a factor of 103-104 (Feigelson & Montmerle 1999). Several authors suggested the possibility of detecting early PMS objects through their hard X-ray emission escaping the highly obscured regions (see Kamata et al. 1997; Hofner & Churchwell 1997; Hofner et al. 2002; Walter 1992). Recently, X-ray surveys have been successful in identifying the young and pre-MS population in star-forming regions, including: i- deeply embedded Class I young stellar objects (YSOs); ii- low-mass T Tauri PMS stars; iii- intermediate-mass Herbig Ae/Be PMS stars; and iv- zero-age MS stars. Moreover, X-ray emission from low-mass PMS stars usually exhibits a strong variability that helps to confirm membership.

In the last decade, X-ray observations of young stars on star-forming regions were intensified thanks to the high spatial resolution and the improved broad-band ([0.2-12.0] and [0.5-10.0] keV) effective area of the XMM-Newton and Chandra satellites. A first X-ray survey in the Carina region by Albacete Colombo et al. (2003) was performed on the basis of two early XMM-Newton observations (rev #115 and #116) centered on $\eta $ Carinae. Because of the spatial resolution of the EPIC[*] camera and relatively short exposure time of the observations ($\sim $35 ks), they detected only 80 X-ray sources, most of them related to the massive OB-type stars with $L_{\rm x}$ $\sim $ 1032-1034  \ensuremath{{\rm erg~s^{-1}}}. Before the observation used here, three Chandra observations were obtained in this region, two (obsId 50 and 1249) in the timed exposure mode, and the third (obsId 51) in the continuous clocking mode, which produces no image. Using only observation obsId 1249, Evans et al. (2003) presented luminosities and hardness ratios of the hot stars in Tr 16, and part of Tr 14. Low-resolution X-ray spectra of luminous sources were discussed by Evans et al. (2004); however, the short exposure time of such an observation ($\sim $9.5 ks) was a serious limitation for the study of intermediate- and low mass stellar population in the region. This limitation exists even if obsId. 50 and obsId. 1249 are combined (Sanchawala et al. 2007), reaching completeness just at X-ray luminosity ($L_{\rm x}$) of $\sim $$\times $ 1031  \ensuremath{{\rm erg~s^{-1}}}, i.e. the X-ray emission level typical of single O- and early B-type stars.

In this paper we present results of the analysis of the deepest X-ray observation ever done in this region ($\sim $90 ks). Section 2 gives details on the observation and data-reduction procedures. Section 3 explains the method used to detect the sources, photon extraction and the construction of the catalog. In Sect. 4 we present results of the cross-correlation with existing near-IR and optical catalogs of objects and their characterization based on their color-color (CC) and color-magnitude (CM) diagrams. Section 5 presents a statistical study of variability in the X-ray domain. Section 6 deals with results of the analysis of extracted X-ray spectra. In Sect. 7 we discuss X-ray luminosities of stars and compare them statistically with the X-ray source population of ONC and Cygnus OB2 star-forming regions. In Sect. 8 we discuss X-ray and stellar parameter of O- and early B-type stars. Finally, in Sect. 9 we give a summary and draw conclusions of the paper.

   
2 The X-ray observations

Trumpler 16 was observed with the ACIS detector on board the Chandra X-ray Observatory (CXO) (Weisskopf et al. 2002) on 2006 August 31[*] (obsId 6402), as part of the Guaranteed Time Observation (GTO) Chandra program. The total effective exposure time was 88.4 ks. The data were acquired in VERY FAINT mode, to ease filtering of non-X-ray events, with six CCD turned on, the four comprising the ACIS-I array [0-3], plus CCDs 6 and 7, part of ACIS-S. However, data from the last two CCDs not are used in the following because of the degraded point spread function (PSF) and reduced effective area resulting from their large distance from the optical axis. The ACIS-I 17' $\times $ 17' field of view (FOV) is covered by 4 chips, each with 1024 $\times $ 1024 pixels (scale 0.49'' px-1). The observation was pointed toward RA = 10$^{\rm h}$44$^{\rm m}$47.93$^{\rm s}$ and Dec = -59$^\circ$43'54.21'', chosen to maximize the number of stars in the FOV and close to the optical axis, but also including most of the OB stars of the cluster. Figure 1-left shows Trumpler 16 as seen in X-rays by our ACIS-I observation.


  \begin{figure}
\par\includegraphics[width=8.25cm,clip]{9999f1A.ps}\includegraphics[width=8.15cm,clip]{9999f1B.ps}\end{figure} Figure 1: Left: color-coded ACIS-I image of the 17' $\times $ 17' field in Trumpler 16 (see color version in the electronic edition). Kernel smoothing was applied to highlight point sources. Energy bands for the RGB image are [0.5:1.5], [1.5:2.2], and [2.2:8.0] keV for the red, green, and blue colors, respectively. Contour lines show the spatial distribution of the dark V-shaped dust lane. Right: ACIS-I image composed with a H$_\alpha $ image. This image clearly show the sharp cut, along the V-shape structure, in the spatial distribution of X-ray sources.
Open with DEXTER

Most of the observed X-ray sources in Fig. 1-left are very likely located toward the central part of Trumpler 16, spatially constrained by the dark V-shaped dust lane of the Carina region (Brooks et al. 1998). However, note the small X-ray source concentration inside this dark structure, located towards the south east part of Fig. 1. Deep near-IR observation are needed to reveal counterparts of such population. We defer to a forthcoming paper the use of some public HST - WFPCII[*] observations to find signatures of star-environment interactions for some of the X-ray sources.

2.1 Data reduction

Data reduction, starting with the Level 1 event list provided by the pipeline processing at the CXO, was performed using CIAO 3.3.0.1[*] and the CALDB 3.1.0 set of calibration files. We produced a level-2 event file using the ACIS_PROCESS_EVENT CIAO task, taking advantage of the VF-mode enhanced background filtering and only retaining events with grades = 0, 2-4, 6 and status = 0. Photon energies were corrected for the time dependence of the energy gain using the CORR_TGAIN CIAO task. Intervals of background flaring were searched for, but none found. We hereafter assume a non-variable background. To improve the sensitivity to faint sources, given the spectrum of the background and that of typical sources, we filtered out events outside the [500:8000] eV energy band.

   
3 Data analysis

We built images in the three band-passes 0.5-1.5 keV (soft), 1.5-2.2 keV (medium), and 2.2-8.0 keV (hard). Before the color-coded image combination, we corrected single band-images by variations in exposure, sensitivity, and vignetting, by computing and dividing with its respective exposure maps. We constructed a color-coded X-ray image of the region by composition of the three soft (red), medium (green), and hard (blue) images (see Fig. 1). This image is a 17 $\times $ 17 arcmin field that comprises the center of Tr 16 and marginally the south east part of Tr 14. In addition to revealing a huge number of X-ray point-like sources with different color-energies, inspection of Fig. 1 suggests the presence of soft diffuse X-ray emission in the region. In particular, Townsley (2006) use a 57 ks ACIS-I observation centered on Tr 14 to study the diffuse X-ray emission in such a region and in the north west part of Tr 16. They explain the diffuse X-ray emission as likely arising from the fast O-star winds that shock and thermalize the surrounding medium. Because homogeneous data analysis techniques are required to do justice to a comparison of diffuse X-ray emission in this dataset with that of Tr 14, this study will be presented in a forthcoming paper.

   
3.1 X-ray source detection

Source detection was performed with the Palermo Wavelet Detection code, PWDetect[*] (Damiani et al. 1997b). It analyzes the data at different spatial scales, allowing the detection of both point-like and moderately extended sources, and efficiently resolving close source pairs. The most important input parameter is the detection threshold (SNR), which we established from the relationship between background level of the observation and expected number of spurious detections due to Poisson noise[*].

Diffuse, soft ($\sim $0.2-1.2 keV) X-ray emission has already been identified in this region (Townsley 2006; Corcoran et al. 1995) and causes different (non-uniform) X-ray background levels across the FOV of our observation. The best way of reducing the diffuse emission contribution is to discard soft photons in the detection procedure. We computed background (BKG) levels in four different energy ranges: 0.5-8.0 keV, 0.7-8.0 keV, 0.9-8.0 keV, and 1.1-8.0 keV. Since exposure maps are needed by the source detection algorithm, we used event files in these energy bands as input to the CIAO tool MKEXPMAP and assumed a monochromatic spectrum (kT=2.0 keV)[*]. The background level for each energy band was computed as the average of values measured in three circular regions free of sources. If we accept 10 spurious detections in the ACIS FOV[*], mean computed background levels can be translated into different SNR thresholds for each energy band, i.e.: 4.7, 4.65, 4.6, and 4.55. These input parameters for PWDetect result in different numbers of detected sources in each energy band, for instance, 1266 sources in the 0.9-8.0 keV band, more than detected in the ``canonical'' 0.5-8.0 keV band (1214 sources). This is consistent with the effect of diffuse soft X-ray emission in masking weak sources.

A careful visual inspection was performed on the 0.5-8.0 keV and 0.9-8.0 keV source lists. We manually rejected 242 and 271 detections respectively, considered spurious either because they were produced by different instrumental artifacts (e.g. CCD gaps, detector edges, false detections along the readout trails), or since they resulted from the ``fragmentation'' into discrete sources of the extended emission making a toroidal ring around the LBV star $\eta $ Carinae (Fig. 1). Furthermore, examination of afterglow[*] contamination led us to discard 29 and 20 false detections. The final source list was constructed by merging both catalogs: the list of 943 sources detected in the 0.5-8.0 keV band and the 1004 sources detected in the 0.9-8.0 keV band. Merging was performed using a criterion of maximum detection significance. The two catalogs coincide in 798 sources, consisting of 322 sources with Sig $_{\rm0.5-8.0}$ > Sig $_{\rm0.9-8.0}$, while for the remaining 476 sources Sig $_{\rm0.9-8.0}$ > Sig $_{\rm0.5-8.0}$ (where Sigx-y is the detection significance in the x-y energy band). Sixty-one sources were detected only in the 0.5-8.0 keV band, while 176 are exclusive to the 0.9-8.0 keV band and would have otherwise remained hidden because of the diffuse emission. This procedure led to a total number of 1035 X-ray sources in the entire field, which we analyze here.

   
3.2 Photon extraction

Table 1: Trumpler 16 X-ray source catalog (see the electronic version for the complete table).

Even with the high spatial resolution of the Chandra ACIS-I camera, the high source density in Trumpler 16, source photon extraction is not an easy task. Although circular regions would contain a relatively large fraction of the PSF for almost all source photons, the extended wings of the PSF mean that very large regions would be needed, incurring in the risk of contamination from nearby sources. Moreover, the resulting inclusion of a large number of background events would reduce the signal-to-noise of weak sources. On the other hand, extraction from regions that are too small may reduce the photon statistics for further spectral and timing analysis. To address these issues, we decided to use ACIS E XTRACT (AE) v3.79 (Broos et al. 2002), an IDL-based package that makes use of TARA[*], CIAO, and FTOOLS[*] software.

This task reduces the problem of accounting for non-Gaussian shapes of the local PSF by calculating the shape of the PSF model at each individual source's position. For some sources, the background level is affected by the extended PSF wings of the bright sources ($\eta $ Carinae and WR25) in the FOV. AE computes source background locally, by defining background extraction regions as circular annuli with inner radii 1.1 times the maximum distance between the source and the 99% PSF contour, and outer radii defined so that the regions contain more than 100 ``background'' events. To exclude contamination of the regions by nearby sources, background events are taken from a ``Swiss cheese'' image that excludes events within the inner annuli radii of all the 1035 sources.

The AE source extraction was performed using a PSF model that contains a specified fraction of source events ( $f_{\rm PSF}$). Generally, we choose $f_{\rm PSF}$ = 90%, and computed the contours from the PSF for a mono-energetic source with E = 1.49 keV. For 9.9% of the sources in the denser parts of the Trumpler 16 field this fraction was reduced to avoid contamination with other nearby sources, in the most extreme cases down to $f_{\rm PSF}$ $\sim $ 50% (just 3 sources).

Following AE science hints, we then refined the initial source positions computed by PWDetect[*] by correlating the source images with the model of local PSF computed by AE libraries. This procedure was only used for those sources lying at off-axis larger than 5 arcmin (316 sources), while for the rest of the source (719 sources) we simply adopted mean photon positions[*]. AE also estimates local background spectra, computes redistribution matrix files (RMFs) and auxiliary response files (ARFs), constructs light curves, performs Kolmogorov-Smirnov variability tests, and computes photometry in 14 different energy bands. Results of AE procedure appears in Table 1, which lists the source number in Col. (1); name according to CXC naming convention[*] (2); sky position (RA and Dec J2000) (3, 4) with relative uncertainty (5); off-axis angle ($\theta$) (6); significance of the detection (Sig.) from PWDetect analysis (7); the source extraction area (8); the PSF fraction within the extraction area, assuming E=1.49 keV (9); the background-corrected extracted source counts in the 0.5-8.0 keV band (NetCnts) (10); the count rates (CR = NetCnts/Exptime/PSF $_{\rm frac}$) in three spectral bands: 0.5-8.0 keV, 0.5-2.0 keV and 2.0-8.0 keV, (11-13); source photon quantiles at 25, 50, and 75% percent in Cols. 14 to 16 (see Sect. 3.3); the median photon energy ( $\overline E_{x}$) in (17). Column 18 is the log ( $P_{\rm ks}$) Kolmogorov-Smirnov probability of non-variability (see Sect. 5), and in (19) there are flags from PWDetect detection code.

   
3.3 X-ray hardness ratios

A commonly used tool for exploring the spectral properties of sources with low photon statistics is the hardness ratio (e.g. Prestwich et al. 2003; Schulz et al. 1989). In this conventional method, the full energy range is divided into two or three sub-bands and the detected source photons are counted separately in each band. Most popular definitions for a single hardness ratio (HR) exist on the basis of only two energy sub-bands: i- HR = H/S or ii- HR = (H-S)/(H+S). By these definitions, HR is very sensitive to small changes (i.e. statistical fluctuations) in the number of photons falling in each band. The requirement of total counts in the full energy band is at least 40 photons (just 36% of our sources satisfy this constraint). Above this limit, HR becomes a ``reliable method'' of estimating the real hardness of sources (Albacete Colombo et al. 2007).

An improved method to resolve this limitation is based on the quantile analysis (Hong et al. 2004). Instead of working with predetermined energy bands, we determine the energy Ex below which the net counts is x% of the total counts of the source. We define quantile Qx as $E^x-E_{\rm min} \over E_{\rm max}-E_{\rm min}$, where in our study $E_{\rm min}$ = 0.5 keV and $E_{\rm max}$ = 8.0 keV. We computed the median Q50 values and quartiles Q25 and Q75 and give values in Table 1. A minor inconvenience of this method is that for a given spectrum, various quantiles cannot be considered independent variables, unlike the counts in different energy bands. However, Hong et al. (2004) overcome this problem by considering the log (Q50/(1-Q50)) vs. 3(Q25/Q75) plane. Based on an extensive set of simulated spectra, they predict the loci of models in this plane. We used a set of absorbed thermal models with plasma temperatures of 0.2, 0.5, 1, 2, 4, 10 keV, and $N_{\rm H}$ equal to 1020, 1021, 0.4 $\times $ 1022, 1022, 4 $\times $ 1022, and 1023 cm-2. Note that the spectrum changes from soft to hard as one goes from left to right in the diagram (see Fig. 2).


  \begin{figure}
\par\includegraphics[width=7.4cm,clip]{9999f2.ps}\end{figure} Figure 2: Quantile color-color diagram (QCCD). The energy scale in the top X-axis shows the median energy values (Q50). The grid pattern represents the location of simulated spectra in the diagram (Hong et al. 2004). Note: filled circles and open boxes refer to X-ray sources with and without 2MASS counterparts (see next section), respectively. Small dots (red) correspond to sources with unreliable quantile values affected by poor photon statistics ($\leq $10 photons).
Open with DEXTER

In a statistical sense, a comparison between sources in the grid models suggest typical $N_{\rm H}$ and kT values distributed around 0.8 $\times $ 1022 cm-2 and 1.5 keV, but dispersed within grid models of $N_{\rm H}$ = 0.4-1.0 $\times $ 1022 cm-2 and kT= 1-3 keV, respectively. The absence of a systematic difference between the position of sources with and without a 2MASS identification on the hardness-ratio plane is probably caused by a combination of two observational bias: (i) related to the limiting magnitude of the 2MASS photometry ( $K_{\rm s}^{\rm limit}$ $\sim $ 14.3 mag); and (ii) the X-ray sensitivity of the Chandra data ( $f_{\rm x}^{\rm limit}$ $\sim $ 10-14.5 erg s-1 cm-2). However, it is not clear from the observational point of view, i.e. throughout the near-IR and X-ray flux source distributions, that a fraction of faint stellar X-rays sources would not have 2MASS counterparts. Deeper near-IR and X-ray observations are needed to unveil this issue.

   
4 Optical and near-IR analysis

   
4.1 Counterparts

The large amounts of gas, dust, and selective extinctions of the region, combined with the absence of available deep (V > 19 mag) optical photometry, prevent finding optical counterparts for most of our X-ray sources. We consider it appropriate to only use optical data to identify the massive star population (typical V < 12, Massey & Johnson 1993) of the region (i.e. Wolf-Rayet, O- and early B-types). We used a recent compendium of massive stars in the Carina Nebula presented by Smith (2006). Out of a total of 60 stars, 44 lie within the 0.0823 deg2 FOV of our X-ray observation. We added the binary (O5.5V + O9.5V) FO 15 (Niemela et al. 2006), not included in the list of Smith (2006). Our final list of 45 hot massive stars is comprised of 1 LBV, 1 Wolf-Rayet, 21 O-type stars, and 22 early B-type stars. With a cross-identification radius of 3 arcsec, a total of 28 X-ray sources were identified. All the O-type stars were identified as X-ray emitters, while just 7/22 ($\sim $31%) of early B-type (SpT. between B0V to B1.5V) have detectable X-ray emission. The list of massive stars with X-ray counterparts and the discussion of their X-ray properties is presented in Sect. 8.

Apart from the massive stars, our X-ray population is essentially composed of low-mass stars. In such a young cluster as Tr 16, low-mass stars are often optically invisible, being embedded and/or obscured by high amounts of gas and dust. We partially solve this problem by using the near-IR observations, on which the impact of dust extinction is reduced. We adopt J (1.25 $\mu$m), H (1.65 $\mu$m), and $K_{\rm s}$ (2.17 $\mu$m) photometry from the Two Micron All Sky Survey (2MASS) Point Source Catalog (PSC)[*]. The 2MASS is complete to magnitudes of 15.8, 15.1, and 14.3 mag in the J, H, and $K_{\rm s}$ bands, respectively. We restricted our photometry analysis to 2MASS sources with flag quality A-C, or D in at least one of the three magnitudes (see explanation in the 2MASS All-Sky Data Release User's Guide). With this restriction, 20 sources were removed from our initial list of 5938 sources in the ACIS FOV of our observation. This leaves a total of 5918 2MASS objects[*].

We began by cross identifying our X-ray source list with the 2MASS catalog. Identification radii, $R_{\rm id}$, were chosen to limit the number of spurious identifications due to chance alignments, $N_{\rm chance}$, and at the same time to include a large number of the true physical associations, $N_{\rm true}$. Identification radii used were 1.0, 1.5, 2.1, and 2.7 arcsec adopted for the four different off-axis angle ranges: [0-2), [2-4), [4-7), and >7 arcmin, respectively[*]. Results of the final identification are presented in Sect. 4.3 and shown in Table 2. The first seven columns are: X-ray source number; 2MASS nomenclature for identified sources; the offset between the two positions; J, H, $K_{\rm s}$ magnitudes; photometry quality flag (Ph.Q); confusion flag (Cont). A total of 660 X-ray sources out of the 1035 in our list were identified with 665 2MASS objects. Five X-ray sources (#96, #382, #401, #816, and #1034) were identified with two 2MASS counterparts each. After a visual inspection, we kept only the closer counterparts. Thus the final list of near-IR counterparts consists of 660 entries.

Table 2: Near-IR counterparts of Trumpler 16 X-ray sources. The complete version is available in the electronic version.

We estimate the expected number of extragalactic sources in our detection list by following the Flaccomio et al. (2006) procedure. We consider the ACIS count-rates of non-stellar sources in the Chandra Deep Field North (CDFN, Barger et al. 2003; Alexander et al. 2003) and estimate absorption corrected count-rates assuming $N_{\rm H}$ = 5 $\times $ 1021 cm-2 (from $A_{\rm v}$ $\sim $ 3.6, see Sect. 4.3) using PIMMS and assuming power-law spectra with index 1 and 2 (Giacconi et al. 2001). We then compare these count rates with upper limits taken at random positions in the ACIS FOV. For $\Gamma$ between 1 and 2 we obtain 72 to 95 expected extragalactic sources. Given the intrinsic near-IR fluxes of these sources and the absorption toward Trumpler 16, they are expected to be among the 385 without NIR counterparts (cf. Flaccomio et al. 2006). This means that no more than 18 to 24% of the unidentified X-ray population is related to extragalactic sources.

   
4.2 Unidentified X-ray sources

A large population of young stars, proto-stars, deeply embedded in dense circumstellar gas and dust should be present in the Carina Nebula (Smith et al. 2003). However optical and near-IR counterparts of young stellar objetcs (YSOs) are difficult to detect. Fortunately, X-ray emission is expected in YSOs along all their initial phases (Montmerle et al. 2000). X-ray photons easily escape from dense circumstellar material, where absorption process becomes important mainly for energies below 1.2 keV (Morrison & McCammon 1983b). This makes hard X-ray energies the most appropriate ``window'' to detect counterparts of deeply embedded young sources.

Of the 375 X-rays sources without near-IR counterparts, just a small fraction ($\sim $20%) is expected to be extra-galactic contamination (see Sect. 4.1). We are thus dealing with about $\sim $300 candidate young (first stage) low-mass stars, highly obscured by circumstellar material. In our data these sources typically have lower X-ray photon statistics than those with a near-IR counterpart, i.e. $\sim $8 vs. 28 average photons, respectively. In X-rays, no quantitative differences in the median energy and spectral quantiles were found for the X-ray sources with and without near-IR counterparts.

   
4.3 Near-IR properties of identified X-ray sources

We now investigate the near-IR properties of the X-ray sources. For this purpose we restrict our analysis to sources with high-quality photometry (Ph.Q = AAA) and no confusion (Cont. = 000). With these requirements the total number of IR sources in the ACIS FOV is reduced from 5918 to 2178. We have also set a further requirement on near-IR counterparts of X-ray sources, that their J, H, and $K_{\rm s}$ magnitude errors be all <0.1 mag. All these requirements yield 367 X-ray sources with good near-IR counterparts, out of the original 660.


  \begin{figure}
\par\includegraphics[width=8.25cm,clip]{9999f3.ps}\end{figure} Figure 3: JHK$_{\rm s}$ color-color diagram from high-quality 2MASS photometry. F ILLED CIRCLES and DOTS refer to X-ray detected and undetected 2MASS sources, respectively. A STERISKS show the source population with intrinsic $K_{\rm s}$ excess. Identified 2MASS source known as O- and early B-type stars are indicated by SQUARES. The main sequence is shown for reference. We also show the classical T Tauri stars (CTTS) locus of Meyer et al. (1997) and reddening vectors ( dashed lines) with length corresponding to $A_{\rm v}$ = 13 mag. Note: (i) the 2MASS photometry of $\eta $ Carinae is severely affected, with Phot.Qual. = ``DDD''; (ii) the peculiar position of the red box with $H-K_{\rm s}<$ -0.1 corresponds to the massive binary CPD-592628AB (SpT: O9.5V+B0.3V).
Open with DEXTER

Figure 3 shows the J-H vs. $H-K_{\rm s}$ color-color (CC) diagram for these AAA-flagged sources. We also plot for comparison the MS (Kenyon & Hartmann 1995), the classical T Tauri stars (CTTS) locus of Meyer et al. (1997), and reddening vectors starting from these loci and with slope ( $A_{K{\rm s}}/E(H-K_{\rm s})$ = 0.125) corresponding to the extinction law given by Hanson (2003). Trumpler 16 members with purely photospheric emission should lie in this reddening band. Otherwise, YSOs, such as CTTSs, and Herbig Ae/Be stars, because of the NIR excess emission originating in the inner parts of their circumstellar disks, are often found to the right of this band, i.e. in the CTTS locus. Fifty-one (out of 367) X-ray sources, i.e. likely Trumpler 16 members, have colors consistent with the (reddened) CTTS locus. This means a fraction of 51/339 (28 OB stars were discarded) $\sim $15% of all (low-mass) identified X-ray sources in the CC diagram. Of all 51 sources with intrinsic K-excess i.e. disk-star systems, eleven (Src-Id: 36, 41, 209, 230, 773, 966, 993, 996, 1002, 1003, and 1009) appear below the CTTS vectors, but with intrinsic bright $K_{\rm s}$ magnitudes, as is shown in both panels of Fig. 4. They are probably intermediate- to high-mass young stars with an intrinsic K-band excess that would be produced by massive accretion disks and/or extended envelopes surrounding massive YSOs. If confirmed, they will contribute about 40% of the total massive star population of the Trumpler 16 region.

Figure 4-left shows the $K_{\rm s}$ vs. $J-K_{\rm s}$ color magnitude (CM) diagram for the same stars plotted in Fig. 3. We also show for reference the expected cluster locus: the intrinsic $K_{\rm s}$ magnitudes and $J-K_{\rm s}$ colors for stars earlier than B5V were taken from the MS calibration of Knödlseder (2000) and Bessell & Brett (1989). For later spectral types (masses between 0.1 and 7 $M_\odot $), we adopted the 3 Myr isochrone from Siess et al. (2000), converted to the observational plane using the calibration given by Kenyon & Hartmann (1995). The adopted MS and 3 Myr isochrone overlap satisfactorily.

In order to estimate the typical visual absorption of cluster members, we computed the distance of each X-ray source from the cluster locus along the reddening direction on the $K_{\rm s}$ vs. J-H plane (see Fig. 4-right). Resulting $A_{\rm v}$ values for individual sources[*] are listed in Col. 9 of Table 2. Note that for 11.2 < $K_{\rm s}<$ 12.8 the absorption cannot be constrained because the reddening vector intersects the cluster locus more than once. In both CM diagrams, around 20 X-ray sources lie to the left of or relatively close to the reddened cluster loci. These are likely to be foreground MS stars, thus their $A_{\rm v}$ values, either cannot be computed or are close to zero. In a statistical sense, the median $A_{\rm v}$ value of OB stars ($A_{\rm v}$ = 2.0  $\pm$ 0.8 mag) is lower than what is computed for low-mass stars ($A_{\rm v}$ = 3.6  $\pm$ 2.4 mag, considered to be the typical absorption of the cluster). Obviously, the above estimates depend on the reliability of the assumed cluster locus and on the assumption that disk-induced excesses do not significantly affect the J and H magnitudes. As a caution against possible contamination and/or anomalous J magnitude, but in particular for the H-band, the median $A_{\rm v}$ for the low-mass star population was estimated by discarding sources with intrinsic $K_{\rm s}$ excesses (i.e. labeled with asterisks in Figs. 3 and 4). We must note that: (i) the dispersion along the $J-K_{\rm s}$ axis (see Fig. 4-left) indicates differential absorption of the region, and translates into a broad $A_{\rm v}$ distribution, with a spread 1 $\sigma \sim$ 2.4 mag; and (ii) the difference between median $A_{\rm v}$ values of low-mass and OB stars suggests a clearing effect of strong winds and the radiation field of massive stars on their surrounding environment. This conclusion has also been reached by Albacete Colombo et al. (2007) for massive stars of the Cyg OB2 region. However, compared to the Cyg OB2 region, the fraction of disk-stars members in the Trumpler 16 region is about four times greater.

Finally, we used the 2MASS J-band magnitudes to obtain an estimate of stellar masses for 510 counterparts (of a total 660 identified stars) with J-band Phot.Qual. = ``A'' to ``D''. We computed the mass vs. J mag relationship appropriate for the cluster mean age (3 Myr), distance (DM = 11.78 mag), and extinction ($A_{\rm v}$ = 3.6 mag), is with the last obtained as described above for the cluster locus in the CM diagrams[*]. We used Siess et al. (2000) models to compute PMS tracks of low- and intermediate-mass stars (masses $\leq $$M_\odot $). We interpolated the J mag - mass relation, using the J magnitude vs. mass relationship at 3 Myr (see Col. 10 of Table 2), to compute individual masses of stars[*]. Unfortunately, the relation mass-J mag degenerates in the mass ranges 0.2-0.65 $M_\odot $ (3 sources) and 2.7-4.52 $M_\odot $ (21 sources). We indicate these sources with ``Mass-deg.'' flag in Col. 11 of Table 2, and give mean mass values for sources lying in these two ranges, i.e. 0.42 and 3.55, respectively. Computed masses over 7 $M_\odot $ are potentially affected by large uncertainties in the extrapolated J mag - mass relation, and these values were excluded from further analysis. In Table 2 we give masses for a total of 510 stars, of which 410 range between 0.65 and 2.52 $M_\odot $.

   
5 X-ray variability


  \begin{figure}
\par\includegraphics[width=8.3cm,clip]{9999f4A.ps}\hspace*{3mm}
\includegraphics[width=8.3cm,clip]{9999f4B.ps}\end{figure} Figure 4: CM diagrams of the Trumpler 16 region. Symbols as in Fig. 4. The two parallel curves indicate the expected cluster loci for the assumed distance and no reddening, and for the mean reddening $A_{\rm v}$ = 3.6 mag. Masses are indicated on the left side of the unreddened isochrone. Identified OB stars show $K_{\rm s}$ magnitudes below 11. Stars with $K_{\rm s}$ excess show typically $J-K_{\rm s}>1$. Note: the limitation imposed by the 2MASS photometry ( $K_{\rm s}\sim $ 14.5 mag) corresponds to a mass limit of 1.1 $M_\odot $. Labels near OB stars refer to Simbad names ( left) and respective Spectral types ( right). The fact that stars Tr16-244 (O4If) and HD 93129AB (WN6ha+O4f) appear above the 120 $M_\odot $ track does not mean that they are more massive, since their evolutionary stages cannot be compared with a MS isochrone. $\eta $ Carinae was discarded because of the bad 2MASS photometry (Phot.qual. = ``DDD'', see text). Sources labeled with numbers indicate stars (with masses above 12 $M_\odot $) probably surrounded by disks.
Open with DEXTER


  \begin{figure}
\par\hspace*{1.5mm}\includegraphics[width=5.5cm,clip]{9999f5A.ps}...
...5D.ps}\hspace*{1.5mm}
\includegraphics[width=5.5cm,clip]{9999f5C.ps}\end{figure} Figure 5: Light curves (in the 0.5-8.0 keV band) showing different variability scenarios occurring among our 77 variable sources. Bin size is 900 s. The source number is indicated in each panel. Upper panels show a flare-like behavior with very quick rises and decay phases of hours, while bottom panels are slowly modulated variability. A probable flare+modulation variability is occurring in the last panel.
Open with DEXTER

PMS stars have high levels of X-ray activity that are commonly attributed to a ``scaled up'' solar-like corona formed by active regions. X-ray variability over a wide range of time scales is common in all magnetically active stars (e.g. Güdel 2004; Favata & Micela 2003; Feigelson & Montmerle 1999). On long time scales, this includes rotational modulation of active regions, their emergence and evolution, and magnetic cycles (e.g. Flaccomio et al. 2005; Marino et al. 2003). Most of the observed variations have short time scales ($\sim $hours), however, and can be attributed to small-scale flares triggered by magnetic reconnection events.

We first investigated X-ray variability in our sources using the non binned one-sample Kolmogorov-Smirnov (KS) test (Press et al. 1992). This test compares the distribution of photon arrival times with what is expected for a constant source. The test was applied to photons in the source's extraction regions, which also contain background photons. Given that the background was found to be constant with time (Sect. 2.1), the results, i.e. the confidence with which we can reject the hypothesis that the flux was constant during our observation, can be attributed to the source photons. Column 18 of Table 1, reports the logarithm of the KS-test significance with values <-4 truncated at that value: sources with $\log~(P_{\rm KS}$) < -3.0 can be considered almost definitely variable because we expect at most one of the 1035 sources (i.e. $\leq $0.1%) to be erroneously classified as variable. Seventy-seven X-ray sources ($\approx$7.4% of the total) fall in this category. Fifty-five sources with -2.0 < $\log~(P_{\rm KS}$) < -3.0 can be considered as likely to be variable, although about half of them suffer from low photon statistics. These numbers of sources are lower limits to the total number of variable sources in the region for several reasons: (i) most of the observed variability is in the form of flares, i.e. events that are shorter than our observation and with a duty cycle that may be considerably longer (Wolk et al. 2005); (ii) the sensitivity of statistical tests to time variability of a given relative amplitude depends critically on photon statistics (see Albacete Colombo et al. 2007). Hereafter, we consider those 77 sources with $\log~(P_{\rm ks}$) < -3 as variable.

To get a more accurate description of the detected variability, we extracted binned light-curves for each of the 77 variable sources in the region. We adopted a bin length of 900 s, a compromise between bins that are long enough to reach a good signal-to-noise ratio per bin for most sources and short enough to resolve the decay phase of typical flares. Since the background of our observation is both low (negligible for many sources) and constant in time, we did not apply any background subtraction to the presented light curves. In Fig. 5 we show examples of the different behaviors among the light-curves of variable sources. Source #254 like others (#35, #270, #438, #468, #480, #596, #761, and #975) experience ``impulsive'' flares with very quick rises and decay phases of only a few hours. Others (sources #503, #36, #41, #136, #164, #260, #489, #523, #731, and #811) show longer (2 to 10 h) flares. In several instances a second impulsive event is visible during the exponential decay of a previous flare (e.g. sources #36, #52, #87, #196, #564, #600, #620, and #623). The case of source #696 is a combination of both variability types, with two consecutive flares. Other sources like #714 and #651 (e.g. #27, #251, #271, #793, #839, #890, #904 and #980) show light curves that, instead of showing typical flares, are characterized by slow continuous rises or decays that might be explained by rotational modulation of non-homogeneously distributed plasma (Flaccomio et al. 2005). Finally, light-curves, like those observed for sources #524, #71, and #933, seem to be related to a combination of flare like activity and rotational modulation.

5.1 Variability in massive stars

Because X-ray emission from O stars, which is believed to be unrelated to solar-like magnetic activity, comes from the integrated emission from many small shocks randomly occurring in their strong winds (Feldmeier et al. 1997; Owocki & Cohen 1999), on average, global X-ray variability is not expected to occur. However, it is surprising that three (out of 28) massive stars, namely Tr 16-11 (B1.5 V, source #136), Tr 16-5 (B1 V, source #489), and the binary HD 93205 (O3.5V((f))+O8V, source #242), are significantly variable, with $\log~(P_{\rm KS}$) values lower than -3. The origin of the observed flare-like variability in the first two sources (B-type stars) is probably coronal activity of unresolved late-type companions. This hypothesis was also proposed for Tr16-11 itself by Evans et al. (2003).

Figure 6 shows the X-ray light-curve of the HD 93205 binary system. Two different processes may be acting simultaneously to explain the observed variability:

We note that the origin of X-ray variability in high mass stars is beyond of the scope of this work, because it needs a larger X-ray data set. We leave this subject to a further investigation.

   
6 Spectral analysis

To characterize the hot plasma responsible for the X-ray emission of Trumpler 16 stars and to estimate their intrinsic X-ray luminosities, we analyzed the ACIS spectra of the 615 (out of all 1035) sources with more than 20 net photons (NET_CNTS), i.e. corrected for local background. Spectral parameters for sources with less than 20 net photons are much too ill-constrained (Albacete Colombo et al. 2007), and thus were not determined. Moreover, high local background could affect the reliability of computed spectral parameters. We defined the fraction  $f_{\rm cont}$ as the ratio between BKG_CNTS (number of background counts in the source extraction region) and NET_CNTS. We accepted spectral fits for sources with $f_{\rm cont}$ < 1 and NET_CNTS $\geq$ 20 photons. Of all 1035 sources, only 563 satisfy both conditions above, while 119 show $f_{\rm cont}$ < 1, and 353 lie in the low statistics regime (i.e. NET_CNTS < 20 ph).

Source and background spectra in the 0.5-8.0 keV band were produced with AE (see Sect. 3.2), along with individual ``redistribution matrices files'' (RMF) and ``ancillary response files'' (ARF). For model-fitting, spectra were grouped to have a specified number of events in each energy bin. Grouping was tuned to the source statistics, and we chose 2, 5, 7, 10, and 60 counts per channel for sources with net-counts in the following ranges: [20-40], [40-100], [100-200], [200-500], and [500-**]. Spectral fitting of background-subtracted spectra was performed with XSPEC v12.0 (Arnaud 2004) and our own shell and TCL scripts to automate the process as described in Flaccomio et al. (2006). Because background-corrected spectra are not appropriately handled by C-statistics (Getman et al. 2005), best-fit parameters for the chosen models were computed by chi-squared ($\chi^2$) minimization.

We fit our spectra by assuming emission by a thermal plasma, in collisional ionization equilibrium, as modeled by the APEC code (Smith et al. 2001). Elemental abundances are not easily constrained with low-statistics spectra and were fixed at Z=0.3 $Z_\odot$, with solar abundance ratios taken from Anders & Grevesse (1989). The choice of sub-solar abundances is suggested by several X-ray studies of star-forming regions (e.g. Preibisch 2003; Feigelson et al. 2002). Absorption was accounted for using the WABS model, parameterized by the hydrogen column density, $N_{\rm H}$ (Morrison & McCammon 1983a). In Table 3 we give best-fit parameters ( $\chi_\nu^2$, $N_{\rm H}$, kT and $L_{\rm x}$) of the sources.

Table 3: X-ray spectral fits of Trumpler 16 sources. The complete version is available in the electronic version.

Except for 28 massive O- and early B-type stars, we fit source spectra with one-temperature (1T) plasma models using an automated procedure. To reduce the risk of finding a local minimum in the $\chi^2$ spaces, our procedure chooses the best fit among several obtained starting from a grid of initial values of the model parameters: $\log~(N_{\rm H})$ = 21.0, 21.7, 22.0, 22.4, 22.7 and 23.0 cm-2 and kT = 0.5, 0.75, 1.0, 2.0, 5.0 keV. The best-fit values of $\log~(N_{\rm H})$ < 20.3 cm-2 were truncated at 20.3 for two cases (sources #150 and #944) because, in the 0.5-8.0 keV energy range, ACIS spectra are insensitive to lower column densities. In a similar way, above 10 keV Chandra is not able to distinguish between such high temperatures. Therefore, 74 best-fit values of kT above 8 keV were truncated to that value. They are indicated with a flag HARD-TAIL in Table 3.

Figure 7 shows the distribution of best-fit $\log~(N_{\rm H})$ values for the 563 fitted sources. They appear to be normally distributed with a median $\log N_{\rm H}$ $\sim $ 21.73 ($N_{\rm H}$ = 5.37 $\times $ 1021 cm-2) and an FWHM of $\sim $0.4 dex. The log-normal distribution of the sources is indicated with the Gaussian curve[*]. The computed median of $N_{\rm H}$ ($\sim $5.37 $\times $ 1021 cm-2) is converted to a median $A_{\rm v}$ = 3.35 by use of the Vuong et al. (2003) relation: $N_{\rm H}$/$A_{\rm v}$ = 1.6 $\times $ 1021 atoms cm-2 mag-1. We also tested the relation $A_{\rm v}$ = 0.56 $N_{\rm H}$ + 0.23 [$N_{\rm H}$ in 1021] (Predehl & Schmitt 1995), for which the median $A_{\rm v}$ is 3.23 mag. Both of these values are in good agreement with the median $A_{\rm v}$ = 3.6 mag computed from our near-IR analysis (see Sect. 4.3). The 1$\sigma $ dispersion of the $N_{\rm H}$ distribution is 0.4 dex. It is translated into typical a $A_{\rm v}$ range between 1.3 to 6.7 mag of visual extinction. The X-ray sources without near-IR counterparts seem to be distributed towards higher absorption values (median $\log~(N_{\rm H})$ $\sim $ 21.9 cm-2) with respect to those with near-IR counterparts. Unfortunately, they generally have poor X-ray photon statistics and consequently a less reliable estimation of their X-ray spectral parameters.


  \begin{figure}
\par\includegraphics[width=8.5cm,angle=0]{9999f7A.ps}\includegraphics[width=8.5cm,angle=0]{9999f7B.ps}\end{figure} Figure 7: Left: the solid histogram shows the $N_{\rm H}$ ( $A_{\rm v}^{\rm xspec}$) distribution of 415 spectral fitted sources with 2MASS counterparts. The upper scale was converted using an $N_{\rm H}$/$A_{\rm v}$ ratio of 1.6 $\times $ 1021. The thick dashed histogram shows the $A_{\rm v}^{\rm near-IR}$ mag distribution for all near-IR counterparts as was computed in Sect. 4. Both distributions peak at $A_{\rm v}^{\rm xspec}$ = 3.35 (1$\sigma $ = 0.4) mag and $A_{\rm v}^{\rm near-IR}$ = 3.26 mag (1$\sigma $ = 0.3 dex). The shadowed histogram shows the $N_{\rm H}$ distribution of 148 sources with no near-IR counterpart. Right: same as the left panel for plasma temperatures (kT) with, in addition, the distribution for ``flaring'' sources (black-filled histogram). The peak of the overal distribution is at $\sim $1.6 keV.
Open with DEXTER

Unlike the $N_{\rm H}$ distribution, the kT distribution of plasma temperatures is not log-normal. It peaks at $\sim $1.6 keV, has a median $\sim $1.95 keV, and shows an extended hard tail attributed both to variable sources ( $\log~P_{\rm KS}<-3$) with harder spectra (median kT = 3.25 keV, as expected from coronal heating processes involved in flare-like activity) and to highly absorbed sources (i.e. $A_{\rm v}$ $\geq$ 6 mag) showing a median kT $\sim $ 2.6 keV (while those in the range 1 $\leq $ $A_{\rm v}$ $\leq $ 6 mag are softer and distributed with a median kT=1.75 keV).

Finally, suspected single massive stars show typically soft spectra with median kT=0.62 keV, while for known massive binaries, this value rises to kT=2.1 keV, no doubt due to hard X-ray photons being produced in the colliding wind region (CWR) of the massive O+OB binaries. Details of X-ray spectral characteristics of massive stars in the region are presented in Sect. 8.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{9999f8.ps}\end{figure} Figure 8: The X-ray luminosity distribution for 592 low-mass stars computed from a single $L_{\rm x}$ count-rate conversion factor is shown as a solid line. The grey filled histogram shows the $L_{\rm x}$ distribution for 354 X-ray sources without near-IR counterpart. The 77 variable sources are shown with the black-filled histogram. Upper-right inset (Hot stars) shows the unabsorbed $L_{\rm x}$ from X-ray spectral fits (see Sect. 8). Error bars are 1$\sigma $ Poisson errors.
Open with DEXTER

   
7 X-ray luminosity of stars

Unabsorbed X-ray luminosities were computed for those sources with available spectral fits for the [0.5-8.0] keV energy range. For sources with no available and/or reliable spectral fit (119+353 out of 1035) $L_{\rm x}$ were computed using a single count rate to $L_{\rm X}$ conversion factor (CF)[*]. CF in Trumpler 16 is 8.3 10.35.2 $\times $ 1033 erg/ph. Upper and lower values correspond to $\pm$ 1$\sigma $ uncertainties, respectively.

In Fig. 8 we show the X-ray luminosity distribution for low-mass[*] stars, while known OB stars are plotted separately in the upper inset histogram. The $L_{\rm x}$ distribution of sources has been plotted separately for (i) 592 X-ray sources with near-IR counterpart, not including 28 OB stars and variable sources with a median log ($L_{\rm x}$$\sim $ 30.5  \ensuremath{{\rm erg~s^{-1}}}. The peak of the distribution indeed marks the completeness limit of our X-ray observation. (ii) There are 354 unidentified sources, except for variable ones, which appear systematically less luminous than those with a 2MASS counterpart, with median log ($L_{\rm x}$$\sim $ 30.2  \ensuremath{{\rm erg~s^{-1}}}. (iii) Seventy seven variable sources, showing a median log ($L_{\rm x}$$\sim $ 31.1  \ensuremath{{\rm erg~s^{-1}}}, i.e. about 4 times higher than the observed $L_{\rm x}$ for similar stars in a quiescent (non-flaring) phase. The upper inset in Fig. 8 shows the unabsorbed $L_{\rm x}$ for 28 massive stars (including 13 binary systems) in the region. By the following results and discussion presented in Sect. 8, binaries appear with typical $L_{\rm x}$ over 10+32  \ensuremath{{\rm erg~s^{-1}}}, higher than observed for those suspected single B-type stars. As we discuss in Sect. 8, massive binaries have typical X-ray luminosities $L_{\rm x}$ $\geq$ 10+32  \ensuremath{{\rm erg~s^{-1}}}, higher than those of single OB stars.

We examined how X-ray activity depends on stellar mass for Trumpler 16 low-mass stars, and compares the results with those already known for the ONC and Cyg OB2 SFRs. To increase statistics, we used $L_{\rm x}$ values computed by means of the average CF. While $L_{\rm x}$ and mass for Cyg OB2 stars were computed following the same procedures used here (Albacete Colombo et al. 2007), for ONC stars they have been computed differently in the literature[*]. For the sake of homogeneity, we re-computed $L_{\rm x}$ and masses of ONC stars by using a single count rate to $L_{\rm x}$ conversion factor (CF $_{\rm
onc}$ = 7.52 $\times $ 10+32 erg/ph) and 2MASS photometry, respectively.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{9999f9.ps}\end{figure} Figure 9: X-ray luminosity vs. stellar mass for stars in the near-IR sample with masses determined from the PMS models of Siess et al. (2000). Tr16 stars are plotted by open and filled circles, corresponding to $L_{\rm x}$ from CF and spectral fits, respectively. Crosses indicate Cyg OB2 sources (Albacete Colombo et al. 2007), while small boxes are ONC sources (Preibisch et al. 2005). The thick gray line shows the linear regression fit to the Tr16 low-mass (0.7-2.5 $M_\odot $) stars.
Open with DEXTER

In Fig. 9 we show a plot of $L_{\rm x}$ vs. star mass. We used circles for all 510 Tr16 sources with estimated masses. Filled (329) and open (181) circles indicate $L_{\rm x}$ values computed from spectral fits and using the CF, respectively. We performed a linear regression for all sources in the 0.7-2.5 $M_\odot $ range: log ($L_{\rm x}$) = 30.26($\pm$0.11) + 1.0($\pm$0.09)log ($M/M_\odot $) with a standard deviation in the residuals of 0.38 dex. The power-law slope we find here agrees with the one found by Albacete Colombo et al. (2007) for the Cyg OB2 region: log ($L_{\rm x}$) = 30.33($\pm$0.16) + 0.71($\pm$0.13)log ($M/M_\odot $) for masses in the 0.5-3.0 $M_\odot $ range. This slope also agrees with what we find for the ONC stars in the mass range 0.5-3.0 $M_\odot $, namely 0.82  $\pm$  0.09.

Changes in the X-ray activity of stars with different ages of the SFRs has been previously reported by Preibisch & Feigelson (2005). To address this issue, we computed $L_{\rm x}$ detection limits for low-mass stars in the Tr16 (3 Myr), CygOB2 (2 Myr), and ONC (1 Myr) observations, as log ($L_{\rm x}$$\sim $ 30.5, 30.3 and 28.5  \ensuremath{{\rm erg~s^{-1}}}, respectively, above which a source is detectable anywhere in the FOV. A mass-dependent completeness fraction ( $f_{\rm comp}$) for our survey of Trumpler 16 stars was then computed by adopting the X-ray luminosity functions given by Preibisch & Feigelson (2005): $f_{\rm comp}$ is $\sim $5% for 0.1-0.5 $M_\odot $, $\sim $40% for 0.5-0.9 $M_\odot $, and $\sim $55% in the 0.9-1.2 $M_\odot $ range. Our survey of Tr 16 is statistically complete for masses $\geq$1.5 $M_\odot $. In Fig. 10 we present $L_{\rm x}$ cumulative distributions for Trumpler 16, Cyg OB2, and ONC stars in the mass range 1.5-2.5 $M_\odot $ where all three data sets are complete. While Trumpler 16 and Cyg OB2 show very similar distributions, the ONC looks quite different from to Tr16. The distance between two distributions (D=0.31) was computed using a two-sample Kolmogorov-Smirnov test (Press et al. 1992). We are able to confirm, with a probability $\geq$99.9% ( $\log~(P_{\rm ks})$ = -4.15), that young ONC (1 Myr) stars, with masses in the range 1.5-2.5 $M_\odot $, are intrinsically more luminous in X-rays than their Tr16 (3 Myr) counterparts.


  \begin{figure}
\par\includegraphics[width=8.4cm,clip]{9999f10.ps}\end{figure} Figure 10: Cumulative distributions of X-ray luminosities for Trumpler 16 low-mass stars with masses 1.5 $\leq $ $M/M_\odot $ $\leq $ 2.5 (thick solid line). Dotted and dashed lines represent the $L_{\rm x}$ cumulative distributions for ONC and CygOB2, respectively, in the same mass range.
Open with DEXTER

   
8 X-rays from massive stars

Table 4: Stellar parameters and X-ray spectral results for massive stars.

A variety of different physical mechanisms are responsible for the observed X-ray emission in OB-type stars. The most widely accepted explanation invokes multiple small-scale shocks in the inner layers of their radiation-driven stellar winds (e.g. Feldmeier et al. 1997). In recent years has gained importance a plasma heating model known as magnetically channeled wind shock (MCWS) (Schulz et al. 2003; Owocki et al. 2005; ). Moreover, WR+OB and/or O+OB interacting wind binary systems produce an excess of X-ray emission from the CWR (Zhekov & Skinner 2000; Stevens et al. 1992; Pittard & Stevens 2002).

Given the relatively large number of massive OB-type stars in the region, it is most relevant to show in Table 4 most appropiate stellar parameters of the massive stars ( \ensuremath{\mathit{L}_{\rm bol}}, \ensuremath{\dot{M}}, and $L_{\rm wind}$ computed as ${1 \over
2}\ensuremath{\dot{M}} $ $v_\infty^2$)[*]. We computed $L_{\rm x}$ by means of an absorbed ( WABS) thermal plasma model ( APEC). Metal abundance was fixed at Z = 0.3 $Z_\odot$ in fitting faint sources (<100 ph.), while it was left as a free parameter for the remaining cases (see Table 4). In Fig. 11 we show the $L_{\rm x}$/ $L_{\rm bol}$ relation. In spite of the observed scatter, one sees that the median $L_{\rm x}$/ \ensuremath{\mathit{L}_{\rm bol}} for binaries is about 7 times larger ($\sim $8.3 $\times $ 10-7) than that of suspected single stars ($\sim $1.1 $\times $ 10-7)[*]. Four of the 13 known binaries in the analysis are well-separated systems: the LBV $\eta $ Carinae, a probable long period binary with $P\sim$ 2026  $\pm$ 2 days (Daminelli et al. 2007), the Wolf-Rayet star WR25 (HD 93162), a 208-days period binary system (Gamen et al. 2006), and HDE 303308, resolved as a binary system with a component separation of about $\sim $38 AU projected along the fine guidance sensor (FGS) y-axis[*]. The stars Tr16 #9 (O9.5 V+?) and Tr16 #23 (O7 V+?) show photometric variability probably related to a secondary component (Nelan et al. 2004), we also consider them to be binaries.

Besides the several works (e.g. Evans et al. 2003; Sanchawala et al. 2007; Seward & Chlebowski 1982; Albacete Colombo et al. 2003; Corcoran et al. 1995; Seward et al. 1979; Evans et al. 2004) about the origin of the observed X-ray emission on massive stars of this region, this goal is far beyond the scope of our paper and so extensively discussed here.

   
9 Summary and conclusion

We report here results of a deep Chandra X-ray observation pointed toward the $\sim $3 Myr old star-forming region Trumpler 16. Source detection was performed using the PWDdetect code, identifying 1035 X-ray sources in the 17' $\times $ 17' ACIS-I FOV. Most of these seem to be outside the obscured V-shaped region of dust and gas. Star formation in this part of the masked region has probably been disrupted and/or diminished as the stellar winds are blocked inside the cloud, due to the efficiency of the photo-evaporation processes caused by interactions of the nearby hot massive stars with the dense dust and gas structures.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{9999f11.ps}\end{figure} Figure 11: Log ($L_{\rm x}$/ $L_{\rm bol}$) versus log ( $L_{\rm bol}$) relation for massive stars in Tr16. Open boxes refer to suspected single stars and filled circles to known binaries. Notes: i- because of flare-like variability (see Sect. 5.1) observed in Tr16-#5 (B1.5V) and Tr16-#11 (B1V), their X-ray emission is probably dominated by a unknown lowmass star companion. ii- The star Tr16-#22 (O8.5V+?) has the highest log ($L_{\rm x}$/ $L_{\rm bol}$) ratio ($\sim $-5.56) among all MS OB-type stars of the region.
Open with DEXTER

Data extraction was performed using the semi-automated IDL-based ACIS E XTRACT package, which is well-suited to the analysis of observations of crowded fields such as ours.

The X-ray source list was cross-identified with optical and near-IR (2MASS) catalogs: 28 X-ray sources (of 44 within the FOV) were identified with optically characterized OB members of Trumpler 16 and 760 with 2MASS sources. Among these sources, almost all are believed to be Trumpler 16 members. About 90 X-ray sources without optical/NIR counterparts are estimated to have of extragalactic nature (AGNs), while the remaining X-ray sources with no counterpart are likely to be associated with members that are fainter than the 2MASS completeness limit.

To characterize the previously unidentified likely cluster members with NIR counterparts, we placed them on NIR color-magnitude ($K_{\rm s}$ vs. $H-K_{\rm s}$) and color-color ( $H-K_{\rm s}$ vs. J-H) diagrams. A first estimate of interstellar extinction was obtained by adopting a 3 Myr isochrone for the low- and intermediate-mass stars and assuming that O- and early B-type stars lie on the MS. We find a median visual absorption for OB stars of $A_{\rm v}$ $\sim $ 2.0 mag, while low mass likely members seem to be slightly more absorbed, $A_{\rm v}$ $\sim $ 3.6 mag. We also use the 3 Myr isochrone and the J magnitude to estimate masses of likely members assuming that they share the same distance and absorption. Our sample of X-ray selected members with near-IR counterparts reaches down to M = 0.5-0.6 $M_\odot $, and is very likely complete down to $\sim $1.5 $M_\odot $. From the $H-K_{\rm s}$ vs. J-H diagram we estimate that $\sim $15% (51/339) of low-mass stars have NIR excesses, finding it to be quite a high percentage with respect to the 2 Myr old Cyg OB2. We believe that the disk fraction in young SFRs is more dependent on the spatial morphology of gas and dust around massive stars, which may enhance photo-evaporation, and thus shorten disk lifetimes, than on the total number of massive stars in the region.

At least 77 sources, i.e. $\sim $7.4%, were found to be variable within our observation with a confidence level greater than 99.9%. Only three of the 28 detected O- and early B-type stars were detected as variable during our 90-ks observation, in spite of the high statistics of the OB stars' light curves. These exceptions are the known binary O3.5V+O8V star HD 93205 (our source #242), showing a rather linear decay of the count rate during the observation plus a short-term variability, the B1.5 V star Tr16-11 (source #136), and the B1V star Tr16-5 (source #489). The last two show a flare-like variability probably related to unresolved low-mass companions.

We modeled the ACIS X-ray spectra of sources with more than 20 photons and $f_{\rm cont}$ < 1. We assumed an absorbed single-component thermal emission model. The median $\log~(N_{\rm H})$ of the sources is 21.73 (cm-2). This value agrees well with the median $A_{\rm v}$ computed from the near-IR diagram. The median kT of low-mass stars is 2.6 keV. Sources associated with O- and early B-type stars are instead quite soft (median kT: 0.60 keV). Absorption-corrected X-ray luminosities of OB stars were calculated from the best-fit spectral models. O and B-type stars are the most luminous, with $L_{\rm x}$ = 2.5 $\times $ 1030-6.3 $\times $ 1033 erg s-1. Their X-ray and bolometric luminosities roughly agree with the relation $L_{\rm x}$/ $L_{\rm bol}=10^{-7}$, albeit with an order of magnitude dispersion. Low-mass stars have $L_{\rm x}$ ranging between 1030and 1031  \ensuremath{{\rm erg~s^{-1}}} (median $L_{\rm x}$ = 2.8 $\times $ 1030). Variable low mass stars are on average 0.5 dex brighter (log ($L_{\rm x}$$\sim $ 31.0  \ensuremath{{\rm erg~s^{-1}}}). These X-ray luminosities are consistent with those of similar mass stars in the slightly younger (2 Myr) Cyg OB2 region. However, in the mass range 1.5-2.5 $M_\odot $, the ONC (1 Myr) shows a higher X-ray activity level than observed in Trumpler 16 stars in the same mass range. We believe that the age-$L_{\rm x}$ activity connection is an acceptable explanation of this result.

Acknowledgements
We thank the referee, Marc Gagne, for his time and many useful comments that improved this work. This publication makes use of 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 the National Aeronautics and Space Administration and the National Science Foundation. J.F.A.C. acknowledges support by the Marie Curie Fellowship Contract No. MTKD-CT-2004-002769 of the project ``The Influence of Stellar High Energy Radiation on Planetary Atmospheres'', and the host institution INAF - Osservatorio Astronomico di Palermo (OAPA). J.F.A.C. is a researcher member of the Consejo Nacional de Investigaciones Científicas y Tecnológicas (CONICET)-Argentina and acknowledges support from this institution. G.M., F.D. and S.S. acknowledge financial support from the Ministero dell'Università e della Ricerca research grants, and ASI/INAF Contract I/023/05/0.

References

 

  
Online Material

Table 1: Trumpler 16source catalog.

Table 2: Near-IR counterparts of Trumpler 16 X-ray sources.

Table 3: X-ray source spectroscopy.



Copyright ESO 2008