A&A 365, L329-L335 (2001)
M. Audard1, E. Behar2, M. Güdel1, A. J. J. Raassen3, D. Porquet4, R. Mewe3, C. R. Foley5, and G. E. Bromage6
Send offprint request: M. Audard
1 - Paul Scherrer Institut, Würenlingen and Villigen, 5232 Villigen PSI, Switzerland
2 -
Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA
3 -
SRON Laboratory for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands
4 -
CEA/DSM/DAPNIA, Service d'Astrophysique, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France
5 -
Mullard Space Science Lab., University College London, Surrey RH5 6NT, UK
6 -
Centre for Astrophysics, University of Central Lancashire, Preston PR1 2HE, UK
Received 2 October 2000 / Accepted 14 November 2000
Abstract
We present the high-resolution RGS X-ray spectrum of the stellar binary Capella observed by the
XMM-Newton satellite. A multi-thermal approach has first been
applied to fit the data and derive elemental abundances. Using the latter, the
emission measure distribution has been reconstructed using a Chebychev polynomial fit.
Its shape is found to display a sharp peak around 7 MK,
consistent with previous EUVE and ASCA results. A smaller but
significant amount of emission measure is required around 1.8 MK in order to explain the
"O VII He-like triplet and the "C VI Ly
line. We have applied the temperature
diagnostics of dielectronic recombination satellite lines to the He-like "O VII triplet
to constrain the cool plasma temperature, and have obtained a lower limit consistent with
the global reconstruction of the emission measure distribution.
We have used line ratios from the forbidden, intercombination, and resonance lines of the
"O VII triplet to derive an average density for the cool coronal plasma
(
cm-3). Implications for the coronal structure of
Capella are discussed.
Key words: stars: abundances - stars: activity - stars: coronae - stars: individual: Capella - X-rays: stars
Author for correspondance: audard@astro.phys.ethz.ch
Capella (
Aurigae; 13 Aurigae; HD 34029;
HR 1708) is, at a distance of 12.93 pc (Perryman et al.
1997), one of the brightest X-ray
objects visible in the sky. It is composed of G1 III + G8 III
star (Strassmeier & Fekel 1990; Hummel et al.
1994) separated by 56.47 mas (Hummel et al. 1994). The orbital
period
(P=104 d) is not linked to the rotation period of each component,
the G1 giant completing about 12 revolutions in one orbital period (Hummel et al.
1994). Catura et al. (1975)
first detected weak X-ray emission from Capella, quickly confirmed by Mewe
et al. (1975). Subsequent extreme ultraviolet (EUV) and X-ray observations
have been prolific with most satellites (e.g., Cash et al.
1978; Holt et al. 1979; Mewe et al.
1982; Vedder & Canizares 1983; Lemen et al.
1989; Dupree et al. 1993;
Schrijver et al. 1995; Favata et al.
1997; Brickhouse et al. 2000; Behar et al. 2000;
Brinkman et al. 2000; Canizares
et al. 2000), however still leaving unresolved problems
in the interpretation of Capella's coronal spectrum.
Linsky et al. (1998) found that the contribution of both stars
to the total flux of the coronal "Fe XXI line in the ultraviolet regime
was similar.
Dupree et al. (1993) showed that iron
("Fe XV-XXIV) was dominating the
EUVE spectrum. Their emission measure (EM)
distribution ranged from 0.1 to 63 MK, with minimum EM
around 1 MK and a sharp peak around 6 MK. Based on lines of highly
ionised "Fe XXI, they derived an electron density of
-1013 cm-3. However, their EM distribution was not in
agreement with later BeppoSAX results (Favata et al. 1997). Brickhouse et al.
(2000) studied simultaneous EUVE and
ASCA observations of Capella. They found that the low
first-ionization-potential (FIP) elements Mg, Si, S, and Fe have
coronal abundances consistent with solar photospheric values, while
the high-FIP element Ne appears to be underabundant by a factor of
3 to 4. However, they were not able to constrain the O abundance, while Brickhouse
(1996) derived a subsolar O abundance from the EUVE data.
Dupree & Brickhouse (1996) found a long-term variability in the EUV intensities
of "Fe XXI to "Fe XXIV by up to a factor of 4.
Recently, first results on Capella were obtained from Chandra HETG/LETG
(Behar et al. 2000; Brinkman et al. 2000;
Canizares et al. 2000; Mewe et al. 2001b;
Ness et al. 2001), confirming the dominance of highly ionised
Fe lines in the X-ray spectrum of Capella. Density diagnostics
applied to the "C V, "N VI and "O VII triplets for LETG data implied a low density
regime (
,
and <5 times 109 cm-3, respectively).
In HETG data, the "O VII triplet, formed at low T, gave a slightly higher density
(0.8-
cm-3), while the
"Mg XI and "Si XIII triplets, formed at higher T, gave upper limits near
and
cm-3. Little (or no) evidence for opacity effects
in the 15.014 Å "Fe XVII line has been seen.
Behar et al. (2000) reproduced fairly well the HETG Fe L-shell spectrum by assuming a single
electron temperature of 0.6 keV. The present paper presents first
results from the observation of Capella with XMM-Newton.
Capella was observed several times by XMM-Newton (Jansen et al. 2001) for calibration purposes. In this Letter, we present the Reflection Grating Spectrometer (RGS; den Herder et al. 2001) data of the on-axis observation (2000-03-25, 11:36:59 UT until 2000-03-26, 02:53:49 UT) which had a low instrumental background. The effective exposure times were 52.3 and 52.4 ksec for RGS1 and RGS2, respectively. The data from the other X-ray instruments onboard XMM-Newton could not be used, because they were severely piled-up and optically contaminated by the brightness of Capella.
The data were analysed with the official ESA XMM Science Analysis System (SAS) software,
version 4.1, and an update of several RGS tasks together with the latest calibration files available at the time of the
analysis. The metatask RGSPROC 0.77 was used to process the RGS data. Spectra were
extracted along the dispersion direction using a spatial mask together with a cut in the plane
of dispersion angle vs. CCD energy. The satellite pointing was stable, except for a short (600 s)
deviation of the attitude pointing that had no significant influence on the RGS spectra of Capella.
The RGS response matrices were generated by RGSRMFGEN 0.29. Although the detailed
description of the response of both spectrometers is expected to
evolve over time, we should note that the description of the RGS2 response is more advanced
than for RGS1. Some systematic errors may be introduced for RGS1 results. However, any future
RGS analysis with a more advanced response is not expected to give results that deviate more than
10-25 % from the current results.
A number of individual line fluxes have been measured in the RGS spectra (Table 1).
Ion | ![]() |
RGS | HETG </I>a | LETG </I>b |
"Mg XII | 8.421 |
![]() |
1.5 |
![]() |
"Mg XI | 9.169 |
![]() |
3.5 |
![]() |
"Mg XI | 9.231 |
![]() |
0.6 |
![]() |
"Mg XI | 9.314 |
![]() |
1.9 |
![]() |
"Fe XVII | 15.014 |
![]() |
30.4 |
![]() |
"Fe XVII | 16.775 |
![]() |
20.0 |
![]() |
"Fe XVII | 17.051 |
![]() |
26.4 |
![]() |
"Fe XVII | 17.100 |
![]() |
24.4 |
![]() |
"O VIII | 18.969 |
![]() |
26.3 |
![]() |
"O VII | 21.602 |
![]() |
9.7 |
![]() |
"O VII | 21.804 |
![]() |
2.6 |
![]() |
"O VII | 22.101 |
![]() |
7.4 |
![]() |
"N VII | 24.781 |
![]() |
5.5 |
![]() |
"C VI | 33.734 |
![]() |
![]() |
![]() |
Figure 1: RGS1 data with the overlaid best-fit model (red thick line) derived from a reconstructed emission measure distribution (Fig. 2). The labels identify major lines. Note the different scales in each panel |
Open with DEXTER |
In this Letter, we have used a "global fitting'' approach.
This approach allows us to obtain a self-consistent solution that best fits the data and
that takes into account simultaneously the contributions of the continuum components
and of all emission lines. Line blends and the overlap of line wings are therefore fully accounted for, and
the line fluxes are correctly reproduced within the accuracy of the code. However, as will be
discussed in the following
sections, the fit depends on the completeness of the atomic database used by the spectral fitting code.
We have used the publicly available Utrecht software SPEX 2.0 (Kaastra et al.
1996). It contains a collisional ionization equilibrium (CIE) model that is
equivalent to the MEKAL (Mewe et al. 1995) code available in the XSPEC software (Arnaud
1996), except that the former includes a significant update of the wavelengths of Fe L-shell
lines between 10-18 Å, based on the solar data by Phillips et al. (1999).
A spectral model using three CIE components with a column density
cm-2
(Linsky et al. 1993)
allowed us to obtain an initial representation of
the temperature structure in Capella. This model was used to derive coronal elemental abundances
from a variety of bright and weak lines. We added further components (up to ten) in the multi-T
approach to investigate whether they are required for a better description of the fit parameters.
This was not the case. We therefore interpret our abundances as sufficiently accurate within the
limitations of the spectral code. The best-fit results are given in Table 2 for RGS2
data, and for completeness also for RGS1.
Parameter | RGS1 | RGS2 |
![]() |
=18.255 | =18.255 |
kT1 [keV] | 0.159+0.006-0.006 | 0.13+0.01-0.01 |
kT2 [keV] | 0.593+0.002-0.002 | 0.58+0.003-0.004 |
kT3 [keV] | 1.14+0.07-0.07 | 0.93+0.07-0.07 |
![]() |
51.90+0.03-0.02 | 51.65+0.08-0.09 |
![]() |
52.86+0.009-0.01 | 52.98+0.01-0.01 |
![]() |
51.81+0.05-0.05 | 51.85+0.08-0.09 |
C | 0.50+0.04-0.04 | 0.52+0.04-0.04 |
N | 0.82+0.07-0.07 | 1.01+0.08-0.08 |
O | 0.44+0.02-0.01 | 0.39+0.01-0.01 |
Ne | 0.51+0.03-0.03 | 0.34+0.02-0.02 |
Mg | 1.13+0.06-0.05 | 0.74+0.04-0.04 |
Si | 0.41+0.06-0.06 | 0.48+0.06-0.06 |
S | 0.14+0.02-0.02 | 0.09+0.02-0.02 |
Ar | 0.23+0.1-0.1 | 0.17+0.14-0.14 |
Ca | 0.27+0.08-0.08 | 0.27+0.11-0.11 |
Fe | 0.62+0.02-0.01 | 0.51+0.01-0.01 |
Ni | 0.88+0.07-0.07 | 0.54+0.05-0.05 |
We estimate, in the following, the influence of the fitting discrepancies on our results.
The most important discrepancies occur for the "Si XIII
lines at 6.6 Å, for the "Fe XVII-XX lines at 11.0-11.4 Å,
at 12.8 Å, at 13.5 Å, and for the "Fe XVIII lines around 16 Å. The spectral code
underestimates the flux around 9.6-10.6 Å, which has been interpreted by Brickhouse et al.
(2000) as being due to missing high excitation (n>5) lines of "Fe XVII-XIX.
Furthermore, the energy dependency of the collision strengths of Fe L-shell lines needs
to be updated in SPEX/MEKAL. Additionally, we note that the spectral code fails to pick up some of the weak lines observed
in the long-wavelength part of the spectrum. HULLAC (Bar-Shalom et al. 1998)
calculations indicate that these should mostly be attributed to L-shell
emission from Si, S, Ar, and Ca. The mentioned discrepancies are related to the incompleteness
of the atomic database of the CIE model in SPEX, and similarly of the MEKAL model in XSPEC.
We conclude that L-shell lines from several elements (e.g., S, Si, Ca, Ar, Ni) and high excitation Fe
L-shell lines are insufficiently described in these codes.
To test the robustness of the derived elemental abundances, we have
iteratively eliminated emission lines and parts of the continuum
that showed poor fit results, thus moving from global fitting towards a "single line analysis''
approach. As an aside, we note here that the overlapping line wings as well as several line blends in the
present RGS spectrum introduce considerable uncertainty if single line fluxes in the
Fe L-shell region are measured without modeling; we therefore kept the "global'' approach
for this test, even when eventually only a few bright line lines contributed
to the results. Despite the large reduction of the spectral information, the
abundances turned out to be quite robust. In all test runs,
the Fe abundance was confined to within 0.50-0.68 (times the solar photospheric value),
O within 0.24-0.39, and Ne within 0.75-1.05.
The S abundance is basically derived from weak L-shell lines in the RGS band. The determination of its
abundance is mostly influenced by the "S XIV line at 24.2 Å in the model
that is not as strong in the data. When the data around this line are removed, the abundance increases to
,
and the model better fits the remaining weak S lines.
Similarly, by removing L-shell Si ions (in the long wavelength band),
the bright Si He-like triplet is correctly fitted, with an abundance of 0.7-1.0 times the solar
photospheric value.
Using coronal elemental abundances derived from the 3-T model, the EM distribution has been
self-consistently reconstructed. Figure 2 shows one realisation of the EM distribution
using a Chebychev polynomial of order 5.
![]() |
Figure 2: Realization of the EM distribution using Chebychev polynomials of order 5. The distribution does not show any significant EM above 10 MK |
Open with DEXTER |
Using the line ratios R=f/i and G = (i+f)/r (see, e.g., Gabriel & Jordan 1969;
Pradhan 1982) from the fluxes of the resonance line (r), the
intercombination line (i) and the forbidden line (f) of the "O VII and,
tentatively, "Mg XI triplets,
we have derived average coronal densities. The Si triplet is strongly blended because
of the decreasing spectral resolution of the RGS at short wavelengths, while the
Ne triplet is heavily blended by Fe and Ni lines. For the weak "N VI triplet, the line intensities
were not unambiguously determined. Note that due to the spatially unresolved nature
of stellar coronal X-ray emission, any derived value of a density should be considered
as a weighted average of the densities in the various regions of both
coronae of Capella.
We used theoretical calculations (Mewe et al. 2001a) that take into account
the radiative and dielectronic recombinations, and the electronic collisional excitations (see Porquet
& Dubau 2000 for the atomic data); they also take into account the influence of the
radiation field (photo-excitation) which is important for low-Z ions.
We refer to Ness et al. (2001) for additional details.
From the measured RGS line fluxes (Table 1), the ratios for "O VII are
and
,
implying an average electron temperature
of
MK and an upper limit for the average density of
cm-3.
For the Mg triplet, we tentatively get
and
,
leading
to
MK and
cm-3. However, due to
the low spectral resolution at short wavelengths in the RGS (e.g., compared to HETG), the
Mg triplet line fluxes are difficult to measure, therefore the derived density
should be taken with caution.
The "O VII spectral region has been separately investigated for the presence or absence of
dielectronic recombination (DR) satellite lines. DR satellite lines of He-like spectral lines in
hot collisional plasmas are excellent additional tracers of the cool (<1 MK) plasma
otherwise not sufficiently constrained by the available lines. It is, to our knowledge, the first
time that such diagnostics is applied to the X-ray spectrum of a stellar corona.
The SPEX (and MEKAL) spectral code does presently not include a sufficient
description of the DR satellite lines for low-Z ions such O. Therefore, the lower temperature
component of the EM
distribution is not well constrained below 1 MK. Using the HULLAC code, we have calculated the emitted spectrum
of the He-like lines of "O VII
including the "O VI DR satellite lines in the low-density limit. All of the
1snl and 1s2
(
= 2 to 4) levels are included in the computations. The
resulting theoretical O spectra, as a function of electron temperature, are depicted in Fig. 3
together with the data.
![]() |
Figure 3: Observed "O VII triplet with calculated spectra of the "O VII and "O VI DR satellite lines for six different electron temperatures. The plots are normalized to the strongest line in each spectrum. Note the steep temperature sensitivity of the DR satellites (see text) |
Open with DEXTER |
Acknowledgements
MA acknowledges support from the Swiss National Science Foundation (grants 2100-049343 and 2000-058827), from the Swiss Academy of Sciences and from the Swiss Commission for Space Research. The Space Research Organization Netherlands (SRON) is supported financially by NWO. CRF acknowledges financial support from the UK Particle Physics and Astronomy Research Council.