J. Sanz-Forcada1 - F. Favata2 - G. Micela1
1 - INAF - Osservatorio Astronomico di Palermo
G. S. Vaiana, Piazza del Parlamento, 1, Palermo, 90134, Italy
2 - Astrophysics Division - Research and Science Support Department
of ESA, ESTEC, Postbus 299, 2200 AG Noordwijk, The Netherlands
Received 7 October 2003 / Accepted 6 November 2003
Abstract
We report a detailed analysis of the coronal abundance of 4 stars
with varying levels of activity and with accurately known
photospheric abundances. The coronal abundance is determined here
using a line flux analysis and a full determination of the
differential emission measure. Photospheric abundance values are
taken from the literature. Previous coronal abundance
determinations have generally been compared to solar photospheric
abundances; from this a number of general properties have been
inferred, such as the presence of a coronal metal depletion with an
inverse First Ionization Potential correlated with activity level.
Here we show that, when coronal abundances are
compared with real photospheric values for the individual stars, the
resulting pattern can be very different. Also, we present evidence
that, in some cases, the coronal metal abundance may not be uniform
in the corona; in particular it can vary with the
temperature of the emitting plasma.
Key words: stars: coronae - stars: abundances - stars: late-type - X-rays: stars - line: identification
Following the initial results from the ASCA-SIS instrument (e.g. White et al. 1994; Drake 1996), which showed that the low-resolution X-ray spectra of active coronal sources were often best fit with a low metal abundance (typically about 0.3 the solar photospheric abundance), the issue of the metal abundance of the coronal plasma has become one of the main areas of interest for X-ray spectroscopy of active stars. Due to the limited spectral resolution of the ASCA spectra, the determination of abundance for elements other than Fe was often not accurate (Favata & Micela 2003). The availability of the XMM-Newton and Chandra high-resolution X-ray spectra of coronal sources has allowed observing individual, resolved spectral lines of a number of elements in the corona of active stars, and thus to address the issue of the abundance of a number of elements. Unfortunately up to now the results have not been as conclusive as initially expected. One reason is that often the same data analyzed with different approaches result in significantly different values for the element abundances (see Favata & Micela 2003, for a number of examples).
One of the issues most often addressed by coronal abundance studies is
whether a bias of the First Ionization Potential (FIP) of the
element is present. Such an "FIP effect'' is found in the solar
corona, where the abundance of elements with low FIP (<10 eV; e.g.
Mg, Si, Fe) is observed to be enhanced with respect to
their photospheric values, while elements with high FIP (
10 eV;
e.g. O, Ne, Ar) remain at abundances similar to the photospheric
values. Such a behavior is observed in the upper solar atmosphere at
temperatures higher than
1 MK, both in the non-active regions
and in full-disk spectra (Laming et al. 1995; Feldman & Laming 2000). EUVE observations
showed the same FIP bias to be present in stars with low activity
levels, such as
Cen AB (G2V+K1V, Drake et al. 1997), and
possibly (only limited knowledge of their photospheric abundance was
available at that time)
Boo A (G8V, Laming & Drake 1999) and
Eri (K2V, Laming et al. 1996). EUVE spectra of more active
stars did not show any FIP effect, with however a rather low coronal
Fe abundance with respect to the solar photosphere, which suggested
some metal depletion in the corona of the active stars. Feldman & Laming (2000)
review some of the results on solar and stellar coronal abundances
previous to Chandra and XMM-Newton.
Table 1: Observation log.
The new generation of X-ray satellites XMM-Newton and
Chandra offers access to high resolution spectra in a
wavelength region where numerous spectral lines from different
elements are formed at coronal temperatures. This permits one to measure
the fluxes of individual lines, in principle allowing a much better
determination of the coronal thermal structure and abundances.
Initial results on coronal abundances determined from data from the
two observatories have been reviewed by Drake (2002) and
Favata & Micela (2003). These new results seem to point towards a scenario with
stars with low activity levels showing, when their coronal and
photospheric abundances are compared, a solar-like FIP effect
(case of
Cen, Raassen et al. 2003) or no FIP bias (for
Procyon, Raassen et al. 2002). Active stars (like HR 1099, II Peg, AR Lac or
AB Dor), with a quite different coronal thermal structure, dominated
by material emitting at temperatures of
10 MK, show a pattern
characterized by a general depletion of elements with low FIP (when
they are compared with solar photospheric values) while the elements
with high FIP remain at solar (or sometimes higher) photospheric
values (Huenemoerder et al. 2003; Sanz-Forcada et al. 2003b; Brinkman et al. 2001; Huenemoerder et al. 2001; Drake et al. 2001) although photospheric
abundances of elements other than Fe are not known for HR 1099.
Audard et al. (2003) analyzed XMM-Newton observations of 5 active
binary systems, concluding that low-FIP elements are subject to greater
depletion in more active stars, while high-FIP elements remain at
similar levels (although photospheric abundances were available only
for one of the objects in their sample, namely
And). In
general, the comparison of stellar coronal abundances with solar
photospheric values is of little use, given the wide range of
abundances found among late type stars (e.g. Cayrel de Strobel et al. 2001). Most of
the stars for which coronal abundances have been determined have no
published photospheric abundances, so that a true test of
FIP-related effects cannot be done. This is especially true for active
stars, for which the large rotational velocity complicates the
measurements of photospheric abundances. In most cases the available
photospheric abundances are limited to Fe (e.g. Randich et al. 1993), and
even these in many cases are likely to be affected by systematic
errors. Moreover, different methods to determine the distribution of
material with temperature in the corona (the so called "emission
measure distribution'', EMD) sometimes produce contradictory results,
depending on the analysis technique employed (Favata & Micela 2003; Sanz-Forcada et al. 2003b.
In the present work we determine the coronal thermal structure and
abundances of four stars with different activity levels (in terms of
)
observed with relatively high statistics with
Chandra and XMM-Newton (see Table 1), and for
which there are accurate photospheric abundance measurements
available, to explore the relation between activity
levels and coronal relative abundance patterns. The chosen stars are
Procyon (one of the stars with lowest
values), the intermediate activity star
Eri, the
long-period active binary
And and the active binary
V851 Cen. These four objects are among the few stars for which both
high-resolution coronal spectra with good statistics and measurements
of photospheric abundances for several elements are available.
The observations are described in Sect. 2, a brief description
of the analysis techniques is presented in Sect. 3.
The results for each star are discussed in
Sect. 4, followed by the discussion and conclusions in
Sect. 5.
![]() |
Figure 1:
Light curves of Procyon (LETGS), |
| Open with DEXTER | |
We have analyzed observations conducted using different X-ray
high-resolution spectrographs; however, the spectral lines observed in
the different spectra are formed in very similar temperature ranges
(
-7.8), hence yielding information of the
same region of the corona (see Sect. 3). The X-ray
luminosity of the targets has been calculated by modeling the observed
spectra, as listed in Table 1.
An XMM-Newton observation of V851 Cen (P.I.
F. Favata) was executed in March 2003
(Table 1). XMM-Newton carries out simultaneous
observation with the EPIC (European Imaging Photon Camera) PN and MOS
detectors (sensitivity range 0.15-15 keV and 0.2-10 keV
respectively), and with the RGS (Reflection Grating
Spectrometer, den Herder et al. 2001) (
-38 Å,
/
-500). The RGS data have been
reduced employing the standard SAS (Science Analysis Software) version 5.4.1 package, removing in the RGS spectra the time intervals when
the background was higher than 0.4 cts/s in CCD #9, to
ensure a "clean'' spectrum. Light curves were obtained by selecting
a circle centered on the source in the EPIC-pn and EPIC-MOS images,
and subtracting the background count rate taken from a nearby area.
The light curve shows a slow decline (
20%) of the source flux
(Fig. 1). High resolution spectra corresponding to the
first and second orders of the RGS have been used to measure the line
fluxes simultaneously, improving the quality of the measurements.
Fluxes measured for second order lines are in good agreement with
those in the first order (Sanz-Forcada et al. 2003b), hence they can be used to
improve the quality of the measurements due to their better spectral
resolution.
An observation with the Chandra High Energy Transmission
Grating Spectrograph (HETG, Weisskopf et al. 2002) of
And was
retrieved from the Chandra archive (Table 1). The
HETGS is made of two gratings, HEG (High Energy Grating,
-15,
-1070),
and MEG (Medium Energy Grating
-27,
-970) operating simultaneously.
Standard reduction tasks from the CIAO v2.3 package have been employed
in the reduction of data and the extraction of the HEG and MEG
spectra. The light curve of
And, obtained from the MEG and HEG first order spectra, shows no significant events, with variations
of less than 20% of the total flux (Fig. 1).
Observations with the Chandra Low Energy Transmission Grating
Spectrograph (LETGS,
-175,
-1000 Weisskopf et al. 2002) and the High Resolution Camera
(HRC-S) of Procyon and
Eri were retrieved from the
Chandra archive as listed in Table 1. Data were
reduced using standard reduction tasks present in the CIAO v2.3
package. The positive and negative spectral orders were summed for
the flux measurements, although at long wavelengths (
40 Å)
separate measurements were made in the plus and minus spectra due to
problems with wavelength calibration or detector gaps. Lines formed in
the first dispersion order but contaminated with contribution from
higher dispersion orders were not employed in the analysis. The two
observations available for Procyon were summed to achieve
better statistics, resulting in an exposure time of 139 ks.
Light curves were obtained from the LEG spectra (first and higher
orders) of Procyon and
Eri as shown in Fig. 1.
![]() |
Figure 2:
Element abundances in the corona of Procyon,
|
| Open with DEXTER | |
![]() |
Figure 3:
Coronal to photospheric abundance ratio for Procyon,
|
| Open with DEXTER | |
X-ray spectra from stellar coronae have numerous lines from different
elements formed at different temperatures. This variety of
temperatures in the lines formation makes it necessary to determine in
detail the thermal structure of the target star, commonly
parameterized through the emission measure distribution (EMD) as a
function of temperature, defined as
[cm-3]. In the present work we apply a line-based method
to reconstruct the EMD of the stellar coronae (this method
is preferred to the global fitting techniques, see Sanz-Forcada et al. 2003b), with a
step in temperature of 0.1 dex, the same as employed in the
Astrophysical Plasma Emission Database (APED v1.3, Smith et al. 2001),
which is used for comparison of the observed with the predicted line
fluxes as explained below.
The EMD reconstruction was conducted by measuring the fluxes from
spectral lines present in the high resolution X-ray spectra, following
Sanz-Forcada et al. (2003b). The Interactive Spectral Interpretation System
(ISIS, Houck & Denicola 2000) software package, provided by the MIT/CXC, was
employed to measure the line fluxes, through convolution of the
spectral response. These measured line fluxes were then corrected for
interstellar medium absorption (ISM) (although such correction is
important only for very few lines).
Theoretical fluxes can be predicted by assuming an initial EMD which
is combined with the emissivity functions of the different spectral
lines. These fluxes are then compared with the observed fluxes, and
the EMD is changed to determine a solution that shows a better
agreement between predicted and observed fluxes. An iterative process,
as explained in Sanz-Forcada et al. (2003b), results in the "best'' solution for
both the EMD and the abundance pattern. Although this solution is not
unique, the use of a large number of spectral lines with a good
temperature coverage constrains the interval of possible solutions to
a small range of values around the solution proposed. Error bars for
the EMD can be calculated with a Monte Carlo method, consisting
of the variation of the observed fluxes by ![]()
and the
determination of the best solution among 1000 possible variations
around the initial EMD solution. Such an exercise is performed for 1000 sets
of possible line fluxes, providing error bars for the solution. The
abundances of the different elements are determined during the
calculation of the EMD, with error bars resulting from the quadratic
summation of the errors on the measurement of observed fluxes (
,
and the dispersion of the
predicted-to-observed line ratios [
]
of lines of the same element (the latter permits us to evaluate
indirectly the errors related to atomic models). The abundances of
the elements are determined relative to Fe, with the [Fe/H]
abundance derived afterwards from the best fit to the global spectrum
for the given EMD and abundances pattern. The coronal abundances
determined, with their photospheric counterparts are shown in
Figs. 2 and 3, and are listed in
Table 2. EMD distributions (Table 3)
are shown in Fig. 4. The line fluxes measured in the four
targets are listed in
Tables 4-7
.
![]() |
Figure 4:
Emission Measure Distribution for Procyon, |
| Open with DEXTER | |
Procyon (
CMi A, HR 2943, HD 61421) is a well-studied nearby
F4IV star (3.497 pc Perryman et al. 1997). A rotational period of
9.1 d can be deduced from
km s-1 (de Medeiros & Mayor 1999)
and i=32
(Irwin et al. 1992). Photospheric abundances of different
elements have been determined by many authors, resulting in a
metallicity close to solar. In this work we will use the values
calculated by Edvardsson et al. (1993), with C abundance from Varenne & Monier (1999).
First studies on the coronal abundances were carried out using EUVE
observations: Drake et al. (1995) calculated, using a line-based analysis, the EMD for low and high FIP elements separately, but they did not find
any "FIP'' or "inverse FIP'' effect. Global fits to
Chandra/LETG and XMM-Newton high-resolution spectra were
employed by Raassen et al. (2002) to determine the abundances of Procyon, also
concluding that no FIP-related effects are present in the corona of
Procyon (although no comparison with photospheric values is present in
their study).
In the present work line fluxes were corrected from the ISM assuming a
value of
(Sanz-Forcada et al. 2003a).
Several problems affect the analysis of the Chandra Procyon
spectra. First, the shape of the continuum of a relatively cool
corona like that of Procyon, in the wavelength range covered here,
does not allow us to establish the [Fe/H] coronal abundance accurately (a
wide range of values show similar line-to-continuum ratios), thus we
fixed it at its photospheric value (note that the rest of coronal
abundances are calculated relative to Fe, and the general level of the EMD depends on [Fe/H] too). Also, the determination of coronal
abundances is affected, in the LETG wavelength range, by the lack of a
complete temperature coverage provided by Fe (only two lines are
formed at
), while many lines from other elements
have little temperature overlap with Fe lines. However it is possible
to include in the analysis the Fe lines observed with EUVE
(Fe XI- XVI), obtaining a better temperature
coverage for one element, useful for the determination of abundances
of other elements. Measurements of EUVE Fe line fluxes in the range
170-370 Å were made by Sanz-Forcada et al. (2003a) from observations taken in
1993, 1994 and 1999 (the latter is partially simultaneous with the
LETG campaign, and no substantial changes were detected in the EUV
flux level). The inclusion of the EUVE lines allows us to cover the range
-6.8, where almost all lines of other elements are
also formed.
Table 2: Coronal and photospheric abundances of the elements ([X/H], solar units) in the target stars.
The abundances determined for Procyon do not show substantial
differences between photospheric and coronal relative abundances (note
that no error bars are published for the photospheric abundances),
with the most notable difference being a relative underabundance of O
in the corona. Nevertheless, an interesting feature arises when
fitting the EMD and abundances. Assuming that the EMD (see
Fig. 4) is correct for all elements, higher (relative to Fe) abundance of S, O, Ne, Si and Mg at lower temperatures would
improve the fit of their respective lines (the case of the S VII lines is very remarkable), while elements like Ni, Ca,
Al, C or N cover a gradient of temperatures too small to reveal any
trend. This effect could be also caused by the Fe
abundance increasing
with temperature, from the minimum of the EMD at
to
the maximum at
.
This might indicate a trend of a
mild FIP effect which increases with temperature in the corona of
Procyon (Drake et al. 1995 also discuss a possible temperature
dependence of the coronal abundances in Procyon).
In the solar corona, the FIP effect is indeed larger at
than at
(Feldman & Laming 2000).
Although suggestive, this hypothesis cannot presently be confirmed
with high confidence for Procyon given the present uncertainties in
atomic models and in the determination of the EMD of this star.
Similarly, an increase in the [Ne/Fe] abundance with temperature
between
and
would explain
the discrepancies found in the determination of the coronal Ne
abundance of AD Leo from Ne VIII lines (Sanz-Forcada & Micela 2002) and
from Ne IX and Ne X lines (Maggio et al. 2002).
As one of the nearest cool stars (d = 3.22 pc),
Eri (HD 22049, HR 1084) is frequently studied for its activity,
abundances, and more recently, stellar planetary companions (see
e.g., Hatzes et al. 2000). With
d (Baliunas et al. 1983),
Eri (K2V) has an intermediate activity level.
Photospheric abundances derived by Zhao et al. (2002) are used here for
comparison with the coronal values, and a value of
(Sanz-Forcada et al. 2003a) is adopted to correct for
ISM absorption. The first high-resolution analysis on coronal
abundances of
Eri was made by Laming et al. (1996) who analyzed
EUVE spectra, calculating separate line-based EMDs for low- and
high-FIP elements, concluding that a solar-like FIP effect was present
in the corona of
Eri. Coronal electron density diagnostics
were reported by Ness et al. (2002) from the Chandra/LETG
observations of
Eri.
In the present work we derive the coronal abundance of 11 elements
from the Chandra/LETG spectrum of
Eri. The corona is
dominated by emitting material in the range
-7.2,
with a distribution compatible with the presence of a peak at
as EUVE data suggest (Sanz-Forcada et al. 2003a), but also with two
peaks at
and 6.8, as the present data, also
compatible with EUVE fluxes, suggest. The [Fe/H] coronal abundance is
better constrained than for Procyon, with a 0.2 dex uncertainty; the
agreement between predicted and observed fluxes is also better than
for Procyon. A mild FIP-like bias is detected in the corona of
Eri, affecting only Ca and Ni (Fig. 3), with
a progressive lower coronal abundance for elements with higher FIP -
while in the Sun the FIP effect consists of a enhancement by more
than a factor of 3 of the elements with
eV.
And (HD 222107, HR 8961) is an active binary system with a
G8 IV-III primary and an unseen companion in an orbit with
d (Walker 1944).
And is one of the few
active stars for which measurements of photospheric abundances for a
good number of elements (Donati et al. 1995; Savanov & Berdyugina 1994) are available, also thanks
to its low rotation rate (
d, Perryman et al. 1997).
Coronal abundances from high resolution spectra were not obtained from
EUVE data because the ISM absorption did not allow the measurement of fluxes
of lines for elements other than Fe , see][]paperland, and the
global fit of the XMM/RGS spectra gave no clear trend with FIP
(Audard et al. 2003). A preliminary analysis of the Chandra/HETG and
LETG spectra was reported by Dupree et al. (2000) and Brickhouse et al. (2003)
respectively.
Table 3: Emission Measure Distribution of the target stars.
The low ISM absorption
= 18.45, ][]wood96 makes the correction of fluxes unnecessary
(affecting the line fluxes at less than the 1% level). The high
spectral resolution and good statistics achieved in the observation
allow us to obtain the best determination of the element abundances of
the four stars, also supported by a good agreement between
observed and predicted fluxes for most lines. The fluxes measured in
the EUV range (Sanz-Forcada et al. 2001) are also well predicted using the EMD
determined here. The coronal abundances determined here are remarkably
consistent with the photospheric values calculated by Donati et al. (1995),
and they are also close to the values calculated by Savanov & Berdyugina (1994).
V851 Cen (HD 119285) is an active binary system with a K2IV-III
primary and an unseen companion, with a photometric period of 12.05 d
(Lloyd Evans & Koen 1987; Cutispoto et al. 2001), slightly longer than its orbital period
(
d, Saar et al. 1990). V851 Cen has a small
projected rotational velocity (
km s-1, Saar et al. 1990)
which allowed Katz et al. (2003) and Morel et al. (2003) to derive accurate
photospheric abundances, revising the value of
determined by Randich et al. (1993) to a value of
.
In the
present work the coronal abundances and thermal structure of this star
is determined for the first time. A correction of the observed fluxes
due to the absorption of the ISM was applied, assuming
,
the value determined for the nearby
star
Cen (Fruscione et al. 1994)
. We use the photospheric abundances
provided by the "Method 1'' of Morel et al. (2003), as shown in
Fig. 2. When the coronal abundances are compared to solar
photospheric ones, V851 Cen shows a behavior very similar to that of
stars like AR Lac or AB Dor, with a decrease of abundance with FIP
until Si, and a progressive increment for elements with high FIP.
However, when these abundances are compared with the stellar
photospheric values, the conclusions are different, with only a small
depletion of coronal Fe, Mg and Si, and a possible enhancement of the
elements with lowest FIP (Ca and Ni), although still consistent with
the photospheric values (Fig. 3). The determination of
the photospheric abundance of elements with larger FIP would be
desirable to confirm any trend.
The results found in the present sample show that the issue of coronal
abundances in stars other than the Sun is complex and still poorly
understood. Solar coronal abundance is by itself a complex issue, with
the actual values depending on which part of the corona is
observed (with coronal holes, active regions and flares showing quite
different abundance patterns). Results obtained from full disk spectra
(thus comparable to stellar spectra) show enhancements by at least a
factor of 3 in elements with
eV (Laming et al. 1995).
Young solar active regions show coronal abundance compatible with the
photospheric values, but they soon develop a FIP bias that reaches
average coronal values in two to three days (Sheeley 1996; Widing & Feldman 2001). The
FIP effect that is observed in the solar corona seems to disappear with
increasing activity level of the star, as it has been shown in cases
like II Peg (Huenemoerder et al. 2001), AR Lac (Huenemoerder et al. 2003), AB Dor
(Sanz-Forcada et al. 2003b),
And and V851 Cen. A solar-like FIP effect
appears to be present in
Cen A (Raassen et al. 2003), a star very
similar to the Sun in many physical parameters, as well as in other
solar analogs with low levels of activity, such as
UMa and
Ori (Güdel et al. 2002), although only the Fe photospheric
abundance is known for these two objects (so that abundance patterns
assume a solar-like mix). A mid-activity level star,
Eri,
could be showing the progressive disappearance of the FIP effect, with
elements like Mg and Fe already showing photospheric abundances in the
corona, and leaving a marginal coronal abundance enhancement only for Ca and Ni. Finally, Procyon seems to be quite different to the Sun,
with no FIP-related trends evident, although there could be an
increment of Fe abundance with the coronal temperature relative to
other elements with higher FIP (see Sect. 4.1). The FIP
effect detected in the solar corona above quiet regions grows with
temperature, in a range similar to that covered here for Procyon
(Feldman & Laming 2000).
The interpretation of the metal depletion in the corona of active
stars that has been given in the past from initial results is somewhat
in contradiction with the present results. Many of the stars for
which a "metal abundance deficiency syndrome'' (MAD) has been claimed
do not have measurements of the photospheric abundance available (or
they are limited to Fe abundance determinations), and published work
is limited to the comparison of stellar coronal patterns with solar
photospheric abundances (a rather arbitrary practice). The case of
And is quite remarkable, with compatible coronal and
photospheric abundances for all elements tested, and a definitely
non-solar photospheric abundance pattern. V851 Cen, with a similar
coronal thermal structure (Fig. 4) shows a marginal coronal
abundance depletion. Other cases like II Peg (Huenemoerder et al. 2001) or AR Lac
(Huenemoerder et al. 2003) also show moderate depletion when they are compared
to the respective photospheric values, but strong metal depletion
seems to be clear only for AB Dor (Sanz-Forcada et al. 2003b, and references therein). However AB Dor is a fast rotator, and the only
available determination of photospheric abundances, carried out by
Vilhu et al. (1987), needs further confirmation given the difficulty of the
measurements in a spectrum so highly affected by rotational broadening.
While the absence of a FIP effect in active stars, at least at the solar
level, is well supported at this point, the presence of a general
"MAD effect'' or an "inverse FIP effect'' increasing with activity
is not yet clearly established, given the discrepancies found between
stars with quite similar levels of activity and thermal coronal
structure like AB Dor,
And and V851 Cen. Moreover, the
abundance pattern could also be changing in different parts of the
corona, depending on the coronal temperature, as present in Procyon.
A larger sample of stars with both photospheric and coronal
abundances, determined in a homogeneous way, would be necessary
to establish a general trend between coronal abundances and
activity levels.
Acknowledgements
We acknowledge support by the Marie Curie Fellowships Contract No. HPMD-CT-2000-00013. We have made use of data obtained through the Chandra data archive, operated by the Smithsonian Astrophysical Observatory for NASA. This research has also made use of NASA's Astrophysics Data System Abstract Service. We are grateful to the referee, Dr. J. M. Laming, for the careful reading and useful comments brought to the manuscript.