Free Access
Issue
A&A
Volume 505, Number 2, October II 2009
Page(s) 755 - 762
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/200811374
Published online 11 August 2009

The X-ray spectrum of the classical TTauri star TW Hydrae observed by LETGS aboard Chandra

A. J. J. Raassen1,2

1 - SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands
2 - Astronomical Institute ``Anton Pannekoek'', Kruislaan 403, 1098 SJ Amsterdam, The Netherlands

Received 18 November 2008 / Accepted 27 July 2009

Abstract
Aims. I present the analysis of the X-ray spectrum of the classical TTauri star TW Hydrae, and determine emission measure distributions, elemental abundances, electron temperatures, densities, and line fluxes.
Methods. TW Hya was observed in the wavelength range 5-175 Å by CHANDRA-LETGS. I performed a multi-temperature fit, based on a collisional ionisation equilibrium (CIE) model for optically thin plasmas. I measured line fluxes by fitting a Gaussian line profile, folded through the response matrix, to the observed line features and I applied these values to establish plasma properties.
Results. No line features were established in the spectrum above 40 Å. Due to the absence of line features above 40 Å $N_{\rm H}$ is calculated to be somewhat higher than assumed so far. The observed ``coronal'' temperature covers a range from about 1.5 to 20 MK. The abundance pattern shows a very strong enhancement of Ne and some enhancement of N. From forbidden-intercombination line ratios of O VII and of Ne IX I obtain a value for the electron density of $n_{\rm e}\simeq 10^{12}~{\rm cm}^{-3}$. This is higher than the value for electron densities in coronal plasmas of older late type main-sequence stars. A blackbody component was added to describe the shape of the continuum around 15 Å better. The emitted X-ray luminosity is in the order of 1030 erg/s.

Key words: techniques: spectroscopic - stars: individual: TW Hya - stars: coronae

1 Introduction

TTauri stars (TTS) are young low-mass pre-main sequence stars, especially present in star-forming regions. They can be divided into two types: classical TTauri stars (cTTS) and weak-line TTauri stars (wTTS). Different from the cTTS, the weak-line TTauri stars do not show strong H$_\alpha$ emission. Classical TTauri stars are surrounded by a disc of gas and dust, while the weak-line TTauri stars have no surrounding disc. Both types of TTauri stars are X-ray emitters. The X-rays were observed by the Einstein Satellite (e.g. Feigelson & DeCampli 1981) and by ROSAT (e.g. Feigelson et al. 1993). The X-rays from weak-line TTauri stars are commonly interpreted as being produced by a scaled up coronal structure. The same is true for classical TTauri stars. However, in cTTS additional mechanisms might be present e.g., jets, winds, or accretion from the disc to the stellar surface. In the last case it is expected that the material funnels from the disc to the stellar surface along magnetic field lines, resulting in hot X-ray emitting spots. TW Hya belongs to the group of classical TTauri stars. It is spectral type K8V (SIMBAD) and nearby, at a distance of 56.4 parsec (Perryman et al. 1997a, b), recently corrected to 53.7 parsec (van Leeuwen 2007). It is surrounded by a circumstellar disc. A wide range of periods has been given over the years (e.g., Herbst & Koret 1988; Mekkaden 1998; Rucinski et al. 2008). Recently, Setiawan et al. (2008) have concluded that a Jupiter-like planet is present close to the star, inside the disc and with a period of 3.56 days. However, Huélamo et al. (2008) state that a cool spot on the stellar surface is more likely to explain the observations.

Since the era of high-resolution X-ray spectroscopy started after the launches of Chandra and XMM-Newton, many papers on the high-resolution X-ray spectrum of TW Hya have appeared. The RGS aboard XMM-Newton was used by Stelzer & Schmitt (2004) and Robrade & Schmitt (2006), while HETGS aboard Chandra was applied by Kastner et al. (2002) and Ness & Schmitt (2005). Robrade & Schmitt (2006) and Stelzer & Schmitt (2004) made a fit to the spectrum with three temperature bins, while Stelzer & Schmitt (2004), Kastner et al. (2002), and Ness & Schmitt (2005) investigated individual lines and line fluxes. Güdel & Telleschi (2007) studied the O VII/O VIII line ratio in an ensemble of wTTS and cTTS.

For RGS and HETGS the wavelength ranges are at the higher end limited to about 35 Å. Applying LETGS, however, I may extend the wavelength range to far higher wavelengths and check the effects of the models above the RGS and HETGS limits.

The strategy in this paper is to perform a multi-temperature fit to the total spectrum (Sect. 3). In this fitting procedure firstly two different models have been applied (CIE and NEI) and the optimal number of T-components has been determined. Results of the global fitting are compared with results from the analysis based on individual lines, described in Sect. 4. In Sect. 5 I discuss the agreements and differences in comparison with the papers mentioned above (Robrade & Schmitt 2006; Stelzer & Schmitt 2004; Kastner et al. 2002; and Ness & Schmitt 2005).

2 Observations

The spectrum of TW Hya was obtained on August 22, 2006 during a 151.64 ks observation using HRC-S in combination with LETG aboard CHANDRA. The log of the observations is shown in Table 1. The data were pipelined with CIAO3.4 with CALDB version 3.4.2 and prepared to be implemented in SPEX by the programs extract and crele developed at SRON.

 \begin{figure}
\par\includegraphics[angle=-90,width=8.8cm,clip]{11374f1.ps}
\end{figure} Figure 1:

The zeroth order lightcurve of the LETGS observation of TW Hya. The time bin size is 1800 s.

Open with DEXTER

Table 1:   Observation log of the data of TW Hya.

In the lightcurve based on the total detector area some time bins showed saturation, i.e., they registered over 184 counts/s. Although this saturation will only sparely influence the spectral counts, the saturated time bins were excluded from the further procedure. Figure 1 shows the zeroth order lightcurve. The time bin size is 1800 s. The average count rate is 0.07 counts/s. This lightcurve shows variability up to about 50% over the total exposure time. However, no specific time intervals could be selected to obtain a reliable spectrum of a ``flare'' state. No periodicity was established during the observation.

The observed LETGS spectrum covers the range from 1 to 176 Å. The final spectrum is the result of the sum of the spectrum of the +1 order and the -1 order. The background, which is high at longer wavelengths, was subtracted. No significant line features (apart from higher order lines) are established above 40 Å.

3 Multi-temperature fitting

3.1 A 10-T CIE and NEI model as a starter

I started the analysis of the LETGS-spectrum of TW Hya by means of two different spectral models: a model for optically thin plasma in collisional ionisation equilibrium (CIE) and a model for a non-equilibrium plasma (NEI). This approach was chosen to investigate the possibility of establishing different mechanisms for hot X-ray radiation i.e., a coronal part as for the wTTS and an accretion part from in-falling disc material. The two different models resulted in a solution with three significant temperature components only and with comparable $\chi^2$ values. However, the parameter U (= $n_{\rm e} \times t)$ in the NEI model obtained a value that high that equilibrium is reached and the NEI model mimics the CIE model.

From this point on I use a 3-T model for optically thin plasma in collisional ionisation equilibrium (CIE) as implemented in SPEX version 2.0[*] (Kaastra et al. 1996a) in combination with MEKAL (Mewe et al. 1995; Kaastra et al. 1996b). The ionisation equilibrium is based on calculations by Arnaud & Rothenflug (1985).

3.2 A 3-T model

A 3-T model was fitted to the data over the total wavelength range from 5 to 175 Å. The T-parameters, emission measures and abundances were freely varied. Two temperature components, describing most of the line features, are established at 0.127 keV and 0.264 keV. They describe very well the observed line structure above 10 Å. Based on the Mg XI, Mg XII and Si XIII, Si XIV line features below 10 Å a third temperature component at 1.45 keV is determined. A fourth T-component of 0.65 keV appeared to model the continuum. These components are in agreement with the temperatures found by Robrade & Schmitt (2006) and Stelzer & Schmitt (2004) in a different way and based on observations with RGS on board of XMM-Newton.

However, fitting the spectrum with these four CIE components results in some deviations, i.e., the model continuum is a bit lower than the data around 15 Å and a bit higher below 10 Å, while weak line features below 10 Å are not recognised as lines by the model. I notice the same deviations when applying the parameter set given by Robrade & Schmitt (2006). To solve this deficiency in the shape of the underlying continuum around 15 Å, a blackbody component was introduced, replacing the T = 0.65 keV component. The blackbody component has a well determined temperature of 0.187 keV, just in between the two cool temperature components. This component is strongly related to the 0.60 keV temperature component of Robrade & Schmitt (2006), but doesn't affect the other T-components. From Table 2 of Günther et al. (2007) we notice that the CIE component at 0.7 keV doesn't improve their fit very significantly. The blackbody component improved the shape of the continuum and it lowered the $\chi^2$ of the fit significantly. Due to the lower continuum below 10 Å the weak line features appear more pronounced. It results in a higher value of the abundance of Si, compared to the value given by Robrade & Schmitt (2006). It is possible that my blackbody component takes into account many weak lines, that are not included in the models and produce a pseudo-continuum.

The spectrum and the model are shown in Fig. 2. The grey line (green in the electronic version) is the model fit to the spectrum.

The parameter values of the fit are collected in Table 2. This table shows the interstellar column density $N_{\rm H}$, the determined temperature values, the emission measures and the emitted luminosities of the different temperature bins, the electron densities, and the abundances. The emission measure is defined as $EM=n_{\rm e}n_{\rm H}V$, in which $N_{\rm H} = 0.85~\times~n_{\rm e}$. Except for the density, broadening, and w-shift, the parameters were all fitted simultaneously.

 \begin{figure}
\par\psfig{figure=11374f2a.ps,angle=-90,height=5.5truecm,width=8....
...igure=11374f2b.ps,angle=-90,height=5.5truecm,width=8.8truecm,clip=}
\end{figure} Figure 2:

The LETGS spectrum of TW Hya together with the fitted CIE model, grey line (green in the electronic version). The top panel shows the total spectrum and the bottom panel shows the area between 0 and 40 Å.

Open with DEXTER

Table 2:   Parameters of the 3-T fit to the spectrum of TW Hya.

The interstellar column density $N_{\rm H}$ was free to vary. Only one $N_{\rm H}$ parameter was used for all four components. However, it is not necessarily true that all plasma undergoes the same absorption. The continuum at higher wavelength is weak and not very reliable due to the high value of the subtracted background. However, line features or the absence of line features are indicative for the interstellar column density. The absence of the Ne VIII lines around 88 and 98 Å indicates a value for $N_{\rm H}$  $\geq 2~\times 10^{20}~{\rm cm}^{-2}$. My fitted value ( $N_{\rm H}= 5.2~\times 10^{20}~{\rm cm}^{-2}$) is higher than that value and also slightly higher than the values used by Robrade & Schmitt (2006) and Stelzer & Schmitt (2004). This higher value is necessary to suppress Si XII line features near 44 Å that appear in the model, but are absent in the observed spectrum. These lines are outside the wavelength range observed by Robrade & Schmitt (2006) and Stelzer & Schmitt (2004).

A slight shift in $\lambda$ ($\approx$1.2e-4) was applied to optimise the positions of the model lines with respect to the data. This shift corresponds to about 35 km s-1. The value of $v_{\rm rad}$ of TW Hya is 12.5 km s-1 (Reid 2003), while the satellite and earth orbital motion result in about -12 km s-1 at time of observation. My fitted value is higher than the result of these geometrical motions. An instrumental wavelength calibration bias is likely on forehand and cannot be ruled out.

A line broadening parameter, which corresponds to a velocity broadening of 287 km s-1, was added to improve the fit. This parameter is added to the instrumental line spread function, which is taken into account by means of folding the model with the response matrix. Ness & Schmitt (2005) obtain the best fit results by introducing an ${\it FWHM}$ of 0.016 Å for all lines of their HETGS observation. This is equivalent to about 200 km s-1. They do not mention whether this broadening is included with the instrumental broadening or not.

The densities given in Table 2 are obtained by fitting to the line triplets of He-like ions (see also Sect. 4.3) in a narrow selected wavelength range around 22 Å (O VII) and around 13.5 Å (Ne IX) only. During this fit procedure the oxygen and neon abundances and electron density ($n_{\rm e}$) were free to vary, while all other parameters were fixed at values obtained in the overall fitting procedure. The abundance served as normalisation parameter. The obtained electron densities are $n_{\rm e} = 6.3(2.6)~\times~10^{11}$ cm-3 based on the O VII line triplet and $n_{\rm e} = 2.5(0.4)~\times~10^{12}$ cm-3 based on Ne IX.

The emitted luminosities are calculated over the energy range from 0.3 to 10 keV. They are obtained for all given temperature bins. The uncertainty in the luminosities is assumed to be the same as the relative uncertainties of the emission measures. Robrade & Schmitt (2006) and Stelzer & Schmitt (2004) give luminosity values in slightly different energy ranges. In their energy ranges, 0.2-10 keV and 0.45-2.25 keV, respectively, they obtain $2.0 \times 10^{30}$ erg/s and $1.4 \times 10^{30}$ erg/s. The luminosities in these ranges obtained in the model fit of this study (corrected for the distance they use) are 2.05 $\times10^{30}$ erg/s and $1.47 \times 10^{30}$ erg/s, respectively. They are in very good agreement with those of Robrade & Schmitt (2006) and Stelzer & Schmitt (2004). It shows that the luminosity is a very robust quantity. About 55% of the emitted luminosity is due to line radiation, while the resulting 45% is radiated in the continuum or a pseudo-continuum produced by many weak indistinguishable lines not present in the models. These percentages are comparable with values (60% versus 40%) obtained using the parameter set of Robrade & Schmitt (2004).

The abundances are coupled for the three temperature components. This is done to stabilise the fit procedure and is not based on physical grounds. The abundances are relative to solar photospheric values (Anders & Grevesse 1989), except for iron. For Fe the value of Grevesse & Sauval (1999) has been applied. It differs by a factor 1.45. For the obtained abundances the over-abundance of neon (a factor of ten higher than oxygen) is most remarkable, while the abundance of nitrogen is a factor 3 higher in comparison with oxygen. The other abundances are close to solar. Iron seems to be depleted. My absolute abundance values are far higher than values given by Stelzer & Schmitt (2004). This can be explained by the fact that the abundances are strongly anti-correlated with the emission measures. Therefore in many papers an abundance ratio related to the abundance of an arbitrarily chosen, generally dominantly present element is given. The abundance pattern, i.e., abundances relative to oxygen (A/O), agrees very well with the results of Stelzer & Schmitt (2004). My carbon abundance is somewhat higher as a result of my higher $N_{\rm H}$ value. The quantity $A \times EM$ is robust. My values for the emission measures are far lower than those in the paper by Stelzer & Schmitt (2004). However, the ratio between the EM2 and EM1, given by them as 2.35, is the same as in my fit i.e., 2.36.

The errors given in the paper are quite different for different quantities. The T parameters, emitted luminosities, abundance ratios (A/O), the observed line fluxes, and $A \times EM$, given in Sect. 4 are robust. Their statistical 1$\sigma$ errors given in Tables 2-4 are robust. For the abundances and the emission measures, however, the errors given in Table 2 are statistical errors that belong to the given fit-parameters. The absolute values of the abundance and emission measures might change considerably more as long as the abundance ratios (A/O) and $A \times EM$ remain unchanged.

4 Individual line fluxes

The individual observed line photon fluxes of the LETGS spectrum of TW Hya have been measured. Many line features are located between 10 and 25 Å. Unfortunately, the number of line features above 40 Å are scarce. The most obvious features above 40 Å are higher order lines of oxygen and neon. Lines below 10 Å are weak, but definitely present, as they are well positioned in relation to experimental laboratory wavelengths in databases (Kelly 1987, and NIST[*]).

For each observed line a Gaussian profile was folded through the response matrix and fitted to the line feature together with a power-law taking into account for the continuum. For lines used for the temperature and density determination in Sect. 4.3 also the following (maybe more sophisticated) method was applied. This was done to avoid contributions by blends. After the final overall multi-temperature fit (Sect. 3.2), the three wavelength ranges that belong to He-like ions (28.5-29.8 Å for N VI), (21.4-22.4 Å for O VII, and 13.2-14 Å for Ne IX) were selected. Then the ions N VI, O VII, and Ne IX were ignored from the model. Now the model describes the total spectrum except the He-like line triplet of the ignored ions. These lines were measured by folding a Gaussian through the response matrix and by fitting to the lines, as described in Sect. 4, but without the additional power-law, since the continuum is taken into account by the multi-temperature model. During this fit all parameters of the model were fixed, except the Gaussians. This way the positions, observed photon fluxes, and widths ( ${\it FWHM}$) of the lines were established. Thanks to the folding of the Gaussian profile through the response matrix and the Line-Spread-Function (LSF) the given ${\it FWHM}$ are net values, in addition to the instrumental width. The measured wavelengths, widths and line photon fluxes are collected in Table 3 together with the theoretical wavelengths from laboratory measurements (Kelly 1987).

Table 3:   The observed line photon fluxes.

Especially at higher wavelengths the lines are slightly shifted to longer wavelength. Unfortunately the line features do not extent to wavelengths above 40 Å. An underlying instrumental wavelength imperfection, cannot be ruled out. However, all lines given by Stelzer & Schmitt (2004) in their Table 2, based on observations with RGS, are (even more pronounced) shifted to longer wavelengths, albeit they do not report this shift and they do not give error estimates for their wavelength measurements. Moreover, Günther & Schmitt (2008) measure positive line shifts of the same order of magnitude, based on observations with FUSE and HST (see their Tables 4 and 5). It should be kept in mind that these observations concern another wavelength area and lower stages of ionisation. Therefore they might be related to other physical processes. The observed wavelength shift in the X-ray regime seems to be real. If not from physical origin, the observed shift would imply a calibration issue for both instruments (LETGS and RGS).

For a number of lines an ${\it FWHM}$ could be obtained. However, the line widths are in many cases upper limits. The ${\it FWHM}$ of a line can be linked to a velocity (v) by the formula: $FWHM(\AA)/\lambda_{\rm obs} = 2\sqrt{\ln 2}v/c$. Applying this formula to the stronger isolated lines at 13.447 Å, 18.969 Å, 21.602 Å, and 24.781 Å, I obtain values for v of 321(94), 256(57), $\mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$ }}}\hbox{$<$ }}}$358, and 334(130) km s-1, respectively. In Sect. 3.2 I have obtained an overall broadening parameter (v=287 km s-1) to improve the multi-temperature fit; a value comparable with these values. ${\it FWHM}$'s of the order of a few hundreds of km s-1 are also found by Günther & Schmitt based on FUSE observations, but in a wavelength area around 1000 Å and based on cooler ions. Their broadening might originate from other processes and not directly be linked to the broadening in the X-ray spectrum.

The last column shows the observed line photon fluxes as measured at Earth in units 10-6 photons cm-2 s-1. The lines of H-like neon, oxygen, nitrogen, and carbon are strong, as well as He-like lines of neon and oxygen. For these lines the uncertainties of the flux values are of the order of 10%. The iron lines are moderately strong. Especially for neon and oxygen higher Rydberg states are visible up to n = 5. Their relative uncertainties are higher.

4.1 Emission measures from individual line fluxes

Line fluxes (F), emission measure (EM), and abundances (A) are related to each other by means of the formula:

\begin{displaymath}F = (4.1807 \times 10^{-54} \times \lambda \times A \times EM \times 10^{-s}) / {\rm d}p^{2}.
\end{displaymath}

In this formula F is the line flux in units of photons cm-2 s-1, $\lambda$ is the wavelength of the line in Å, A is the abundance relative to solar photospheric values, EM is the emission measure in units cm-3, dp is the distance of the source in parsecs and s is the emissivity factor given by Mewe et al. (1985), applied for the temperature ($T_{\rm m}$) of maximum emissivity. In Table 4 the values of the quantity $A \times EM$ obtained from individual line fluxes are given. The given uncertainties in this quantity are based on the uncertainties of the measured line fluxes and an assumed average uncertainty of 15% in the theoretical transition probabilities of the lines. The latter value might be a bit low for weak lines and a bit too high for strong resonance transitions. Since the forbidden and intercombination lines in He-like ions behave as communicating tubes under influence of the density I have used the sum (f+i) in this table.

Table 4:   Calculated values of $A \times EM$, based on individual line fluxes.

From the table we notice that the value of $A \times EM$ differ up to a factor of 10. For comparable values of $T_{\rm m}$ the values of EM, however, should be the same. The difference in value of $A \times EM$in these cases is caused by deviations from solar abundance values for the involved ions. As an example: the ions O VII and N VII, although appearing in the same temperature range, show a difference of a factor 3 in the values of $A \times EM$. This indicates an overabundance of nitrogen compared to oxygen by a factor of about 3. The same way (based on O VIII and Ne IX), an overabundance of neon compared to oxygen is established by a factor of about 10. Absolute values of the abundances are difficult to determine on the basis of line fluxes, since they are anti-correlated with the emission measure. $A \times EM$ is a robust quantity.

4.2 Comparison with observed line photon fluxes from other papers

In this subsection I compare the observed line photon fluxes obtained in this work with the most recent results from high-resolution X-ray spectroscopy, i.e., Stelzer & Schmitt (2004) using RGS aboard XMM-Newton, Kastner et al. (2002), and Ness & Schmitt (2005), applying HETGS aboard Chandra (see Table 5). The values given by Ness & Schmitt (2005) have been recalculated from energy fluxes to photon fluxes. For quite a number of lines the comparison is satisfactory, taking into account the uncertainties in the obtained fluxes. However, some severe deviations that exceed the given uncertainties are noticed. This is especially true for the comparison between the work of Stelzer & Schmitt (2004) and Kastner et al. (2002).

Table 5:   Comparison of observed line photon fluxes.

The deviations are not limited to weak lines with high uncertainties in their measured flux values, but they appear also for the strong Ne IX lines around 13.5 Å and the strong O VIII line at 19.0 Å (see Fig. 3). These figures show the observed line photon fluxes measured by Stelzer & Schmitt (2004), Kastner et al. (2002), and this work versus wavelength and versus temperature of maximum emissivity of the line. From these figures it is difficult to conclude whether the differences are wavelength dependent, and therefore due to the use of three different instruments, or temperature dependent, and therefore related to variability of emission measures in the source TW Hya. The line ratios of the Ne IX lines, however, are comparable for the three spectra (see also the next section).

 \begin{figure}
\par\psfig{figure=11374f3a.ps,angle=-90,height=5.5truecm,width=8....
...igure=11374f3b.ps,angle=-90,height=5.5truecm,width=8.8truecm,clip=}
\end{figure} Figure 3:

Line fluxes versus wavelength (top) and versus log T (bottom) from: - Stelzer & Schmitt (2004), * Kastner et al. (2002) and $^{\circ }$ This work (black, red, and blue in the electronic version, respectively). Fluxes are given in $10^{-6}~{\rm photons~cm}^{-2}~{\rm s}^{-1}$.

Open with DEXTER

4.3 Temperature and density determination based on individual line flux ratios

The ratios between the line fluxes of the resonance lines from two ionisation stages of the same elements provide valuable information on the temperature regime of the plasma. Apart from that the line triplets in the He-like ions, consisting of a resonance line (r), an intercombination line (i) and a forbidden line (f), are strong diagnostic tools for temperature as well as for density (Gabriel & Jordan 1969).

4.3.1 Temperature

The temperature can be established, based on individual line fluxes. The ratio between the summed fluxes of the line triplets in the He-like ions and the resonance lines of H-like ions as well as the ratio (f+i)/r in the He-like line-triplets are temperature sensitive. As values for the ratios between the He-like line triplets and H-like resonance lines I obtain 0.79(.23) for N , 1.04(.09) for O , 4.74(.45) for Ne , and 2.2(1.6) for Si , resulting in temperatures of 0.180(.022) keV, 0.244(.007) keV, 0.296(.007) keV, and 1.13(0.50) keV, respectively. These results are based on tables of line-emissivities given by Mewe et al. (1985) and MEKAL. My ratio for O , as well as my $N_{\rm H}$ value are a bit higher than the values used by Güdel & Telleschi (2007) in their Fig. 2.

As ratios (f+i)/r in the He-like ions N VI, O VII, and Ne IX the values of 1.06(.55), 0.994(.146), and 0.978(.074) were established, respectively. These values correspond to temperatures of $\mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$ }}}\hbox{$<$ }}}$0.28 keV, 0.129(.043) keV, and 0.18(.05) keV, as derived from tables given by Porquet et al. (2001). They all cover nicely the values obtained in the multi-T fit (Sect. 3.2).

4.3.2 Density

By definition (see Sect. 3.2), knowledge of the electron density $n_{\rm e}$ offers the possibility to derive the emitting volume, assuming a source with spatially constant $n_{\rm e}$.

The ratios (f/i) of C V, N VI, O VII, Ne IX, Mg XI, and Si XIII serve as density diagnostic. The line triplet of S XV is not resolved and given as one feature in Table 3. Their lines are located in the range from 5 to 45 Å covered by LETGS. Here, I will focus on the ions N VI, O VII, and Ne IX. They are the stronger of the line triplets and well resolved.

Considering the measured line fluxes (see Table 3), the f/i ratios of N VI, O VII, and Ne IX are $\mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$ }}}\hbox{$<$ }}}$.46, 0.15(.06) and 0.63(.07), respectively. These ratios result in densities of $\mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$ }}}\hbox{$>$ }}}$ $6.5\times 10^{10}$ cm-3 for N VI, $n_{\rm e}= 8.2(-2.4,+~6.4)\times 10^{11}$ cm-3 for O VII, and of $n_{\rm e}=3.5(1.0)\times 10^{12}$ cm-3 for Ne IX (Porquet et al. 2001). This is in fair agreement with the values obtained from the multi-temperature fit (see Table 2).

5 Discussion

In this section I will discuss my results and pay special attention to differences and agreements between my results and results given in former publications. It is difficult to conclude where differences come from: whether they originate in the use of three different instruments or different models, in erroneous interpretations, or in actual variation of the source.

The three obtained CIE-temperatures (0.127(.012), 0.264(.007), and 1.45(.43) keV) are in very good agreement with values established by Stelzer & Schmitt (2004) (0.172, 0.272, and 1.72 keV, respectively) and by Robrade & Schmitt (2006) (0.24, 0.60, and 1.37 keV, respectively). To improve the shape op the continuum a blackbody component of 0.187(.016) keV was introduced that replaces the temperature bin of 0.60 keV found by Robrade & Schmitt. From Table 2 in Günther et al. (2007) we notice that this CIE temperature component (0.72 keV) is not very significant. The introduction of my blackbody component models the shape of the continuum from 5 to 15 Å better. Thanks to that it offers the possibility to fit some weak line features of silicon and sulphur properly. It reduces $\chi^2$ significantly; from 537 (applying parameters given by Robrade & Schmitt 2006) to 334 (this work) for 299 degrees of freedom and just covering the wavelength range from 5 to 35 Å. It is possible that my ``blackbody'' continuum describes a continuum formed by many weak lines not present in the databases.

My fitted value of $N_{\rm H}$ is higher than the values given by Stelzer & Schmitt (2004) and Robrade & Schmitt (2006). This higher value is necessary to suppress the appearance of lines above 40 Å in the model that are absent in the data. A fit with a Gaussian at the theoretical position of the Si XII lines at 44 Å (where no feature is visible) yields an upper limit of 7$\times$10-6 photons cm-2 s-1. These lines also appear in the model when using the parameter sets given by Stelzer & Schmitt (2004) and Robrade & Schmitt (2006), based on RGS, to describe the LETGS spectrum. Their parameters, however, were based on data sets limited to 35 Å, not covering the higher wavelength range where these lines are located.

The abundance pattern (A/O) is quite well the same as for the two other papers, especially concerning the Model B in Stelzer & Schmitt (2004). Neon and nitrogen are enhanced. My carbon abundance relative to oxygen is a bit higher than theirs, but that is due to my higher $N_{\rm H}$ value. My absolute abundance values, however, are far higher (by a factor 2-3) than those given in the two other papers. This results in correspondingly lower values of the emission measures, as the quantity $A \times EM$ is robust. Despite the big difference in absolute values for the emission measure the overall emission measure shape agrees with Stelzer & Schmitt (2004) and the quantity $A \times EM$ is comparable for the three papers (Stelzer & Schmitt 2004; Robrade & Schmitt 2006; and this work).

The electron density $n_{\rm e}$, derived from O VII and Ne IX is high compared to electron densities observed in stellar coronae. My obtained $n_{\rm e}$ value of about $ 10^{12}~{\rm cm}^{-3}$ corresponds to values given in the other papers. It is related to the temperature regime around .2 keV and considered to be an indication for a process of accretion.

Although the overall pattern of the line fluxes and line flux ratios agree with Stelzer & Schmitt (2004) and Kastner et al. (2002) it is wondering that even for strong features the absolute line flux values differ considerably. This is especially true for the line fluxes given by Stelzer & Schmitt (2004) and Kastner et al. (2002). Thanks to the fact that the ratios agree better than the absolute values the conclusions drawn concerning the density, temperatures, and abundances are the same for the three papers. The underlying reason for the difference in line flux might be related to the use of different instruments or to the variability of the source, due to higher volumes of accreted material. The latter is suggested by the difference in line fluxes of Ne X and Ne IX. The ratio (f+i)/r in He-like ions is close to 1, indicating a plasma in collisional equilibrium.

The measured line positions are slightly shifted to higher wavelength. The corresponding velocity of about 30 km s-1 is higher than the satellite motion or $v_{\rm rad}$ of TW Hya. The lines are broadened, in addition to the instrumental LSF, by an ${\it FWHM}$ of $\sim$0.025 Å. This is broader than comparable lines in the coronal spectrum of Capella. Shift and broadening might be fingerprints of the velocities in the funneling accreted plasma from the disc to the stellar surface, although a calibration bias cannot be ruled out. In the latter case it might affect LETGS as well as RGS.

6 Summary

The overall conclusions based on the LETGS spectrum of TW Hya confirm the statements made by Robrade & Schmitt (2006), Stelzer & Schmitt (2004) and Kastner et al. (2002), concerning the dominance of a high density plasma ( $n_{\rm e}\simeq 10^{12}~{\rm cm}^{-3}$) at temperatures of about 0.2 keV. This is different from coronal cool star spectra for which higher temperatures and lower densities are determined (e.g., Raassen et al. 2003).

The abundances of neon (10$\times$) and nitrogen (3$\times$) are strongly enhanced in comparison with the oxygen abundance.

Thanks to the extended wavelength region of LETGS more information about the spectrum above 40 Å was obtained, resulting in a somewhat higher value of the $N_{\rm H}$, i.e., the sum of the interstellar absorption and the circumstellar absorption just near TW Hya, compared to values used or obtained by Stelzer & Schmitt (2004) and Robrade & Schmitt (2006).

Differences in line fluxes between the papers Stelzer & Schmitt (2004), Kastner et al. (2002), Ness & Schmitt (2005) and this work are an indication for variability in mass accretion rate, but might also be due to instrumental differences.

An additional blackbody component is added to improve the multi-temperature fit and to describe a deficiency in continuum radiation between 5 and 15 Å in the model.

Acknowledgements
The SRON National Institute for Space Research is supported financially by NWO.

References

Footnotes

... version 2.0[*]
http://www.sron.nl/~divisions/hea/spex/version2.0
... NIST[*]
http://physics.nist.gov/PhysRefData/ASD/index.html

All Tables

Table 1:   Observation log of the data of TW Hya.

Table 2:   Parameters of the 3-T fit to the spectrum of TW Hya.

Table 3:   The observed line photon fluxes.

Table 4:   Calculated values of $A \times EM$, based on individual line fluxes.

Table 5:   Comparison of observed line photon fluxes.

All Figures

  \begin{figure}
\par\includegraphics[angle=-90,width=8.8cm,clip]{11374f1.ps}
\end{figure} Figure 1:

The zeroth order lightcurve of the LETGS observation of TW Hya. The time bin size is 1800 s.

Open with DEXTER
In the text

  \begin{figure}
\par\psfig{figure=11374f2a.ps,angle=-90,height=5.5truecm,width=8....
...igure=11374f2b.ps,angle=-90,height=5.5truecm,width=8.8truecm,clip=}
\end{figure} Figure 2:

The LETGS spectrum of TW Hya together with the fitted CIE model, grey line (green in the electronic version). The top panel shows the total spectrum and the bottom panel shows the area between 0 and 40 Å.

Open with DEXTER
In the text

  \begin{figure}
\par\psfig{figure=11374f3a.ps,angle=-90,height=5.5truecm,width=8....
...igure=11374f3b.ps,angle=-90,height=5.5truecm,width=8.8truecm,clip=}
\end{figure} Figure 3:

Line fluxes versus wavelength (top) and versus log T (bottom) from: - Stelzer & Schmitt (2004), * Kastner et al. (2002) and $^{\circ }$ This work (black, red, and blue in the electronic version, respectively). Fluxes are given in $10^{-6}~{\rm photons~cm}^{-2}~{\rm s}^{-1}$.

Open with DEXTER
In the text


Copyright ESO 2009

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.