© ESO, 2017

1. Introduction

The present paper addresses the occurrence and properties of magnetic activity cycles in young solar-like stars with ages spanning the range 4−95 Myr. The study of activity cycles requires the availability of long-term time series of one or more magnetic activity indexes.

The solar magnetic activity has been well monitored in time by recording the measurements of different activity indexes such as the total solar irradiance (TSI), the sunspots number, the flare occurrence rate, and the intensity of specific spectral lines. The analysis of long-term time series of these activity proxies revealed that the Sun exhibits periodic or quasi-periodic variability phenomena occurring over different time scales ranging from 30 d to 80 yr. The 30-d periodicity is induced by the solar rotation that modulates the visibility of spots and faculae over the solar disk. An 11-yr period, the so-called Schwabe cycle, is related to the evolution of the solar magnetic field and is associated with a cyclic variation of the sunspots number and the average latitude at which spots and faculae occur.

The level of the solar magnetic activity is a complex function of time and other signals are superimposed over the 11-yr cycle.

A 154-d periodicity, for instance, was detected in γ-ray activity (Rieger et al. 1984), in the Mt. Wilson sunspot index (Ballester et al. 2002), and in the sunspot area (Lean 1990; Oliver et al. 1998). These 154-d cycles are usually named Rieger cycles because they were detected, for the first time, by Rieger et al. (1984).

The expression quasi-biennial variations is instead used to refer to periodic or quasi-periodic phenomena that occur with timescales ranging from 1 yr to 3 yr. Gnevyshev (1967, 1977), for instance, analyzed variations in sunspots numbers,the coccurence rate of large flares and coronal green-line emission and found that each solar cycle exhibits two maxima separated byintervals of 2–3 yr. Variability phenomena occurring in a 2-yr timescale have also been detected in the 35-yr TSI time series analyzed by Ferreira Lopes et al. (2015). A review of the quasi-biennial variations detected in solar activity and the physical processes invoked to explain such variations can be found in Hathaway (2015). Finally, on longer timescales, the Sun activity is characterized by the 22-yr Hale cycle (Hale et al. 1919) that is associated with a reverse in the polarity of sunspots and the 80-yr Gleisseberg cycle (Gleissberg 1958) that is a cyclic modulation of the amplitude of the 11-yr cycles.

Stars with a spectral type later than F5 have a magnetic activity similar to that observed in the Sun and these stars are therefore characterized by activity cycles similar to those occurring in our star. The analysis of long-term photometric time series has allowed the detection of activity cycles in a wide sample of late-type stars. The first studies on stellar activity cycles were based on the results of the Mount Wilson Ca II H&K survey (Wilson 1968) and were conducted by Wilson (1978), Durney et al. (1981), Baliunas & Vaughan (1985), and Baliunas et al. (1995). These authors measured the rotation period P_rot and cycle length P_cyc for tens of stars by analyzing the temporal variations of the photometric fluxes in two 0.1 nm passbands centered on the core of the Ca II H&K emission lines. Recently, Oláh et al. (2016) analyzed the full Ca II H&K time series collected at Mount Wilson for a sample of 29 stars. Their work exploits 36-yr time series that are 10 yr longer than those analyzed by Baliunas et al. (1995) and allow a better frequency resolution.

Messina & Guinan (2002) exploited long-term photometry in the Johnson V band to measure the rotation period P_rot and the length of activity cycle P_cyc in a sample of six young solar analogues. Oláh et al. (2009) analyzed long-term spectroscopic and photometric observations to study the activity cycles in a sample of 20 late-type stars.

Lehtinen et al. (2016) studied activity trends in a sample of 21 young solar-type stars by analyzing long-term photometry collected in Johnson B and V bands. Suárez Mascareño et al. (2016) detected long-term activity cycles in a sample of 49 nearby main-sequence stars by analyzing the V-band time series collected by the All Sky Automatic Survey (ASAS, Pojmanski 1997). Recently, high precision photometric time series acquired by space-borne telescopes were used to study stellar cycles. Vida et al. (2014), for instance, exploited Kepler data to study activity cycles in 39 fast-rotating late-type stars and Ferreira Lopes et al. (2015) studied activity cycles in 16 late-type stars by analyzing CoRoT time series.

Since stellar activity cycles have been detected, a link has been searched for between the cycle length P_cyc, magnetic activity indexes, and global stellar parameters such as the rotation period P_rot, Rossby Number Ro and age (see e.g. See et al. 2016; Oláh et al. 2016; García et al. 2014, and references therein). These relationships are important to probe the theoretical models developed to simulate the magnetic dynamos. The data currently available are still incomplete especially for young solar-like stars and the picture of the relationships between magnetic activity and global stellar parameters is far from being clear.

The relationship between the parameters P_cyc and P_rot is, for instance, not fully understood. Several authors found that P_cyc and P_rot are correlated with each other (see e.g. Baliunas et al. 1996; Böhm-Vitense 2007; Oláh et al. 2009, 2016). Other authors found no correlation between the two parameters (Lehtinen et al. 2016). Other works found a correlation in main-sequence F, G, and K stars and no correlation in M-type stars (Savanov 2012; Vida et al. 2014; Suárez Mascareño et al. 2016). This last feature has been interpreted as a difference in the dynamo mechanisms acting in FGK stars and M stars.

The study of the relationship between the ratio and the Rossby number showed that stars seem to lie in different activity branches (see e.g. Brandenburg et al. 1998; Lehtinen et al. 2016) possibly connected to different kinds of dynamos. Böhm-Vitense (2007) remarked that the Sun has an anomalous position with respect to these branches. It is unclear if the Sun is a special case or if other stars exhibit similar behavior.

The evolution of the ratio with the stellar age has been studied by Soon et al. (1993) and, more recently, by Oláh et al. (2016) in the age range 200–6200 Myr.

In the present work we searched for activity cycles in a sample of 90 late-type stars belonging to young loose stellar associations. In Distefano et al. (2016, hereafter Paper I) we investigated the rotational properties of these stars and determined a lower limit for their surface differential rotation (SDR). In the present paper we focused on their magnetic activity. As remarked in Paper I, the ages and spectral types of these stars span the ranges 4–95 Myr and G5-M4, respectively. This means that our sample comprises both fully convective stars and stars that have a radiative core surrounded by a convective envelope. Hence the data used here allow the characterization of activity cycles in stars with different structures and ages and extend the age range investigated by Soon et al. (1993) and Oláh et al. (2016).

The data used to perform our analysis and the procedure followed to measure the length of activity cycles are described in Sects. 2 and 3, respectively. The results of our analysis are reported in Sect. 4 and discussed in Sect. 5. In Sect. 6, the conclusions are drawn.

2. Data

2.1. Targets

In this paper we exploited long-term photometric time series collected by the ASAS survey and we searched for activity cycles in a sample of 90 solar-like stars belonging to young loose stellar associations. These targets are listed in Table 2 together with their main astrophysical parameters. In Table 1 we reported, for each association, the number of stars investigated and the number of stars in which we were able to detect one or more activity cycles. The same ASAS time series have been exploited in Paper I to investigate the rotational properties of our target stars.

In Paper I we also investigated 22 stars for which the SuperWASP (Butters et al. 2010) time series were available. These data have a photometric accuracy better than ASAS data. However, a typical SuperWASP dataset spans only a 4-yr interval and is characterized by observation gaps of several months. The ASAS time series, despite their lower photometric precision, are best suited to study long-term variability. Indeed these time series span an interval of about 9 yr and are characterized by shorter observation gaps due to the use of two different observing stations located at Las Campanas Observatory, Chile and at Haleakala, Maui, Hawaii, respectively. For this reason we decided to analyze the ASAS data exclusively.

2.2. Magnetic activity indexes

All works cited in the introduction measured the lengths of activity cycles by searching for a periodicity in the long-term time series of different activity indexes.

In the present paper, we exploited the long-term time series of three activity indexes that have been obtained from the analysis of the photometric data collected by the ASAS survey. In Paper I, we processed the ASAS data with a sliding window algorithm that divided each photometric time series in segments of a given length T. For each segment, we computed three activity indexes: the rotation period P_rot, inferred by the rotational modulation of the light curve by magnetically active regions (ARs); the median magnitude (V_med); and the interquantile range (IQR), defined as the difference between the 75th and the 25th percentiles of V magnitude values. Different examples of these activity index time series can be seen in Figs. 2–5 of Paper I and in Figs. 1–4, 12–13 of the present work.

Fig. 1 Vmed time series for the star ASAS J011315-6411.6. A cycle with length Pcyc = 1801 d was detected by the period-search algorithms. The black line is the sinusoid best-fitting the data.

Fig. 2 IQR time series for the star ASAS J070153-4227.9. A cycle with length Pcyc = 1818 d (red line) was detected by the period search algorithms. A visual inspection reveals also two shorter cycles with lengths 220 and 290 d, respectively. The dark and gray lines were obtained by fitting the data with smoothing cubic splines and were plotted to highlight the cycles detected by eye.

Fig. 3 Prot time series for the star ASAS J083656-7856.8. A visual inspection revealed two cycles with lengths 170 d and 130 d, respectively. Two smoothing cubic splines are overplotted to highlight the cycles with dark gray and light gray lines, respectively.

In Paper I we processed the ASAS photometric time series using three different values of the sliding-window length, i.e. T = 50 d, 100 d and 150 d, producing three different time series for each activity index. The results given in Paper I have been inferred by the time series obtained with T = 50 d. As discussed there, the use of sliding-window lengths of T = 100 d and T = 150 d tends to filter out the variability phenomena occurring over a timescale shorter than 100 d, flattening the amplitude variations of the different activity indexes (see. Figs. 2–5 of Paper I). However, the activity index time series obtained with T = 100 and 150 d are more regular and smoother than those obtained with T = 50 d and are more suitable to detect long-term variations due to activity cycles. Therefore, while T = 50 d is an optimal choice for analyzing SDR as in Paper I, 100-d segments, as used here, are better suited to search for activity cycles. This difference in T obviously does not introduce any inconsistencies in, for example, correlation analysis between activity cycles and SDR.

Fig. 4 Vmed time series for the star ASAS J082406-6334.1. The time series shows a long-term trend suggesting a cycle longer than 3000 d. A short cycle of 180 d is also visible and highlighted with a smoothing cubic spline in black.

In Paper I, we showed that the rotation period detected in different segments slightly changes in time. This trend could be induced by the combined effect of stellar SDR and migration of Ars in latitude. Indeed, ARs placed at different latitudes rotate with different frequencies in case of SDR. The migration of ARs, due to a Schwabe-like cycle, induces a temporal variation in the detected photometric period P_rot. Such a parameter can be therefore regarded as a tracer of the spot latitude migration.

This activity index has to be treated with caution. Lehtinen et al. (2016) and Lanza et al. (2016) noticed that the different rotation periods detected in different segments could also be a consequence of the ARs growth and decay occurring simultaneously at various latitudes and longitudes on the star. However, as remarked by Lehtinen et al. (2016), the use of segments with a short duration should limit this effect; see the discussion on the sliding-window length in Sect. 3.1 of Paper I and in Sect. 4.1 of Lehtinen et al. (2016) for more details. This activity index has recently been used by Vida et al. (2014) to investigate stellar cycles on Kepler data.

The V_med index is related to the axisymmetric part of the spot distribution on the stellar surface. In fact, the higher the fractional area covered by spots evenly distributed in longitude, the lower the average photometric flux in the V band. This is one of the most used activity index to investigate activity cycles (see e.g. Rodonò et al. 2000; Messina & Guinan 2002; Lehtinen et al. 2016)

The IQR index is proportional to the stellar variability amplitude. The variability amplitude is also an index widely used to study activity cycles (see e.g. Rodonò et al. 2000; Ferreira Lopes et al. 2015; Lehtinen et al. 2016). The variability of solar-like stars is due to different phenomena acting in different timescales such as the intrinsic evolution of spots and faculae, the intrinsic evolution of ARs complexes, and flux rotational modulation induced by ARs. In this paper, the IQR index was computed in 100-daylong segments. In these intervals the main source of variability is the rotational modulation (see e.g. Messina & Guinan 2002; Lanza et al. 2003). Hence the IQR index is proportional to the amplitude of the rotational modulation signal and is related to the non-axisymmetric part of the spot distribution.

Finally, for each star we also computed a global magnetic index that is given by the average value ⟨ IQR ⟩ of the IQR time series. The ⟨ IQR ⟩ values are reported in Table 2 for all the stars investigated here.

Fig. 5 Vmed time series for the AB Dor A. The gray squares are used to indicate the Vmed values computed in the 50-d segments.

2.3. Computation of the global parameters for the target stars

In this paper we investigate the relationships between the cycles lengths P_cyc, activity index ⟨ IQR ⟩, and global parameters of our targets. Our analysis was focused on the following parameters: the photometric shear ΔΩ_phot, effective temperature T_eff, convective turnover timescale τ_C, Rossby number Ro, radius at the bottom of the convection zone R_C, dynamo number D_N, and magnetic diffusivity timescale τ_diff. The values ΔΩ_phot, T_eff, τ_C, and Ro were computed in Paper I whereas R_C, D_N, and τ_diff were computed in this work for the first time.

The photometric shear of a given star is defined as (1)where P_min and P_max are the minimum and maximum values of the P_rot index and are interpreted as the rotation periods of two ARs occurring at two different latitudes. This parameter can be regarded as a lower limit for the SDR that is defined as (2)where Ω_Eq and Ω_pole are the stellar rotation frequencies at the equator and at the poles.

The effective temperature T_eff and convective turnover timescale τ_C have been derived in Paper I using Spada et al. (2013).

The Rossby number Ro is defined as (3)where τ_C is derived by the theoretical isochrones of Spada et al. (2013). Other works (Saar & Brandenburg 1999; Lehtinen et al. 2016) use the definition given by Brandenburg et al. (1998): (4)In these works, the convective turnover timescale is computed through the semi-empirical equation given by Noyes et al. (1984) that expresses τ_C as a function of the color index B−V. In order to make a comparison between our results and those of the abovementioned works, we computed the Ro values according to the Eq. (4) and using the τ_C values as given by the formula of Noyes et al. (1984). In the rest of the paper, the values computed by means of Eq. (4) are indicated with the symbol Ro_Br, where Br stands for Brandenburg.

The dynamo number is defined as (5)with (6)where D_CZ is the depth of the convective zone, l_m is the mixing length and H_D is the density scale height, η_t is the turbulent diffusivity, and ΔΩ the rotational shear. The parameters D_CZ, η_t, l_m, and H_D occurring in Eqs. (5) and (6) were derived by the theoretical models of Spada et al. (2013). The parameter ΔΩ was replaced with the parameter ΔΩ_phot computed in Paper I.

Finally, the magnetic diffusivity timescale is defined as (7)where η_t is the turbulent diffusivity; τ_diff is the timescale for the turbulent dissipation of magnetic energy over a length scale D_CZ. All the parameters computed for the different targets are reported in Table 2.

3. Method

We searched for activity cycles by running the Lomb-Scargle (Lomb 1976; Scargle 1982) and the phase dispersion minimization (PDM; Jurkevich 1971; Stellingwerf 1978) algorithms on the P_rot, V_med, and the IQR time series. The ASAS time series span, on average, an interval T = 3000 d. We assumed, as in Distefano et al. (2012), that a period P can be detected if P ≤ 0.75T and we performed our period search in the range (100–2000 d). We computed the False Alarm Probability (FAP) associated with a given period P_cyc according to the procedure described in Sect. 3 and we flagged a period as valid if it satisfied the requirement FAP < 0.1%. In some time series, visual inspection allowed the identification of cyclical patterns that the Lomb-Scargle and PDM algorithms were not able to detect. The failure of the period search algorithms is due to the quasi-periodic nature of the stellar cycles and to a variety of reasons: in some cases a cyclical pattern occurs only in a limited interval of the whole time series; in other cases a cyclical pattern is visible in the complete time series, but the length of the cycle changes in time. Finally, in some cases, the visual inspection suggests that the star is characterized by a cycle longer than the whole time series length. Since the size of our sample is relatively small, we decided to visually inspect each time series aiming at detecting and measuring the lengths of the cycles missed by the Lomb-Scargle and PDM algorithms. In this case, we flagged a P_cyc value as valid if at least one complete oscillation is visible (i.e., if at least two consecutive maxima or minima are clearly visible).

The uncertainty on P_cyc is estimated through the equation (8)where δν_cyc is the uncertainty associated with the frequency of the cycle ν_cyc = 1 /P_cyc. This uncertainty is given by the equation: (9)where δν_samp is the uncertainty due to the limited and discrete sampling of time series and δν_noise is the uncertainty due to data noise; δν_samp and δν_noise are computed according to Kovacs (1981) through the equations where T is the interval time spanned by the time series, σ is the standard deviation of the data before subtracting any periodic signal, N the number of time series points, and A the amplitude of the signal. In cases in which P_cyc values were determined by visual inspection we did not give a δP_cyc estimate since there is not an objective way to obtain it.

False alarm probability computation

The typical Monte Carlo method used to assess the FAP of a given period P₀, which has a power z₀ in the Lomb-Scargle periodogram, consists in simulating N synthetic time series with the same sampling of the original time series. The synthetic data points are usually generated by means of the equation (12)where ⟨ x ⟩ is the mean value of the original time series and R(0,σ) is a Gaussian random variate with a zero mean and a dispersion σ given by the standard deviation associated with the real data points. The Lomb-Scargle periodogram is computed for each time series and the highest peak Z_max is retained for each periodogram. The FAP associated with P₀ is then computed as the fraction of synthetic time series for which Z_max ≥ z₀.

The previous approach is valid only if the points of the original time series are uncorrelated, i.e., if two consecutive data points are independent of each other (see for example the discussion in Herbst & Wittenmyer 1996; Stassun et al. 1999; Rebull 2001). This is not the case for our time series, for which two consecutive elements of the activity indexes have been computed in segments that partially overlap with each other.

For this reason, we estimated the FAP by means of the Monte Carlo approach suggested by Brown et al. (1996) and also adopted by Rebull (2001) and Parihar et al. (2009). In this approach, the synthetic data points are generated by means of a recursive procedure (13)where x₀ and x_i are the first element and the ith element of the simulated time series and R(0,σ) denotes a Gaussian random variate centered on zero with dispersion σ (equal to the standard deviation of the real data). The parameters α and β are defined as (14)where Δt = t_i−t_i−1 is the time interval between the two consecutive points x_i−1 and x_i, and L_corr is the correlation timescale. The time series defined by Eqs. (13) is such that two points x_i and x_j satisfy the property (15)where C₂(i,j) is the two-points correlation function and Δt_ij = t_i−t_j the time between x_i and x_j (Brown et al. 1996). In our case, we adopted L_corr = T/ 3 were T is the length of the sliding window used to generate the activity index time series.

In order to assess the FAP probability associated with the period detected in a given time series we followed the Monte Carlo procedure described above by generating 10 000 synthetic time series. If the period was found by means of the Lomb-Scargle algorithm, the FAP was computed by taking the fraction of synthetic time series for which Z_max>z₀. If the period P₀ was found by means of the PDM algorithm, we ran the PDM on each synthetic time series and we retained the minimum value θ_min of the PDM periodogram for each time series. The FAP was then computed by taking the fraction of the synthetic time series for which θ₀<θ_min, where θ₀ is the PDM estimator associated with P₀.

4. Results

4.1. Results of time series analysis

We detected activity cycles in 67 stars. We found two distinct cycles in 32 of the stars and more than two cycles in 16 stars. The results of our analysis are reported in Table 3. In this table we listed the P_cyc lengths found for each star. For each cycle we reported a flag indicating the activity index from which P_cyc was inferred and a second flag indicating the method used to detect the period. If P_cyc was inferred by means of the Lomb-Scargle or the PDM algorithm, we also reported the associated FAP.

In Figs. 1–4 we report some examples of our analysis. In Fig. 1 we show the V_med time series for the star ASAS J011315-6411.6. The period search algorithms detected in this case a cycle with P_cyc = 1801 d. The sinusoid best fitting the data is overplotted (red continuous line).

In Fig. 2 we plot the IQR time series for the star ASAS J070153-4227.9. The Lomb-Scargle periodgram of this time series gives a highly significant peak at P = 1818 d. However, visual inspection also reveals the two shortest cycles with lengths of 220 and 290 d, respectively. These cycles could be secondary cycles analoguous to the Rieger cycles or to the quasi-biennal oscillations superimposed on the 11-yr solar cycle.

In Fig. 3 we report the P_rot time series for the star ASAS 083656-7856.8. The Lomb-Scargle and PDM algorithms were not able to detect any significant period in this star. However two cycles with lengths 170 d and 130 d are clearly visible in the data. These features are very common in the Sun and in solar-like stars. For instance, the Rieger cycles and the quasi-biennal oscillations have been well detected in the Sun but, as remarked by Oláh et al. (2016), they are not continuously present in the solar data registered until now. Oláh et al. (2016) noticed in the Mount Wilson time series that some of the detected cycles are only temporarily seen and that young stars are characterized by cycles whose duration changes in time. Rieger-like cycles were also detected by Lanza et al. (2009) in CoRoT-2 and by Bonomo & Lanza (2012) in Kepler-17.

Finally in Fig. 4, we report the V_med time series of the star ASAS J082406-6334.1. The visual inspection suggests the existence of a cycle longer than 3000 d. A shorter cycle with P_cyc = 180 d is also detected by eye and is highlighted with the black continuous line. Some of the P_cyc values detected by visual inspection are shorter than the sliding-window length used to process the ASAS time series. In fact, the 100-d sliding window attenuates the periodic signals shorter than 100 d but, in some cases, does not completely suppress them. So if these signal have a sufficient amplitude, they can still be detected after the segmentation procedure (see Appendix A for details).

4.2. Case of AB Dor A

To the best of our knowledge, the cycles reported in Table 3 have been detected here for the first time. The only exception in our sample is given by AB Dor A (HD 36705) which is identified with the ASAS ID J052845-6526.9. AB Dor A is a fast rotating (P_rot = 0.51 d) K1V star that is well studied in the literature and is the primary component of the quadruple system AB Dor, which also comprises the stars AB Dor Ba, AB Dor Bb, and AB Dor C. AB Dor Ba and AB Dor Bb are the components of the binary system AB Dor B, which was resolved for the first time by Janson et al. (2007) and is located at 8.9′′ ± 0.1′′ from AB Dor A (Martin & Brandner 1995). AB Dor C is a close companion of AB Dor A. It was detected by Guirado et al. (1997) and it is located at about 0.16′′ from AB Dor A. Guirado et al. (2010) made a dynamical estimate of the AB Dor A mass and obtained M = 0.86 ± 0.09 M_⊙, which is in good agreement with the photometric estimate M = 0.93 M_⊙ made in Paper I and based on the theoretical isochrones of Spada et al. (2013; see Paper I for details).

The magnetic activity of AB Dor A has been widely studied in the literature by means of spectroscopic and photometric data collected at different wavelengths (see e.g. Drake et al. 2015; Lalitha & Schmitt 2013; Budding et al. 2009; Lalitha & Schmitt 2013; Jeffers et al. 2007; Messina et al. 2006; Järvinen et al. 2005, and references therein)

Järvinen et al. (2005) merged the photometric data collected in several works and obtained a V-band time-serie spanning 1978–2000. These authors analyzed these data with the inversion technique developed by Berdyugina (1998) to study the temporal evolution of the spots longitude distribution. They analyzed the temporal variations of the spots longitudes and mean brightness by means of a Fourier analysis and detected a primary cycle with length P_cyc1 = 21 ± 3 yr and a secondary cycle with length P_cyc2 ~ 5.5 yr. The primary cycle is mainly associated with the mean magnitude variations, while the secondary cycle is a flip-flop cycle, i.e., a periodic switch of the longitude where the dominant spot concentration occurs. Järvinen et al. (2005) remarked that their estimate of the primary cycle length is not very accurate because of the sparseness of data in the years 1978–1985. We merged the photometric data reported in literature and summarized in Järvinen et al. (2005) with the ASAS data and we obtained a 27-yr V-band time series. We processed this longterm dataset with the sliding window algorithm described in Paper I and we obtained three 27-yr time series for the activity indexes V_med, P_rot, and IQR. In this specific case, we used a 50-d sliding window because the duration of the different observing seasons usually do not exceed two months and therefore, the use of a 100-d sliding window should be meaningless. The analysis of the V_med time series with the Lomb-Scargle algorithm revealed an activity cycle of length P_cyc = 16.78 ± 2 yr. In Fig. 5 we plot the long-term photometric time series and the sinusoid that best fits the data.

The use of the activity indexes time series also allowed the detection of secondary cycles with lengths of 400 d, 190 d, and 90 d, respectively.

Fig. 6 Pcyc length versus stellar rotation period Prot. The filled circles indicate the primary cycles and the empty circles the secondary cycles. The arrows indicate the stars for which the visual inspection of time series revealed Pcyc ≥ 3000 d. The continuous line, dashed line, and dotted line represent the Aa, Ab, and I sequences identified by Böhm-Vitense (2007), respectively. The dash-dotted line represents the locus where the ratio is equal to the ratio between the typical length of the Rieger cycle and the solar rotation period. The red asterisk indicates the solar Schwabe cycle; the green asterisk indicates the solar Rieger cycle.

5. Discussion

5.1. Relationships between cycle length, rotation period and Rossby number

Since stellar activity cycles were first detected, a relationship has been searched between the cycle period P_cyc and t stellar rotation period P_rot. Böhm-Vitense (2007) analyzed a sample of stars obtained by selecting the best quality spectrophotometric data collected at Mt. Wilson. She identified three different stellar sequences that have an almost constant ratio . The Aa and Ab sequences (where A stands for active) consist of young and active stars. The ratio is between 400 and 500 for the Aa stars and is about 300 for the Ab stars. The I sequence (where I stands for inactive) comprises old and less active stars for which the ratio is about 90. Böhm-Vitense (2007) speculates that the three sequences are due to different kind of dynamos excited by different phenomena. According to her interpretation, the dynamo should be generated by SDR in A sequence stars and by a vertical shear in I sequence stars.

In Fig. 6 we plot the cycle lengths versus stellar rotation periods as in Böhm-Vitense (2007). The filled circles indicate the primary cycles and the empty circles indicate the secondary cycles. We also overplotted the three sequences identified by Böhm-Vitense (2007) and a fourth line corresponding to the P_cyc/P_rot ≃ 5 ratio of the solar Rieger cycle. The figure shows that the primary cycles do not follow any particular trend and that their lengths seem to be uncorrelated with the rotation periods. The secondary cycles are also uniformly distributed between the Ab and I sequences and do not follow any particular pattern. Some of the secondary cycles lie beneath the I sequence and their P_cyc/P_rot ratio is very close to that of the solar Rieger-cycle.

Fig. 7 ratio versus inverse of the stellar rotation period. The filled circles indicate the primary cycles and the empty circles indicate the secondary cycles. The slope m = 1.02 ± 0.06 of the best-fitting straight line indicates a lack of correlation between Pcyc and Prot. The black asterisk indicates solar Schwabe cycle; the green asterisk indicates the solar Rieger cycle.

Other authors (see e.g. Baliunas et al. 1996; Vida et al. 2014; Oláh et al. 2016) searched for a correlation between and the parameter that, according to Baliunas et al. (1996), should be proportional to the dynamo number D. In these works a linear fit is performed investigating the relationship between the logarithms of the two quantities. The slope m extracted from the linear fit is the exponent of the power law P_cyc/P_rot = (1 /P_rot)^m. If m is close to 1 then no correlation exists between P_cyc and P_rot. Baliunas et al. (1996) found m = 0.74, which is in agreement with the value m = 0.76 recently found by Oláh et al. (2016).

In Fig. 7 we plotted versus for our targets. We performed a linear fit between the two quantities and we found m = 1.02 ± 0.06 that indicates no correlation between P_cyc and P_rot. Our results are slightly higher than those found by Baliunas et al. (1996) and Oláh et al. (2016) but are in close agreement with those found by Savanov (2012) and Lehtinen et al. (2016).

The lack of correlation between P_cyc and P_rot seen in Fig. 6 and in Fig. 7 suggests that the dynamo mechanism occurring in our young targets stars is different from that acting in the older Mt. Wilson stars selected by Böhm-Vitense (2007).

Fig. 8 Ratio of ratio versus RoBr number. The filled bullets indicate the primary cycles and the empty bullets indicate the secondary cycles. The black points and the black crosses indicate the results by Saar & Brandenburg (1999) and Lehtinen et al. (2016), respectively. The red continuous line depicts the inactive sequence. The black dotted line depicts the curve inferred by Lehtinen et al. (2016) and connects the active, the transitional and the super active sequence. The RoBr ranges corresponding to the three sequences are indicated with the black arrows and labelled A, T, and S, respectively.

Fig. 9 ⟨ IQR ⟩ versus parameters τdiff, Teff, DN, RC, and τC, respectively. In each panel the Spearman coefficient rS and the corresponding two-sided p-value are shown.

The lack of correlation between P_cyc and P_rot seen in our targets could be related to their low Rossby number values. In fact, most of our targets are fast rotating stars (P_rot< 10 d) with Ro values falling in the range (0.004–0.17) and Ro_Br values in the range (0.002–0.06). These ranges are very different from those covered by the Mt. Wilson stars. Saar & Brandenburg (1999) extended the sample of the Mt. Wilson stars with an ensemble of young and fast rotating stars and showed that these stars formed a third branch in the (, ) plane. This branch was called the super active branch and was clearly distinct from the active and inactive branches and was attributed to a different dynamo mechanism. Recently, Lehtinen et al. (2016) exploited long-term photometry to measure the cycle length P_cyc of 21 young active stars and found that the active and the super active branch are connected by a transitional branch. In Fig. 8 we plotted the values versus the values inferred for our stars. In the plot we also reported the data from Saar & Brandenburg (1999) and from Lehtinen et al. (2016). The primary cycles of our stars are in good agreement with the transitional branch as defined by Lehtinen et al. (2016). Also Lehtinen et al. (2016) did not find any correlation between P_cyc and P_rot in their targets. Hence, we can conclude that in stars belonging to the transitional branch, the cycle length P_cyc and the rotation period P_rot are uncorrelated.

5.2. Relationship between P_cyc and global stellar parameters

We investigated how the cycle lengths P_cyc are related to the global stellar parameters of our targets. We evaluated the degree of the correlation between P_cyc and a given parameter by computing the Spearman rank-order correlation coefficient (r_S) (Press et al. 1992). This coefficient is a nonparametric measure of the monotonicity of the relationship between two datasets.

A value of r_S close to 0 implies that the two datasets are poorly correlated, whereas the closer r_S is to ±1 the stronger the monotonic relationship between the two variables.

The statistical significance of a given value r_S = r₀ can be evaluated by computing the two-sided p-value, i.e., the probability (P( | r_S | ≥ | r₀ |) under the null hypothesys that two invesigated datasets are non-monotonically correlated (Press et al. 1992).

In Table 2 we reported the Spearman correlation coefficients between the length P_cyc of the primary cycles and the various stellar parameters. The corresponding p-values are also reported. The values of | r_S | are very close to 0 in all the cases indicating a poor correlation between P_cyc and the different variables.

5.3. Relationship between ⟨ IQR ⟩ and global stellar parameters

As remarked in Sect. 2, the IQR index is proportional to the amplitude of the rotational modulation signature and is related to the non-axisymmetric part of the spot distribution. The amplitude of rotational modulation is a widely used activity index and can be regarded as a robust proxy of the surface magnetic activity. García et al. (2013) demonstrated that the variance of the TSI, computed on a 60-d sliding window, closely mimics the variations of the 10.7 cm radio flux that in turn is a good indicator of the solar magnetic activity (see Bruevich et al. 2014, for details). We computed the Spearman coefficient between ⟨ IQR ⟩ and different stellar parameters to investigate how the mean level of the stellar surface magnetic activity is linked to the stellar properties. In Table 3 we reported the values of r_S for the different parameters. In this case, r_S values are significantly higher then those computed for P_cyc. A good correlation (| r_S | ≃ 0.5) is seen between ⟨ IQR ⟩ and T_eff, τ_C, τ_diff, R_C and D_N. The sign of correlation is positive for τ_C, τ_diff and D_N and negative for T_eff and R_C.

In Fig. 9 we reported ⟨ IQR ⟩ versus the various stellar parameters for which a significant correlation is seen (r_S ≃ 0.5). Despite the scatter in the data, the different plots clearly show that ⟨ IQR ⟩ tends to increase with τ_diff, τ_C, and T_eff and has a negative correlation with T_eff and R_C.

A weaker but still significant correlation is also seen between ⟨ IQR ⟩ and the photometric shear ΔΩ_phot. In Fig. 10 we plotted ⟨ IQR ⟩ versus ΔΩ_phot. The various stars are color-coded according to their effective temperature T_eff. The plot shows that

• ΔΩ_phottends to increase with temperature (as demonstrated in Paper I);

• the sign of the correlation between ⟨ IQR ⟩ and ΔΩ_phot and ⟨ IQR ⟩ is negative (r_S = −0.43). This means that the lower ΔΩ_phot the higher ⟨ IQR ⟩.

These results are in very good agreement with the theoretical models developed by Kitchatinov & Olemskoy (2011). In fact, these models predict that even a small SDR is very efficient for dynamos in M-type stars, it is less efficient in K and G type and even a strong SDR is inefficient in F-type stars.

Fig. 10 ⟨ IQR ⟩ versus ΔΩphot. The points are color coded according the stellar effective temperature Teff.

5.4. Evolution of P_cyc and ⟨ IQR ⟩ with the stellar age

Oláh et al. (2016) investigated how the ratio depends on the stellar age. They found that exhibits a large scatter in the younger and more active stars of their sample whereas the older stars have about the same value. The age at which becomes constant is, according to their work, at about 2.2 Gyr, i.e., the age of the Vaughan-Preston gap (Vaughan & Preston 1980). The sample of stars analyzed by Oláh et al. (2016) covers the age range 200–6200 Myr. Our work permits us to extend the age range studied by Oláh et al. (2016) and to investigate the transition age between the pre-main sequence (PMS) and the main sequence (MS) phase. Our analysis is slightly different from that performed by Oláh et al. (2016). We decided to investigate the relationship between P_cyc and the stellar age instead of the relationship between and the stellar age. In fact, the trend of the between ratio, in the age range covered by our targets, could be dominated by the complex P_rot evolution seen in the PMS stars (see e.g. Lanzafame & Spada 2015, and references therein). In the top panel of Fig. 11, we reported the length of the primary cycles P_cyc versus the stellar age. The blue filled symbols indicate the cycles identified in the present work and the red empty symbols mark those found by Oláh et al. (2016). The picture shows that P_cyc is about constant in the age range 4−300 Myr. After 300 Myr, the cycle lengths are more scattered and seem to increase with the stellar age. After 2.2 Gyr, that is the age corresponding to the Vaughan-Preston gap, the scatter in P_cyc decreases as described by Oláh et al. (2016). Although our data cover an interval of about 3000 d, which could prevent the detection of longer cycles, we observed long-term trends suggesting cycles longer than 3000 d in just 11 of the 90 investigated targets. For this reason, we are confident that our data are not seriously affected by a selection bias and that the age range investigated here is really, on average, characterized by cycles shorter than those observed in the older Mt. Wilson stars. In the bottom panel of Fig. 11 we plotted ⟨ IQR ⟩ versus the stellar age. The picture shows that ⟨ IQR ⟩ tends to decrease with the stellar age as also indicated by the Spearman coefficient r_S = 0.39. This result is in agreement with all the works that studied the evolution of the level of magnetic activity versus age (see e.g. Žerjal et al. 2017; Oláh et al. 2016, and references therein)

Fig. 11 Top panel: Pcyc versus the stellar age. The blue filled circles refer to the primary cycles detected in the present work and the red empty circles to the data reported by Oláh et al. (2016). The blue arrows indicate the stars for which the visual inspection of time series revealed Pcyc ≥ 3000 d. Bottom panel: activity index ⟨ IQR ⟩ versus the stellar age is shown.

5.5. Butterfly diagrams

Fig. 12 Time series for star ASAS J070030-7941.8. Top panel: Vmed time series for the star ASAS J070030-7941.8 is shown. Middle panel: Prot time series is shown. Bottom panel: IQR time series is shown. A primary cycle with length Pcyc = 1834 d was detected in the Vmed time series and the sinusoid best fitting the data was overplotted on it (red dashed line). The Prot and the IQR time series exhibit cycles (dark continuous lines) shorter than the primary cycle. The black dotted lines indicate the times at which maxima occur in the Prot time series.

In the Sun, the latitudes at which ARs occur change versus time according to the so-called “butterfly-diagram”. At the beginning of a solar cycle, ARs emerge in belts located at intermediate latitudes (± 35deg). As the cycle progresses, the ARs formation regions migrate toward the equator. At the beginning of the next cycle, the ARs form again at intermediate latitudes. The solar SDR combined with the ARs migration determines a decrease in the rotation period detectable along the 11-yr cycle because the ARs located at intermediate latitudes rotate more slowly than those located at the equator. As remarked in Paper I and in Sect. 2.2, the activity index P_rot of a given star can therefore be regarded as a tracer of the ARs migration in latitude. This activity index cannot give information about the exact latitudes at which ARs occur but can give insight into the nature of stellar magnetic field. In fact, in a star with a solar-like dynamo, we expect that the trend of P_rot with time mimics that of the solar butterfly diagram, i.e., that during a given cycle P_rot decreases in time and that it rises rapidly to higher values at the beginning of the next cycle. Messina & Guinan (2003), for instance, analyzed the trend of P_rot in six young solar analogues. They noticed that three of these stars follow a solar behavior whereas the other three show an anti-solar behaviour with P_rot increasing during the cycle. The P_rot time series of our target stars, in general, have fewer points than the V_med and IQR time series. Indeed, as remarked in Paper I, P_rot is detectable in a given time interval only if ARs maintain a stable configuration. In Figs. 12 and 13 we reported the P_rot time series for the stars ASAS J070030-7941.8 and ASAS J055329-8156.9 which are the targets with the highest number of P_rot determinations. The V_med time series and IQR time series of the two stars are also reported for comparison. The stars have two primary cycles with length P_cyc = 1834 d and P_cyc = 1290 d, respectively. The two sinusoids that best fit the data and correspond to the primary cycles are overplotted on the V_med time series. The visual inspection of the pictures shows that in both stars:

• the P_rot index oscillates with a cycle shorter than the primary cycle;hence the ARs latitude migration occurs in a timescale shorterthan the length of the primary cycle;

• in some intervals, the trend of P_rot vs time mimics the trends of V_med and IQR and the maxima and minima of P_rot correspond to local maxima and minima of the V_med and IQR time series; in other intervals the maxima and minima of P_rot are uncorrelated with those of the V_med and IQR time series. This implies that in some time intervals the changes in AR latitudes are correlated with the variations in the ARs areas and in the mean level of magnetic activity whereas in some other intervals there is no correlation;

• the P_rot variations with time are very different from those occurring in the Sun. Rise patterns are followed by decreasing patterns. This implies that ARs first migrate from regions with shorter rotation periods to regions with higher rotation periods and then migrate in the opposite sense. In the Sun the migration occurs in only one sense i.e. from the intermediate latitudes to the equator.

Similar trends are observed in all the targets analyzed here.

Fig. 13 Time series for star ASAS J055329-8156.9. Top panel: the Vmed time series for the star ASAS J055329-8156.9 is shown. Middle panel: the Prot time series is shown. Bottom panel: the IQR time series is shown. A primary cycle with length Pcyc = 1290 d was detected in the Vmed time series and the sinusoid best fitting the data was overplotted on it (red dashed line). The Prot and IQR time series exhibit cycles (dark continuous lines) shorter than the primary cycle. The black dotted lines indicate the times at which maxima and minima occur in the Prot time series.

The patterns seen in the P_rot time series of our targets can be regarded as a further evidence that the dynamo acting in these stars is very different from that acting in the Sun.

6. Conclusions

In the present work we searched for activity cycles in stars belonging to young loose stellar associations covering the age range 4−95 Myr. We investigated the correlation between cycle properties and global stellar parameters and how cycle properties evolve with the stellar age. In particular, our work extends the age range covered by the Mt Wilson stars whose properties have recently been reanalyzed by Oláh et al. (2016).

We analyzed the long-term time series of three different activity indexes and we were able to detect activity cycles in 67 stars and to measure their length P_cyc. About half of these stars show multiple and complex cycles according to the results found by Oláh et al. (2016) for young and active stars. Some of the detected secondary cycles have a ratio similar to that of the solar Rieger-cycles .

We investigated how P_cyc is correlated with the photometric shear ΔΩ_phot, effective temperature T_eff, convective turnover timescale τ_C, Rossby number Ro, radius at the bottom of the convection zone R_C, dynamo number D_N, and magnetic diffusivity timescale τ_diff. We found that P_cyc is essentially uncorrelated with the abovementioned parameters in our targets. The lack of correlation between P_cyc and P_rot suggests that the dynamo acting in our targets is different from that acting in the active and inactive branches determined by Böhm-Vitense (2007). Indeed, the location of our targets in the plane is in good agreement with the transitional branch as determined by Lehtinen et al. (2016). The activity index ⟨ IQR ⟩, which can be regarded as a proxy of the magnetic surface activity level, is positively correlated with τ_C, τ_diff, D_N and negatively correlated with T_eff and R_C. The analysis of the relationship between ⟨ IQR ⟩, SDR, and T_eff supports the theoretical models developed by Kitchatinov & Olemskoy (2011). Indeed, in agreement with these models, even a small differential rotation is efficient for dynamos in M-type stars, but it becomes less efficient or completely inefficient at increasing T_eff (cf. Fig. 10).

The analysis of the butterfly diagrams of our target stars shows two main features. The first is that the ARs migration in latitude occurs over a timescale shorter than the primary cycle; the second is that the latitudes at which ARs emerge seem to oscillate from a maximum to a minimum latitude and vice-versa. This is very different from the solar behaviour where the spots migrate only from intermediate latitudes to the equator. We merged our data with those analyzed by Oláh et al. (2016) and we studied how P_cyc evolve with the stellar age. We found that P_cyc is about constant and does not show significant correlation with the stellar age in the range (4−300 Myr); after 300 Myr P_cyc values are quite scattered and tend to increase with the stellar age; and after 2.2 Gyr, P_cyc values tend to be less scattered and seem to converge to the solar value as described by Oláh et al. (2016). The activity index ⟨ IQR ⟩ decreases with the stellar age in the age range 4–95 Myr. Finally, we merged the ASAS time series of the star AB Dor A with the photometric data from previous works and we detected a cycle with length P_cyc = 16.78 ± 2 yr and shorter secondary cycles with lengths of 400 d, 190 d, and 90 d.

Acknowledgments

The authors are grateful to the referee Lauri Jetsu for helpful comments and suggestions.

References

Appendix A: Effect of the sliding-window algorithm on a periodic signal

The segmentation procedure used to derive activity indexes is equivalent to a low-pass filter. This preserves the signals with a frequency lower than a certain cutoff frequency and attenuates signals with frequencies higher than the cutoff frequency. In our case, the use of a sliding-window with length T = 100 d attenuates signals with periods shorter than 100-d but, in some cases, does not completely suppress them. We performed different tests by simulating sinusoidal signals with different amplitudes and periods and by processing them with our segmentation algorithm. In the top panel of Fig. A.1 we plotted a simulated time series obtained by combining two sinusoidal signals with periods P₁ = 300 d, P₂ = 45 d and amplitudes A₁ = 0.4, A₂ = 0.6, respectively. A white gaussian noise with variance σ = 0.1 was added to the simulated data. In the bottom panel, we plotted the filtered signal obtained after processing the simulated time series with our segmentation algorithm. The filter has attenuated but not suppressed the 45-d signal. In Fig. A.2 we reported a similar test obtained by combining two period signals P₁ = 300 d, P₂ = 45 d and amplitudes A₁ = 0.4, A₂ = 0.1. In this second case the signal with P = 45 d is completely suppressed.

Fig. A.1 Top panel: a simulated time series obtained by combining two sinusoidal signals with periods P1 = 300 d, P2 = 45 d and amplitudes A1 = 0.4, A2 = 0.6, respectively. Bottom panel: filtered time series. The segmentation algorithm attenuates but does not suppress the signal with P = 45 d.

Fig. A.2 Top panel: a simulated time series obtained by combining two period signals with periods P1 = 300 d, P2 = 45 d and amplitudes A1 = 0.4, A2 = 0.1, respectively. Bottom panel: filtered time series. The 45-d is completely suppressed in this case.

All Tables

All Figures

	Fig. 1 Vmed time series for the star ASAS J011315-6411.6. A cycle with length Pcyc = 1801 d was detected by the period-search algorithms. The black line is the sinusoid best-fitting the data.
In the text

	Fig. 2 IQR time series for the star ASAS J070153-4227.9. A cycle with length Pcyc = 1818 d (red line) was detected by the period search algorithms. A visual inspection reveals also two shorter cycles with lengths 220 and 290 d, respectively. The dark and gray lines were obtained by fitting the data with smoothing cubic splines and were plotted to highlight the cycles detected by eye.
In the text

	Fig. 3 Prot time series for the star ASAS J083656-7856.8. A visual inspection revealed two cycles with lengths 170 d and 130 d, respectively. Two smoothing cubic splines are overplotted to highlight the cycles with dark gray and light gray lines, respectively.
In the text

	Fig. 4 Vmed time series for the star ASAS J082406-6334.1. The time series shows a long-term trend suggesting a cycle longer than 3000 d. A short cycle of 180 d is also visible and highlighted with a smoothing cubic spline in black.
In the text

	Fig. 5 Vmed time series for the AB Dor A. The gray squares are used to indicate the Vmed values computed in the 50-d segments.
In the text

Fig. 6 Pcyc length versus stellar rotation period Prot. The filled circles indicate the primary cycles and the empty circles the secondary cycles. The arrows indicate the stars for which the visual inspection of time series revealed Pcyc ≥ 3000 d. The continuous line, dashed line, and dotted line represent the Aa, Ab, and I sequences identified by Böhm-Vitense (2007), respectively. The dash-dotted line represents the locus where the ratio is equal to the ratio between the typical length of the Rieger cycle and the solar rotation period. The red asterisk indicates the solar Schwabe cycle; the green asterisk indicates the solar Rieger cycle.

In the text

	Fig. 7 ratio versus inverse of the stellar rotation period. The filled circles indicate the primary cycles and the empty circles indicate the secondary cycles. The slope m = 1.02 ± 0.06 of the best-fitting straight line indicates a lack of correlation between Pcyc and Prot. The black asterisk indicates solar Schwabe cycle; the green asterisk indicates the solar Rieger cycle.
In the text

Fig. 8 Ratio of ratio versus RoBr number. The filled bullets indicate the primary cycles and the empty bullets indicate the secondary cycles. The black points and the black crosses indicate the results by Saar & Brandenburg (1999) and Lehtinen et al. (2016), respectively. The red continuous line depicts the inactive sequence. The black dotted line depicts the curve inferred by Lehtinen et al. (2016) and connects the active, the transitional and the super active sequence. The RoBr ranges corresponding to the three sequences are indicated with the black arrows and labelled A, T, and S, respectively.

In the text

	Fig. 9 ⟨ IQR ⟩ versus parameters τdiff, Teff, DN, RC, and τC, respectively. In each panel the Spearman coefficient rS and the corresponding two-sided p-value are shown.
In the text

	Fig. 10 ⟨ IQR ⟩ versus ΔΩphot. The points are color coded according the stellar effective temperature Teff.
In the text

	Fig. 11 Top panel: Pcyc versus the stellar age. The blue filled circles refer to the primary cycles detected in the present work and the red empty circles to the data reported by Oláh et al. (2016). The blue arrows indicate the stars for which the visual inspection of time series revealed Pcyc ≥ 3000 d. Bottom panel: activity index ⟨ IQR ⟩ versus the stellar age is shown.
In the text

Fig. 12 Time series for star ASAS J070030-7941.8. Top panel: Vmed time series for the star ASAS J070030-7941.8 is shown. Middle panel: Prot time series is shown. Bottom panel: IQR time series is shown. A primary cycle with length Pcyc = 1834 d was detected in the Vmed time series and the sinusoid best fitting the data was overplotted on it (red dashed line). The Prot and the IQR time series exhibit cycles (dark continuous lines) shorter than the primary cycle. The black dotted lines indicate the times at which maxima occur in the Prot time series.

In the text

Fig. 13 Time series for star ASAS J055329-8156.9. Top panel: the Vmed time series for the star ASAS J055329-8156.9 is shown. Middle panel: the Prot time series is shown. Bottom panel: the IQR time series is shown. A primary cycle with length Pcyc = 1290 d was detected in the Vmed time series and the sinusoid best fitting the data was overplotted on it (red dashed line). The Prot and IQR time series exhibit cycles (dark continuous lines) shorter than the primary cycle. The black dotted lines indicate the times at which maxima and minima occur in the Prot time series.

In the text

	Fig. A.1 Top panel: a simulated time series obtained by combining two sinusoidal signals with periods P1 = 300 d, P2 = 45 d and amplitudes A1 = 0.4, A2 = 0.6, respectively. Bottom panel: filtered time series. The segmentation algorithm attenuates but does not suppress the signal with P = 45 d.
In the text

	Fig. A.2 Top panel: a simulated time series obtained by combining two period signals with periods P1 = 300 d, P2 = 45 d and amplitudes A1 = 0.4, A2 = 0.1, respectively. Bottom panel: filtered time series. The 45-d is completely suppressed in this case.
In the text