Issue 
A&A
Volume 600, April 2017



Article Number  A120  
Number of page(s)  5  
Section  Stellar atmospheres  
DOI  https://doi.org/10.1051/00046361/201629429  
Published online  12 April 2017 
Further evidence for a subyear magnetic chromospheric activity cycle and activity phase jumps in the planet host τ Boötis
Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany
email: jschmitt@hs.unihamburg.de
Received: 29 July 2016
Accepted: 4 October 2016
We examine the Sindex data, obtained in the context of the Mount Wilson H&K project for the nearby Ftype star τ Boo, for the presence of possible cyclic variations on timescales below one year and “phase jump” episodes in the observed Sindex activity levels, to determine whether such features are persistent properties of the chromospheric activity of τ Boo and possibly other latetype stars. Within the Mount Wilson H&K project τ Boo was observed during 1278 individual nights, albeit with a very inhomogeneous coverage ranging from 2 to 137 observations per year. Our analysis shows that periodical variations with timescales on the order of 110–120 days are a persistent feature of the Mount Wilson data set. Furthermore we provide further examples of “phase jump” episodes, when the observed Sindex activity drops from maximum to minimum levels on timescales of one to two weeks, hence such features also appear to occur on a more or less regular basis in τ Boo.
Key words: stars: activity / stars: chromospheres / stars: latetype
© ESO, 2017
1. Introduction
Recently the data obtained in the context of the socalled Mount Wilson H&K project have been publicly released and become available to the astronomical community at large. The hitherto most complete description of this Mount Wilson monitoring data set has been presented by Baliunas et al. (1995), who discuss the Sindex “light” curves of 112 stars; one of these is τ Boo (=HD 120136), which shows, according to Baliunas et al. (1995), a chromospheric cycle of 11.6 ± 0.5 yr. τ Boo is a bright (m_{V} = 4.5 mag) star of spectral type F7 with an M2 type companion, which was observed several thousand times during the Mount Wilson program. In addition to the 11.6 yr cycle described by Baliunas et al. (1995) there is evidence for cyclic variability on substantially shorter timescales.
τ Boo is now one of the brightest known planet hosts. τ Boo b, a socalled hot Jupiter, was discovered by Butler et al. (1997) with an orbit period of 3.3 days, which is very close to the rotation period of τ Boo A. For these reasons and its substantial magnetic activity, the τ Boo system is very interesting also from the point of view of starplanetinteraction (SPI). Lanza (2013) describes a theoretical framework linking a magnetic interaction between star and planet and shows that such effects could significantly enhance evaporative effects due to EUV and Xray radiation from the host. As far as τ Boo is concerned, Shkolnik et al. (2008) present Ca II K data and argue that these data may indicate SPI, however, it is probably fair to say that actual observational clearcut demonstrations of SPI are still very elusive.
Mittag et al. (2017) describe an Sindex variability of τ Boo with a timescale of about 120 days observed in the years 2013–2016 with the TIGRE facility (Schmitt et al. 2014) and show that the same period is also consistent with cyclic variability at Xray wavelengths. Similarly, Mengel et al. (2016), using their Sindex time series taken with the NARVAL instrument in the years 2007–2015, deduce a 117day period; we identify this period with the 116day period previously mentioned, but not described by Baliunas et al. (1997). Furthermore, based on spectropolarimetric observations with ESPaDOnS and NARVAL, subsequent Zeeman Doppler Imaging by Donati et al. (2008) provides evidence for magnetic field reversals on τ Boo; further observations of polarity changes with the same instrumental setup have been reported by Fares et al. (2009), Fares et al. (2013), and by Mengel et al. (2016), and the observed field reversals suggest a possible magnetic cycle of 240 or 740 days.
Using their TIGRE monitoring data, Mittag et al. (2017) show that the maxima of the Sindex data might well be associated with polarity reversals of τ Boo as is observed for the Sun, so that the “magnetic” cycle period would be twice as large. Finally, Mittag et al. (2017) report a “phase jump” of the activity cycle of τ Boo, i.e., an episode in the spring of 2016, when the activity more or less suddenly dropped from maximum to minimum values. The purpose of this research note is to investigate to what extent similar phenomena are present in the Mount Wilson time series of τ Boo.
Fig. 1 Mount Wilson Sindex time series of τ Boo. 
2. Observations and data analysis
2.1. Mount Wilson data
As already mentioned, the data obtained in the context of the Mount Wilson Project have recently been publicly released and can be downloaded^{1}. A detailed description of the data is also provided at this web site. The available data specifically include the star identification; the calibrated S index, which we use in this paper as the basis for our analysis; an instrument code indicating with which instrument the data was taken; and the date of the observation as well as other material.
Naturally, it suggests itself to investigate whether the star τ Boo also shows similar phenomena as observed in recent years and reported by Mittag et al. (2017) in the Mount Wilson time series. In the context of the Mount Wilson Program, τ Boo has been intensively observed for more than 30 yr and the data base lists more than 4000 entries of individual observations. However, a closer inspection of these data shows that, normally, τ Boo was observed three times per night, in a few cases even more often, and in one single night more than 100 observations had been obtained. On the other hand, especially in the early years, there is typically only one observation per night available. All observations available within a single night were averaged to obtain a more homogeneous data set and so we end up in this fashion with observations of τ Boo in 1278 individual nights spread over the years 1967–2001. In Table 1 we provide a list with the number of nights per year when τ Boo was observed. In every year between 1967 and 2001 H&K observations of τ Boo have been carried out, the number of nights varies from 2 (in 1980) to 137 (in 1984). The “light curve” produced from this data is shown in Fig. 1. One recognizes a few outliers with Sindex values above 0.25, which are ignored in our subsequent analysis.
Overview of available τ Boo Ca II H&K observations obtained in the Mount Wilson (MW) and TIGRE programs.
2.2. Discrete autocorrelation function
One of the methods to determine recurrent structures and periods in light curves is autocorrelation analysis. Normally one assumes data with equidistant temporal spacing, in which case the calculation of the autocorrelation for a given time lag τ is straightforward and described in every text book. In data collected in the context of groundbased astronomy (or geophysics for that matter) the recorded time series are usually irregularly sampled and hence the computation of, for example, an autocorrelation function becomes a bit more cumbersome.
While it is possible to resample a given light curve onto a regularly spaced time grid, in our opinion any application of such methods is precluded by the way astronomical data are usually sampled (for example, seasonal variations). In the context of AGN work, Edelson & Krolik (1988) introduced the socalled discrete correlation function as a method to compute cross or autocorrelations without the need to interpolate any data. Briefly, their method works as follows: we let s(t) be some function that is sampled at some (nonequidistant) times t_{i}, i = 1,N with values s_{i} = s(t_{i}). We can then consider all pairs of data points, which are characterized by some pairwise lag Δt_{i,j} = t_{i} – t_{j}; since we are dealing with autocorrelation functions (rather than crosscorrelation functions) we need to consider only positive lags (i.e., i>j). We let $\overline{)\mathit{s}}$ and σ_{s} denote the first two moments of the distribution function of the function s(t), the unbinned discrete autocorrelation function UBAF is defined as $\begin{array}{ccc}{\mathrm{UBAF}}_{\mathit{i,j}}\mathrm{=}\frac{\mathrm{(}{\mathit{s}}_{\mathit{i}}\mathrm{}\overline{)\mathit{s}}\mathrm{)}\mathrm{\ast}\mathrm{(}{\mathit{s}}_{\mathit{j}}\mathrm{}\overline{)\mathit{s}}\mathrm{)}}{{\mathit{\sigma}}_{\mathit{s}}^{\mathrm{2}}}\mathrm{\xb7}& & \end{array}$(1)Edelson & Krolik (1988) point out that in the case of data with significant errors the denominator in Eq. (1) has to be modified to ${\mathit{\sigma}}_{\mathrm{2}}^{\mathrm{2}}\mathrm{}{\mathit{e}}_{\mathit{s}}^{\mathrm{2}}$, where e_{s} is a typical error; since in our case the intrinsic scatter of the data is much larger than the measurement errors, we use Eq. (1) as is. Each time lag pair contributes toward the estimate of the discrete autocorrelation function (DAF), which can be constructed (in binned form) as follows: to determine DAF at some lag τ, we identify all time lag pairs, satisfying τ−Δτ/ 2 < Δt_{i,j}<τ + Δτ/ 2, where Δτ is the spacing of the time lags where the DAF is calculated. We let there be M such pairs, then $\begin{array}{ccc}\mathrm{DAF}\mathrm{\left(}\mathit{\tau}\mathrm{\right)}\mathrm{=}\frac{\mathrm{\u27e8}{\mathrm{UBAF}}_{\mathit{i,j}}\mathrm{\u27e9}}{\mathit{M}}\mathrm{\xb7}& & \end{array}$(2)Edelson & Krolik (1988) proceed with providing some error estimates for the soconstructed function DAF(τ). We compute errors by bootstrapping the data using permuted data sets with the same temporal sampling as the original data. Performing these calculations on the Mount Wilson data of τ Boo (shown in Fig. 1), we obtain the DAF shown in Fig. 2 (solid histogram); the dashed histograms show the extent of the 10 000 bootstrapped realizations comprising 67.5%, 90%, 99%, and 99.9% of the autocorrelation values at each time lag considered. As is clear from Fig. 2, there are highly significant autocorrelations in the DAF of the Mount Wilson τ Boo data. We recognize a peak at small time lags, while periods of ≈50 days have no measurable autocorrelation. For periods of ≈120 days and ≈240 days, however, we find a nonzero autocorrelation at a significance of more than 99.99%. The autocorrelation error distribution is asymmetric and attains maximal values at about 190 days or half a year. This can be attributed to the observing pattern of τ Boo, which can be observed well from about January to August.
Fig. 2 Discrete autocorrelation function of Mount Wilson Sindex time series (solid histogram). The dashed histograms denote the maximum values of the autocorrelation functions at the 0.67%, 90%, 99%, and 99.9% level; see text for details. 
As a sanity check, we perform the same computations as carried out for τ Boo for the daily solar Sun spot record, available online^{2}. That Sun spot record dates back to 1818, prior to that date the available daily observations are too sparse. We further note that the 198 yr since then comprise about 18 solar cycles, while the approximately 36 yr of Mount Wilson monitoring of τ Boo correspond to about 144 cycles. In order to simulate the seasonal Mount Wilson sampling we sample the solar data with the same sampling pattern as encountered for τ Boo and obtain the DACF for the Sun spot record shown in Fig. 3, which clearly shows the solar 11yr cycle and its harmonics. So clearly, the DACF retrieves the known solar cycle from the daily sun spot data when observed with a typical astronomical sampling pattern.
2.3. Lomb Scargle analysis
In order to characterize any periodicities in the activity of τ Boo further and to identify possible changes in the periodicity pattern, we investigate the changes on shorter timescales as follows. Since τ Boo disappears behind the Sun for a few months, it is natural to group the available data into observation seasons. Following our TIGRE approach, we then consider three adjacent seasons and up in this fashion with a total of 11 periods each comprising three observation seasons. We list the dates of these observational periods and the number of nights with observations of τ Boo in Table 2. Since we are interested only in variations on short timescales, we rectified the data for each season by subtracting the respective mean for that season.
LS period analysis in Mount Wilson data.
Zechmeister & Kürster (2009) derive a generalized LombScargle (LS) periodogram and show that the normalized power spectrum p(ω) at some frequency ω is given by the expression $\begin{array}{ccc}\mathit{p}\mathrm{\left(}\mathit{\omega}\mathrm{\right)}\mathrm{=}\frac{{\mathit{\chi}}_{\mathrm{0}}^{\mathrm{2}}\mathrm{}{\mathit{\chi}}^{\mathrm{2}}\mathrm{\left(}\mathit{\omega}\mathrm{\right)}}{{\mathit{\chi}}_{\mathrm{0}}^{\mathrm{2}}}\mathit{,}& & \end{array}$(3)where ${\mathit{\chi}}_{\mathrm{0}}^{\mathrm{2}}$ and χ^{2}(ω) are the values of the χ^{2} statistics at frequency values of zero and ω, respectively. Given some observed data o_{i} (i = 1...N) at times t_{i} with errors σ_{i} the χ^{2} statistics is calculated as usual through $\begin{array}{ccc}{\mathit{\chi}}^{\mathrm{2}}\mathrm{=}\sum _{\mathit{i}\mathrm{=}\mathrm{1}}^{\mathit{N}}{\left(\frac{{\mathit{m}}_{\mathit{i}}\mathrm{}{\mathit{o}}_{\mathit{i}}}{{\mathit{\sigma}}_{\mathit{i}}}\right)}^{\mathrm{2}}\mathit{,}& & \end{array}$(4)and the model values m_{i} at times t_{i} are given by $\begin{array}{ccc}{\mathit{m}}_{\mathit{i}}\mathrm{=}\mathit{m}\mathrm{\left(}{\mathit{t}}_{\mathit{i}}\mathrm{\right)}\mathrm{=}\mathit{A}\mathrm{sin}\mathrm{\left(}\mathit{\omega}{\mathit{t}}_{\mathit{i}}\mathrm{\right)}\mathrm{+}\mathit{B}\mathrm{cos}\mathrm{\left(}\mathit{\omega}{\mathit{t}}_{\mathit{i}}\mathrm{\right)}\mathrm{+}\mathit{C,}& & \end{array}$(5)and Zechmeister & Kürster (2009) present an analytical solution for the bestfit coefficients A–C and, hence, the normalized power p(ω).
For the calculation of the significance of the power recorded at some frequency, we do not follow the approach by Zechmeister & Kürster (2009); rather we perform repeated permutations of the measurements while keeping the observed times the same and carry out a LS analysis of the sopermuted data sets. In this fashion all biases introduced by the temporal sampling pattern of the observations are maintained for the error analysis. To compute the significance we therefore compute predetermined percentiles of the computed normalized power distributions of the permutation at each frequency under investigation.
Fig. 4 Lomb Scargle diagram for the Mount Wilson τ Boo data taken between 1983–1985 (solid line); we also indicate the confidence levels of 0.675, 0.954, 0.9973, 0.999, and 0.9999, respectively. 
In Figs. 4 and 5 we show two examples of our analysis. In each case we performed 1 000 000 random permutations to compute confidence levels of the derived periods. In Fig. 4 one clearly sees a highly significant signal at a period of about four months in addition to a signal at around 85 days, both of which are highly significant well above the 99.99% level. On the other hand, the peak near 160 days is probably a feature introduced by the data sampling, since also the bootstrapped periodograms show a peak near that period. In Fig. 5 again a signal near a period of about four months is apparent, however, the significance is below the formal “3σ” level. The periodogram in Fig. 4 was computed from 283 data points and that in Fig. 5 was computed from 53 data points, indicating the need for dense sampling to reach reliable conclusions on the presence of periodicities in the data.
Fig. 5 Lomb Scargle diagram for the Mount Wilson τ Boo data taken between 1989–1991 (solid line); confidence levels are given as in Fig. 4. 
2.4. Phase jumps
A “phase jump” episode was recorded in the TIGRE S index. In the spring of 2016 the cycle changed its phase by almost 180 degrees, jumping almost from maximum to minimum levels and thereupon rising again. It is obvious to check whether in the Mount Wilson time series there are also signatures of such phase jumps. Clearly, phase jumps can be recognized only in rather densely sampled time series. In the Mount Wilson series these are the years 1983–1986, with a particularly dense sampling available for the year 1984 with 137 different nights. The available data (in 1984) cover six months from January to June and hence cover almost two of the τ Boo activity cycles. In Fig. 6 we show a closeup view of the Mount Wilson data for the observing season 1984 (upper panel) and 1985 (lower panel); for each season we removed the trends and brought all data to the same mean value per season. For each observing season we fit a sine wave with a period of 120 days (solid lines in Fig. 6) to the data, we also propagated the best fit from observing season 1984 into observing season 1985 (dashed line in lower panel of Fig. 6). It appears that some time in the second half of 1984 a “phase jump episode” occurred that was very similar to that observed in the spring of 2016.
Fig. 6 Mount Wilson Sindex time series in observing season 1984 (upper panel) and 1985 (lower panel). The solid lines denote bestfit sine waves; the dashed line in the lower panel denotes the extrapolation of the 1984 bestfit solution into 1985; see text for details. 
Fig. 7 Mount Wilson Sindex time series in the year 1993 (upper panel) and TIGRE Sindex time series in the year 2006 (lower panel); see text for details. 
A second possible “phase jump episode” is visible in the Mount Wilson data taken in 1993. In upper panel of Fig. 7 we plot the Mount Wilson τ Boo Sindex time series; it appears that between days 5105 and 5115, i.e., a tenday period for which there is unfortunately no data, a jump from close to maximum to close to minimum values occurred. For comparison, in the lower panel we show the TIGRE light curve from 2016 (taken from Mittag et al. 2017), where a similar phase jump between day 13 460 and day 13 465 occurred. Therefore activity phase jumps appear to occur on τ Boo, however, the available data coverage precludes us from determining a frequency or occurrence pattern of such events.
3. Discussion and conclusion
We have inspected the data obtained for τ Boo in the context of the Mount Wilson monitoring program for the occurrence of shorter term periodicities. Specifically we investigated to what extent evidence for the fourmonth cycle recently described by Mittag et al. (2017) can be found in this data series. Using both the technique of discrete autocorrelation functions and the technique of Lomb Scargle periodograms, we find strong evidence for the occurrence of periodicities at a period of about 120 days or four months in these data; other shorter periods, especially at 90 days, may also be present in the data. Therefore a cycle with a period of about 120 days seems to be consistently present in τ Boo between at least 1967, when the Mount Wilson observations started and 2016 with the TIGRE observations reported by Mittag et al. (2017). These 49 yr in the life of τ Boo correspond to approximately 147 cycles; in the case of the Sun this number would correspond to more than 1600 yr of continuous cyclic activity. The Mount Wilson time series also provides evidence for phase jumps in the observed Ca activity, however, in general, the temporal sampling of the data is usually too coarse for definite conclusions. There is, naturally, substantial scatter in the data and almost nightly sampling is required to distinguish
the intrinsic scatter from deterministic variations in the data. It is obvious that the dense monitoring of the Mount Wilson S index is a good and observationally very simple method to study the activity cycle of a star such as τ Boo. It would clearly be interesting to accompany these measurements with simultaneous magnetic field measurements and Zeeman Doppler images and check to what extent magnetic topologies change during phase jumps of the activity cycle.
For example, at the web site http://www.sidc.be/silso/datafiles
Acknowledgments
The HK_Project_v1995_NSO data derive from the Mount Wilson Observatory HK Project, which was supported by both public and private funds through the Carnegie Observatories, the Mount Wilson Institute, and the HarvardSmithsonian Center for Astrophysics starting in 1966 and continuing for over 36 yr. These data are the result of the dedicated work of O. Wilson, A. Vaughan, G. Preston, D. Duncan, S. Baliunas, and many others.
References
 Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269 [NASA ADS] [CrossRef] [Google Scholar]
 Baliunas, S. L., Henry, G. W., Donahue, R. A., Fekel, F. C., & Soon, W. H. 1997, ApJ, 474, L119 [NASA ADS] [CrossRef] [Google Scholar]
 Butler, R. P., Marcy, G. W., Williams, E., Hauser, H., & Shirts, P. 1997, ApJ, 474, L115 [NASA ADS] [CrossRef] [Google Scholar]
 Donati, J.F., Moutou, C., Farès, R., et al. 2008, MNRAS, 385, 1179 [NASA ADS] [CrossRef] [Google Scholar]
 Edelson, R. A., & Krolik, J. H. 1988, ApJ, 333, 646 [NASA ADS] [CrossRef] [Google Scholar]
 Fares, R., Donati, J.F., Moutou, C., et al. 2009, MNRAS, 398, 1383 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Fares, R., Moutou, C., Donati, J.F., et al. 2013, MNRAS, 435, 1451 [NASA ADS] [CrossRef] [Google Scholar]
 Lanza, A. F. 2013, A&A, 557, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mengel, M. W., Fares, R., Marsden, S. C., et al. 2016, MNRAS, 459, 4325 [NASA ADS] [CrossRef] [Google Scholar]
 Mittag, M., Robrade, J., Schmitt, J. H. M. M., et al. 2017, A&A, 600, A119 [Google Scholar]
 Schmitt, J. H. M. M., Schröder, K.P., Rauw, G., et al. 2014, Astron. Nachr., 335, 787 [NASA ADS] [CrossRef] [Google Scholar]
 Shkolnik, E., Bohlender, D. A., Walker, G. A. H., & Collier Cameron, A. 2008, ApJ, 676, 628 [NASA ADS] [CrossRef] [Google Scholar]
 Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Overview of available τ Boo Ca II H&K observations obtained in the Mount Wilson (MW) and TIGRE programs.
All Figures
Fig. 1 Mount Wilson Sindex time series of τ Boo. 

In the text 
Fig. 2 Discrete autocorrelation function of Mount Wilson Sindex time series (solid histogram). The dashed histograms denote the maximum values of the autocorrelation functions at the 0.67%, 90%, 99%, and 99.9% level; see text for details. 

In the text 
Fig. 3 Same as Fig. 2, but for solar sun spot record. 

In the text 
Fig. 4 Lomb Scargle diagram for the Mount Wilson τ Boo data taken between 1983–1985 (solid line); we also indicate the confidence levels of 0.675, 0.954, 0.9973, 0.999, and 0.9999, respectively. 

In the text 
Fig. 5 Lomb Scargle diagram for the Mount Wilson τ Boo data taken between 1989–1991 (solid line); confidence levels are given as in Fig. 4. 

In the text 
Fig. 6 Mount Wilson Sindex time series in observing season 1984 (upper panel) and 1985 (lower panel). The solid lines denote bestfit sine waves; the dashed line in the lower panel denotes the extrapolation of the 1984 bestfit solution into 1985; see text for details. 

In the text 
Fig. 7 Mount Wilson Sindex time series in the year 1993 (upper panel) and TIGRE Sindex time series in the year 2006 (lower panel); see text for details. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.