A&A 368, 888-900 (2001)
DOI: 10.1051/0004-6361:20010026
R. Mewe1 - A. J. J. Raassen1,2 - J. J. Drake3 - J. S. Kaastra1 - R. L. J. van der Meer1 - D. Porquet4
1 - Space Research Organization Netherlands (SRON),
Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands
2 -
Astronomical Institute "Anton Pannekoek'', Kruislaan 403,
1098 SJ Amsterdam, The Netherlands
3 -
Smithsonian Astrophysical Observatory, 60 Garden Street,
Cambridge, MA 02138, USA
4 -
CEA/DSM/DAPNIA, Service d'Astrophysique, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France
Received 12 October 2000 / Accepted 15 December 2000
Abstract
We report an analysis of the X-ray spectrum of Capella from 6 to 175
Å obtained with the Low Energy Transmission Grating Spectrometer (LETGS) on board
of the X-ray space observatory CHANDRA. Many emission line features
appear that can be resolved much better as compared to former
instruments (EUVE and ASCA). Coronal electron densities ()
and
temperatures (T) of brightly emitting regions are constrained by an
analysis of ratios of density- and temperature-sensitive lines of
helium-like ions and highly ionized iron atoms. Lines emitted by
e.g., O VII & VIII, Mg X-XII, Si XII-XIV, Fe IX, X & XV-XXIII are
used to derive T. Line ratios in the helium-like triplets of C V, N VI,
O VII, Mg XI, and Si XIII yield T in the range 0.5-10 MK, and
in the range 109-1013 cm-3. The Fe IX/X ratio
yields
0.9 MK, while lines from Fe XVIII to XXII give
6-10 MK. Flux ratios of Fe XX-XXII lines indicate for the
electron densities an upper limit in the range
(2-5) 1012 cm-3.
From line ratios of Fe XVII and Fe XVIII we derive constraints on the
optical depth
of
1-1.5.
An emission measure distribution is derived from Fe line intensities.
Results for element abundances (relative to Fe) from a 4-T model are:
O and Ne/Fe about solar, N, Mg and Si/Fe
2
solar.
The results for Tand
are described in terms of quasi-static coronal loop models
and it is shown that the X-ray emission originates from compact
structures much smaller than the stellar radii.
Key words: stars: individual: Capella - stars: coronae - stars: late-type - stars: activity -
X-rays: stars -
atomic processes vspace-1mm
The G1 III + G8/K0 III binary Capella with a 104.0233-day orbit at
a distance of 13 pc (Strassmeier et al. 1993) is one of the brightest
non-degenerate stellar X-ray sources. Observations by Linsky et al. (1998) with the GHRS spectrometer on the Hubble Space Telescope
(HST) on Sept. 1995 and on April 1996 indicate that the optically
brightest (by
)
and more fastly rotating (v sin i=36 kms-1) G1
component (mass 2.6
,
radius 9.2
)
is about 3-10
times brighter in the chromospheric and transition-region UV lines
than the slowly rotating (v sin i=3 kms-1) G8/K0 star
(2.7
,
12.2
). The Sept. 1995 observations showed that
the coronae around the G1 and G8 stars were equally bright in the hot
(107 K) coronal Fe XXI 1354 Å line, but recent observations with
the HST STIS by Johnson et al. (2001) on October 1999 indicate that the
G8 corona was at least 5 times fainter than the G1 corona in the Fe XXI
line, the latter having changed relative to the Sept. 1995 observations
by only 50% or less. These results suggest that the hot corona of
the G8 star is quite strongly variable whereas that of the G1 is not.
Due to the coronal structure of Capella
its spectrum is rich in
emission lines in the soft X-ray region and has been well studied by many
earlier instruments with lower spectral resolution
(for references see, e.g., Brinkman et al. 2000).
The spectra of Capella obtained with the LETGS are of high
quality and agree with the pre-flight calculated instrumental spectral
resolution (cf. Brinkman et al. 2000).
The wavelength region of the LETGS between 6-42 Å contains a number of
helium-like line "triplets'' from Si XIII to C V which provide density
as well as temperature diagnostics for plasmas with temperatures in
the range 1-10 MK and densities in the range
108-1013 cm-3. The long-wavelength region between
90-140 Å contains a series of lines from -
transitions in Fe XVIII-Fe XXIII which provide density and
temperature diagnostics for relatively hot (
5 MK) and dense
(
1012 cm-3) plasmas (e.g., Mewe et al. 1985, 1991;
Brickhouse et al. 1995).
From EUVE spectra of Capella, Dupree et al. (1993), Schrijver et al. (1995), and Brickhouse (1996) have derived electron densities
in
the range 1012-1013 cm-3 at a temperature
MK. More recently, Brickhouse et al. (2000) have analyzed
the coronal temperature and abundance structure of the coronae of
Capella based on EUVE and ASCA observations covering a broad wavelength
range. Our work presented here
covers the ASCA region and part of the EUVE region,
though in more detail due to the higher spectral resolution
(
0.06 Å) of the LETG spectrometer. In the present
paper, the long-wavelength iron lines have been studied for the
purposes of providing further insights into the coronal
temperature structure as well as for density diagnostics.
Section 2 discusses the observations, the data reduction, the counting
rate/photon flux conversion, and the correction for interstellar
absorption. Section 3 discusses branching ratios for
-independent
lines decaying from the same upper level or within the same ion in
relation to theoretical values and the consequences for optical depth
and source size. In Sect. 4 the analysis of the temperature
diagnostic for helium-like triplets and iron lines is presented and
Sect. 5 discusses briefly the determination of the electron density
from the iron lines and the helium-like triplets.
Section 6 describes the emission measure distribution and
Sect. 7 discusses the results in terms of quasi-static loop models.
![]() |
Figure 1: CHANDRA-LETGS spectrum of Capella observed on 9 Nov. 1999 in +1 order (top), in -1 order (middle), and detail of -1 order spectrum between 12-25 Å (bottom). Spectra are corrected for background. Several prominent lines are labelled with the emitting ions |
Open with DEXTER |
![]() | I+1b |
![]() | I-1b | Fc |
![]() | iond | IDd |
6.20 | 15 ![]() | 6.18 | 13 ![]() | 0.7 | 6.180 | Si XIV | HLy![]() |
6.66 | 43 ![]() | 6.64 | 43 ![]() | 2.0 | 6.648 | Si XIII | Her |
6.70 | 13 ![]() | 6.68 | 19 ![]() | 0.7 | 6.688 | Si XIII | Hei |
6.75 | 21 ![]() | 6.73 | 22 ![]() | 1.0 | 6.740 | Si XIII | Hef |
8.43 | 33 ![]() | 8.41 | 36 ![]() | 2.0 | 8.421 | Mg XII | HLy![]() |
9.18 | 44 ![]() | 9.16 | 50 ![]() | 3.2 | 9.170 | Mg XI | Her |
9.23 | 12 ![]() | 9.22 | 9 ![]() | 0.7 | 9.231 | Mg XI | Hei |
9.32 | 24 ![]() | 9.31 | 23 ![]() | 1.6 | 9.315 | Mg XI | Hef |
12.44 | 33 ![]() | 12.42 | 25 ![]() | 2.3 | 12.430 | Ni XIX | Ne5 |
13.46 | 98 ![]() | 13.45 | 89 ![]() | 7.4 | 13.448 | Ne IX | Her |
14.07 | 39 ![]() | 14.06 | 40 ![]() | 3.1 | 14.045 | Ni XIX | Ne8A |
15.21 | 39 ![]() | 15.20 | 40 ![]() | 3.0 | 15.176 | O VIII | HLy![]() |
18.64 | 20 ![]() | 18.62 | 18 ![]() | 1.7 | 18.627 | O VII | He3A |
18.97 | 275 ![]() | 18.96 | 309 ![]() | 26 | 18.969 | O VIII | HLy![]() |
21.61 | 66 ![]() | 21.61 | 71 ![]() | 9.9 | 21.602 | O VII | Her |
21.81 | 12 ![]() | 21.81 | 15 ![]() | 1.9 | 21.804 | O VII | Hei |
22.11 | 46 ![]() | 22.10 | 49 ![]() | 6.9 | 22.101 | O VII | Hef |
24.79 | 46 ![]() | 24.79 | 47 ![]() | 6.9 | 24.781 | N VII | HLy![]() |
27.00 | 14 ![]() | 27.01 | 17 ![]() | 2.4 | 26.990 | C VI | HLy![]() |
28.44 | 16 ![]() | 28.43 | 14 ![]() | 2.5 | 28.446 | C VI | HLy![]() |
28.79 | 10 ![]() | 28.77 | 13 ![]() | 1.9 | 28.787 | N VI | Her |
29.10 | 5 ![]() | 29.08 | 6 ![]() | 0.9 | 29.084 | N VI | Hei |
29.54 | 9 ![]() | 29.52 | 11 ![]() | 1.8 | 29.534 | N VI | Hef |
33.74 | 43 ![]() | 33.74 | 41 ![]() | 8.7 | 33.736 | C VI | HLy![]() |
40.26 | 20 ![]() | (sum | +1&-1) | 5.0 | 40.268 | C V | Her |
40.73 | 5 ![]() | (sum | +1&-1) | 1.9 | 40.731 | C V | Hei |
41.47 | 5 ![]() | (sum | +1&-1) | 2.2 | 41.472 | C V | Hef |
43.78 | 22 ![]() | 43.78 | 36 ![]() | 3.2 | 43.763 | Si XI | Be5A |
44.03 | 23 ![]() | 44.02 | 29 ![]() | 2.9 | 44.021 | Si XII | Li6A |
44.17 | 41 ![]() | 44.17 | 49 ![]() | 5.1 | 44.165 | Si XII | Li6B |
45.52 | 9 ![]() | 45.51 | 11 ![]() | 1.1 | 45.520 | Si XII | Li7B |
45.69 | 16 ![]() | 45.68 | 20 ![]() | 2.1 | 45.692 | Si XII | Li7A |
49.21 | 14 ![]() | 49.20 | 14 ![]() | 2.1 | 49.222 | Si XI | Be8 |
57.89 | 18 ![]() | 57.89 | 17 ![]() | 3.4 | 57.876 | Mg X | Li5A |
gap | 63.29 | 21 ![]() | 2.4 | 63.295 | Mg X | Li6B |
![]() | I+1b |
![]() | I-1b | Fc |
![]() | iond | IDd |
12.14 | 142 ![]() | 12.13 | 142 ![]() | 11 | 12.124 | Fe XVII | Ne4C |
12.134 | Ne IX | HLy![]() | |||||
12.28 | 57 ![]() | 12.27 | 61 ![]() | 4.7 | 12.264 | Fe XVII | Ne4D |
12.84 | 46 ![]() | 12.83 | 49 ![]() | 3.7 | 12.831 | Fe XX | N16 |
13.83 | 90 ![]() | 13.81 | 85 ![]() | 6.8 | 13.826 | Fe XVII | Ne3A |
14.23 | 163 ![]() | 14.22 | 186 ![]() | 14 | 14.208 | Fe XVIII | F14 |
14.39 | 43 ![]() | 14.38 | 39 ![]() | 3.2 | 14.374 | Fe XVIII | F12 |
14.56 | 36 ![]() | 12.54 | 33 ![]() | 2.7 | 14.540 | Fe XVIII | F10 |
15.03 | 426 ![]() | 15.02 | 444 ![]() | 34 | 15.014 | Fe XVII | Ne3C |
15.28 | 145 ![]() | 15.26 | 147 ![]() | 11 | 15.265 | Fe XVII | Ne3D |
15.47 | 32 ![]() | 15.46 | 37 ![]() | 2.7 | 15.474 | Fe XVII | Ne3E |
15.64 | 41 ![]() | 15.62 | 53 ![]() | 3.6 | 15.628 | Fe XVIII | F7 |
15.84 | 35 ![]() | 15.82 | 29 ![]() | 2.5 | 15.831 | Fe XVIII | F6 |
15.88 | 37 ![]() | 15.86 | 48 ![]() | 3.3 | 15.873 | Fe XVIII | F5 |
16.02 | 107 ![]() | 16.00 | 127 ![]() | 9.1 | 16.002 | Fe XVIII | F4e |
16.006 | O VIII | HLy![]() | |||||
16.09 | 135 ![]() | 16.07 | 151 ![]() | 11 | 16.078 | Fe XVIII | F3 |
16.78 | 259 ![]() | 16.77 | 279 ![]() | 20 | 16.780 | Fe XVII | Ne3F |
17.06 | 268 ![]() | 17.05 | 446 ![]() | 31 | 17.055 | Fe XVII | Ne3G |
17.10 | 309 ![]() | 17.09 | 195 ![]() | 24 | 17.100 | Fe XVII | NeM2 |
17.63 | 38 ![]() | 17.60 | 37 ![]() | 3.6 | 17.626 | Fe XVIII | F1 |
50.34 | 62 ![]() | 50.34 | 38 ![]() | 9.0 | 50.350 | Fe XVI | Na6A |
50.56 | 23 ![]() | 50.56 | 11 ![]() | 3.2 | 50.555 | Fe XVI | Na6B |
52.90 | 4 ![]() | gap | 0.6 | 52.911 | Fe XV | Mg1 | |
54.14 | 26 ![]() | gap | 3.9 | 54.142 | Fe XVI | Na7B | |
54.72 | 37 ![]() | gap | 5.8 | 54.728 | Fe XVI | Na7A | |
gap | 62.87 | 18 ![]() | 2.0 | 62.879 | Fe XVI | Na8B | |
gap | 63.71 | 36 ![]() | 4.2 | 63.719 | Fe XVI | Na8A | |
66.40 | 16 ![]() | 66.36 | 31 ![]() | 3.3 | 66.368 | Fe XVI | Na9 |
91.06 | 14 ![]() | 91.01 | 15 ![]() | 3.2 | 91.020 | Fe XIX | O6F |
93.95 | 160 ![]() | 93.92 | 159 ![]() | 37 | 93.923 | Fe XVIII | F4A |
97.90 | 5 | 97.85 | 2 ![]() | 0.9 | 97.880 | Fe XXI | C6C |
101.58 | 28 ![]() | 101.55 | 27 ![]() | 7.6 | 101.550 | Fe XIX | O6B |
102.26 | 10 ![]() | 102.22 | 9 ![]() | 2.6 | 102.220 | Fe XXI | C6D |
103.97 | 54 ![]() | 103.95 | 54 ![]() | 15 | 103.937 | Fe XVIII | F4B |
108.38 | 87 ![]() | 108.36 | 94 ![]() | 25 | 108.370 | Fe XIX | O6A |
109.98 | 12 ![]() | 109.98 | 9 ![]() | 2.9 | 109.970 | Fe XIX | O6E |
110.67 | 3 ![]() | 110.64 | 4 ![]() | 0.9 | 110.630 | Fe XX | N6D |
111.72 | 9 ![]() | 111.70 | 7 ![]() | 2.2 | 111.700 | Fe XIX | O6C |
114.41 | 3 ![]() | 114.41 | 4 ![]() | 1.1 | 114.410 | Fe XXII | B10A |
117.14 | 24 ![]() | 117.16 | 22 ![]() | 6.4 | 117.170 | Fe XXII | B11 |
117.50 | 8 ![]() | 117.53 | 4 ![]() | 1.7 | 117.510 | Fe XXI | C7B |
118.70 | 17 ![]() | 118.70 | 23 ![]() | 5.7 | 118.660 | Fe XX | N6C |
120.00 | 22 ![]() | 120.01 | 24 ![]() | 6.7 | 120.000 | Fe XIX | O6D |
121.87 | 36 ![]() | 121.88 | 33 ![]() | 11 | 121.830 | Fe XX | N6B |
128.76 | 22 ![]() | 128.74 | 22 ![]() | 11 | 128.730 | Fe XXI | C6A |
132.87 | 37 ![]() | 132.87 | 51 ![]() | 20 | 132.850 | Fe XX | N6A |
135.81 | 11 ![]() | 135.82 | 12 ![]() | 5.4 | 135.780 | Fe XXII | B12A |
171.06 | 39 ![]() | off | detector | 44 | 171.075 | Fe IX | A4 |
174.51 | 6 ![]() | off | detector | 18 | 174.535 | Fe X | Cl4A |
Spectra of Capella were obtained in the calibration phase during
6-12 September 1999 (cf. Brinkman et al. 2000)
and subsequently on 9 November 1999, using the High Resolution
Camera (HRC-S) and the Low Energy Transmission Grating (LETG) on board
the X-ray observatory CHANDRA.
In the present paper we study the second spectrum
(exposure time 85.4 ksec).
The HRC-S is equipped with three microchannel plate detectors (10 cm each) in
one row. The 0th order diffraction signal is aimed at the central detector,
such, that for the +1 and -1 diffraction orders, the wavelength regions that
fall in the gaps between the detector plates do not overlap (see
Fig. 1).
This way the
central plate covers nominally the range of -50 (-1 order) to +60 (+1 order) Å, the outer
plates -55 to -165 Å and +65 to +175 Å. The dispersion is 1.15 Å/mm.
Figure 1 shows the LETGS spectra of Capella in the
+1 & -1 order with the gaps between 62-65 Å and 52-56 Å (+1 and -1 order, respectively)
and a detail of the -1 order spectrum between 12-25 Å.
The spectra were rebinned to intervals of 0.02 Å and corrected for
the background.
Around the dispersion axis a source box is defined. The edge of the box is a
function of the wavelength, small at short wavelengths and wide at long
wavelengths, like a bowtie. The background box is defined as the area from 10
to 39.5 arcsec above and under the dispersion axis. The ratio of the source to
background area ranges from 1:25 at the short wavelengths to 1:6 at the edge of
the detector.
For each bin the number of counts is integrated in the cross dispersion
direction inside the boxes. The net number of source counts is calculated as
follows:
,
,
,
where
counts in background box,
counts in source box, and
net source counts in source box with error ds.
By visual inspection we selected the prominent lines and fitted them with
Gaussian profiles. This has yielded values for the
wavelength, FWHM, and the number of photon counts (including
uncertainties) for each line. By comparing wavelengths and line counts with expected values
we have identified all lines given in Tables 1 and 2.
To represent the background due to continuum radiation we introduce a constant shift.
The intensities of the emission line features in the
spectrum are based on integration over the line profile to obtain
the total counts. For isolated lines
the Gaussian line width turns out to be consistent with the average
theoretical instrumental width, which implies that Doppler broadening
can be neglected compared to the instrumental width.
The measured FWHM of the instrumental profile is about 0.068 Å for wavelengths
Å (and increasing with
at long wavelengths). This
is in good agreement with the results of the pre-flight ground calibration.
Blended lines were fitted by multiple Gauss profiles with the proper instrumental width.
The LETGS instrumental profile is
actually closer to a Lorentzian function than to a Gaussian; we have
investigated the effect of using Lorentzian
profiles in the analysis but usually found no significant differences in
derived line parameters.
Counting rates I and
uncertainties are given in Tables 1 and 2.
The observed counting rates I (10-4 c/s) have been converted to photon fluxes F (10-4 phot/cm2/s) using the
preliminaryin-flight effective area model for the spectrometer (for details of the calibration, see van der Meer
et al. 2001).
The observed flux at Earth was further corrected for interstellar absorption, assuming an interstellar column density
cm-2 (see Schrijver et al. 1995), using the cross sections of Morrison & McCammon (1983).
In Tables 1 and 2 we give the photon flux F averaged over the +1 & -1 order and corrected for interstellar absorption.
In the overlapping wavelength regions ( <35 Å) the obtained fluxes are in agreement with values obtained with the CHANDRA-HETGS by
Canizares et al. (2000) and with data from the RGS on board XMM-Newton (see Audard et al. 2001; cf. their Table 1 for a comparison of fluxes).
For the long-wavelength region (
>90 Å) we have compared the LETGS fluxes with the EUVE fluxes obtained by Brickhouse et al. (2000).
The LETGS/EUVE flux ratio was in the range
0.5 to 0.8 which is explained by a different activity level of Capella (see e.g.,
comparison by Linsky et al. 1998).
![]() |
+1a | -1a | av.b | Nc | Md | Ce | Bf |
Si XII | |||||||
![]() |
0.57 | 0.52 | 0.54(.15) | 0.51 | 0.50 | 0.51 | -- |
Fe XVI | |||||||
![]() |
gap | 0.50 | 0.49(.08) | 0.48 | 0.48 | 0.52 | -- |
Fe XVIII | |||||||
![]() |
3.0 | 2.9 | 2.5(.2) | 2.76 | 2.43 | 2.70 | 2.71 |
Fe XIX | |||||||
![]() |
4.0 | 4.0 | 3.5(.4) | 3.75 | 3.51 | 3.64 | 3.75 |
![]() |
2.3 | 3.0 | 2.6(.4) | 1.98 | 1.80 | 1.95 | 2.00 |
]
![]() |
3.3 | 3.8 | 3.5(1.1) | 2.52 | 2.34 | 2.51 | 2.54 |
![]() |
1.42 | 1.28 | 1.4(.4) | 1.27 | 1.27 | 1.29 | 1.27 |
Fe XXI | |||||||
![]() |
2.0 | 4.0 | 2.8(1.7) | 2.42 | 2.37 | 2.39 | 2.42 |
In Table 3, photon flux ratios of two lines decaying from the same upper level are compared to the theoretical branching ratios. The theoretical values are from the NIST compilation, from MEKAL (Mewe et al. 1985, 1995a) as implemented in the SPEX code (Kaastra et al. 1996), from recent calculations at SRON using the Cowan code programs, and from work by Brickhouse et al. (1995). From this table we notice that the theoretical ratios from MEKAL for the Fe XVIII pair and three of the Fe XIX pairs are somewhat lower than those from the three other theories, which are in good agreement with each other. For all except one the theoretical line ratios in Table 3 are within the experimental statistical errors of the Capella line ratios. For some of the ratios the statistical error is also very large, e.g., the ratio 102.22/97.88 for which the deviation between the +1 and -1 order is apparently large, though the flux for 97.88 Å should only be treated as an upper limit. The only disagreement concerns the ratio 101.55/109.97. However, taking into account an uncertainty of 25% for the theoretical values as indicated by the NIST database (classification "C'') all ratios agree with the theoretical values. Brickhouse et al. (1995) observed a deviation for the ratio 101.55/111.7. She suggests a blend for the line at 101.55 Å. Although in our case the ratios related to the 101.55 Å line are about 25% larger than the theoretical value (which is in agreement with some blending) they are all within the combined statistical and NIST error bars.
Schrijver et al. (1995) and Mewe et al. (1995b) have investigated the
possibility that strong resonance photons are scattered out of the
line of sight in late-type stellar coronae. Branching ratios can be
used to check the assumption of the optical thin model because effects
of resonance scattering would affect the measured branching ratio.
From the fact that the intensities of, e.g., the strong resonance
lines Fe XVIII 93.92 and Fe XIX
108.37 are in good
agreement with the intensities of other lines sharing the same upper
levels, we can derive a constraint on the optical depth taking into
account the systematic uncertainties of the theoretical transition
probabilities A (typical 25% for each A, hence 35% for the branching ratio)
which dominate over the statistical errors (typically 10%). If we
rule out a reduction in the resonance line intensity larger than
about 30%, then on the basis of a simple "escape-factor'' model
(escape factor
for a homogeneous
mixture of emitters and absorbers in a slab geometry (e.g., Kaastra
& Mewe 1995)) we can put a constraint on the optical depth:
yields
for the Fe XVIII
93.92 line.
Line transitions of Ne-like Fe XVII yield the strongest emission lines
seen in the LETGS spectrum in the range 15-17 Å (cf. Fig. 1 and Table 2).
For the very strong Fe XVII 15.014 resonance line we can also
derive a constraint on the optical depth as follows.
In principle, the ratio
of the strong Fe XVII 15.014 Å resonance line (with a large oscillator strength
f=2.66) to a presumably optically thin Fe XVII line with a relatively small oscillator strength
(e.g., 16.780 Å, f=0.101, or 15.265 Å, f=0.593)
yields the optical depth in the 15.014 Å line, applying the escape-factor model.
Solar physicists have used this technique
to derive the density in active regions on the Sun from the optical depth and the estimated
plasma dimension
(e.g., Saba et al. 1999).
Though the 16.780 Å line has the smallest f-value we prefer to use the 15.265 Å line
because the 16.780 Å
line can be affected by radiative cascades.
From the measured Fe XVII 15.014/15.265 photon flux ratio of 2.850.14 we derive for a plasma of
7 MK (as derived from the Fe XVII 15.265/Fe XVIII 16.078 ratio, see Sect. 4.1) an upper limit to the optical depth
of the Fe XVII 15.014 Å line by comparing to the corresponding line ratios
given by experiment (2.8-3.2) (Livermore Electron Beam Ion Trap (EBIT)) or theory (3.3-4.7) (cf. Brown
et al. 1998; Bhatia & Doschek 1992).
Taking the EBIT ratio we obtain from the condition
an upper limit of
,
whereas the theoretical
ratio would give a more conservative upper limit of
from
.
![]() |
= | ![]() |
|
![]() |
![]() |
(1) |
![]() |
(2) |
The above analysis, however, is strictly valid only if all flux is concentrated in one single loop. If one observes a collection of loops distributed more or less homogeneously over the stellar surface, the effects of optical depth tend to average out. There may be significant resonance scattering of photons into the line of sight. In this situation it is possible that there is indeed significant optical depth in certain structures, but that this is compensated for by the azimuthal average over the stellar disk whereby some lines of sight pick up photons that are scattered and others lose photons to scattering (e.g., Wood & Raymond 2000). Nevertheless, the above analysis still serves to demonstrate the usefulness of investigating optical depths to provide further clues on the emitting geometry.
![]() |
+1a | -1a | av.b | Mc | Bd |
Fe XVI | |||||
54.14/54.728 | 0.70 | gap | 0.67(.15) | 0.56 | -- |
Fe XIX | |||||
91.02/101.55 | 0.49 | 0.56 | 0.42(.07) | 0.156 | 0.241 |
91.02/108.37 | 0.16 | 0.16 | 0.130(.014) | 0.057 | 0.082 |
91.02/109.97 | 1.1 | 1.7 | 1.09(.19) | 0.287 | 0.476 |
108.37/109.97 | 7.0 | 10.4 | 8.4(1.4) | 5.00 | 5.48 |
108.37/111.70 | 10.0 | 13.6 | 11.4(2.6) | 6.34 | 6.95 |
Fe XX | |||||
118.66/121.83 | 0.48 | 0.68 | 0.55(.07) | 0.517 | 0.517 |
118.66/132.85 | 0.47 | 0.44 | 0.28(.03) | 0.355 | 0.386 |
121.83/132.85 | 0.97 | 0.65 | 0.55(.06) | 0.679 | 0.754 |
Fe XXII | |||||
117.17/135.78 | 2.1 | 1.8 | 1.17(.2) | 0.853 | 0.991 |
In Table 4 ratios are given for lines that belong to the same Fe ion, but that do not decay from the same upper level. In this case they cannot be compared to the ratios of the transition probabilities, the branching ratio as given by NIST and COWAN. However, comparison with MEKAL and Brickhouse et al. (1995) is still possible, as they take into account the populations of the levels due to collisional excitations and radiative decays. Only density-independent lines are used. In Table 4 all ratios in Fe XIX that include the line at 91.02 Å are far too high. This implies blending by an unknown feature or a theoretical flux value far too low for that line. There is some indication for this, given the considerable difference between the MEKAL and Brickhouse values, suggesting some uncertainties in the calculation of the proper flux value. This situation renders this line at 91.02 Å problematic for density and temperature diagnostics. The ratios including the Fe XX line at 132.85 Å also deviate from the theoretical values. However, 132.85 Å is well known to be blended with a line of Fe XXIII; the relative contributions of Fe XX and Fe XXIII depend on the temperature regime and for Capella are expected to be dominated by Fe XX. The deviant behaviour of the lines at 91.02 Å and 132.85 Å appears again in Sects. 4 and 5 where we discuss the temperature and density diagnostics. For completeness, these lines were included in the tables though the results based on their ratios should be considered as less reliable.
![]() |
Figure 2: Temperature dependence of photon flux ratios between Fe XVIII to Fe XXI lines. The observed ratios are indicated |
Open with DEXTER |
![]() |
+1a | -1a | av.b | Tc |
Fe IX/Fe X | ||||
171.07/174.53 | 6.36 | -- | ![]() |
![]() |
Fe XVII/Fe XVIII | ||||
15.265/16.078 | 1.07 | 0.97 | 1.01(.07) | ![]() |
16.780/15.628 | 6.36 | 5.30 | 5.7(.6) | ![]() |
Fe XVIII/Fe XVIII | ||||
15.628/93.92 | 0.25 | 0.33 | 0.10(.01) |
![]() |
Fe XVIII/Fe XIX | ||||
103.94/108.37 | 0.63 | 0.58 | 0.61(.05) | ![]() |
103.94/101.55 | 1.9 | 2.0 | 2.0(.2) | ![]() |
93.92/91.02d | 11.5 | 10.6 | 11.6(1.6) | ![]() |
93.92/108.37 | 1.84 | 1.70 | 1.51(.15) |
![]() |
103.94/111.70 | 6.26 | 7.72 | 7.0(1.5) |
![]() |
Fe XIX/Fe XX | ||||
120.00/118.66 | 1.3 | 1.0 | 1.2(.2) |
![]() |
120.00/121.83 | 0.61 | 0.70 | 0.63(.07) | ![]() |
120.00/132.85d | 0.59 | 0.46 | 0.35(.04) |
![]() |
108.37/118.66 | 5.0 | 4.1 | 4.3(.5) | ![]() |
108.37/121.83 | 2.4 | 2.8 | 2.2(.2) | ![]() |
108.37/132.85d | 2.4 | 1.8 | 1.21(.09) | ![]() |
111.70/118.66 | 0.50 | 0.31 | 0.38(.12) | 7.0 +0.8-0.6 |
111.70/121.83 | 0.24 | 0.21 | 0.21(.06) | 6.8 +0.9-0.5 |
101.55/121.83 | 0.79 | 0.79 | 0.73(.07) | ![]() |
109.97/121.83 | 0.35 | 0.27 | 0.28(.05) | ![]() |
91.02d/121.83 | 0.39 | 0.45 | 0.31(.03) | ![]() |
Fe XX/Fe XXI | ||||
118.66/117.51 | 2.2 | 5.2 | 3.3(1.0) | 8.4 +1.0-0.7 |
121.83/117.51 | 4.6 | 7.6 | 6.5(2.0) | 8.4 +0.9-0.7 |
132.85d/117.51 | 4.8 | 11.7 | 11.9(3.9) | 7.7 +0.9-0.5 |
Fe XXI/Fe XXII | ||||
117.51/117.17 | 0.32 | 0.20 | 0.27(.08) | 9.7 +1.2-0.8 |
![]() |
+1a | -1a | av.b | Tc | Td |
C V | |||||
![]() |
sum | 0.80(.26) |
![]() |
1.0 | |
N VI | |||||
![]() |
1.44 | 1.28 | 1.44(.46) | 0.5+0.5-0.2 | 0.35 |
O VII | |||||
![]() |
0.88 | 0.90 | 0.89(.09) | ![]() |
1.9 |
Mg XI | |||||
![]() |
0.83 | 0.67 | 0.74(.11) | 4.6 +1.4-1.0 | 5.7 |
Si XIII | |||||
![]() |
0.78 | 0.94 | 0.87(.17) | 5 +3-2 | 6 |
C V/C VI | |||||
![]() |
sum | 0.59(.24) |
![]() |
||
N VI/VII | |||||
![]() |
0.213 | 0.269 | 0.27(.04) |
![]() |
|
O VII/VIII | |||||
![]() |
0.239 | 0.229 | 0.37(.02) |
![]() |
|
Mg XI/XII | |||||
![]() |
1.34 | 1.38 | 1.61(.18) | ![]() |
|
Si XIII/XIV | |||||
![]() |
2.79 | 3.33 | 2.82(.44) | ![]() |
In Fig. 3 and Table 5 theoretical ratios of temperature-sensitive lines of higher charge state ions of iron are shown and compared with the observations (Brickhouse et al. 1995). The line ratios from the lower charge state Fe ions Fe IX and Fe X indicate that these ions are formed in a plasma with a temperature of about 1 MK. For the higher Fe ions Fe XVIII, Fe XIX, and Fe XX temperatures of 6 MK or more are observed, but for ions higher than Fe XX the temperature increases to 8 or 10 MK. These results are consistent with those obtained by Brickhouse et al. (2000) on the basis of EUVE observations.
We notice some deviations including the Fe XIX lines at 108 and 120 Å and those at 101, 109, and 111 Å. The first two both decay from the same upper level (922.908 cm-1), while the three others all decay from the level at 984.760 cm-1. In the ratio between the Fe XIX/Fe XX this results in a lower temperature when the two lines at 108 and 120 Å are involved and a somewhat higher temperature when the three other lines are used. Errors induced in the computed populations of the levels in the theoretical model resulting from inaccurate collisional excitation rates are likely to be the underlying reason. We notice from Table 5 a fairly regular increase of the temperature from 6 MK up to 10 MK going from Fe XVIII to Fe XXII. This favours a multi-temperature model; this comes as no surprise based on the earlier analyses using EUVE and ASCA observations that arrived at similar conclusions (e.g., Dupree et al. 1993; Brickhouse et al. 2000).
The wide dynamic wavelength range of the LETGS allows to study
simultaneously the
and
iron line transitions
in the L-shell and to compare the derived temperatures. In
the short wavelength range, around 16 Å the observed features
correspond to
transitions,
while in the wavelength range above 80 Å the
lines originate from
transitions.
From the measured photon flux ratio Fe XVII 15.265Å/Fe XVIII
we derive a temperature of 7.0
0.2 MK using
MEKAL, but after correcting (enhancing) the Fe XVIII flux by a factor 2.14
as indicated by the results from a benchmark study of the MEKAL code
with solar flare spectra by Phillips et al. (1999; their Table A8).
The two Fe XVII/Fe XVIII
ratios result in higher temperatures than the
ratios based on
transitions. As a test case the ratio Fe
XVIII/Fe XVIII between an
and an
transition was calculated. It results in an even higher temperature
of 10 MK. However, this ratio covers a large wavelength range and is
more sensitive to flux calibration errors.
Higher temperature regions (T = 10 MK) are present in the coronae of Capella (e.g., Brickhouse et al. 2000), though with emission measures more than an order of magnitude below those of the dominant temperature of several MK. Lines of higher Fe ionization stages (e.g., Fe XXIV) that could constrain higher temperatures are outside the range in which the LETGS has either insufficient spectral resolution or sensitivity to yield accurate measurents of their fluxes.
The LETGS spectra of Capella contain a number of interesting emission
line features in the short wavelength region between 6 and 45 Å that
originate from the He-like ions Si XIII, Mg XI, Ne IX, O VII, N VI, and C V.
It concerns the resonance line 1sS0 - 1s2p 1P1(r),
the intercombination line 1s
S0 - 1s2p 3P1,2(i)and the forbidden line 1s
S0 - 1s2s 3S1(f). There
are combinations of these
three transitions that are strongly temperature- and density-sensitive. The
resonance line intensity is comparable to the sum of the intensities
of the two other lines and increases at higher temperatures, while for
high density the intercombination line becomes intense compared to the
forbidden line. Of the six ions mentioned above,
the Ne IX triplet could not be used because of severe blending
of the r and i lines by Fe XIX lines and partly blending of the f line by an Fe XVII line.
For C V the i and f lines are rather weak, but for N VI and O VII three clearly separated components
were measured. For Mg XI and Si XIII the intercombination line was
blended by the resonance line, but a
double Gaussian fit enabled us to deconvolve these lines, resulting in
measured wavelengths in very good agreement (within about 0.005-0.01
Å) with the
theoretical values. In the case of Si, the measurement
should be considered less certain due to the lower spectral resolution
at these wavelengths.
From the ratio G = (i + f)/r the temperature is determined for each He-like ion and is given in Table 6. We observe a wide range of temperatures from about 0.5 MK for N VI up to 5 MK for the higher stages of ionization (Si XIII). However, the N VI temperature appears somewhat discrepant in comparison to those derived for C and O and we suspect it might be affected by line blending. The temperatures were all derived using an updated version of the plasma radiative loss and line fitting program SPEX developed at SRON and from recent calculations by Porquet & Dubau (2000) and by Mewe et al. (2001) including all recombinations. Another possibility for obtaining the coronal temperature from He-like ions is to use the ratio between the resonance lines of these ions and the resonance lines of the one higher ionized stages, the H-like ions, of the same element. These ratios and the derived temperatures are collected also in Table 6. The single He-like ion ratios result in a lower temperature than the ratios from He-like versus H-like ions. The conclusion is that the H-like ions are formed in another (hotter) area than the He-like ions. Therefore the obtained temperature is a mixture and average of the two areas.
![]() |
+1a | -1a | av.b |
![]() |
Fe XX |
||||
![]() |
0.15 | 0.16 | 0.15(.06) |
![]() |
![]() |
0.07 | 0.11 | 0.08(.03) |
![]() |
![]() |
0.07 | 0.07 | 0.04(.02) |
![]() |
Fe XXI | ||||
![]() |
0.22 | 0.10 | 0.09(.03) |
![]() |
![]() |
0.34 | 0.20 | 0.17(.05) |
![]() |
Fe XXII | ||||
![]() |
0.14 | 0.18 | 0.17(.10) |
![]() |
![]() |
0.30 | 0.33 | 0.19(.10) |
![]() |
![]() |
+1a | -1a | av.b |
![]() |
![]() |
C V | |||||
![]() |
sum | 0.89(.6) | ![]() |
4 109 | |
N VI | |||||
![]() |
0.52 | 0.53 | 0.49(.17) | ![]() |
6 109 |
O VII | |||||
![]() |
0.26 | 0.31 | 0.28(.05) |
![]() |
3 109 |
Mg XI | |||||
![]() |
0.51 | 0.40 | 0.45(.13) |
![]() |
3 1012 |
Si XIII | |||||
![]() |
0.62 | 0.82 | 0.72(.24) | ![]() |
4 1013 |
From a number of density-sensitive iron lines at wavelengths above 90
Å (e.g., Fe XX to Fe XXII) the electron density can be investigated for
plasmas with relatively high temperatures (
MK). For the
determination of the electron density in the coronal plasma of Capella
ratios of density-dependent lines of Fe XX, Fe XXI, and Fe XXII were
used. The data used and density values obtained are listed in Table 7. The
theoretical density dependence for these lines is taken from the work
by Brickhouse et al. (1995). The theoretical flux ratio-
curves
are very flat for low densities and start to deviate from the
low-density plasma situation only for values of
cm-3.
This makes the method insensitive for the lower limit.
The
values obtained
for the ions Fe XX to Fe XXII, however, are consistently around
1012 cm-3and all with an upper limit of
-
cm-3.
Two of the three transitions of the He-like ions described above are also strongly density-dependent (e.g., Gabriel & Jordan 1969), namely the ratio of the intercombination line (i) and the forbidden line (f). Line ratios (i/f) were derived from measured values and the corresponding plasma densities were derived; the results are listed in Table 8. The uncertainty of the counting rate of the O VII intercombination line justifies an upper limit only. As for the temperature diagnostics, the obtained values are from SPEX and from Porquet & Dubau (2000) and Mewe et al. (2001).
If we compare the results for i/f and
with those obtained by
Ness et al. (2001) in the case of C, N, and O, we note that the values
are slightly different but in agreement within the error bars. Ness
et al. (2001) find for O VII a somewhat different value for i/f(
)
based on a method in which the total
spectrum is fitted with a model for the continuum plus background,
whereas we fit the background-subtracted spectrum. Because below
about 109 cm-3 the dependence of i/f on
becomes very
flat (Fig. 5 of Ness et al.) this gives a much lower upper limit for
(only 5 108 cm-3) than we find but the difference is
still acceptable within the statistical uncertainties.
From the observations of Capella with the High-Energy Transmission
Grating on CHANDRA Canizares et al. (2000) obtain a slightly larger
i/f ratio for O VII (
,
but consistent with our ratio
of
)
implying a higher electron density within the
range (8-20) 109 cm-3. Again the difference in i/f may be
due to a different treatment of the continuum plus background.
For the lower ions C V, N VI, and O VII we obtain densities typically in the
range 109-1010 cm-3 for temperatures between about 0.5
and 2 MK. Such densities are typical of those for active coronal
regions on the Sun. In the case of Capella, similar densities for the
lower temperature (1-2 MK) plasma were suggested by Schrijver et al. (1995)
and by Brickhouse (1996) in analyses of EUVE spectra.
The higher ions Mg XI and Si XIII associated with higher temperatures (
MK)
suggest much higher densities, i.e. 1012-1013 cm-3, that
are consistent with the results obtained for the highly ionized Fe
ions. Such high densities have been observed on the Sun only during
flares, but were derived earlier from EUVE spectra of Capella (cf. Sect. 1).
We note, however, that the density obtained for Si XIII is significantly higher
than the upper limit found from Mg XI.
For the same line complex observed with the CHANDRA HETG (albeit at a different
time), Canizares et al. (2000) obtained a ratio i/f of
,
which is to be compared with our value of
.
It is clear from the recent analysis of HST STIS spectra by
Johnson et al. (2001) that the hot 10 MK component of the Capella coronae
can be highly variable, and so it is therefore
possible that the higher density we derived based on observations made
at a different time is real. However, since the spectral resolution
of the LETGS is fairly low at the wavelength of the Si He-like
complex, rendering line deconvolution more difficult, this
density should perhaps be treated with caution at this time.
We performed an emission measure (EM) analysis of the Capella LETGS data by comparing measured line fluxes with the theoretical fluxes as calculated by MEKAL for a given emission measure and formation temperature for a particular line, assuming solar abundances (Anders & Grevesse 1989). Figure 4 (top) presents the results of the EM distribution based on the analysis of the iron lines. For all Fe-ions used the error bars are given. These error bars contain two components. When having many lines (and values) for one particular ion, the error is dominated by the statistical error of these values. In case of one single transition or two closely lying values for an ion, the error is dominated by the uncertainty in the flux as given in Tables 1 and 2. The larger error bar for Fe XVI is probably due to the fact that all lines used are near the gap between the detectors where the efficiency declines rapidly and the calibration is uncertain.
![]() |
Figure 3: EM distribution derived from Fe IX to Fe XXII line intensities for the Capella CHANDRA-LETGS observations (top) compared to the results from a 4-T fitting with SPEX and from observations with EUVE (Brickhouse et al. 2000) and EXOSAT (Lemen et al. 1989) (bottom) |
Open with DEXTER |
For comparison we give in the bottom panel of Fig. 3 the results from two other instruments: EXOSAT (Lemen et al. 1989) and EUVE (Brickhouse et al. 2000). These distributions are comparable to our results given the fact that the results are obtained at different times and with different instruments.
Using the SPEX code (Kaastra et al. 1996) we have fitted the Chandra Capella spectrum with a 4-temperature model with variable abundances. The result scaled to the solar abundance of iron is shown in the bottom panel of Fig. 3 and agrees satisfactorily with the emission measure distribution given on top. The derived ratios of and Ne to Fe are close to solar or mildly subsolar, while the N, Mg, and Si ratios are about a factor of 2 larger than compared to that taken from Anders & Grevesse (1989). For Mg and Si these results agree with those from Brickhouse et al. (2000) based on an analysis of EUVE and ASCA spectra, but our Ne/Fe ratio is a factor of 3 larger.
Ion |
![]() |
![]() |
Lc |
![]() |
EMe | ff | fg |
C V | 0.03 | 1.0 | 4.3 | 1.0 | 0.74 | 73 | 6.5 |
N VI | 0.05 | 1.4 | 5.1 | 1.4 | 1.10 | 34 | 3.0 |
O VII | 0.07h | 2.0 | 7.4 | 2.0 | 0.71 | 7.4 | 0.68 |
Mg XI | 40h | 6.3 | 0.13 | 6.2 | 6.6 | 1.1-2 | 1.0-3 |
Si XIII | 300i | 10 | 0.043 | 9.8 | 5.3 | 5.2-4 | 4.8-5 |
Fe XX | 30h | 8.9 | 0.34 | 8.7 | 1.2 | 1.5-3 | 1.4-4 |
Fe XXI | 30h | 10 | 0.43 | 9.8 | 0.64 | 6.3-4 | 5.8-4 |
Fe XXII | 30h | 12 | 0.62 | 12 | 0.38 | 2.5-4 | 2.3-5 |
Because the two binary components both contribute to the observed X-ray emission (see Sect. 1), any
interpretation of the composite spectrum can only be given in terms of average stellar properties.
We derive from our observations global properties of the X-ray emitting
coronal structures using a standard model of a quasi-static, semi-circular,
magnetic loop with constant cross section as developed by Rosner et al. (1978; hereafter RTV).
In the RTV loop model
the basic loop parameters are related through
a simple scaling law:
![]() |
(3) |
![]() |
(4) |
Generalizations of the basic RTV model have been developed by Vesecky et al. (1979; hereafter
VAU), including cross sections that increase with height in a way representing loop field structures by a
magnetic line dipole below the chromosphere. The scaling law is only slightly modified by the VAU results which
can be conveniently approximated by introducing into the rhs of Eqs. (3 or 4) an additional multiplication factor
or
,
respectively. Here
is the ratio of the loop cross section at the apex to that at the
loop footpoint, the expansion factor.
It turns out that a large variety of solar loops ranging from small bright points to active-region loops and
large-scale structures can be modeled by these scaling laws (e.g., discussion by Mewe 1991).
Applying Eq. (4) for the derived densities (Sect. 5) we obtain estimates for the loop
semilengths L (Table 9), taking
equal to the temperature of maximum line formation (from MEKAL).
For a given density the emission measure places an additional constraint on the emitting volume (i.e. the "filling factor'').
From the measured line fluxes and the known emissivities at temperature of maximum line formation we derive the emission measure
,
where V is the volume contributing to the emission (for cosmic abundances the hydrogen density
). In Table 9 we give EM for various lines.
Assuming for simplicity that the X-ray corona consists of one class of
N identical loops with parameters
,
we derive for the emission measure
![]() |
(5) |
![]() |
(6) |
The results appear to show the existence of two main classes of loop structures in the coronae of the binary system Capella. On the one hand, the analysis of the lower ions N VI and O VII indicates temperatures, densities, and loop lengths that are more or less representative of active coronal regions on the Sun (with relatively large filling factors).
On the other hand, the higher He-like ions
Mg XI and Si XIII and the highly ionized Fe ions (Fe XX-XXII) indicate
the presence of a class of much more compact and dense, hot loops (with very
small filling factors). On the Sun temperatures and densities this high - as observed in non-flaring active regions in the
coronae of Capella -
are seen only during flares. For both classes of loop systems
and R,
where the pressure scale height
(in cm) is given by
![]() |
(7) |
In all cases the loop sizes are much smaller than the stellar radius so that it is clear that the X-ray emission originates from compact regions that must be supported by strong magnetic fields (e.g., of the order of several hundred Gauss as was suggested by Dupree et al. 1993). Similar conclusions were drawn for other active binaries based on a review of EUVE observations by Drake (1996).
Echo the salient features of earlier EUVE analyses, we have
established different temperature regimes in the Capella coronal
plasmas: from
-2 MK for Fe IX & X, N VI, and O VII,
5-6. MK for Fe XVIII-XX up to
9 MK for Fe XXI and
Fe XXII. The upper limit to the densities obtained from Fe XX-XXII
lines is about 3 1012 cm-3. The density of the plasma
obtained from density sensitive lines of He-like triplets N VI and O
VII is in the range 109-1010 cm-3 for
0.5-2 MK and from Mg XI and Si XIII in the range
1012-1013 cm-3 for
5-9 MK.
Abundances derived for various elements are about solar (O, Ne) or twice solar (N, Mg, Si).
The results for the densities and emission measures are interpreted in terms of quasi-static loop models. They appear to show the existence of two main classes of loop structures in the coronae of Capella. The analysis of the lower ions N VI and O VII indicates temperatures, densities, and loop lengths that are more or less representative for active coronal regions on the Sun (with relatively large filling factors). On the other hand, the higher He-like ions Mg XI and Si XIII and the highly ionized Fe ions (Fe XX-XXII) indicate the presence of a class of much more compact and dense, hot loops (with very small filling factors). On the Sun temperatures and densities this high are seen only during flares. For both classes of loops the size is much smaller than the pressure scale height and the stellar radius.
Comparisons between ratios of lines decaying from the same upper level and branching ratios give confidence in
the in-flight total area calibration and the detector sensitivity as well as the reliability of the calculated
transition probabilities.
Based on photon fluxes of lines, sensitive to resonance scattering, we have derived conservative upper limits to
the optical depth of the lines Fe XVIII 93.92 and Fe XVII
15.014 of 1 and 1.5, respectively.
Acknowledgements
The Space Research Organization Netherlands (SRON) is supported financially by NWO. We thank Nancy Brickhouse for her contribution to the effort to obtain for Capella a long observation time with the LETGS during the initial in-flight calibration of the CHANDRA transmission gratings in the context of the "Emission Line Project'' headed by her and Jeremy Drake. JJD was supported by NASA contract NAS8-39073 to the Chandra X-ray Center during the course of this research. Finally, we wish to thank the referee for helpful comments that allowed us to improve the paper.