A&A 383, 197-201 (2002)
DOI: 10.1051/0004-6361:20011725
Research Note
J. Lyytinen1 - P. Johansson1 - L. Jetsu1 - E. Esko1 - T. Hackman1 - D. S. Hall2 - G. W. Henry3 - S. Kontinen1 - P. Könönen1,4 - S. Maisala1 - A. Palviainen1 - K. Ryynänen1
1 - Observatory, PO Box 14, 00014 University of Helsinki,
Finland
2 - Dyer Observatory, Vanderbilt University, Nashville, TN 37235, USA
3 - Center of Excellence in Information Systems,
Tennessee State University,
330 10th Ave. North, Nashville,
TN 37203-3401, USA
4 - Helsinki University of Technology, Metsähovi Radio
Research Station, Metsähovintie, 02540
Kylmälä, Finland
Received 25 August 2000 / Accepted 26 November 2001
Abstract
Our time series analysis of sixteen BV
light curves of the chromospherically active binary V511 Lyr confirmed
the
rotation period uniquely.
The seasonal periodicity changes of 3.8% indicated
the presence of detectable surface differential rotation.
The significant
periodicity in
the light curve minimum epochs suggested that
the strongest starspot formation in V511 Lyr was concentrated
on one stable active longitude.
Key words: stars: individual: V511 Lyr (HDE 337518) - stars: variable: general - stars: activity - techniques: photometric
Bidelman (1983) noticed
the chromospheric CaII K emission of V511 Lyr (HDE 337518).
Due to its strong coronal emission,
V511 Lyr appears in both catalogues of the bright extreme ultraviolet
sources discovered by the ROSAT-satellite
(Pounds et al. 1993; Pye et al. 1995).
It has been identified as the optical counterpart of
the ROSAT source 2RE 190624+274308 (Mason et al. 1995) and
the EUVE source J 1906+274 (Lampton et al. 1997).
The spectroscopy by Mulliss & Bopp (1994)
revealed that V511 Lyr is a double-lined binary with
chromospheric H
and CaII 8542 Å emission.
Jeffries et al. (1995) determined the orbital period
and concluded that this system probably consists of
two mid-K main sequence stars with nearly circular orbits.
The photometric rotation period,
,
obtained by Henry et al. (1995)
indicated synchronous rotation.
They estimated the spectral-types of the primary (K0V)
and the secondary (K3V), and also measured the respective rotation
velocities
and
.
Fekel (1997) determined the values of
and
.
The most recent estimates are
(G8V) and
(K3V) in Osten & Saar (1998).
We performed a detailed time series analysis of two years of V511 Lyr photometry. Our modelling results are compiled into Table 1, while the data are published only in electronic form in Table 2.
![]() |
Figure 1:
The differential V magnitudes V511 Lyr minus HDE 337653 with
the ephemerides
![]() |
Open with DEXTER |
The differential BV photometry of V511 Lyr was made between
May 1994 and June 1996 with the Vanderbilt-Tennessee State
0.4m Automated Photoelectric Telescope (APT).
The comparison and check stars were
HDE 337653 (C1) and HD 179422 (C2).
Our Table 2 contains seventeen subsets (1st column: SET),
the heliocentric Julian dates (2nd column: HJD)
and the differential BV magnitudes
V511 Lyr minus HDE 337653 (3rd and 4th columns).
Note that we analysed sixteen subsets,
because
contained only two nights of observations.
A suitable subset length was below one month,
because the light curves remained stable for about ten rotations
(Fig. 1:
).
The time distribution of the data naturally contained gaps
without observations, which had to be taken into account when
dividing these data into subsets.
The selected subsets contain a comparable amount of data,
e.g. the number of observing nights has a range
between 12 and 19 in all analysed subsets (Fig. 1).
The photometry was corrected for differential extinction and
transformed to the standard Johnson system.
Each observation was a mean of three measurements.
All observations with a standard deviation greater than
were automatically rejected.
These APT data have already been described by Henry et al. (1995),
who studied the
photometry.
Their external precision estimate,
,
relied on observations of "selected pairs of especially stable stars''.
This is comparable to our standard deviation of
between the seasonal
and
values.
But the standard deviations within the seasons varied
from
to
in
,
and
from
to
in
.
This implies weak variability in C1 or C2.
Considering the spectral-type (Harlan 1969: F5V)
and the absence of CaII H&K emission (Wilson 1966),
photometric variability in HD 179422 is improbable.
No variability has been observed in this star
(e.g. Breger 1969;
Fernie 1982)
and it has been proved to be a reliable photometric comparison
star in studies of the cepheid variable HR 7308
(e.g. Burki et al. 1986).
The mean differences in all data were
and
.
Combining these with V=6.353 and B-V=0.416 of
HD 179422 (Arellano Ferro 1984) gave
V=8.789 and B-V=0.585 for HDE 337653.
Except for the spectral-type (Nesterov et al. 1995: F8),
the available information about HDE 337653 was limited.
Assuming no interstellar extinction,
our apparent V magnitude and the Hipparcos satellite 67.6 pc
distance estimate (Perryman et al. 1997) indicated
a main sequence absolute magnitude (
).
Our B-V for HDE 337653 would then fit the G0V spectral-type.
In conclusion, if there is real photometric variability in C1 or C2,
it is more probable in the former.
But this
variability
is much weaker than that in V511 Lyr,
and does not affect our analysis.
The TSPA-method from Jetsu & Pelt (1999: hereafter Paper I)
was applied to search for periodicity in the normalized magnitudes
of each subset.
Except for the short
,
the B and V magnitudes of each subset
were normalized (Paper I: Eq. (17)).
Our nonlinear second order TSPA model was
The seasonal V light curves revealed
that the spot distribution of V511 Lyr sometimes evolved
significantly during one month (Fig. 1).
The mean (M) reflects the total level of magnetic activity.
The amplitude (A) tells how nonuniformly this activity
is distributed in longitude.
A nearly constant seasonal brightness of V511 Lyr should not be
interpreted as absence of magnetic activity,
because low A occurred at very different M levels
(e.g. Fig. 1:
and 12).
The long-term MB and MV changes were smooth
(Figs. 2a-c).
The effective surface temperature decreased when the star became fainter,
i.e. the
MB-MV colour index correlated
with the mean luminosity (Fig. 2d).
The light curve amplitudes AB and AV
changed more abruptly than the mean,
the most extreme case occurring between
and 12
(Figs. 2e and f).
Because the time span of our photometry is only about two years,
the long-term changes in the mean brightness or the amplitude of the
light curve did not give any definite information about the presence
of an activity cycle in V511 Lyr.
Except for the aforementioned correlation between the colour index
and the mean brightness changes,
there were no obvious correlations between the
different light curve parameters displayed in Fig. 2.
However, the time series analysis of the two years of
photometry of V511 Lyr did allow us to study two important
phenomena already detected in several other chromospherically active stars,
namely the differential rotation (P changes) and
the active longitudes (
changes).
Henry et al. (1995) discovered
the
periodicity in the
photometry,
but noted that the
periodicity
was "nearly as strong''.
The one sidereal day window period (
)
connects these two alternatives and also
introduces other spurious periodicities (Paper I: Eq. (16)).
The obvious combination for V511 Lyr is
,
,
,
,
and
.
Our analysis with the TSPA-method phase residual- and
-criteria
(Paper I: Sect. 5)
confirmed uniquely that
is the real periodicity.
![]() |
Figure 2:
a) The differential V magnitudes
V511 Lyr minus HDE 337653 (Table 2: n=256).
b)-g) The
MB, MV,
MB-MV,
AB, AV
and
P estimates with
![]() |
Open with DEXTER |
If the latitude of the strongest concentrated active region determines
the seasonal light curve period,
the P changes measure the stellar surface differential rotation.
The weighted mean of the reliable P was
with an error of
(Fig. 2g: closed squares).
If no significant period changes occur,
these reliable period estimates should fit a model
of a constant period P'.
The best alternative for such a constant period model is
,
because the definition of
ensures that no other
constant period model can yield a lower
value.
We had to reject the hypothesis,
"The rotation period of V511 Lyr is constant
'',
because the reliable Pgave
having a significance
This supported the presence of differential rotation.
The ratio
can be used to measure the statistical
upper limit for the observed period changes.
For example,
Jetsu et al. (2000) list the values of
,
4.6% and 3.3% for
V1794 Cyg (G5III-IV),
V815 Her (G2V)
and
V368 Cep (K1V), respectively.
They also studied the relation between Z and the differential
rotation coefficients k of 15 other stars in Hall & Busby (1990).
Our reliable P estimates from Table 1 gave
for V511 Lyr.
It is also very probable that the full range of
the P changes in V511 Lyr has not been observed,
because the time span of our photometry is only about two years.
Yet, of all 15 stars studied in Hall & Busby (1990),
only the period changes in BY Dra exceeded those observed here in V511 Lyr.
All reliable P estimates were above the three unreliable ones
connected to low amplitude light curves
(Fig. 2g: open squares).
But there was no real correlation between A and P,
because four other low A light curves had
larger and reliable P estimates (,
6, 12 and 14).
A detailed TSPA analysis example of similar uncertainties encountered
with low amplitude light curves was presented in Sect. 6.2 of Paper I.
Although no correlation between the A and P
was found in the light curves of V511 Lyr,
such a correlation could be present
in some chromospherically active stars.
For example, in the case of a strong solar-type differential rotation,
where the rotation periods in the equatorial regions are shorter
than in the polar regions,
the higher A values could be more frequent
for the lower P values.
![]() |
Figure 3:
The primary and secondary minimum phases with the ephemeris
![]() ![]() ![]() |
Open with DEXTER |
Normalized magnitudes | Original BV magnitudes | |||||||||
SET |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|||
1 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
2 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
3 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
4 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||
5 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
6 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||
7 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
8 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
9 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
10 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
11 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
12 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||
14 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
15 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
16 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
||||
17 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
The light curve primary minimum (
)
occurs when the
strongest active region crosses the central meridian of the
visible stellar disk.
A secondary minimum (
)
represents
another weaker active region.
We analysed all reliable
with the Kuiper-method described in
Jetsu & Pelt (2000, hereafter Paper II).
This method was used to test the hypothesis (H0):
"There is no periodicity in the reliable
epochs of V511 Lyr''.
To cover all P in Table 1,
the periodicities between
and
were tested.
The best periodicity,
,
reached the critical level
(Paper II: Eq. (A.3)).
This critical level
represents the probability that the Kuiper-method
test statistic reaches the observed,
or an even more extreme, value.
Hence the distribution for the phases of the reliable
with the best period
was extremely improbable.
The significance level to reject H0 was fixed to
,
as in Jetsu (1996),
where active longitudes were detected in four RS CVn stars.
Although H0 could not be rejected in V511 Lyr with this
,
we emphasize that the critical level of
was
extremely improbable under H0.
Furthermore, our Rayleigh-method analysis gave the same best periodicity
with a critical level of
(Paper II: Eq. (A.4)).
The seasonal primary and secondary minimum phases
with
are displayed in Fig. 3.
The reliable
concentrated
inside an interval of 0.3 in phase.
The unreliable
and
appeared randomly distributed.
These results suggested that the strongest
starspot formation on V511 Lyr was concentrated
on one active longitude that rotated with a constant period
of
during the whole time span of our photometry.
Our time series analysis of the photometry has arrived at two
apparently contradictory results.
The first result was that the reliable seasonal P determined with
the TSPA-method exhibited very significant changes.
In other words,
the hypothesis of a constant rotation period in V511 Lyr had to rejected.
This was interpreted as a sign of the presence of differential rotation.
The weighted mean of all these seasonal periodicities was
.
The second result was that the Kuiper-method analysis of
a series of time points,
the reliable
estimates,
indicated the presence of an active longitude rotating with a constant
period of
.
The apparent contradiction is that
.
Furthermore, differential rotation indicates variable rotation periods,
but the active longitude indicates a constant rotation period.
Nevertheless, similar results have already been obtained in several other
chromospherically stars, e.g. in FK Com, V 1794 Cyg, V 815 Her and even
in the Sun (Jetsu et al. 1993, 1997, 1999,
2000).
In all these cases, the interpretation has been the same.
There may exist some long-lived structure(s) below the visible solar
or stellar surface that rotate with a constant angular velocity.
The formation of the sun- or starspots
is more frequent close to these structures.
When the sun- or starspots are disconnected from these structures,
the observed rotation periods depend on their latitudinal location,
i.e. on the law of surface differential rotation.
We studied one signature of magnetic activity in V511 Lyr,
the rotational and long-term luminosity changes caused by starspots.
The 3.8% changes in the reliable seasonal periodicities (P) determined
with the TSPA-method indicated the presence of surface differential rotation.
These seasonal periodicity changes were so large
that the hypothesis of a constant period had to be rejected.
The TSPA-method
also gave the light curve primary and secondary minimum epochs
when the active regions crossed the central meridian of the visible
stellar disk.
Our Kuiper- and Rayleigh-method analysis of these epochs indicated
that the starspot formation in V511 Lyr was concentrated on one
active longitude rotating with a constant period of
.
These results suggest the coexistence of differential rotation
(variable seasonal P)
and
active longitudes
(constant
for all
).
One possible interpretation could be a stable and strong magnetic
field structure deeper in the convective zone rotating with a constant period.
The variable seasonal periodicities could then represent
the latitude of the short-lived starspots emerging
to the surface above this structure.
The mean brightness (M) changes related to the total starspot area
were smooth and
correlated with the effective surface temperature.
The abrupt changes in light curve amplitudes (A) implied that
the longitudinally nonuniform part of the starspot distribution
sometimes evolved rapidly, even within one month.
The photometry also provided some new constraints for the binary system V511 Lyr.
Without any interstellar extinction,
the Hipparcos distance 50.5 pc (Perryman et al. 1997) and
the maximum apparent brightness V=8.87 implied an absolute magnitude of
MV1,2=5.35.
The most recent estimates for V511 Lyr were
(G8V primary) and
(K3V secondary)
in Osten & Saar (1998).
These particular spectral-types had
MV1=5.6 with
and
MV2=6.8 with
in Gray (1992).
These predicted
MV1,2=5.29 and
a 3.0 flux ratio between the primary and the secondary in V.
This fits our
MV1,2=5.35 nicely.
Due to the 3.0 flux ratio, the rotation period
obtained from
the light curves most probably represented that of the G8V primary
in a binary with
.
Assuming that the rotation of the less massive K3V secondary
is also synchronized,
our
and the
from
Osten & Saar (1998) yielded the projected radii
and
.
Furthermore,
combining these to the radii in Gray (1992)
gave consistent inclinations of
and
.
Acknowledgements
The research at the Observatory of Helsinki University was performed during spring 2000 by the participants of the lectures "Time Series Analysis in Astronomy''. The work by T.H. was supported by a grant from Wihuri foundation, and by Helsinki University research funding for "Time Series Analysis of Stellar Magnetic Activity'' (No. 974/62/98). Astronomy with automated telescopes at Tennessee State University is supported through NASA grants NCC5-96 and NCC5-228, which funds TSU's Center for Automated Space Science, and NSF grant HRD-9706268, which funds TSU's Center for Systems Science Research. We have made use of the Simbad-database operated at CDS, Strasbourg, France. We wish thank the anonymous referee for valuable comments on our manuscript.