A&A 487, 709-716 (2008)
DOI: 10.1051/0004-6361:200809688
L.-Y. Zhang1,2 - S.-H. Gu1
1 - National Astronomical Observatories/Yunnan Observatory, Joint laboratory for Optical Astronomy,
Chinese Academy of Sciences, Kunming 650011, PR China
2 - Graduate
School of Chinese Academy of Sciences, Beijing 100039, PR China
Received 29 February 2008 / Accepted 10 June 2008
Abstract
Aims. We present the new high-resolution echelle spectra of SZ Psc, obtained in Nov. 2004 and Sep.-Dec. 2006, and study its chromospheric activity.
Methods. By means of the spectral subtraction technique, we analyze our spectroscopic observations including several optical chromospheric activity indicators (the
D3,
D1, D2, H
,
and
infrared triplet lines).
Results. All indicators show that the chromospheric activity of the system is associated with the cooler component. We find that the values of
EW8542/EW8498 are in the range 1-3, which indicates optically thick emission in plage-like regions. The 2006 data suggest the presence of active longitude phenomena. For the
8542 and 8662 and the H
lines, it seems that the excess emission is stronger near the two quadratures of system. This may be anti-correlated with the behavior of the
D1 line. The absorption features are detected in the subtracted H
lines, which could be explained by prominence-like extended material seen on the stellar disk or by mass transfer from the cooler component to the hotter one.
Key words: stars: binaries: eclipsing - stars: binaries: spectroscopic - stars: chromospheres - stars: activity - stars: late-type - stars: individual
![]() |
Figure 1:
The
![]() ![]() ![]() ![]() |
Open with DEXTER |
The first variability of the lightcurve of SZ Psc was found by Jensch (1934). Since then, the light-curve variations have been analyzed by many authors (Jakate et al. 1976; Eaton 1977; Catalano et al. 1978; Eaton et al. 1982; Tunca 1984; Doyle et al. 1994a; Antonopoulou et al. 1995; Lanza et al. 2001; Kang et al. 2003; Eaton & Henry 2007). For the distortion wave on the lightcurves of SZ Psc, this has been explained by the starspot model (Eaton & Hall 1979; Lanza et al. 2001; Kang et al. 2003; Eaton & Henry 2007).
The spectroscopic observations have been carried out by many groups (Jakate
et al. 1976; Bopp & Talcott 1978; Weiler 1978; Bopp 1981; Ramsey &
Nations 1981; Huenemoerder & Ramsey 1984; Fernández-Figueroa et al. 1986; Popper 1988; Fernández-Figueroa et al. 1994; Frasca &
Catalano 1994; Montes et al. 1995; Eaton & Henry 2007).
Jakate et al. (1976) and Popper (1988) derived the spectroscopic
orbital elements on the basis of the radial velocity curves, and
they also found strong chromospheric emission in the
H & K lines from the cooler component. A flare-like
phenomenon was detected in the
H & K lines by
Fernández-Figueroa et al. (1986). Recently, Eaton & Henry
(2007) have revised the orbital elements of the system, but so far
there has been no spectroscopic study of the
IRT lines for SZ Psc.
SZ Psc is very active, because of its extremely variable
H profiles. In 1978, Weiler (1978) found that SZ Psc only
exhibited significant variation in its H
emission lines
near orbital phase 0.8. Moreover, Bopp & Talcott (1978) confirmed
this result and found H
emission only in the phase range 0.6-0.8. However, Ramsey & Nations (1981) found H
emission is in the phase range 0-0.4 and interpreted the
H
behavior as a flare event that locally injects material
into a transient disk or shell structure around the K1 IV star.
Afterward, Bopp (1981) discussed an unusual H
profile,
particularly its large width and double-peaked appearance. Beyond
that, Huenemoerder & Ramsey (1984) found that SZ Psc was highly
active with the emission from a cooler component and that it had
undergone a large H
outburst during and after which the
profiles were suggestive of a circumstellar origin. The extremely
variable H
profiles suggest there are two active
longitudes. Recently, Eaton & Henry (2007) have found that
H
emission seems independent of the orbital phase on
average, but it can increase markedly within
a few orbital cycles.
Up to now, the spectral subtraction technique has been used
widely for chromospherically active binary systems in several lines:
H & K,
b,
D3,
D1, D2, H
,
H
,
IRT (Gunn & Doyle 1997; Gunn et al. 1997; Lázaro & Arévalo 1997; Montes et al. 1997;
Arévalo & Lázaro 1999; Frasca et al. 2000; Montes et al.
2000; Shan et al. 2006; Gálvez et al. 2007, etc.). The technique
was described in detail for obtaining the chromospheric contribution
by Barden (1984) and Montes et al. (1995). Also, it was discussed
and applied in the search for prominence-like events (Gunn
& Doyle 1997; Gunn et al. 1997).
In this paper, new spectroscopic observations of SZ Psc are analyzed using the spectral subtraction technique. Some new results are derived for the chromospherically active activity of SZ Psc.
The new spectroscopic observations of SZ Psc were carried out with
the 2.16 m telescope at the Xinglong station of the National
Astronomical Observatories, China, in five observing runs:
Nov. 20-27, 2004, Sep. 1-6, 2006, Oct. 28-29, 2006, Nov. 28-30, 2006,
and Dec. 8-10, 2006. The Coudé echelle spectrograph with a
spectral resolution about 37 000 and a
pixel
Tektronix CCD detector were used during our observations. The
reciprocal dispersions are 0.079
/pixel for the
D3,
D1, D2spectral region, 0.089
/pixel for the H
spectral
region, 0.113
/pixel for the
8498
spectral region, 0.115
/pixel for the
8542 spectral region, and 0.117
/pixel for the
8662 spectral region. Correspondingly, the
spectral resolution determined as the FWHM of the arc comparison
lines is 0.150, 0.160, 0.207, 0.218, and 0.221
,
respectively.
Table 1: The observational log of SZ Psc.
The spectra were reduced using the IRAF package. The standard procedures were used, which include image trimming, bias subtraction, flat-field division, background subtraction, cosmic ray removal, and spectrum extraction. The wavelength calibration was obtained by taking the spectra of a Th-Ar lamp. Finally, the spectra were normalized by a polynomial fit to the observed continuum. The normalized spectral profiles are displayed in Fig. 1. During our observations, the signal-to-noise ratio (S/N) varied in quality, and in some instances it suffered from poor seeing and intermittent clouds. For most of the observations, the S/N is more than 100. The observing log is listed in Table 1, which includes the observing date, the Heliocentric Julian Date of observation (HJD), orbital phase, exposure time, and the S/N. The orbital phases were calculated using the revised ephemeris derived by Eaton & Henry (2007).
On some nights of our observing seasons when the telluric
lines were heavy, two rapidly rotating early-type stars, HR 7894
(B5IV,
km s-1) and HR 989 (B5V,
km s-1),
were observed. Their observed spectra were fitted using the
high-order spline3 function to derive the telluric line templates.
For the spectra of HR 989, interstellar
lines
were present and removed. The telluric lines in the spectra of SZ Psc were removed by means of these templates obtained on the same
nights. The examples of this procedure in different spectral regions were given by Gu et al. (2002).
![]() |
Figure 2:
Samples of the observed, synthesized, and subtracted
spectra for the
![]() ![]() ![]() |
Open with DEXTER |
The normalized spectra of SZ Psc are analyzed using the spectral subtraction technique. In this method, the synthesized spectra is constructed from artificially rotationally broadened, radial-velocity shifted, and weighted spectra of two inactive stars with the same spectral type and luminosity class as the two components of the active system. For our situation, the synthesized spectrum is constructed by means of the program STARMOD (Barden 1985). We observed several stars with spectral types and luminosity classes similar to those of the individual components of SZ Psc. By comparison, we find that two inactive stars HR 7948 (K1IV) and HR 6669 (F8V) are much better templates for SZ Psc. Thus, we chose them as template stars to construct synthesized spectra for all spectral data of SZ Psc. For the hotter component, Strassmeier et al. (1993) gave F8 IV for the spectral type. However, we find that HR 6669 is a better template for the hotter component. Moreover, HR 6669 does have the advantage of being an extremely inactive star (Wright et al. 2004), and Gray et al. (2001) gave F8 V for its spectral type, so the spectral type of the hotter component should be F8V.
In the course of the analysis, the v sin i values (0/78 km s-1) of the components of SZ Psc were used, which are
determined from the spectra spanning the wavelength ranges 6389-6477
and 6615-6706
with many photospheric lines. These
values are close to 80 km s-1 for the cooler component and less
than 5 km s-1 for the hotter component, as published by Eaton
& Henry (2007). The intensity weight for each of the two components
is obtained from high-quality spectra around phases 0.25 and 0.75 where the two components are separated well. Consequently, the
adopted intensity weight ratios are 0.25/0.75 for
D3 and
D1, D2 spectral
region, 0.20/0.80 for H
spectral region, 0.18/0.82 for
the
8498 spectral region, 0.175/0.825 for the
8542 spectral region, and 0.17/0.83 for the
8662 spectral region.
Samples of the synthesized and subtracted spectra (the
observed spectra minus the synthesized one) obtained in the
D3,
D1, D2,
and
IRT line regions are displayed in Fig. 2.
Since the H
features are bizarre, some samples are
displayed in Fig. 3. For each sample, we plot the observed (lower
solid line), the synthesized (dashed line), and the subtracted
spectra (upper solid line), which is offset for better display. In
these figures, ``H'' and ``C'' represent the hotter and cooler components, respectively.
![]() |
Figure 3:
Samples of the profiles of the H![]() |
Open with DEXTER |
Table 2:
The measurements for the excess emissions
of the
D1, D2, H
,
and
IRT lines.
The new spectroscopic observations allow us to simultaneously analyze the behavior of the different chromospheric activity indicators formed at different atmospheric heights.
The equivalent widths (EWs) of the excess emissions were measured using the IRAF SPLOT task. When the profiles of two components are easily separated, the EWs were evaluated on the subtracted spectra by integrating over the emission profile. In addition, we measured the EWs using the Gaussian fit (For asymmetric profiles, we used several Gaussian profiles to fit), If the profiles are blended, we measured the EWs using several Gaussian fits. For comparison, we also measured the EWs with the Lorenzian fit. The errors of the EWs were estimated by using the difference between the results of these different methods. The net EWs of the excess emissions are listed in Table 2.
Because of nearly a 4 integral-day period for SZ Psc, it was very difficult to observe effectively, and the phase coverage interval of observation was short per night at a single site, so we only had a few spectra per night. To discuss the rotational modulation of chromospheric activity, we used all of the 2006 data because we only have 3 to 4 data points per rotation in different epochs.
The
D1, D2 lines are formed in the
upper photosphere and lower chromosphere. The filled-in absorption
of the core of the
D1, D2 lines
could
be used as chromospheric activity indicators (Andretta et al. 1997; Monte et al. 1997).
The
D lines are characterized by deep
absorption (see Fig. 1). The core of the lines does not show any
reversal. The application of the spectral subtraction technique
reveals that the cores of the
D1, D2lines are weak filled-in absorption for the cooler component. In
some cases, the
D1, D2 lines exhibit
obvious excess
emission from the cooler component.
The measured EWs of the excess emissions in the
D1, line are plotted vs. the orbital phase
in Fig. 4. For
D1, it seems that there are
two minima around the quadratures. That suggests an active longitude
phenomenon during our
observations.
The
D3 line is formed in the upper
chromosphere and is an activity indicator for the Sun and late-type
stars. The emission of the line is a probe for detecting flare-like
events (Zirin 1988). In our set of 20 spectra
in two different observing seasons, we observed no flare-like episodes.
The
IRT lines are important chromospheric
activity indicators for the Sun and late-type stars (Shine & Linsky
1972; Shine & Linsky 1974; Linsky et al. 1979; Foing et al. 1989;
Dempsey et al. 1993; Gunn & Doyle 1997; Montes et al. 1997, 2000; Andretta et al. 2005). They are formed in the lower
chromosphere making them sensitive probes of the temperature minimum
region (Montes et al. 1997).
![]() |
Figure 4:
The EWs of the excess emissions vs. orbital phase for
![]() ![]() |
Open with DEXTER |
For all the observed spectra (Fig. 1), we can see that the
IRT absorption lines are from both components
and that the cooler component shows strong filled-in absorption or
sometimes even small self-reversal core emission. For the spectra
during Nov. 28-30, 2006 and Dec. 8-10, 2006, we note that the
8498 line exhibits weak self-reversal emission
in the line core, while the
8542 and 8662 lines show little or no such
self-reversal core emission. For the spectra during Oct. 28, 29, 2006,
especially the spectra on Oct. 28, 2006, all
IRT lines exhibit obvious self-reversal in the line core. This is
the strongest for the
8498 line, a little
weaker for the
8542 line, and distinctly the
weakest for the
8662 line.
The spectral subtraction reveals that the cooler component
shows clear
IRT excess emission. For the
subtracted spectra, we measured the EWs of the excess emission and
plotted vs. orbital phase in Fig. 4. From this figure, we find no
significant trend in the
8498 data. For the
8542 and 8662 lines, especially for the
8542 line, it seems that there
are two active longitudes near phases 0.25 and 0.75. In addition, the excess emission is also stronger (within the errors) at phase 0.09.
The H line is a very useful indicator of chromospheric
activity, and it is formed in the middle chromosphere (Montes et al. 1997, etc.).
The H line exhibits variation from absorption to
weak emission (see Fig. 1). Also, it is usually filled in with broad
emission. All the subtracted H
spectra indicate that the
cooler component is active. Some excess emission profiles are
broader and centered on the cooler component, for example, the
spectra at phases 0.941 and 0.780. Furthermore, there are obvious
absorption features in the subtracted spectra, at orbital phases
0.171, 0.424, 0.345, 0.663, 0.677, and 0.687. The EWs of the excess
emissions for the 2006 data are plotted against orbital phase in
Fig. 5. The orbital modulation of the H
excess emission is
similar to
8542 and 8662 data. It seems that
there are also two active longitudes near phases 0.35 and 0.75.
![]() |
Figure 5:
As Fig. 4, but for the EWs of H![]() |
Open with DEXTER |
In this paper, we analyze simultaneous spectroscopic observations of several optical chromospheric activity indicators for the RS CVn eclipsing binary system SZ Psc using the spectral subtraction technique. We investigate the details of the excess emission and study the chromospheric emission variation with orbital phase.
For the
IRT lines, we find the cooler
component of SZ Psc shows filled-in absorption or sometimes even
weak self-reversal core emission. This might be due to the low-level
plage and high rotational velocity of the cooler component of SZ
Psc. For some other chromospherically active binary systems, Montes
et al. (2000) also detected a similar feature in the
IRT lines.
For the most obvious self-reversal core emission feature, it
is normal that the reversal is the strongest for the
8498 line, a little weaker for the
8542
line, and weakest of all for the
8662
line. This is consistent with the observation of the
IRT lines of solar plages with different degrees of activity
(Shine & Linsky 1972; Shine 1973; Shine & Linsky 1974). The
increase in self-reversal core emission in the
IRT lines results from the higher temperatures and higher electron
densities at a given line optical depth produced by the steeper
temperature
gradient (Shine & Linsky 1974).
For the ratio between EW8542 and EW8498, we find the values are generally low, between 1-3, which indicates optically thick emission in plage-like regions. In solar plage, the values of EW8542/EW8498 are in the range 1.5-3 (Chester 1991). These low values are also found in other chromospherically active binaries by other authors: Lázaro & Arévalo (1997); Arévalo & Lázaro (1999); Montes et al. (2000); and Gu et al. (2002).
All the H spectra indicate that the cooler
component is active. In the subtracted spectra, the excess emission
profiles exhibit broad wings. The broad emission feature coincides
with the results derived by other authors: Huenemoerder & Ramsey
(1984), Eaton & Henry (2007), etc. There are a
couple of possible explanations for this. First, the broad component
could be interpreted as arising from microflaring. Such broad
H
components are also detected in some other
chromospherically active stars (Hatzes 1995; Montes et al. 2000;
etc). Second, it might be caused by instantaneous mass transfer from
the cooler component to the hotter one (Bopp 1981). Because the
cooler component is filling 82
of its Roche lobe (Lanza et al.
2001),
it is very possible that the mass transfer can happen due to the strong chromospheric activity.
In summary, all the analyzed activity indicators show that
the cooler component is active. This is consistent with the result
that the
h&k emission lines are mainly
enhanced by the cooler component of the system (Doyle et al. 1994a,b; Kang et al. 2003). Moreover, for the
H and K
lines of SZ Psc derived by Fernáandez-Figueroa et al. (1986) and
Fernáandez-Figueroa et al. (1994), there is obvious
emission from the cooler component and no emission from the hotter one.
For the non-eclipsed configuration, the modulation of excess emission with phase indicates the appearance and disappearance of discrete active regions as the star is presenting different hemispheres to us.
For the cooler component, our result suggests the presence
of an active longitude phenomenon during the 2006 epoch. To
demonstrate the possible rotational modulation well, we used
polynomial function to fit the EWs of
D1,
IRT, and H
excess emission profiles,
which are displayed in Figs. 4 and 5. For the
8542 and 8662 and H
lines, the excess emissions (with the
orbital phase) may be correlated basically, especially around two
quadratures. It seems that the emissions are stronger around two
quadratures of the system (phases 0.25 and 0.75). However, for the
D1 line, it seems that the emission is
weaker around the two quadratures, so the
D1 line may be anti-correlated with the
8542 and 8662 and the H
lines. The reason might be that
the
lines could be affected by the presence of
active spot regions, because the reactions of neutral atoms in their
ground states are very sensitive to the photospheric temperature
structure (Barradoy Navascués et al. 2001). Moreover, there is
some evidence that two variable active regions appear around
longitudes
and
by analyzing a sequence of
V-band light curves of SZ Psc (Lanza et al. 2001). In addition, both
active longitudes were also suggested by
Huenemoerder & Ramsey (1984).
For another chromospherically active binary UX Ari, the
favorite active longitudes are also around two quadratures of the
system based on the analysis of several chromospheric activity
indicators
D3,
D1, D2, H
,
and
IRT lines
(Gu 2005). These active longitude features have also been found
using photometry on many other active RS CVn binary systems, such
as EI Eri,
Gem, HK Lac, and RT And
(Olah Hall & Henry 1991; Henry et al. 1996; Berdyugina & Tuominen 1998; and Zhang & Gu 2007).
Table 3: Observed properties of the absorption features identified by Gaussian fitting.
The interesting result from our
high-resolution H subtracted spectra is the detection of
absorption components at some orbital phases in our observations.
The physical properties of the absorptions are listed in Table 3,
which includes orbital phase, the radial velocity measured with
respect to the cooler component's frame of rest, and the EW of the
absorption. Similar
peculiarities can be found in the differential spectra of Huenemoerder & Ramsey (1984).
The absorption features may be explained by the stellar prominence. Stellar prominence results in relative absorption while passing between the stellar disk and the observer, and then scatters photons arising in the chromosphere out of the line of sight. This phenomenon was first reported by Collier Cameron & Robinson (1989) for the rapidly rotating star AB Dor. Since then, stellar prominence has been observed as transient absorption (usually) or emission crossing the profiles of the Balmer lines, not only in single stars but also in binary systems (Hall & Ramsey 1994; Cameron et al. 2002).
If the numbered features in Fig. 3 arise from prominences
(and are not simply due to less active stellar surface areas or the
velocity separation of the components.), we calculate the distances
of the prominences from the rotation axis of cooler component, using
the simple model (Collier Cameron & Robinson 1989). During
Oct. 28, 29, 2006, the distances of the prominences1, 2, and 3 (see
Fig. 3) are 3.8 R, 2.5 R, and
(R is the radius of the
cooler component), respectively. For the prominence3, if we let it
go back to the phase 0.432, the velocity of prominence3 should be
km s-1. On the other hand, the velocity of the
prominence1 detected at phase 0.432 is just inside the range. Thus,
prominence1 and prominence3 might came from the same place. During
the 2004 observing season, the velocities of the two absorptions are -111 km s-1 blueward and 97 km s-1 redward of the central
velocity of the cooler component at phase 0.171, and the velocity of
the absorption3 feature is about 57 km s-1 blueward of the
cooler component at phase 0.424. If absorption1 and 3 have the same
prominence, the distance of the prominence is 1.7 R. The EW of the
absorption3 is larger than one of the absorption1. This indicates
that the prominence absorption should be much more visible against
the disk of the cooler component from the phase 0.171 to 0.424, due
to the much greater projected area. For absorption2 at phase 0.171,
if the distance of the promience2 is 2.7 R (the average value of the
distances of our observed prominences), when moving to 0.424, the
radial velocity should be about 185 km s-1 redward of the
cooler component. It is just the critical point of the excess
emission, and absorption2 should not be visible against the
stellar disk at phase 0.424. This is consistent with our result that there is no absorption feature at the critical point of the redward excess emission at phase 0.424.
In the course of the spectral analysis, it was found that there are some absorption components around the primary profiles at some phases such as 0.663, 0.667, and 0.687. These absorptions may also be a mass transfer effect since the cooler component is very close to its Roche lobe and the mass transfer occurs more easily when considering the strong magnetic activity. These phenomena are also found in other RS CVn stars by other authors: Huenemoerder et al. (1989), Hall & Ramsey (1992), Gu et al. (2002).
Unfortunately, our data do not cover the H
line.
We cannot tell whether this is optically thick or thin structure via
the ratio
(Huenemoerder et al.
1989). It is impossible to define the nature of the absorption precisely at present.
The results from the analysis of our spectral data of SZ Psc can be summarized as follows.
Acknowledgements
The authors would like to thank the observing assistants of the 2.16 m telescope of Xinglong station for their help and support during our observations. We are very grateful to Dr. Montes for providing a copy of the STARMOD program. We would also like to thank the anonymous referee for his or her valuable suggestions and comments, which have led to significant improvement in our manuscript. We would also thank Prof. Chris Sterken and the team at the SWYA conference in Blankenberge, Belgium, which led to many improvements, including in the language of the revised manuscript. This work is supported by the NSFC under grants No. 10373023 and 10773027.