A&A 490, 1055-1070 (2008)
DOI: 10.1051/0004-6361:200809999
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 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 values. The cluster median extinction is
= 3.6 mag, while OB-type stars appear less absorbed, having
= 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
5.37
1021 cm-2 and
1.95 keV. The OB stars appear, softer with a median
0.65 keV. The median X-ray luminosity is 6.3
1030
,
while variable sources show a larger median
value of 13
1030
.
OB-stars have an even higher median
of 80
1030
,
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
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
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
105 solar masses (Grabelsky et al. 1988). The concentration of massive stars (i.e.
20
)
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) 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=
/E(B-V) values (3.48 and 4.16 for Tr 16 and Tr 14 regions). A more reliable distance (2250
180 pc), was derived from proper motion and Doppler velocities of the expanding
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
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 Carinae. Because of the spatial resolution of the
EPIC
camera and relatively short exposure time of the observations (
35 ks), they detected only 80 X-ray sources, most of them related to the massive OB-type
stars with
1032-1034
.
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 (
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 (
)
of
7
1031
,
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 (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.
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'
17' field of view (FOV) is covered by 4 chips, each with 1024
1024 pixels (scale 0.49'' px-1). The observation was pointed toward RA = 10
44
47.93
and Dec = -59
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.
![]() |
Figure 1:
Left: color-coded ACIS-I image of the 17' ![]() ![]() |
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.
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.
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
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.
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 (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 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
> Sig
,
while for the remaining 476 sources Sig
> Sig
(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.
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 ( 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 (
). Generally, we choose
= 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
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 (
) (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
)
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 (
)
in (17). Column 18 is the log (
)
Kolmogorov-Smirnov probability of non-variability
(see Sect. 5), and in (19) there are flags from PWDetect detection code.
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
,
where in our study
= 0.5 keV and
= 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
equal to 1020, 1021, 0.4
1022, 1022, 4
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).
![]() |
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 (![]() |
Open with DEXTER |
In a statistical sense, a comparison between sources in the grid models suggest typical
and kT values distributed around 0.8
1022 cm-2 and 1.5 keV, but dispersed within grid models of
= 0.4-1.0
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 (
14.3 mag); and (ii) the X-ray sensitivity of the Chandra data (
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.
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 (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 m), H (1.65
m), and
(2.17
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
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,
,
were chosen to limit the number of spurious identifications due to chance alignments,
,
and at the same time to include a large number of the true
physical associations,
.
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,
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 = 5
1021 cm-2 (from
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
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.
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 (20%) is expected to be extra-galactic contamination (see Sect. 4.1). We are thus dealing with about
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.
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.
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 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.
![]() |
Figure 3:
JHK![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Figure 3 shows the J-H vs.
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 (
= 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)
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
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
vs.
color magnitude (CM) diagram for the same stars plotted in Fig. 3. We also show for reference the expected cluster locus:
the intrinsic
magnitudes and
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
), 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 vs. J-H plane (see Fig. 4-right). Resulting
values for individual sources
are listed in Col. 9 of
Table 2. Note that for 11.2 <
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
values, either cannot be computed or are close to zero. In a statistical sense, the median
value of OB stars (
= 2.0
0.8 mag) is lower than what is computed for low-mass stars (
= 3.6
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
for the low-mass star population was estimated by discarding sources with intrinsic
excesses (i.e. labeled with asterisks in Figs. 3 and 4). We must note that: (i) the dispersion along the
axis (see Fig. 4-left) indicates differential absorption of the
region, and translates into a broad
distribution, with a spread 1
2.4 mag; and (ii) the difference between median
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 ( = 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
7
). 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
(3 sources) and 2.7-4.52
(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
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
.
![]() |
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 ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
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 (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
) < -3.0 can be considered almost definitely variable because we expect at most one of the 1035 sources (i.e.
0.1%) to be erroneously classified as variable. Seventy-seven
X-ray sources (
7.4% of the total) fall in this category. Fifty-five sources with -2.0 <
) < -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
) < -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.
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
) 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:
![]() |
Figure 6:
X-ray light-curve of the massive binary (O3.5V((f))+O8V) HD 93205. This light curve shows: i- a uniform decay probably related to the orbital motion of the system, with changing absorption ![]() |
Open with DEXTER |
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
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
< 1 and NET_CNTS
20 photons. Of all 1035 sources, only 563 satisfy both conditions above, while 119 show
< 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 () 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 ,
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,
(Morrison & McCammon 1983a). In Table 3 we give best-fit parameters (
,
,
kT and
)
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 spaces,
our procedure chooses the best fit among several obtained starting
from a grid of initial values of the model parameters:
= 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
< 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
values for the 563 fitted sources. They appear to be normally distributed with a median
21.73 (
= 5.37
1021 cm-2) and an FWHM of
0.4 dex. The log-normal distribution of the sources is indicated with the Gaussian curve
. The computed median of
(
5.37
1021 cm-2) is converted to a median
= 3.35 by use of the Vuong et al. (2003) relation:
/
= 1.6
1021 atoms cm-2 mag-1. We also
tested the relation
= 0.56
+ 0.23 [
in 1021] (Predehl & Schmitt 1995), for which the median
is 3.23 mag. Both of these values are in good agreement with the median
= 3.6 mag computed from our near-IR analysis (see Sect. 4.3). The
1
dispersion of the
distribution is 0.4 dex. It is translated into typical a
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
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.
![]() |
Figure 7:
Left: the solid histogram shows the ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Unlike the
distribution, the kT distribution of plasma temperatures is not log-normal. It peaks at
1.6 keV, has a median
1.95 keV, and shows an extended hard tail attributed both to variable sources (
)
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.
6 mag) showing a median kT
2.6 keV (while those in the range
1
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.
![]() |
Figure 8:
The X-ray luminosity distribution for 592 low-mass stars computed from a single ![]() ![]() ![]() ![]() |
Open with DEXTER |
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)
were computed using a single count rate to
conversion factor (CF)
. CF in Trumpler 16 is 8.3
10.35.2
1033 erg/ph. Upper and lower values correspond to
1
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
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 (
)
30.5
.
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 (
)
30.2
.
(iii) Seventy seven variable sources, showing a median
log (
)
31.1
,
i.e. about 4 times higher than the observed
for similar stars in a quiescent (non-flaring) phase. The upper inset in Fig. 8 shows the unabsorbed
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
over 10+32
,
higher than observed for those suspected single B-type stars. As we discuss in Sect. 8, massive binaries have typical X-ray luminosities
10+32
,
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
values
computed by means of the average CF. While
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
and masses of ONC stars by
using a single count rate to
conversion factor (CF
= 7.52
10+32 erg/ph) and 2MASS photometry, respectively.
![]() |
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 ![]() ![]() |
Open with DEXTER |
In Fig. 9 we show a plot of vs. star mass. We used
circles for all 510 Tr16 sources with estimated masses. Filled (329)
and open (181) circles indicate
values computed from spectral fits
and using the CF, respectively. We performed a linear regression for
all sources in the 0.7-2.5
range: log (
) = 30.26(
0.11) + 1.0(
0.09)log (
)
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 (
) = 30.33(
0.16) + 0.71(
0.13)log (
)
for masses
in the 0.5-3.0
range. This slope also agrees with what we find for the ONC stars in the mass range 0.5-3.0
,
namely 0.82
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
detection limits for low-mass stars in the
Tr16 (3 Myr), CygOB2 (2 Myr), and ONC (1 Myr) observations, as log (
)
30.5, 30.3 and 28.5
,
respectively, above which a source is detectable anywhere in the FOV. A mass-dependent
completeness fraction (
)
for our survey of Trumpler 16 stars was then computed by adopting the X-ray luminosity functions given by Preibisch & Feigelson (2005):
is
5% for 0.1-0.5
,
40% for 0.5-0.9
,
and
55% in the
0.9-1.2
range. Our survey of Tr 16 is statistically complete for masses
1.5
.
In Fig. 10 we present
cumulative distributions for Trumpler 16, Cyg OB2, and ONC stars in the mass range 1.5-2.5
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
99.9% (
= -4.15), that young ONC (1 Myr) stars, with masses in the range 1.5-2.5
,
are intrinsically more luminous in X-rays than their Tr16 (3 Myr) counterparts.
![]() |
Figure 10:
Cumulative distributions of X-ray luminosities for Trumpler 16 low-mass stars with masses 1.5 ![]() ![]() ![]() ![]() |
Open with DEXTER |
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 (
,
,
and
computed as
)
. We computed
by means of an absorbed ( WABS) thermal plasma model ( APEC). Metal
abundance was fixed at Z = 0.3
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
/
relation. In spite of the observed scatter, one sees
that the median
/
for binaries is about 7 times larger (
8.3
10-7) than that of suspected single stars (
1.1
10-7)
. Four of the
13 known binaries in the analysis are well-separated systems: the LBV
Carinae, a probable long period binary with
2026
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
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.
We report here results of a deep Chandra X-ray observation
pointed toward the 3 Myr old star-forming region Trumpler 16. Source
detection was performed using the PWDdetect code, identifying 1035 X-ray sources in the 17'
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.
![]() |
Figure 11:
Log (![]() ![]() ![]() ![]() ![]() ![]() |
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 ( vs.
)
and color-color (
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
2.0 mag, while low
mass likely members seem to be slightly more absorbed,
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
,
and is very likely
complete down to
1.5
.
From the
vs. J-H diagram we estimate that
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. 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
< 1. We assumed an absorbed single-component
thermal emission model. The median
of the sources is
21.73 (cm-2). This value agrees well with the median
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
= 2.5
1030-6.3
1033 erg s-1. Their X-ray and bolometric luminosities roughly agree with the relation
/
,
albeit with an order of magnitude
dispersion. Low-mass stars have
ranging between 1030and 1031
(median
= 2.8
1030). Variable low
mass stars are on average 0.5 dex brighter (log (
)
31.0
).
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
,
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-
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.
Table 1: Trumpler 16source catalog.
Table 2: Near-IR counterparts of Trumpler 16 X-ray sources.
Table 3: X-ray source spectroscopy.