A&A 468, 529-540 (2007)
DOI: 10.1051/0004-6361:20066318
The XMM-Newton extended survey of the Taurus molecular cloud
M. Güdel1 -
S. L. Skinner2 -
S. Yu. Mel'nikov3 -
M. Audard4,
-
A. Telleschi1 -
K. R. Briggs1
1 - Paul Scherrer Institut, Würenlingen and Villigen,
5232 Villigen PSI, Switzerland
2 -
CASA, 389,
University of Colorado,
Boulder, CO 80309-0389,
USA
3 -
Ulugh Beg Astronomical Institute,
Academy of Sciences of Uzbekistan,
Astronomicheskaya ul. 33,
Tashkent 700052, Uzbekistan
4 -
Columbia Astrophysics Laboratory,
Mail Code 5247, 550 West 120th Street,
New York, NY 10027,
USA
Received 29 August 2006 / Accepted 8 November 2006
Abstract
Context. The generation of X-rays in accreting T Tauri stars (TTS) is thought to be predominantly due to energy dissipation in magnetic fields, but alternative X-ray generation mechanisms have been proposed, such as heating in accretion shocks near the stellar surface, or in shocks forming in jets.
Aims. We test models and trends discussed in the literature using X-ray data from the classical TTS T Tau.
Methods. High-resolution spectroscopy from the Reflection Grating Spectrometers on XMM-Newton is used to infer electron densities, element abundances and the thermal structure of the X-ray source. We also discuss the ultraviolet light curve obtained by the Optical Monitor, and complementary ground-based photometry. A high-resolution image from Chandra constrains contributions from the two companions of T Tau N.
Results. The X-ray grating spectrum is rich in emission lines, but shows an unusual mixture of features from very hot (30 MK) and very cool (1-3 MK) plasma, both emitted by similar amounts of emission measure. The cool plasma confirms the picture of a soft excess in the form of an enhanced O VII/O VIII Ly
flux ratio, similar to that previously reported for other accreting TTS. Diagnostics from lines formed by this plasma indicate low electron densities (
1010 cm-3). The Ne/Fe abundance ratio is consistent with a trend in pre-main sequence stars in which this ratio depends on spectral type, but not on accretion.
Conclusions. On the basis of line density diagnostics, we conclude that the density of the cool "soft-excess'' plasma is orders of magnitude below that predicted for an accretion shock, assuming previously determined accretion rates of (3-6)
yr-1. We argue that loading of magnetic field lines with infalling material suppresses the heating process in a part of the corona. We thus suggest that the X-ray production of T Tau is influenced by the accretion process although the X-rays may not form in the bulk of the accretion footpoints.
Key words: stars: coronae - stars: formation - stars: pre-main sequence - X-rays: stars - stars: individual: T Tauri
Low-mass pre-main sequence stars maintain magnetic fields that are thought to provide the energy for high-energy particles and radiation as shown by their strong X-ray and radio emission (for reviews, see, e.g., Güdel 2004; Feigelson & Montmerle 1999; Güdel 2002; Feigelson et al. 2007). Vigorous sources of X-ray radiation have been found among class I protostars (e.g., Imanishi et al. 2001; Koyama et al. 1994), accreting classical T Tau stars (CTTS, e.g., Feigelson & DeCampli 1981; Walter & Kuhi 1981), and weak-line T Tau stars (WTTS; e.g., Neuhäuser et al. 1995; Walter 1986) that are no longer surrounded by a thick, accreting disk. A number of observational features have supported a picture in which X-ray emission is due to magnetic coronal processes similar to mechanisms seen in the solar corona: rotational modulation points to a concentration of hot plasma in various magnetic active regions (Flaccomio et al. 2005); flares heat plasma on time scales of minutes to hours, reaching temperatures of a few tens of MK (Wolk et al. 2005); and electron temperatures continuously exceed 10 MK (e.g., Skinner et al. 2003; Imanishi et al. 2001); such temperatures cannot be achieved by liberating gravitational energy from infall, nor by shocks in moving gas (e.g., jets) commonly observed around young stars at other wavelengths. Also, high-energy electrons observed through their gyrosynchrotron emission give direct evidence for particle acceleration, for which magnetic reconnection processes provide an ideal basis.
Nevertheless, accreting pre-main sequence stars add circumstellar disks and accretion flows to the picture of the stellar environment; stellar magnetic fields may reach out to the inner border of the disk from where they guide lightly ionized mass to the stellar surface (Calvet & Gullbring 1998). WTTS and CTTS are therefore ideal objects to test to what extent the stellar environment controls coronal magnetic fields and energy release. Except for disks and mass accretion streams, the two types of stars appear to be identical in nature, occupying the same region in the Hertzsprung-Russell Diagram (HRD).
Evidence for some sort of influence of the disk environment onto X-ray production
is broad: i) star-disk magnetic fields may brake stellar rotation, thus damping the
magnetic dynamo and inducing weaker overall magnetic activity. In the Taurus Molecular
Cloud, for example, CTTS rotate on average less rapidly than WTTS (characteristic
rotation periods of
d and 2-4 d, respectively; Bouvier et al. 1993).
Many TTS may, however, emit at the empirical saturation limit,
,
for rotation periods below 10 d (Preibisch et al. 2005), so that
a difference between the dynamo operation of CTTS and WTTS may be difficult to identify; ii) despite possible saturation effects, it has now been clearly established that both the X-ray
luminosity,
,
and
are, for otherwise similar stellar characteristics,
lower by a factor of
2 for CTTS than for WTTS
(Neuhäuser et al. 1995; Preibisch et al. 2005; Telleschi et al. 2007a). It is unclear in what sense accretion is
responsible for this result; iii) on the other hand, the electron temperature
of the X-ray emitting plasma may be higher in CTTS than in WTTS (Tsujimoto et al. 2002; Imanishi et al. 2001; Telleschi et al. 2007a);
iv) specific changes in the time behavior of X-rays were also observed in outburst phenomena
ascribed to strong accretion events in FU Ori-type and EX Lupi-type variables (e.g., a strong
flux increase of hard emission during outburst, Kastner et al. 2006,2004a; or a softening of
the X-rays during outburst, Audard et al. 2005).
Recently, a more direct role of accretion has been suggested for the generation
of X-rays in CTTS, namely emission from gas that has been shock-heated in the
accretion flow near the stellar surface. While free-fall velocities of a few 100 km s-1
on CTTS are sufficient to heat shocks to a few MK, a problem is X-ray absorption by the gas
stream itself, and by the surrounding photosphere if the shock is buried below the surface
(Drake 2005; Calvet & Gullbring 1998; Lamzin 1999). However, the temperatures attained in such shocks
are clearly not sufficient to explain the extremely hot X-ray sources of CTTS with temperatures
exceeding 10 MK.
On the other hand, in the exceptional case of the CTTS TW Hya, X-ray spectroscopic
evidence for a predominantly cool plasma led Kastner et al. (2002) to suggest
that X-ray production by accretion may be relevant in this star. This conjecture has
been further supported by evidence for unusually high electron densities of
cm-3 derived from density-sensitive line flux ratios of O VII and Ne IX in high-resolution X-ray spectra of the CTTS TW Hya, BP Tau, and V4046 Sgr
(Kastner et al. 2002; Robrade & Schmitt 2006; Günther et al. 2006; Schmitt et al. 2005; Stelzer & Schmitt 2004), supported by flux ratios
of Fe XVII lines (Ness et al. 2005). Stelzer & Schmitt (2004) further
argued that high Ne/Fe and N/Fe abundance ratios may point to accreting gas that is
depleted of Fe, presumably because Fe (and several other elements) easily condense onto
dust grains further out in the disk, leaving an accreting gas phase enriched by Ne and N.
Two large survey projects of X-ray emission in star-forming regions have addressed the
above questions anew. Preibisch et al. (2005) found no support for an accretion scenario in
in the Chandra Orion Ultradeep Project (COUP); the bulk of the
observed X-ray emission is related to very hot, coronal plasma. Accretion may in fact suppress
the X-ray luminosity by mass loading of magnetic fields. High-resolution X-ray spectroscopy
collected in the framework of the XMM-Newton Extended Survey of the Taurus Molecular Cloud
(XEST) has shown a series of new facets (Telleschi et al. 2007b):
i) very high (1011 cm-3) electron densities are not recorded in all
accreting pre-main sequence stars, posing problems for the accretion scenario even in the most
favorable case of 100% accretion surface filling; ii) abundance "anomalies'', in particular a high Ne/Fe abundance ratio, are common to many TTS, regardless of whether or not
they are accreting. On the other hand, the Ne/Fe abundance ratio seems to depend on the stellar
spectral type; iii) while the bulk plasma observed in CTTS shows very similar thermal
properties as plasma in WTTS, there is soft excess emission in accreting stars; this
excess is visible as (and is defined by) an enhanced ratio of the O VII/O VIII
Ly
line fluxes compared with WTTS or zero-age main-sequence (ZAMS) stars. Because
the maximum formation temperatures of O VII and O VIII are 2 MK and 3 MK, respectively (and O VIII is formed up to considerably higher temperatures),
an enhanced ratio indicates the presence of additional plasma below 3 MK.
Observation of this excess requires high-resolution spectroscopy that resolves emission
lines formed at 1-4 MK such as lines from O VII or N VI and N VII.
To further our understanding of these issues, detailed case studies
are needed. Ideally, a bright source will produce detailed
high-resolution grating spectra of the line-dominated range below 1 keV.
This allows us to assess anomalies in the thermal structure and electron densities from line ratios.
To access lines formed at low (1-4 MK) temperatures (located predominantly at wavelengths
longward of 17 Å), the photoelectric absorption should be moderate (hydrogen column density
of no more than a few times
1021 cm-2). Only few CTTS fulfill these requirements. One of the most outstanding
among them is the CTTS T Tau N. We report detailed X-ray spectroscopic
and imaging studies of this object, also searching for emission related to its two embedded
companions, T Tau Sab.
The specific goal of this paper is to use T Tau as a test case for hypotheses of X-ray production related to magnetic and accretion processes. We add complementary information from optical and ultraviolet observations. The plan of the paper is as follows. In Sect. 2, we introduce the T Tau system, and describe our observations in Sect. 3. Section 4 presents our imaging studies, Sect. 5 introduces light curves, and Sect. 6 presents results from our spectroscopic analysis. Section 7 contains a discussion and our conclusions.
Historically, the optically revealed northern component, T Tau N, was the prototype
of accreting low-mass pre-main sequence stars. Subsequently, a strongly variable IR companion was
found
to its south, now designated as T Tau S.
It was soon recognized to be a protostar probably co-eval with
T Tau N (Dyck et al. 1982) and dominating the system flux at
m.
The southern component was later found to be double, with a separation of
only
or 13 AU (Duchêne et al. 2002; Koresko 2000). An orbital analysis of
astrometric radio and infrared data infers masses of
and
for Tau Sa and Sb, respectively (Duchêne et al. 2006),
suggesting that the former is an embedded Herbig Ae star. The mass of T Tau N is also higher than
the masses of typical TTS, namely
(Güdel et al. 2007a),
and its rotation period is P = 2.8 d (Herbst et al. 1986). The distance to T Tau is
140 pc
(Loinard et al. 2005).
The T Tau triple may be at a critical transition stage between protostars and
fully optically revealed TTS; some of its observational properties appear to be
related to its particular orientation that places us into an optimal position
for observations. van Langevelde et al. (1994a) and van Langevelde et al. (1994b) found evidence
for infall of molecular material and a collimated outflow, pointing at 13
to the observer. Solf & Böhm (1999) suggest a pole-on view for T Tau N, which is also clearly supported
by mm observations of a well-developed, near-face-on massive accretion disk
(Akeson et al. 1998). The latter authors point out that the disk is similar in mass
to the solar nebula in which our planetary system formed. The pole-on view
of the N component results in a relatively low extinction, while from larger angles
all T Tau components would appear to be embedded accreting young stars. A CO emission ring, identified with the outer walls of a bi-conical outflow region, provides
evidence that we are looking down into a low-extinction cavity which the stellar
wind has nearly swept free of molecular gas (Momose et al. 1996).
T Tau Sa appears to be surrounded by an opaque disk seen nearly edge-on
(Duchêne et al. 2005,2006). The disk then is no more than 3 AU in radius given the
close orbit of the Sb companion (with a semi-major axis of 11.6 AU and a periastron distance
of
6 AU).
The T Tau system has been detected as a radio source
with properties of an extended wind but also of non-thermal radio emission
(Phillips et al. 1993; Skinner & Brown 1994), including gyrosynchrotron emission from an outflow (Ray et al. 1997) and strongly polarized emission from magnetic fields inferred to be
of kGauss strength (Smith et al. 2003; Loinard et al. 2005).
The T Tau system was well detected in X-rays with the ROSAT PSPC
(Neuhäuser et al. 1995; Stelzer et al. 2000), with a quiescent count rate of (0.03
0.009) cts s-1, and multiple energetic and rapid
flares developing on
1 h timescales and reaching up to 0.2 cts s-1.
Table 1 summarizes basic parameters of our XMM-Newton and Chandra X-ray observations, including the centroid coordinates of the X-ray images of T Tau. Both exposures comprise data of four stellar/substellar systems that are thought to be members of the Taurus Molecular Cloud: these are, apart from the T Tau system, the M0 CTTS IRAS 04187+1927, the WTTS RX J0422.1+1934, and the brown dwarf 2M J04221332+1934. All have been detected as X-ray sources. Further details on the former two stars are summarized in Güdel et al. (2007a), while the brown dwarf is described in Grosso et al. (2007).
Table 1: X-ray observing log and coordinates of T Tau.
We obtained an 83 ks exposure of the T Tau system on
August 15-16, 2005, with XMM-Newton (Jansen et al. 2001) using standard observing
setups
.
Specifically, the European Photon Imaging Cameras (EPIC)
obtained images with an angular resolution of
(FWHM), using
the full window with a diameter of 30
for the PN camera (Strüder et al. 2001), while
the MOS1 and MOS2 cameras (Turner et al. 2001) observed in the large-window (W3) and small-window (W2)
configurations that constrain the central CCD to a field of 5.5
5.5
and
1.8
1.8
,
respectively. These constrained window modes were
chosen to prevent the detectors from pile-up effects in case of strong flares (but no strong
flares that would induce pile-up were actually observed). All EPIC cameras
used the medium blocking filter.
The two Reflection Grating Spectrometers (RGS, den Herder et al. 2001) were operating in spectroscopy mode, providing X-ray spectra with a resolution of approximately 60-70 mÅ (FWHM) in the wavelength range of 5-35 Å. Two chip failures earlier in the mission left a gap in the RGS1 spectrum in the range of 10.6-13.8 Å and in the RGS2 spectrum in the range of 20.0-24.1 Å.
The Optical Monitor (OM, Mason et al. 2001) observed in imaging and high-time resolution
mode with the UVW1 ultraviolet filter inserted. The latter transmits light roughly in the
2500-3500 Å range. The imaging field of view is a square with a size of 17
17
.
The high-time resolution recording applies, however, only to a central window of
10.5
10.5
in which T Tau was positioned.
The X-ray data analysis of the XMM-Newton observations is described in detail by Güdel et al. (2007a). For details on the OM data reduction, see Audard et al. (2007). The data reduction made use of standard SAS v.6.1.0 software (i.e., the epchain, emchain, rgsproc, omfchain, and omichain tasks). Spectra and images were extracted with the evselect task. For the EPIC data, we used canned response matrices appropriate for SAS v.6.1.0, whereas ancillary response files (ARFs) were generated using the SAS task arfgen. The RGS response matrix was generated with the rgsrmfgen task. Background radiation fluctuated between a high and a low level during the observation, each occupying approximately half of the time. The high particle radiation can affect accurate background subtraction in both RGS, in particular in long-wavelength regions where the observed, true continuum is low and suppression of the source spectrum by photoelectric absorption may be severe. We have experimented with either extracting only low-background data, or with constraining the width of the source extraction region in the RGS detector plane. We found that extracting 70% (down from standard 90%) of the PSF width was optimum, and this choice was preferred to selecting in time while keeping the extraction width at the standard value.
The nearby WTTS RX J0422.1+1934 was bright enough (amounting to 41% of the count rate of T Tau in the EPIC detectors, see Güdel et al. 2007a) that it slightly contaminated the background region on the RGS detector. The source region of this star was therefore excluded from the background region of T Tau. It does not directly contaminate the T Tau spectrum, however, as it is located at the edge of the RGS detector, far outside the standard spectral extraction width for T Tau.
We further obtained a Chandra X-Ray Observatory (Weisskopf et al. 2002) exposure of 26 ks duration, using
the High Resolution Camera (HRC-I, Murray et al. 2000) with a spatial resolution of
but no
energy resolution in the X-ray range. The HRC field of view is about 30
in diameter.
The data were reduced using standard CIAO v.3.0.2a software
tools
.
New level 1 events were created with the hrc_process_events task, and level 2 events
after applying filters and good time interval (GTI) selection. A light curve
analysis revealed no significant variability. A formal source detection was performed
with the wavdetect algorithm in order to obtain precise centroid coordinates for
the T Tau source.
We obtained ground-based optical observations of T Tau during three nights around the Chandra observations. Photometric monitoring was performed in the UBVR bands with the 48-cm AZT-14 reflector at the Mt. Maidanak Observatory during the nights of 2003 October 22/23, 23/24, and 24/25. For these observations we used a single-channel photoelectric photometer in pulse-counting mode, and the photometric data were reduced to Johnson's international system. We observed HD 27560 as a comparison star and HD 284415 as check star. Other observational information can be found in Mel'nikov et al. (2005).
The T Tau triple with a separation of 0
7 between the N and Sab components cannot
be resolved by the XMM-Newton EPIC cameras. The Chandra HRC image, on the
other hand, could marginally reveal binary structure if the N and S components showed
comparable flux levels. The centroid coordinate of the stellar X-ray image
(Table 1) is very close to T Tau N.
Figure 1 shows the HRC-I image of the T Tau source. We modeled
a point-spread function (PSF) using the CIAO task mkpsf. For this, we used a circle
with a radius of 15 pixels (with a pixel size of 0.13
)
around the centroid.
The PSF was then
fitted to T Tau image in the Chandra sherpa software (Freeman et al. 2001)
(the version implemented in CIAO version 3.3.0.1 was used). The model also contained
a spatially constant background. We found the centroid to be slightly
but insignificantly offset from the true stellar position
(offset "X-ray minus true'' of
in RA and
in declination
).
Next, two PSFs were fitted, now assuming two sources with fixed (offset-corrected)
positions, the correct relative coordinates, and the same PSF width.
Although a very faint extension toward the S is suggested in the Chandra X-ray image, the best-fit amplitude at this position is only
% of the amplitude of the T Tau N model,
and is therefore at best a marginal detection (at the 2.4
level).
In summary, the Chandra HRC image (Fig. 1) shows that T Tau N is the
principal contributor to the X-ray source. This finding is
supported by spectral analysis below that shows low photoelectric absorption compatible
only with T Tau N. A small contribution from T Tau S cannot, however, be excluded. Because its visual extinction is
mag (for T Tau Sb, Duchêne et al. 2002), we expect a hydrogen column density of
1022 cm-2, assuming standard
gas-to-dust ratios (Vuong et al. 2003). Some counts could therefore be expected at energies above
1.5-2 keV only. We estimated an X-ray luminosity of T Tau S
with the PIMMS software
,
adopting
1022 cm-2, an electron temperature of either
7 MK or 28 MK (corresponding to either of the higher temperatures found for T Tau N), and 6.1% of the total counts. We found
1030 erg s-1 and
1030 erg s-1, respectively, for the two temperatures.
The true values may be higher if
is higher, but lower if we overestimated
the count rate contributions in the marginal fit. In any case, given the small contribution compared to T Tau N and lack of strong significance, we will assume that the observed XMM-Newton spectra are from T Tau N.
![]() |
Figure 1:
The T Tau source on the Chandra HRC-I detector (pixel size:
![]() ![]() ![]() ![]() |
Open with DEXTER |
Figure 2 shows the X-ray and ultraviolet light curves from the XMM-Newton observations.
![]() |
Figure 2: Light curves of the T Tau source, extracted from the XMM-Newton observations. From top to bottom: combined EPIC PN+MOS1 light curve for soft (0.4-1 keV) counts; hard (1-7.3 keV) counts; hardness ratio, i.e., ratio between hard and soft count rates; total light curve (0.4-7.3 keV; the arrows mark the times between which both EPIC cameras were operational); and the OM UVW1 light curve (dots illustrate bin count rates before manual correction; see text for details). Horizontal bars give time intervals for each OM exposure. |
Open with DEXTER |
We extracted the combined EPIC PN and MOS1 background-subtracted light curves for the softer and harder portion of
the spectrum (energies of 0.4-1.0 keV and 1.0-7.3 keV, respectively; MOS2 was not used because of difficulties in defining a background curve in the small-window mode).
These light curves are shown
in the upper two panels of Fig. 2. Both count rates are slowly decaying, probably from a flare with a peak that appears to have been recorded at the beginning of the observation. The harder flux decays more rapidly but then also levels off earlier, leading to a hardness increase in the second half of the observation, illustrated in the third panel that shows the ratio "hard:soft''. The fourth panel shows the light curve for the entire spectral range. Note that during the first 1350 s, the PN camera did not record data (interval before the time marked by the left arrow). We scaled the MOS1 light curve by an average factor to
the estimated PN+MOS1 count rate. Similarly, only the PN observed during the last 2830 s (time after the second arrow).
The light curve has again been scaled from the PN count rate to the expected total count rate.
The OM UVW1 light curve is shown in the bottom panel. It was recorded in a series of exposures, each with a duration of 3540 s, followed by a gap of 340 s. A few exposures were not successfully recorded.
The high count rate induced large coincidence losses which led to inaccurate calibration
of the incident count rate in some exposures, visible as obvious discontinuities with respect to the count rates in previous or following exposures. The affected exposures were manually corrected by multiplying their count rates
to produce a smoothly connected light curve (Fig. 2; the uncorrected fluxes
are also shown; the count rates in the first exposure were not corrected and are uncertain given the long time gap before the second exposure).
Three features are to be mentioned: i) there is no detailed correlation between the X-ray and the UV light curves; ii) both curves show, however, a general decay in the first half of the observation, perhaps as a consequence of
the weak flare seen in X-rays. The time of the UV peak is uncertain, given the unknown correction factor for the first exposure. It may have occurred before the start of the observations; iii) the UV light curve increases in the second half of the observation, in contrast to the X-rays. However, the X-ray hardness also increases,
and both UV flux and hardness show a peak around t = 6
104 s. It is possible that both features relate
to a flare-like heating event that produced both hot coronal plasma (seen in the harder X-ray band) and enhanced emission in the chromosphere and transition region (seen in the UVW1 band).
![]() |
Figure 3: Optical UBVR light curves of T Tau, collected during three consecutive nights (2003 October 22-25). The abscissae give JD-2 452 900.0, and the ordinates are in magnitudes. Both day-to-day variations and short-term fluctuations (in the U band) can be seen. |
Open with DEXTER |
The light curves in UBVR obtained at Mt. Maidanak Observatory during the ground-based monitoring period in 2003 around the observing campaign with Chandra are presented in Fig. 3. No overlap with the Chandra observing interval was achieved. These optical observations nearly cover the 2.8-day photometric period of T Tau (Herbst et al. 1986) which corresponds to the rotational period of the star. Our data show day-by-day variability of T Tau: while on October 22/23 the brightness of the star is rather stable, the brightness must have risen in all bands between the first and the second observation; prominent slopes can be seen in the light curves during the next two days, in opposite directions, first decreasing in brightness, then increasing again. Given that the rotation period of the star is 2.8 d, at least two minima and two maxima must have occurred during the stellar rotation covered by these observations. Magnetically active stars indeed often show two longitude intervals with preferred magnetic activity (e.g., Berdyugina et al. 1998). The amplitude of this variability seems to increase with decreasing wavelength and therefore reaches its maximum for the U band. This is particularly evident during the October 24/25 observation when the U-band curve approaches the B-band brightness. This day-by-day variability is similar to the long-term variations in the XMM-Newton UV light curve.
![]() |
Figure 4: Statistical analysis of variability in the U band light curve of T Tau ( top) and the check star HD 284415 ( bottom). The solid histograms give the number of measurements with a given magnitude offset from the average (after removing a slow trend for T Tau). The dashed lines mark the maximum rms error for a measurement at the average magnitude. |
Open with DEXTER |
![]() |
Figure 5: Fluxed, combined RGS1+2 spectrum. The background has not been subtracted, but its spectrum is overplotted as a red, dotted histogram, and also separately shown in the narrow, lower panels on the same flux scale. The bin width is 58 mÅ. The brightest line features contributing to the model fit are labeled. A few important features that have not been detected are marked with dotted lines. |
Open with DEXTER |
Another type of variability which is suggested in the U-band curve is defined by fast small-amplitude variations. In order to estimate the accuracy of our observations, we calculated histograms of the brightness distribution for both the check star and T Tau itself (Fig. 4). For the second and third T Tau observation, we first removed a slow trend in the brightness, using a polynomial fit that describes the nightly trend but no shorter features. The x-axis gives the deviation of points from the mean magnitude after removal of the trend. The maximum rms error of a single measurement in the U band calculated from the observations is 0.011 mag for the check star and 0.021 mag for T Tau (in the U band the check star was roughly three times brighter than T Tau). We plot these errors with dashed lines in the histograms. The histogram of T Tau suggests that some of the variability is due to real fluctuations beyond statistical errors, while this is not the case for the check star.
We present the combined and fluxed RGS1+2 spectrum in Fig. 5.
Given the significant background radiation level, we prefer to first show the full spectrum
before background subtraction and the background spectrum itself separately. The
potential influence of background features on the spectrum at the source
position can thus be examined in detail. The unsubtracted spectrum also retains maximum
signal-to-noise ratios to reveal weak line features. (We also inspected the spectrum extracted for low-background episodes
only. Apart from a lower S/N, it very closely resembles the spectrum shown here.)
For the subsequent quantitative spectral analysis, we will use background-subtracted RGS spectra
and provide separate illustrations (see below). We mark all spectral
line features that contained (3-4)
10-6 photons cm-2 s-1, identifying the dominant ions contributing to the flux (a few of these lines may not be individually recognized in the spectrum). Dotted vertical lines point to the location of a few
important fainter lines that remained undetected.
The spectrum is rich in lines, but also reveals a strong continuum shortward of 17 Å.
The most unusual aspect of the T Tau X-ray spectrum is its
mixture of features associated with extremely hot plasma, such as the continuum
and several lines of highly-ionized Fe formed above 10 MK,
and lines attributed
to very cool plasma such as those from N VII and O VII.
The ratio between the fluxes in the O VII He-like line triplet at 21.6-22.1 Å and the O VIII Ly line at 18.97 Å is near unity. This is unusual
for magnetically active main-sequence stars and WTTS, while other CTTS have shown similar
ratios, as noted in a survey
of TTS X-ray grating spectra presented by Telleschi et al. (2007b). Among main-sequence solar analogs,
only inactive, old stars reveal such low flux ratios (Telleschi et al. 2005). We also note
the unusually clear appearance of two He
lines of O VII, at 17.7 Å and at 18.6 Å and the presence of the Ly
line of N VII
at 24.7 Å. On the other hand, several features in the unsubtracted spectrum
coincident with the He
lines of N VI at 28.5-29.5 Å and the He
lines of N VI at 23.3 Å and at 23.8 Å are not significant
because the background spectrum contains similar features at the respective wavelengths.
We emphasize that the identified features of O VII do not suffer from irregular background.
We have performed spectral fits in XSPEC (Arnaud 1996), using the vapec collisional
ionization equilibrium thermal line+continuum model combined with a common
photoelectric absorption component, parameterized by the neutral hydrogen column density .
We use background-subtracted spectra for this analysis.
The source spectrum is sufficiently rich in spectral features that
a determination of several element abundances is warranted. We note that alternative
methods such as reconstruction of emission measure distributions from extracted line
fluxes could in principle be applied; however, most of the important line features
required for such an analysis are heavily blended, and from no element (including Fe)
do we see lines from more than two ionization stages at sufficient signal-to-noise ratios.
We therefore prefer to use template spectra in XSPEC that take into account all
unresolved blends and faint lines present in the code, while the principal
thermal components are determined by the strongest lines and the continuum. A comparison between
these two approaches for similar but brighter X-ray spectra from solar analogs has
shown little difference (Telleschi et al. 2005).
We simultaneously fitted spectra of RGS1 and RGS2 between 7 Å and 25 Å, complemented with a portion of the MOS1 spectrum in the range of 1.5-9.35 Å.
We use only one MOS spectrum confined to the harder range because we prefer to
give full weight to the high-resolution RGS spectra wherever possible. The MOS spectrum complements information on the Mg, Si, S, and Fe abundances and constrains the
hottest plasma components. All spectra were rebinned so that they contained at least 15 counts per
bin before background subtraction.
Three thermal components were required for a good fit, defined by their temperatures T1,2,3 and
the associated emission measures EM1,2,3.
The results are shown in Table 2, together with 68% and 90% confidence ranges.
The reduced
is excellent (1.09 for 382 degrees of freedom = d.o.f.). We also list
the "average temperature''
,
which has been defined as the logarithmic
average of the three temperatures, weighted with the respective emission measures;
its 68% range is an approximation calculated using Gaussian error propagation.
This temperature is considerably lower than
MK found in the XEST survey
from EPIC data (Güdel et al. 2007a). This is because the EPIC spectral fit does not
recognize the softest emission for which the cool component at 1.8 MK is responsible.
If we calculate
for the two hotter components, we obtain
MK,
in agreement with the EPIC survey value.
The fits are illustrated in Fig. 6 for two spectral ranges
of the RGS detectors. The RGS data are shown at a higher resolution (fewer counts per bin) than actually used for the fit, in order to better illustrate the individual line features.
Table 2: Results from the 3T model fit.
![]() |
Figure 6: Spectral fit to the background-subtracted RGS1 (black) and RGS2 (red) data in the 11-18 Å ( top) and 18-23 Å ( bottom) ranges. Note the missing data in the RGS1 spectrum below 13.8 Å and in the RGS2 spectrum above 20 Å. These gaps are due to chip failures earlier in the mission. The lower panels in each plot show the residuals in terms of sigmas. |
Open with DEXTER |
The unusual mixture of very hot (30 MK) and very cool (
1.8 MK) plasma
is confirmed. In the intrinsic, unabsorbed model spectrum, the O VII resonance
line is in fact the strongest line in the entire X-ray spectrum, exceeding even the O VIII Ly
line by a factor of 2.0. Such ratios are otherwise typical only for extremely inactive stars with cool coronae, such as Procyon (Raassen et al. 2002).
The amount of cool emission measure is similar to the total
emission measure of the hotter two components although the precise ratio strongly depends on ;
the latter sensitively affects the soft part of the spectrum to which the cool plasma mostly
contributes. The value of
1021 cm-2
is somewhat higher than found from EPIC data
alone (Güdel et al. 2007a),
1021 cm-2 (from a 2-T analysis). Part of this discrepancy could be induced by slight, systematic
calibration inaccuracies at the long-wavelength end of the RGS. Further, for spectra
of limited S/N, an acceptable spectral fit could overestimate the cool EM if
is slightly too high. However, the basic result is evident in the flux ratios of the
O VII and O VIII lines at similar wavelengths.
We find a ratio between the observed fluxes in O VII and O VIII of 1.06
0.29. From the model, we infer that about half of the O VIII Ly
line flux is due
to the two hotter components. The other half must be due to a cooler component and therefore
amounts to about 50% of the total O VII flux at 21.6-22.1 Å. Keeping the hotter
two components fixed in the model but enforcing
1021 cm-2 as in
the XEST survey (Güdel et al. 2007a, noting that
of this magnitude principally influences temperature and emission measure of the coolest component), we calculated conditions for T1 and EM1 to produce the observed flux ratio, given the contributions from the hotter
components. There is a family of solutions for (T1, EM1)
below T1 = 2.1 MK, with the coolest emission measure still being 17% larger than EM2,
and this solution is obtained for
T1 = 1.59 MK. We conclude that even for the lower
,
the coolest emission measure is comparable to the intermediate and hotter EMs.
The element abundances (given with respect to solar photospheric abundances of
Anders & Grevesse 1989, except for Fe for which we used the value of Grevesse & Sauval 1999 as a basis)
are common among TTS (Telleschi et al. 2007b). We note
a relatively large abundance ratio for Ne/Fe (=). The abundance of N is very poorly
determined, given the very low S/N of any potentially detected line. The usual trend toward an inverse first ionization effect in magnetically active stars (Brinkman et al. 2001; Güdel et al. 2001)
is marginally recovered.
We used the O VII He-like line triplet at 21.6-22.1 Å to infer electron
densities in the coolest plasma component. In high-density environments, the forbidden transition
(
)
at 22.1 Å can become suppressed by electron collisions
in favor of the intercombination transitions (
)
at 21.8 Å (Gabriel & Jordan 1969; Blumenthal et al. 1972, see Telleschi et al. 2007b
for applications to high-resolution spectra of other pre-main sequence stars). The
density-sensitive range approximately covers
.
The limited S/N makes individual
extraction of the three line fluxes difficult. The i line is not significantly
detected, while the r and f lines are present (at the 2.3 and 3.6
level, respectively), suggesting low densities. Formally, the unabsorbed fluxes at the
distance of the Earth are
10-4 photons cm-2 s-1in the (r, i, f) lines, resulting in a best-fit f/i ratio of 4.0. However, we note that there is
an additional constraint on the relative line fluxes. The so-called G ratio, (r+i)/f, where r, i, and f denote the fluxes in the resonance, intercombination, and forbidden lines, is
determined by the thermal structure of the plasma. We therefore fixed our best-fit
thermal model and made use of the density calculations implemented in the vmekal
code in XSPEC to infer the densities. We only fitted
and, for optimum adjustment of
the line fluxes, the emission measure of the coolest component that predominantly
forms the lines. The high f/i ratio requires low densities. The best-fit values are below the density-sensitive
range (in agreement with the above f/i estimate) so that we quote only the
68% and 90% upper limits in Table 2. These are approximately 2.6
1010 cm-3 (68% confidence)
and
8
1010 cm-3 (90% confidence) using Gehrels weights on our spectral data
.
The X-ray luminosities, ,
in Table 2 are given for the
0.3-10 keV range (adopted for the
XEST survey, Güdel et al. 2007a) and also for the 0.1-10 keV range to include
important contributions from the coolest component. The latter indeed doubles
compared to the EPIC analysis in Güdel et al. (2007a). Not considering
the contributions from the cool plasma, we find
1031 erg s-1for the 0.3-10 keV range, nearly identical to the XEST survey value of
0.94
1031 erg s-1 based on a two-component spectral fit to the EPIC PN data.
We now discuss the X-ray properties of T Tau N as found in this paper in the context of X-ray features that have previously been associated with accreting pre-main sequence stars, and with properties seen in the U-band and optical light curves. We address, in turn, the following characteristics: i) flux ratios in He-like triplets, indicative of electron densities; ii) cool plasma; iii) abundance anomalies; and iv) optical and ultraviolet variability. We refer to Telleschi et al. (2007b) and Telleschi et al. (2007c) for further detailed discussions of high-resolution X-ray spectra from pre-main sequence stars.
He-like line flux ratios and electron densities. High electron densities in cool
plasma could be a consequence of accretion shocks. The post-shock temperature
in the accretion column near the stellar surface is
,
where the velocity vis about 0.9 times the free-fall velocity
(Calvet & Gullbring 1998),
is the proton mass, k is the Boltzmann constant, and the mean molecular weight
for a fully ionized gas. For T Tau, we use
and
as listed in Güdel et al. (2007a). Then,
km s-1 and
106 K. The soft excess is thus compatible with expected post-shock temperatures of the free-falling gas.
The accretion rate of T Tau is
yr-1 (Calvet et al. 2004; White & Ghez 2001).
For the accretion luminosity, we use the expression given
by Telleschi et al. (2007c) but for a somewhat larger corotation radius:
,
where
erg s-1),
,
yr-1 and
.
We thus obtain
,
i.e., at least two orders of magnitude more than
required to produce the radiative output in the softest X-rays. Some energy will
be conducted to cooler layers, and the majority of the X-rays may be entirely absorbed
(Lamzin 1999).
Lastly, we can estimate the shock density, using the strong-shock condition n2 = 4n1, where n1 and n2 are the pre-shock and post-shock densities, the former relating to the mass accretion rate and the accreting surface fraction F (filling factor - assuming that all flows produce similar shocks). Telleschi et al. (2007b) gave the expression
We also inspected the Ne IX density-sensitive triplet formed at somewhat higher temperatures (4 MK).
The signal-to-noise ratio is low and the Ne IX triplet is strongly blended
with several lines of highly ionized Fe. Nevertheless, there is a clear indication of a strong f line in that triplet as well, again pointing to low densities (the Ne triplet is density-sensitive in the range of approximately
cm-3).
In summary, then, standard accretion shock physics does not seem to be compatible with a shock-heating scenario for the soft spectral component forming the bulk of the O VII (and Ne IX) lines in T Tau.
We emphasize that T Tau is now the second accreting pre-main sequence star that shows a high f/i flux ratio, after similar reports on the Herbig star AB Aur (Telleschi et al. 2007c). This is to be compared with
three CTTS previously reported to show high densities,
cm-3, namely TW Hya (Kastner et al. 2002), BP Tau (Schmitt et al. 2005), and V4046 Sgr (Günther et al. 2006).
Abundance anomalies. High abundances of the volatile elements N and Ne in TW Hya with respect to the Fe abundance have been suggested to be due to stronger condensation of Fe and some other elements to dust grains in the circumstellar disk (Stelzer & Schmitt 2004). If the star preferentially accretes from the gas phase, then an elevated Ne/Fe ratio would be measured in the coronal or shock-heated plasma.
Anomalously high Ne abundances have been found in the X-ray source of TW Hya (Kastner et al. 2002; Stelzer & Schmitt 2004), but also in non-accreting evolved stars such as members of the RS CVn class (Ne/Fe abundance ratios up to 13.4, Audard et al. 2003), and in pre-main sequence stars thought to be non-accreting (weak-line) TTS, such as TWA 5 (Ne/Fe
10; Argiroffi et al. 2005) and HD 98800 (Ne/Fe
5; Kastner et al. 2004b). Telleschi et al. (2007b) have systematically analyzed RGS spectra of
pre-main sequence stars and compared Fe/Ne and Fe/O abundance ratios for them. They found that i) both ratios tend to be higher for young G-type stars than for K- or M-type stars; this trend still holds for near-zero-age main-sequence
stars; and ii) there is no difference between CTTS and WTTS with regard to this trend.
Figure 7 shows the Fe/Ne ratio derived from high-resolution spectroscopy as a function of spectral type for TTS and ZAMS stars (for references, see Telleschi et al. 2007b). The trend for an increasing ratio toward earlier spectral types is clearly seen, independently of the evolutionary status (CTTS, WTTS, ZAMS). The ratio for T Tau fits in well for a K0 star. The same holds true for the Fe/O ratio (not shown; see Telleschi et al. 2007b). We conclude that the abundances observed in T Tau follow a trend that holds for a larger sample of young stars, and this trend - as judged from a still small sample of stars - does not seem to be related to accretion properties.
![]() |
Figure 7: The Fe/Ne abundance ratio in the X-ray sources of various pre-main sequence stars and ZAMS stars as a function of spectral class. The abundance measurements for stars other than T Tau were taken from the previous literature and from Telleschi et al. (2007b) (see that paper for a comprehensive description and references; values were converted, if necessary, to refer to the solar photospheric abundances of Anders & Grevesse 1989 and, for Fe, Grevesse & Sauval 1999). The value determined in the present work for T Tau is shown in bold face. Symbols mark different types of stars: squares: ZAMS stars; triangles: WTTS; diamonds: CTTS. |
Open with DEXTER |
Soft excess. The most outstanding property of the old CTTS TW Hya is its extremely soft X-ray emission (Kastner et al. 2002; Stelzer & Schmitt 2004).
However, no such X-ray sources have been identified in large samples of CTTS in Orion (Preibisch et al. 2005) or in Taurus
(Güdel et al. 2007a). On the other hand, Telleschi et al. (2007b) have suggested that, apart from a dominant hot plasma, there
is a soft excess in accreting pre-main sequence stars, best revealed by an unusually large ratio between the (summed) photon fluxes
in the O VII He-like triplet and the O VIII Ly line, found to be around unity for CTTS but typically less than 0.5 for WTTS. The intrinsic ratio between the luminosities of these lines is of course
altered by photoelectric absorption, i.e., the O VII/O VIII ratio tends to decrease for increasing
.
The
value for T Tau compares with the largest in the Telleschi et al. (2007b) sample, yet the
observed flux ratio is also among the largest. As reported above, we found a ratio of 1.06
0.29. In the unabsorbed model spectrum,
the O VII r line is even a factor of two stronger than the O VIII Ly
line, and in fact is the strongest line in the
X-ray spectrum.
T Tau thus confirms the conjecture that CTTS reveal a soft excess, despite the dominance of an extremely hard component in the overall X-ray spectrum. T Tau is the most extreme case for a soft excess in any TTS studied by high-resolution spectroscopy, except for TW Hya in which the soft component entirely dominates.
Optical and ultraviolet variability. Accretion shocks induce excess emission at short optical and ultraviolet wavelengths.
Variability in these bands seen in accreting TTS has therefore been attributed to short-term variations in the mass accretion,
or simply inhomogeneities at certain locations in the streams, rather than to flare-like energy release events (Gullbring et al. 1996).
Variability amplitudes are typically modest (a few tenths of a magnitude) and last for 0.6 h to several hours. The gas heats
only moderately, again in contrast to chromospheric gas in flares (Gullbring et al. 1996). The induced variability may be explained
by accretion rate changes of the order of 10%, which will then require a scale length of the inhomogeneities of 5
1010 cm for
the CTTS BP Tau and, combined with the free fall velocity, correctly explains the variability time scales (see Gullbring et al. 1996, for
further details). It is possible that these inhomogeneities reflect in variations of the softest X-ray flux if X-rays are
produced in the accretion shock, but a dedicated monitoring study of BP Tau with ROSAT and ground-based UBVRI photometry revealed
an absence of correlated events; X-ray variability may thus not be related to accretion-induced optical and UV variability (Gullbring et al. 1997).
New evidence was provided by the comprehensive XEST survey: a comparison of the X-ray light curves of many TTS with simultaneous light curves obtained in the U band or a shorter-wavelength ultraviolet band showed little
correlated behavior, except for a few examples of flares that were observed in both bands (Audard et al. 2007).
There is mixed evidence in our data. First, both the XMM-Newton ultraviolet data and the ground-based U-band data show clear evidence
for both day-to-day variations and short-term variability. The former shows time scales of at least half a day, compatible with modulation
of accretion hot spots due to rotation. The short-term events (e.g., at 6
104 s in Fig. 2, or at 37.35 d in the bottom panel
of Fig. 3) are similar to those reported by Gullbring et al. (1996). Conclusions are ambiguous, however: During the XMM-Newton
observations, there is a long-term trend in X-ray hardness that appears to correlate with the UV flux; in particular, there is a decay
from the beginning of the observation that may be related to flare emission. The harder X-rays seen in our spectra cannot be produced
by accretion shocks as the infall velocity is too low (see above). Some of the variations of the UV flux may therefore be due
to an excess in chromospheric and transition region heating (producing excess UV emission) in concert with coronal heating (producing
harder X-rays). Also, the short event at 6
104 s in the UV appears
to be reflected in the X-ray hardness curve. The temporal relation is typical for flares: The ultraviolet excess, a signature of
impulsive heating of the chromosphere, precedes the X-ray signal, a consequence of filling of magnetic loops with hot plasma.
Thus, overall a significant contribution to UV variability may come from coronal energy release events, although a clear distinction
is difficult. This view is supported by the large survey presented by Audard et al. (2007) that shows several obvious X-ray flares with
correlated U-band or ultraviolet features.
Modulation effects due to rotation are also possible. The UBVR observations in Fig. 3 show some evidence (e.g., within
the middle or the lower panels, or between the upper and middle panels) for a correlation between B or V magnitude and the U-B color, indicating that the source is bluer when brighter and redder when fainter. This is expected for both
heated regions and cool spots. Heated regions could be due to either
accretion hot spots or magnetically active chromospheric regions. However, the inclination angle of T Tau is small,
i = 19-23 deg (Eislöffel & Mundt 1998; Herbst et al. 1997). Rotationally modulated surface features must therefore be located at latitudes
of approximately 20 deg around the equator. The standard magnetic accretion scenario (e.g., Calvet & Gullbring 1998) predicts that accretion
hot spots form at high latitudes. It is therefore more likely that, at least during our UBVR observations, the modulation
was due to active regions (unless the accretion also occurs along more complicated magnetic field structures toward
equatorial regions). The indication of (at least) two minima and two maxima in Fig. 3 could be due to two preferred active longitudes producing chromospherically active regions near the equator (e.g., Berdyugina et al. 1998).
Our discussion thus points to the following conclusions: i) if mass accretion streams induce X-rays in shocks near the surface, then densities should be much higher than inferred from the observed density-sensitive line flux ratios. It is possible that some of the observed X-ray emission originates in accretion shocks, but only a fraction of the accretion stream, subject to very low densities, would be responsible; ii) abundance anomalies do not seem to reflect accretion properties, and T Tau is no exception; iii) T Tau is, however, the most outstanding example, apart from TW Hya, of an accreting TTS with a soft excess in the X-ray spectrum. Such excess is not found in non-accreting young stars. We thus suggest that the X-ray production of T Tau is influenced by the accretion process although the X-rays may not form in the bulk of the accretion footpoints.
Shocks could form in outflows and jets that are related to the accretion process (see Güdel et al. 2005, for the example of DG Tau A). In that case, low temperatures and low densities would be expected. Since T Tau does drive outflows and is surrounded by Herbig-Haro objects (see references in Sect. 1), this possibility cannot be excluded, but explicit evidence is lacking. The high-resolution Chandra HRC image (Fig. 1) does not show any evidence for emission that might be related to a jet, in contrast to what has been found for DG Tau A (Güdel et al. 2005). The entire X-ray spectrum can be explained by a multi-thermal plasma subject to an overall photoelectric absorption component, again in contrast to the spectra of proposed jet-related X-ray sources (Güdel et al. 2007b).
An alternative location for the interaction between the plasma heating process and accretion is the corona or magnetosphere itself. Accreting material is guided along magnetic field lines and can increase the density in these regions considerably (similar to the estimates presented for the footpoints of the accretion flows, see above). Depending on the stream density and the strength of the magnetic fields near the accretion streams, the magnetic geometry may be considerably shaped by the mass stream itself; magnetic field lines be straightened by the flow, perhaps suppressing coronal magnetic energy release so that coronal heating may have little or no effect on infalling gas.
If, however, a small fraction of the mass loads onto magnetic fields of active coronal regions,
it will lower the overall temperature and slightly increase
the density there; cooling thus also becomes more efficient. A fraction of the coronal
volume affected by a small amount of the infalling material could thus produce a cool coronal component, i.e., an observable soft excess.
A suggestion along those lines was made by Preibisch et al. (2005) who proposed that the lower
ratio of CTTS compared to WTTS in Orion is due to accretion flows completely suppressing coronal heating in the respective magnetospheric regions.
We estimate what fraction of the accreting material would have to be heated to T = 1.8 MK as observed in the softest component,
assuming that the entire soft component is emitted by such material (rather than by plasma evaporating from the stellar surface).
We provide this estimate for the limiting case that any parcel of inflowing material that becomes subject to coronal heating is heated to T at once, after which
it is permanently cooling by radiation, without further heating. The energy rate to heat accreting, ionized gas to T = 1.8 MK is
for hydrogen gas, where
is the rate of new, cool electrons flowing
into the heating region; there is an equal rate of protons.
Adopting
yr-1 for the total accretion rate, we find
1042z s-1, where z is the fraction of the infalling gas stream that is heated. Therefore
1032z erg s-1. We assume that all energy is radiatively
lost by the plasma at 1.8 MK. Modeling the energy losses of the cool plasma component in XSPEC, we find that only 17% of the
radiative power is emitted above 0.3 keV, or 28% above 0.1 keV. Because the cool component emits 5.2
1030 erg s-1
into the 0.3-10 keV band (Sect. 6), the total radiative energy loss rate of this plasma is
3
1031 erg s-1. We conclude that
z = 1.9-3.8% of the accretion flow
would need to be heated to low coronal temperatures. This fraction will be smaller if the gas is
continually heated while radiating.
We conclude, then, that T Tau implies an important role of accretion in the production of high-energy emission. We suggest that part of the accreting material adds a cool coronal component by being heated to modest temperatures while flowing into active regions. We are not able to predict why the X-ray excess is found at temperatures of about 2 MK, but we suggest that the additional loading of magnetic fields by cool, accreting mass will suppress heating to the otherwise common coronal temperatures of >10 MK.
Acknowledgements
We thank Svetlana Berdyugina for helpful comments on photometry and active longitudes. We warmly acknowledge financial support by the International Space Science Institute (ISSI) in Bern to the XMM-Newton XEST team. X-ray astronomy research at PSI has been supported by the Swiss National Science Foundation (grants 20-66875.01 and 20-109255/1). Part of 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 the USA (NASA). M.A. acknowledges support by NASA grants NNG05GQ13G and NNG05GF92G for the XMM-Newton project. In addition, he acknowledges support from a Swiss National Science Foundation Professorship (PP002-110504). Support for this work was provided by the National Aeronautics and Space Administration through Chandra Award Number SAO GO3-4004X issued by the Chandra X-Ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics and Space Administration under contract NAS8-03060.