A&A 468, 391-403 (2007)
DOI: 10.1051/0004-6361:20065559
The XMM-Newton extended survey of the Taurus molecular cloud
N. Grosso1 - K. R. Briggs2 - M. Güdel2 - S. Guieu1 - E. Franciosini3 - F. Palla4 - C. Dougados1 - J.-L. Monin1,5 - F. Ménard1 - J. Bouvier1 - M. Audard6 - A. Telleschi2
1 - Laboratoire d'Astrophysique de Grenoble,
Université Joseph-Fourier,
38041 Grenoble Cedex 9, France
2 - Paul Scherrer Institut, 5232 Villigen und Würenlingen,
Switzerland
3 - INAF - Osservatorio Astronomico di Palermo, Piazza del
Parlamento 1, 90134 Palermo, Italy
4 - INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico
Fermi 5, 50125 Firenze, Italy
5 - Institut Universitaire de France
6 - Columbia Astrophysics Laboratory, Columbia University,
550 West 120th Street, New York, NY 10027, USA
Received 5 May 2006 / Accepted 5 August 2006
Abstract
Aims. We report the X-ray properties of young (3 Myr) bona fide brown dwarfs of the Taurus Molecular Cloud (TMC).
Methods. The XMM-Newton Extended Survey of the TMC (XEST) is a large program designed to systematically investigate the X-ray properties of young stellar/substellar objects in the TMC. In particular, the area surveyed by 15 XMM-Newton pointings (of which three are archival observations), supplemented with one archival Chandra observation, allows us to study 17 brown dwarfs with M spectral types.
Results. Half of this sample (9 out of 17 brown dwarfs) is detected; 7 brown dwarfs are detected here for the first time in X-rays. We observed a flare from one brown dwarf. We confirm several previous findings on brown dwarf X-ray activity: a log-log relation between X-ray and bolometric luminosity for stars (with
)
and brown dwarfs detected in X-rays, which is consistent with a mean X-ray fractional luminosity
;
for the XEST brown dwarfs, the median of
(including upper limits) is -4.0; a shallow log-log relation between X-ray fractional luminosity and mass; a log-log relation between X-ray fractional luminosity and effective temperature; a log-log relation between X-ray surface flux and effective temperature. We find no significant log-log correlation between the X-ray fractional luminosity and
.
Accreting and nonaccreting brown dwarfs have a similar X-ray fractional luminosity. The median X-ray fractional luminosity of nonaccreting brown dwarfs is about 4 times lower than the mean saturation value for rapidly rotating low-mass field stars. Our TMC brown dwarfs have higher X-ray fractional luminosity than brown dwarfs in the Chandra Orion Ultradeep Project.
Conclusions. The X-ray fractional luminosity declines from low-mass stars to M-type brown dwarfs, and as a sample, the brown dwarfs are less efficient X-ray emitters than low-mass stars. We thus conclude that while the brown dwarf atmospheres observed here are mostly warm enough to sustain coronal activity, a trend is seen that may indicate its gradual decline due to the drop in photospheric ionization degree.
Key words: stars: low-mass, brown dwarfs - X-rays: stars - ISM: individual objects: Taurus molecular cloud
Pre-main-sequence low-mass stars, i.e. T Tauri stars, show a high
level of X-ray emission, which is generally attributed to an active
corona, an enhanced version of the magnetic corona on the contemporary
Sun (see for a review Feigelson & Montmerle 1999). The dynamo mechanism
producing the magnetic field in these fully convective stars is still
discussed (Preibisch et al. 2005a). X-ray observations of brown
dwarfs (BDs) allow investigation of the magnetic activity in the
substellar regime, where masses are lower than about 0.075 .
The first X-ray detection of a BD was reported by Neuhäuser & Comerón (1998), who identified in the Cha I dark cloud an X-ray source detected by ROSAT - in early 1991, before the near-infrared observations of the first BDs were reported (Rebolo et al. 1995; Nakajima et al. 1995) - as the counterpart of a young BD. ROSAT detected only a few BDs or very low-mass stars, all of them young and located in star-forming regions (Mokler & Stelzer 2002; Neuhäuser et al. 1999; Comerón et al. 2000).
The new generation of X-ray satellites, XMM-Newton and Chandra, allows now
detection of more BDs thanks to their increase of sensitivity. The
number of young or intermediate-age BDs detected in X-rays
(e.g., Stelzer et al. 2004; Briggs & Pye 2004; Preibisch & Zinnecker 2001; Tsuboi et al. 2003; Imanishi et al. 2001; Stelzer et al. 2006; Preibisch & Zinnecker 2002; Ozawa et al. 2005)
exceeds largely the one of (older) field BDs, which have been detected
so far mainly during X-ray flares
(Rutledge et al. 2000; Stelzer 2004). Recently, the X-ray properties of 33
young BDs (Slesnick et al. 2004) in the Orion Nebula Cluster (ONC) located
at 450 pc were studied by Preibisch et al. (2005b) in the Chandra Orion Ultradeep Project (COUP).
We investigate in this paper the X-ray properties of a sample of 17 BDs in the nearby (140 pc) Taurus Molecular Cloud (TMC), surveyed in the XMM-Newton Extended Survey of the TMC (XEST; Güdel et al. 2007). First, we define in Sect. 2 the TMC BD sample, review previous X-ray results on it, and report the BDs detected in the XEST. Analysis of the X-ray variability of the TMC BDs is presented in Sect. 3, and X-ray spectral properties are determined in Sect. 4. In Sect. 6, we compare the X-ray luminosities of the BDs and low-mass stars observed in the XEST. Comparison of the X-ray fractional luminosities of BDs in the XEST and COUP is made in Sect. 7. We discuss the origin of the BD X-ray emission in the broader context of the X-ray emission of cool stars in Sect. 8. Finally, we summarize our results in Sect. 9.
For a typical TMC member, with an age of 3 Myr, the
stellar/substellar boundary at 0.075
is between M6V and
M6.5V according to pre-main-sequence tracks from Baraffe et al. (1998). Therefore, we
compiled from the literature objects with spectral type later than
M6V to build our bona fide BD sample
.
The first optical imaging and spectroscopic surveys of BDs in the TMC
concentrated on high stellar density regions, while the majority of
the volume occupied by the molecular clouds was left unexplored
(Luhman 2004; Briceño et al. 2002; Luhman 2000; Briceño et al. 1998). Recently, Guieu et al. (2006)
performed a large scale optical imaging survey of TMC with the
Canada-France-Hawaii telescope with the CFH12k and MEGACAM large-scale
optical cameras, covering a total area of 28 deg2, and
encompassing the densest parts of the cloud as well as their
surroundings. By employing all-sky catalogs,
Luhman (2006) considered a
area, large
enough to encompass all of the TMC. These optical surveys of the TMC,
conducted since nearly 10 years, yielded the identification of
42 BDs.
![]() |
Figure 1:
The XMM-Newton Extended Survey of the Taurus Molecular
Cloud (XEST). Contours show the 12CO emission (Dame et al. 1987) of
the Taurus Molecular Cloud (TMC), overlaid on a visual
extinction map (linear colour scale, with black colour
corresponding to
![]() ![]() |
Open with DEXTER |
Table 1:
The sample of the 17 young BDs surveyed in the
XEST. Column (2) gives the reference of the discovery paper. Col. (3)
give the 2MASS counterparts. Column (7) gives the references for the
spectral type (Col. 4), the optical extinction (when not
available derived from
using Rieke & Lebofsky 1985, Col. 5), and the
reference of the equivalent width of H
(
,
given in Col. 6) if it is different. Negligible optical extinction in
Col. (5) are indicated by 0. The effective temperature in Col. (8) has
been computed from the spectral type using the temperature scale
,
which is valid for young
stars with M spectral type (Guieu et al. 2006). The visual extinction in
Col. (5) is taken when available from the literature; for 2MASS J0421
and 2MASS J0422, we averaged the visual extinction computed from
J-H,
,
and the spectral type, using a dwarf sequence
(compiled from the literature; e.g., Leggett et al. 1998).
The bolometric luminosity of the substellar photosphere in
Col. (9) has been computed from I, J-band magnitudes and
(see Guieu et al. 2006); for 2MASS J0414 and MHO 4 (without
I-band magnitude available) the reference for the luminosity is
Luhman (2004) and Briceño et al. (2002), respectively. Column (10)
indicates accreting sources based mainly on
(see Sect. 6.2). The last column indicates detection in X-rays
(this work; see Table 2). References:
B98 = Briceño et al. (1998); B02 = Briceño et al. (2002); G06 = Guieu et al. (2006);
L04 = Luhman (2004); L06 = Luhman (2006); M01 = Martín et al. (2001);
M05 = Muzerolle et al. (2005a).
Two BDs of this sample were previously detected in X-rays by the ROSAT/Position Sensitive Proportional Counter (PSPC) - which was sensitive only to soft X-rays from 0.1 to 2.4 keV - but without being identified as substellar objects.
The BD MHO 4 (Luhman 2000; Briceño et al. 1998) was first
detected during a 4.0 ks ROSAT/PSPC
pointed observation of V826 Tau in August 1992 with an X-ray count rate of
4 counts ks-1 and proposed as a new weak-line T Tauri star of
Taurus-Auriga (source RXJ0431.3+1800 of Wichmann et al. 1996). MHO 4
was also detected during a 7.7 ks ROSAT/PSPC pointed observation of the dark cloud
Lynds 1551 in February 1993 with an X-ray count rate of 3.1 counts ks-1 and
proposed as a new member of the TMC (source L1551X 15 of Carkner et al. 1996,
associated with a "faint, very red
star''). Neuhäuser et al. (1999) reported these two X-ray
detections of MHO 4, considered at that time as a BD
candidate.
CFHT-BD-Tau 4, one of the BDs of Martín et al. (2001), was
detected during a 10.6 ks ROSAT/PSPC pointed observation of the
Heiles dark cloud 2 NW in March 1993 with an X-ray count rate of
4.5 counts ks-1 corresponding to a luminosity of
1029 erg s-1 (source HCL 2 NW-7a of Briceño et al. 1999, with no
optical counterpart but a NIR counterpart). Mokler & Stelzer (2002)
reported this X-ray source as the counterpart of this BD.
The higher sensitivity of the new generation of X-ray satellites,
XMM-Newton (Jansen et al. 2001) and Chandra (Weisskopf et al. 2002), allows detection
of BDs which are from 10 to
100 times fainter in
X-rays.
The XMM-Newton Extended Survey of the TMC (XEST) is a large program designed to systematically investigate the high-energy properties of young stellar and substellar objects in this nearest star-forming region. It is composed of 19 dedicated XMM-Newton pointings plus 8 archival XMM-Newton observations. The total XMM-Newton exposure time is about 300 h shared among these 27 pointings (Güdel et al. 2007). We will focus here on the X-ray detection of BDs with the three X-ray telescopes of XMM-Newton - two European Photon Imaging Cameras (EPIC) equipped with MOS CCD arrays (Turner et al. 2001), and the third carrying a pn CCD camera (Strüder et al. 2001) - simultaneous U-band observations obtained with the XMM-Newton optical/UV monitor are discussed in Grosso et al. (2007).
Figure 1 shows the XEST area and the locations of the 42 TMC BDs. Fifteen XMM-Newtonpointings (of which three are archival observations) surveyed 16 BDs. We looked also for serendipitous Chandraobservations at the Chandra X-ray Center's archives (see online Appendix A, Table A.1). We supplemented the XEST with the ACIS-I observation of L1527 which allow us to survey and detect one more BD (CFHT-BD-Tau 4); and with the ACIS-I observation of L1551 which provides a second epoch observation of MHO 4. Therefore, we surveyed in X-rays 17 out of the 42 TMC BDs, which allows us to survey the X-ray emission of 40% of the known BDs in the TMC.
Table 1 lists the properties of these 17 BDs.
Figure 2 gives their location in
the HR diagram. The coolest object of the surveyed BD sample has
a spectral type M9.5V, corresponding to about 0.02 (Baraffe et al. 2002). Their luminosities range from
0.075
to 0.001
.
The optical extinctions of
this sample range from a negligible value to about 9 mag with a mean
and a median of 1.5 and 0.5 mag, respectively.
The average low extinction of this sample is mainly due to the initial
selection process for TMC candidate BDs based on optical photometry,
and therefore biased towards low extinction BDs.
![]() |
Figure 2:
HR diagram of the 17 BDs of the TMC surveyed in the XEST. The
references for the computation of the effective temperatures and the
luminosities are given in Table 1. The
pre-main-sequence tracks from
Baraffe et al. (1998) are shown for comparison. Continuous lines
show mass tracks from 0.1 down to 0.02 ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3: Light curve of the X-ray flare from CFHT-BD-Tau 1. The bin size is 1000 s. The grey thick line shows the fit of the non-zero count rates using a quiescent level plus an exponential rise and decay. |
Open with DEXTER |
The list of X-ray sources in the XEST (Güdel et al. 2005) was completed with X-ray sources of the two archival Chandra observations (see online Appendix A for details on the data reduction of Chandra data). We cross-correlated the BD positions with the positions of all X-ray sources. We detected 9 BDs in X-rays (see Table 2). The two BDs previously detected with ROSAT were redetected with XMM-Newton and Chandra. The detection of MHO 4 in L1551 was previously published by Favata et al. (2003) with XMM-Newton, and Bally et al. (2003) with Chandra. Seven BDs are detected here for the first time in X-rays. The detection rate in X-rays of these 17 BDs is therefore 53%.
We detected mainly the BDs with luminosities greater than
0.01
,
i.e. the BDs with spectral type
earlier or equal than M8V (Fig. 2).
![]() |
Figure 4: X-ray spectra of TMC BDs. From top to bottom and left to right: 2MASS J0422 (pn spectrum), MHO 4 ( XMM-Newton and Chandra spectra), CFHT-Tau 5 ( XMM-Newton spectra), CFHT-BD-Tau 1 ( XMM-Newton spectra), CFHT-BD-Tau 4 ( Chandra spectra). Black, red, and green code for XMM-Newton/EPIC pn, MOS1, and MOS2 spectra, respectively. The continuous lines show our best fits obtained with an absorbed single-temperature or two-temperature (MHO 4 spectra) plasma model (see Table 3). |
Open with DEXTER |
Obvious variability was observed only from CFHT-BD-Tau 1, a
nonaccreting M7V BD, which displayed an X-ray (coronal) flare during
30% of the observation period (Fig. 3). The shape
of this light curve, a fast rise followed by exponential decay, is
typical of X-ray flares from
young stellar objects (e.g., Imanishi et al. 2003; Favata et al. 2005). We fit the
non-zero count rates with a quiescent level plus an exponential rise and decay
(Fig. 3). We estimate from this light curve fit
that the quiescent level contributes to 54% of the total counts
(quiescent+flare) detected from CFHT-BD-Tau 1.
The quiescent emission level is
counts ks-1, the flare amplitude is
counts ks-1, i.e. about 8 times the quiescent
level. The time scales of the rise and decay phase are
min and
min, respectively. Both
values are in the lower range of time scales observed in
the systematic study of X-ray flares from low-mass young stellar objects
in the
Ophiuchi star-forming region with
Chandra (Imanishi et al. 2003). In this latter star-forming region,
Ozawa et al. (2005) observed with XMM-Newton from GY310, a young M8.5V BD
(Wilking et al. 1999), an X-ray flare with a similar amplitude (peaking
6 times above the quiescent level) and decay time scale, but with a
rise phase about 5 times longer.
Assuming that the count rate scales linearly with the X-ray luminosity
(e.g., Mitra-Kraev et al. 2005), we compute, from the average X-ray
luminosity of CFHT-BD-Tau 1 (see Table 3),
that the total energy released by the CFHT-BD-Tau 1 flare is about
erg.
Table 2: Detection list of young BDs in the XEST. Column (1) numbers correspond to labels in Fig. 1. Columns (2)-(4) give the BD, satellite and target names, respectively. The naming of XMM-Newton sources in Col. (5) follows the convention of Güdel et al. (2007), where the two and three digits code for the field and the source number in this field, respectively. Columns (6)-(9) give X-ray source positions, total positional uncertainties, and distance to the 2MASS position, respectively. References: [FGMSD03] = Favata et al. (2003), [BFR03] = Bally et al. (2003). There is only pn data for X-ray source #1 as it falls in the gap of the MOS CCDs in window mode in this observation. X-ray source #4b is affected by pn gap. For X-ray source #9 only MOS data are available.
The X-ray light curves of the other BDs show no large variability (Fig. B.1). Two BDs, observed at two differents epochs, show variations in X-ray luminosity lower than a factor of 2 (Fig. B.2; see below Tables 3 and 4). We looked for a possible modulation of the EPIC light curve of MHO 4 (see Fig. B.2) using a Lomb-Scargle periodogram (e.g., Flaccomio et al. 2005), but we found no statistically significant periodic signal. The BD count rates observed in the XEST are mainly representative of a "quiescent'' emission, however (micro)flaring activity could be present even though unresolved in the light curves as observed in active stars (Güdel et al. 2003).
In total, we have 11 X-ray observations of the 9 BDs detected in the XEST
(Table 2). We performed spectral fits for
the two Chandra sources, and for these XMM-Newton sources with more than 100
counts (Fig. 4). The computation of XMM-Newton source
spectra with associated
redistribution matrix files (RMFs) and auxiliary response files (ARFs)
is detailed in Güdel et al. (2007), for Chandra observations see online
Appendix A. Spectral fitting was performed with XSPEC
(version 11.3; Arnaud 1996) using one or two optically thin
thermal plasma models (MEKAL; Mewe et al. 1995) with 0.3 times
the solar elemental abundances, combined with an X-ray absorption
model (WABS; Morrison & McCammon 1983), i.e. the same model as used
in COUP to allow us to directly compare with the result from
Preibisch et al. (2005b). For XMM-Newton spectra, we fitted simultaneously
the EPIC pn, MOS1, and MOS2 spectra.
Table 3 gives our best fit parameters obtained
for 4 BDs in 3 XMM-Newton and 2 Chandra observations. The XMM-Newton and Chandra spectra
of the brightest source, MHO 4, are better described with a
two-temperature plasma. The two plasma temperatures found separately
from XMM-Newton and Chandra spectra are consistent with each other. For CFHT-Tau 5,
the hydrogen column density, ,
determined by one
temperature spectral fitting is 3 (with an acceptable range of 2-4)
times lower than the value
predicted from
mag (Guieu et al. 2006) using the
ratio of Vuong et al. (2003) combined with
the
ratio of Cardelli et al. (1989) for
,
which leads to
cm-2 mag-1. This may point out a non-canonical
dust-to-gas ratio on the line of sight towards this object.
Table 3:
Spectral properties of BDs obtained from spectral fitting. To
fit the spectra (Fig. 4), we used a WABS
absorption model (Morrison & McCammon 1983) combined with one or two
MEKAL optically thin thermal plasma model (Mewe et al. 1995) with 0.3 times
the solar elemental abundances. MEKAL plasma models were computed
rather than
interpolated from a pre-calculated table. We used
statistics with
standard weighting. Column (4) gives the net source counts collected by
the instrument given in Col. (3). Confidence ranges at the 68% level
(
;
corresponding to
for Gaussian
statistics) are given in parentheses. The value of reduced
and
,
the degrees of freedom, are indicated in Col. (10). The emission
measures in Cols. (8) and (9) and the X-ray luminosity in the
0.5-8 keV energy range corrected for absorption in Col. (11) were
computed assuming a distance of 140 pc for the TMC.
The X-ray fractional luminosity,
,
is
given in the last column. For CFHT-BD-Tau 1, the second line gives an
estimate of the quiescent X-ray luminosity derived for the light curve
fit (see Sect. 3).
Table 4:
Spectral properties of BDs obtained from quantile analysis. Column (3)
gives the instrument name (MOS stands for MOS1+MOS2). Column (4) and (5)
indicate the exposure and the
net source counts collected by this instrument (i.e. for MOS the
exposure average and the sum of net source counts). Column (6) gives the
energies below which the net counts are 25%, 50%, and 75% of the
source net counts in the 0.5-7.3 keV energy range. The resulting
position in the quantile diagram (Fig. 5) is given in
Col. (7). The observed optical extinctions (see Table 1) were used to
disentangle temperature double solutions in the quantile diagram (#2 on pn), and
to estimate the hydrogen column density (Col. 8) for sources with low
constraint in the quantile diagram (#2, #6 and #7 on pn, and #9
on MOS), using the relation
cm-2 mag-1 (Cardelli et al. 1989; Vuong et al. 2003);
these values are between brackets in Col. 8. Where negligible optical
extinction was measured, we adopted
cm-2. When no constraint on the plasma temperature was
obtained, we adopted 1 keV (value between brackets in Col. 9). The
X-ray luminosity in the 0.5-8 keV energy range corrected for
absorption in Col. (10) was computed assuming a distance of 140 pc
for the TMC. The X-ray fractional luminosity,
,
is given in the last column.
![]() |
Figure 5: Quantile diagram for pn ( top) and MOS1+MOS2 ( bottom). The grid of hydrogen column density and plasma temperature were computed following the prescriptions of Hong et al. (2004) using in SHERPA an absorption model ( WABS; Morrison & McCammon 1983) multiplied with a single optically thin thermal plasma model ( MEKAL; Mewe et al. 1995) with 0.3 times the solar elemental abundances. The background map in grey levels indicates the degeneracy level of grid parameters (see online Appendix C). The MOS1+MOS2 quantile diagram is used only for X-ray sources without available pn data. In the PN quantile diagram, asterisks and diamond mark the plasma parameters obtained from spectral fitting; the dashed line indicates the X-ray colour locus corresponding to the parameter uncertainties of source #4a (diamond; Table 3). The spectrum of #3 requires a two-temperature plasma model (see Table 3). |
Open with DEXTER |
For the faintest X-ray sources, we used a quantile analysis, the new
spectral classification technique for X-ray sources proposed by
Hong et al. (2004). This technique uses a quantile diagram based on the
X-ray colours:
,
and
,
where
,
with
the energy where the
net counts are x% of the source net counts in the
-
energy range. We take here
keV and
keV.
With the spectral model of Sect. 4.1, we calculated for
pn and MOS1+MOS2, from the Perl script provided by
Hong et al. (2004), a grid of hydrogen
column densities and plasma temperatures (
,
kT) in SHERPA
and the corresponding X-ray colours. We adopted on-axis RMF and ARF;
and for MOS1+MOS2, the MOS1 and MOS2 RMFs were averaged and the ARFs
added. Figure 5 shows the resulting quantile diagrams. Grid
lines indicate loci of iso-column densities and iso-temperatures.
We noted that the emission lines of the spectral model produce folds
in the parameter grid, introducing some degeneracy of the X-ray
colours; i.e. several different values of (,
kT) have the same
,
,
,
and hence the same
X-ray colours. The degeneracy of the quantile diagram is studied in details in
online Appendix C. The result of this analysis is the
degeneracy map plotted in grey levels in Fig. 5, which
indicates the number of different (
,
kT) values at each position in the
quantile diagram.
We calculated for all XMM-Newton sources pn (or MOS1+MOS2) X-ray colours with background
subtraction (Table 4) and errors following the
prescriptions of Hong et al. (2004). The MOS1+MOS2 quantile diagram is
used only for X-ray sources without available pn data.
We checked for sources #1, #4a and #5
that the plasma parameters derived from the pn X-ray colours are
consistent with the plasma parameters obtained from spectral fitting
marked with asterisks and diamond in Fig. 5. For comparison, we indicate in Fig. 5 the X-ray colour locus
corresponding to the parameter uncertainties found in the fitting of
the XMM-Newton/EPIC spectra of source #4a, the brightest X-ray source
(Table 3).
The spectrum of source #3 requires to be fitted by a two-temperature
plasma model, and has X-ray colours which are outside the grid
computed with a single-temperature plasma model.
Optical extinctions were used to disentangle temperature double
solutions (source #2 on pn; see online Appendix C), and to
estimate the hydrogen column density for sources with low constrained
X-ray colours (sources #2, #6 and #7 on pn, and #9 on MOS) using
the relation
cm-2 mag-1 (Cardelli et al. 1989; Vuong et al. 2003). Where negligible
optical extinction was measured, we adopted
cm-2. When there was no constraint on the plasma
temperature, we adopted kT=1 keV (Table 4).
Then, in XSPEC we
computed from these plasma parameters and the source's ARF and RMF the
X-ray luminosity, corrected for absorption in the 0.5-8 keV energy
range needed to reproduce the observed net count rate in the
0.5-7.3 keV energy range. For comparison, the X-ray luminosities
obtained by spectral fitting of sources #5a and #6 and the ones
obtained by quantile analysis show only a difference of 0.1 dex.
In the following, for BDs observed and detected twice, we use the logarithmic average of their X-ray luminosity. For CFHT-BD-Tau 1, we use the quiescent X-ray luminosity.
We calculated an upper limit to the X-ray luminosity of each
undetected BD in the 0.5-8 keV band (see Table 5).
Counts within 10
of the optical/near infrared position were
extracted from the summed EPIC soft-band (0.5-2 keV) image. The
expected number of counts in the absence of emission from the BD was
determined from the identical region of the corresponding
"reconstructed image'', which is composed of background and detected
sources output by the XMM-Newton Science Analysing System (SAS) source detection algorithm, EMLDETECT
(Güdel et al. 2007). An upper limit, at the 95% confidence level, to
the number of counts in the region from the BD was computed using the
prescription of Kraft et al. (1991), which accounts for Poissonian
fluctuations in the counts from the background, and those from the
BD
.
This results in an upper limit to the sum of counts collected by three different
detectors. The expected number of counts collected by
detector i, Ni, from a source of flux
is dependent on the
fraction of the total source counts that were collected within the
extraction region,
,
the exposure time,
,
and the source energy flux required to produce a count rate
of 1 count s-1, Ki, such that
.
Therefore, the upper limit to the source energy flux is
computed from the upper limit to the sum of source counts in all the
detectors,
,
as
.
A model source,
constructed at the position of the BD from the calibration point
spread function used by EMLDETECT, was used to calculate
,
and convolved with exposure images generated by the SAS command EEXPMAP to compute the effective on-axis exposure
time at the BD position,
,
on each detector. The
on-axis response of each detector was generated using the SAS
command ARFGEN and the appropriate canned response matrix.
Ki was calculated in XSPEC assuming an isothermal
coronal plasma of temperature 1 keV (11.6 MK) with metallicity 0.3 times
that of the solar corona (Anders & Grevesse 1989), and a hydrogen column
density obtained from the optical extinction using the relation
cm-2 mag-1(Cardelli et al. 1989; Vuong et al. 2003). Upper limits to X-ray luminosities in the 0.5-8 keV
energy range were calculated assuming a distance of 140 pc to the
TMC.
The median of the upper limits of our sample is about
erg s-1. However, the deep archival observation
of SU Aur helped to detect 2MASS J0455 down to half this
luminosity level.
To put the X-ray luminosity properties of the XEST BDs into context,
we compare them to the X-ray luminosity of low-mass stars in the XEST
(see Güdel et al. 2007).
For objects observed and detected twice, we use the
logarithmic average of their X-ray luminosity.
For TMC members that are unresolved multiple stars
in X-rays, we considered the total bolometric
luminosities. Figure 6 shows, for the XEST low-mass
stars (
,
which corresponds at
3 Myr to
;
Siess et al. 2000), protostars (Class I
sources) and BDs, the X-ray luminosities
compared to the bolometric luminosities. Considering only the X-ray
detected BDs and low-mass (proto)stars, we determined with the parametric
EM (Expectation-Maximization) algorithm implemented in ASURV
(Feigelson & Nelson 1985) the following linear regression fit:
,
with a standard deviation of 0.4 dex in
for the low-luminosity objects
(L*
10
). This relation is very similar to the
relations found for low-mass stars in other young clusters
(Preibisch et al. 2005a; Feigelson & Montmerle 1999). It is consistent with a linear
relation between X-ray and bolometric luminosity characterized by
which is valid from the
low-mass stars to the substellar regime. However, the bulk of the upper
limits of undetected BDs are below this average X-ray fractional
luminosity, which suggests a lower X-ray fractional
luminosity for the BDs with about
.
Taking into account upper limits of undetected BDs, the median of
for the XEST BDs is -4.0 (see
Sect. 7). The X-ray fractional luminosity of XEST BDs is hence
lower than the one of XEST low-mass stars.
To investigate the relation between X-ray fractional luminosity
and physical parameters when one moves from low-mass stars to the
substellar regime, we focus on objects with spectral type M0 or later, which
corresponds for an age of 3 Myr to masses and luminosities lower than
about 1
and 0.7
,
respectively. We used the temperature scale given in
Table 1 to convert M spectral types to effective
temperatures.
We determined source masses from effective temperature
and luminosities, using the pre-main-sequence tracks of
Baraffe et al. (1998). Masses were interpolated linearly along the
isochrones. For sources located above the 1 Myr isochrone, we
prolonged the convective tracks vertically. We attributed
0.02
to KPNO-Tau 4, i.e. the only BD above the
0.02
track. The source IRAS S04301+261 (spectral type
M0), which is located below the main-sequence was discarded.
We restricted our sample to stars which are not
multiple in order not to introduce extra assumptions when distributing the
X-ray flux to unresolved components. Figure 7 shows
the X-ray fractional luminosity versus mass for this
sample. Spearman's rank correlation coefficient computed with ASURV indicates a correlation with the probability of the null
hypothesis (i.e. no correlation)
.
The EM
algorithm yields the shallow linear regression fit:
,
with a
standard deviation of 0.5 dex in
.
This relation is consistent with the relation found by
Preibisch et al. (2005a) for the T Tauri stars (
)
of
the ONC. This relation implies that the X-ray fractional luminosity
decreases by a factor of about 3 from 1
stars to
0.03
BDs.
Table 5:
Upper limits to the BD X-ray luminosities. The hydrogen column
density in unit of 1021 cm-2 given in Col. (3) is obtained
from the optical extinction using the relation
cm-2 mag-1 (Cardelli et al. 1989; Vuong et al. 2003). Exposure times in
Col. (4) are for summed EPIC (pn+M1+M2) data, in units equivalent for
a pn-on-axis observation. Columns (5) and (6) give the upper limits at
the 95% confidence level for net counts in the 0.5-2 keV energy
range and the X-ray luminosity in the 0.5-8 keV energy range
corrected for absorption, respectively. The upper limit to the
X-ray fractional luminosity,
,
is given in
the last column.
![]() |
Figure 6:
X-ray luminosity vs. bolometric luminosity for the TMC
BDs (diamonds; and arrows for upper limits) and
the TMC members detected in the XEST (white and black dots
represent low-mass stars and protostars, respectively; Güdel et al. 2007). The
dotted lines indicate from bottom to top an X-ray fractional
luminosity,
![]() ![]() |
Open with DEXTER |
![]() |
Figure 7: X-ray fractional luminosity vs. mass for the young BDs of the TMC and single TMC members of the XEST. The symbols are as in Fig. 6. The vertical dotted lines indicates the stellar/substellar boundary. The solid and dashed lines show a linear regression fit and standard deviation, respectively. |
Open with DEXTER |
![]() |
Figure 8: X-ray fractional luminosity vs. effective temperature for the young BDs of the TMC and single TMC members of the XEST. The symbols are as in Fig. 6. The vertical dotted lines indicates the stellar/substellar boundary. The solid and dashed lines show a linear regression fit and standard deviation, respectively. |
Open with DEXTER |
![]() |
Figure 9: X-ray surface flux vs. effective temperature for the young BDs of the TMC and single TMC members of the XEST with spectral type M0 or later. The symbols are as in Fig. 6. The solid and dashed lines show a linear regression fit and standard deviation, respectively. |
Open with DEXTER |
Figure 8 shows the
X-ray fractional luminosity versus the effective temperature for our
sample ranging from 3840 K to
2500 K. Spearman's rank
correlation coefficient computed with ASURV indicates a
correlation with the probability of the null hypothesis (i.e. no
correlation)
.
The EM algorithm yields the linear
regression fit:
,
with a standard deviation of 0.5 dex in
.
This
relation implies that the X-ray fractional luminosity decreases by a
factor of about 3 from hot coronae of solar-mass stars to cooler
atmospheres of M9V BDs.
We computed the X-ray surface flux, i.e. the X-ray
luminosity divided by the surface area of the source, which is
computed from its bolometric luminosity and effective
temperature:
,
where
is the Stefan-Boltzmann constant.
Combined with the decrease of the
ratio with the
effective temperature, this formula implies that the X-ray surface
flux decreases with the effective temperature with a power-law slope
greater than 4. Figure 9 shows the X-ray surface flux
versus the effective temperature for our sample. The X-ray surface
fluxes range from
105 to
erg s-1 cm-2. For comparison, the
X-ray surface flux of the solar corona during the solar cycle varies
from
to
erg s-1 cm-2(Peres et al. 2004). Spearman's rank correlation coefficient indicates a
log-log correlation with a high confidence level (
P(0)=10-4). The
EM algorithm yields:
,
with a standard deviation of 0.5 dex in
.
The slope of this correlation is well steeper than 4 due to the
decline of the
ratio with the effective
temperature. This correlation predicts that a M9V BD should have an
X-ray surface flux in its corona
10 times weaker than in an M0 star. However, this is still 4 times higher than the X-ray surface
flux of the solar corona at the solar cycle maximum (Peres et al. 2004). A
similar trend was found by Preibisch et al. (2005b) for the BDs of ONC (see
Sect. 7 for further comparison between the XEST and the COUP
BD sample).
Studying a BD sample including young BDs (in the ONC, Ophiuchi, and IC 348) plus one intermediate BD (TWA 5B)
and one old (LP 944-20) field BD, Tsuboi et al. (2003) found a log-log correlation
between the X-ray fractional luminosity (value independent of the
distance assumption), and
(
), and noted that BDs
with
greater than 100
were not detected in
X-rays. We note that including the corresponding (six) X-ray upper
limits, this correlation vanishes. Moreover, this proposed correlation
is strongly biased by the low X-ray activity of LP 944-20 - detected
only during an X-ray flare (Rutledge et al. 2000) - which is likely
explained by the cooling of this old (500 Myr) field BD
(Stelzer et al. 2006).
![]() |
Figure 10:
X-ray fractional luminosity of TMC BDs vs.
![]() |
Open with DEXTER |
![]() |
Figure 11:
Equivalent width of H![]() ![]() |
Open with DEXTER |
measurements are available for all the XEST BDs but 2MASS J0455
(Table 1), allowing us to investigate the relationship
between X-ray fractional luminosity and H
emission from a
well-defined sample of young BDs. Figure 10 shows the
X-ray fractional luminosity versus
.
We detected in
X-rays BDs both with low and high
.
We compute Spearman's rank
correlation coefficient, and find no significant log-log correlation
between the X-ray fractional luminosity and
.
![]() |
Figure 12: Cumulative distributions of the X-ray fractional luminosities for nonaccreting (continuous line) and accreting (dashed line) BDs. In the top panel, diamonds and arrows indicate detections and upper limits, respectively of the two samples. Only Kaplan-Meier estimator error bars of the nonaccreting sample are shown to clarify the plot. Two-population statistical methods show that both samples are drawn from the same underlying distribution. |
Open with DEXTER |
The H
emission line cannot be used as a tracer of accretion in
low-mass stars and BDs without a priori knowledge on the limit of pure
chromospheric H
emission in these objects. The canonical limit
of 10
for
was first used to disentangle
pure chromospheric activity in Weak-line T Tauri Stars (WTTSs) and
H
emission line excess in Classical T Tauri stars (CTTSs)
produced by accretion. Martín (1998) proposed
limits depending on the spectral type.
Barrado y Navascués & Martín (2003) improved this criterion and extended it into
the substellar domain. This empirical criterion is the saturation at
,
based on observations of
nonaccreting stars in young open clusters, which corresponds
physically to the maximum amount of energy that can be released in
nonthermal processes by the chromosphere, i.e.
of the total emitted energy. Low-mass stars or BDs
exceeding this limit are accreting. Figure 11 shows
versus spectral type for the BDs of XEST. The
dashed line shows the saturation at
which increases from
24
at M 6 to
148
at L0 (Barrado y Navascués & Martín 2003).
MHO 4 and CFHT-BD-Tau 3 are just above this limit, but
high-resolution spectrum showed no indication of accretion
(Mohanty et al. 2005). Therefore, we classify these BD as
nonaccreting. For the other BDs, our classification based on
is in agreement with detailed high-resolution accretion
identification (Mohanty et al. 2005; Muzerolle et al. 2005b). We find 6 accreting and
8 nonaccreting BDs.
The X-ray detection rates of
accreting and nonaccreting BDs are 33%
27% (2/6) and
63%
36% (5/8), respectively
.
We compute with ASURV the cumulative
distributions of the X-ray fractional luminosities - including upper
limits - with the Kaplan-Meier estimator for accreting and nonaccreting
BDs (Fig. 12).
Two-population statistical methods provided by ASURV show that the
probability for the null hypothesis that both samples are drawn from
the same underlying distribution is 0.75. Therefore, the accreting and
nonaccreting BDs have similar X-ray fractional luminosities. The
median value of the X-ray fractional luminosities is -3.9 for the
nonaccreting BDs.
A large difference is found between accreting and nonaccreting
low-mass stars. In the TMC, the CTTSs are 2.2 times less luminous
than the WTTS (Güdel et al. 2006). A similar result was found in the ONC
where the CTTSs (selected with 8542
Ca II line) are 2.1
times less luminous than the WTTSs
(Flaccomio et al. 2003; Preibisch et al. 2005a). Preibisch et al. (2005a) proposed that
in accreting objects magnetic reconnection cannot heat the dense
plasma in mass-loaded accreting field lines to X-ray
temperatures. If true, this implies that coronal activity is less
affected by the accretion in BDs than in low-mass stars.
The median X-ray fractional luminosity of nonaccreting BDs in the XEST
is 4 times lower than the mean saturation value for rapidly rotating
low-mass (
)
field stars
(
;
Pizzolato et al. 2003), whereas in
ONC the median X-ray fractional luminosity of nonaccreting low-mass
stars is consistent with this saturation level (Preibisch et al. 2005a).
![]() |
Figure 13: Cumulative distributions of the X-ray fractional luminosities of XEST (continuous line) and COUP (dashed line) BDs. In the top panel, diamonds and arrows indicate detections and upper limits, respectively of the two samples. Only Kaplan-Meier estimator error bars of the nonaccreting sample are plotted to clarify the plot. Two-population statistical methods show that XEST BDs are more active in X-rays than COUP BDs. |
Open with DEXTER |
![]() |
Figure 14: X-ray fractional luminosity vs. spectral type for objects of type M 5 and later. Detections of late M field stars from Fleming et al. (1993) are shown as asterisks. The circles show low-mass stars of the TMC detected in X-rays (Güdel et al. 2007). Diamonds and thick arrows show BDs in the TMC. The other X-ray detected BDs (see Preibisch et al. 2005b, and references therein) are shown by gray filled squares. For very cool objects with strong flares, the values at flare peak are shown by triangles, connected by dotted lines to the quiescent emission. Some symbols have been slightly moved in spectral type to avoid overlaps. |
Open with DEXTER |
Preibisch et al. (2005b) in COUP studied the X-ray properties of young BDs
spectroscopically identified in the near-IR (Slesnick et al. 2004) having
spectral types between M6V and M9V. Eight out of the 33
BDs
were clearly detected as X-ray sources down to a detection limit of
.
The near-IR selection helped to find extincted
BDs (Slesnick et al. 2004). Consequently, the apparently low detection
rate of BDs in COUP is in many cases related to the substantial
extinction of these BDs. Considering only the ONC BDs with
mag reduces the median visual extinction from 5.6 to 2.2 mag,
and leads to an X-ray detection rate of 40% (6 out of 15 BDs),
similar to the X-ray detection rate observed for TMC BDs.
We compute the cumulative distributions of the X-ray fractional luminosities
of the whole sample of XEST BDs, and of the COUP BDs with
mag (Fig. 13). Two-population statistical
methods show that the probability for the null hypothesis that both
samples are drawn from the same underlying distribution is lower than
0.05. The former and the latter samples have median X-ray
fractional luminosity equal to -4.0 and -4.2, respectively.
Therefore, the XEST BDs are 1.6 times more active in X-rays than the COUP
BDs. The origin of this difference is not yet understood.
We discuss now the origin of the BD X-ray emission in the broader context of the X-ray emission of cool stars. Following Preibisch et al. (2005b), we compare the X-ray fractional luminosities of TMC BDs with other cool objects.
Figure 14 shows the X-ray fractional
luminosities of TMC objects with spectral type equal or later than M 5,
X-ray detected BDs (Preibisch et al. 2005b, and references
therein; Stelzer et al. 2006), the evolved BD Gl 569
Ba,b (Stelzer 2004), and the field BD LP 944-20
(Martín & Bouy 2002; Rutledge et al. 2000). Cool field stars with spectral type M 5
or later are also indicated: M field stars from
Fleming et al. (1993); very-low mass field stars VB 10 (Fleming et al. 2003),
LHS 2065 (Schmitt & Liefke 2002), HC 722
(Preibisch et al. 2005b; Slesnick et al. 2004,2005). All these cool objects show
a similar level of X-ray fractional luminosities ranging from
10-3 to
10-4 (note that, compared to
Fig. 8, the spectral type range has been reduced
by a factor of two in Fig. 14).
Based on Baraffe et al. (1998) evolutionary tracks, objects of spectral
type M 7 with an age of 1 Gyr are not BDs, but low-mass stars twice
as massive as a typical TMC BD having an M 7 spectral type and an age of 3 Myr.
Moreover, such very cool stars also have surface gravities about 40
times higher than in a typical TMC BD.
This shows, as found by Preibisch et al. (2005b) from the COUP BDs, that
the X-ray activity of BD coronae is not strongly dependent of the BD
mass and the BD surface gravity. Therefore, this implies that the
relation that we have obtained, for the young objects of TMC with
spectral type M, between X-ray fractional luminosity and mass
(Fig. 7) is in fact the consequence of the more
fundamental relation between X-ray fractional luminosity and effective
temperature (Fig. 8). This latter relation agrees
with the overall result (field dwarfs and young BDs) in
Fig. 14: of the 15 sources shown with spectral
types M8.5V or later, only 4 have any detected quiescent emission; the
rest are either not detected at all or (in 3 cases) detected only
during strong flares.
By considering a subsample of
evolved BDs (from 1 Myr to
1 Gyr) with a small
range of masses (0.05-0.07
)
but with effective
temperatures ranging from 3000 K down to 1000 K, Stelzer et al. (2006) found
evidence for a similar decline of the fractional X-ray luminosity with the
effective temperature; showing that the atmospheric temperature plays
a crucial role in determining the level of X-ray activity.
The growing evidence for a decline of the coronal activity with the
effective temperature is analogous to the decline in chromospheric
(H)
activity seen in field dwarfs at similar spectral types,
which show a slow decline from mid- to late M in
and a sharp drop-off around M9V (Mohanty & Basri 2003).
Although a definitive turn-over of the fractional X-ray luminosity
toward L-type BDs is not proven for the present sample
due to sparse statistics, we note that the fractional luminosity
declines from low-mass stars to M-type BDs, and as a sample, the
BDs are less efficient X-ray emitters than low-mass stars.
We thus conclude that while the BD atmospheres observed here are
mostly warm enough to sustain coronal activity, a trend is seen
that may indicate its gradual decline due to the drop in
photospheric ionization degree (Mohanty & Basri 2003).
With the XEST we detected 9 (out of 17) young BDs with spectral
type ranging from M6.25V to M8V; 7 BDs are detected here for the
first time in X-rays. This BD sample surveyed in X-rays allowed us to
investigate the magnetic activity in the substellar regime. We recovered a
well-known relation between X-ray and bolometric luminosity for stars (with
L* 10
)
and BDs detected in X-rays, which is
consistent with a mean X-ray fractional luminosity
.
For the XEST BDs, the median of
(including upper limits) is -4.0. The X-ray
fractional luminosity of XEST BDs is hence lower than the one of XEST
low-mass stars.
A shallow relation is found between X-ray fractional luminosity and mass. We show that the X-ray fractional luminosity decline by a factor of about 3 from hot coronae of solar-mass stars to cooler atmospheres of M9V BDs. Consequently, a relation is found between X-ray surface flux and effective temperature, which implies a decrease of about one magnitude in the X-ray surface flux from a M0V star to a M9V BD.
No significant relation is found between the X-ray
fractional luminosity and
.
We used
to identify accreting and nonaccreting BDs.
Accreting and nonaccreting BDs have a similar X-ray fractional
luminosities. The median X-ray fractional luminosity of
nonaccreting BDs are about 4 times lower than the mean
saturation value for rapidly rotating low-mass field stars. BDs have higher X-ray
fractional luminosities in the TMC than in Orion.
We confirm, as previously observed in the Cha I (Stelzer et al. 2004) and in the ONC (Preibisch et al. 2005b), that there is no dramatic change of the magnetic activity at the stellar/substellar boundary. Young BDs of spectral type M are sufficiently warm to sustain an active corona. The young BDs in the TMC, with a median spectral type of M 7.5, have on average an X-ray surface flux which is 7 times higher than the one observed in the solar corona at the solar cycle maximum.
Deeper X-ray observations of the coolest M-type BDs in the TMC are needed to investigate a possible turn-over of the fractional X-ray luminosity of TMC BDs around spectral type M9V.
Note added in proof
2MASS J04335245+2612548, a TMC BD with spectral type M8.5V (Luhman et al. 2006) which was reported after the acceptance of this paper, is neither detected in X-rays by XMM-Newton/EPIC or in the U-band by the XMM-Newton optical/UV monitor; it is located on the X-ray PSF wings and the U-band "smoke ring'' of IT Tau (XEST-18-030). Therefore, adding this newly identified BD to our TMC BD sample wouldn't change our conclusions.
Acknowledgements
We thank the anonymous referee for his constructive comments that helped to improve this paper; Beate Stelzer for a careful reading of the manuscript; and the International Space Science Institute (ISSI) in Bern for significant financial support of the project team. This research is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. X-ray astronomy research at PSI has been supported by the Swiss National Science Foundation (grants 20-66875.01 and 20-109255/1). M.A. acknowledges support from NASA grant NNG05GF92G.
Table A.1:
Archival Chandra observations which surveyed
serendipitously the TMC BDs. Column (1) gives the name of the
instrument used: "I'' (
field of view) and "S''
(
field of view) for imaging and spectroscopy ACIS
CCD (Garmire et al. 2003), respectively. Column (8) give the name of the
TMC BD whitin the Chandra field of view. Column (9) indicates whether we use this archival
data to supplement the XEST or not (see following notes). Notes:
in the ACIS-S observation #3364, the TMC BD KPNO-Tau 2 was located
13
off-axis on ACIS-S2, and it is not detected; we obtained here a better
constraint on the X-ray luminosity of this BD using the sum of the
two XMM-Newton exposures, rather than this ACIS-S observation. In the ACIS-S
observation #4488, the TMC BD 2MASS J0421 was located
5.7
off-axis on ACIS-S3, the pipeline detection algorithm found
no source at this location.
Table A.1 gives the log of Chandra observations which surveyed serendipitously the TMC BDs, and used to supplement the XEST.
Data reduction of archival Chandra data were performed with CIAO
software package version 3.2.1, XSPEC version 11.3, the Penn State charge transfer
inefficiency (CTI) corrector
version 1.45, and the acis_extract package version 3.67.
We followed Penn State University (PSU) procedures
(e.g., Getman et al. 2005) starting from data reduction from Level 1
event files provided by the Chandra X-Ray Center (CXC). We used
the PSU CTI Corrector (Townsley et al. 2002) to correct partially the data
for CCD CTI caused mainly by radiation damage at the beginning of the
Chandra mission. Source detections were performed with the CIAO's task wavdetect. 2MASS counterparts were used to
correct small boresight errors. The acis_extract package was
then used to extract source photons, estimate local background,
construct source and background spectra, compute redistribution matrix
files (RMFs) and auxiliary response files (ARFs), construct light
curves and time-energy diagrams, and perform automated spectral
grouping and fitting.
Figure B.1 shows the X-ray light curves of the brightest TMC BDs, which are consistent with quiescent activity. Figure B.1 shows X-ray light curves of BDs observed at two different epochs.
In contrast to the continuum models used as illustration of quantile diagram by Hong et al. (2004), which connect one single couple of physical parameters to unique X-ray colours, we found that the spectra of optically thin plasma (i.e. a continuum bremsstrahlung plus emission lines) may have several couple of physical parameters producing the same X-ray colours. To use the quantile analysis safely, we need to reckon for each physical parameter the number of possible values per X-ray colour.
We divided first the quantile diagram in a regular grid with a
resolution. Then, the uneven grid of X-ray colours of
each physical parameter was projected to this regular grid using
spline interpolation. Each physical parameter therefore generates a
3-dimensional surface regularly sampled in the X-ray colour
plane. Finally, we found numerically the number of intersections
between this surface and a line of constant colour, i.e. the
degeneracy level.
Figures C.1 and C.2 show the
result for the hydrogen column density and the plasma temperature grids,
respectively, computed with WABS MEKAL plasma model
with 0.3 times the solar elemental abundances and XMM-Newton EPIC pn RMF and ARF. The left
pannels show the degeneracy maps, which are nearly identical when one
considers the uneven sampling of the model parameter; in
Fig. 5 we plotted conservatively the maximum of these two
degeneracy maps. The right pannels show for each parameter the
corresponding 3-dimensional surface in the X-ray colour plane with
folds producing multiple solutions. The quantile diagram is mainly composed
of areas where X-ray colours correspond to a unique couple of
parameters, but areas with multiple solutions - from 2 to 5 - are also
present. However, the area surface decreases with the number of
multiple solutions.
Figure C.3 shows an illustration of degenerated
solutions if y=1.3. For example if
there is
a unique solution for
and kT. If
there is a double solution: a first solution with low
and high
kT, and a second solution with high
and low
kT. Generally speaking, if there are N solutions of
and kT ordered by increasing values, the solution number n(with n ranging from 1 to N) is defined as the
-value
of
and the
-value of kT. Degeneracy on
both parameter can then be fixed using extra knowledge on only one
parameter.
We used the knowledge of the optical extinction to derive the value
of the hydrogen column density using
the relation
cm-2 mag-1 (Cardelli et al. 1989; Vuong et al. 2003). The best estimate of the
temperature is then the temperature solution corresponding to the
hydrogen column density solution which is the closest to this value.
![]() |
Figure C.3:
Illustration of degenerated solutions. Left: horizontal cut of
the ![]() |