Issue |
A&A
Volume 519, September 2010
|
|
---|---|---|
Article Number | A34 | |
Number of page(s) | 27 | |
Section | Catalogs and data | |
DOI | https://doi.org/10.1051/0004-6361/200911873 | |
Published online | 09 September 2010 |
Results from DROXO
III. Observation, source list, and X-ray properties of sources detected in the ``Deep Rho Ophiuchi XMM-Newton Observation''
I. Pillitteri1,2,9 - S. Sciortino2 - E. Flaccomio2 - B. Stelzer2 - G. Micela2 - F. Damiani2 - L. Testi3 - T. Montmerle4 - N. Grosso5,6 - F. Favata7 - G. Giardino8
1 - DSFA, Università degli Studi di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy
2 - INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy
3 - ESO - Karl-Scharzschild-Strasse 2, 85748 Garching bei München, Germany
4 - Laboratoire d'Astrophysique de Grenoble Université Joseph-Fourier, Grenoble, France
5 - Université de Strasbourg, Observatoire Astronomique de Strasbourg, 11 rue de l'université, 67000 Strasbourg, France
6 - CNRS, UMR 7550, 11 rue de l'université, 67000 Strasbourg, France
7 - ESA - Planning and Community Coordination Office, Science Programme, Paris, France
8 - Astrophysics Division - RSSD ESA, ESTEC, Noordwijk, The Netherlands
9 - SAO-Harvard Center for Astrophysics, Cambridge MA, USA
Received 18 February 2009 / Accepted 26 April 2010
Abstract
Context. X-rays from very young stars are powerful probes to
investigate the mechanisms at work in the very first stages of the star
formation and the origin of X-ray emission in very young stars.
Aims. We present results from a 500 ks long observation of
the Rho Ophiuchi cloud with a XMM-Newton large program named DROXO,
aiming at studying the X-ray emission of deeply embedded young stellar
objects (YSOs).
Methods. The data acquired during the DROXO program were reduced
with SAS software, and filtered in time and energy to improve the
signal to noise of detected sources; light curves and spectra were
obtained.
Results. We detected 111 sources, 61 of them associated with Ophiuchi
YSOs as identified from infrared observations with ISOCAM.
Specifically, we detected 9 out of 11 Class I objects, 31 out
of 48 Class II and 15 out 16 Class III objects. Six objects
out of 21 classified Class III candidates are also detected. At
the same time we suggest that 15 Class III candidates that
remain undetected at
are not related to the cloud population. The global detection rate is
64%. We have achieved a flux sensitivity of
erg s-1 cm-2. The
to
ratio shows saturation at a value of
-3.5 for stars with
K or 0.7
as observed in the Orion Nebula. The plasma temperatures and the
spectrum absorption show a decline with YSO class, with Class I
YSOs being hotter and more absorbed than Class II and III
YSOs. In one star (GY 266) with infrared counterpart in 2MASS and
Spitzer catalogs we have detected a soft excess in the X-ray spectrum,
which is best fitted by a cold thermal component less absorbed than the
main thermal component of the plasma. This soft component hints at
plasma heated by shocks due to jets outside the dense circumstellar
material.
Key words: stars: coronae - stars: formation - X-rays: stars - Galaxy: formation
1 Introduction
Low-mass stars in the pre-main sequence (PMS) phase are characterized by intense X-ray emission. X-ray activity has not been yet firmly confirmed in the very initial protostellar cores that accrete material from a surrounding envelope (Class 0 objects), because X-rays are either absent or completely absorbed by the dense circumstellar material. In Class I, II, and III objects and until the Zero Age Main Sequence stage the X-ray emission is very strong when compared to solar-age stars and also characterized by impulsive variability (Feigelson & Montmerle 1999). In the case of Class III objects all X-rays are thought to originate from a magnetized stellar corona, which likely resembles a scaled-up version of the corona of late-type Main Sequence stars and the Sun. Additional X-ray production mechanisms are possibly at work in Class I and II objects. At these young evolutionary stages the magnetic field drives the material in-falling from the disk to the stellar surface where the matter becomes shocked. Furthermore, outflows along the disk axis may interact with the circumstellar material and form shocks. The plasma in the shocks in both in- and outflows may become heated up to a few million degrees, thus emitting X-rays (Favata et al. 2002; Sacco et al. 2008; Giardino et al. 2007; Kastner et al. 2005; Pravdo et al. 2001).
Significant effort has been devoted to the X-ray observation of star-forming regions,
which provide natural laboratories where to study the X-ray emission from young stars and its
implications for the mechanisms of star formation.
The Ophiuchi cloud is among the nearest star-forming regions (120 pc, Loinard et al. 2008)
and has been extensively studied from the infrared (IR) to the X-rays bands.
The
Ophiuchi cloud is shaped as a dense multicore structure,
hosting more than 200 members comprising young stellar objects
(YSOs) in all evolutionary stages from Class 0 to Class III (Wilking et al. 2008).
While a dispersed population of optically visible stars associated with the Upper Scorpius OB association with an age of
5 Myr
is present in the region, the studies in the IR band have revealed a
younger population of Pre Main Sequence (PMS) stars and protostellar
objects with ages of only 0.3-1 Myr (Luhman & Rieke 1999),
which are thus younger than YSOs in other star-forming regions like the Taurus Molecular Cloud (
1-5 Myr) and the Orion Nebula (
1 Myr, Hillenbrand 1997).
The mid-IR survey with ISOCAM on board the ISO satellite reported by
Bontemps et al. (2001, hereafter Bo01)
allowed the detection of a population of 16 Class I, 123
Class II and 77 Class III YSOs, adding 71 previously
unknown members of Classes I and II. X-ray studies carried
out with the Einstein and ROSAT satellites had revealed several tens of embedded Class I, II,
and III stars that are highly active in X-rays and are characterized by a very strong time
variability (Casanova et al. 1995; Montmerle et al. 1983).
Casanova et al. (1995) tentatively identified several X-ray sources with Class I protostars with
ROSAT/PSPC observations. However, X-ray observations with the
ROSAT/High-Resolution Imager were needed to confirm X-ray emission
from Class I protostars (Grosso 2001; Grosso et al. 1997).
Recent observations with the XMM-Newton (Ozawa et al. 2005) and Chandra satellites (Imanishi et al. 2001; Flaccomio et al. 2003)
have increased the number of candidate PMS members of the region.
These studies suggested that accreting YSOs have increasing X-ray
activity going from Class I to Class III and decreasing
plasma temperatures. The aim of the Deep Rho Ophiuchi XMM-Newton Observation (DROXO) was to obtain a high-sensitivity survey in the Core F of
Ophiuchi
Cloud by means of a long, almost continuous observation. We report the
data analysis and the X-ray properties of the PMS stars observed
during the DROXO program. The structure of the paper is as follows:
Sects. 2 and 3 describe the observation and the data
analysis, Sects. 4 and 5 describe the sensitivity of the
survey and the nature of the sources detected in DROXO. In Sect. 6
we discuss the X-ray emission of the sample of classified YSOs, in
particular we explore the behavior of plasma temperatures, absorption,
X-ray luminosity and ratio
with respect to
effective temperatures, masses and evolutionary stages.
In Sect. 7 we give a summary. Appendices A and B contain
tables with the list of detected sources, the best-fit parameters of
spectra, the upper limits to rate for undetected YSOs in the ISOCAM
sample, the Spitzer counterparts to DROXO sources, and an example of
one page of the atlas (published electronically only) with the EPIC
spectrum and light curve of each source.
2 The observation
The Deep Rho Ophiuchi XMM Observation is a large program carried with the XMM-Newton
satellite pointed toward core F in the
Ophiuchi Cloud for almost eight consecutive days. The nominal pointing was at RA = 16
27
19.5
and Dec = -24
41
40.9
(J2000), and the net exposure time was
515 ks.
The pointing approximately coincides with that of the 33 ks XMM-Newton
exposure studied by Ozawa et al. (2005). The observation has been carried out in five subsequent satellite orbits (0961
0965), keeping the same position angle in all orbits.
During the first orbit chip no. 6 of the MOS 1 camera was damaged apparently by a micrometeorite impact
and has stopped functioning, so that only a
28.5 ks exposure is
available by that chip in DROXO. No damage was registered on the other chips of MOS 1. During orbit no 0964 the
Ophiuchi exposure was split in two segments, separated by
25 ks, due to a triggered TOO observation.
During the first three orbits the drift of the satellite was larger
than 6
,
and this influenced the computation of exposure maps and
the subsequent analysis of summed data as explained below.
![]() |
Figure 1: Merged EPIC image of events recorded in the time-filtered data to enhance the signal-to-noise of faint sources. Colors encode the following bands: 0.3-1.0 keV (red), 1.0-2.5 keV (green), 2.5-8.0 keV (blue), respectively. The MOS 1, 2, and PN images are normalized by effective area and exposure time to reduce instrumental artifacts like CCD gaps. |
Open with DEXTER |
In Fig. 1 we show the EPIC MOS1, MOS2 and PN images added in a pseudo-color Red-Green-Blue image. The bands chosen for red, green and blue components are: 0.3-1.0 keV, 1.0-2.5 keV, 2.5-8.0 keV. Sources with a softer/less absorbed spectrum are redder than sources with a hard/highly absorbed spectrum. Each CCD image has been divided by the proper exposure map and by the average effective area in the energy band to normalize the different efficiencies of the three cameras.
3 Data analysis
The observation data files (ODF) were processed with the SAS software
(version 6.5) to produce full field-of-view event lists calibrated
in both energy and position. We subsequently filtered these event files
and retained only the events with energy
in the 0.3-10 keV band and those that triggered simultaneously at
most
two nearby pixels. This step was executed for each detector
(MOS 1, MOS 2, PN) and for each of the exposure segments of
the five orbits (hereafter step 1). The large satellite drift is
not taken into account automatically by the SAS software which results
in wrong exposure maps when they are created with the default values.
We needed
to reduce the ATTREBIN parameter of task EEXPMAP to 0.5
and to choose the most accurate algorithm (of which the default was the
fastest and less precise one) to obtain correctly computed exposure
maps.
3.1 Source detection
We performed source detection with the PWXDETECT code developed at INAF-Osservatorio Astronomico di Palermo (Damiani et al. 1997a,b). The code allows the detection of sources starting from unbinned photon positions recorded in several datasets from MOS 1, MOS 2 and PN cameras, through a multiscale analysis of mexican-hat wavelet convolved images. All data were properly scaled by time and effective area of each CCD detector, obtaining a flux image of the sum of all exposures. At the end of the process the count rates of detected sources were re-scaled to a reference instrument which, in our case, was the EPIC MOS 1. We refer to the count rates derived in this way as MOS1 equivalent count rates.
In order to improve the detection of faint sources, we filtered the data obtained at step 1, excluding those time intervals with a high background count rate. In fact, the background during DROXO was highly variable for a significant fraction of the exposure time (see Fig. 2). We excluded all events registered in intervals when the rate on the whole image, shown in Fig. 2, was higher than a given threshold. The rate threshold was chosen in a way that maximizes the signal-to-noise-ratio (SNR) for faint sources (cf. Damiani et al. 2003; Sciortino et al. 2001) and improves the source detection process toward faint sources as expected. The net cleaned exposure times after this filtering are 198.1 ks (38%) 273.2 ks (53%) and 213.2 ks (41%) for MOS 1, 2 and PN, respectively.
![]() |
Figure 2: Top panel: light curve of all events recorded with MOS 1. Labels for the individual satellite orbits are given on top of the plot. The gray shaded area represents the time interval when MOS 1 was turned off after the micro-meteorite impact. The horizontal line is the threshold count rate that maximizes the SNR of the full image as described in the text. Bottom panel: signal-to-noise function vs. cumulative time. |
Open with DEXTER |
We detected 111 point sources with a significance threshold corresponding to two spurious detections in the whole field-of-view. The threshold for the detection of faint sources was determined from the analys is of a large set of simulations of background-only images and then running the detection code on these images. The simulations provide a value for the significance threshold to retain at most 1-2 spurious sources per field in real data.
Source positions, off-axis distance, exposure times, and
X-ray count rates are listed in Table A.1 for all detected
sources. We also report the 2 MASS designation and other names
from the literature for the optical counterparts identified in
Sect. 5.1. The last column
of Table A.1 contains a flag pointing to the source of the
previous identifications. We built an atlas of spectra and light curves
for each source which is available online; in Appendix B we show a page of this atlas as an example.
![]() |
Figure 3: Two examples of choices for the source and background extraction regions. The left panel shows two nearby but distinct sources. The right panel shows an extreme case in which the region for the faint source is shrunk to minimize the influence of the bright source. The background regions are chosen as near as possible to the sources to which they refer and at the same distance of the readout node for PN. |
Open with DEXTER |
3.2 X-ray spectra
We produced both light curves and spectra of all sources by selecting the photons in circular regions around the source positions. The radii of the regions depend on the source intensity, on possible of nearby sources and the geometry of the CCD. We used regions on the same CCD for source and background; for the background, we avoided to include pixels and columns poorly calibrated in energy. Figure 3 shows two examples of the choice of source and background extraction regions. For PN spectra a further constraint was to have both source and background extraction regions at nearly the same distance from the CCD readout node.
To produce the spectra we filtered the photons with respect to FLAG
(chosen to be equal to zero) and PATTERN
(less or equal to four) as indicated in the SAS User's Guide. Given
that the choice of
GTIs was tailored toward the faintest sources, at this stage we further
improved the choice of GTIs for bright sources. The time-filtering we
adopted for spectra starts from the Good Time Intervals
(GTIs) defined initially to perform source detection and iteratively adds temporal bins from the light-curve that contribute to increase the total SNR of the spectrum. The procedure starts by adding the time bin with the largest
individual SNR and is then iterated considering time-bins of
decreasing SNR until no gain in total SNR is obtained.
As a result, the net exposure time of each spectrum is different.
3.2.1 Background subtraction
We noticed that scaling the background by the geometric areas of the source and background extraction regions leads to incorrect estimates of the background, especially for faint sources and/or during times of high background. The effect is that background-corrected light curves of faint sources are either directly or inversely correlated with the background light curves and often lead to negative net count-rates. We understood this as the effect of a vignetted background component that is not properly taken into account by a purely geometric scaling factor. In order to correctly estimate the background contribution to the photons extracted in the source region we constructed background maps for each instrument (MOS 1, 2 and PN) and for each orbit. These were built by removing large regions around detected sources from the images, and by subsequently smoothing and interpolating the maps over the source extraction areas. Although the resulting scaling factors differ from the purely geometric ones by only a few percent, the difference is relevant in cases when the background dominates the count rate in the source regions and mitigates the above mentioned spurious effects on the light curves.
3.2.2 Model fitting of spectra
We analyzed the spectra of all X-ray sources in the 0.3-10 keV band. For a given EPIC detector (MOS 1, MOS 2 or PN) spectra from all orbits were summed up; analogously, background spectra were obtained; the response matrices and ancillary response files for each spectrum were multiplied and then summed by weighting by the exposure time. The background was scaled according to the procedure described in Sect. 3.2.1.
The spectra were grouped prior to fitting with XSPEC
v. 12.3 to obtain at least a minimum SNR in each bin, by
considering both source and background photons. In order to obtain the
largest number of meaningful grouped spectra we adopted two schemes of
grouping procedure based on high and low SNR of the final spectrum. For
this purpose we used the same procedure as in the ACIS_EXTRACT package for the analysis of Chandra ACIS data, adapted to our EPIC spectra to take into account the background.
We grouped the spectra using two thresholds for the minimum SNR to be
obtained in each bin, i.e. we generated two sets of spectra, one set
with a high and one with a low SNR per bin. The minimum SNR in each
spectral bin was imposed on the basis of the source count statistics,
defining our low SNR binning and
defining our high SNR binning. Where possible, we tried to obtain binned spectra with at least eight bins.
The spectra of all detectors were fitted simultaneously. In some cases we had to discard the data of one or two of the three EPIC cameras. These cases occur when a source is on CCDs gaps of one or two cameras, thus leading to a wrong estimate of the point spread function fraction contained in the extraction area.
The spectra were fitted
with one-temperature (1-T), two-temperature (2-T) and,
in some cases, three temperature (3-T) APEC models (Smith et al. 2001)
plus absorption (WABS) (Morrison & McCammon 1983),
the free parameters were the absorption column ,
the temperatures, and the
emission measures. The abundance pattern was fixed
to that found in PMS coronae from the Orion Nebula in the COUP survey (Maggio et al. 2007). Only in four cases described below (Elias 29, SR 12A, IRS 42, GY 266) a more complex model was required.
For Elias 29 (src. 38) and SR 12A (src 53) the global abundance
scaling was left as a free parameter to achieve a better fit.
For Elias 29 we confirm an unusually high abundance (
)
with respect to those found in PMS coronae (
),
as already reported by Favata et al. (2005). The source IRS 42/GY 252 (src. 54) is located in the wings of the much brighter X-ray source
corresponding to SR 12A (see right panel of Fig. 3).
To take account of the contribution by SR 12A we added the best fitting 3-T APEC model multiplied
by a constant factor representing the amount of contamination to the model of IRS 42/GY252.
This scaling factor was a free fitting parameter.
The spectrum of IRS 42/GY 252 itself can be described with a 1-T APEC model. The spectrum and the best-fit model are shown in Fig. 4.
As can be seen, IRS 42/GY 252 is highly absorbed and the soft
emission is attributed entirely to the contamination by SR 12A.
For the source no. 61 (GY 266) we used a model given by the
sum of two differently absorbed APEC models to account for the soft excess visible below 1.0 keV, as discussed in Sect. 6.4 and Fig. 11.
![]() |
Figure 4: Spectrum and best-fit model of IRS42/GY252. The stepped curves are the contributions of SR 12A (light line) and that of IRS42 (dark line). The spectrum of IRS42/GY252 is modeled by an absorbed single temperature APEC model. |
Open with DEXTER |
Generally, the best fit was chosen on basis of the probability
of obtaining a
higher than the observed one
.
Our threshold for acceptable fits was
.
Whenever this criterion was satisfied by the spectra with high SNR per bin
(typically for medium-high statistics spectra) we chose these;
otherwise we selected the fits obtained from the spectra with low SNR per bin.
We always selected the best-fit model with the smallest number of free parameters.
In seven cases, even for some bright sources, no formally acceptable fit was found and we
allowed for a lower
.
However, we accepted the best fit results
by visually checking that the overall shape of the spectrum is well reproduced by the model.
With this procedure we obtained spectral fits for 91/111 sources.
Table A.2
summarizes the results. The columns report the source number, the
available EPIC datasets, the data sets we chose for fitting, a flag
indicating the binning of the spectra we used (high or low SNR), the
type of model, the model best-fit parameters, the unabsorbed flux and
luminosity in the 0.3-10 keV band, the
statistics, degrees of freedom and the probability
.
For those weak sources without a spectral analysis, we calculated
fluxes and X-ray luminosities in the 0.3-10.0 keV band by using
their count rates and PIMMS software
assuming a 1-T model with kT and
equal to the median of values derived from the best-fit procedure of spectra (see Sect. 4) and a distance of 120 pc.
3.3 Photospheric parameters
In Sects. 5 and 6 we will focus our discussion on the X-ray properties of the sample of PMS stars classified by Bo01.
The photospheric stellar parameters (
,
and masses) for the Class II and III objects from the list of
Bo01 were estimated from the near IR (2MASS) photometry. The procedure
we used closely follows that adopted by Bo01 and improved by Natta et al. (2006). We assumed that the J-band
emission from these sources is dominated by the stellar photosphere and
is only marginally contaminated by the emission from circumstellar
material and also that the IR colors of Class II sources can
be described as the emission from a passive circumstellar disk as
described by Meyer et al. (1997).
These assumptions obviously do not apply to Class I sources and
for this reason photospheric parameters for these objects were not
estimated.
We used the Cardelli et al. (1989) extinction law with Rv=4.4, which we think is
appropriate for Ophiuchus. A small number of sources (15%)
have colors slightly bluer than those of reddened main sequence stars,
presumably due to photometric uncertainties, which are on the order of
0.1 mag, while the offsets of these objects with respect to the
reddened sequence range between 0 and 0.15 mag.
Dereddening these sources by extrapolating the colors
of Class II and III sources would produce an overestimate of
the extinction. For these objects we assigned the colors of the closest
photosphere model on the reddened main sequence. Table A.5 lists the effective temperatures, masses and bolometric luminosities of ISOCAM objects.
The values of the J-band extinction we derive are very similar to ones from Natta et al. (2006), with only the significant exception of WL 16, for which our procedure produces a significantly higher extinction.
4 Sensitivity of the survey
![]() |
Figure 5: Sensitivity map of the DROXO field in units of ct s-1 per point source (logarithmic scale). The map refers to the sum of PN and MOS instruments. The intensity scale is indicated in the vertical strip on the right side of the map. Crosses indicate the positions of detected sources. |
Open with DEXTER |
![]() |
Figure 6:
Left panel: MOS 1 rate from Ozawa et al. (2005) vs. DROXO MOS 1 rate. The lines mark the region for rate variations of factors 2 and 5. Right panel: the same for the ACIS rate from Flaccomio et al. (2003).
Upper limits to count rate of sources undetected in the compared
samples are indicated by vertical and horizontal arrows. Errors are
quoted at 1 |
Open with DEXTER |
The faintest detected source in DROXO has a net count rate of
ct ks-1 in MOS equivalent units.
In Fig. 5
we plot a sensitivity map in units of count rate per point source of
the DROXO field of view. It is built starting from a smoothed
background map, after removing the contributions of the sources, and
taking into account the vignetted exposure map and the threshold used
for detection. The map shows that the sensitivity varies by a
factor 2.5 in the area covered
by the three EPIC cameras, and it is quite constant across a 60%
fraction of the field of view because of the characteristics of the
detector. To translate the limit count rate in a limit flux we needed a
conversion factor that
depends on the spectrum temperature and, critically, on the absorption.
From the analysis of the spectra we know that the absorption is also
strongly variable from source to source by more than a factor 100,
which can be due to local material around the sources or dense cloud
material in front of the objects, or both. The plasma temperatures that
we find from the spectra vary from
keV to a few keV. The combined action of local absorption and plasma temperatures of undetected sources
causes the limit flux to vary more than the range observed in the sensitivity map showed in Fig. 5.
For example, we can estimate a conversion factor from count rates to fluxes by using the median of the column absorption (
cm-2) and of the plasma temperatures (kT = 3.1 keV) from the spectral fits, and this gives us
erg cts-1 cm-2. This yields a limiting unabsorbed flux in the 0.3-10 keV band of
erg s-1 cm-2 and a luminosity of
erg s-1. The cfs derived from all sources are in the range
erg cm-2 cts-1, yielding a limit flux comprised in the range
erg cm-2 s-1 and luminosities
erg s-1,
respectively. These values are indicative of the sensitivity we
achieved across the field of view, but local strong absorption can
significantly lower the actual sensitivity in that position.
4.1 Comparison with other X-ray surveys of
Ophiuchi
In Fig. 6 we show the scatter plots of the MOS 1 count rate from Ozawa et al. (2005) and Chandra-ACIS rate (which has an effective area comparable to the EPIC MOS) from Flaccomio et al. (2003)
vs. the DROXO MOS equivalent rate. In DROXO we detected 80 of the X-ray sources reported by Ozawa et al. (2005), while 7 of the sources found by Ozawa and collaborators remain undetected in DROXO; analogously we identified 62 sources
from Flaccomio et al. (2003),
while 30 of those sources have not been detected. We computed an upper
limit to the count rates for those sources that are detected by Flaccomio et al. (2003) and Ozawa et al. (2005) but are undetected in DROXO (indicated as arrows in Fig. 6).
We also re-analyzed the previous XMM-Newton observation reported by Ozawa et al. (2005)
with the same detection procedure as in DROXO, and we calculated the
upper limits to count rates for DROXO sources undetected in the Ozawa
survey. Analogously we calculated upper limits to count
rates for two sources detected in DROXO but undetected by Flaccomio et al. (2003). The count rates of sources detected in DROXO and in the two other surveys globally agree within a factor 5.
Because there is no systematic trend in the scatter we conclude that it
is likely to be attributed to variability. Although a detailed study of
the time variability is beyond the scope of this paper,
we report here that 52% of the sources have variable rates at more than
90% significance level when compared with the Ozawa et al. survey,
and 79% are found variable to be compared with the survey of Flaccomio
et al.
In the 33 ks XMM-Newton survey of Ozawa et al. (2005) the faintest detected source has a count rate of 0.5 ct ks-1. Using the conversion factor derived by us for DROXO the corresponding luminosity is
erg s-1,
which scales reasonably well with the exposure times of both surveys.
Flaccomio et al.'s survey reaches a limiting rate similar to
DROXO, likely due to the lower Chandra-ACIS background with respect to
that found in our EPIC observation. Most of the ACIS sources undetected
in DROXO are near the sensitivity limit of DROXO, and a variability of
factor 2 can easily explain their missed detections.
Sources 33 (GY 304) and 43 are brighter in DROXO than in Ozawa et al. survey by more than a factor five. DoAr 25 suffered of strong pile-up in the Ozawa et al. survey. These authors do not report a direct measurement of the count rate for this object, but they derived a luminosity from the spectral analysis. From its X-ray luminosity we inferred a count rate greater than 200 ct/ks. We checked that in DROXO DoAr 25 and the other two brightest sources in the field, SR 12A and IRS 55, are not affected by pile-up. The variability of DoAr 25 with respect to the XMM observation detailed in Ozawa et al. (2005) is at least a factor 4.
Among the most variable sources, WL 2 (labeled in Fig. 6) shows a big flare in the DROXO light curve (Fig. 7). The quiescent rate of WL 2 is 6 ct ks-1 and the peak rate is
300 ct ks-1 explaining its large offset position in the scatter plot of the right panel in Fig. 6. The quiescent rate is consistent with the one measured by Ozawa et al. (2005). The long duration of this flare on WL 2 (
35 ks) underlines the need for long observations of YSOs to properly assess their quiescent emission.
![]() |
Figure 7:
PN Light curve of WL 2 (source nr. 15).
The bold light curve is the net rate, while the thin line is the total rate (source + background). A large flare lasting |
Open with DEXTER |
5 The nature of the X-ray sources
5.1 Optical and IR counterparts of the DROXO sources
We searched for optical and IR counterparts of DROXO sources in the 2MASS, Spitzer C2D (Evans et al. 2003a) and
ISOCAM (Bo01) surveys as well as in the optical and IR surveys of Natta et al. (2006), Wilking et al. (2005),
Barsony et al. (2005) and Luhman & Rieke (1999).
The match radius takes into account the uncertainties on the X-ray
source position in DROXO, which is generally on the order of
, and other catalogs. The
mean value for the match radius is
with a 0.1-0.9 quantile range of
.
As anticipated in Sect. 3.1, Table A.1
reports the literature names of DROXO counterparts and a flag
indicating in which catalogs the source is identified. Sources with
counterparts in the ISOCAM survey are indicated with letter ``I''
in the ``Flag'' column; sources identified only in Spitzer are
indicated with letter ``S''; sources which have only a X-ray
identification in the surveys by Flaccomio et al. (2006) and/or Ozawa et al. (2005)
are indicated with letter ``X''. Among the 111 X-ray
sources 78 have a counterpart in 2 MASS and 61 have a
counterpart in ISOCAM data. Two objects have no 2 MASS and no
ISOCAM counterpart, but have been identified in the literature (IRS50
and WSB46).
An additional 16 X-ray sources have been newly identified with Spitzer. This leaves 16 DROXO sources without known optical/IR counterpart. Of these, 6 were detected in previous X-ray surveys and 10 are presented here for the first time. We calculated upper limits to the X-ray fluxes and luminosities for the sample of undetected ISOCAM objects falling in the field of view of DROXO using the conversion factor derived in Sect. 4.
5.2 Unidentified X-ray sources
The inspection of light curves and spectra can give some
indication on the nature of sources with unknown optical/IR
counterpart, designated ``U'' or ``X'' in Table A.1.
Five unidentified sources (Src. 5, 7, 16, 26, 58)
are too faint for spectral analysis, but the other five unidentified
X-ray sources (Src. # 12, 19, 88, 96, 110) have spectral
parameters compatible with the expectation for
a YSO in Oph, i.e. kT of a few keV and
.
Furthermore, Src. 110
shows impulsive time variability similar to flares typical for PMS
stars.
There are also six sources (Src. # 20, 29, 41, 45, 48, 95) without
optical/IR counterpart that have been detected in previous X-ray
surveys. For Src. # 20, 29 and 41 we had too few counts to obtain
meaningful spectral fits, and their light curves show some variability,
but no clearly identifiable flare. For Src. # 45, 48 and 95 the
spectral analysis gives a high absorption (above 1022 cm-2).
While the temperature of Src. # 95 is not constrained, Src. # 45 and 48
have plasma temperatures of 5.4 and 4.5 keV. For these two sources
a power law best fit to their spectra is also acceptable. These three
objects could have characteristics consistent with those of highly
embedded YSOs, although they are undetected on the millimetric surveys
of Jørgensen et al. (2008); Johnstone et al. (2000); Motte et al. (1998).
The lack of 2MASS counterparts suggests that these sources could have
an extragalactic nature. On the other hand, it is not ruled out tthat
they are very low mass PMS stars or even brown dwarfs.
5.3 Mid-IR photometry from Spitzer
The Spitzer photometry for the X-ray sources identified in the C2D
survey are listed in Table A.4. We derived magnitudes in the 3.6,
4.5, 5.8, 8.0 m IRAC bands and constructed
the color-color diagram shown in Fig. 8
(color index
[3.6]-[4.5] vs. color index
[5.8]-[8.0]).
This diagram provides a rough classification of YSOs according to Allen et al. (2004) and Hartmann et al. (2005).
Normal, unreddened stars or Class III/Weak T Tauri stars
with very low reddening should occupy a region centered
around (0,0). Reddening due to matter along the line of sight
tends to disperse vertically data points along the [3.6]-[4.5] IRAC
color index, while the [5.8]-[8.0] color index is not affected by this
source of reddening (see Flaherty et al. 2007).
Infrared emission from a circumstellar disk in Class I and II
sources moves the objects both vertically and horizontally toward
redder values of both indexes. We expect to find Class II YSOs in
the region marked with the black box or above it and very embedded
protostars (Class 0/I) are expected to lie at
and at
[5.8]-[8.0] > 1 (Allen et al. 2004).
![]() |
Figure 8:
IRAC color-color diagrams of objects in the DROXO field with information on the reddening from J or V band photometry or from |
Open with DEXTER |
We plotted only the objects for which we have information on the reddening from J and V band photometry or from .
In the left panel of Fig. 8
the symbols represents the IRAC colors before correction for reddening.
Objects detected in DROXO are marked with squares, diamonds are
Class I YSOs as
classified by Bo01, triangles are Class II YSOs, asterisks are
Class III YSOs. The points with ISOCAM classification follow the
distribution outlined above, with Class III YSOs around or
dispersed vertically above the (0,0) point, and Class II YSOs
located in the rectangle or dispersed vertically above this region.
In order to investigate the nature of these latter
``dispersed'' objects we dereddened the colors of those objects for
which we had information on the extinction. If available
we used the reddening
or
from the literature, otherwise the column absorption
from the X-ray spectra with the conversion factor of
(Mathis 1990).
To convert these extinctions to the Spitzer bands we followed
the calibrations by Flaherty et al. (2007) and Carpenter (2001).
The relation between reddening in 2MASS
and IRAC bands provided by Flaherty et al. (2007) for
Ophiuchi indicates that no reddening is present in the
[5.8]-[8.0] color index.
Therefore, dereddening shifts the objects vertically downwards in Fig. 8. The right hand panel of Fig. 8 shows the dereddened IRAC colors. Evidently, most objects with
[3.6]-[4.5] > 0.7in the left hand side of Fig. 8, fall in the ``canonical'' IRAC
Class II area after dereddening. Therefore, the majority of objects
with apparently protostellar IRAC colors are probably strongly reddened
Class II sources. We note also that the ISOCAM sub-sample without X-ray counterparts has the same color distribution
as the ISOCAM/DROXO sample. This indicates that no bias is present against X-ray properties.
Prisinzano et al. (2008) found that
Class I protostars in the Orion Nebula can be separated into two
distinct subclasses: the first, Ia, with a rising SED from K to 8 m and lower X-ray emission level than the second one, Ib, characterized by rising SED up to 4.5
m.
This latter group shows IRAC colors more similar to those of
Class II objects. The first subclass, Class Ia, is instead
populated by more embedded protostars with lower and, perhaps, more
absorbed X-ray emission. They find 23 Class 0-Ia YSOs, 22
Class Ib YSOs
and 148 Class II YSOs. In COUP the fraction of Class Ia on
Class II YSOs is
16%, the fraction of Class I versus Class II YSOs in Orion is
30%. In DROXO the fraction of Class I objects with [3.6]-[4.5] > 1 (after dereddening) is
7% with respect to Class II objects.
This sample should be composed by very embedded objects similar to
Class Ia defined by Prisinzano et al. (2008). The number of very embedded protostars in Core F of
Ophiuchi is thus very low with respect to Class II objects. The same evidence has been obtained by Jørgensen et al. (2008). In that paper the fraction of Class I to Class II is reported to be lower in
Ophiuchi than that present in other star-forming regions like Perseus (10% in
Ophiuchi
vs. 90% in Perseus, respectively), hinting for a different age or
formation time-scales for these two star-forming regions.
Three objects (GY 92, GY 289, GY 203) which have been
classified as Class III by Bo01 are found with Spitzer colors
similar to Class II YSOs. Furthermore, we find three objects
classified as Class II YSOs (GY 146, GY 240,
GY 450) that are bluer than 0.4 mag in ([5.8]-[8.0]) IRAC
color. These objects may have transition disks (see Kim et al. 2009, and references therein). The star with the reddest IRAC colors located in the top right corner of Fig. 8
is WL 22/GY 174. It is classified as a Class II object
and probably it suffers of strong foreground extinction as pointed out
by Wilking et al. (2001). This is supported also from our spectral fit to its X-ray spectrum: we find a
absorption value of
cm-2, a factor of 5 higher than the average of the DROXO sample.
6 Results
6.1 Coronal temperatures and absorption
Table A.2 lists the temperatures obtained through a thermal model
fit to the X-ray spectra as described in Sect. 3.2.2.
Most of the spectra are reasonably well described by a model with a
single thermal component; usually we find temperatures higher than
1.5 keV and absorption column
higher than 1022 cm-2.
In a few cases two or three thermal components are needed to fit the
spectra, depending on the characteristics of the spectrum and its count
statistics. For these cases we calculated the average of the two or
three plasma temperatures weighted by their emission measures to obtain
a representative mean temperature for these spectra. The median of the
representative plasma temperatures is
3.1 keV with a 1
range between 1.6 and 8 keV (the 10%-90% quantile range is
0.8-13 keV). The 10%-90% quantile range of the
column is
cm-2 and the median is
cm-2.
![]() |
Figure 9:
Top left panel: X-ray luminosities (0.3-10 keV) versus
effective temperatures of YSOs identified with ISOCAM
objects (Bo01). Symbols are: squares Class II YSOs, triangles
Class III YSOs. Upper limits to luminosities were calculated for
ISOCAM objects not detected in DROXO (vertical arrows). Top right panel:
|
Open with DEXTER |
6.2 X-ray luminosities and stellar parameters
We examine the relation between X-ray luminosity and X-ray to
bolometric luminosity ratio and mass and effective temperature for all
X-ray sources for which stellar parameters could be determined in
Sect. 3.3 (Fig. 9).
The X-ray luminosities for the sample with known stellar parameters increase from
erg s-1 up to
erg s-1 in the 3000-5000 K range (Fig. 9, top left panel) or
0.1-0.7
(left bottom panel). Most of the YSOs have effective temperatures
around 3000-5000 K (K type stars), with the exception of the
three hot stars WL 5, WL 16 and WL 19.
The fraction of X-ray detections among ISOCAM YSOs increases from
to
.
The boxes in the bottom left panel of Fig. 9 are the 10-90% quantiles obtained with ASURV software (Schmitt 1985; Lavalley et al. 1992) for
values in four ranges of mass (<0.2, 0.2-0.5,
). For these mass ranges the median values of
are
28.3, 29.0, 30.0 and 30.0, respectively.
The COUP sample contains very few upper limits in each mass range whereas in DROXO the fraction of upper limits is
30% reducing the
medians in each mass range. As discussed below in Sect. 6.3,
we suggest that a fraction of Class III YSOS are likely spurious
cloud members. We fitted the relation between X-ray luminosity and mass
excluding these suspect members in the Class III sample obtaining
with the same procedure used by Preibisch et al. (2005)
for COUP sample. We also fitted X-ray luminosity and mass relation for
Class II and Class III samples separately, finding similar
slopes but different normalizations. The relations are

and

for Class II and Class III YSOs respectively. The slopes of these relations are very similar to those found for ONC in COUP program (1.44, Preibisch et al. 2005). The lower normalizations that we find in the





The
ratio increases with decreasing stellar masses (Fig. 9, bottom and top right panels) saturating at
-3.5 for stars cooler than 5000 K and with masses of
0.7
,
slightly below the ``canonical'' value of -3 observed in MS young stars, but very similar to the value reported by Preibisch et al. (2005) for Orion YSOs (-3.6).
From the relation between mass and X-ray luminosity and the hypothesis that the ratio
for PMS stars is saturated at a level of -3.5,
Telleschi et al. (2007) derived an empirical mass - bolometric luminosity relation for PMS stars
which is shallower than the relation that yields for MS stars. By considering that the saturation limit of
for our sample of YSOs of
Ophiuchi is similar to that of TMC and ONC (
-3.5, see Fig. 9) and that the slope of the relation between mass and
is quite similar for
Ophiuchi, ONC and TMC, it is suggested that the relation between mass and
found from Telleschi et al. (2007)
for coeval PMS stars should apply also for our sample.
Among the three hottest stars, characterized by effective temperatures
around 10 000 K, we detected the Class III YSOs
WL 19 and WL 5, while WL 16 (Class II) remains
undetected. X-ray spectra of WL 19 and WL 5 show quite hot
plasma temperatures (kT = 3.7 keV and 4.5 keV, respectively) and high absorption (
and
cm-2, respectively). The undetected Class II YSO WL 16 is a peculiar object: it consists of a massive star (
,
)
that illuminates a circumstellar disk visible only at mid-IR wavelengths (Ressler & Barsony 2003). It suffers of strong absorption (
)
likely due to a foreground screen of cloud material.
Taking into account this high absorption, the upper limit to luminosity
is inversely correlated with the plasma temperature: for kT = 4 keV we obtain
,
for kT = 1 keV
,
for kT = 0.5 keV
.
According to their position in the HR diagram these three objects are
intermediate-mass pre-MS stars and their X-ray luminosities and
ratios are in the range typically observed for Herbig Ae/Be stars (Stelzer et al. 2009).
The coolest object is the binary system WL 2/GY 128 (Barsony et al. 2005).
A discrepancy between its spectral type and the effective temperature
is present in the literature. While its spectral type is comprised
between K and M as reported by Luhman & Rieke (1999), its temperature, estimated by Natta et al. (2006), is very low (2300 K). As discussed in Sect. 4.1, WL 2 has undergone a huge flare during DROXO. Its quiescent X-ray luminosity is
1029 erg s-1
which is typical for young K-M type stars but unexpectedly high
for a low mass brown dwarf. Likely the photospheric temperature of this
object is more similar to that of late K or M-type stars.
Only two bona fide brown dwarfs are in the field of view. We detected GY 310 (log
erg s-1), but not GY 141, which was detected by Ozawa et al. (2005) during a flare with a flux 90 times higher than in a previous Chandra observation.
6.3 X-ray emission of different YSO classes
We examined the X-ray detection rates for YSOs in different evolutionary states referring to the YSO classification of Bo01.
Their catalog comprises 16 Class I, 123 Class II, and 38 bona-fide
Class III. The latter are classified as Class III YSOs
on the basis of absence of IR excess and detection in X-ray images (ROSAT) and/or radio
band (VLA). Given the low ROSAT sensitivity at



![[*]](/icons/foot_motif.png)


Previous studies of star-forming regions have shown that for each
given mass range, Class I and II YSOs show lower X-ray
luminosities than Class III YSOs (Telleschi et al. 2007; Neuhäuser et al. 1995; Preibisch et al. 2005; Flaccomio et al. 2003; Prisinzano et al. 2008).
Figure 10
(top left panel) shows the Kaplan-Meier estimators of X-ray luminosity
function (XLF) of Class I, II, and III objects following the
Bo01 classification. Upper limits are mostly concentrated below the
lowest detection values for Class II and III YSOs. We plot also
the XLF of Class III YSOs (dots and dashed line) after excluding
15 upper limits of Class III YSO candidates from Bo01. In
this way we try to correct the bias introduced by spurious members
that very likely contaminate the Class III YSO sample as discussed
above.
This corrected Class III XLF is similar to the XLF of Class I YSOs.
The corrected Class III XLF shows higher luminosity levels with respect to Class II XLF. Two sample tests yield a probability of 99% for the two distributions of being different,
thus suggesting that in
Ophiuchi
Class III YSOs are more luminous than Class II YSOs. We
observe also that Class I YSOs have similar X-ray luminosity
levels compared with Class III objects.
The high X-ray detection rate and high X-ray luminosities among
Class I objects are surprising. For comparison, in ONC Prisinzano
et al. (2008) have found that the X-ray luminosities of
Class I are lower than those of more evolved Class II YSOs.
Our X-ray bright Class I YSOs could be explained as an effect of
different mass distributions in different samples but, given that we
cannot estimate masses of Class I objects with our method, we
cannot test this hypothesis.
Table 1: X-raydetection rates and median of logarithm of X-ray luminosities and plasma temperatures for the different YSO groups following the classification given by Bontemps et al. (2001).
![]() |
Figure 10:
Top left panel: Kaplan-Meier estimators of cumulative distributions of |
Open with DEXTER |
The distribution of mean plasma temperatures in Fig. 10
(top right panel) shows that on average the emitting plasma in
Class I YSOs is hotter than in Class II and Class III
YSOs. The median plasma temperatures for Class I, II and III
are reported in Table 1. Two-sample tests give a probability of 99% to reject the null hypothesis that the distributions
of temperature for Class III and Class I YSOs are drawn from the same distribution, the probability is 98% when comparing Class II and Class I temperatures, and
90% when comparing Class II and Class III YSOs temperature distributions. As expected, the distributions of
show that the absorption is higher on average in Class I YSOs than in Class II and III (see Fig. 10 bottom left and right panels), indicating circumstellar matter in Class I objects.
6.4 Source with a soft X-ray excess
We discuss the case of the source nr. 61, which shows
an excess of soft photons in the 0.3-1.0 keV band significantly
different
from a typical coronal singly-absorbed multi-temperature plasma. This
object has no counterpart in the ISOCAM survey, although in the
literature is identified with GY 266 and is classified as a
variable star by Alves de Oliveira & Casali (2008). Furthermore, it has an IR counterpart in 2MASS and Spitzer catalogs. From IRAC and MIPS photometry this object is cataloged as a normal star (Evans et al. 2003b).
The DROXO spectrum of this object is shown in Fig. 11.
We modeled the spectrum with two thermal components, which were differently absorbed. The soft component has
kT = 0.470.230.67 keV (5.5 MK) and is
absorbed by a
column of
cm-2. This is one of the softest thermal components we have determined in the whole sample of DROXO spectra. The hot component has
kT = 1.81.42.4 keV absorbed by
cm-2. The emission measure ratio of soft to hot component is
1%. The hot component has the typical temperature and column density of a YSO in
Oph, whereas the weakly absorbed component is unusually soft.
Güdel et al. (2007) have reported similar examples of X-ray spectra
with soft excess from the XEST survey. They modeled them with a 2-T plasma
with individual absorptions and hypothesized that the soft component be
associated to the emission from shocks in jets. Bonito et al. (2007)
have studied the soft X-ray emission that could arise from bipolar jets
in YSOs through extensive MHD simulations. They find that a jet less
dense than the ambient medium in which it propagates can emit soft
X-rays with plasma temperatures of 2-3 MK, very similar to that of
the soft component we found. Moreover, Bally et al. (2003) and Favata et al. (2006)
observed that the X-ray emission of HH 154 jet is located at the
base of the shock and that morphological changes are detected on a time
scale of four years. The X-ray luminosity of that jet is
erg s-1. In our case the unabsorbed
luminosity of the soft component in the spectrum of GY 266 is much lower than that of HH 154 (
erg s-1).
Given the angular resolution of EPIC we cannot spatially resolve the
soft component from the hot component. Thanks to angular resolution of
Chandra (Güdel et al. 2008) resolved separated soft emission from the coronal emission in the star DG Tau belonging to the Taurus Molecular Cloud.
The coincidence of the elongated soft X-ray emission with the direction of
the optical jet represents a clear case for the jet scenario in DG Tau.
In our case the lack of IR excess is difficult to reconcile with
the circumstellar disk, which is expected along with the jet.
Explaining the origin of the X-ray soft component of GY 266 requires further investigation of this star.
![]() |
Figure 11: X-ray spectrum of source # 61 (GY 266) with a soft component in the 0.3-1 keV energy range, less absorbed than the hot component. |
Open with DEXTER |
7 Discussion and conclusions
We described the Deep Rho Ophiuchi XMM-Newton Observation (DROXO),
aimed at exploring at high sensitivity the X-ray emission of YSOs in the
Ophiuchi core F region. We detected 111 X-rays sources, and for 91 of
them we obtained a model fit to their X-ray spectra. By using optical
and IR data we estimated the photospheric parameters of most of the
sample of X-ray detected YSOs.
We find 10 unidentified sources and a further 6 sources already detected in previous X-ray surveys, but without optical/IR counterpart. Three of them (src. num. 45, 48 and 95) show light curves with some impulsive variability. The spectra of src. num. 45 and 48 have a good fit with an absorbed power law, while src. num. 95 has a spectrum compatible with a thermal model. The lack of 2MASS counterparts could be a hint that they are of extragalactic nature, but, given the sensitivity of 2MASS catalog, it is not possible to rule out they are very low mass PMS stars or even brown dwarfs, especially src. num. 95.
The sensitivity of the survey is
erg s-1 cm-2 in flux and
(erg s
in
luminosity, but it strongly depends on the local absorption, which is
largely variable in the field-of-view. The sample of 96 classified
YSOs given by Bo01 in the DROXO field of view has allowed us to explore
the X-ray emission from Class I, II, and III YSOs of
Ophiuchi. The X-ray properties of
Ophiuchi
PMS stars obtained from DROXO were compared with those of the Orion
Nebula Cloud PMS stars studied in the COUP survey.
When fitting the relation between mass and X-ray luminosity with a
power law, the index that we find for DROXO sample is quite similar to
that found for ONC and TMC (Telleschi et al. 2007; Preibisch et al. 2005).
With respect to COUP we find a lower normalization, suggesting that on
average the PMS stars in the range of masses between 0.1-2
in our sample are less luminous than their analogs in COUP. Young stellar objects in
Ophiuchi exhibit a saturation of the ratio
/
near the ``canonical'' value of
10-3 for masses between 0.5 and 1
,
while stars with masses below 0.5
show a lower limit of saturation
/
.
This is consistent with what was found by Preibisch et al. (2005) for the COUP sample.
We detected a large fraction of YSOs in the field of view, 94% of which
are bona fide Class III stars, 77% of Class II and 82% of
Class I YSOs, respectively. We confirm the high detection rate
among Class I YSOs found by Imanishi et al. (2001). The detection rate in our Class I sample is higher than the analog rate found in the COUP survey by Prisinzano et al. (2008) (82% vs. 62%) despite the higher sensitivity of
COUP with respect to DROXO.
Prisinzano et al. (2008) make a distinction between
Class Ia (characterized by a rising SED from K to 8 m) and Class Ib (rising SED up to 4.5
m).
The detection rates in these two subclasses are 44% and 82%
respectively. X-ray luminosities of Class Ib objects are higher
than those of Class Ia objects.
The rate of X-ray detection of Class Ib objects in Orion is
identical to that of our sample of Class I objects in
Ophiuchi.
With support of Spitzer photometry we suggest that only a small
fraction of our Class I stars are deeply embedded objects. From
X-ray data the high detection rate suggests that our Class I
sample could be mainly formed by Class Ib YSOs as defined by
Prisinzano et al. (2008).
These objects in Orion have X-ray luminosities more similar to
Class II objects. In
Oph we find that our nine detected Class I stars are on average more luminous than Class II stars.
Bontemps et al. (2001) list
39 stars defined as Class III candidates, 21 of which are
surveyed in DROXO. Six out of 21 are detected in DROXO, providing
support for their PMS nature. The other 15 are undetected below
.
Likely they are not PMS stars
and thus not members of the
Oph
cloud. By means of XLFs we evaluated the level of emission of
Class I, II, and III YSOs. After excluding the suspect
Class III YSOs, we find that Class I YSOs emit at the same
level of Class III YSOs (median
and 29.7, respectively) and both samples are more luminous than
Class II YSOs (median
).
In COUP Class III objects are more luminous than Class I and II
for masses between 0.5 and 1.2
,
while for masses in 0.1-0.5
Class Ib sources are slightly more luminous than Class II and III (see Fig. 8 in Prisinzano et al. 2008).
The analysis of X-ray spectra leads us to conclude that the mean plasma temperatures of Class I protostars are higher than in Class II and Class III YSOs. We find a clear trend of decreasing plasma temperatures passing from Class I to Class III objects, differently from what is observed in ONC where Prisinzano et al. (2008) do not find a significant evolution of temperatures from Class I to Class III objects. Four Class III stars and one Class II star (GY 112, SR 12A, GY 296, GY 380 and GY 3) have spectra noticeably softer than the rest of YSOs with mean temperatures around 0.7-0.8 keV.
The absorption derived from X-ray spectra is higher in Class I
than in Class II and III. Given the high absorption among
Class I YSOs, it is impossible to determine if a soft thermal
component in their spectra is completely absorbed or is not present at
all. This soft component would reduce the intrinsic kT
for Class I objects, making them more similar to the ones in the
ONC. However, higher absorption should lead to an underestimate of the
luminosities of heavily absorbed spectra like those of Class I
objects, while we observe
a luminosity for them higher than that of in Class II and
Class III YSOs.
Furthermore, Prisinzano et al. (2008) have found in extensive
simulations that varying
does not significantly influence the XLF.
The star GY 266 shows a peculiar X-ray spectrum composed by two thermal components that are differently absorbed. The hot and heavily absorbed component is similar to that found in other PMS stars, while the soft one, less absorbed, could arise from unresolved jet shocks.
AcknowledgementsThe authors acknowledge an anonymous referee for the useful comments that improved the paper. I.P., S.S., E.F., B.S., G.M., and F.D. acknowledge financial support from ASI/INAF contract nr. I/023/050. LT acknowledges support from ASI-INAF I/016/07/0.
Appendix A: Tables
Table A.1:
List of sources detected in the DROXO EPIC images.
The columns refer to X-ray positions, position accuracy,
summed exposure times from the three EPIC cameras, count rates scaled to MOS detector units (see Sect. 3.1) with 1
errors, 2MASS identifier, ISOCAM identifier, SED classification given by Bontemps et al. (2001),
literature names, flag column for identifications (see footnote).
Coordinate errors take into account position uncertainties calculated
by the detection code and the best-match radius from the optical/IR
objects. Class III candidates from Bo01 are indicated by ``III?''.
Table A.2: X-ray Parameters from model fit to source spectra. The errors are quoted at the 90% confidence region of best fit. Unabsorbed X-ray fluxes and luminosities are calculated in the 0.3-8 keV energy band.
Table A.3:
Upper limits to count rate, X-ray flux and luminosity in the 0.3-10 keV for the ISOCAM YSOs (Bontemps et al. 2001) undetected in the DROXO field of view. Count rates are scaled to MOS detector units as in Table A.1 (see Sect. 3.1).
Class III candidates from Bo01 are indicated by ``III?''. Fluxes
and luminosities are calculated by assuming a conversion factor from
count rates of
and a distance of 120 pc. For ISOCAM nr. 92 (WL 16) we have taken into account the large absorption toward this star (
mag, cf. Sect. 6.2).
Table A.4:
Spitzer IRAC fluxes of IR counterparts to
DROXO X-ray sources from the catalog by Evans et al. (2003b).
The
index is fitted from K to Spitzer MIPS1 band (
m). The classification given in the last column is that given in Evans et al. (2003b).
Table A.5: Effective temperatures and bolometric luminosities of ISOCAM objects in DROXO (see Sect. 3.3).
Appendix B: Atlas of DROXO sources
We show here an example of the atlas that we have produced for each DROXO source. It is reported the number of the source, the literature name as in Bontemps et al. (2001), the EPIC image, the J-band 2MASS image, the EPIC spectrum, and the light curve. The atlas is available only online at the following web address: http:www.astropa.unipa.it/ pilli/atlas_droxo_sources.pdf
![]() |
Figure B.1:
|
Open with DEXTER |
References
- Allen, L. E., Calvet, N., D'Alessio, P., et al. 2004, ApJS, 154, 363 [NASA ADS] [CrossRef] [Google Scholar]
- Alves de Oliveira, C., & Casali, M. 2008, A&A, 485, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bally, J., Feigelson, E., & Reipurth, B. 2003, ApJ, 584, 843 [NASA ADS] [CrossRef] [Google Scholar]
- Barsony, M., Ressler, M. E., & Marsh, K. A. 2005, ApJ, 630, 381 [NASA ADS] [CrossRef] [Google Scholar]
- Bonito, R., Orlando, S., Peres, G., Favata, F., & Rosner, R. 2007, A&A, 462, 645 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bontemps, S., André, P., Kaas, A. A., et al. 2001, A&A, 372, 173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Carpenter, J. M. 2001, AJ, 121, 2851 [NASA ADS] [CrossRef] [Google Scholar]
- Casanova, S., Montmerle, T., Feigelson, E. D., & Andre, P. 1995, ApJ, 439, 752 [Google Scholar]
- Damiani, F., Flaccomio, E., Micela, G., et al. 2003, ApJ, 588, 1009 [NASA ADS] [CrossRef] [Google Scholar]
- Damiani, F., Maggio, A., Micela, G., & Sciortino, S. 1997a, ApJ, 483, 350 [NASA ADS] [CrossRef] [Google Scholar]
- Damiani, F., Maggio, A., Micela, G., & Sciortino, S. 1997b, ApJ, 483, 370 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, N. J., Dunham, M. M., Jørgensen, J. K., et al. 2009, ApJS, 181, 321 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, II, N. J., Allen, L. E., Blake, G. A., et al. 2003a, PASP, 115, 965 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, II, N. J., Allen, L. E., Blake, G. A., et al. 2003b, PASP, 115, 965 [NASA ADS] [CrossRef] [Google Scholar]
- Favata, F., Fridlund, C. V. M., Micela, G., Sciortino, S., & Kaas, A. A. 2002, A&A, 386, 204 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Favata, F., Micela, G., Silva, B., Sciortino, S., & Tsujimoto, M. 2005, A&A, 433, 1047 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Favata, F., Bonito, R., Micela, G., et al. 2006, A&A, 450, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Feigelson, E. D., & Montmerle, T. 1999, ARA&A, 37, 363 [NASA ADS] [CrossRef] [Google Scholar]
- Flaccomio, E., Micela, G., & Sciortino, S. 2003, A&A, 402, 277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Flaccomio, E., Micela, G., & Sciortino, S. 2006, A&A, 455, 903 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Flaherty, K. M., Pipher, J. L., Megeath, S. T., et al. 2007, ApJ, 663, 1069 [CrossRef] [Google Scholar]
- Giardino, G., Favata, F., Pillitteri, I., et al. 2007, A&A, 475, 891 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Grosso, N. 2001, A&A, 370, L22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Grosso, N., Montmerle, T., Feigelson, E. D., et al. 1997, Nature, 387, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Güdel, M., Telleschi, A., Audard, M., et al. 2007, A&A, 468, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Güdel, M., Skinner, S. L., Audard, M., Briggs, K. R., & Cabrit, S. 2008, A&A, 478, 797 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hartmann, L., Megeath, S. T., Allen, L., et al. 2005, ApJ, 629, 881 [NASA ADS] [CrossRef] [Google Scholar]
- Hillenbrand, L. A. 1997, AJ, 113, 1733 [NASA ADS] [CrossRef] [Google Scholar]
- Imanishi, K., Tsujimoto, M., & Koyama, K. 2001, ApJ, 563, 361 [NASA ADS] [CrossRef] [Google Scholar]
- Johnstone, D., Wilson, C. D., Moriarty-Schieven, G., et al. 2000, ApJ, 545, 327 [NASA ADS] [CrossRef] [Google Scholar]
- Jørgensen, J. K., Johnstone, D., Kirk, H., et al. 2008, ApJ, 683, 822 [NASA ADS] [CrossRef] [Google Scholar]
- Kastner, J. H., Franz, G., Grosso, N., et al. 2005, ApJS, 160, 511 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, K. H., Watson, D. M., Manoj, P., et al. 2009, ApJ, 700, 1017 [NASA ADS] [CrossRef] [Google Scholar]
- Lavalley, M. P., Isobe, T., & Feigelson, E. D. 1992, BAAS, 839, 24 [Google Scholar]
- Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJ, 675, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Luhman, K. L., & Rieke, G. H. 1999, ApJ, 525, 440 [NASA ADS] [CrossRef] [Google Scholar]
- Maggio, A., Flaccomio, E., Favata, F., et al. 2007, ApJ, 660, 1462 [NASA ADS] [CrossRef] [Google Scholar]
- Mathis, J. S. 1990, ARA&A, 28, 37 [NASA ADS] [CrossRef] [Google Scholar]
- Mayne, N. J., & Naylor, T. 2008, MNRAS, 386, 261 [NASA ADS] [CrossRef] [Google Scholar]
- Meyer, M. R., Calvet, N., & Hillenbrand, L. A. 1997, AJ, 114, 288 [NASA ADS] [CrossRef] [Google Scholar]
- Montmerle, T., Koch-Miramond, L., Falgarone, E., & Grindlay, J. E. 1983, ApJ, 269, 182 [NASA ADS] [CrossRef] [Google Scholar]
- Morrison, R., & McCammon, D. 1983, ApJ, 270, 119 [NASA ADS] [CrossRef] [Google Scholar]
- Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 [NASA ADS] [Google Scholar]
- Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Neuhäuser, R., Sterzik, M. F., Schmitt, J. H. M. M., Wichmann, R., & Krautter, J. 1995, A&A, 297, 391 [NASA ADS] [Google Scholar]
- Ozawa, H., Grosso, N., & Montmerle, T. 2005, A&A, 438, 661 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Padgett, D. L., Rebull, L. M., Stapelfeldt, K. R., et al. 2008, ApJ, 672, 1013 [NASA ADS] [CrossRef] [Google Scholar]
- Pravdo, S. H., Feigelson, E. D., Garmire, G., et al. 2001, Nature, 413, 708 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Preibisch, T., Kim, Y.-C., Favata, F., et al. 2005, ApJS, 160, 401 [NASA ADS] [CrossRef] [Google Scholar]
- Prisinzano, L., Micela, G., Flaccomio, E., et al. 2008, ApJ, 677, 401 [NASA ADS] [CrossRef] [Google Scholar]
- Ressler, M. E., & Barsony, M. 2003, ApJ, 584, 832 [NASA ADS] [CrossRef] [Google Scholar]
- Sacco, G. G., Argiroffi, C., Orlando, S., et al. 2008, A&A, 491, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schmitt, J. H. M. M. 1985, ApJ, 293, 178 [Google Scholar]
- Sciortino, S., Micela, G., Damiani, F., et al. 2001, A&A, 365, L259 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, ApJ, 556, L91 [NASA ADS] [CrossRef] [Google Scholar]
- Stelzer, B., Robrade, J., Schmitt, J. H. M. M., & Bouvier, J. 2009, A&A, 493, 1109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Telleschi, A., Güdel, M., Briggs, K. R., Audard, M., & Palla, F. 2007, A&A, 468, 425 [Google Scholar]
- Wilking, B. A., Bontemps, S., Schuler, R. E., Greene, T. P., & André, P. 2001, ApJ, 551, 357 [NASA ADS] [CrossRef] [Google Scholar]
- Wilking, B. A., Meyer, M. R., Robinson, J. G., & Greene, T. P. 2005, AJ, 130, 1733 [NASA ADS] [CrossRef] [Google Scholar]
- Wilking, B. A., Gagné, M., & Allen, L. E. 2008, Star Formation in the Ophiuchi Molecular Cloud, ed. B. Reipurth, 351 [Google Scholar]
Footnotes
- ... impact
- http://xmm.vilspa.esa.es/external/xmm_news/items/MOS1-CCD6/index.shtml
- ... software
- See http://xmm.vilspa.esa.es/sas
- ... online
- See http:www.astropa.unipa.it/ pilli/atlas_droxo_sources.pdf
- ... light-curve
- The binning of the light curve for each source is chosen in
an iterative way, starting from
ks and progressively enlarging
until the bin with maximum net counts has at least 50 counts or
becomes larger than 10 Ks.
- ...
software
- See http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html
- ... data
- See http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html for the details of the algorithm.
- ...(Evans et al. 2003a)
- See also http://ssc.spitzer.caltech.edu/legacy/c2dhistory.html
- ... objects
- Four of six are classified as ``YSO'' and two of six are classified as ``star'' in the Spitzer catalog, Evans et al. (2003b); Padgett et al. (2008); Evans et al. (2009) for this scheme of classification.
All Tables
Table 1: X-raydetection rates and median of logarithm of X-ray luminosities and plasma temperatures for the different YSO groups following the classification given by Bontemps et al. (2001).
Table A.1:
List of sources detected in the DROXO EPIC images.
The columns refer to X-ray positions, position accuracy,
summed exposure times from the three EPIC cameras, count rates scaled to MOS detector units (see Sect. 3.1) with 1
errors, 2MASS identifier, ISOCAM identifier, SED classification given by Bontemps et al. (2001),
literature names, flag column for identifications (see footnote).
Coordinate errors take into account position uncertainties calculated
by the detection code and the best-match radius from the optical/IR
objects. Class III candidates from Bo01 are indicated by ``III?''.
Table A.2: X-ray Parameters from model fit to source spectra. The errors are quoted at the 90% confidence region of best fit. Unabsorbed X-ray fluxes and luminosities are calculated in the 0.3-8 keV energy band.
Table A.3:
Upper limits to count rate, X-ray flux and luminosity in the 0.3-10 keV for the ISOCAM YSOs (Bontemps et al. 2001) undetected in the DROXO field of view. Count rates are scaled to MOS detector units as in Table A.1 (see Sect. 3.1).
Class III candidates from Bo01 are indicated by ``III?''. Fluxes
and luminosities are calculated by assuming a conversion factor from
count rates of
and a distance of 120 pc. For ISOCAM nr. 92 (WL 16) we have taken into account the large absorption toward this star (
mag, cf. Sect. 6.2).
Table A.4:
Spitzer IRAC fluxes of IR counterparts to
DROXO X-ray sources from the catalog by Evans et al. (2003b).
The
index is fitted from K to Spitzer MIPS1 band (
m). The classification given in the last column is that given in Evans et al. (2003b).
Table A.5: Effective temperatures and bolometric luminosities of ISOCAM objects in DROXO (see Sect. 3.3).
All Figures
![]() |
Figure 1: Merged EPIC image of events recorded in the time-filtered data to enhance the signal-to-noise of faint sources. Colors encode the following bands: 0.3-1.0 keV (red), 1.0-2.5 keV (green), 2.5-8.0 keV (blue), respectively. The MOS 1, 2, and PN images are normalized by effective area and exposure time to reduce instrumental artifacts like CCD gaps. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Top panel: light curve of all events recorded with MOS 1. Labels for the individual satellite orbits are given on top of the plot. The gray shaded area represents the time interval when MOS 1 was turned off after the micro-meteorite impact. The horizontal line is the threshold count rate that maximizes the SNR of the full image as described in the text. Bottom panel: signal-to-noise function vs. cumulative time. |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Two examples of choices for the source and background extraction regions. The left panel shows two nearby but distinct sources. The right panel shows an extreme case in which the region for the faint source is shrunk to minimize the influence of the bright source. The background regions are chosen as near as possible to the sources to which they refer and at the same distance of the readout node for PN. |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Spectrum and best-fit model of IRS42/GY252. The stepped curves are the contributions of SR 12A (light line) and that of IRS42 (dark line). The spectrum of IRS42/GY252 is modeled by an absorbed single temperature APEC model. |
Open with DEXTER | |
In the text |
![]() |
Figure 5: Sensitivity map of the DROXO field in units of ct s-1 per point source (logarithmic scale). The map refers to the sum of PN and MOS instruments. The intensity scale is indicated in the vertical strip on the right side of the map. Crosses indicate the positions of detected sources. |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Left panel: MOS 1 rate from Ozawa et al. (2005) vs. DROXO MOS 1 rate. The lines mark the region for rate variations of factors 2 and 5. Right panel: the same for the ACIS rate from Flaccomio et al. (2003).
Upper limits to count rate of sources undetected in the compared
samples are indicated by vertical and horizontal arrows. Errors are
quoted at 1 |
Open with DEXTER | |
In the text |
![]() |
Figure 7:
PN Light curve of WL 2 (source nr. 15).
The bold light curve is the net rate, while the thin line is the total rate (source + background). A large flare lasting |
Open with DEXTER | |
In the text |
![]() |
Figure 8:
IRAC color-color diagrams of objects in the DROXO field with information on the reddening from J or V band photometry or from |
Open with DEXTER | |
In the text |
![]() |
Figure 9:
Top left panel: X-ray luminosities (0.3-10 keV) versus
effective temperatures of YSOs identified with ISOCAM
objects (Bo01). Symbols are: squares Class II YSOs, triangles
Class III YSOs. Upper limits to luminosities were calculated for
ISOCAM objects not detected in DROXO (vertical arrows). Top right panel:
|
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Top left panel: Kaplan-Meier estimators of cumulative distributions of |
Open with DEXTER | |
In the text |
![]() |
Figure 11: X-ray spectrum of source # 61 (GY 266) with a soft component in the 0.3-1 keV energy range, less absorbed than the hot component. |
Open with DEXTER | |
In the text |
![]() |
Figure B.1:
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.