A&A 365, L344-L352 (2001)
M. Güdel 1 - M. Audard 1 - H. Magee 2 - E. Franciosini 3 - N. Grosso 4 - F. A. Cordova 5 -
R. Pallavicini
3 - R. Mewe 6
Send offprint request: M. Güdel
1 - Paul Scherrer Institut, Würenlingen & Villigen, 5232 Villigen PSI,
Switzerland
2 -
Mullard Space Science Laboratory, University
College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK
3 -
Osservatorio Astronomico di Palermo "G. S. Vaiana'', Piazza del Parlamento 1,
90134 Palermo, Italy
4 -
Max-Planck-Institut für Extraterrestrische Physik, PO Box 1312, 85741 Garching,
Germany
5 -
University of California, Santa Barbara, CA 93106, USA
6 -
SRON Laboratory for Space Research, Sorbonnelaan 2, 3584 CA
Utrecht, The Netherlands
Received 2 October 2000 / Accepted 20 October 2000
Abstract
We present results of a comprehensive X-ray investigation of the multiple
stellar system Castor AB + YY Gem. The observations were obtained with the
XMM-Newton observatory. We report the first spatially separated
detection and identification of both Castor A and B as frequently
flaring X-ray sources. Spectral modeling of the YY Gem and Castor coronal
sources indicates a broad plasma distribution between 2-15 MK.
Elemental abundances are found to be below solar photospheric values,
with the abundance of Ne being highest (0.7 times solar photospheric).
The abundances increase during a large flare, most notably
for Fe. The "O VII He-like line triplet indicates average coronal densities
of the cooler plasma of a few times 1010 cm-3 for YY Gem and Castor.
The YY Gem light curve shows three deep eclipses, indicating that the coronae
of both binary components are similarly active and are relatively compact.
A 3-D light curve inversion shows that coronal active regions follow
density scale heights compatible with spectroscopically
measured temperatures. The dominant densities also agree with
spectroscopic values. Most active regions are located between
degrees
latitude, in agreement with findings from previous Doppler imaging. We also
report the tentative detection of line broadening due to orbital motion.
Key words: stars: activity - stars: coronae - stars: flare - stars: abundances - stars: binaries: eclipsing - stars: individual: YY Gem, Castor
Author for correspondance: guedel@astro.phys.ethz.ch
Close binary systems are of great value for the study of atmospheric structuring. Although short-period binaries are often tidally locked, the mutual influence of the components on the internal structure and dynamo is poorly understood; it appears that the principal effect of binarity is the enforced rapid rotation of the components, which generates strong magnetic fields through a dynamo (Guinan & Giménez 1993).
Comparison between coronal properties of binary components and single, active field stars potentially provides information on the generation of the dynamo. The ubiquitous magnetic activity in close binary stars suggests the presence of a strong magnetic dynamo in the stellar interiors which, according to current theory, is generated by convection and differential rotation. Given their strong tidal effects, however, close binary stars are expected to show little or no differential rotation. It thus seems that differential rotation plays a minor role in the generation of the dynamo, or that the stellar interiors rotate differentially despite the observed tidal locking of the photospheres (Guinan & Giménez 1993). Eclipsing binaries further offer a unique advantage for the study of the spatial extent, the distribution, and the location of emitting material. Such systems may provide information on the location of the underlying dynamo. Some stars or binaries indicate the presence of solar-like coronal structures that are comparable with solar active regions (Schmitt & Kürster 1993), while observations of other systems suggest extended structures in particular between the components, perhaps related to intrabinary magnetic fields (Uchida & Sakurai 1983; Siarkowski et al. 1996). Giant magnetic structures have been identified on RS CVn binaries through eclipse mapping (White et al. 1990). Spectral differences at different phases of the eclipse may indicate radial structuring, but reported results are ambiguous and contradictory (Walter et al. 1983; White et al. 1990; Ottmann et al. 1993).
Strong polar magnetic fields have been proposed for magnetically active stars (Schüssler & Solanki 1992) based on arguments involving Coriolis forces in the stellar interior. Doppler imaging of active stars indeed suggests active regions concentrated at high latitudes (e.g., Strassmeier & Rice 1998). Corresponding coronal geometries have also been reported from radio observations (Benz et al. 1998; Mutel et al. 1998) and from X-ray observations of a flaring region that was occulted during an eclipse of Algol (Schmitt & Favata 1999). Radio Very Long Baseline Interferometry (VLBI) observations provide clues for very extended magnetospheres, including structures of the size of the intrabinary distance in binary systems (e.g., Mutel et al. 1985, 1998; Alef et al. 1997; Benz et al. 1998). The question then is whether such features are generally found on active stars, whether they contain hot plasma, whether they are preferred on single stars or on binaries, and whether there is additional strong activity at low latitudes. To address these questions, we have obtained a long observation of the Castor triplet of stellar binaries with the XMM-Newton observatory.
Castor is a hierarchical sextuplet system at a distance of
pc
(ESA 1997). The principal optical components
are Castor A (A1 V) and Castor B (A5Vm) in an orbit with a period of
467 yrs and a separation of presently
.
Each component
is accompanied by an invisible late-type star in a close orbit. Castor A's
companion is most likely to be a late K star (Güdel & Schmitt 1996) in
a 9.21 d eccentric orbit, while Castor B's companion appears to be an
early M dwarf (Güdel & Schmitt 1996) in a 2.93 d circular orbit (orbit
information from Batten et al. 1978; Heintz 1988; Haisch et al. 1990).
YY Gem (= Castor C) is a short-period spectroscopic binary and flare
star (e.g., Haisch et al. 1990) located 71
south of the
Castor system. It is gravitationally bound to Castor but only marginally so
(Anosova & Orlov 1991). It consists of two almost identical dMe stars in a circular 0.814 d
orbit. Due to the fortunate orbit inclination angle of
(Pettersen 1976),
the system produces nearly total primary and secondary eclipses, where the
primary has a slightly larger radius (0.66 vs. 0.58
)
and mass
(0.62 vs. 0.57
)
than the secondary (Haisch et al. 1990).
Doppler imaging shows magnetic photospheric spot
concentrations in mid-latitude bands (around
), with a weaker spot
at the equator (Hatzes 1995).
YY Gem has been extensively observed in X-rays. Pallavicini et al. (1990)
report quiescent and flaring emission not only from YY Gem itself, but they
detect the (unresolved) Castor AB system as a somewhat weaker X-ray source as
well. This led to some speculation about possible X-ray emission from A stars.
Later observations obtained by ROSAT (PSPC & HRI) and by ASCA
addressed this problem further (Schmitt et al. 1994; Gotthelf et al. 1994).
Castor AB was identified again, revealing a surprisingly high rate of
flares. Positional analysis indicated that the emission is most
probably coming from Castor A.
The observed X-ray luminosity requires a mass
if a single star is responsible (Schmitt et al. 1994).
Schmitt et al. (1994) present the first high-resolution radio
observations of the Castor system, revealing radio emission not only
from YY Gem as reported earlier (Linsky & Gary 1983), but also from
Castor A, where the flux ratio
is
similar to the X-ray flux ratio. The discovery of the third radio
source, coincident with Castor B, showed that
all three binary systems are "magnetically active'' (Güdel & Schmitt 1996).
VLBI studies revealed a coronal size of
cm
for YY Gem, corresponding to the intrabinary distance (Alef et al. 1997).
The absence of a detected radio eclipse supports the large size.
However, indications of deep eclipses were reported from X-ray
observations (Haisch et al. 1990), clearly suggesting different structures
for radio and X-ray emitting plasma.
Instrument | UT range | JD range 2451659.0 + |
RGS | 18:20:42 - 00:49:19 | 0.26438 - 0.53425 |
02:54:14 - 21:15:19 | 0.62100 - 1.38564 | |
MOS | 19:22:12 - 00:45:25 | 0.30708 - 0.53154 |
03:55:40 - 20:52:16 | 0.66366 - 1.36963 | |
pn | 18:42:50 - 00:47:50 | 0.27975 - 0.53322 |
03:16:26 - 20:54:41 | 0.63641 - 1.37131 |
The observations reported here were obtained by XMM-Newton
(Jansen et al. 2001) in two closely separated intervals (Table 1).
Two EPIC MOS cameras (Turner et al. 2001) provide an angular resolution
of 4
,
just comparable with the separation
between Castor A and B. Both MOS were operated in the "partial window mode 2'', thus
including both Castor and YY Gem.
The EPIC pn camera (Strüder et al. 2001) was operated in the small
window mode, with Castor AB outside the field of view.
The two RGS detectors (den Herder et al. 2001) were in spectroscopy mode.
The optical monitor (Mason et al. 2001) was in a closed position,
due to the strong optical light of Castor AB.
All data were reduced with the XMM-Newton Science Analysis System (SAS) software version 4.1, using several updates of individual tasks. For the RGS data, standard processing was performed using the RGSPROC task, followed by spectral extraction and response generation. The RGS fully separates Castor from YY Gem. The RGS spectral products were analysed in the Utrecht spectral software SPEX using the Collisional Ionization Equilibrium model CIE (Kaastra et al. 1996).
The EPIC MOS data were processed using EMPROC. For the spectral analysis, we
used a response matrix provided by Turner et al. (2001).
The EPIC pn data were reduced
with EPPROC. A preliminary calibration was applied
to the small window data. We used EVSELECT for source extraction.
Spectral analysis was performed in XSPEC (Arnaud 1996) using multi-T
models for optically thin plasma with variable abundances (VMEKAL)
that are compatible with the CIE model.
![]() |
Figure 1: EPIC MOS1 image of the Castor system. The components are labeled. Insets: contour plots of Castor AB data subsets. Contours are normalized relative to peak. Top: only flares on Castor A. Middle: low-level emission. Bottom: only flares on Castor B (see Fig. 2b) |
Open with DEXTER |
![]() |
Figure 2:
Top: Co-added light curve of YY Gem
from MOS1&2 and pn binned to 300 s
(the first ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3: RGS1&2 fluxed spectra of YY Gem and Castor AB (bin resolution: 0.04375 Å). Note the different vertical scales in the two panels |
Open with DEXTER |
The EPIC MOS show for the first time spatially separated images of Castor A and B and thus prove that both Castor binaries (i.e., most probably their low-mass companions) are X-ray sources (Fig. 1), compatible with their radio detections (Güdel & Schmitt 1996). Their separation and position angle coincide with the optically determined values (Heintz 1988). Overall, Castor A is X-ray brighter. Both stars flare frequently (Fig. 2b), in fact at a surprisingly high rate, so that no part of the light curve is constant, with Castor B's flare rate being somewhat higher. A closer look at Fig. 2b shows that flares on Castor B reveal rapid decays, while on Castor A they decay less rapidly. The light curve of YY Gem will be discussed further below.
Figure 3 shows the coadded and fluxed spectra of YY Gem and Castor AB
from all RGS observations. A series of lines of highly
ionized Fe and several lines of the Ly and He series are visible, most notably from
C, N, O, and Ne. The significant continuum
indicates the presence of hot plasma in both sources.
The flux ratio between the forbidden (f) and the intercombination (i) lines of the
He-like "O VII triplet at 22 Å is density sensitive in the range
109-1013 cm-3
(Mewe et al. 1985).
We have fitted (Fig. 4) the triplet for the quiescent RGS1 spectrum
of YY Gem in SPEX and found
cm-3,
where the errors signify 90% confidence limits.
To confirm this value, we derived the individual f and i fluxes
by fitting the line spread function to the lines, to obtain
a flux ratio
.
Given that the emission is most probably
originating from plasma with
MK (Sect. 4.3), the new calculations
by Mewe et al. (2001) for the "O VII triplet indicate
cm-3, in agreement with SPEX.
For Castor, we find a similar value from SPEX,
cm-3,
although it is not constrained
toward low densities.
![]() |
Figure 4: Density-sensitive He-like "O VII triplet of quiescent YY Gem |
Open with DEXTER |
---
![]() |
Castor AB | |||
Parameter | RGS2 | MOS1 | pn | MOS1 |
kT1 | 0.29+0.05-0.04 | 0.37+0.03-0.05 | 0.25+0.02-0.02 | 0.17+0.8-0.07 |
kT2 | 0.66+0.04-0.03 | 0.74+0.05-0.06 | 0.60+0.03-0.02 | 0.67+0.02-0.02 |
kT3 |
![]() |
1.82+0.42-0.30 | 1.12+0.21-0.10 |
![]() |
logEM1 b | 51.98+0.12-0.18 | 52.07+0.06-0.16 | 51.80+0.10-0.08 |
![]() |
logEM2 b | 52.26+0.11-0.13 | 52.21+0.07-0.08 | 51.98+0.13-0.09 | 51.75+0.13-0.14 |
logEM3 b | 51.67+0.36-0.70 | 51.65+0.17-0.13 | 51.85+0.10-0.11 | 51.02+0.19-0.41 |
C | 0.48+0.30-0.16 | = 1 | = 1 | = 1 |
N | 0.60+0.35-0.19 | = 1 | = 1 | = 1 |
O | 0.29+0.13-0.06 | 0.35+0.04-0.04 | 0.38+0.04-0.04 | 0.22+0.17-0.11 |
Ne | 0.75+0.38-0.15 | 0.63+0.09-0.10 | 0.74+0.34-0.29 | 0.46+0.19-0.14 |
Mg | 0.17+0.27-0.17 | 0.17+0.07-0.07 | 0.30+0.10-0.09 | 0.26+0.15-0.10 |
Si | 0.44+0.73-0.44 | 0.29+0.09-0.08 | 0.42+0.09-0.08 | 0.25+0.11-0.09 |
S | 0.14+0.23-0.14 | 0.27+0.17-0.16 | 0.36+0.15-0.15 | 0.47+0.32-0.28 |
Ca | = 1 | = 1 | = 1 | = 1 |
Fe | 0.14+0.05-0.03 | 0.15+0.03-0.02 | 0.24+0.05-0.04 | 0.12+0.06-0.04 |
Ni | 0.21+0.23-0.21 | 0.00+0.16-0.00 | = 0.21 | 0.13+0.35-0.13 |
![]() |
183/321 | 186/140 | 372/360 | 160/130 |
---
![]() |
||||
Parameter | MOS1 | MOS2 | pn | |
kT1 | - | - | 0.77+0.07-0.06 | |
kT2 | 3.22+0.32-0.32 | 3.47+0.39-0.35 | 3.21+0.27-0.22 | |
logEM1 b | - | - | 51.59+0.11-0.08 | |
logEM2 b | 52.58+0.03-0.04 | 52.54+0.04-0.05 | 52.49+0.03-0.03 | |
O | 0.15+1.04-0.15 | 1.80+1.54-1.24 | 1.53+0.61-0.52 | |
Ca | 3.50+3.24-3.26 | 0.00+2.24-0.00 | 4.83+2.69-2.71 | |
Fe | 0.74+0.25-0.22 | 0.85+0.36-0.29 | 0.74+0.16-0.15 | |
![]() |
84/124 | 123/123 | 369/354 |
a Errors give 90![]() ![]() |
were held fixed ("= x''). "...'' indicate unconstrained parameters. |
Abundances are relative to solar photospheric abundances (Anders & Grevesse 1989). |
b Due to the present lack of Good Time Intervals for pn data, we normalized the count |
rates with the length of the time interval; this implies a lower limit for the EM values for pn. |
c These flare components are additional to the fixed quiescent 3-T model. |
EMs are averages over rise and flare peak. |
We fitted separately RGS2, MOS and pn data for the quiescent YY Gem before the large flare at 1.35 d with 3-T coronal models (Table 2). For the MOS, the time interval was t=1.157-1.351 d (Fig. 2a). For the RGS we included data between t=0.85-1.19 d after which time a new RGS observing segment started. Both intervals show very similar hardness characteristics. We modeled the flare at t>1.35 d for MOS and pn by keeping the quiescent model found before fixed, and adding one or two further thermal models, for which we fitted abundances, T, and EMs. Multi-T fit results are reported in Table 2.
All three fits to the quiescent spectra agree quite well.
The most notable trends are low abundances for elements
with a low First Ionization Potential (FIP), most clearly for Fe,
while high-FIP elements, in particular Ne, show higher abundances.
The poor result for S from the RGS is due to weak lines with poor atomic
data (Audard et al. 2001). Figure 5 shows an EM distribution recovered
from the RGS2 spectrum, using Chebychev polynomials of degree 5. A broad distribution
from 2-15 MK compatible with the 3-T fit results is evident. There is little plasma with
T < 2 MK: a number of lines of "C VI, "N VI, "N VII, and
"O VII with formation temperatures in this range constrain this regime.
During the flare, some elemental abundances increase, in particular Fe which increases by a
factor of 3-5. Note, however, that the significance of the flare abundance values
for the other elements are clearly at the limit of the present calibration and
the S/N available.
![]() |
Figure 5: EM distribution for the quiescent YY Gem, derived from RGS2. A polynomial degree of 5 was used |
Open with DEXTER |
![]() |
Figure 6: YY Gem EPIC pn spectra (top: flaring; bottom:quiescent) |
Open with DEXTER |
Figure 6 illustrates spectral fits to the YY Gem pn data for quiescence
and for the large flare. More detailed time-resolved spectroscopy was performed with MOS,
with five segments selected as follows (Fig. 2a):
(1) secondary eclipse -
(2) flare peak at 0.8 d - (3) quiescent emission (around 1 d) -
(4) primary eclipse - (5) flare (including rise) at 1.37 d. Two T components
were sufficient for the given S/N ratio. The results are shown graphically
in Fig. 7. The heating effects during the flares are evident. Note also the
trend for higher metal abundances (in particular of Fe and Ne) during the large flare,
while no significant changes are seen during the eclipses.
How different is Castor AB from YY Gem? Close inspection of the RGS spectra in
Fig. 3 reveal higher flux ratios of
"O VIII/"O VII and "Ne X/"Ne IX for
YY Gem than for Castor, indicating that YY Gem's corona is on average hotter despite
the ongoing flaring on Castor. Some support for this result comes from multi-T fits
reported in Table 2 (from the integrated MOS1 observation).
The ratio between the EMs of the hottest to the intermediate component is larger for YY Gem.
The abundance pattern is, within the errors, rather similar.
![]() |
Figure 7:
Time resolved MOS1 (diamonds) and MOS2 (triangles) spectroscopy of YY Gem based
on 2-T spectral fits, showing the parameters kT1, kT2, abundances of O, Ne, Mg, Si,
S, Fe, and Ni. The numbers on the x-axis identify time
intervals (see text and Fig. 2a). Error bars indicate 90% confidence ranges based on
![]() |
Open with DEXTER |
![]() |
Figure 8:
Line broadening of the "O VIII Ly![]() |
Open with DEXTER |
Radial velocity variations due to the orbital motion of the two stars in YY Gem may be evident as broadening of the spectral lines during the orbital period. Eclipse times and orbital phases were taken from Haisch et al. (1990) who predict a primary eclipse at JD = 2451660.11676, consistent with the three eclipses seen in the XMM-Newton observations (Fig. 2a).
The RGS spectra for the second observation were subdivided into time
segments for which maximum or minimum line broadening would be expected (i.e.
around quadrature
or 90
and conjunction 0
and
). Two
segmentations were used, firstly the observation was divided into 8 segments
with segment 6 centered on the primary eclipse. Each segment, with the exception
of segment 1, was approximately 8800 s in duration. A second
set of time segments, shifted by about 2200 s relative to the first, was
also constructed which optimized the counts in the first segment.
Line broadening and centroids for the strongest emission line in the spectrum,
that of "O VIII at 18.969 Å, were determined firstly with SPEX
using the calibrated
line spread function in the response matrix for the RGS and secondly within the software
package DIPSO (Howarth et al. 1998) by correcting for the
effective area variation and fitting a combination of two Gaussian functions, which
was found to give an acceptable representation of the line spread function.
The variation of line broadening for RGS2 (average for the two segmentations)
for each of the methods used
is shown in Fig. 8. It shows the expected sinusoidal variation with orbital
phase, with minimum broadening occurring near conjunction (phases of 0
and
).
It should be noted that the value for the 8th segment was determined
without including the flare. The amplitude of the variation is approximately
12 mÅ, corresponding to a velocity of 190 km s-1 which is in reasonable
agreement with the expected value of up to 160 km s-1, determined from the orbital
motion (including the stellar sizes). No anomalous broadening, resulting from the
inhomogeneity seen in the light curve modeling (Sect. 5), was detected.
The sinusoidal variation is not clearly evident in the results obtained from the RGS1 spectra. Although there is some variation during the observation it is not clearly correlated with the orbital period. We note, however, that although the detailed description of the response of both spectrometers is expected to evolve over time, we have preferred RGS2 for the present analysis as this description is more advanced than for RGS1.
The long light curve available for YY Gem includes three eclipses (two primary
and one secondary eclipse). We have applied a method first presented by
Siarkowski (1992) and Siarkowski et al. (1996) for a 3-D deconvolution
of X-ray light curves. In short, the stellar environment is divided into
regular volume bins that are occupied by variable amounts of EM.
After each iteration step, the model and the observed light curves
are compared. For each phase, the EMs of all visible
volume elements are increased by a factor equal to the ratio of observed
flux to modeled flux.
We have slightly modified the original procedure by applying volume bins
defined in spherical coordinates.
Bins are located at positions
.
In the radial direction,
we subdivide the interval from R* to 1.99 R* for the primary
(2.15 R* for the secondary) into 16 regular bins. In
(co-latitude)
direction, we define 36 equidistant bins from 0 (north pole) to
(south pole).
Finally, the width of the volume bins in
(longitude) direction is
such that all volumes at a given radius are equal, i.e.,
,
with 72 equidistant bins at the equator. To prevent ambiguities at
the pole, we keep at least four
bins at each co-latitude.
Several words of caution are in order. i) The light curve clearly shows flares
which must not enter into the modeling process. We therefore removed obvious
flares at 0.48 d, 0.7 d, 0.8 d, and 1.37 d (Fig. 2a), and linearly interpolated from the
pre-flare to the post-flare count rate. The gap around 0.6 d was also
linearly interpolated. This is a serious limitation of the present modeling
also in the light of the possible presence of further small flares.
We selected the interval 0.55-1.36 d for modeling, corresponding to
one orbit period. ii) Light curve inversion is an ill-posed problem and generally
satisfies multiple solutions (Siarkowski 1992; Siarkowski et al. 1996). Physically reasonable spatial restrictions
should be introduced to confine the problem to sensible solutions. We have done
so by restricting the radial extent of the coronae to 1 R*above the photosphere but, as we will see, this constraint will automatically be
fulfilled for most of the EM. iii) There is much correlation between neighboring bins, depending
on the time resolution and the signal-to-noise ratio of the measurements.
The finest structures in the solutions are ambiguous or not
significant. iv) Since the observer is located close to the orbital
plane, some north-south ambiguity could be introduced,
i.e., features are mirrored at the equator.
Some features may also be found on both stars at similar coordinates relative
to the star centers.
![]() |
Figure 9:
Model for the coronal structure in the dMe binary YY Gem derived from light curve
inversion. The primary star is in the center and is eclipsed by the secondary at phase 0.
The panels show the aspect in phase steps of 45![]() |
Open with DEXTER |
![]() |
Figure 10: Light curve modeling results for YY Gem. Primary and secondary components are designated with "1'' and "2'' and shown by solid and dashed lines, respectively. a) Distribution of EM in stellar latitude. b) Average radial density profiles; vertical lines indicate photospheric levels. c) Distribution of density values in all bins around each star. d) Relative luminosity contributions from plasma of given density (assuming volume filling factor = 1 in each volume bin). The spectroscopically measured density ranges (two methods, see Sect. 4.2) are indicated at the bottom |
Open with DEXTER |
With these caveats in mind, our goal is not to reconstruct a real coronal
image but to investigate statistical properties of physically reasonable
solutions that are compatible
with the light curve. While the most stringent restrictions are due to
the eclipses, the most difficult part of the modeling is the underlying
quasi-steady emission. Any
axisymmetric distribution with matching total luminosity can reproduce this
emission if it is sufficiently extended or lies close to the uneclipsed north
poles of the stars.
We have studied three special but suggestive cases that are
specified as starting conditions: i) a diffuse underlying corona with
;
ii) a diffuse underlying corona with
const. iii) Addition of a strong,
uneclipsed active region at each of the two north poles.
Case (iii) suppresses diffuse axisymmetric components during the iteration, essentially leaving the inhomogeneous structure for modeling. The latter was found identically from all three starting conditions. Figure 9 illustrates a possible model for one orbital revolution of the two stars as seen from the Earth (solution for case [i]; the primary star is held in the center). Note that
![]() |
(1) |
Figure 10c shows the statistical distribution of densities in all bins considered.
While most bins are of low density (i.e., the high-density plasma has a
small filling factor), the emissivity of a plasma scales with
so that most of the observed emission stems from plasma with densities
of log
(Fig. 10d). This is again in gratifying agreement with
the explicit density measurements in the RGS spectra. Finally, the total luminosities
of the two stars are very similar: the non-axisymmetric portions modeled
above contribute
1.3 1029 erg s-1 and
1.5 1029 erg s-1 from
the primary and the secondary, respectively.
The total system luminosity is
3.7 1029 erg s-1 (as determined from
the MOS data).
We have addressed a number of issues relevant for stellar coronal structure by making
use of a high-quality XMM-Newton observation of the Castor/YY Gem system.
The eclipse observations allow us to assess the
statistical properties of the coronal plasma, such as its
density scale height and the location of strongly emitting features.
Our 3-D light curve inversion shows excellent agreement with spectroscopic
results: i) the average density scale height implies a temperature that
is in good agreement with the range measured from
the available spectroscopy; ii) densities that dominate the X-ray emissivity
agree well with spectroscopically determined average densities
of a few times 1010 cm-3. iii) The active regions are mainly distributed
between degrees in latitude, in agreement with previous findings from
Doppler imaging (Hatzes 1995). The preference of mid-latitude regions is
at variance with reports on polar magnetic spots on active stars. We
may only speculate here whether the tidal interaction in this close system
of almost fully convective M stars plays a role for this scenario. On the
other hand, eclipse modeling does not exclude the presence of active regions
near the north polar regions of the system.
At any rate, there is clear evidence
that the coronae are highly structured and that both components
are similarly active:
i) Deep eclipses require similar amounts of EM on both stars, and
relatively compact active regions. We find no indications for
very extended magnetospheres as inferred from radio observations (Alef et al. 1997).
This is, in retrospect, not surprising: radio emission is due to high-energy electrons
trapped in magnetic fields; they radiate more efficiently in a low-density
environment where they suffer fewer collisions. On the other hand, efficient X-ray emission
originates preferentially in high-density environments since the emissivity scales
with
.
Such places are naturally found closer to the stellar surfaces.
ii) The asymmetries seen in the eclipse light curves clearly require structured
coronae. And iii) we see indications for periodic line broadening, compatible
in phase and amplitude with two orbiting stars of similar brightness.
We have investigated the thermal and elemental abundance structure both during quiescence and during flares. We find a broad distribution of EM above 2 MK, with little material below. The EM distribution peaks at about 8 MK, with a tail that may reach 15-20 MK as derived from EPIC data. There is no appreciable change during the eclipses, but we do see significant increases of some abundances during the larger flare that reaches temperatures around 40 MK.
Finally, we report the definitive detection and identification of both Castor A and B as flaring X-ray sources. Although the Castor AB source was detected in X-rays before (Pallavicini et al. 1990; Schmitt et al. 1994; Gotthelf et al. 1994), the individual coronal components have hitherto been detected only as radio sources (Güdel & Schmitt 1996). The new detections support the view that both Castor A and B have low-mass companions that are magnetically very active, as judged from the high flare rate.
Acknowledgements
We thank Slava Zavlin (MPE) who provided the CTI model for the PN small window mode. Frank Haberl (MPE) processed the pn data and made events files available to us. M. A. 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. H. M. acknowledges financial support from PPARC. E. F. and R. P. acknowledge support from the Italian Space Agency. N. G. is supported by the European Union (Marie Curie Individual grant; HPMF-CT-1999-00228). F. A. C. has been supported for data analysis by NASA grant NAS5-97119. SRON is supported financially by NWO.