J.-U. Ness1 - M. Güdel2 - J. H. M. M. Schmitt1 - M. Audard3 - A. Telleschi2
1 - Universität Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany
2 - Paul Scherrer Institut, Würenlingen & Villingen, 5232 Villingen PSI, Switzerland
3 - Columbia Astrophysics Laboratory, 550 West 120th Street, New York, NY 10027, USA
Received 24 March 2004 / Accepted 5 July 2004
Spatial information from stellar X-ray coronae cannot be assessed directly, but scaling laws from the solar corona make it possible to estimate sizes of stellar coronae from the physical parameters temperature and density. While coronal plasma temperatures have long been available, we concentrate on the newly available density measurements from line fluxes of X-ray lines measured for a large sample of stellar coronae with the Chandra and XMM-Newton gratings. We compiled a set of 64 grating spectra of 42 stellar coronae. Line counts of strong H-like and He-like ions and Fe XXI lines were measured with the CORA single-purpose line fitting tool by Ness & Wichmann (2002). Densities are estimated from He-like f/i flux ratios of O VII and Ne IX representing the cooler (1-6 MK) plasma components. The densities scatter between from the O VII triplet and between from the Ne IX triplet, but we caution that the latter triplet may be biased by contamination from Fe XIX and Fe XXI lines. We find that low-activity stars (as parameterized by the characteristic temperature derived from H- and He-like line flux ratios) tend to show densities derived from O VII of no more than a few times 1010 cm-3, whereas no definitive trend is found for the more active stars. Investigating the densities of the hotter plasma with various Fe XXI line ratios, we found that none of the spectra consistently indicates the presence of very high densities. We argue that our measurements are compatible with the low-density limit for the respective ratios ( cm-3). These upper limits are in line with constant pressure in the emitting active regions.
We focus on the commonly used Rosner et al. (1978) scaling law to derive loop lengths from temperatures and densities assuming loop-like structures as identical building blocks. We derive the emitting volumes from direct measurements of ion-specific emission measures and densities. Available volumes are calculated from the loop-lengths and stellar radii, and are compared with the emitting volumes to infer filling factors. For all stages of activity we find similar filling factors up to 0.1.
Key words: X-rays: stars - stars: coronae - stars: late-type - stars: activity - techniques: spectroscopic
The magnetic outer atmosphere of the Sun, the corona, was recognized in radio
and X-ray emission. While the radio emission is associated with bremsstrahlung
and cyclotron emission from free electrons in the hot plasma, the X-ray
emission is produced by bremsstrahlung and line emission. Stellar coronal
activity is therefore investigated primarily in these two bands of the
electromagnetic spectrum. We focus on the potentials offered by X-ray
spectroscopy of stellar coronae. Systematic measurements with Einstein, ROSAT,
ASCA, and other satellite-based X-ray missions revealed that all late-type
main sequence stars from type M to F have coronal X-ray emission. Schmitt (1997)
found X-ray surface fluxes covering four orders of magnitude. Coronal activity
appears thus to be a universal process for cool stars of spectral types F-M,
but no correlation between spectral type and degree of activity could be
established. The only fundamental stellar parameter found to be correlated
with the X-ray luminosity is the rotational velocity
age, e.g., Pallavicini et al. 1981), suggesting that a magnetic dynamo process is involved
in producing the X-ray coronae. Some information on the spatial
distribution of coronal plasma was inferred by indirect means such as modeling
of eclipses and rotational modulation
(e.g., Schmitt & Kürster 1993; Güdel et al. 1995,2003; Siarkowski et al. 1996; White et al. 1990), but such analyses can
only be carried out for very special systems with advantageous geometries. X-ray
spectra allow us to gather a more general insight into the physical properties
of a large variety of stellar coronae. The analysis of X-ray coronae has in the
past been possible only with very limited spectral resolution or with low
sensitivity. Plasma densities (as the subject of this work) could not be
measured from these spectra, because information from spectral lines was not
available. Nevertheless, temperature distributions (or emission measure
distributions EMD) and coronal abundances could be estimatedfrom low-resolution
spectra by the application of global fit approaches. A model spectrum
composed of a continuum and all known emission line fluxes formed under assumed
temperature conditions and with assumed elemental abundances (mostly only scaled
to solar abundances) is iterated using one or two temperature components left
free to vary. The information on bremsstrahlung continuum and emission lines is
extracted from atomic databases containing line emissivities as a function of
plasma temperature under assumptions of solar elemental composition and
collisionally ionized plasma. A spectral model is thus composed as the sum of
bremsstrahlung and all lines formed under the assumption of thermal equilibrium,
and model parameters are the equilibrium temperature and elemental abundances.
The first detailed survey of low resolution X-ray spectra for a large sample of
130 late-type stars (B-V colors redder than 0.0) was presented by
Schmitt et al. (1990) using Einstein data. For each spectrum the temperature structure
was obtained with global spectral models. Different approaches were tested
ranging from isothermal plasma with one or two (absorbed) temperature
components to continuous emission measure distributions.
A much smaller sample concentrating on a sample covering the Sun in time was
analyzed by Güdel et al. (1997) using ROSAT and ASCA data. Their sample
consisted of nine G stars in different stages of evolution.
From MEKAL and Raymond-Smith models they found that the older stars (with slow
rotation) contained only a single cool temperature component in the
emission measure distribution while the younger, more active stars had a
bimodal emission measure distribution with a similar cool component
and an additional hot component apparently independent of the cooler component
suggesting an additional heating mechanism for the more active stars. The
hotter temperature component
was found to scale with the X-ray
|Figure 1: Density-sensitive Fe XXI line flux ratios as predicted by the APEC database for the temperature range between log(T)=6.6( lower borders) and 7.6 ( upper borders). In the right panel we show the term diagram explaining the formation of the 102.22 Å line as an example.|
Despite relatively good spatial resolution with previous X-ray satellites, the technology of their photon counting detectors reached only moderate spectral resolution in the X-ray regime. Very few attempts were made to use dispersive gratings converting spatial resolution into spectral resolution with the potential to resolve individual emission lines. For grating spectroscopy sufficient light is needed and it is therefore only feasible for the brightest sources with long exposure times. In the extreme ultraviolet range this technique has been successfully applied, e.g., with the EUVE mission. Many low-temperature Fe lines from stages Fe X to Fe XVI can be measured with EUVE and are sensitive to densities (Schmitt et al. 1994). Also, at higher temperatures density-sensitive Fe lines can be measured with EUVE (Dupree et al. 1993), but they are only sensitive for relatively high ranges of density cm-3. A summary of results from EUVE measurements is presented by Bowyer et al. (2000). Since the density information is inferred from mostly weak lines the results tend to be ambiguous. The hotter plasma regions of Capella and AB Dor were investigated by Sanz-Forcada et al. (2003a); Dupree et al. (1993) using Fe XIX-XXII line ratios. Densities as high as 1013 cm-3 were reported, suggesting very compact emission regions. However, later Chandra LETGS measurements of Capella contradict these results (Mewe et al. 2001). Because of low sensitivity, EUVE data could only be obtained for some bright sources such as Capella or AB Dor. The apertures of the new missions Chandra and XMM-Newton are large enough to allow grating spectroscopy for many stellar coronae, and X-ray spectra of unprecedented spectral resolution are available. We are able now to measure densities and temperatures from line flux ratios. In this paper we describe the formalism for calculating densities from selected emission line fluxes in Sect. 2. We then present the results from density-sensitive ratios representative for O VII and Ne IX plasma as well as carbon-like Fe XXI lines in Sect. 4. In Sect. 6 we discuss our results.
Spectroscopic information on coronal plasma densities for stars other than the Sun first became possible with the advent of high resolution EUVE spectra ( ; Abbott et al. 1996) that allowed the separation of individual spectral lines. Even with this resolution the available diagnostics often tended to be ambiguous because of the poor signal-to-noise ratio of observed spectra or due to blended lines. Studies of density-sensitive lines of Fe XIX to Fe XXII in EUVE spectra of the RS CVn system Capella revealed some evidence for high densities of to 1013 cm-3at coronal temperatures near 107 K reported by Dupree et al. (1993), but often the densities derived from different lines or ions varied greatly despite the similar formation temperatures. More recent analyses of EUVE Fe XIX-XXII lines (Sanz-Forcada et al. 2003a) and Chandra HETGS Fe XXI lines for AB Dor (Sanz-Forcada et al. 2003b) also returned high densities , but again, not consistently for all considered lines. Since many of the lines used in the analyses are intrinsically faint, even the HETGS data suffer from rather low signal-to-noise. Also the consequences of unidentified blends can be more severe for faint lines, which could be the reason for the discrepant densities derived from Fe lines in the same ionization stage.
The spectrometers on board the X-ray telescopes Chandra and XMM-Newton make it possible to measure emission lines at wavelengths shorter than 95 Å with far better signal to noise and spectral resolution. In particular, the Fe L-shell and M-shell lines and lines of the He-like ions from carbon to silicon are available for density measurements. A few of the density-sensitive Fe lines measurable with EUVE in the 120 Å range can also be measured with the Low Energy Transmission Grating Spectrometer (LETGS) on board Chandra but only upper limits were found by Mewe et al. (2001) ( cm-3). They point out that the Fe XIX to Fe XXII line ratios are only sensitive above 1011 cm-3, so that no tracer for low densities for the hotter plasma component is available, neither with EUVE nor with Chandra or XMM-Newton.
The XMM-Newton and Chandra grating spectra allow high precision measurements of individual line fluxes and line flux ratios. Line fluxes are used to compute emission measures in specific lines, which can be combined to differential emission measure distributions (e.g., Schmitt & Ness 2004; Ness et al. 2003a). Ratios of certain line fluxes are sensitive to temperatures or densities and can be used to describe local conditions in stellar coronae describing the physical conditions in the line-forming regions. In this paper we analyse the He-like triplets of oxygen and neon that probe the low-temperature component of the plasma, and four density-sensitive Fe XXI lines that probe the hotter component of the plasma in a larger sample of stars.
With six electrons (1s22s22p2) the Fe XXI ion is carbon-like. Its ground configuration is split into the states 3P0 (ground state), 3P1, 3P2, 1D2, and 1S0. Transitions within the ground configuration, which are naturally forbidden by definition, do occur; for example, the 3P1-3P0 transition is located at 1354 Å observable with HST (e.g., Linsky et al. 1998; Robinson et al. 1994). Fe XXI at UV and X-ray wavelengths involve transitions from the 3S1 and 3D1 states of the excited configuration 1s22s2p3 (see term diagram in right panel of Fig. 1).
Collisional excitations from the ground state 3P0 predominantly occur into the excited 3D1 state, producing the strong extreme ultraviolet (XUV) Fe XXI at 128.73 Å, which is essentially independent of density and can be used as a reference line reflecting the Fe abundance and the ionization balance for the stage Fe XXI. Collisional excitation from the ground state into the 3S1 state of the excited configuration is much less likely than collisional excitation from the "excited'' ground state 3P2-3S1. The population of the state 3S1 of the excited configuration therefore depends critically on the population of the state 3P2. In a low-density plasma the population of the latter will be small and eventually decay radiatively into the 3P0 level; in a high density plasma the 3P2 state will be collisionally depopulated into the 3S1 level of the excited configuration, which in turn decays radiatively into 3S1-3P2 (102.22 Å) or 3S1-3P1 (97.88 Å). The latter two transitions are therefore density-sensitive, because they depend on the density-sensitive population of the 3P2 state. Similar considerations apply to the ground- and excited levels 3P1 and 3P2 (Mason et al. 1979) leading to transitions at 3P1-3P1 (117.5 Å), and 3P2-3P2 (121.21 Å). All these lines are in the band pass of the Chandra LETGS and can be used to estimate densities of the Fe XXI emitting plasma at MK. In Fig. 1 we show predicted line flux ratios as a function of density (theoretical emissivities were taken from the APEC database, e.g., Smith et al. 2001). Although these lines are all in the same ionization stage, the line flux ratios depend slightly on the plasma temperature. This is illustrated by the associated curves for a low temperature log(T)=6.6 and a high temperature log(T)=7.6 for each ratio in Fig. 1. Clearly, temperature primarily affects the low-density limit. To be conservative we will use the high-temperature theoretical ratios for comparison with our measured ratios, yielding higher theoretical line flux ratios.
The derivation of densities with He-like triplets originated
in solar observations (Gabriel & Jordan 1969). The excited state transitions
1P1, 3P1, and 3S1 to the ground state
1S0 are by convention named resonance line (r), intercombination line
(i), and forbidden line (f), respectively. The ratio f/i is density-sensitive
due to collisional excitations 3S
P1 in high-density
plasmas. These transitions compete with radiative transitions induced by
possible external radiation sources, namely the stellar surface. An analytical
description was given by, e.g., Gabriel & Jordan (1969); Blumenthal et al. (1972):
Since the gratings were to a large extent designed to measure the He-like
triplets of N VI up to Si XIII it is not surprising that such
analyses have been carried out for quite a few individual sources
(e.g., Ness et al. 2001; Stelzer & Schmitt 2004; Güdel et al. 2001; Audard et al. 2001; Ness et al. 2003a,2002c). Especially the
O VII triplet has been analyzed, because the lines are strongest and
least blended. From the measured f/i-ratios densities were calculated. For
Capella and Procyon the densities were found to be at the lower end of the
sensitivity range (Ness et al. 2001). The first study of f/i-ratios in a
sample of stellar coronae was carried out by Ness et al. (2002d), who measured
f/i-ratios for all He-like ions measurable with the LETGS for a sample of
ten stellar coronae. For inactive stars (with low
only low density limits were found, while for the active stars in their
sample higher densities were encountered, although a little surprisingly
only low-density limits were measured for some RS CVn stars.
Ness et al. (2002d) concluded that for the
high-temperature plasma LETGS data offer no conclusive tracer for densities
because of blending problems with the Ne IX triplet, which
is better measured with the HETGS (see also Ness et al. 2003a).
Especially for the more active stars only a very small fraction of the X-ray emitting
plasma is produced in the temperature range where the O VII triplet
|Figure 2: Comparison of theoretical predictions for density sensitive f/i ratios for O VII ( top panel) and Ne IX ( bottom panel) from Eq. (2), the CHIANTI database, Porquet et al. (2001), and the APEC database. Filled areas represent the varying electron temperatures as given in the legends (hatched: CHIANTI, dark-grey shaded: Porquet et al. (2001), light-shaded: APEC). Lower temperatures yield lower f/i ratios.|
The purpose of this paper is the analysis of He-like f/i ratios for a large sample of X-ray spectra obtained with the Reflection Grating Spectrometer (RGS) on board XMM-Newton and the Low Energy Transmission Grating (LETGS) and the High Energy Transmission Grating Spectrometers (HETGS; consisting of the Medium Energy Grating MEG and the High Energy Grating HEG) on board Chandra. The LETGS spectra are also used to measure density-sensitive Fe XXI lines. The aspect of resonant line scattering in stellar coronae has been addressed with a large sample of grating spectra (Ness et al. 2003b) and strong evidence is found that opacity effects can generally be neglected in coronal plasmas. We apply the same procedures for data reduction as in Ness et al. (2003b). For internal consistency, the LETGS spectra presented by Ness et al. (2002d) are re-analyzed and included in our sample. We focus on density measurements with the O VII and Ne IX triplets and Fe XXI ratios, and estimate systematic emitting volumes and filling factors for the different coronae.
We prefer to measure line counts with a program developed specifically for this task named named CORA (Ness & Wichmann 2002). The lheasoft package XSPEC can also do the job, but for a large number of different spectra the CORA program is more efficient.
The CORA program measures line counts from the raw spectrum, i.e., the instrumental background is not subtracted and the background spectrum is instead added to the iterated model spectrum. While the measurement of line counts for the Chandra gratings is straightforward (Ness et al. 2003b), more difficulties arise for RGS spectra. For Capella and AB Dor we measured RGS line counts with both CORA and XSPEC (convolving a -function profile with the line spread function) and we found consistent results within the errors. However, line fluxes obtained with the CORA program were found to be systematically lower than the results obtained with XSPEC. A possible reason for these systematic discrepancies could be that XSPEC uses a wavelength-dependent response matrix, which provides a more accurate instrument description than CORA (which uses only approximate analytical line profiles). We tested the effects arising from different treatments of the instrumental line profiles by exporting the line profiles used for the XSPEC fits into CORA and found similar results compared to what we obtained using the analytical Lorentzian profile function, without any systematic trend. Obviously the Lorentzian used by CORA is an adequate representation for the RGS line profiles for our purposes.
Another source of systematic errors is the estimation of source background.
As source background we consider the combination of continuum emission and
the sum of unresolved weak lines above the measured instrumental background.
This problem has been extensively described by Ness et al. (2003b) who developed
a modified median routine providing a parameterized method to determine a source
continuum value. We present an alternative approach, viz. a optimization of the source continuum.
In the CORA program the continuum is assumed to be constant in sufficiently small wavelength regions, which is justified for individual line fitting with nearby lines. However, use of the median value as a source background value as applied in the CORA program is only valid as long as more than 50% of the bins belong to the desired background. A significant difference between the LETGS, for which the CORA program was originally developed, and the RGS spectra in our sample is the profile function. The wings of the RGS line profile are wider and contribute to the continuum value obtained with the median function and CORA will therefore return an overestimated continuum and thus underestimated line counts. Our modified method uses the median value as a start value and we minimize the value iterating only the source background value. The other parameters of the model spectrum (wavelengths and line widths) are kept fixed to the given initial values but the line counts are optimized with the implemented likelihood method in each iteration step.
We tested the new procedure and found it stable and particularly useful for wavelength regions with few emission lines. For RGS spectra this method returned systematically lower source background values, while for LETGS spectra these values were consistent with the median values. The discrepancies with the XSPEC results are significantly reduced. We point out that even without the new procedure the discrepancies are not significant within the errors. We also tested the fit procedure for the Fe XVII lines measured by Ness et al. (2003b), but found that it did not work well. We attribute these difficulties to the large number of lines in the 15 Å region. The disadvantage of our approach is that all line features not selected to be measured as emission lines increase the source background value in the attempt to minimize the value in those wavelength regions. We thus conclude that the parameterized median value must be used for line-crowded regions, while the approach represents a non-parameterized procedure for measuring in wavelength regions where all line features are selected to be fitted.
For our study of stellar coronae we selected a sample of coronal sources
as large as possible. We gathered 64 grating spectra of 42 stellar coronae;
we specifically discuss 22 RGS spectra, 16 LETGS spectra, and 26 HETGS spectra.
The reduction of the spectra has been carried out with standard SAS and CIAO
routines and is described by Ness et al. (2003b). The details of the observations
with exposure times and derived luminosities are summarized in Table 1
of Ness et al. (2003b). X-ray luminosities, averaged over the complete observations
(i.e., including flares), are obtained by summing up all
dispersed photons in the wavelength range 5.15-38.19 Å after
In order to
have the largest possible data sample for our systematic analysis of
densities we extracted all additional stellar spectra that were publicly
available by 31 January 2004 from the Chandra archive. These additional
observations are listed in Table 1 using the same format as in
Ness et al. (2003b). For internal consistency the LETGS observations of Algol and
Capella were reanalyzed for this work and the results are also listed in
Table 1 for comparison. We extracted effective areas for flux
conversion from the Capella observations
and used these areas for all observations. Since this paper deals only with
line ratios (and therefore uses only ratios of effective areas), this procedure
is sufficiently accurate. Nevertheless we compared these effective
areas with individually extracted areas and found sufficient agreement for all
instruments. The complete stellar sample used for this work is described in
Sect. 3.2 and all stellar properties relevant for this paper are
listed in Table 2.
|Star||Exposure time [ks]||(1028 erg/s)|
|MEG||1 order||2 order||1 order||2 order|
a5.15-38.19 Å. b5.15-21.5 Å.
For this paper we measure the O VII triplet (21.6/21.8/22.1 Å, log ) and the Ly line of O VIII (18.97 Å, log ) with the RGS1, the LETGS, and the MEG. The RGS2 cannot measure O VII because of chip failure, and the HEG does not cover the O VII triplet. The Ne IX triplet (13.44/13.55/13.7 Å, log ) is severely blended with highly ionized Fe lines (Ness et al. 2003a) and the blends can only be resolved with the HEG and the MEG. Although the RGS and LETGS cover the 13.5 Å region we do not analyze those data here, because de-blending is too complicated and must be done in a future paper. The measured counts and derived f/i ratios (using effective areas as described in Ness et al. 2003b) are listed in Tables A.2 and A.5. The densities derived from Eq. (2) are also listed with additional description in Sect. 4.2.2. For the hot sources for which LETGS spectra are available, we measured the lines at 128.73 Å, 117.5 Å, 102.22 Å, and 97.87 Å. The measured line counts are listed in Table 3 and discussed in Sect. 4.2.1.
|Star||HD/Gl||Spectr. Typea||Distancea||Va||B.C.b||log( )|
aFrom Simbad. bFrom Kaler (1989).
dFrom (used for analysis). eMeasured with ROSAT (5.2-124 Å).
In Table 2 we list all relevant stellar parameters. The spectral type information has been taken from the Simbad database. It can be seen that a broad range of coronae is included in the sample. The sample covers stars with extremely high flare activity, RS CVn systems, and other double systems. The distances are also from Simbad and are based on Hipparcos parallaxes. An important parameter for our analysis is the stellar radius which is needed for scaling the derived coronal loop sizes to typical geometries. For some stars in the sample radii are not available in the literature, and we adopted a procedure to calculate stellar radii from the apparent visual magnitudes V and distances (and thus absolute magnitudes) and spectral types, listed in Table 2. We estimate bolometric corrections and effective temperatures from the spectral type interpolating tables from Kaler (1989). We use the bolometric luminosity , calculated from the absolute luminosity and the bolometric correction to calculate the stellar radius for each star from with the Stefan-Boltzmann constant. In Table 2 we list the effective temperatures thus derived, bolometric luminosities, and stellar radii in comparison to what we found in the literature. Good agreement with most values from the literature is found and we use our derived radii for further analysis. In the last column we list the X-ray luminosity from ROSAT (e.g., Hünsch et al. 1999, wavelength range 5.2-124 Å). In Fig. 3 we compare these values with X-ray luminosities obtained from the new spectra (wavelength range 5.15-38.2 Å) listed in Tables A.1 and A.4. The regression fit has a slope slightly lower than one, indicating that the ROSAT fluxes for the more active stars have been underestimated. With the best-fit regression parameters the discrepancies are no higher than a factor of two.
|128.73 Å||117.5 Å||102.22 Å||98.87 Å|
|Cet||587.4 31.6||178.8 23.1||198.9 24.7||98.20 17.6|
|And||471.4 29.4||139.5 20.7||144.3 24.9||109.7 18.9|
|AD Leo||61.17 14.2||15.48 11.7||21.64 12.9||14.89 9.14|
|Algol||354.0 25.4||120.3 19.4||145.3 19.9||77.86 14.3|
|Capella||1086. 43.1||501.7 33.9||362.0 36.6||69.97 21.1|
|EK Dra||52.02 10.4||28.27 11.5||9.360 12.1||-|
|HR1099||386.8 25.6||156.4 19.6||142.3 22.2||63.43 15.4|
|UX Ari||167.6 19.3||71.58 17.6||35.31 17.0||55.07 13.0|
|YY Gem||46.10 10.5||29.74 10.7||8.531 9.01||17.86 7.87|
In the past, coronal temperatures were estimated via global spectral fits from low-resolution spectra. With high-resolution spectra now available we can additionally determine average coronal temperatures by considering individual line fluxes, that yield a temperature characterizing the formation of, e.g., the O VII triplet. The ratios of line fluxes originating from adjacent ionization stages of the same chemical element allow the determination of temperatures independent of the respective elemental abundances. Since for the conversion of measured line fluxes to plasma temperatures and for the construction of synthetic spectra from global plasma emission models the same databases are used, the systematic errors are basically the same. However, a careful choice of strong lines for the calculation of line ratios can reduce systematic errors, since strong lines are mostly not as much affected by uncertainties as weaker lines. Also, line blends from unidentified weaker lines cannot harm strong lines as much as weaker lines. A detailed discussion of the interpretation of line ratios and global fit approaches is given in Güdel (2004).
The strongest lines in the grating spectra are the H-like and He-like lines of carbon, nitrogen, oxygen, neon, magnesium, and silicon. We measured the line fluxes of these lines in order to calculate temperature-sensitive ratios for these ions for each star in our sample. In Fig. 4 we show our results plotted vs. the total X-ray luminosities (cf. Tables A.1 and A.4 for oxygen); a clear trend can be recognized for all ions. We overplotted best-fit linear regressions and find decreasing slopes for increasing Z, thus increasing formation temperatures of the H-like and He-like states. A correlation of plasma temperature with the degree of activity has long been known (e.g., Schmitt et al. 1990; Schrijver et al. 1984). In terms of a putative emission measure distribution, this trend would indicate a predominance of higher emission measure at higher temperature in stars with a higher activity level. When considering H- and He-like lines of individual elements, this trend should be reflected in a larger ratio of H-like to He-like line fluxes in the more active stars, and this is what we now recover from our analysis.
A continuous temperature distribution suggests a mixture of temperatures making the adjacent lines an ideal means for defining interpolation points fixing the shape of the temperature distribution. Such an approach has been introduced by, e.g., Schmitt & Ness (2004) for finding abundance-independent emission measure distributions. For our purposes we can conclude that the temperatures derived from the line ratios of H-like and He-like lines characterize the plasma conditions around the plasma that produces the He-like lines used for the density analysis.
We measured line fluxes from both density-sensitive He-like lines and carbon-like Fe XXI lines. With the He-like lines we probe only the "cool'' plasma component, while with the Fe XXI densities we probe the "hotter'' coronal components.
|Figure 3: Comparison of ROSAT luminosities (5.2-124 Å) listed in Table 2 with measurements for RGS1, RGS2, LETGS, MEG ( 5.15-38.2 Å), and HEG ( 5.15-21.5 Å) from Ness et al. (2003b) and Table 1. The regression curve (ROSAT) = (Chandra/XMM) with parameters R1 and R2 as given in the legend is overplotted. The slope R2 is near unity, indicating good cross-calibration between instruments.|
As described in Sect. 2.2.1 the Fe XXI density diagnostics are
essentially based on the appearance of certain lines in high-density plasmas.
Our search for densities exceeding
is therefore based on
detections of these density-sensitive lines. In Table 3 we list the
results for our Fe XXI line count measurements. Table 3
reveals that for all stars studied the reference line at 128.73 Å is the
strongest line, yielding flux ratios below one
for all the stars in our sample. In plasmas with densities exceeding
cm-3, the 121.21 Å is expected to be the strongest
Fe XXI line, yet it is detected in none of our sample stars. The "best''
cases for detection are the LETGS spectra of Algol and Capella (see
Fig. 5 left panel), but the statistical significance of the
"features'' appearing at 121.21 Å is very low. Even if these features are
taken as real, the measured line fluxes do not imply densities higher than
1012 cm-3. We also investigated the Fe XXI lines at
142.16 Å and at 145.65 Å for Capella, which were used by Dupree et al. (1993).
They measured flux ratios with EUVE
12.8, but we are unable to detect any significant flux at these wavelengths
in the LETGS spectrum of Capella. The grey shaded areas in
Fig. 5 show the expected spectrum using the line counts for the
128.73 Å line and the EUVE flux ratios; it can be seen that such high
densities would have been measurable with the LETGS.
|Figure 4: Line flux ratios of H-like (Ly) and He-like (r) lines vs. activity indicator for the respective ions of oxygen, neon, magnesium, and silicon. Linear fits yield good approximations with decreasing slopes for higher-Z ions. The right axes give average temperatures derived from comparison of the measured line ratios with temperature-sensitive predictions from the APEC line database.|
In Fig. 6 we plot the line flux ratios measured for our sample
stars between the detected Fe XXI lines at 117.5 Å, 102.22 Å, and
97.87 Å, all with respect to the Fe XXI line at 128.73 Å. To
convert line counts into fluxes we used the effective areas and ISM
transmissions listed in Table 3. Note that we did not consider
individual ISM transmissions
for each star; since this effect is small and differential; the error is
smaller than the statistical error of our measurements. The grey lines in
Fig. 6 represent the line flux ratios computed from APEC for
the case of a low density plasma. As can be seen from Fig. 6 all
the ratios for 117.5 Å/128.73 Å are above the computed low density
limit, however, all the observed line ratios are consistent with a value of
0.25. Since we consider it unlikely that all observed coronae are above the
low-density limit at the same density, a far more plausible explanation
is that all coronae are in the low-density limit and that a flux ratio of
0.25 is a more appropriate value for the low-density limit than the computed
value of 0.16. Similar conclusions apply to the flux ratios of the
102.2 Å/128.73 Å and 97.87 Å/128.73 Å lines, where the computed
low-density values are 0.17 and 0.07, which has to be compared to the observed
values of 0.25 and 0.10. Again, the most plausible explanation is that all
coronae are in the low density limit, which is consistent with the
non-detection of the Fe XXI lines at 121.21 Å, 142.16 Å, and
145.65 Å. We also point out that in no case all Fe XXI line
ratios yield consistent high densities. Individual deviations from the
low-density limit could be due to unidentified blending or other uncertainties
in the atomic data bases. In particular, inclompleteness of atomic databases
result in a bias towards higher densities. Unknown emission lines could mimic
high densities when unexpectedly showing up where we expect to see
Our conclusion is that densities above 1013 cm-3 can definitely be
ruled out, and densities above
cm-3 appear highly
|Figure 5: Measurement of Fe XXI lines for Algol and Capella. At most marginal evidence for the presence of Fe XXI lines at 121.21 Å ( left panel) and at 142.16 Å as well as 145.65 Å ( right panel) is present. The 128.73 Å line is shifted due to calibration uncertainties. The grey shaded areas indicate the expected spectrum for high density plasma as found by Dupree et al. (1993).|
|Figure 6: Line flux ratios for measured Fe XXI lines at 97.9 Å, 102.2 Å, 117.5 Å, 121.21 Å, and 128.8 Å. The grey lines represent the theoretical low-density limits. The average ratios, weighted with the measurement errors, are marked with the black lines.|
Our measurements of line fluxes for O VII and Ne IX lines are used to determine f/i ratios which are converted to electron densities with Eq. (2); the derived densities (and 1higher limits) are listed in Tables A.2 and A.5. The radiation term describing the contribution to f/i ratios from radiatively induced transitions is negligible for Ne IX and for O VII for most of our sources. For O VII we calculate values for the stars with the highest effective temperatures (cf. Table 2) from IUE measurements at 1630 Å. The method is described in Ness et al. (2001,2002b). For Algol, Capella, and Procyon we calculate values for of 2.18, 0.003, and 0.01, respectively. For Algol the source of UV radiation is the companion B star (Ness et al. 2002c) and depending on the phase geometry during the observation can be significantly lower. The density values listed in Tables A.2 and A.5 take radiation effects into account.
Specific problems in the measurement of the Ne IX lines arise from
the complicated blending structure, studied by Ness et al. (2003a)
for Capella using the best SNR data available. According to Ness et al. (2003a)
the intercombination line could possibly be blended with an additional line of
Fe XIX. Unfortunately this line is only predicted by the APEC
line database, but it could not even be resolved with the HEG or in any
laboratory measurements. For Capella, Ness et al. (2003a) found
this line to contribute to the measured flux with about a third of the
total flux measured at 13.55 Å, thus pushing
the f/i ratio into the low-density limit (for Capella). When inspecting our
Tables A.2 and A.5 we find systematically higher densities from
Ne IX line ratios than from O VII line ratios. Since this is
critical for the assumption of constant pressure in the X-ray
emitting structures, we will test whether systematically higher densities
are still found, if the blending is accounted for.
We extracted the respective line emissivities from the
APEC database and show their temperature dependence in Fig. 8.
According to the APEC, blending can be significant at higher temperatures
and for larger Fe/Ne abundance ratios (which are rather small for most
coronal sources). To assess the amount of expected Fe XIX
contamination we adopt a scaling procedure, measuring the line flux of
supposedly isolated strong Fe XIX lines and scale with the theoretical
ratio of line fluxes.
We extracted emissivities for the Fe XIX line at 13.462 Å and found a value for the ratio of the emissivities and the blending line at 13.551 Å of 8. We then measured line fluxes with the HEG, scaled these with the emissivity ratio 8, and in Table 4 we list line counts thus predicted for the 13.551 Å line. These are contrasted with the original measurements and for most sources higher f/i ratios, yielding lower densities, are indeed found. A similar behavior was found when the 13.518 Å line was used. In Table 5 only sources with particularly high densities are listed, and for Ne IX we calculate densities from the corrected f/i ratios. From Table 5 we conclude that the densities derived from the Ne IX triplet are still systematically higher than those derived from the O VII triplet even if blending is taken into account. Hence the O VII- and Ne IX-emitting layers cannot be at the same pressure. We caution, however, that the blending is purely theoretical and all conclusions rely on the accuracy of the APEC database. Additional blending is predicted to occur from an Fe XX line (see Fig. 8), but we did not make estimates for the contribution from this line, because no Fe XX lines are available for a scaling procedure. For further analysis we use the non-corrected f/i ratios for Ne IX.
|Figure 7: Comparison of total X-ray luminosity and O VII ( left) and Ne IX-specific luminosities. The dashed lines indicate lines of equal luminosity.|
|Figure 8: Iron lines blending the Ne IX intercombination line in high-temperature plasmas predicted by the APEC database assuming solar elemental abundances.|
|Star||ctsa||i [cts]||f [cts]||new f/i||old f/i|
|O VII||Ne IX|
As an activity indicator specific only for the O VII and Ne IX emitting layers
we calculate an ion-specific X-ray
luminosity from the sum of the three He-like line fluxes r+i+f as described
by Ness et al. (2003b). These luminosities are plotted in Fig. 7 in
comparison with the total X-ray luminosities (also taken from Ness et al. 2003b
and listed in Tables A.1 and A.4). Clearly, the O VII and
Ne IX luminosities strongly correlate with the total X-ray luminosity
and we conclude that these luminosities represent the overall degree of
magnetic activity of the coronae at least as well as the total X-ray luminosities.
This is supported by our finding that the average temperatures derived from
the respective ions correlate with the overall X-ray luminosity. We find
that for the least active stars O VII emission contributes on
average less than 10% to the overall luminosity and Ne IX emission
at most 7% for stars of intermediate activity. For the active
stars the percentage drops to below 3% for oxygen and neon.
From the ion-specific luminosities we can also calculate ion-specific emission
measures, but a temperature structure has to be assumed for this procedure. For
simplicity we assume an isothermal plasma at the peak formation temperature for
each ion (
Note that the emission peak around the He-like ions is very narrow.
The ion-specific emission measure is then calculated from
Figure 11 (for O VII) and Fig. 12 (for Ne IX) show the
central results of our measurements in graphical form.
In the upper panels of these figures we show the measured f/iratios and the densities derived with these measurements with Eq. (2)
vs. the ion-specific luminosities (
). The low-density limit R0 is marked with a vertical
dotted line in the upper left panels, only measurements with f/i-ratios <R0yield actual density measurements. Therefore all f/i measurements resulting in
low-density limits have been marked by light colors. In the bottom left panels
we plot the densities versus the emission measure obtained from the ion-specific
luminosities calculated with Eq. (3). In the plot we include lines of
equal emitting (coronal) volumes
For each ion we estimate available volumes
potentially be filled with coronal plasma. We use stellar radii (cf.
Table 2) and a scale height H:
and H3 are neglected) and compare
with the ion specific emitting volumes.
For an estimate of the height we assume the plasma to be confined in a
uniform distribution of loop structures obeying the loop scaling law by
Rosner et al. (1978, RTV):
|Figure 9: Comparison of coronal volumes obtained from O VII and Ne IX (only MEG measurements). The solid line marks the line of equal volumes.|
|Figure 10: Comparison of temperatures derived from Eq. (6) and from Si H-like and He-like line ratios representing the hot component for those stars where HETGS and LETGS spectra are available.|
From these considerations we derive the available volume
|Figure 11: Measured f/i ratios and derived densities versus the O VII specific X-ray luminosities ( upper two panels). The low-density limit is marked with a vertical dotted line and measurements with only upper limits are marked by light symbols. Lower left panel: densities versus O VII specific emission measure at T=2 MK. Lower right: emitting volume versus available volume (Eq. (7) with from Eq. (6); numbers are listed in Tables A.2/A.5, A.1/A.4, and A.3/A.6).|
|Figure 12: Same as Fig. 11 for Ne IX.|
|Figure 13: Filling factors obtained from ratios of coronal volumes (derived from O VII and Ne IX densities) and assumed available volumes. The latter depend on assumed loop-top temperatures, which we derive from the H-like and He-like line ratios (marked with red bullets) and from (Eq. (6)). As activity indicator we use here the surface fluxes from .|
For inference of sizes in stellar coronae the scaling laws derived for the solar corona can be used to link physical properties with spatial sizes under the assumption of loop-like geometries as identical building blocks. The physical parameters required for the application of the scaling laws are the plasma temperature, emission measure, and density. While good estimates of plasma temperatures and emission measures are available from low-resolution spectroscopy, estimates of plasma density do require high-resolution spectroscopy. The He-like triplets provide a direct way to measure densities by using collisionally induced reduction of the so-called f/i line flux ratio with increasing plasma density. The density-sensitive regimes of the different He-like ions increase with atomic number and the best measurable density range is provided by the O VII triplet and the Ne IX triplet. The N VI and C V triplets are also good density tracers; however, the lines of both are quite faint. C V can only be measured with the LETGS and can be blended with Ne IX/Fe XIX third order lines for the more active stars (Ness et al. 2001). The ions of C V and N VI are produced at rather low temperatures and therefore probe only the cooler coronal plasma. The O VII triplet is very prominent and can be measured by all grating instruments with high precision, but again it represents only the low-temperature regions of a multi-temperature plasma. The Ne IX triplet is sensitive to somewhat higher densities and is produced at higher temperatures than O VII. However, the measurement of Ne IX is difficult because of Fe XIX lines blending with the Ne IX lines, particularly the important intercombination line (Ness et al. 2003a). The most reliable results on the Ne IX lines can therefore be derived from the HETGS data for the Ne IX f/i ratios.
The results of our density measurements can be summarized as follows: first, for O VII the measured f/i-ratios are in the range 1 to the low density limit of 3.95. The coronal source with the hitherto lowest measured O VII f/i-ratio is Algol, which yields for RGS, LETGS, and MEG data, i.e., in three independent measurements. While Algol is indeed very active, it is in our opinion very likely that these low f/i-values are affected by the radiation field of Algol's primary. This is supported by the lack of such low f/i-ratios for Ne IX (for higher-Z He-like ions the effects from UV radiation fields become lower; cf. Fig. 8 in Ness et al. 2002d). Second, for Ne IX the measured f/i-ratios are in the range 1 to the low density limit of 3.4. A handful of stars like 44 Boo, Algol, AR Lac, and EV Lac have the lowest values, but discrepancies appear between measurements with different instruments and even between simultaneous measurements (in MEG and HEG). Third, in no case do we have significant density measurements from Fe XXI. Although all measured line ratios should yield consistent densities, we found that none of our spectra returned consistently high densities. We further found no detections for lines that typically appear in high density plasmas. From the upper limits of our Fe XXI density estimates the typically encountered densities in coronal plasmas are definitely not higher than cm-3. Therefore the Si XIII and the Mg XI triplets will yield only low-density limits in coronal plasmas. Recent systematic studies of Si and Mg He-like f/i ratios by Testa et al. (2004) indeed revealed only low-density limits for Si XIII for all stars in their sample, but some density measurements for Mg XI are also reported. Deviations from the Si f/i low-density limit (systematically higher f/i values as expected) were argued to imply too low a theoretical low-density f/i value. Testa et al. (2004) point out that densities in stellar coronae do not exceed as reported by, e.g., Sanz-Forcada et al. (2003b) for AB Dor. The deviations found for Mg f/i ratios were found to be related to the ratio of X-ray luminosity and bolometric luminosity, but no discernible trend with the X-ray surface flux was found. A particular difficulty was that the Mg XI lines are blended with lines of the Ne Ly series (n>5), increasing the formal measurement errors. Testa et al. (2004) successfully disentangled the lines, but admit that residual Ne blending might still be present.
It is instructive to inspect the "low'' f/i-ratios for O VII and Ne IX where measurements with good SNR and high resolution (i.e., MEG) are available. In Table 5 we list only thosemeasurements where the oxygen f/i-ratios are below 2.2 (within the errors) and the neon f/i-ratios are below 2.0. The peculiar role of Algol becomes apparent; its low O VII f/i-ratio stands out, while the Ne IX f/i-ratio does not. The densities derived from O VII and Ne IX usually differ, the Ne IX densities being higher than the corresponding O VII densities. However, using the MEG values the densities for 44 Boo, EV Lac, and II Peg are consistent, while they are inconsistent for the RS CVn stars Capella and CrB. The case of Capella appears especially striking: While three different measurements (with RGS, LETGS, MEG) yield consistent "high'' values of the O VII f/i-ratio, both MEG and HEG yield consistent "low'' values for the Ne IX f/i-ratio. The Ne IX measurement for Capella has been discussed in great detail by Ness et al. (2003a) and it was found that the intercombination line could be blended with an additional Fe XIX line that cannot be resolved with MEG and HEG. In the case of Capella it turned out that accounting for the predicted amount of blending leads to the low-density limit, thus to densities consistent with the densities obtained from O VII. We note that the densities (for Capella) derived from Ne IX are fully consistent with the upper limits derived from Fe XXI (even without accounting for theoretically predicted blending); at any rate, the case of Capella (and possibly that of CrB) appears somewhat peculiar.
With the measured densities of the cool and hot plasma component we computed the emitting volumes of these plasma components (cf. Tables A.3 and A.6). Surprisingly, these volumes are rather small, for the "best'' data sets with the smallest errors one finds volumes between 10 29 - 1030.5 cm3for O VII and somewhat smaller values for Ne IX. For example, for EV Lac we find cm3 for oxygen and neon; assuming a filling factor of unity, one would obtain a coronal scale height of 100 km, which appears pathologically small. We therefore conclude that the filling factor is far away from unity, and that conclusion is substantiated by more sophisticated calculations of the coronal volume. We first calculated the volume of a maximum corona consisting of hydrostatic loops with isothermal temperatures derived from the Ly/He-like line ratios for the respective ions and the measured density, and second, computed maximal coronal volumes from the relationship between the temperature of the hotter plasma and the total stellar X-ray luminosity reported by Güdel et al. (1997) and, again, find small overall filling factors.
So far, the discussion assumed isolated temperature components. Under reasonable physical conditions one expects the pressure to stay approximately constant along any magnetic field line since the high observed temperatures imply large pressure scale heights. Under this assumption the high-density Ne IX emission regions cannot be magnetically connected with the high-density O VII emission regions because the Ne densities would have to be below the O densities, the opposite of what is observed. On the other hand, assuming isobaric loops, any loop emitting in Ne IX will also emit in the O VII lines; also, the density of these O VII layers will be even larger than those of the Ne IX emission layers. Therefore the observed O VII emission would have to be composed of at least two components, a high-temperature, high-density component, which contributes rather little to the observed flux in the forbidden O VII line, and a lower-temperature, lower-density component, which contributes to the bulk of the forbidden O VII line. A detailed modeling in terms of physically consistent loops is beyond the scope of this paper; here, we just consider the following numerical example of EV Lac. The observed ratios of the fluxes of the He-r line of Ne IX and the Ly lines imply temperatures of about log , the observed Ne f/i-ratio of 2.92 implies densities of about 1011.1 cm-3. The O VII emission would then be located at a density of 1011.6 cm-3, which in turn would lead to an O VII f/i-ratio of 0.29 for this material. Decomposing the observed O VII f and i emission into a "hot'' component (with f/i = 0.29) and a "cool'' component (with an assumed low-density limit of f/i = 3.95), so that the overall f/i-ratio is equal to the observed f/i-ratio of 1.52, leads to the following numbers: , , , , i.e., a situation where the f-line is dominated by the low-pressure component, but the i-line by the high-pressure component. The filling factor of the low-pressure component is undetermined and could be large, while the filling factor of the high-pressure component, which contributes most of the flux, would definitely be quite small and dominant by virtue of its high density.
What, then, do we learn about coronal structure? First, an inactive corona is generally dominated by cool plasma (1-4 MK), and this plasma never covers a large fraction of the surface as is the case for the non-flaring Sun. For active stars, cool plasma might actually cover a larger fraction of the stellar surface, but more strikingly, a hot component appears. This hot component has been known since the earliest stellar X-ray observations, but its nature has been debated. Its characteristic temperature seems to be correlated with the activity level, and in extremely active stars it reaches temperatures that on the Sun are known exclusively during flares. In our investigation, we cannot confine the extent of the hot plasma because we do not find conclusive indications for a definitive density of this plasma component. We rather argue that we measure upper limits. Yet we have added two further pieces of information: first, the cooler plasma component cannot cover a large fraction of the stellar surface. And second, there must be spatially separate plasma components at different temperatures (e.g., those detected by the O VII and the Ne IX density analysis). We tentatively argue that the hotter plasma loops fill the space between cooler loops until much of the corona is dominated by the hot plasma. Why, then, do hot loops become progressively more important as the stellar activity increases? As the magnetic activity level and consequently the surface magnetic filling factor increases, the coronal magnetic fields become denser, leading to increased interactions between neighboring field lines, which leads to increased heating (Güdel et al. 1997). The increased heating rate inevitably drives chromospheric material into the loops until an equilibrium is attained. The X-ray luminosity of the hot plasma thus rapidly increases as we move to more active stars. In their most extreme form, such interactions lead to increased levels of flaring, again resulting in increased amounts of hot, luminous plasma. In this picture, the cooler loops are post-flare loops that are still over-dense while returning to their equilibrium state.
Whatever the cause for the increased heating, a relatively cool component, corresponding to typical active regions as seen on the Sun, appears to be present in all stars at a similar level of surface coverage. Our survey has shown that this component reveals densities that may exceed 1010 cm-3 but a clear systematic trend with the overall activity level does not seem to be present.
This work is based on observations obtained with Chandra and XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and the USA (NASA). This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. J.-U.N. acknowledges financial support from Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) under 50OR98010. AT and MG acknowledge support from the Swiss National Science Foundation (grant No. 2000-066875).
|Line flux (10-13ergcm-2s-1)||log(EM/cm-3)|
|star||Instr.||[1028 erg/s]||O VIII||O VII(r)||[1028 erg/s]||[1028 erg/s]||O VIIa||Ne IXb|
|O VII||Ne IX|
|Star||Instr.||i [cts]||f [cts]||f/i||log()||i [cts]||f [cts]||f/i||log()|
|O VII||Ne IX|
aEmitting coronal volumes from Eq. (4). bAvailable volumes from Eq. (7).
cFilling factor .
|Line flux (10-13ergcm-2s-1)||log(EM/cm-3)|
|Star||Instr.||[1028 erg/s]||O VIII||O VII(r)||[1028 erg/s]||[1028 erg/s]||O VIIa||Ne IXb|
|O VII||Ne IX|
|star||Instr.||i [cts]||f [cts]||f/i||log()||i [cts]||f [cts]||f/i||log()|
|O VII||Ne IX|
aEmitting coronal volumes from Eq. (4). bAvailable volumes from Eq. (7).
cFilling factor .