Issue 
A&A
Volume 607, November 2017



Article Number  A87  
Number of page(s)  39  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201630262  
Published online  20 November 2017 
Stellar rotation periods determined from simultaneously measured Ca II H&K and Ca II IRT lines
^{1} Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany
email: mmittag@hs.unihamburg.de
^{2} Department of Astronomy, University of Guanajuato, 36000 Guanajuato, Mexico
Received: 15 December 2016
Accepted: 10 August 2017
Aims. Previous studies have shown that, for latetype stars, activity indicators derived from the Ca ii infraredtriplet (IRT) lines are correlated with the indicators derived from the Ca ii H&K lines. Therefore, the Ca ii IRT lines are in principle usable for activity studies, but they may be less sensitive when measuring the rotation period. Our goal is to determine whether the Ca ii IRT lines are sufficiently sensitive to measure rotation periods and how any Ca ii IRT derived rotation periods compare with periods derived from the “classical” Mount Wilson Sindex.
Methods. To analyse the Ca II IRT lines’ sensitivity and to measure rotation periods, we define an activity index for each of the Ca II IRT lines similar to the Mount Wilson Sindex and perform a period analysis for the lines separately and jointly.
Results. For eleven latetype stars we can measure the rotation periods using the Ca ii IRT indices similar to those found in the Mount Wilson Sindex time series and find that a period derived from all four indices gives the most probable rotation period; we find good agreement for stars with already existing literature values. In a few cases the computed periodograms show a complicated structure with multiple peaks, meaning that formally different periods are derived in different indices. We show that in one case, this is due to data sampling effects and argue that denser cadence sampling is necessary to provide credible evidence for differential rotation. However, our TIGRE data for HD 101501 shows good evidence for the presence of differential rotation.
Key words: stars: activity / stars: chromospheres / stars: latetype / stars: rotation
© ESO, 2017
1. Introduction
In addition to mass, the rotation period of a latetype main sequence star is among the most important stellar parameters. Due to the rotationageactivity connection one can estimate the star’s chromospheric activity level (e.g. Noyes et al. 1984) and its age (Barnes 2007), given its rotation period. Furthermore, the measurement of the rotation period at different epochs can – at least in principle – provide information on the differential rotation of the underlying star (Donahue et al. 1996).
The rotation period of a star can be measured via photometry or by using spectral lines that are sufficiently sensitive to chromospheric activity. Both methods rely on the presence of spots or plage regions located on and thus rotating with the star, imprinting a rotational modulation on the photospheric continuum or chromospheric emission line flux. It is clear that the method only works with an inhomogeneous distribution of spotted and plage regions over the star. However, recent results from the Kepler mission show the success of the photometric method through the measurement of the rotation periods of many thousands of stars (e.g. Reinhold et al. 2013; McQuillan et al. 2014); we highlight that all the stars measured by Kepler are quite faint and that no such measurements for brighter stars are available.
For the spectroscopic method, the Ca ii H&K lines are normally used, since these lines are the most sensitive ones to chromospheric activity of all spectral lines in the optical waveband, and many examples of Ca ii determined rotation periods can be found in the literature (e.g. Noyes et al. 1984; Donahue et al. 1996; Hempelmann et al. 2016). For the very coolest stars, however, the Ca ii H&K lines are harder and harder to observe (at least with smaller telescopes), because these stars become quite faint at the wavelengths of the Ca ii H&K lines.
A possible alternative to the Ca ii H&K lines is the use of the Ca ii infrared triplet (with lines located at 8498, 8542, and 8662 Å, thereafter Ca II IRT). In earlier studies the general correlation between the Ca ii H&K and the Ca II IRT lines was shown (see, e.g. MartínezArnáiz et al. 2011) and previously obtained results show that the Ca II IRT lines are generally usable to also study chromospheric activity. Whether or not it is possible to measure rotation periods using the Ca II IRT lines and, if so, how the periods derived from the Ca II IRT lines compare to the periods derived from the Ca ii H&K lines, remain open questions. Since the chromospheric emission measured in the Ca II IRT and Ca ii H&K lines originates from different chromospheric layers (Montes et al. 2000), periods derived from these lines need not necessarily coincide.
Given that these lines, used as chromospheric proxies, may vary in flux on short time scales, it is clear that in any meaningful comparison study of the Ca ii H&K and Ca II IRT lines one must use strictly simultaneously obtained observations. In this paper we therefore use a sample of such strictly simultaneously carried out monitoring observations. We first present our definition of new Ca ii IRT activity indices, which we use to measure the rotation period. Furthermore, we present a method to create a periodogram from a number of different, but simultaneously measured activity indices. The results of the rotation measurements of the Ca II IRT lines are then compared with the results obtained using the Mount Wilson Sindex (Vaughan et al. 1978; Baliunas et al. 1995) only; a mean rotation period is estimated and compared with literature values. Finally we discuss to what extent there is evidence of differential rotation in our data and present our conclusions.
2. Observation, data reduction, and synthetic spectra
For this investigation we chose twelve stars from our ongoing activity study of stars with spectral types F to K with the TIGRE telescope; here we use only data from the first monitoring program from August 2013 to July 2014, selecting those stars with the highest significance in their periodograms derived from their Ca ii H&K indices alone (cf., Hempelmann et al. 2016). In Table 1, the investigated stars are listed with the number of observations, their visual magnitude and their colour index B−V (taken from Hipparcos catalogue ESA 1997).
List of the observed objects.
The stars were observed with the TIGRE telescope, located at the La Luz Observatory (located in central Mexico). TIGRE is a 1.2 m fully robotic telescope with the fibrefed Échelle spectrograph HEROS (Heidelberg Extended Range Optical Spectrograph). This spectrograph is equipped with two spectral channels (blue & red), a spectral resolution of R ≈ 20 000, and a spectral range from ≈3800 Å to ≈8800 Å with a small gap of 100 Å centred at 5800 Å; a more detailed description of the TIGRE facility has been provided by Schmitt et al. (2014).
The data are reduced with TIGRE/HEROS standard reduction pipeline, which is implemented in IDL and based on the IDL reduction package REDUCE by Piskunov & Valenti (2002). The TIGRE/HEROS standard reduction pipeline is designed as a fully automatic pipeline including wavelength calibration. It contains all necessary reduction steps to reduce Échelle spectra (bias subtraction, order definition, wavelength calibration, spectra extraction, flat fielding); additional information can be found in Hempelmann et al. (2016) and Mittag et al. (2016).
One of the main scientific goals of TIGRE is the longterm monitoring of cool stars. Therefore, an additional part of the TIGRE/HEROS reduction pipeline is the estimation of the instrumental Sindex (thereafter, S_{TIGRE}). This instrumental Sindex can be converted into the “classical” S_{MWO}index with the following linear conversion relation derived by Mittag et al. (2016): (1)During the course of our analysis we also make use of synthetic spectra computed with the radiation transfer code PHOENIX version 16 (Hauschildt et al. 1999). The PHOENIX code provides atmosphere calculations for novae and supernovae as well as for stars and brown dwarfs. Assuming local thermal equilibrium (LTE), Husser et al. (2013) compute the most recent stellar atmosphere library based on the PHOENIX code. The stellar model spectra used in this work were calculated for the mean stellar parameters listed in Mittag et al. (2013, Table 2); they are based on the models computed by Husser et al. (2013), but the elements H, He, Ca i–iii, Na i–II, and Mg i–iii are treated in nonLTE, which results in improved descriptions of these lines.
3. Ca II IRT based activity indices
3.1. Index definitions
To characterise the activity observed in the Ca II IRT lines, the subtraction method (e.g. MartínezArnáiz et al. 2011) could be used, where the reference spectrum of an inactive star is subtracted from the spectrum of the star of interest. The reference spectrum must come from a star as similar in stellar parameters, such as B−V, log g, and metallicity, as is possible to ensure minimal differences in the photospheric spectra. Naturally, finding such a suitable reference spectrum may pose a problem.
A simpler and more practical method is the definition of an activity index in analogy to the “classical” Mount Wilson Sindex by determining the ratio between the counts in the line core of each of the Ca II IRT lines and the counts in a reference spectral range. The great advantage of this purely empirical approach is its simplicity with an index directly obtained from the data without any reference to a synthetic or comparison inactive spectrum.
Fig. 1
A TIGRE spectrum of the Sun (moon light) of the Ca II IRT lines with the bandpasses chosen for the definition of the S_{IRT}indices: the top panel shows the spectral regions used for the S_{IRT_8498}index, the middle panel those used for the S_{IRT_8542}index and the bottom panel those used for the S_{IRT_8662}index. 
We then define our S_{IRT}indices for the Ca II IRT lines using the following equations: and (4)where N is the number of counts in the bandpasses explained below and illustrated in Fig. 1. For our S_{IRT}indices we always use an 1 Å bandpass centred in the line core to estimate the counts in the Ca II IRT lines. In Fig. 1, the three Ca II IRT lines with the bandpasses of the respective S_{IRT}index are shown. Some care in the definition of the locations of the reference bands is required since in the spectral range of the Ca II IRT some telluric lines (mostly H_{2}O lines) are located. The line intensity of the telluric lines is timevariable, therefore this must be considered in the definition of the used line continua, and we have selected spectral ranges for the line continua without any or very low telluric line contamination. Since the considered spectral range contains a large number of stronger telluric lines, this resulted in only 2 Å wide bandpasses used for the reference continua for the Ca II IRT lines at 8498 and 8542 Å; the defined continua for the Ca II IRT line at 8498 Å are centred at 8490 and 8509 Å and for the Ca II IRT line at 8542, 8530, and 8566 Å. For the Ca II IRT line at 8662 Å, the line continua are defined as a 10 Å wide bandpass centred at 8641 and 8690 Å, respectively. The three lines in the bandpass centred at 8690 Å are stellar lines of the elements Si i and Fe i. We note that the continuum reference bands are not symmetric with respect to the respective lines.
3.2. Index renormalisation
Because of the asymmetric placement of the quasicontinuum reference bands we must use the “true” spectral trend and slope for the considered spectral range (8450–8720 Å) to obtain S_{IRT}indices consistent with the respective flux ratios. We adapt the method described by Hall & Lockwood (1995) to renormalise the Ca II H&K region for measuring the instrumental Mount Wilson Svalue for the spectral range of the Ca II IRT. To estimate the spectral slope in the Ca II IRT region, we use synthetic PHOENIX spectra (see Sect. 2), calculated for main sequence stars in the B −V range from 0.44 to 1.6, as follows: First, we normalise the computed PHOENIX spectrum by its maximal flux value, then we select seven supporting points (at 8456, 8484, 8523, 8570, 8630, 8697, and 8715 Å) for the continuum (x,y) and fit a linear function with the approach y = A + Bx, where y is the normalised flux and x the wavelength. The resulting fit coefficients – the constant A and the slope B – are then functions of the colour B−V and are plotted against B−V in Figs. 2 and 3.
Fig. 2
Upper panel: best fit coefficients of the constant A of the linear fit vs. B−V. The solid line represents the determined relation. Lower panel: the residuals. 
Fig. 3
Upper panel: best fit coefficients of the slope B of the linear fit multiplied with 10^{5} vs. B−V. The solid line represents the determined relation. Lower panel: the residuals. 
The mean Sindices of the Ca II H&K and Ca II IRT lines and the standard deviation (σ).
Finally, we empirically determine a relation describing the dependence of the constant A and the slope B on the B−V–colour to calculate a correction function corr(λ,B−V) to renormalise the whole Ca II IRT range. For the correction function in the spectral range 8450–8720 Å, we obtain the following relations: (5)where x = λ−8450 Å and the relations for the constant and slope are: where z = (B−V)−0.4.
We now briefly describe our estimation of the S_{IRT}indices using the Ca II IRT line at 8498 Å as an example. First, the wavelength is corrected with the barycentric velocity shift. Next, the Ca II IRT line region is continuumnormalised. Then we measure the radial velocity by a crosscorrelation with a synthetic PHOENIX spectrum, taking into account the rotational linebroadening by vsin(i) and the spectral resolution of the spectrograph (R ≈ 20 000). After the radial velocity correction, the “true” tilt of the spectrum is calculated with the correction function Eq. (5) and multiplied by the continuumnormalised spectrum. From this spectrum we calculate the S_{IRT_8498}index using Eq. (2). The counts in the bandpasses are integrated and the ratio is calculated; to estimate the errors of the individual bandpasses, we perform a MonteCarlo simulation for each bandpass, taking into account the signaltonoise ratio in each single pixel. The standard deviation of the integral variations is used as the error of the corresponding bandpass. Finally the error of the S_{IRT_8498}index is estimated by error propagation of the errors of the single bandpasses.
3.3. Ca II based indices and correlations
For the sample stars considered in this paper, the mean S_{MWO} and S_{IRT} indices and the corresponding standard deviations are listed in Table 2. Furthermore, the S_{IRT}indices vs. the S_{MWO}indices are plotted for each sample star in Figs. A.1–A.12. An example is shown in Fig. 4 for the S_{IRT_8498} vs. S_{MWO} indices measured by TIGRE for the star HD 131977. Obviously, the two indices appear highly correlated, and in order to quantify the correlation between the S_{MWO} and S_{IRT}indices we calculate Spearsman’s ρ correlation coefficients for all of our sample stars and list the results in Table 3. As is clear from Table 3, for most stars, the S_{IRT} and S_{MWO}indices are clearly correlated, and only for a few stars, do we find little or no correlation between the various indices. In Fig. 5 we plot the derived values of ρ between the three Ca II IRT lines and Ca II H&K lines as a function of the measured dispersion in the S_{MWO}indices and find that stars with a larger spread in the S_{MWO}index tend to show a larger correlation coefficient. The significance of the correlations between σ(S_{MWO}) and ρ(S_{MWO}S_{IRT_8498}), and ρ(S_{MWO}S_{IRT_8542}) and ρ(S_{MWO}S_{IRT_8662}) are 99.9%, 99.7%, and 99.3%. There are two possible explanations for the trends. First, the Ca II IRT lines may have a lower sensitivity to stellar activity. Second, noise may be a problem, hiding, especially for the lowest dispersion values, any underlying correlation.
The correlation coefficient (Spearsman’s ρ) between S_{MWO} and S_{IRT}indices time series.
Fig. 4
Correlation between S_{IRT_8498} and S_{MWO} for HD 131977 with a Spearman’s ρ of 0.924. 
Fig. 5
Spearman’s ρ correlation coefficients of the correlation between S_{MWO} and S_{IRT}indices vs. the standard deviation (σ) of the S_{MWO} time series. 
4. Period and significance estimation
4.1. Detrending
Before analysing our time series for periodic variations, we have to determine and remove possible trends in the data; such trends could be caused, for example, by longterm trends in the activity of the sample stars. For our study, the precise nature of these trends is not relevant, however, detrending is important for the period estimation, since an unaccounted trend can attenuate the power derived in the periodogram. The common methods to search for periodicities are the LombScargle periodogram (Lomb 1976; Scargle 1982; Horne & Baliunas 1986) or the generalized LombScargle periodogram (thereafter GLS; Zechmeister & Kürster 2009), which assume a purely sinusoidal wave form, and the GLS by Zechmeister & Kürster (2009), which additionally takes into account a constant offset, but no systematic trend.
Fig. 6
Example for the estimation of the polynomial order for the detrending of the S_{MWO} time series of HD 43162: the upper panel shows the maximal power p(ω), the panel below the results of the Ftest, where the significance of the scatter reduction is compared to the residuals of the nondetrended time series, the third panel shows the significance of the scatter reduction compared to detrending with the polynomial fit of order N1 and in the lowest panel we plot the estimated periods; the dotted lines show the formal 1σ limit. 
We perform detrending by using polynomials of varying order. To decide which order of polynomial is appropriate, we detrend the time series with polynomials of different order, analyse the periodic variations via GLS, and estimate the period with the largest peak height. The thus determined period is then removed from the data via a sinusoidal wave fit. To estimate the significance of the trend, the Ftest is used to estimate the significance of the reduction of the scatter in the residuals. The criterion for finding the most suitable order of the polynomial fit is the height of the largest peak (maximal power (p(ω))) and hence the significance of the corresponding period. Furthermore, the results of the Ftest are used to check the significance of the used polynomial fit and so to find the optimal order of the polynomial fit. Here, the significance limit of the scatter reduction is set to “1σ” because the intrinsic scattering is relatively large, meaning that a higher limit is inappropriate.
As an example of our procedure we consider the case of HD 43162 (see Figs. 6 and A.4). In the upper panel, the maximal peak heights (Max p(ω)) are shown as a function of the polynomial order and it is apparent that the peak height increases significantly with a detrending polynomial of first order. The following panel shows the significance of the scatter reduction compared to the residuals of the nondetrended time series. In the third panel of Fig. 6, the significance of the scattering reduction compared to the detrending with the polynomial fit of order N−1 are shown to check whether the significance changes exceed the 1σ level. It is obvious that a fit with a polynomial of order larger than 1 is not warranted. In the lower panel, the periods vs. polynomial order are depicted to show the dependence of the derived period on the order of the polynomial.
Finally, we point out that the detrending of the investigated time series is a somewhat controversial issue, because the details of the applied detrending procedures can influence the resulting periodograms. To quantify the influence of the detrending polynomial on the period, we calculate the ratio between the standard deviation and mean of the periods for the different polynomial fits (see the listing in Table 4). The periods do not show any significant dependencies on the used polynomial orders and are stable to within a few percent, except for the case of HD 75732, where no significant rotation period was found in the Ca II IRT lines. Furthermore, the periods vs. polynomial orders are shown in the corresponding panels in Figs. A.1–A.12 to visualise the general problem of detrending.
Ratio between the standard deviations of the periods and the mean period obtained from the different polynomial orders for the detrending.
4.2. Period determination
For the search for periodic variations in the time series of the Ca ii H&K and Ca ii IRT lines, we carry out a Fourier analysis using the GLS method by Zechmeister & Kürster (2009). Given M data values y_{j} at times t_{j} and uncertainties σ_{j}, the formalism computes a χ^{2} value with the approach y_{j} = Asinωt_{j} + Bcosωt_{j} + C, with constants A, B and C determined from the data, to generate a periodogram as a function of the chosen frequency ω. The periodogram is computed from the following expression (Zechmeister & Kürster 2009); (8)where is a value for the dispersion around the mean and is used for the normalisation of χ^{2}(ω); the detailed formulae to compute p(ω) can be found in Zechmeister & Kürster (2009).
In our case, we are dealing not only with a single series of data y_{j} at times t_{j}, but rather we are dealing with a total of four index time series y_{i,j}, where i = 1,4 and j = 1,M, each representing one of our Ca ii index time series. If one is not interested in the phase of the different time series we can “combine” the periodograms in the following way: Zechmeister & Kürster (2009) show that the normalised power spectrum p(ω) (with the cross terms set to zero by suitably shifting the times t_{j}) can be written as (Eq. (5) in their paper) (9)where the nomenclature is the same as in Zechmeister & Kürster (2009). The two sums in Eq. (9) extend over all the data points in the given data set at some frequency ω. In our case, several data sets provide information on the power at frequency ω. We therefore consider the mean power (averaged over the N data sets) through (10)where the quantities subscripted with the index i refer to the ith data set. Since the data sets are independent (the phases do not need to agree) we can use the formulae given in the Appendix of Zechmeister & Kürster (2009) for each data set and find (11)which is a simple averaged periodogram.
In this fashion we can “average” the periodograms from our four line indices into a single, mean periodogram, an example of which is shown in Fig. 7. Here, the averaged periodogram and the original periodograms of the four time series for HD 43162 are shown.
Fig. 7
Averaged periodogram and individual periodograms from the S_{MWO} and the S_{IRT} time series for HD 43162; the solid line depicts the averaged periodogram, the dotted line the S_{MWO}, the dashed line the S_{IRT_8498}, the dasheddotted line the S_{IRT_8542} and the dashedtripledotted line the S_{IRT_8662} periodograms. 
4.3. Period errors, significance, and FAP estimation
To estimate the error on the determined period we use the equation given by Baliunas et al. (1995), (12)where P is the period, σ_{n} the standard deviation of the residuals (O–C), T the total length of the observation interval, A the amplitude of the signal, and N the total number of data points.
Fig. 8
Significance estimation of the S_{MWO} time series for HD 43162: Upper panel: logarithm of the probability (Prob(p>p_{0})) vs. logarithm of one minus significance level (1−p_{0}). The solid line shows the linear fit and the dotted line the quadratic fit. Lower panel: the residuals of the linear fit are depicted as crosses and the quadratic fit as asterisks. 
For the most probable rotation period (P_{M}) obtained from the averaged periodogram, the error is estimated in the same fashion. Here, the amplitude and the σ_{n} for each time series are estimated with the period found in the averaged periodogram and the values of ΔP for each time series are estimated. Finally, the ΔP values of each time series (ΔP_{i}) are combined to the error of the most probable rotation period (ΔP_{M}) with the formalism to estimate the error of the weighted mean: (13)To estimate the probability levels in the periodogram, a bootstrap method is used, permuting the recorded data 10^{5} times, while maintaining the times of the measurements. With this method, we cannot estimate the exact probability of a peak in the periodogram, but we can estimate probability limits. Given the probabilities and the corresponding significance levels, we can derive a general relation between the peak height and the probability in order to estimate the peak probability. Therefore, the approach for the estimation of the probabilities is used from Zechmeister & Kürster (2009): (14)where p_{0} is the peak height and N the number of data points. The exponent is a normalisation factor for the periodogram to convert the power into a probability. This exponent can be substituted so that we obtain (15)or (16)The latter equation describes a linear relation between the logarithmic probability and the logarithm of the peak height with the normalisation factor as the slope.
Fig. 9
Significance estimation of the signal in the averaged periodogram for HD 111395: Upper panel: logarithm of the probability (Prob(p>p_{0})) vs. logarithm of one minus significance level (1−p_{0}). The solid line shows the linear fit and the dotted line the quadratic fit. Lower panel: the residuals of the linear fit are depicted as crosses and the quadratic fit as asterisks. 
Periods found using the S_{MWO} and S_{IRT} time series with the corresponding FAP.
Fig. 10
Logarithm of the estimated FAPs of the derived periods in the individual time series vs. objects. The red bars show the log(FAP) values for the period derived from the S_{MWO} time series, the blue bars that from the S_{IRT_8498} time series, the yellow bars that from the S_{IRT_8542} and the green bars that from the S_{IRT_8662} time series. 
In Fig. 8 we present an example of our procedures and plot the logarithmic probabilities against the logarithm of the peak heights for the bootstrapped GLS periodograms for HD 43162; we find that a linear relation describes the general trend of the relation between the logarithmic probability and logarithm of the peak height. In this particular case we obtain a slope of 11.7 which differs slightly from the calculated value for the exponent of 12.5 (using N = 28, see Eq. (14)). Furthermore, a visual inspection of the residuals shows a systematic deviation from a linear curve. We therefore add a quadratic correction term to Eq. (16) and write (17)which provides a better description of the relation between the logarithmic probability and logarithm of the peak height. The standard deviation of the residuals for the linear fit is 0.019 and for the quadratic fit is 0.006. The Ftest gives a significance of ≈99.5% for the scatter reduction with the quadratic fit compared to the linear fit. Another example for the quadratic fit is shown in Fig. 9, where we show the significance estimation for the averaged periodogram of the four time series for HD 111395; it is obvious that the quadratic fit describes this relationship better than the linear fit: the standard deviation of the residuals is ≈0.11 for the linear fit and 0.02 for the quadratic fit; and with the Ftest, we obtain a significance >99.9% for the scatter reduction with the quadratic fit compared to the linear fit. Needless to say, in some cases, a linear fit is possible, but for consistency we use the quadratic fit.
In comparison to Eq. (14), the normalisation factor is estimated independently from the number of data points and we obtain a more accurate peak probability.
With Eq. (17), the probability of a peak can be calculated, but to estimate the false alarm probability (thereafter, FAP) of the peak one needs the number of independent frequencies. To estimate these frequencies we use the method described by Zechmeister & Kürster (2009) and assume the number of independent frequencies , where Δf (f_{max}−f_{min}) is the frequency range and δf the frequency resolution. The frequency resolution is estimated as the full width at half maximum (FWHM) of a peak in the periodogram. The frequency range corresponds to the highest reasonable frequency (f_{max}) because the time series are long enough to assume f_{max} ≫ f_{min} and so that applies Δf = f_{max} (Zechmeister & Kürster 2009). To estimate the highest reasonable frequency, the time spans between consecutive observations (ΔT) are calculated. and the most frequent value of ΔT is multiplied by two to obtain the smallest reasonable period and so the highest reasonable frequency. This definition is based on the definition of the Nyquist frequency for evenly sampled time series.
Finally, following Zechmeister & Kürster (2009) the FAP is estimated as (18)where M is the number of independent frequencies and Prob(p>p_{0}) is estimated with Eq. (17).
5. Results of the rotation period measurement
In the following, we present our results of the measurement of the stellar rotation period for the Ca ii H&K lines and the Ca ii IRT lines as described in Sect. 4.
5.1. Time series for individual lines
Fig. 11
Comparison of the periods derived from the S_{IRT}indices time series vs. the periods derived from the S_{MWO} time series (in the upper panels) with the identity as solid line. The lower panels show the relative difference (in percent) between the periods derived in the S_{IRT}indices and S_{MWO} time series vs. the period obtained from the S_{MWO} time series with a 3σ error for the period ratios. The periods with a FAP less than 4.55 × 10^{2} (2σ) are labelled with a plus symbol and periods with a FAP greater than 4.55 × 10^{2} with an asterisk. 
The results from the analysis of the time series using individual lines are listed in Table 5. Each step in the process of rotation period determination for each of the 12 studied stars is documented in different graphs found in Figs. A.1 to A.12, including the raw time series of the indices with the used fit for the detrending, determined periods as a function of the order of the detrending polynomial, periodograms, and phasefolded time series of the detrended indices.
For 11 out of 12 sample objects, we find a significant rotation period with a FAP lower than 2.7 × 10^{3} (3σ) in at least one of the Ca ii IRT lines and a correlation between the indices; see Figs. A.1 to A.12. Only for HD 75732 do we find no period using the Ca ii IRT lines and also no correlation between the indices; see Fig. A.6. Moreover, the derived period in the S_{MWO} time series of HD 75732 also strongly depends on the chosen detrending procedure; see Fig. A.6. The resulting period range is 40.5 to 37.4 days for detrending with a polynomial of the first, second, and third order. The reason for this behaviour is a longterm trend and a periodic variation with a period of ≈100 days. With the used polynomial fit of the third order, the longterm trend and the 100 days variation is removed and so we obtain a 37.4 day period for HD 75732.
To compare the FAPs of the periods found in the S_{MWO} and in the S_{IRT} time series, the logarithm of the FAPs for each time series are shown as a bar plot vs. the objects in Fig. 10, except for HD 75732. Figure 10 shows that in only four cases is the FAP of the period derived from the S_{MWO} time series the smallest FAP of all time series considered.
In Fig. 11 the periods obtained from the S_{MWO} and the S_{IRT}indices are plotted. These periods follow the identity with a small deviation for some objects. To quantify this deviation, we calculate the relative difference between the period found from the S_{IRT} and S_{MWO}indices (see Fig. 11); the difference between the periods is lower than 9%. The majority of the periods are identical to within the 3σ errors. However, a few periods are clearly different. These variations of the periods are of the order of the variations of the minimal and maximal rotation periods to the mean rotation periods listed in Donahue et al. (1996, Table 1); for a possible reason see our discussion in Sect. 6.
5.2. Most probable rotation period
For eleven stars, we find a rotation period in the S_{MWO} time series and in the individual S_{IRT} time series. To estimate the most probable rotation period and the corresponding FAP, the periodograms are combined with the method described in Sect. 4 and shown in Figs. A.1 to A.12. Here, we use all four time series even if one time series shows a weak signal in the periodogram. The results for the eleven stars are listed in Table 6. For completeness, the determined period for HD 75732 is also listed in Table 6 although
Literature period, most probable rotation period and the corresponding FAP.
this period is only apparent in the S_{MWO} time series.
Finally, we compare our most probable periods with periods found in the literature (cf., see Fig. 12 and Table 6), which is the case for eight stars. For six stars, our periods agree with the literature values to within 5%. For HD 72905 and for HD 75732 the difference between our period and the literature value is ≈7% and ≈11%, respectively. The differences between ours and the literature values may be explained with the possible differential rotation of the star, as is shown by stars listed in Donahue et al. (1996, Table 1). Furthermore, we have no information about the uncertainty of the periods for the stars with a literature period from Baliunas et al. (1998) and Wright et al. (2011). Therefore, we cannot assess the true difference between these values. For HD 75732, the large discrepancy between ours and the literature value may be explained by a detrending effect; see Sect. 5.1.
Fig. 12
Upper plot: comparison of our most probable rotation periods with the literature values with the identity as solid line. Lower panel: the relative difference in percent between our most probable rotation periods and the literature values with 3σ error for the period ratios. 
6. Evidence for differential rotation?
The measurement of the rotation period using spectral line variations is only possible if the star exhibits a stable distribution of spots and/or active regions on a timescale comparable with the length of the timeseries with a lifetime longer than the rotation period itself. These stellar active regions can evolve with time so that one may measure different periods at different epochs if the star rotates differentially and the active regions dominating the overall activity signal change their latitudes. This effect can be directly observed for the Sun (e.g. Hempelmann & Donahue 1997), while for stars the situation is far more elusive, since their surfaces cannot be spatially resolved with activity indices measurements. Nevertheless, there have been several studies of differential rotation for stars using this method; for example, Donahue et al. (1996) interpret the variations of the rotation periods measured with their S_{MWO} time series as a measurement of differential rotation. We note that in our sample of twelve stars (see Table 5), four stars are also included in the study of Donahue et al. (1996); namely HD 35296, HD 39587, HD 72905 and HD 101501.
In our case, we measure the rotation period only at one epoch, which should lead to the same rotation period measured in the S_{MWO} and in the S_{IRT} time series. However, for the stars HD 35296, HD 39587, HD 97334, and HD 101501, we actually find different rotation periods using the different indices. An inspection of the periodograms of the stars HD 35296 and HD 39587 (shown in Figs. A.1 and A.3) shows a complicated structure with multiple peaks, but very similar for the different activity indices – the only difference being the location of the highest peak. Another case, possibly indicative of differential rotation, is that of HD 101501 (Fig. A.8), which shows very broad peaks in the periodograms of the different indices, some of which even show some structure; for example, the S_{IRT_8662} time series shows a small bump at a period of ≈17 days, which corresponds to the main peak found in the S_{MWO} time series. Our rotation periods for the stars HD 35296 and HD 39587 determined with the individual indices fall within the minimum and maximum period found by Donahue et al. (1996) for their different data epochs, while for the star HD 101501 we find a slightly lower period than the minimum period by Donahue et al. (1996).
To check the significance of the different periods we carry out a “prewhitening” analysis by removing our most probable period (see Table 6, Col. 3) from the time series and repeating the GLS analysis on these prewhitened data for HD 35296, HD 39587, HD 97334 and HD 101501. Here, we test whether the peak structures are caused by aliasing or not. In the following we discuss the results for these four stars.
Fig. 13
Periodograms for HD 39587 of the S_{MWO} and S_{IRT_8542} time series without the removal of the 5.36 day period (depicted as solid line) and after prewhitening (dashed line). 
Fig. 14
Periodograms for HD 39587 of the S_{MWO} and S_{IRT_8542} time series without the removal of the 5.13 day period (depicted as solid line) and after prewhitening (dashed line). 
Fig. 15
Periodograms for HD 39587 of the S_{MWO} and S_{IRT_8542} time series without the removal of the 5.6 day period (depicted as solid line) and after prewhitening (dashed line). 
HD 35296: the results for the star HD 35296 are shown in Figs. A.13 and A.14 for the four indices. In the periodograms of the prewhitened time series, no significant power peaks remain, which suggests that the complicated power distribution of the “raw” time series can be explained by aliasing effects.
HD 39587: in Fig. A.15, the periodograms of the prewhitening analysis of the 5.36 days removal are shown for the four indices and additionally in Fig. 13 for the S_{MWO} and S_{IRT_8542} time series. These figures show that the peaks are generally removed, except for the S_{IRT_8542} time series. Here, the peak at 5.6 day is significantly, but not completely reduced. At this point, it is worth mentioning that the most probable period is the mean of the 5.13 and 5.6 days and therefore, these signals could be averaged out so that in the prewhitening analysis the peaks are also removed. Hence, the prewhitening analysis is also performed for the periods found in the S_{MWO} time series (5.13 days; see Figs. 14 and A.16) and in the S_{IRT_8542} time series (5.6 days; see Figs. 15 and A.17). We find in the case of the 5.13 days removal that the 5.6 days peak is still present in the periodograms for the O–C data especially for the S_{MWO} and S_{IRT_8542} time series; see Fig. 14. In the other case, the 5.6 day removal, the peak at the period of 5.13 days in the (O–C) periodogram for the S_{MWO} time series is clearly and significantly visible. Therefore, we conclude that the periods 5.13 and 5.6 days should be real but we cannot completely exclude aliasing.
Fig. 16
Periodograms for HD 97334 of the S_{MWO} and S_{IRT_8498} time series without removal of the 7.92day period (depicted as solid line) and after prewhitening (dashed line). 
Fig. 17
Periodograms for HD 97334 of the S_{MWO} and S_{IRT_8498} time series without removal of the 8.26day period (depicted as solid line) and after prewhitening (dashed line). 
HD 97334: for the prewhitening analysis, only the periods found in the S_{MWO} and S_{IRT_8498} time series are used because the most probable period is similar to the period found in the S_{MWO} time series. The results for all four time series are shown in Figs. A.18 and A.19. First, the prewhitening analysis is performed for the period found in the S_{MWO} times series (7.92 day; see Figs. A.18 and 16). In this case, we find that the peak at ≈8.3 days is still visible in the periodogram for the O–C data for the S_{IRT_8542} time series; see Fig. 16. For the prewhitening analysis with the 8.26 days period (see Figs. A.19 and 17) we find a peak at ≈7.8 days in the periodogram for the O–C data for the S_{MWO} time series, which is close to the 7.92day period. Hence, we conclude that the difference between the periods 7.92 and 8.26 days could be real, but just like the case of HD 39587, we cannot fully exclude aliasing.
HD 101501: the periodograms of the prewhitening analysis are shown for the four indices in Fig. A.20. For the S_{MWO} time series a clear peak is visible in the periodogram after the removal of the most probable period. This peak at ≈17.5 days is slightly shifted to the peak before the removal but it is clearly there. Also in the S_{IRT_8662} time series this peak (at ≈18 days) is visible; see Fig. 18.
Fig. 18
Periodograms for HD 101501 of the S_{MWO} and S_{IRT_8662} time series without the removal of the 16.06 day period is depicted with the solid line and the data after the removal (O–C) with a dashed line. 
We conclude that this period difference between S_{MWO} time series and S_{IRT_8542} and S_{IRT_8662} time series cannot be a sampling effect and therefore tentatively interpret this difference as evidence of differential rotation.
We note in this context that the excess fluxes for the Ca II H&K and the Ca II IRT lines are formed in different chromospheric layers (Montes et al. 2000). Therefore, one cannot exclude the possibility that in some cases the location of the dominant active regions, where the Ca II H&K and the Ca II IRT lines are formed, are different and hence lead to different periods in the different lines, as suggested by the result of the prewhitening analysis for HD 101501. On the other hand, our analysis for HD 35296 (shown in Fig. A.13) demonstrates that the period differences appearing in the GLS periodograms are a sampling effect and we conclude that here more observations with a denser sampling cadence are required to arrive at conclusive results in this area of study.
7. Discussion and conclusions
In this study we measure rotation periods of twelve stars using the three Ca II IRT lines and the simultaneously measured Ca II H&K lines. To do so, we introduce S_{IRT}indices to characterise the excess emission in the Ca ii IRT line cores and use the common S_{MWO}index for the Ca II H&K lines. For eleven stars we find a significant rotation period at least in one of the three S_{IRT} indices; the significance of these periods determined from their S_{IRT} time series is comparable to the significance of the periods found in the simultaneously measured S_{MWO} time series, and in a few cases is even larger. Thus we conclude that the Ca II IRT lines can indeed be used for the determination of the rotation period, and for the most robust period determinations one should use all the individual lines; we present a procedure to properly “add” the LombScargle periodograms obtained for different activity indices. For those stars with previous period measurements available in the literature (eight stars) we find agreement between our newly measured values and the literature values to within better than 7%, except HD 75732 with ≈11%.
In general, the rotation periods measured at the same epoch but in different lines are expected to be identical to within the errors and this behaviour is found in our current and previous study of HD 131977 (Mittag et al. 2015) and for six additional stars in the program stars studied here. Five stars show a somewhat different behaviour: one object shows a clear period in the Ca II H&K lines, but in all Ca II IRT lines no periodic signal is visible in the periodogram; interestingly, this is the object with the lowest activity in our sample. In the other four cases, the derived GLS periodograms show a rather complicated structure. A prewhitening for HD 35296 at the most probable period and the 3.502day period for all four lines reduces the remaining periodograms to noise, indicating that the structures in the periodogram are due to data sampling, which we intend to improve in the future. On the other hand the result of prewhitening the data for HD 101501 shows that the determined periods are clearly different and are not likely caused by aliasing.
In summary, we conclude that the Ca II IRT lines are usable as activity indicators and that the Ca II IRT lines are sufficiently sensitive to measure the rotation period of a star. With the possibility to use the Ca II IRT lines as activity indicators, an additional chromospheric layer (to that where the Ca II H&K lines are formed) becomes accessible to observations.
Acknowledgments
We are grateful to our referee Dr. M. Kürster for his helpful and very constructive referee report, which led to substantial improvement of this paper. This research has made valuable use of the VizieR catalogue access tool, CDS, Strasbourg, France. The original description of the VizieR service was published in A&AS 143, 23. We are grateful for the financial support received through the DFG project No. SCHM 1032/491 and we further acknowledge travel support by bilateral grants CONACyTDFG (No.192334) and CONACyTDAAD (No. 207772) in the past 2 years, as well as by the CONACyT mobility grant (No. 207662).
References
 Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269 [NASA ADS] [CrossRef] [Google Scholar]
 Baliunas, S. L., Donahue, R. A., Soon, W., & Henry, G. W. 1998, in Cool Stars, Stellar Systems, and the Sun, eds. R. A. Donahue, & J. A. Bookbinder, ASP Conf. Ser., 154, 153 [Google Scholar]
 Barnes, S. A. 2007, ApJ, 669, 1167 [NASA ADS] [CrossRef] [Google Scholar]
 Donahue, R. A., Saar, S. H., & Baliunas, S. L. 1996, ApJ, 466, 384 [NASA ADS] [CrossRef] [Google Scholar]
 ESA, ed. 1997, The HIPPARCOS and TYCHO catalogues. Astrometric and photometric star catalogues derived from the ESA HIPPARCOS Space Astrometry Mission, ESA SP, 1200 [Google Scholar]
 Hall, J. C., & Lockwood, G. W. 1995, ApJ, 438, 404 [NASA ADS] [CrossRef] [Google Scholar]
 Hauschildt, P. H., Allard, F., & Baron, E. 1999, ApJ, 512, 377 [NASA ADS] [CrossRef] [Google Scholar]
 Hempelmann, A., & Donahue, R. A. 1997, A&A, 322, 835 [NASA ADS] [Google Scholar]
 Hempelmann, A., Mittag, M., GonzalezPerez, J. N., et al. 2016, A&A, 586, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Horne, J. H., & Baliunas, S. L. 1986, ApJ, 302, 757 [NASA ADS] [CrossRef] [Google Scholar]
 Husser, T.O., Wendevon Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lomb, N. R. 1976, Ap&SS, 39, 447 [NASA ADS] [CrossRef] [Google Scholar]
 MartínezArnáiz, R., LópezSantiago, J., CrespoChacón, I., & Montes, D. 2011, MNRAS, 414, 2629 [NASA ADS] [CrossRef] [Google Scholar]
 McQuillan, A., Mazeh, T., & Aigrain, S. 2014, ApJS, 211, 24 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Mittag, M., Schmitt, J. H. M. M., & Schröder, K.P. 2013, A&A, 549, A117 [Google Scholar]
 Mittag, M., Hempelmann, A., GonzalezPerez, J. N., & Schmitt, J. H. M. M. 2015, in Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, 18, eds. G. T. van Belle, & H. C. Harris, 549 [Google Scholar]
 Mittag, M., Schröder, K.P., Hempelmann, A., GonzálezPérez, J. N., & Schmitt, J. H. M. M. 2016, A&A, 591, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Montes, D., FernándezFigueroa, M. J., De Castro, E., et al. 2000, A&AS, 146, 103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763 [NASA ADS] [CrossRef] [Google Scholar]
 Piskunov, N. E., & Valenti, J. A. 2002, A&A, 385, 1095 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reinhold, T., Reiners, A., & Basri, G. 2013, A&A, 560, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Scargle, J. D. 1982, ApJ, 263, 835 [NASA ADS] [CrossRef] [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]
 Vaughan, A. H., Preston, G. W., & Wilson, O. C. 1978, PASP, 90, 267 [NASA ADS] [CrossRef] [Google Scholar]
 Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Single time series
To better illustrate our time series analysis of the different indices we show in Figs. A.1 to A.12 some graphs that represent different steps and concepts of our analysis procedure for each star. The upper four panels show the power of the highest peak, the significance of the scatter reduction compared to the nondetrended time series, the significance of the scatter reduction compared to the detrended time series with the polynomial fit of order N−1, and the period as a function of the order of the polynomial fit for all four used indices. The significance of the scatter reduction is estimated with the Ftest. A detail description can be found in Sect. 4.1. The following six middle panels show the correlation betweenthe new S_{IRT}indices and the S_{MWO}index for the original data (upper three panels) and for the data detrended with the used polynomial fit (lower three panels). For HD 101501, all time series are not detrended so that only the correlations of the original data show. After a page break, the upper four panels show the original time series of all four indices, with the black line denoting the best polynomial fit used for detrending. The middle two panels show: (left) the obtained periodogram for all four indices with the mean probability level (Prob, see Sect. 4.3) 1σ, 2σ, 3σ and 4σ; and (right) detrended and phasefolded time series for all four indices with the black line denoting the bestfitting sinusoidal. More information about the period estimation can be found in Sect. 4.2. The lower two panels show the window function of the time series on the left side and the mean periodogram of the four time series on the right side.
Fig. A.1
Plots for the time series analysis of HD 35296 and, for the detrending, the polynomial orders of 3, 1, 1, and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542}, and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.1
continued. 
Fig. A.2
Plots for the time series analysis of HD 37394, and, for the detrending, the polynomial orders of 1, 0, 0 and 1 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.2
continued. 
Fig. A.3
Plots for the time series analysis of HD 39587 and, for the detrending, the polynomial orders of 2, 0, 1 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.3
continued. 
Fig. A.4
Plots for the time series analysis of HD 43162 and, for the detrending, the polynomial orders of 1, 1, 1 and 1 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.4
continued. 
Fig. A.5
Plots for the time series analysis of HD 72905 and, for the detrending, the polynomial orders of 0, 0, 1 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.5
continued. 
Fig. A.6
Plots for the time series analysis of HD 75732 and, for the detrending, the polynomial orders of 3, 0, 0 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.6
continued. 
Fig. A.7
Plots for the time series analysis of HD 97334 and, for the detrending, the polynomial orders of 2, 2, 0 and 1 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.7
continued. 
Fig. A.8
Plots for the time series analysis of HD 101501 and, for the detrending, the polynomial orders of 0, 0, 0 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.8
continued. 
Fig. A.9
Plots for the time series analysis of HD 111395 and, for the detrending, the polynomial orders of 2, 2, 0 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.9
continued. 
Fig. A.10
Plots for the time series analysis of HD 115043 and, for the detrending, the polynomial orders of 2, 2, 1 and 2 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.10
continued. 
Fig. A.11
Plots for the time series analysis of HD 118972 and, for the detrending, the polynomial orders of 0, 0, 1 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.11
continued. 
Fig. A.12
Plots for the time series analysis of HD 131977 and, for the detrending, the polynomial orders of 2, 2, 2 and 2 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 
Fig. A.12
continued. 
Fig. A.13
Periodograms for HD 35296 of the four time series without the removal of the most probable period (3.604 days) is depicted with the solid line and the data after the removal (O–C) with a dashed line. 
Fig. A.14
Periodograms for HD 35296 of the four time series without the removal of the 3.502 day period is depicted with the solid line and the data after the removal (O–C) with a dashed line. 
Fig. A.15
Periodograms for HD 39587 of the four time series without the removal of the most probable period (5.36 day) is depicted with the solid line and the data after the removal (O–C) with a dashed line. 
Fig. A.16
Periodograms for HD 39587 of the four time series without the removal of the 5.13 day period is depicted with the solid line and the data after removal (O–C) with a dashed line. 
Fig. A.17
Periodograms for HD 39587 of the four time series without the removal of the 5.6 day period is depicted with the solid line and the data after the removal (O–C) with a dashed line. 
Fig. A.18
Periodograms for HD 97334 of the four time series without the removal of the 7.92 day period is depicted with the solid line and the data after removal (O–C) with a dashed line. 
Fig. A.19
Periodograms for HD 97334 of the four time series without the removal of the 8.26 day period is depicted with the solid line and the data after removal (O–C) with a dashed line. 
Fig. A.20
Periodograms for HD 101501 of the four time series without the removal of the most probable period (16.06 day) is depicted with the solid line and the data after the removal (O–C) with a dashed line. 
All Tables
The mean Sindices of the Ca II H&K and Ca II IRT lines and the standard deviation (σ).
The correlation coefficient (Spearsman’s ρ) between S_{MWO} and S_{IRT}indices time series.
Ratio between the standard deviations of the periods and the mean period obtained from the different polynomial orders for the detrending.
Periods found using the S_{MWO} and S_{IRT} time series with the corresponding FAP.
All Figures
Fig. 1
A TIGRE spectrum of the Sun (moon light) of the Ca II IRT lines with the bandpasses chosen for the definition of the S_{IRT}indices: the top panel shows the spectral regions used for the S_{IRT_8498}index, the middle panel those used for the S_{IRT_8542}index and the bottom panel those used for the S_{IRT_8662}index. 

In the text 
Fig. 2
Upper panel: best fit coefficients of the constant A of the linear fit vs. B−V. The solid line represents the determined relation. Lower panel: the residuals. 

In the text 
Fig. 3
Upper panel: best fit coefficients of the slope B of the linear fit multiplied with 10^{5} vs. B−V. The solid line represents the determined relation. Lower panel: the residuals. 

In the text 
Fig. 4
Correlation between S_{IRT_8498} and S_{MWO} for HD 131977 with a Spearman’s ρ of 0.924. 

In the text 
Fig. 5
Spearman’s ρ correlation coefficients of the correlation between S_{MWO} and S_{IRT}indices vs. the standard deviation (σ) of the S_{MWO} time series. 

In the text 
Fig. 6
Example for the estimation of the polynomial order for the detrending of the S_{MWO} time series of HD 43162: the upper panel shows the maximal power p(ω), the panel below the results of the Ftest, where the significance of the scatter reduction is compared to the residuals of the nondetrended time series, the third panel shows the significance of the scatter reduction compared to detrending with the polynomial fit of order N1 and in the lowest panel we plot the estimated periods; the dotted lines show the formal 1σ limit. 

In the text 
Fig. 7
Averaged periodogram and individual periodograms from the S_{MWO} and the S_{IRT} time series for HD 43162; the solid line depicts the averaged periodogram, the dotted line the S_{MWO}, the dashed line the S_{IRT_8498}, the dasheddotted line the S_{IRT_8542} and the dashedtripledotted line the S_{IRT_8662} periodograms. 

In the text 
Fig. 8
Significance estimation of the S_{MWO} time series for HD 43162: Upper panel: logarithm of the probability (Prob(p>p_{0})) vs. logarithm of one minus significance level (1−p_{0}). The solid line shows the linear fit and the dotted line the quadratic fit. Lower panel: the residuals of the linear fit are depicted as crosses and the quadratic fit as asterisks. 

In the text 
Fig. 9
Significance estimation of the signal in the averaged periodogram for HD 111395: Upper panel: logarithm of the probability (Prob(p>p_{0})) vs. logarithm of one minus significance level (1−p_{0}). The solid line shows the linear fit and the dotted line the quadratic fit. Lower panel: the residuals of the linear fit are depicted as crosses and the quadratic fit as asterisks. 

In the text 
Fig. 10
Logarithm of the estimated FAPs of the derived periods in the individual time series vs. objects. The red bars show the log(FAP) values for the period derived from the S_{MWO} time series, the blue bars that from the S_{IRT_8498} time series, the yellow bars that from the S_{IRT_8542} and the green bars that from the S_{IRT_8662} time series. 

In the text 
Fig. 11
Comparison of the periods derived from the S_{IRT}indices time series vs. the periods derived from the S_{MWO} time series (in the upper panels) with the identity as solid line. The lower panels show the relative difference (in percent) between the periods derived in the S_{IRT}indices and S_{MWO} time series vs. the period obtained from the S_{MWO} time series with a 3σ error for the period ratios. The periods with a FAP less than 4.55 × 10^{2} (2σ) are labelled with a plus symbol and periods with a FAP greater than 4.55 × 10^{2} with an asterisk. 

In the text 
Fig. 12
Upper plot: comparison of our most probable rotation periods with the literature values with the identity as solid line. Lower panel: the relative difference in percent between our most probable rotation periods and the literature values with 3σ error for the period ratios. 

In the text 
Fig. 13
Periodograms for HD 39587 of the S_{MWO} and S_{IRT_8542} time series without the removal of the 5.36 day period (depicted as solid line) and after prewhitening (dashed line). 

In the text 
Fig. 14
Periodograms for HD 39587 of the S_{MWO} and S_{IRT_8542} time series without the removal of the 5.13 day period (depicted as solid line) and after prewhitening (dashed line). 

In the text 
Fig. 15
Periodograms for HD 39587 of the S_{MWO} and S_{IRT_8542} time series without the removal of the 5.6 day period (depicted as solid line) and after prewhitening (dashed line). 

In the text 
Fig. 16
Periodograms for HD 97334 of the S_{MWO} and S_{IRT_8498} time series without removal of the 7.92day period (depicted as solid line) and after prewhitening (dashed line). 

In the text 
Fig. 17
Periodograms for HD 97334 of the S_{MWO} and S_{IRT_8498} time series without removal of the 8.26day period (depicted as solid line) and after prewhitening (dashed line). 

In the text 
Fig. 18
Periodograms for HD 101501 of the S_{MWO} and S_{IRT_8662} time series without the removal of the 16.06 day period is depicted with the solid line and the data after the removal (O–C) with a dashed line. 

In the text 
Fig. A.1
Plots for the time series analysis of HD 35296 and, for the detrending, the polynomial orders of 3, 1, 1, and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542}, and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.1
continued. 

In the text 
Fig. A.2
Plots for the time series analysis of HD 37394, and, for the detrending, the polynomial orders of 1, 0, 0 and 1 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.2
continued. 

In the text 
Fig. A.3
Plots for the time series analysis of HD 39587 and, for the detrending, the polynomial orders of 2, 0, 1 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.3
continued. 

In the text 
Fig. A.4
Plots for the time series analysis of HD 43162 and, for the detrending, the polynomial orders of 1, 1, 1 and 1 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.4
continued. 

In the text 
Fig. A.5
Plots for the time series analysis of HD 72905 and, for the detrending, the polynomial orders of 0, 0, 1 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.5
continued. 

In the text 
Fig. A.6
Plots for the time series analysis of HD 75732 and, for the detrending, the polynomial orders of 3, 0, 0 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.6
continued. 

In the text 
Fig. A.7
Plots for the time series analysis of HD 97334 and, for the detrending, the polynomial orders of 2, 2, 0 and 1 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.7
continued. 

In the text 
Fig. A.8
Plots for the time series analysis of HD 101501 and, for the detrending, the polynomial orders of 0, 0, 0 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.8
continued. 

In the text 
Fig. A.9
Plots for the time series analysis of HD 111395 and, for the detrending, the polynomial orders of 2, 2, 0 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.9
continued. 

In the text 
Fig. A.10
Plots for the time series analysis of HD 115043 and, for the detrending, the polynomial orders of 2, 2, 1 and 2 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.10
continued. 

In the text 
Fig. A.11
Plots for the time series analysis of HD 118972 and, for the detrending, the polynomial orders of 0, 0, 1 and 0 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.11
continued. 

In the text 
Fig. A.12
Plots for the time series analysis of HD 131977 and, for the detrending, the polynomial orders of 2, 2, 2 and 2 are used for S_{MWO}, S_{IRT_8498}, S_{IRT_8542} and S_{IRT_8662}. A description of these figures can be found at the beginning of this Appendix. 

In the text 
Fig. A.12
continued. 

In the text 
Fig. A.13
Periodograms for HD 35296 of the four time series without the removal of the most probable period (3.604 days) is depicted with the solid line and the data after the removal (O–C) with a dashed line. 

In the text 
Fig. A.14
Periodograms for HD 35296 of the four time series without the removal of the 3.502 day period is depicted with the solid line and the data after the removal (O–C) with a dashed line. 

In the text 
Fig. A.15
Periodograms for HD 39587 of the four time series without the removal of the most probable period (5.36 day) is depicted with the solid line and the data after the removal (O–C) with a dashed line. 

In the text 
Fig. A.16
Periodograms for HD 39587 of the four time series without the removal of the 5.13 day period is depicted with the solid line and the data after removal (O–C) with a dashed line. 

In the text 
Fig. A.17
Periodograms for HD 39587 of the four time series without the removal of the 5.6 day period is depicted with the solid line and the data after the removal (O–C) with a dashed line. 

In the text 
Fig. A.18
Periodograms for HD 97334 of the four time series without the removal of the 7.92 day period is depicted with the solid line and the data after removal (O–C) with a dashed line. 

In the text 
Fig. A.19
Periodograms for HD 97334 of the four time series without the removal of the 8.26 day period is depicted with the solid line and the data after removal (O–C) with a dashed line. 

In the text 
Fig. A.20
Periodograms for HD 101501 of the four time series without the removal of the most probable period (16.06 day) is depicted with the solid line and the data after the removal (O–C) with a dashed line. 

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.